The European Regional Seas Ecosystem Model (ERSEM) is one of the most established ecosystem models for the lower trophic levels of the marine food web in the scientific literature. Since its original development in the early nineties it has evolved significantly from a coastal ecosystem model for the North Sea to a generic tool for ecosystem simulations from shelf seas to the global ocean. The current model release contains all essential elements for the pelagic and benthic parts of the marine ecosystem, including the microbial food web, the carbonate system, and calcification. Its distribution is accompanied by a testing framework enabling the analysis of individual parts of the model. Here we provide a detailed mathematical description of all ERSEM components along with case studies of mesocosm-type simulations, water column implementations, and a brief example of a full-scale application for the north-western European shelf. Validation against in situ data demonstrates the capability of the model to represent the marine ecosystem in contrasting environments.

This paper presents the first decadal reanalysis simulation of the biogeochemistry of the North West European shelf, along with a full evaluation of its skill, confidence, and value. An error‐characterized satellite product for chlorophyll was assimilated into a physical‐biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The results showed that the reanalysis improved the model simulation of assimilated chlorophyll in 60% of the study region. Model validation metrics showed that the reanalysis had skill in matching a large data set of in situ observations for 10 ecosystem variables. Spearman rank correlations were significant and higher than 0.7 for physical‐chemical variables (temperature, salinity, and oxygen), ∼0.6 for chlorophyll and nutrients (phosphate, nitrate, and silicate), and significant, though lower in value, for partial pressure of dissolved carbon dioxide (∼0.4). The reanalysis captured the magnitude of pH and ammonia observations, but not their variability. The value of the reanalysis for assessing environmental status and variability has been exemplified in two case studies. The first shows that between 325,000 and 365,000 km² of shelf bottom waters were vulnerable to oxygen deficiency potentially threatening bottom fishes and benthos. The second application confirmed that the shelf is a net sink of atmospheric carbon dioxide, but the total amount of uptake varies between 36 and 46 Tg C yr⁻¹ at a 90% confidence level. These results indicate that the reanalysis output data set can inform the management of the North West European shelf ecosystem, in relation to eutrophication, fishery, and variability of the carbon cycle.

A reanalysis is a physically consistent set of optimally merged simulated model states and historical observational data, using data assimilation. High computational costs for modeled processes and assimilation algorithms has led to Earth system specific reanalysis products for the atmosphere, the ocean and the land separately. Recent developments include the advanced uncertainty quantification and the generation of biogeochemical reanalysis for land and ocean. Here, we review atmospheric and oceanic reanalyzes, and more in detail biogeochemical ocean and terrestrial reanalyzes. In particular, we identify land surface, hydrologic and carbon cycle reanalyzes which are nowadays produced in targeted projects for very specific purposes. Although a future joint reanalysis of land surface, hydrologic, and carbon processes represents an analysis of important ecosystem variables, biotic ecosystem variables are assimilated only to a very limited extent. Continuous data sets of ecosystem variables are needed to explore biotic‐abiotic interactions and the response of ecosystems to global change. Based on the review of existing achievements, we identify five major steps required to develop terrestrial ecosystem reanalysis to deliver continuous data streams on ecosystem dynamics.

We assimilated phytoplankton functional types (PFTs) derived from ocean color into a marine ecosystem model, to improve the simulation of biogeochemical indicators and emerging properties in a shelf sea. Error‐characterized chlorophyll concentrations of four PFTs (diatoms, dinoflagellates, nanoplankton, and picoplankton), as well as total chlorophyll for comparison, were assimilated into a physical‐biogeochemical model of the North East Atlantic, applying a localized Ensemble Kalman filter. The reanalysis simulations spanned the years 1998–2003. The skill of the reference and reanalysis simulations in estimating ocean color and in situ biogeochemical data were compared by using robust statistics. The reanalysis outperformed both the reference and the assimilation of total chlorophyll in estimating the ocean‐color PFTs (except nanoplankton), as well as the not‐assimilated total chlorophyll, leading the model to simulate better the plankton community structure. Crucially, the reanalysis improved the estimates of not‐assimilated in situ data of PFTs, as well as of phosphate and pCO₂, impacting the simulation of the air‐sea carbon flux. However, the reanalysis increased further the model overestimation of nitrate, in spite of increases in plankton nitrate uptake. The method proposed here is easily adaptable for use with other ecosystem models that simulate PFTs, for, e.g., reanalysis of carbon fluxes in the global ocean and for operational forecasts of biogeochemical indicators in shelf‐sea ecosystems.