Abstract:Numerical models of ocean biogeochemistry are becoming the major tools used to detect and predict the impact of climate change on marine resources and to monitor ocean health. However, with the continuous improvement of model structure and spatial resolution, incorporation of these additional degrees of freedom into fidelity assessment has become increasingly challenging. Here, we propose a new method to provide information on the model predictive skill in a concise way. The method is based on the conjoint use of a k -means clustering technique, assessment metrics, and Biogeochemical-Argo (BGC-Argo) observations. The k -means algorithm and the assessment metrics reduce the number of model data points to be evaluated. The metrics evaluate either the model state accuracy or the skill of the model with respect to capturing emergent properties, such as the deep chlorophyll maximums and oxygen minimum zones. The use of BGC-Argo observations as the sole evaluation data set ensures the accuracy of the data, as it is a homogenous data set with strict sampling methodologies and data quality control procedures. The method is applied to the Global Ocean Biogeochemistry Analysis and Forecast system of the Copernicus Marine Service. The model performance is evaluated using the model efficiency statistical score, which compares the model–observation misfit with the variability in the observations and, thus, objectively quantifies whether the model outperforms the BGC-Argo climatology. We show that, overall, the model surpasses the BGC-Argo climatology in predicting pH, dissolved inorganic carbon, alkalinity, oxygen, nitrate, and phosphate in the mesopelagic and the mixed layers as well as silicate in the mesopelagic layer. However, there are still areas for improvement with respect to reducing the model–data misfit for certain variables such as silicate, pH, and the partial pressure of CO 2 in the mixed layer as well as chlorophyll- a -related, oxygen-minimum-zone-related, and particulate-organic-carbon-related metrics. The method proposed here can also aid in refining the design of the BGC-Argo network, in particular regarding the regions in which BGC-Argo observations should be enhanced to improve the model accuracy via the assimilation of BGC-Argo data or process-oriented assessment studies. We strongly recommend increasing the number of observations in the Arctic region while maintaining the existing high-density of observations in the Southern Oceans. The model error in these regions is only slightly less than the variability observed in BGC-Argo measurements. Our study illustrates how the synergic use of modeling and BGC-Argo data can both provide information about the performance of models and improve the design of observing systems.

A causality-based method for multi-model comparison: Application to relationships between atmospheric and marine biogeochemical variables

Responses of Summer Phytoplankton Biomass to Changes in Top-Down Forcing: Insights from Comparative Modelling

A comparison of two causal methods in the context of climate analyses

Structure of the Atlantic Meridional Overturning Circulation in Three Generations of Climate Models

On the ability of OMIP models to simulate the ocean mixed layer depth and its seasonal cycle in the Arctic Ocean

Using machine learning and Biogeochemical-Argo (BGC-Argo) floats to assess biogeochemical models and optimize observing system design

Impact Of Coherent Ocean Stratification On Amoc Reconstruction By Coupled Data Assimilation With A Biased Model

Influence of Two Atmospheric Transport Models on Inferring Sources and Sinks of Atmospheric Co2

North Atlantic Ocean-Originated Multicentennial Oscillation of the AMOC: A Coupled Model Study

Centennial‐Scale Variability of the Atlantic Meridional Overturning Circulation in CMIP6 Models Shaped by Arctic–North Atlantic Interactions and Sea Ice Biases

Effects of Model Resolution and Ocean Mixing on Forced Ice-Ocean Physical and Biogeochemical Simulations Using Global and Regional System Models

Diagnosing the causes of AMOC slowdown in a coupled model: a cautionary tale

A study of mesoscale air–sea interaction in the Southern Ocean with a regional coupled model

Inverse Modeling of Annual Atmospheric Co2 Sources and Sinks 2. Sensitivity Study

Centennial-scale variability of the Atlantic Meridional Circulation in CMIP6 models shaped by Arctic-North Atlantic interactions and sea ice biases

A study of capturing Atlantic meridional overturning circulation (AMOC) regime transition through observation-constrained model parameters

Contrasting Spatial Structures of Atlantic Multidecadal Oscillation Between Observations and Slab Ocean Model Simulations

Subpolar Gyre – AMOC – Atmosphere Interactions on Multidecadal Timescales in a Version of the Kiel Climate Model

The Seasonal Cycle of Atmospheric Co2: A Study Based on the Ncar Community Climate Model (Ccm2)

Ecosystem model optimization using in situ flux observations: benefit of Monte Carlo versus variational schemes and analyses of the year-to-year model performances

Large diversity in AMOC internal variability across NEMO-based climate models