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.

Biogenic aerosol precursors from phytoplankton production can affect cloud properties, especially in remote regions such as the Arctic Ocean. Reliable estimates on variability and trend of these precursors are required as extensive measurements in the Arctic are still scarce. We present a setup of the coupled ocean biogeochemical model FESOM2.1-REcoM3 where we integrated dissolved carboxylic acid containing polysaccharides (PCHO) and Transparent Exopolymer Particles (TEP) to describe these precursors in the upper ocean. We define PCHO as one part of the excreted organic carbon, which can then aggregate to form larger particles, TEP. Compared to observations, the simulation provides a valid TEP estimate with mean concentrations of 200-400 µg C L-1 on the continental shelves and 10-50 µg C L-1 in the central basins (0-30 m depth range). Further, the simulation for 1990-2019 reveals a significant positive trend of TEP of 0.5-3 µg C L-1 yr-1 during July-September in the Amerasian Basin (+3.5% yr-1), the Canadian Archipelago (+1.2% yr-1) and the Kara Sea (+0.8% yr-1), in contrast to the eastern Fram Strait (-0.4% yr-1), the Barents Sea (-0.3% yr-1), and parts of the Eurasian Basin with a significant decrease of -0.5-2 µg C L-1 yr-1. Our study provides for the first time an integration of TEP formation, aggregation and remineralization processes into a global ocean biogeochemical model. This simulation assembles valuable data on biogenic aerosol precursors, and as such, fills a gap on which Earth System Models can greatly benefit to improve the understanding of aerosol feedbacks within the Arctic climate.

We assess the Southern Ocean CO2 uptake (1985-2018) using data sets gathered in the REgional Carbon Cycle Assessment and Processes Project phase 2 (RECCAP2). The Southern Ocean acted as a sink for CO2 with close agreement between simulation results from global ocean biogeochemistry models (GOBMs, 0.75±0.28 PgCyr-1) and pCO2-observation-based products (0.73±0.07 PgCyr-1). This sink is only half that reported by RECCAP1. The present-day net uptake is to first order a response to rising atmospheric CO2, driving large amounts of anthropogenic CO2 (Cant) into the ocean, thereby overcompensating the loss of natural CO2 to the atmosphere. An apparent knowledge gap is the increase of the sink since 2000, with pCO2-products suggesting a growth that is more than twice as strong and uncertain as that of GOBMs (0.26±0.06 and 0.11±0.03 PgCyr-1 decade-1 respectively). This is despite nearly identical pCO2 trends in GOBMs and pCO2-products when both products are compared only at the locations where pCO2 was measured. Seasonal analyses revealed agreement in driving processes in winter with uncertainty in the magnitude of outgassing, whereas discrepancies are more fundamental in summer, when GOBMs exhibit difficulties in simulating the effects of the non-thermal processes of biology and mixing/circulation. Ocean interior accumulation of Cant points to an underestimate of Cant uptake and storage in GOBMs. Future work needs to link surface fluxes and interior ocean transport, build long overdue systematic observation networks and push towards better process understanding of drivers of the carbon cycle.

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.

As a contribution to the Regional Carbon Cycle Assessment and Processes phase 2 (RECCAP2) project, we present synthesized estimates of Arctic Ocean sea-air CO2 fluxes and their uncertainties from 8 surface ocean pCO2-observation products, 18 ocean biogeochemical hindcast and data assimilation models and 6 atmospheric inversions. For the period of 1985−2018, the Arctic Ocean was a net sink of CO2 of 116 ± 4 TgC yr−1 in the pCO2 products and 92 ± 30 TgC yr−1 in the models. The CO2 uptake peaks in late summer and early autumn, and is low in winter when sea ice inhibits sea-air fluxes. The long-term mean CO2 uptake in the Arctic Ocean is primarily caused by steady-state fluxes of natural carbon (70 ± 15 %), and enhanced by the atmospheric CO2 increase (19 ± 5 %) and climate change (11 ± 18 %). The annual mean CO2 uptake increased from 1985 to 2018 at a rate of 31 ±13 TgC yr−1dec-1 in the pCO2 products and 10 ± 4 TgC yr−1dec-1 in the models. Moreover, 77 ± 38 % of the trend in the net CO2 uptake over time is caused by climate change, primarily due to rapid sea ice loss in recent years. Both, the mean CO2 uptake and the trend, is substantially weaker in the atmospheric inversions. Uncertainties across all estimates are large, in the pCO2 products because of scarcity of observations and in the models because of missing processes.

The coastal ocean contributes to regulating atmospheric greenhouse gas concentrations by taking up carbon dioxide (CO2) and releasing nitrous oxide (N2O) and methane (CH4). Major advances have improved our understanding of the coastal air-sea exchanges of these three gasses since the first phase of the Regional Carbon Cycle Assessment and Processes (RECCAP in 2013), but a comprehensive view that integrates the three gasses at the global scale is still lacking. In this second phase (RECCAP2), we quantify global coastal ocean fluxes of CO2, N2O and CH4 using an ensemble of global gap-filled observation-based products and ocean biogeochemical models. The global coastal ocean is a net sink of CO2 in both observational products and models, but the magnitude of the median net global coastal uptake is ~60% larger in models (-0.72 vs. -0.44 PgC/yr, 1998-2018, coastal ocean area of 77 million km2). We attribute most of this model-product difference to the seasonality in sea surface CO2 partial pressure at mid- and high-latitudes, where models simulate stronger winter CO2 uptake. The global coastal ocean is a major source of N2O (+0.70 PgCO2-e /yr in observational product and +0.54 PgCO2-e /yr in model median) and of CH4 (+0.21 PgCO2-e /yr in observational product), which offsets a substantial proportion of the net radiative effect of coastal \co uptake (35-58% in CO2-equivalents). Data products and models need improvement to better resolve the spatio-temporal variability and long term trends in CO2, N2O and CH4 in the global coastal ocean.

Fossil fuel combustion, land use change and other human activities have increased the atmospheric carbon dioxide (CO2) abundance by about 50% since the beginning of the industrial age. The atmospheric CO2 growth rates would have been much larger if natural sinks in the land biosphere and ocean had not removed over half of this anthropogenic CO2. As these CO2 emissions grew, uptake by the ocean increased in response to increases in atmospheric CO2 partial pressure (pCO2). On land, gross primary production (GPP) also increased, but the dynamics of other key aspects of the land carbon cycle varied regionally. Over the past three decades, CO2 uptake by intact tropical humid forests declined, but these changes are offset by increased uptake across mid- and high-latitudes. While there have been substantial improvements in our ability to study the carbon cycle, measurement and modeling gaps still limit our understanding of the processes driving its evolution. Continued ship-based observations combined with expanded deployments of autonomous platforms are needed to quantify ocean-atmosphere fluxes and interior ocean carbon storage on policy-relevant spatial and temporal scales. There is also an urgent need for more comprehensive measurements of stocks, fluxes and atmospheric CO2 in humid tropical forests and across the Arctic and boreal regions, which are experiencing rapid change. Here, we review our understanding of the atmosphere, ocean, and land carbon cycles and their interactions, identify emerging measurement and modeling capabilities and gaps and the need for a sustainable, operational framework to ensure a scientific basis for carbon management.

The ocean coastal-shelf-slope ecosystem west of the Antarctic Peninsula (WAP) is a biologically productive region that could potentially act as a large sink of atmospheric carbon dioxide. The duration of the sea-ice season in the WAP shows large interannual variability. However, quantifying the mechanisms by which sea ice impacts biological productivity and surface dissolved inorganic carbon (DIC) remains a challenge due to the lack of data early in the phytoplankton growth season. In this study, we implemented a circulation, sea-ice and biogeochemistry model (MITgcm-REcoM2) to study the effect of sea ice on phytoplankton blooms and surface DIC. Results were compared with satellite sea-ice and ocean color, and research ship surveys from the Palmer Long Term Ecological Research (LTER) program. The simulations suggest that the annual sea-ice cycle has an important role in the seasonal DIC drawdown. In years of early sea-ice retreat there is a longer growth season leading to larger seasonally integrated net primary production (NPP). Part of the biological uptake of DIC by phytoplankton, however, is counteracted by increased oceanic uptake of atmospheric CO2. Despite lower seasonal NPP, years of late sea-ice retreat show larger DIC drawdown, attributed to lower air-sea CO2 fluxes and increased dilution by sea-ice melt. The role of dissolved iron and iron limitation on WAP phytoplankton also remains a challenge due to the lack of data. The model results suggest sediments and glacial meltwater are the main sources in the coastal and shelf regions, with sediments being more influential in the northern coast.