EEMtoolbox: A user-friendly R package for flexible ensemble ecosystem modeling

Abstract

O_LIForecasting 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. C_LIO_LIEnsemble 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 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 (generalized Lotka-Volterra). C_LIO_LITo 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 generalized Lotka-Volterra model, to two additional model structures (the multi-species Gompertz, and the Bimler-Baker model), and additionally allows users to define their own model structures. C_LIO_LIWe demonstrate the usage of EEMtoolbox by modelling the introduction of sihek (extinct-in-the-wild) on Palmyra Atoll in the Pacific Ocean. With its simple interface, our package facilitates straightforward generation of EEM parameter sets, thus unlocks advanced statistical methods supporting conservation decisions using ecosystem network models. C_LI