EEMtoolbox: A User‐friendly R Package for Flexible Ensemble Ecosystem Modelling

Forecasting ecosystem changes due to disturbances or conservation interventions is essential to improve ecosystem management and anticipate unintended consequences of conservation decisions. Mathematical models allow practitioners to understand the potential effects and unintended consequences via simulation. However, calibrating these models is often challenging due to a paucity of appropriate ecological data. Ensemble ecosystem modelling (EEM) is a quantitative method used to parameterize models from theoretical ecosystem features rather than data. Two approaches have been considered to find parameter values satisfying those features: a standard accept–reject algorithm, appropriate for small ecosystem networks, and a sequential Monte Carlo (SMC) algorithm that is more computationally efficient for larger ecosystem networks. In practice, using SMC for EEM generation requires advanced statistical and mathematical knowledge, as well as strong programming skills, which might limit its uptake. In addition, current EEM approaches have been developed for only one model structure (generalised Lotka–Volterra). To facilitate the usage of EEM methods, we introduce EEMtoolbox, an R package for calibrating quantitative ecosystem models. Our package allows the generation of parameter sets satisfying ecosystem features by using either the standard accept–reject algorithm or the novel SMC procedure. Our package extends the existing EEM methodology, originally developed for the generalised Lotka–Volterra model, to two additional model structures (the multispecies Gompertz and the Bimler–Baker model) and additionally allows users to define their own model structures. We demonstrate the usage of EEMtoolbox by simulating changes in species abundance immediately after the release of the sihek ( Todiramphus cinnamominus , extinct‐in‐the‐wild species) on Palmyra Atoll in the Pacific Ocean. With its simple interface, our package facilitates straightforward generation of EEM parameter sets, thus unlocking advanced statistical methods supporting conservation decisions using ecosystem network models.