Search

​​​​​​​Ocean CO2 uptake and acidification in response to human activities are driven primarily by the rise in atmospheric CO2 but are also modulated by climate change. Existing work suggests that this “climate effect” influences the uptake and storage of anthropogenic carbon and acidification via the global increase in ocean temperature, although some regional responses have been attributed to changes in circulation or biological activity. Here, we investigate spatial patterns in the climate effect on surface ocean acidification (and the closely related carbonate chemistry) in an Earth system model under a rapid CO2-increase scenario and identify a different driving process. We show that the amplification of the hydrological cycle, a robustly simulated feature of climate change, is largely responsible for the spatial patterns in this climate effect at the sea surface. This “hydrological effect” can be understood as a subset of the total climate effect, which includes warming, hydrological cycle amplification, circulation, and biological changes. We demonstrate that it acts through two primary mechanisms: (i) directly diluting or concentrating dissolved ions by adding or removing freshwater and (ii) altering the sea surface temperature, which influences the solubility of dissolved inorganic carbon (DIC) and acidity of seawater. The hydrological effect opposes acidification in salinifying regions, most notably the subtropical Atlantic, and enhances acidification in freshening regions such as the western Pacific. Its single strongest effect is to dilute the negative ions that buffer the dissolution of CO2, quantified as alkalinity. The local changes in alkalinity, DIC, and pH linked to the pattern of hydrological cycle amplification are as strong as the (largely uniform) changes due to warming, explaining the weak increase in pH and DIC seen in the climate effect in the subtropical Atlantic Ocean.

Nitrous oxide (N2O) is a long-lived potent greenhouse gas and stratospheric ozone-depleting substance that has been accumulating in the atmosphere since the preindustrial period. The mole fraction of atmospheric N2O has increased by nearly 25 % from 270 ppb (parts per billion) in 1750 to 336 ppb in 2022, with the fastest annual growth rate since 1980 of more than 1.3 ppb yr−1 in both 2020 and 2021. According to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR6), the relative contribution of N2O to the total enhanced effective radiative forcing of greenhouse gases was 6.4 % for 1750–2022. As a core component of our global greenhouse gas assessments coordinated by the Global Carbon Project (GCP), our global N2O budget incorporates both natural and anthropogenic sources and sinks and accounts for the interactions between nitrogen additions and the biogeochemical processes that control N2O emissions. We use bottom-up (BU: inventory, statistical extrapolation of flux measurements, and process-based land and ocean modeling) and top-down (TD: atmospheric measurement-based inversion) approaches. We provide a comprehensive quantification of global N2O sources and sinks in 21 natural and anthropogenic categories in 18 regions between 1980 and 2020. We estimate that total annual anthropogenic N2O emissions have increased 40 % (or 1.9 Tg N yr−1) in the past 4 decades (1980–2020). Direct agricultural emissions in 2020 (3.9 Tg N yr−1, best estimate) represent the large majority of anthropogenic emissions, followed by other direct anthropogenic sources, including fossil fuel and industry, waste and wastewater, and biomass burning (2.1 Tg N yr−1), and indirect anthropogenic sources (1.3 Tg N yr−1) . For the year 2020, our best estimate of total BU emissions for natural and anthropogenic sources was 18.5 (lower–upper bounds: 10.6–27.0) Tg N yr−1, close to our TD estimate of 17.0 (16.6–17.4) Tg N yr−1. For the 2010–2019 period, the annual BU decadal-average emissions for both natural and anthropogenic sources were 18.2 (10.6–25.9) Tg N yr−1 and TD emissions were 17.4 (15.8–19.20) Tg N yr−1. The once top emitter Europe has reduced its emissions by 31 % since the 1980s, while those of emerging economies have grown, making China the top emitter since the 2010s. The observed atmospheric N2O concentrations in recent years have exceeded projected levels under all scenarios in the Coupled Model Intercomparison Project Phase 6 (CMIP6), underscoring the importance of reducing anthropogenic N2O emissions. To evaluate mitigation efforts and contribute to the Global Stocktake of the United Nations Framework Convention on Climate Change, we propose the establishment of a global network for monitoring and modeling N2O from the surface through to the stratosphere. The data presented in this work can be downloaded from https://doi.org/10.18160/RQ8P-2Z4R (Tian et al., 2023).

We present the development and evaluation of MOM6-COBALT-NWA12 version 1.0, a 1/12∘ model of ocean dynamics and biogeochemistry in the northwest Atlantic Ocean. This model is built using the new regional capabilities in the MOM6 ocean model and is coupled with the Carbon, Ocean Biogeochemistry and Lower Trophics (COBALT) biogeochemical model and Sea Ice Simulator version-2 (SIS2) sea ice model. Our goal was to develop a model to provide information to support living-marine-resource applications across management time horizons from seasons to decades. To do this, we struck a balance between a broad, coastwide domain to simulate basin-scale variability and capture cross-boundary issues expected under climate change; a high enough spatial resolution to accurately simulate features like the Gulf Stream separation and advection of water masses through finer-scale coastal features; and the computational economy required to run the long simulations of multiple ensemble members that are needed to quantify prediction uncertainties and produce actionable information. We assess whether MOM6-COBALT-NWA12 is capable of supporting the intended applications by evaluating the model with three categories of metrics: basin-wide indicators of the model's performance, indicators of coastal ecosystem variability and the regional ocean features that drive it, and model run times and computational efficiency. Overall, both the basin-wide and the regional ecosystem-relevant indicators are simulated well by the model. Where notable model biases and errors are present in both types of indicator, they are mainly consistent with the challenges of accurately simulating the Gulf Stream separation, path, and variability: for example, the coastal ocean and shelf north of Cape Hatteras are too warm and salty and have minor biogeochemical biases. During model development, we identified a few model parameters that exerted a notable influence on the model solution, including the horizontal viscosity, mixed-layer restratification, and tidal self-attraction and loading, which we discuss briefly. The computational performance of the model is adequate to support running numerous long simulations, even with the inclusion of coupled biogeochemistry with 40 additional tracers. Overall, these results show that this first version of a regional MOM6 model for the northwest Atlantic Ocean is capable of efficiently and accurately simulating historical basin-wide and regional mean conditions and variability, laying the groundwork for future studies to analyze this variability in detail, develop and improve parameterizations and model components to better capture local ocean features, and develop predictions and projections of future conditions to support living-marine-resource applications across timescales.

The global ocean is losing oxygen with warming. Observations and Earth system model projections, however, suggest that this global ocean deoxygenation does not equate to a simple and systematic expansion of tropical oxygen minimum zones (OMZs). Previous studies have focused on the Pacific Ocean; they showed that the outer OMZ deoxygenates and expands as oxygen supply by advective transport weakens, the OMZ core oxygenates and contracts due to a shift in the composition of the source waters supplied by slow mixing, and in between these two regimes oxygen is redistributed with little effect on OMZ volume. Here, we examine the OMZ response to warming in the Indian Ocean using an ensemble of Earth system model high-emissions scenario experiments from the Coupled Model Intercomparison Project Phase 6. We find a similar expansion–redistribution–contraction response but show that the unique ocean circulation pathways of the Indian Ocean lead to far more prominent OMZ contraction and redistribution regimes than in the Pacific Ocean. As a result, only the outermost volumes (oxygen>180 µmol kg−1) expand. The Indian Ocean experiences a broad oxygenation in the southwest driven by a reduction in waters supplied by the Indonesian Throughflow in favor of high-oxygen waters supplied from the southern Indian Ocean gyre. Models also project a strong localized deoxygenation in the northern Arabian Sea due to the rapid warming and shoaling of marginal sea outflows (Red Sea and Persian Gulf) and increases in local stratification with warming. We extend the existing conceptual framework used to explain the Pacific OMZ response to interpret the response in the Indian Ocean.

A historical dataset of river chemistry and discharge is presented for 140 monitoring sites along the US East Coast, the Gulf of Mexico, and the US West Coast from 1950 to 2022. The dataset, referred to here as River Chemistry for the U.S. Coast (RC4USCoast), is mostly derived from the Water Quality Database of the US Geological Survey (USGS) but also includes river discharge from the USGS's Surface-Water Monthly Statistics for the Nation and the U.S. Army Corps of Engineers. RC4USCoast provides monthly time series as well as long-term averaged monthly climatological patterns for 21 variables including alkalinity and dissolved inorganic carbon concentration. It is mainly intended as a data product for regional ocean biogeochemical models and carbonate chemistry studies in the US coastal regions. Here we present the method to derive RC4USCoast and briefly describe the rivers' carbonate chemistry patterns. The dataset is publicly available at https://doi.org/10.25921/9jfw-ph50 (Gomez et al., 2022).

Marine heat waves (MHWs) are a recurrent phenomenon in the northeast Pacific that impact regional ecosystems and are expected to intensify in the future. Prior work showed that these events, including the 2014–2015 “warm blob”, are associated with widespread surface nutrient declines in the subpolar Alaska Gyre (AG) and the North Pacific Transition Zone (NPTZ) but reduced chlorophyll concentrations in the NPTZ only. Here we explain the contrast between these two regions using a global ocean biogeochemical model (MOM6-COBALT) with Argo float and ship-based observations to investigate how MHWs influence marine productivity. We find that phytoplankton and zooplankton production respond relatively modestly to MHWs in both regions. However, differences in the response to seasonal iron and nitrogen limitation between large (>10 µm) and small (<10 µm) phytoplankton size classes explain the differences in ecosystem response to MHWs across the two biomes. During MHWs, reduced nutrient supply limits large phytoplankton production in the NPTZ (−13 % annually) but has a limited impact on the already iron-limited large phytoplankton population in the AG (−2 %). In contrast, MHWs yield a springtime increase in small phytoplankton in both regions due to shallower mixed layers and weaker light limitation. These modest changes are in apparent contradiction with prior estimates suggesting a collapse in net community production during the warm blob. We show, however, that 70 % of the decline in net community production previously calculated from nitrate Argo data can be attributed to artifacts in the method and that only 30 % can be attributed to interannual variability, in line with our model-based results. Although modest, the primary production anomalies associated with MHWs modify the phytoplankton size distribution, resulting in a significant shift towards small phytoplankton production (i.e., lower large-to-small-phytoplankton ratio) and reduced secondary and export production, especially in the NPTZ.

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 datasets and methodology 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 first time, an approach is shown to reconcile the difference in our ELUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of 10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for 2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLAND was 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. The global atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOS relative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2020, but discrepancies of up to 1 GtC yr−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 changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, 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 datasets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 2021).

The temporal variability of the sea surface partial pressure of CO2 (pCO2) and the underlying processes driving this variability are poorly understood in the coastal ocean. In this study, we tailor an existing method that quantifies the effects of thermal changes, biological activity, ocean circulation and freshwater fluxes to examine seasonal pCO2 changes in highly variable coastal environments. We first use the Modular Ocean Model version 6 (MOM6) and biogeochemical module Carbon Ocean Biogeochemistry And Lower Trophics version 2 (COBALTv2) at a half-degree resolution to simulate coastal CO2 dynamics and evaluate them against pCO2 from the Surface Ocean CO2 Atlas database (SOCAT) and from the continuous coastal pCO2 product generated from SOCAT by a two-step neuronal network interpolation method (coastal Self-Organizing Map Feed-Forward neural Network SOM-FFN, Laruelle et al., 2017). The MOM6-COBALT model reproduces the observed spatiotemporal variability not only in pCO2 but also in sea surface temperature, salinity and nutrients in most coastal environments, except in a few specific regions such as marginal seas. Based on this evaluation, we identify coastal regions of “high” and “medium” agreement between model and coastal SOM-FFN where the drivers of coastal pCO2 seasonal changes can be examined with reasonable confidence. Second, we apply our decomposition method in three contrasted coastal regions: an eastern (US East Coast) and a western (the Californian Current) boundary current and a polar coastal region (the Norwegian Basin). Results show that differences in pCO2 seasonality in the three regions are controlled by the balance between ocean circulation and biological and thermal changes. Circulation controls the pCO2 seasonality in the Californian Current; biological activity controls pCO2 in the Norwegian Basin; and the interplay between biological processes and thermal and circulation changes is key on the US East Coast. The refined approach presented here allows the attribution of pCO2 changes with small residual biases in the coastal ocean, allowing for future work on the mechanisms controlling coastal air–sea CO2 exchanges and how they are likely to be affected by future changes in sea surface temperature, hydrodynamics and biological dynamics.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate – the “global carbon budget” – is important 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 methodology 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) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. 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 last decade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ±  0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budget imbalance BIM of −0.1 GtC yr−1 indicating a near balance between estimated sources and sinks over the last decade. For the year 2019 alone, the growth in EFOS was only about 0.1 % with fossil emissions increasing to 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was 5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminary data for 2020, accounting for the COVID-19-induced changes in emissions, suggest a decrease in EFOS relative to 2019 of about −7 % (median estimate) based on individual estimates from four studies of −6 %, −7 %, −7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. Comparison of estimates from diverse approaches and observations shows (1) no consensus in the mean and trend in land-use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent discrepancy between the different methods for the ocean sink outside the tropics, particularly in the Southern Ocean. 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 (Friedlingstein et al., 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) 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) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. 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 last decade available (2009–2018), EFF was 9.5±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.9±0.02 GtC yr−1 (2.3±0.01 ppm yr−1), SOCEAN 2.5±0.6 GtC yr−1, and SLAND 3.2±0.6 GtC yr−1, with a budget imbalance BIM of 0.4 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in EFF was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC yr−1, reaching 10 GtC yr−1 for the first time in history, ELUC was 1.5±0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5±0.9 GtC yr−1 (42.5±3.3 GtCO2). Also for 2018, GATM was 5.1±0.2 GtC yr−1 (2.4±0.1 ppm yr−1), SOCEAN was 2.6±0.6 GtC yr−1, and SLAND was 3.5±0.7 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in EFF of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. 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 (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013). The data generated by this work are available at https://doi.org/10.18160/gcp-2019 (Friedlingstein et al., 2019).

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use and 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) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. 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 last decade available (2008–2017), EFF was 9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1, SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of 0.5 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2017 alone, the growth in EFF was about 1.6 % and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017, ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 405.0±0.1 ppm averaged over 2017. For 2018, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.7 % (range of 1.8 % to 3.7 %) based on national emission projections for China, the US, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. The analysis presented here shows that the mean and trend in the five components of the global carbon budget are consistently estimated over the period of 1959–2017, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations show (1) no consensus in the mean and trend in land-use change emissions, (2) a persistent low agreement among the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models, originating outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018, 2016, 2015a, b, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2018.

Radium-228 (228Ra), an almost conservative trace isotope in the ocean, supplied from the continental shelves and removed by a known radioactive decay (T1∕2 = 5. 75 years), can be used as a proxy to constrain shelf fluxes of other trace elements, such as nutrients, iron, or rare earth elements. In this study, we perform inverse modeling of a global 228Ra dataset (including GEOSECS, TTO and GEOTRACES programs, and, for the first time, data from the Arctic and around the Kerguelen Islands) to compute the total 228Ra fluxes toward the ocean, using the ocean circulation obtained from the NEMO 3.6 model with a 2° resolution. We optimized the inverse calculation (source regions, cost function) and find a global estimate of the 228Ra fluxes of 8.01–8. 49 × 1023 atoms yr−1, more precise and around 20 % lower than previous estimates. The largest fluxes are in the western North Atlantic, the western Pacific and the Indian Ocean, with roughly two-thirds in the Indo-Pacific Basin. An estimate in the Arctic Ocean is provided for the first time (0.43–0.50 × 1023 atoms yr−1). Local misfits between model and data in the Arctic, the Gulf Stream and the Kuroshio regions could result from flaws of the ocean circulation in these regions (resolution, atmospheric forcing). As radium is enriched in groundwater, a large part of the 228Ra shelf sources comes from submarine groundwater discharge (SGD), a major but poorly known pathway for terrestrial mineral elements, including nutrients, to the ocean. In contrast to the 228Ra budget, the global estimate of SGD is rather unconstrained, between 1.3 and 14. 7 × 1013 m3 yr−1, due to high uncertainties on the other sources of 228Ra, especially diffusion from continental shelf sediments. Better precision on SGD cannot be reached by inverse modeling until a proper way to separate the contributions of SGD and diffusive release from sediments at a global scale is found.

The seasonal upwelling along the west coast of India (WCI) brings nutrient-rich, oxygen-poor subsurface waters to the continental shelf, favoring very low oxygen concentrations in the surface waters during late boreal summer and fall. This yearly-recurring coastal hypoxia is more severe during some years, leading to coastal anoxia that has strong impacts on the living resources. In the present study, we analyze a 1/4° resolution coupled physical–biogeochemical regional oceanic simulation over the 1960–2012 period to investigate the physical processes influencing the oxycline interannual variability off the WCI, that being a proxy for the variability on the shelf in our model. Our analysis indicates a tight relationship between the oxycline and thermocline variations in this region on both seasonal and interannual timescales, thereby revealing a strong physical control of the oxycline variability. As in observations, our model exhibits a shallow oxycline and thermocline during fall that combines with interannual variations to create a window of opportunity for coastal anoxic events. We further demonstrate that the boreal fall oxycline fluctuations off the WCI are strongly related to the Indian Ocean Dipole (IOD), with an asymmetric influence of its positive and negative phases. Positive IODs are associated with easterly wind anomalies near the southern tip of India. These winds force downwelling coastal Kelvin waves that propagate along the WCI and deepen the thermocline and oxycline there, thus preventing the occurrence of coastal anoxia. On the other hand, negative IODs are associated with WCI thermocline and oxycline anomalies of opposite sign but of smaller amplitude, so that the negative or neutral IOD phases are necessary but not the sufficient condition for coastal anoxia. As the IODs generally start developing in summer, these findings suggest some predictability to the occurrence of coastal anoxia off the WCI a couple of months ahead.

Anthropogenic changes in atmosphere–ocean and atmosphere–land CO2 fluxes have been quantified extensively, but few studies have addressed the connection between land and ocean. In this transition zone, the coastal ocean, spatial and temporal data coverage is inadequate to assess its global budget. Thus we use a global ocean biogeochemical model to assess the coastal ocean's global inventory of anthropogenic CO2 and its spatial variability. We used an intermediate resolution, eddying version of the NEMO-PISCES model (ORCA05), varying from 20 to 50 km horizontally, i.e. coarse enough to allow multiple century-scale simulations but finer than coarse-resolution models (∼ 200 km) to better resolve coastal bathymetry and complex coastal currents. Here we define the coastal zone as the continental shelf area, excluding the proximal zone. Evaluation of the simulated air–sea fluxes of total CO2 for 45 coastal regions gave a correlation coefficient R of 0.8 when compared to observation-based estimates. Simulated global uptake of anthropogenic carbon results averaged 2.3 Pg C yr−1 during the years 1993–2012, consistent with previous estimates. Yet only 0.1 Pg C yr−1 of that is absorbed by the global coastal ocean. That represents 4.5 % of the anthropogenic carbon uptake of the global ocean, less than the 7.5 % proportion of coastal-to-global-ocean surface areas. Coastal uptake is weakened due to a bottleneck in offshore transport, which is inadequate to reduce the mean anthropogenic carbon concentration of coastal waters to the mean level found in the open-ocean mixed layer.

During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to-model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.

This study aims to evaluate the potential for impacts of ocean acidification on North Atlantic deep-sea ecosystems in response to IPCC AR5 Representative Concentration Pathways (RCPs). Deep-sea biota is likely highly vulnerable to changes in seawater chemistry and sensitive to moderate excursions in pH. Here we show, from seven fully coupled Earth system models, that for three out of four RCPs over 17% of the seafloor area below 500 m depth in the North Atlantic sector will experience pH reductions exceeding −0.2 units by 2100. Increased stratification in response to climate change partially alleviates the impact of ocean acidification on deep benthic environments. We report on major pH reductions over the deep North Atlantic seafloor (depth >500 m) and at important deep-sea features, such as seamounts and canyons. By 2100, and under the high CO2 scenario RCP8.5, pH reductions exceeding −0.2 (−0.3) units are projected in close to 23% (~15%) of North Atlantic deep-sea canyons and ~8% (3%) of seamounts – including seamounts proposed as sites of marine protected areas. The spatial pattern of impacts reflects the depth of the pH perturbation and does not scale linearly with atmospheric CO2 concentration. Impacts may cause negative changes of the same magnitude or exceeding the current target of 10% of preservation of marine biomes set by the convention on biological diversity, implying that ocean acidification may offset benefits from conservation/management strategies relying on the regulation of resource exploitation.

The considerable uncertainties in the carbon budget of the Southern Ocean are largely attributed to unresolved variability, in particular at a seasonal timescale and small spatial scale (~ 100 km). In this study, the variability of surface pCO2 and dissolved inorganic carbon (DIC) at seasonal and small spatial scales is examined using a data set of surface drifters including ~ 80 000 measurements at high spatiotemporal resolution. On spatial scales of 100 km, we find gradients ranging from 5 to 50 μatm for pCO2 and 2 to 30 μmol kg−1 for DIC, with highest values in energetic and frontal regions. This result is supported by a second estimate obtained with sea surface temperature (SST) satellite images and local DIC–SST relationships derived from drifter observations. We find that dynamical processes drive the variability of DIC at small spatial scale in most regions of the Southern Ocean and the cascade of large-scale gradients down to small spatial scales, leading to gradients up to 15 μmol kg−1 over 100 km. Although the role of biological activity is more localized, it enhances the variability up to 30 μmol kg−1 over 100 km. The seasonal cycle of surface DIC is reconstructed following Mahadevan et al. (2011), using an annual climatology of DIC and a monthly climatology of mixed layer depth. This method is evaluated using drifter observations and proves to be a reasonable first-order estimate of the seasonality in the Southern Ocean that could be used to validate model simulations. We find that small spatial-scale structures are a non-negligible source of variability for DIC, with amplitudes of about a third of the variations associated with the seasonality and up to 10 times the magnitude of large-scale gradients. The amplitude of small-scale variability reported here should be kept in mind when inferring temporal changes (seasonality, interannual variability, decadal trends) of the carbon budget from low-resolution observations and models.

Ocean ecosystems are increasingly stressed by human-induced changes of their physical, chemical and biological environment. Among these changes, warming, acidification, deoxygenation and changes in primary productivity by marine phytoplankton can be considered as four of the major stressors of open ocean ecosystems. Due to rising atmospheric CO2 in the coming decades, these changes will be amplified. Here, we use the most recent simulations performed in the framework of the Coupled Model Intercomparison Project 5 to assess how these stressors may evolve over the course of the 21st century. The 10 Earth system models used here project similar trends in ocean warming, acidification, deoxygenation and reduced primary productivity for each of the IPCC's representative concentration pathways (RCPs) over the 21st century. For the "business-as-usual" scenario RCP8.5, the model-mean changes in the 2090s (compared to the 1990s) for sea surface temperature, sea surface pH, global O2 content and integrated primary productivity amount to &plus;2.73 (±0.72) °C, −0.33 (±0.003) pH unit, −3.45 (±0.44)% and −8.6 (±7.9)%, respectively. For the high mitigation scenario RCP2.6, corresponding changes are +0.71 (±0.45) °C, −0.07 (±0.001) pH unit, −1.81 (±0.31)% and −2.0 (±4.1)%, respectively, illustrating the effectiveness of extreme mitigation strategies. Although these stressors operate globally, they display distinct regional patterns and thus do not change coincidentally. Large decreases in O2 and in pH are simulated in global ocean intermediate and mode waters, whereas large reductions in primary production are simulated in the tropics and in the North Atlantic. Although temperature and pH projections are robust across models, the same does not hold for projections of subsurface O2 concentrations in the tropics and global and regional changes in net primary productivity. These high uncertainties in projections of primary productivity and subsurface oxygen prompt us to continue inter-model comparisons to understand these model differences, while calling for caution when using the CMIP5 models to force regional impact models.

The expansion of OMZs (oxygen minimum zones) due to climate change and their possible evolution and impacts on the ecosystems and the atmosphere are still debated, mostly because of the unability of global climate models to adequatly reproduce the processes governing OMZs. In this study, we examine the factors controlling the oxygen budget, i.e. the equilibrium between oxygen sources and sinks in the northern Arabian Sea OMZ using an eddy-resolving biophysical model. Our model confirms that the biological consumption of oxygen is most intense below the region of highest productivity in the western Arabian Sea. The oxygen drawdown in this region is counterbalanced by the large supply of oxygenated waters originated from the south and advected horizontally by the western boundary current. Although the biological sink and the dynamical sources of oxygen compensate on annual average, we find that the seasonality of the dynamical transport of oxygen is 3 to 5 times larger than the seasonality of the biological sink. In agreement with previous findings, the resulting seasonality of oxygen concentration in the OMZ is relatively weak, with a variability of the order of 15% of the annual mean oxygen concentration in the oxycline and 5% elsewhere. This seasonality primarily arises from the vertical displacement of the OMZ forced by the monsoonal reversal of Ekman pumping across the basin. In coastal areas, the oxygen concentration is also modulated seasonally by lateral advection. Along the western coast of the Arabian Sea, the Somali Current transports oxygen-rich waters originated from the south during summer and oxygen-poor waters from the northeast during winter. Along the eastern coast of the Arabian Sea, we find that the main contributor to lateral advection in the OMZ is the Indian coastal undercurrent that advects southern oxygenated waters during summer and northern low-oxygen waters during winter. In this region, our model indicates that oxygen concentrations are modulated seasonally by coastal Kelvin waves and westward-propagating Rossby waves. Whereas on seasonal time scales the sources and sinks of oxygen are dominated by the mean vertical and lateral advection (Ekman pumping and monsoonal currents), on annual time scales we find that the biological sink is counterbalanced by the supply of oxygen sustained by mesoscale structures (eddies and filaments). Eddy-driven advection hence promotes the vertical supply of oxygen along the western coast of the Arabian Sea and the lateral transport of ventilated waters offshore the coast of Oman and southwest India.

Low-latitude (LL) oceans account for up to half of global net primary production and export 1-5 . It has been argued that the Southern Ocean dominates LL primary production and export 6 , with implications for the response of global primary production and export to climate change 7 . Here we applied observational analyses and sensitivity studies to an individual model to show, instead, that 72% of LL primary production and 55% of export is controlled by local mesopelagic macronutrient cycling. A total of 34% of the LL export is sustained by preformed macronutrients supplied from the Southern Ocean via a deeper overturning cell, with a shallow preformed northward supply, crossing 30° S through subpolar and thermocline water masses, sustaining only 7% of the LL export. Analyses of five Coupled Model Intercomparison Project Phase 6 (CMIP6) models, run under both high-emissions low-mitigation (shared socioeconomic pathway (SSP5-8.5)) and low-emissions high-mitigation (SSP1-2.6) climate scenarios for 1850-2300, revealed significant across-model disparities in their projections of not only the amplitude, but also the sign, of LL primary production. Under the stronger SSP5-8.5 forcing, with more substantial upper-ocean warming, the CMIP6 models that account for temperature-dependent remineralization promoted enhanced LL mesopelagic nutrient retention under warming, with this providing a first-order contribution to stabilizing or increasing, rather than decreasing, LL production under high emissions and low mitigation. This underscores the importance of a mechanistic understanding of mesopelagic remineralization and its sensitivity to ocean warming for predicting future ecosystem changes., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

The Southern Ocean is a major sink of atmospheric CO 2 , but the nature and magnitude of its variability remains uncertain and debated. Estimates based on observations suggest substantial variability that is not reproduced by process-based ocean models, with increasingly divergent estimates over the past decade. We examine potential constraints on the nature and magnitude of climate-driven variability of the Southern Ocean CO 2 sink from observation-based air-sea O 2 fluxes. On interannual time scales, the variability in the air-sea fluxes of CO 2 and O 2 estimated from observations is consistent across the two species and positively correlated with the variability simulated by ocean models. Our analysis suggests that variations in ocean ventilation related to the Southern Annular Mode are responsible for this interannual variability. On decadal time scales, the existence of significant variability in the air-sea CO 2 flux estimated from observations also tends to be supported by observation-based estimates of O 2 flux variability. However, the large decadal variability in air-sea CO 2 flux is absent from ocean models. Our analysis suggests that issues in representing the balance between the thermal and non-thermal components of the CO 2 sink and/or insufficient variability in mode water formation might contribute to the lack of decadal variability in the current generation of ocean models. This article is part of a discussion meeting issue 'Heat and carbon uptake in the Southern Ocean: the state of the art and future priorities'.

Agricultural intensification in India has increased nitrogen pollution, leading to water quality impairments. The fate of reactive nitrogen applied to the land is largely unknown, however. Long-term records of riverine nitrogen fluxes are nonexistent and drivers of variability remain unexamined, limiting the development of nitrogen management strategies. Here, we leverage dissolved inorganic nitrogen (DIN) and discharge data to characterize the seasonal, annual, and regional variability of DIN fluxes and their drivers for seven major river basins from 1981 to 2014. We find large seasonal and interannual variability in nitrogen runoff, with 68% to 94% of DIN fluxes occurring in June through October and with the coefficient of variation across years ranging from 44% to 93% for individual basins. This variability is primarily explained by variability in precipitation, with year- and basin-specific annual precipitation explaining 52% of the combined regional and interannual variability. We find little correlation with rising fertilizer application rates in five of the seven basins, implying that agricultural intensification has thus far primarily impacted groundwater and atmospheric emissions rather than riverine runoff. These findings suggest that riverine nitrogen runoff in India is highly sensitive to projected future increases in precipitation and intensification of the seasonal monsoon, while the impact of projected continued land use intensification is highly uncertain.

Carbon storage by the ocean and by the land is usually quantified separately, and does not fully take into account the land-to-ocean transport of carbon through inland waters, estuaries, tidal wetlands and continental shelf waters-the 'land-to-ocean aquatic continuum' (LOAC). Here we assess LOAC carbon cycling before the industrial period and perturbed by direct human interventions, including climate change. In our view of the global carbon cycle, the traditional 'long-range loop', which carries carbon from terrestrial ecosystems to the open ocean through rivers, is reinforced by two 'short-range loops' that carry carbon from terrestrial ecosystems to inland waters and from tidal wetlands to the open ocean. Using a mass-balance approach, we find that the pre-industrial uptake of atmospheric carbon dioxide by terrestrial ecosystems transferred to the ocean and outgassed back to the atmosphere amounts to 0.65 ± 0.30 petagrams of carbon per year (±2 sigma). Humans have accelerated the cycling of carbon between terrestrial ecosystems, inland waters and the atmosphere, and decreased the uptake of atmospheric carbon dioxide from tidal wetlands and submerged vegetation. Ignoring these changing LOAC carbon fluxes results in an overestimation of carbon storage in terrestrial ecosystems by 0.6 ± 0.4 petagrams of carbon per year, and an underestimation of sedimentary and oceanic carbon storage. We identify knowledge gaps that are key to reduce uncertainties in future assessments of LOAC fluxes., (© 2022. Springer Nature Limited.)

The ocean is the main source of thermal inertia in the climate system. Ocean heat uptake during recent decades has been quantified using ocean temperature measurements. However, these estimates all use the same imperfect ocean dataset and share additional uncertainty due to sparse coverage, especially before 2007. Here, we provide an independent estimate by using measurements of atmospheric oxygen (O 2 ) and carbon dioxide (CO 2 ) - levels of which increase as the ocean warms and releases gases - as a whole ocean thermometer. We show that the ocean gained 1.29 ± 0.79 × 10 22 Joules of heat per year between 1991 and 2016, equivalent to a planetary energy imbalance of 0.80 ± 0.49 W watts per square metre of Earth's surface. We also find that the ocean-warming effect that led to the outgassing of O 2 and CO 2 can be isolated from the direct effects of anthropogenic emissions and CO 2 sinks. Our result - which relies on high-precision O 2 atmospheric measurements dating back to 1991 - leverages an integrative Earth system approach and provides much needed independent confirmation of heat uptake estimated from ocean data.

New estimates of p CO 2 from profiling floats deployed by the Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) project have demonstrated the importance of wintertime outgassing south of the Polar Front, challenging the accepted magnitude of Southern Ocean carbon uptake (Gray et al., 2018, https://doi:10.1029/2018GL078013). Here, we put 3.5 years of SOCCOM observations into broader context with the global surface carbon dioxide database (Surface Ocean CO 2 Atlas, SOCAT) by using the two interpolation methods currently used to assess the ocean models in the Global Carbon Budget (Le Quéré et al., 2018, https://doi:10.5194/essd-10-2141-2018) to create a ship-only, a float-weighted, and a combined estimate of Southern Ocean carbon fluxes (<35°S). In our ship-only estimate, we calculate a mean uptake of -1.14 ± 0.19 Pg C/yr for 2015-2017, consistent with prior studies. The float-weighted estimate yields a significantly lower Southern Ocean uptake of -0.35 ± 0.19 Pg C/yr. Subsampling of high-resolution ocean biogeochemical process models indicates that some of the differences between float and ship-only estimates of the Southern Ocean carbon flux can be explained by spatial and temporal sampling differences. The combined ship and float estimate minimizes the root-mean-square p CO 2 difference between the mapped product and both data sets, giving a new Southern Ocean uptake of -0.75 ± 0.22 Pg C/yr, though with uncertainties that overlap the ship-only estimate. An atmospheric inversion reveals that a shift of this magnitude in the contemporary Southern Ocean carbon flux must be compensated for by ocean or land sinks within the Southern Hemisphere., (©2019. American Geophysical Union. All Rights Reserved.)

The ocean is the main source of thermal inertia in the climate system 1 . During recent decades, ocean heat uptake has been quantified by using hydrographic temperature measurements and data from the Argo float program, which expanded its coverage after 2007 2,3 . However, these estimates all use the same imperfect ocean dataset and share additional uncertainties resulting from sparse coverage, especially before 2007 4,5 . Here we provide an independent estimate by using measurements of atmospheric oxygen (O 2 ) and carbon dioxide (CO 2 )-levels of which increase as the ocean warms and releases gases-as a whole-ocean thermometer. We show that the ocean gained 1.33 ± 0.20  × 10 22  joules of heat per year between 1991 and 2016, equivalent to a planetary energy imbalance of 0.83 ± 0.11 watts per square metre of Earth's surface. We also find that the ocean-warming effect that led to the outgassing of O 2 and CO 2 can be isolated from the direct effects of anthropogenic emissions and CO 2 sinks. Our result-which relies on high-precision O 2 measurements dating back to 1991 6 -suggests that ocean warming is at the high end of previous estimates, with implications for policy-relevant measurements of the Earth response to climate change, such as climate sensitivity to greenhouse gases 7 and the thermal component of sea-level rise 8 .

A decrease in the 13 C/ 12 C ratio of atmospheric CO 2 has been documented by direct observations since 1978 and from ice core measurements since the industrial revolution. This decrease, known as the 13 C-Suess effect, is driven primarily by the input of fossil fuel-derived CO 2 but is also sensitive to land and ocean carbon cycling and uptake. Using updated records, we show that no plausible combination of sources and sinks of CO 2 from fossil fuel, land, and oceans can explain the observed 13 C-Suess effect unless an increase has occurred in the 13 C/ 12 C isotopic discrimination of land photosynthesis. A trend toward greater discrimination under higher CO 2 levels is broadly consistent with tree ring studies over the past century, with field and chamber experiments, and with geological records of C 3 plants at times of altered atmospheric CO 2 , but increasing discrimination has not previously been included in studies of long-term atmospheric 13 C/ 12 C measurements. We further show that the inferred discrimination increase of 0.014 ± 0.007‰ ppm -1 is largely explained by photorespiratory and mesophyll effects. This result implies that, at the global scale, land plants have regulated their stomatal conductance so as to allow the CO 2 partial pressure within stomatal cavities and their intrinsic water use efficiency to increase in nearly constant proportion to the rise in atmospheric CO 2 concentration., Competing Interests: The authors declare no conflict of interest.

All Earth System models project a consistent decrease in the oxygen content of oceans for the coming decades because of ocean warming, reduced ventilation and increased stratification. But large uncertainties for these future projections of ocean deoxygenation remain for the subsurface tropical oceans where the major oxygen minimum zones are located. Here, we combine global warming projections, model-based estimates of natural short-term variability, as well as data and model estimates of the Last Glacial Maximum (LGM) ocean oxygenation to gain some insights into the major mechanisms of oxygenation changes across these different time scales. We show that the primary uncertainty on future ocean deoxygenation in the subsurface tropical oceans is in fact controlled by a robust compensation between decreasing oxygen saturation (O 2sat ) due to warming and decreasing apparent oxygen utilization (AOU) due to increased ventilation of the corresponding water masses. Modelled short-term natural variability in subsurface oxygen levels also reveals a compensation between O 2sat and AOU, controlled by the latter. Finally, using a model simulation of the LGM, reproducing data-based reconstructions of past ocean (de)oxygenation, we show that the deoxygenation trend of the subsurface ocean during deglaciation was controlled by a combination of warming-induced decreasing O 2sat and increasing AOU driven by a reduced ventilation of tropical subsurface waters.This article is part of the themed issue 'Ocean ventilation and deoxygenation in a warming world'., (© 2017 The Author(s).)

Global ocean acidification is caused primarily by the ocean's uptake of CO2 as a consequence of increasing atmospheric CO2 levels. We present observations of the oceanic decrease in pH at the basin scale (50 °S-36 °N) for the Atlantic Ocean over two decades (1993-2013). Changes in pH associated with the uptake of anthropogenic CO2 (ΔpHCant) and with variations caused by biological activity and ocean circulation (ΔpHNat) are evaluated for different water masses. Output from an Institut Pierre Simon Laplace climate model is used to place the results into a longer-term perspective and to elucidate the mechanisms responsible for pH change. The largest decreases in pH (∆pH) were observed in central, mode, and intermediate waters, with a maximum ΔpH value in South Atlantic Central Waters of -0.042 ± 0.003. The ΔpH trended toward zero in deep and bottom waters. Observations and model results show that pH changes generally are dominated by the anthropogenic component, which accounts for rates between -0.0015 and -0.0020/y in the central waters. The anthropogenic and natural components are of the same order of magnitude and reinforce one another in mode and intermediate waters over the time period. Large negative ΔpHNat values observed in mode and intermediate waters are driven primarily by changes in CO2 content and are consistent with (i) a poleward shift of the formation region during the positive phase of the Southern Annular Mode in the South Atlantic and (ii) an increase in the rate of the water mass formation in the North Atlantic.