Abstract. 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.

A carbon budget for the northwest European continental shelf seas (NWES) was synthesized using available estimates for coastal, pelagic and benthic carbon stocks and flows. Key uncertainties were identified and the effect of future impacts on the carbon budget were assessed. The water of the shelf seas contains between 210 and 230 Tmol of carbon and absorbs between 1.3 and 3.3 Tmol from the atmosphere annually. Offshelf transport and burial in the sediments account for 60-100 and 0-40% of carbon outputs from the NWES, respectively. Both of these fluxes remain poorly constrained by observations and resolving their magnitudes and relative importance is a key research priority. Pelagic and benthic carbon stocks are dominated by inorganic carbon. Shelf sediments contain the largest stock of carbon, with between 520 and 1600 Tmol stored in the top 0.1 m of the sea bed. Coastal habitats such as salt marshes and mud flats contain large amounts of carbon per unit area but their total carbon stocks are small compared to pelagic and benthic stocks due to their smaller spatial extent. The large pelagic stock of carbon will continue to increase due to the rising concentration of atmospheric CO 2 , with associated pH decrease. Pelagic carbon stocks and flows are also likely to be significantly affected by increasing acidity and temperature, and circulation changes but the net impact is uncertain. Benthic carbon stocks will be affected by increasing temperature and acidity, and decreasing oxygen concentrations, although the net impact of these interrelated changes on carbon stocks is uncertain and a major knowledge gap. The impact of bottom trawling on benthic carbon stocks Frontiers in Marine Science | www.frontiersin.org 1 March 2020 | Volume 7 | Article 143Legge et al.Carbon on the Northwest European Shelf is unique amongst the impacts we consider in that it is widespread and also directly manageable, although its net effect on the carbon budget is uncertain. Coastal habitats are vulnerable to sea level rise and are strongly impacted by management decisions. Local, national and regional actions have the potential to protect or enhance carbon storage, but ultimately global governance, via controls on emissions, has the greatest potential to influence the long-term fate of carbon stocks in the northwestern European continental shelf.

Abstract. The ERSEM model 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 part 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-West European shelf. Validation against in situ data demonstrates the capability of the model to represent the marine ecosystem in contrasting environments.

The benthic environment is a crucial component of marine systems in the provision of ecosystem services, sustaining biodiversity and in climate regulation, and therefore important to human society. With the contemporary increase in computational power, model resolution and technological improvements in quality and quantity of benthic data, it is necessary to ensure that benthic systems are appropriately represented in coupled benthic-pelagic biogeochemical and ecological modelling studies. In this paper we focus on five topical challenges related to various aspects of modelling benthic environments: organic matter reactivity, dynamics of benthic-pelagic boundary layer, microphytobenthos, biological transport and small-scale heterogeneity, and impacts of episodic events. We discuss current gaps in their understanding and indicate plausible ways ahead. Further, we propose a three-pronged approach for the advancement of benthic and benthic-pelagic modelling, essential for improved understanding, management and prediction of the marine environment. This includes: (A) development of a traceable and hierarchical framework for benthic-pelagic models, which will facilitate integration among models, reduce risk of bias, and clarify model limitations; (B) extended cross-disciplinary approach to promote effective collaboration between modelling and empirical scientists of various backgrounds and better involvement of stakeholders and end-users; (C) a common vocabulary for terminology used in benthic modelling, to promote model development and integration, and also to enhance mutual understanding.

A novel pan-European marine model ensemble was established, covering nearly all seas under the regulation of the Marine Strategy Framework Directive (MSFD), with the aim of providing a consistent assessment of the potential impacts of riverine nutrient reduction scenarios on marine eutrophication indicators. For each sea region, up to five coupled biogeochemical models from institutes all over Europe were brought together for the first time. All model systems followed a harmonised scenario approach and ran two simulations, which varied only in the riverine nutrient inputs. The load reductions were evaluated with the catchment model GREEN and represented the impacts due to improved management of agriculture and wastewater treatment in all European river systems. The model ensemble, comprising 15 members, was used to assess changes to the core eutrophication indicators as defined within MSFD Descriptor 5. In nearly all marine regions, riverine load reductions led to reduced nutrient concentrations in the marine environment. However, regionally the nutrient input reductions led to an increase in the non-limiting nutrient in the water, especially in the case of phosphate concentrations in the Black Sea. Further core eutrophication indicators, such as chlorophyll-a, bottom oxygen and the Trophic Index TRIX, improved nearly everywhere, but the changes were less pronounced than for the inorganic nutrients. The model ensemble displayed strong consistency and robustness, as most if not all models indicated improvements in the same areas. There were substantial differences between the individual seas in the speed of response to the reduced nutrient loads. In the North Sea ensemble, a stable plateau was reached after only three years, while the simulation period of eight years was too short to obtain steady model results in the Baltic Sea. The ensemble exercise confirmed the importance of improved management of agriculture and wastewater treatments in the river catchments to reduce marine eutrophication. Several shortcomings were identified, the outcome of different approaches to compute the mean change was estimated and potential improvements are discussed to enhance policy support. Applying a model ensemble enabled us to obtain highly robust and consistent model results, substantially decreasing uncertainties in the scenario outcome.