ERSEM 15.06: a generic model for marine biogeochemistry and the ecosystem dynamics of the lower trophic levels

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.

Marine regime shifts in ocean biogeochemical models: a case study in the Gulf of Alaska

Regime shifts have been reported in many marine ecosystems, and are often expressed as an abrupt change occurring in multiple physical and biological components of the system. In the Gulf of Alaska, a regime shift in the late 1970s was observed, indicated by an abrupt increase in sea surface temperature and major shifts in the catch of many fish species. This late 1970s regime shift in the Gulf of Alaska was followed by another shift in the late 1980s, not as pervasive as the 1977 shift, but which nevertheless did not return to the prior state. A thorough understanding of the extent and mechanisms leading to such regime shifts is challenged by data paucity in time and space. We investigate the ability of a suite of ocean biogeochemistry models of varying complexity to simulate regime shifts in the Gulf of Alaska by examining the presence of abrupt changes in time series of physical variables (sea surface temperature and mixed layer depth), nutrients and biological variables (chlorophyll, primary productivity and plankton biomass) using change-point analysis. Our study demonstrates that ocean biogeochemical models are capable of simulating the late 1970s shift, indicating an abrupt increase in sea surface temperature forcing followed by an abrupt decrease in nutrients and biological productivity. This predicted shift is consistent among all the models, although some of them exhibit an abrupt transition (i.e. a significant shift from one year to the next), whereas others simulate a smoother transition. Some models further suggest that the late 1980s shift was constrained by changes in mixed layer depth. Our study demonstrates that ocean biogeochemical can successfully simulate regime shifts in the Gulf of Alaska region, thereby providing better understanding of how changes in physical conditions are propagated from lower to upper trophic levels through bottom-up controls.

Defining scenarios of future vectors of change in marine life and associated economic sectors

Addressing the multitude of challenges in marine policy requires an integrated approach that considers the multitude of drivers, pressures, and interests, from several disciplinary angles. Scenarios are needed to harmonise the analyses of different components of the marine system, and to deal with the uncertainty and complexity of the societal and biogeophysical dynamics in the system. This study considers a set of socio-economic scenarios to (1) explore possible futures in relation to marine invasive species, outbreak forming species, and gradual changes in species distribution and productivity; and (2) harmonise the projection modelling performed within associated studies. The exercise demonstrates that developing interdisciplinary scenarios as developed in this study is particularly complicated due to (1) the wide variety in endogeneity or exogeneity of variables in the different analyses involved; (2) the dual role of policy decisions as variables in a scenario or decisions to be evaluated and compared to other decisions; and (3) the substantial difference in time scale between societal and physical drivers.

iMarNet: an ocean biogeochemistry model intercomparison project within a common physical ocean modelling framework

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

Quantitative pathways for Northeast Atlantic fisheries based on climate, ecological–economic and governance modelling scenarios

Here we present quantitative projections of potential futures for ecosystems in the North Atlantic basin generated from coupling a climate change-driven biophysical model (representing ecosystem and fish populations under climate change) and a scenario-driven ecological–economic model (representing fleets and industries under economic globalization). Four contrasting scenarios (Baseline, Fortress, Global Commons, Free Trade) were defined from the perspective of alternative regional management and governance of the oceanic basin, providing pathways for the future of ecosystems in the Northeast Atlantic basin by 2040. Results indicate that in the time frame considered: (1) the effects of governance and trade decisions are more significant in determining outcomes than the effects of climate change alone, (2) climate change is likely to result in a poleward latitudinal shift of species ranges and thus resources, with implications for exploitation patterns, (3) the level of fisheries regulation is the most important factor in determining the long term evolution of the fisheries system, (4) coupling climate change and governance impacts demonstrates the complex interaction between different components of this social–ecological system, (5) an important driver of change for the future of the North Atlantic and the European fishing fleets appears to be the interplay between wild fisheries and aquaculture development, and finally (6) scenarios demonstrate that the viability and profit of fisheries industries is highly volatile. This study highlights the need to explore basin-scale policy that combines medium to long-term environmental and socio-economic considerations, and the importance of defining alternative sustainable pathways.

Modelling the Stoichiometric Regulation of C-Rich Toxins in Marine Dinoflagellates

Toxin production in marine microalgae was previously shown to be tightly coupled with cellular stoichiometry. The highest values of cellular toxin are in fact mainly associated with a high carbon to nutrient cellular ratio. In particular, the cellular accumulation of C-rich toxins (i.e., with C:N > 6.6) can be stimulated by both N and P deficiency. Dinoflagellates are the main producers of C-rich toxins and may represent a serious threat for human health and the marine ecosystem. As such, the development of a numerical model able to predict how toxin production is stimulated by nutrient supply/deficiency is of primary utility for both scientific and management purposes. In this work we have developed a mechanistic model describing the stoichiometric regulation of C-rich toxins in marine dinoflagellates. To this purpose, a new formulation describing toxin production and fate was embedded in the European Regional Seas Ecosystem Model (ERSEM), here simplified to describe a monospecific batch culture. Toxin production was assumed to be composed by two distinct additive terms; the first is a constant fraction of algal production and is assumed to take place at any physiological conditions. The second term is assumed to be dependent on algal biomass and to be stimulated by internal nutrient deficiency. By using these assumptions, the model reproduced the concentrations and temporal evolution of toxins observed in cultures of Ostreopsis cf. ovata, a benthic/epiphytic dinoflagellate producing C-rich toxins named ovatoxins. The analysis of simulations and their comparison with experimental data provided a conceptual model linking toxin production and nutritional status in this species. The model was also qualitatively validated by using independent literature data, and the results indicate that our formulation can be also used to simulate toxin dynamics in other dinoflagellates. Our model represents an important step towards the simulation and prediction of marine algal toxicity.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs); six types of phytoplankton, three types of zooplankton, and heterotrophic bacteria. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing zooplankton, and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean High Nutrient Low Chlorophyll (HNLC) region during summer. When model simulations do not represent crustacean macrozooplankton grazing, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there was no iron deposition from dust. When model simulations included the developments of the zooplankton component, the simulation of phytoplankton biomass improved and the high chlorophyll summer bias in the Southern Ocean HNLC region largely disappeared. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community rather than iron limitation. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.

Global Carbon Budget 2015

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2, and land-cover change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005–2014), EFF was 9.0 ± 0.5 GtC yr−1, ELUC was 0.9 ± 0.5 GtC yr−1, GATM was 4.4 ± 0.1 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 3.0 ± 0.8 GtC yr−1. For the year 2014 alone, EFF grew to 9.8 ± 0.5 GtC yr−1, 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr−1 that took place during 2005–2014. Also, for 2014, ELUC was 1.1 ± 0.5 GtC yr−1, GATM was 3.9 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and SLAND was 4.1 ± 0.9 GtC yr−1. GATM was lower in 2014 compared to the past decade (2005–2014), reflecting a larger SLAND for that year. The global atmospheric CO2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in EFF will be near or slightly below zero, with a projection of −0.6 [range of −1.6 to +0.5] %, based on national emissions projections for China and the USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the global economy for the rest of the world. From this projection of EFF and assumed constant ELUC for 2015, cumulative emissions of CO2 will reach about 555 ± 55 GtC (2035 ± 205 GtCO2) for 1870–2015, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2015).

Impacts of Climate Modes on Air-Sea Heat Exchange in the Red Sea

The impacts of various climate modes on the Red Sea surface heat exchange are investigated using the MERRA reanalysis and the OAFlux satellite reanalysis datasets. Seasonality in the atmospheric forcing is also explored. Mode impacts peak during boreal winter [December–February (DJF)] with average anomalies of 12–18 W m−2 to be found in the northern Red Sea. The North Atlantic Oscillation (NAO), the east Atlantic–west Russia (EAWR) pattern, and the Indian monsoon index (IMI) exhibit the strongest influence on the air–sea heat exchange during the winter. In this season, the largest negative anomalies of about −30 W m−2 are associated with the EAWR pattern over the central part of the Red Sea. In other seasons, mode-related anomalies are considerably lower, especially during spring when the mode impacts are negligible. The mode impacts are strongest over the northern half of the Red Sea during winter and autumn. In summer, the southern half of the basin is strongly influenced by the multivariate ENSO index (MEI). The winter mode–related anomalies are determined mostly by the latent heat flux component, while in summer the shortwave flux is also important. The influence of the modes on the Red Sea is found to be generally weaker than on the neighboring Mediterranean basin.

Coastal upwelling and downwelling forcing of circulation in a semi-enclosed bay: Ria de Vigo

Semi-enclosed bays in upwelling regions are exposed to forcing related to winds, currents and buoyancy over the shelf. The influence of this external forcing is moderated by factors such as connectivity to the open ocean, shelter by surrounding topography, dimensions of the bay, and freshwater outflows. Such bays, preferred locations for ports, mariculture, marine industry, recreational activities and coastal settlement, present a range of characteristics, understanding of which is necessary to their rational management. Observations in such a semi-enclosed bay, the Ria de Vigo in Spain, are used to characterize the influence of upwelling and downwelling pulses on its circulation. In this location, near the northern limit of the Iberian upwelling system, upwelling events dominate during a short summer season and downwelling events the rest of the year. The ria response to the external forcing is central to nutrient supply and resultant plankton productivity that supports its high level of cultured mussel production. Intensive field studies in September 2006 and June 2007 captured a downwelling event and an upwelling event, respectively. Data from eight current profiler moorings and boat-based MiniBat/ADCP surveys provided an unprecedented quasi-synoptic view of the distribution of water masses and circulation patterns in any ria. In the outer ria, circulation was dominated by the introduction of wind-driven alongshore flow from the external continental shelf through the ria entrances and its interaction with the topography. In the middle ria, circulation was primarily related to the upwelling/downwelling cycle, with a cool, salty and dense lower layer penetrating to the inner ria during upwelling over the shelf. A warmer, lower salinity and less dense surface layer of coastal waters flowed inward during downwelling. Without external forcing, the inner ria responded primarily to tides and buoyancy changes related to land runoff. Under both upwelling and downwelling conditions, the flushing of the ria involved shelf responses to wind pulses. Their persistence for a few days was sufficient to allow waters from the continental shelf to penetrate the innermost ria. Longer term observations supported by numerical modeling are required to confirm the generality of such flushing events in the ria and determine their typical frequency, while comparative studies should explore how these scenarios fit into the range of conditions experienced in other semi-enclosed bays.

ERSEM 15.06: a generic model for marine biogeochemistry and the ecosystem dynamics of the lower trophic levels

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.