The carbonate system of the NW European shelf: sensitivity and model validation

The ocean plays an important role in regulating the climate, acting as a sink for carbon dioxide, perturbing the carbonate system and resulting in a slow decrease of seawater pH. Understanding the dynamics of the carbonate system in shelf sea regions is necessary to evaluate the impact of Ocean Acidification (OA) in these societally important ecosystems. Complex hydrodynamic and ecosystem coupled models provide a method of capturing the significant heterogeneity of these areas. However rigorous validation is essential to properly assess the reliability of such models. The coupled model POLCOMS–ERSEM has been implemented in the North Western European shelf with a new parameterization for alkalinity explicitly accounting for riverine inputs and the influence of biological processes. The model has been validated in a like with like comparison with North Sea data from the CANOBA dataset. The model shows good to reasonable agreement for the principal variables, physical (temperature and salinity), biogeochemical (nutrients) and carbonate system (dissolved inorganic carbon and total alkalinity), but simulation of the derived variables, pH and pCO2, are not yet fully satisfactory. This high uncertainty is attributed mostly to riverine forcing and primary production. This study suggests that the model is a useful tool to provide information on Ocean Acidification scenarios, but uncertainty on pH and pCO2 needs to be reduced, particularly when impacts of OA on ecosystem functions are included in the model systems.

Exploring the Red Sea seasonal ecosystem functioning using a three-dimensional biophysical model

The Red Sea exhibits complex hydrodynamic and biogeochemical dynamics, which vary both in time and space. These dynamics have been explored through the development and application of a 3-D ecosystem model. The simulation system comprises two off-line coupled submodels: the MIT General Circulation Model (MITgcm) and the European Regional Seas Ecosystem Model (ERSEM), both adapted for the Red Sea. The results from an annual simulation under climatological forcing are presented. Simulation results are in good agreement with satellite and in situ data illustrating the role of the physical processes in determining the evolution and variability of the Red Sea ecosystem. The model was able to reproduce the main features of the Red Sea ecosystem functioning, including the exchange with the Gulf of Aden, which is a major driving mechanism for the whole Red Sea ecosystem and the winter overturning taking place in the north. Some model limitations, mainly related to the dynamics of the extended reef system located in the southern part of the Red Sea, which is not currently represented in the model, still need to be addressed.

The biogeochemical structuring role of horizontal stirring: Lagrangian perspectives on iron delivery downstream of the Kerguelen Plateau

Field campaigns are instrumental in providing ground truth for understanding and modeling global ocean biogeochemical budgets. A survey however can only inspect a fraction of the global oceans, typically a region hundreds of kilometers wide for a temporal window of the order of (at most) several weeks. This spatiotemporal domain is also the one in which the mesoscale activity induces through horizontal stirring a strong variability in the biogeochemical tracers, with ephemeral, local contrasts which can easily mask the regional and seasonal gradients. Therefore, whenever local in situ measures are used to infer larger-scale budgets, one faces the challenge of identifying the mesoscale structuring effect, if not simply to filter it out. In the case of the KEOPS2 investigation of biogeochemical responses to natural iron fertilization, this problem was tackled by designing an adaptive sampling strategy based on regionally optimized multisatellite products analyzed in real time by specifically designed Lagrangian diagnostics. This strategy identified the different mesoscale and stirring structures present in the region and tracked the dynamical frontiers among them. It also enabled back trajectories for the ship-sampled stations to be estimated, providing important insights into the timing and pathways of iron supply, which were explored further using a model based on first-order iron removal. This context was essential for the interpretation of the field results. The mesoscale circulation-based strategy was also validated post-cruise by comparing the Lagrangian maps derived from satellites with the patterns of more than one hundred drifters, including some adaptively released during KEOPS2 and a subsequent research voyage. The KEOPS2 strategy was adapted to the specific biogeochemical characteristics of the region, but its principles are general and will be useful for future in situ biogeochemical surveys.

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.

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.