Description: The Southern Ocean is the most important area of anthropogenic carbon (Cant) uptake in the world ocean, only rivalled in importance by the North Atlantic Ocean. Significant variability on decadal time-scales in the uptake of Cant in the Southern Ocean has been observed and modelled, likely with consequences for the interior ocean storage of Cant in the region, and implications for the global carbon budget. Here we use eight cruises between 1973 and 2012 to assess decadal variability in Cant storage rates in the southeast Atlantic sector of the Southern Ocean. For this we employed the extended multiple linear regression (eMLR) method. We relate variability in DIC (dissolved inorganic carbon) storage, which is assumed to equal anthropogenic carbon storage, to changes in ventilation as observed from repeat measurements of transient tracers. Within the Antarctic Intermediate Water (AAIW) layer, which is the dominant transport conduit for Cant into the interior ocean, moderate Cant storage rates were found without any clear temporal trend. In Subantarctic Mode Water (SAMW), a less dense water mass found north of the Subantarctic Front and above AAIW, high storage rates of Cant were observed up to about 2005 but lower rates in more recent times. The transient tracer data suggest a significant speed-up of ventilation in the summer warmed upper part of AAIW between 1998 and 2012, which is consistent with the high storage rate of Cant. A shift of more northern Cant storage to more southern storage in near surface waters was detected in the early 2000s. Beneath the AAIW the eMLR method as applied here did not detect significant storage of Cant. However, the presence of the transient tracer CFC-12 all through the water column suggests that some Cant should be present, but at concentrations not reliably quantifiable. The observed temporal variability in the interior ocean seems at a first glance to be out of phase with observed surface ocean Cant fluxes, but this can be explained by the time delay for the surface ocean signal to manifest itself in the interior of the ocean.

Description: We analyzed 214 new core-top samples for their CaCO3 content from shelves all around Antarctica in order to understand their distribution and contribution to the marine carbon cycle. The distribution of sedimentary CaCO3 on the Antarctic shelves is connected to environmental parameters where we considered water depth, width of the shelf, sea-ice coverage and primary production. While CaCO3 contents of surface sediments are usually low, high (〉 15%) CaCO3 contents occur at shallow water depths (150–200 m) on the narrow shelves of the eastern Weddell Sea and at a depth range of 600–900 m on the broader and deeper shelves of the Amundsen, Bellingshausen and western Weddell Seas. Regions with high primary production, such as the Ross Sea and the western Antarctic Peninsula region, have generally low CaCO3 contents in the surface sediments. The predominant mineral phase of CaCO3 on the Antarctic shelves is low-magnesium calcite. With respect to ocean acidification, our findings suggest that dissolution of carbonates in Antarctic shelf sediments may be an important negative feedback only after the onset of calcite undersaturation on the Antarctic shelves. Macrozoobenthic CaCO3 standing stocks do not increase the CaCO3 budget significantly as they are two orders of magnitude lower than the budget of the sediments. This first circumpolar compilation of Antarctic shelf carbonate data does not claim to be complete. Future studies are encouraged and needed to fill data gaps especially in the under-sampled southwest Pacific and Indian Ocean sectors of the Southern Ocean.

Description: The influence of eddy structures on the seasonal depletion of dissolved inorganic carbon (DIC) and carbon dioxide (CO2) disequilibrium was investigated during a trans-Atlantic crossing of the Antarctic Circumpolar Current (ACC) in austral summer 2012. The Georgia Basin, downstream of the island of South Georgia (54-55°S, 36-38°W) is a highly dynamic region due to the mesoscale activity associated with the flow of the Subantarctic Front (SAF) and Polar Front (PF). Satellite sea-surface height and chlorophyll-a anomalies revealed a cyclonic cold core that dominated the northern Georgia Basin that was formed from a large meander of the PF. Warmer waters influenced by the SAF formed a smaller anticyclonic structure to the east of the basin. Both the cold core and warm core eddy structures were hotspots of carbon uptake relative to the rest of the ACC section during austral summer. This was most amplified in the cold core where greatest CO2 undersaturation (−78 μatm) and substantial surface ocean DIC deficit (5.1 mol m−2) occurred. In the presence of high wind speeds, the cold core eddy acted as a strong sink for atmospheric CO2 of 25.5 mmol m−2 day−1. Waters of the warm core displayed characteristics of the Polar Frontal Zone (PFZ), with warmer upper ocean waters and enhanced CO2 undersaturation (−59 μatm) and depletion of DIC (4.9mol m−2). A proposed mechanism for the enhanced carbon uptake across both eddy structures is based on the Ekman eddy pumping theory: (i) the cold core is seeded with productive (high chlorophyll-a) waters from the Antarctic Zone and sustained biological productivity through upwelled nutrient supply that counteracts DIC inputs from deep waters; (ii) horizontal entrainment of low-DIC surface waters (biological uptake) from the PFZ downwell within the warm core and cause relative DIC-depletion in the upper water column. The observations suggest that the formation and northward propagation of cold core eddies in the region of the PF could project low-DIC waters towards the site of Antarctic Intermediate Water formation and enhance CO2 drawdown into the deep ocean.

Description: Harde (2017) proposes an alternative accounting scheme for the modern carbon cycle and concludes that only 4.3% of today's atmospheric CO2 is a result of anthropogenic emissions. As we will show, this alternative scheme is too simple, is based on invalid assumptions, and does not address many of the key processes involved in the global carbon cycle that are important on the timescale of interest. Harde (2017) therefore reaches an incorrect conclusion about the role of anthropogenic CO2 emissions. Harde (2017) tries to explain changes in atmospheric CO2 concentration with a single equation, while the most simple model of the carbon cycle must at minimum contain equations of at least two reservoirs (the atmosphere and the surface ocean), which are solved simultaneously. A single equation is fundamentally at odds with basic theory and observations. In the following we will (i) clarify the difference between CO2 atmospheric residence time and adjustment time, (ii) present recently published information about anthropogenic carbon, (iii) present details about the processes that are missing in Harde (2017), (iv) briefly discuss shortcoming in Harde's generalization to paleo timescales, (v) and comment on deficiencies in some of the literature cited in Harde (2017).

Description: Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1% relative to 2020, with fossil emissions at 10.1±0.5GtCyr-1 (9.9±0.5GtCyr-1 when the cement carbonation sink is included), and ELUC was 1.1±0.7GtCyr-1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 10.9±0.8GtCyr-1 (40.0±2.9GtCO2). Also, for 2021, GATM was 5.2±0.2GtCyr-1 (2.5±0.1ppmyr-1), SOCEAN was 2.9 ±0.4GtCyr-1, and SLAND was 3.5±0.9GtCyr-1, with a BIM of -0.6GtCyr-1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO2 concentration averaged over 2021 reached 414.71±0.1ppm. Preliminary data for 2022 suggest an increase in EFOS relative to 2021 of +1.0% (0.1% to 1.9%) globally and atmospheric CO2 concentration reaching 417.2ppm, more than 50% above pre-industrial levels (around 278ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959-2021, but discrepancies of up to 1GtCyr-1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at 10.18160/GCP-2022 (Friedlingstein et al., 2022b).

Description: The cycling of carbon in the oceans is affected by feedbacks driven by changes in climate and atmospheric CO2. Understanding these feedbacks is therefore an important prerequisite for projecting future climate. Marine biogeochemistry models are a useful tool but, as with any model, are a simplification and need to be continually improved. In this study, we coupled the Finite-volumE Sea ice–Ocean Model (FESOM2.1) to the Regulated Ecosystem Model version 3 (REcoM3). FESOM2.1 is an update of the Finite-Element Sea ice–Ocean Model (FESOM1.4) and operates on unstructured meshes. Unlike standard structured-mesh ocean models, the mesh flexibility allows for a realistic representation of small-scale dynamics in key regions at an affordable computational cost. Compared to the previous coupled model version of FESOM1.4–REcoM2, the model FESOM2.1–REcoM3 utilizes a new dynamical core, based on a finite-volume discretization instead of finite elements, and retains central parts of the biogeochemistry model. As a new feature, carbonate chemistry, including water vapour correction, is computed by mocsy 2.0. Moreover, REcoM3 has an extended food web that includes macrozooplankton and fast-sinking detritus. Dissolved oxygen is also added as a new tracer. In this study, we assess the ocean and biogeochemical state simulated with FESOM2.1–REcoM3 in a global set-up at relatively low spatial resolution forced with JRA55-do (Tsujino et al., 2018) atmospheric reanalysis. The focus is on the recent period (1958–2021) to assess how well the model can be used for present-day and future climate change scenarios on decadal to centennial timescales. A bias in the global ocean–atmosphere preindustrial CO2 flux present in the previous model version (FESOM1.4–REcoM2) could be significantly reduced. In addition, the computational efficiency is 2–3 times higher than that of FESOM1.4–REcoM2. Overall, it is found that FESOM2.1–REcoM3 is a skilful tool for ocean biogeochemical modelling applications.

Description: The partitioning of CO2 between atmosphere and ocean depends to a large degree not only on the amount of dissolved inorganic carbon (DIC) but also on alkalinity in the surface ocean. That is also why ocean alkalinity enhancement (OAE) is discussed as one potential approach in the context of negative emission technologies. Although alkalinity is thus an important variable of the marine carbonate system, little knowledge exists on how its representation in models compares with measurements. We evaluated the large-scale alkalinity distribution in 14 CMIP6 Earth system models (ESMs) against the observational data set GLODAPv2 and show that most models, as well as the multi-model mean, underestimate alkalinity at the surface and in the upper ocean and overestimate it in the deeper ocean. The decomposition of the global mean alkalinity biases into contributions from (i) physical processes (preformed alkalinity), which include the physical redistribution of biased alkalinity originating from the soft tissue and carbonates pumps; (ii) remineralization; and (iii) carbonate formation and dissolution showed that the bias stemming from the physical redistribution of alkalinity is dominant. However, below the upper few hundred meters the bias from carbonate dissolution can gain similar importance to physical biases, while the contribution from remineralization processes is negligible. This highlights the critical need for better understanding and quantification of processes driving calcium carbonate dissolution in microenvironments above the saturation horizons and implementation of these processes into biogeochemical models. For the application of the models to assess the potential of OAE to increase ocean carbon uptake, a back-of-the-envelope calculation was conducted with each model's global mean surface alkalinity, DIC, and partial pressure of CO2 in seawater (pCO2) as input parameters. We evaluate the following two metrics: (1) the initial pCO2 reduction at the surface ocean after alkalinity addition and (2) the uptake efficiency (ηCO2) after air–sea equilibration is reached. The relative biases of alkalinity versus DIC at the surface affect the Revelle factor and therefore the initial pCO2 reduction after alkalinity addition. The global mean surface alkalinity bias relative to GLODAPv2 in the different models ranges from −85 mmol m−3 (−3.6 %) to +50 mmol m−3 (+2.1 %) (mean: −25 mmol m−3 or −1.1 %). For DIC the relative bias ranges from −55 mmol m−3 (−2.6 %) to 53 mmol m−3 (+2.5 %) (mean: −13 mmol m−3 or −0.6 %). All but two of the CMIP6 models evaluated here overestimate the Revelle factor at the surface by up to 3.4 % and thus overestimate the initial pCO2 reduction after alkalinity addition by up to 13 %. The uptake efficiency, ηCO2, then takes into account that a higher Revelle factor and a higher initial pCO2 reduction after alkalinity addition and equilibration mostly compensate for each other, meaning that resulting DIC differences in the models are small (−0.1 % to 1.1 %). The overestimation of the initial pCO2 reduction has to be taken into account when reporting on efficiencies of ocean alkalinity enhancement experiments using CMIP6 models, especially as long as the CO2 equilibrium is not reached.

Description: 〈jats:p〉Abstract. In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14CDIC. 〈/jats:p〉