This study characterized ocean biological carbon pump metrics in the second iteration of the REgional Carbon Cycle Assessment and Processes (RECCAP2) project, a coordinated, international effort to constrain contemporary ocean carbon air-sea fluxes and interior carbon storage trends using a combination of observation-based estimates, inverse models, and global ocean biogeochemical models. The analysis here focused on comparisons of global and biome-scale regional patterns in particulate organic carbon production and sinking flux from the RECCAP2 model ensemble against observational products derived from satellite remote sensing, sediment traps, and geochemical methods. There was generally encouraging model-data agreement in large-scale spatial patterns, though with substantial spread across the model ensemble and observational products. The global-integrated, model ensemble-mean export production, taken as the sinking particulate organic carbon flux at 100 m (6.41 ± 1.52 Pg C yr–1), and export ratio defined as sinking flux divided by net primary production (0.154 ± 0.026) both fell at the lower end of observational estimates. Comparison with observational constraints also suggested that the model ensemble may have underestimated regional biological CO2 drawdown and air-sea CO2 flux in high productivity regions. Reasonable model-data agreement was found for global-integrated, ensemble-mean sinking particulate organic carbon flux into the deep ocean at 1000 m (0.95 ± 0.64 Pg C yr–1) and the transfer efficiency defined as flux at 1000m divided by flux at 100m (0.121 ± 0.035), with both variables exhibiting considerable regional variability. Future modeling studies are needed to improve system-level simulation of interaction between model ocean physics and biogeochemical response.

The ocean is a major carbon sink and takes up 25-30% of the anthropogenically emitted CO2. A state-of-the-art method to quantify this sink are global ocean biogeochemistry models (GOBMs) but their simulated CO2 uptake differs between models and is systematically lower than estimates based on statistical methods using surface ocean pCO2 and interior ocean measurements. Here, we provide an in-depth evaluation of ocean carbon sink estimates from 1980 to 2018 from a GOBM ensemble. As sources of inter-model differences and ensemble-mean biases our study identifies the (i) model set-up, such as the length of the spin-up, the starting date of the simulation, and carbon fluxes from rivers and into sediments, (ii) the ocean circulation, such as Atlantic Meridional Overturning Circulation and Southern Ocean mode and intermediate water formation, and (iii) the oceanic buffer capacity. Our analysis suggests that the late starting date and biases in the ocean circulation cause a too low anthropogenic CO2 uptake across the GOBM ensemble. Surface ocean biogeochemistry biases might also cause simulated anthropogenic fluxes to be too low but the current set-up prevents a robust assessment. For simulations of the ocean carbon sink, we recommend in the short-term to (1) start simulations in 1765, when atmospheric CO2 started to increase, (2) conduct a sufficiently long spin-up such that the GOBMs reach steady-state, and (3) provide key metrics for circulation, biogeochemistry, and the land-ocean interface. In the long-term, we recommend improving the representation of these metrics in the GOBMs.