Description: 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.

Description: In geoscience and other fields, researchers use models as a simplified representation of reality. The models include processes that often rely on uncertain parameters that reduce model performance in reflecting real-world processes. The problem is commonly addressed by adapting parameter values to reach a good match between model simulations and corresponding observations. Different optimization tools have been successfully applied to address this task of model calibration. However, seeking one best value for every single model parameter might not always be optimal. For example, if model equations integrate over multiple real-world processes which cannot be fully resolved, it might be preferable to consider associated model parameters as random parameters. In this paper, a random parameter is drawn from a wide probability distribution for every singe model simulation. We developed an optimization approach that allows us to declare certain parameters random while optimizing those that are assumed to take fixed values. We designed a corresponding variant of the well known Covariance Matrix Adaption Evolution Strategy (CMA-ES). The new algorithm was applied to a global biogeochemical circulation model to quantify the impact of zooplankton mortality on the underlying biogeochemistry. Compared to the deterministic CMA-ES, our new method converges to a solution that better suits the credible range of the corresponding random parameter with less computational effort.