Description: Global biogeochemical ocean models contain a variety of different biogeochemical components and often much simplified representations of complex dynamical interactions, which are described by many (≈10–≈100) parameters. The values of many of these parameters are empirically difficult to constrain, due to the fact that in the models they represent processes for a range of different groups of organisms at the same time, while even for single species parameter values are often difficult to determine in situ. Therefore, these models are subject to a high level of parametric uncertainty. This may be of consequence for their skill with respect to accurately describing the relevant features of the present ocean, as well as their sensitivity to possible environmental changes. We here present a framework for the calibration of global biogeochemical ocean models on short and long time scales. The framework combines an offline approach for transport of biogeochemical tracers with an Estimation of Distribution Algorithm (Covariance Matrix Adaption Evolution Strategy, CMAES). We explore the performance and capability of this framework by five different optimizations of six biogeochemical parameters of a global biogeochemical model. First, a twin experiment explores the feasibility of this approach. Four optimizations against a climatology of observations of annual mean dissolved nutrients and oxygen determine the extent, to which different setups of the optimization influence model's fit and parameter estimates. Because the misfit function applied focuses on the large-scale distribution of inorganic biogeochemical tracers, parameters that act on large spatial and temporal scales are determined earliest, and with the least spread. Parameters more closely tied to surface biology, which act on shorter time scales, are more difficult to determine. In particular the search for optimum zooplankton parameters can benefit from a sound knowledge of maximum and minimum parameter values, leading to a more efficient optimization. It is encouraging that, although the misfit function does not contain any direct information about biogeochemical turnover, the optimized models nevertheless provide a better fit to observed global biogeochemical fluxes.

Description: The Apparent Oxygen Utilisation (AOU) is a classical measure of the amount of oxygen respired in the ocean's interior. We show that AOU systematically overestimates True Oxygen Utilisation (TOU) in 6 coupled circulation-biogeochemical ocean models. This is due to atmosphere–ocean oxygen disequilibria in the subduction regions, consistent with previous work. We develop a simple, new, observationally-based approach which we call Evaluated Oxygen Utilisation (EOU). In this approach, we take into account the impact of the upper ocean oxygen disequilibria into the interior, considering that transport takes place predominantly along isopycnal surfaces. The EOU approximates the TOU with less than half of the bias of AOU in all 6 models despite large differences in the physical and biological components of the models. Applying the EOU approach to a global observational dataset yields an oxygen consumption rate 25% lower than that derived from AOU-based estimates, for a given ventilation rate.

Description: Highlights: • We examine the role of marine particle for regulating trace element distribution. • We review the state of the art for modelling the oceanic distribution of specific tracers: Thorium, Protactinium, Iron, and Aluminium. • We review the state of the art for modelling particle distribution in large scale ocean biogeochemical model. The distribution of trace elements in the ocean is governed by the combined effects of various processes, and by exchanges with external sources. Modelling these represents an opportunity to better understand and quantify the mechanisms that regulate the oceanic tracer cycles. Observations collected during the GEOTRACES program provide an opportunity to improve our knowledge regarding processes that should be considered in biogeochemical models to adequately represent the distributions of trace elements in the ocean. Here we present a synthesis about the state of the art for simulating selected trace elements in biogeochemical models: Protactinium, Thorium, Iron and Aluminium. In this contribution we pay particular attention on the role of particles in the cycling of these tracers and how they may provide additional constraints on the transfer of matter in the ocean.

Description: Although of substantial importance for marine tracer distributions and eventually global carbon, oxygen, and nitrogen fluxes, the interaction between sinking and remineralization of organic matter, benthic fluxes and burial is not always represented consistently in global biogeochemical models. We here aim to investigate the relationships between these processes with a suite of global biogeochemical models, each simulated over millennia, and compared against observed distributions of pelagic tracers and benthic and pelagic fluxes. We concentrate on the representation of sediment–water interactions in common numerical models, and investigate their potential impact on simulated global sediment–water fluxes and nutrient and oxygen distributions. We find that model configurations with benthic burial simulate global oxygen well over a wide range of possible sinking flux parameterizations, making the model more robust with regard to uncertainties about the remineralization length scale. On a global scale, burial mostly affects oxygen in the meso- to bathypelagic zone. While all model types show an almost identical fit to observed pelagic particle flux, and the same sensitivity to particle sinking speed, comparison to observational estimates of benthic fluxes reveals a more complex pattern, but definite interpretation is not straightforward because of heterogeneous data distribution and methodology. Still, evaluating model results against observed pelagic and benthic fluxes of organic matter can complement model assessments based on more traditional tracers such as nutrients or oxygen. Based on a combined metric of dissolved tracers and biogeochemical fluxes, we here identify two model descriptions of burial as suitable candidates for further experiments and eventual model refinements.

Description: We present a suite of experiments with a hierarchy of biogeochemical models of increasing complexity coupled to an offline global ocean circulation model based on the “transport matrix method”. Biogeochemical model structures range from simple nutrient models to more complex nutrient-phytoplankton–zooplankton-detritus-DOP models. The models’ skill is assessed by various misfit functions with respect to observed phosphate and oxygen distributions. While there is generally good agreement between the different metrics employed, an exception is a cost function based on the relative model-data misfit. We show that alterations in parameters and/or structure of the models – especially those that change particle export or remineralization profile – affect subsurface and mesopelagic phosphate and oxygen, particularly in the upwelling regions. Visual inspection of simulated biogeochemical tracer distributions as well as the evaluation of different cost functions suggest that increasing complexity of untuned, unoptimized models, simulated with parameters commonly used in large-scale model studies does not necessarily improve performance. Instead, variations in individual model parameters may be of equal, if not greater, importance.

Description: Amodel is presented that simulates the formation of marine aggregates from particles of different origin inside amodel of pelagic biological processes. Experiments are carried out with parameterizations appropriate for different types of aggregates, using different kinds of physical forcing, and compared to observations of dissolved inorganic nitrogen (DIN), particulate organic nitrogen (PON), marinesnow concentration, and sedimentation. The occurrence of large, macroscopically visible aggregates (marinesnow) can best be simulated with parameterizations that have been derived from in situ observations of marinesnow, but not with aparameterization sufficient for dense particles. The parameterization strongly determines the amount and timing of deep export, as well as the post-bloom development of the food web in the upper layers. Detritus in aggregates plays a role mainly during times when zooplankton are abundant, as e.g. in the western Arabian Sea during Southwest Monsoon. Then the large aggregates as fast sinking vehicles may remove detritus quickly from shallow and mid-water depth, preventing the accumulation of nutrients that are produced via detritus decomposition. In this region, detritus contributes strongly to deep sedimentation. The nitrogenbudget at this location with regard to the observations cannot be closed: depending on model type, either the model simulates too high sedimentation, or too high DIN. Possible causes for this mismatch include undercollection by sediment traps, inaccurate representation of physical processes in the model and the neglect of biological processes, such as production of dissolved organic matter or denitrification.

Description: This study presents results from 46 sensitivity experiments carried out with three structurally simple (2, 3, and 6 biogeochemical state variables, respectively) models of production, export and remineralization of organic phosphorus, coupled to a global ocean circulation model and integrated for 3000 years each. The models’ skill is assessed via different misfit functions with respect to the observed global distributions of phosphate and oxygen. Across the different models, the global root-mean square misfit with respect to observed phosphate and oxygen distributions is found to be particularly sensitive to changes in the remineralization length scale, and also to changes in simulated primary production. For this metric, changes in the production and decay of dissolved organic phosphorus as well as in zooplankton parameters are of lesser importance. For a misfit function accounting for the misfit of upper-ocean tracers, however, production parameters and organic phosphorus dynamics play a larger role. Regional misfit patterns are investigated as indicators of potential model deficiencies, such as missing iron limitation, or deficiencies in the sinking and remineralization length scales. In particular, the gradient between phosphate concentrations in the northern North Pacific and the northern North Atlantic is controlled predominantly by the biogeochemical model parameters related to particle flux. For the combined 46 sensitivity experiments performed here, the global misfit to observed oxygen and phosphate distributions shows no clear relation to either simulated global primary or export production for either misfit metric employed. However, a relatively tight relationship that is very similar for the different model of different structural complexity is found between the model-data misfit in oxygen and phosphate distributions to simulated meso- and bathypelagic particle flux. Best agreement with the observed tracer distributions is obtained for simulated particle fluxes that agree most closely with sediment trap data for a nominal depth of about 1000 m, or deeper.

Description: Ocean biogeochemistry has a large effect on marine biogeochemical tracer distribution, comparable to that of circulation. However, the exact form of many biogeochemical processes, and therefore their parameterisation in global models, is not well known. Using the "Transport Matrix Method" (Khatiwala et al., 2005) as an offline tool to quickly equilibrate global coupled ocean biogeochemical models, we carry out sensitivity analyses to assess the skill of different biogeochemical, physical, and numerical setups. The ultimate aim is to find the "best" (with respect to metrics tailored for specific research questions) biogeochemical model, that may eventually be integrated into e.g. high-resolution models,or models that simulate transient scenarios.

Description: Global biogeochemical ocean models are used to predict the future evolution of so-called oxygen minimum zones (OMZ), and the associated environmental and possible socio-economic impacts. Different models give different results and vary largely in their biogeochemical, physical and numerical setup. In order to assess the ability of the models to describe the present state as a necessary condition for skillful predictions into the future, they are usually compared against observed distributions of oxygen and other variables, such as thickness of oxygen minimum zones, nutrients, tracers for circulation and/or water mass age. We here examine different metrics for skill evaluation particularly of model representations of oxygen (and OMZs), for a wide range of global biogeochemical models. Among the metrics considered are Taylor plots, volume distributions of oxygen, volume of OMZ, preformed oxygen, and metrics that combine various diagnostic biogeochemical tracers. We finally investigate the impact these metrics may have for the “choice” of any best model, and discuss their applicability for different research or societal questions.

Description: In the aftermath of an earthquake and tsunami on 11 March 2011 radioactive 137Cs was discharged from a damaged nuclear power plant to the sea off Fukushima Dai-ichi, Japan. Here we explore its dilution and fate with a state-of-the-art global ocean general circulation model, which is eddy-resolving in the region of interest. We find apparent consistency between our simulated circulation, estimates of 137Cs discharged ranging from 0.94 p Bq (Japanese Government, 2011) to 3.5 ± 0.7 p Bq (Tsumune et al., 2012), and measurements by Japanese authorities and the power plant operator. In contrast, our simulations are apparently inconsistent with the high 27 ± 15 p Bq discharge estimate of Bailly du Bois et al. (2012). Expressed in terms of a diffusivity we diagnose, from our simulations, an initial dilution on the shelf of 60 to 100 m2 s−1. The cross-shelf diffusivity is at 500 ± 300 m2 s−1 significantly higher and variable in time as indicated by its uncertainty. Expressed as an effective residence time of surface water on the shelf, the latter estimate transfers to 43 ± 16 days. As regards the fate of 137Cs, our simulations suggest that activities up to 4 mBq l−1 prevail in the Kuroshio-Oyashio Interfrontal Zone one year after the accident. This allows for low but detectable 0.1 to 0.3 m Bq l−1 entering the North Pacific Intermediate Water before the 137Cs signal is flushed away. The latter estimates concern the direct release to the sea only.