Abstract. 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 to  ≈ 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 timescales. The framework combines an offline approach for transport of biogeochemical tracers with an estimation of distribution algorithm (Covariance Matrix Adaption Evolution Strategy, CMA-ES). We explore the performance and capability of this framework by five different optimizations of six biogeochemical parameters of a global biogeochemical model, simulated over 3000 years. 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 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 timescales, 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.

Abstract. Global biogeochemical ocean models are often tuned to match the observed distributions and fluxes of inorganic and organic quantities. This tuning is typically carried out “by hand”. However, this rather subjective approach might not yield the best fit to observations, is closely linked to the circulation employed and is thus influenced by its specific features and even its faults. We here investigate the effect of model tuning, via objective optimisation, of one biogeochemical model of intermediate complexity when simulated in five different offline circulations. For each circulation, three of six model parameters have been adjusted to characteristic features of the respective circulation. The values of these three parameters – namely, the oxygen utilisation of remineralisation, the particle flux parameter and potential nitrogen fixation rate – correlate significantly with deep mixing and ideal age of North Atlantic Deep Water (NADW) and the outcrop area of Antarctic Intermediate Waters (AAIW) and Subantarctic Mode Water (SAMW) in the Southern Ocean. The clear relationship between these parameters and circulation characteristics, which can be easily diagnosed from global models, can provide guidance when tuning global biogeochemistry within any new circulation model. The results from 20 global cross-validation experiments show that parameter sets optimised for a specific circulation can be transferred between similar circulations without losing too much of the model's fit to observed quantities. When compared to model intercomparisons of subjectively tuned, global coupled biogeochemistry–circulation models, each with different circulation and/or biogeochemistry, our results show a much lower range of oxygen inventory, oxygen minimum zone (OMZ) volume and global biogeochemical fluxes. Export production depends to a large extent on the circulation applied, while deep particle flux is mostly determined by the particle flux parameter. Oxygen inventory, OMZ volume, primary production and fixed-nitrogen turnover depend more or less equally on both factors, with OMZ volume showing the highest sensitivity, and residual variability. These results show a beneficial effect of optimisation, even when a biogeochemical model is first optimised in a relatively coarse circulation and then transferred to a different finer-resolution circulation model.

No abstract

Global biogeochemical ocean models rely on many parameters, which govern the interaction between individual components, and their response to the physical environment. They are often assessed/calibrated against quasi‐synoptic data sets of dissolved inorganic tracers. However, a good fit to one observation might not necessarily imply a good match to another. We investigate whether two different metrics—the root‐mean‐square error to nutrients and oxygen and a metric measuring the overlap between simulated and observed oxygen minimum zones (OMZs)—help to constrain a global biogeochemical model in different aspects of performance. Three global model optimizations are carried out. Two single‐objective optimizations target the root‐mean‐square metric and a sum of both metrics, respectively. We then present and explore multiobjective optimization, which results in a set of compromise solutions. Our results suggest that optimal parameters for denitrification and nitrogen fixation differ when applying different metrics. Optimization against observed OMZs leads to parameters that enhance fixed nitrogen cycling; this causes too low nitrate concentrations and a too high global pelagic denitrification rate. Optimization against nutrient and oxygen concentrations leads to different parameters and a lower global fixed nitrogen turnover; this results in a worse fit to OMZs. Multiobjective optimization resolves this antagonistic effect and provides an ensemble of parameter sets, which help to address different research questions. We finally discuss how systematic model calibration can help to improve models used for projecting climate change and its effect on fisheries and climate gas emissions.