Description: 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)

Description: We present a mapped climatology (GLODAPv2.2016b) of ocean biogeochemical variables based on the new GLODAP version 2 data product (Olsen et al., 2016; Key et al., 2015), which covers all ocean basins over the years 1972 to 2013. The quality-controlled and internally consistent GLODAPv2 was used to create global 1°  ×  1° mapped climatologies of salinity, temperature, oxygen, nitrate, phosphate, silicate, total dissolved inorganic carbon (TCO2), total alkalinity (TAlk), pH, and CaCO3 saturation states using the Data-Interpolating Variational Analysis (DIVA) mapping method. Improving on maps based on an earlier but similar dataset, GLODAPv1.1, this climatology also covers the Arctic Ocean. Climatologies were created for 33 standard depth surfaces. The conceivably confounding temporal trends in TCO2 and pH due to anthropogenic influence were removed prior to mapping by normalizing these data to the year 2002 using first-order calculations of anthropogenic carbon accumulation rates. We additionally provide maps of accumulated anthropogenic carbon in the year 2002 and of preindustrial TCO2. For all parameters, all data from the full 1972–2013 period were used, including data that did not receive full secondary quality control. The GLODAPv2.2016b global 1°  ×  1° mapped climatologies, including error fields and ancillary information, are available at the GLODAPv2 web page at the Carbon Dioxide Information Analysis Center (CDIAC; doi:10.3334/CDIAC/OTG.NDP093_GLODAPv2).

Description: Stratospheric ozone depletion and emission of greenhouse gases lead to a trend of the southern annular mode (SAM) toward its high-index polarity. The positive phase of the SAM is characterized by stronger than usual westerly winds that induce changes in the physical carbon transport. Changes in the natural carbon budget of the upper 100 m of the Southern Ocean in response to a positive SAM phase are explored with a coupled ecosystem-general circulation model and regression analysis. Previously overlooked processes that are important for the upper ocean carbon budget during a positive SAM period are identified, namely, export production and downward transport of carbon north of the polar front (PF) as large as the upwelling in the south. The limiting micronutrient iron is brought into the surface layer by upwelling and stimulates phytoplankton growth and export production but only in summer. This leads to a drawdown of carbon and less summertime outgassing (or more uptake) of natural CO2. In winter, biological mechanisms are inactive, and the surface ocean equilibrates with the atmosphere by releasing CO2. In the annual mean, the upper ocean region south of the PF loses more carbon by additional export production than by the release of CO2 into the atmosphere, highlighting the role of the biological carbon pump in response to a positive SAM event.

Description: Version 2 of the Global Ocean Data Analysis Project (GLODAPv2) data product is composed of data from 724 scientific cruises covering the global ocean. It includes data assembled during the previous efforts GLODAPv1.1 (Global Ocean Data Analysis Project version 1.1) in 2004, CARINA (CARbon IN the Atlantic) in 2009/2010, and PACIFICA (PACIFic ocean Interior CArbon) in 2013, as well as data from an additional 168 cruises. Data for 12 core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl 4) have been subjected to extensive quality control, including systematic evaluation of bias. The data are available in two formats: (i) as submitted but updated to WOCE exchange format and (ii) as a merged and internally consistent data product. In the latter, adjustments have been applied to remove significant biases, respecting occurrences of any known or likely time trends or variations. Adjustments applied by previous efforts were re-evaluated. Hence, GLODAPv2 is not a simple merging of previous products with some new data added but a unique, internally consistent data product. This compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg −1 in dissolved inorganic carbon, 6 µmol kg −1 in total alkalinity, 0.005 in pH, and 5 % for the halogenated transient tracers.

Description: 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 all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2006–2015), EFF was 9.3 ± 0.5 GtC yr−1, ELUC 1.0 ± 0.5 GtC yr−1, GATM 4.5 ± 0.1 GtC yr−1, SOCEAN 2.6 ± 0.5 GtC yr−1, and SLAND 3.1 ± 0.9 GtC yr−1. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1, showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr−1 that took place during 2006–2015. Also, for 2015, ELUC was 1.3 ± 0.5 GtC yr−1, GATM was 6.3 ± 0.2 GtC yr−1, SOCEAN was 3.0 ± 0.5 GtC yr−1, and SLAND was 1.9 ± 0.9 GtC yr−1. GATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4 ± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2 % (range of −1.0 to +1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Niño conditions of 2015–2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565 ± 55 GtC (2075 ± 205 GtCO2) for 1870–2016, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2016).

Description: We combined data sets of measured sedimentary calcium carbonate (CaCO3) and satellite-derived pelagic primary production to parameterize the relation between CaCO3 content on the Antarctic shelves and primary production in the overlying water column. CaCO3 content predicted in this way was in good agreement with the measured data. The parameterization was then used to chart CaCO3 content on the Antarctic shelves all around the Antarctic, using the satellite-derived primary production. The total inventory of CaCO3 in the bioturbated layer of Antarctic shelf sediments was estimated to be 0.5 Pg C. This quantity is comparable to the total CO2 uptake by the Southern Ocean in only one to a few years (dependent on the uptake estimate and area considered), indicating that the dissolution of these carbonates will neither delay ocean acidification in this area nor augment the Southern Ocean CO2 uptake capacity.

Description: We present measurements of pCO2, O2 concentration, biological oxygen saturation (ΔO2/Ar), and N2 saturation (ΔN2) in Southern Ocean surface waters during austral summer, 2010–2011. Phytoplankton biomass varied strongly across distinct hydrographic zones, with high chlorophyll a (Chl a) concentrations in regions of frontal mixing and sea ice melt. pCO2 and ΔO2/Ar exhibited large spatial gradients (range 90 to 450 µatm and −10 to 60%, respectively) and covaried strongly with Chl a. However, the ratio of biological O2 accumulation to dissolved inorganic carbon (DIC) drawdown was significantly lower than expected from photosynthetic stoichiometry, reflecting the differential time scales of O2 and CO2 air-sea equilibration. We measured significant oceanic CO2 uptake, with a mean air-sea flux (~ −10 mmol m−2 d−1) that significantly exceeded regional climatological values. N2 was mostly supersaturated in surface waters (mean ΔN2 of +2.5%), while physical processes resulted in both supersaturation and undersaturation of mixed layer O2 (mean ΔO2phys = 2.1%). Box model calculations were able to reproduce much of the spatial variability of ΔN2 and ΔO2phys along the cruise track, demonstrating significant effects of air-sea exchange processes (e.g., atmospheric pressure changes and bubble injection) and mixed layer entrainment on surface gas disequilibria. Net community production (NCP) derived from entrainment-corrected surface ΔO2/Ar data, ranged from ~ −40 to 〉 300 mmol O2 m−2 d−1 and showed good coherence with independent NCP estimates based on seasonal mixed layer DIC deficits. Elevated NCP was observed in hydrographic frontal zones and stratified regions of sea ice melt, reflecting physical controls on surface water light fields and nutrient availability.

Description: The Surface Ocean CO2 Atlas (SOCAT) is a synthesis of quality-controlled fCO2 (fugacity of carbon dioxide) values for the global surface oceans and coastal seas with regular updates. Version 3 of SOCAT has 14.7 million fCO2 values from 3646 data sets covering the years 1957 to 2014. This latest version has an additional 4.6 million fCO2 values relative to version 2 and extends the record from 2011 to 2014. Version 3 also significantly increases the data availability for 2005 to 2013. SOCAT has an average of approximately 1.2 million surface water fCO2 values per year for the years 2006 to 2012. Quality and documentation of the data has improved. A new feature is the data set quality control (QC) flag of E for data from alternative sensors and platforms. The accuracy of surface water fCO2 has been defined for all data set QC flags. Automated range checking has been carried out for all data sets during their upload into SOCAT. The upgrade of the interactive Data Set Viewer (previously known as the Cruise Data Viewer) allows better interrogation of the SOCAT data collection and rapid creation of high-quality figures for scientific presentations. Automated data upload has been launched for version 4 and will enable more frequent SOCAT releases in the future. High-profile scientific applications of SOCAT include quantification of the ocean sink for atmospheric carbon dioxide and its long-term variation, detection of ocean acidification, as well as evaluation of coupled-climate and ocean-only biogeochemical models. Users of SOCAT data products are urged to acknowledge the contribution of data providers, as stated in the SOCAT Fair Data Use Statement. This ESSD (Earth System Science Data) "living data" publication documents the methods and data sets used for the assembly of this new version of the SOCAT data collection and compares these with those used for earlier versions of the data collection (Pfeil et al., 2013; Sabine et al., 2013; Bakker et al., 2014). Individual data set files, included in the synthesis product, can be downloaded here: doi:10.1594/PANGAEA.849770. The gridded products are available here: doi:10.3334/CDIAC/OTG.SOCAT_V3_GRID.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2020 is an update of the previous version, GLODAPv2.2019. The major changes are data from 106 new cruises added, extension of time coverage to 2019, and the inclusion of available (also for historical cruises) discrete fugacity of CO2 (fCO2) values in the merged product files. GLODAPv2.2020 now includes measurements from more than 1.2 million water samples from the global oceans collected on 946 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 106 new cruises with the data from the 840 quality-controlled cruises of the GLODAPv2.2019 data product using crossover analysis. Comparisons to empirical algorithm estimates provided additional context for adjustment decisions; this is new to this version. The adjustments are intended to remove potential biases from errors related to measurement, calibration, and data-handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 µmol kg−1 in dissolved inorganic carbon, 4 µmol kg−1 in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete fCO2, were not subjected to bias comparison or adjustments.

Description: Anthropogenic emissions of CO2 to the atmosphere have modified the carbon cycle for more than 2 centuries. As the ocean stores most of the carbon on our planet, there is an important task in unraveling the natural and anthropogenic processes that drive the carbon cycle at different spatial and temporal scales. We contribute to this by designing a global monthly climatology of total dissolved inorganic carbon (TCO2), which offers a robust basis in carbon cycle modeling but also for other studies related to this cycle. A feedforward neural network (dubbed NNGv2LDEO) was configured to extract from the Global Ocean Data Analysis Project version 2.2019 (GLODAPv2.2019) and the Lamont–Doherty Earth Observatory (LDEO) datasets the relations between TCO2 and a set of variables related to the former's variability. The global root mean square error (RMSE) of mapping TCO2 is relatively low for the two datasets (GLODAPv2.2019: 7.2 µmol kg−1; LDEO: 11.4 µmol kg−1) and also for independent data, suggesting that the network does not overfit possible errors in data. The ability of NNGv2LDEO to capture the monthly variability of TCO2 was testified through the good reproduction of the seasonal cycle in 10 time series stations spread over different regions of the ocean (RMSE: 3.6 to 13.2 µmol kg−1). The climatology was obtained by passing through NNGv2LDEO the monthly climatological fields of temperature, salinity, and oxygen from the World Ocean Atlas 2013 and phosphate, nitrate, and silicate computed from a neural network fed with the previous fields. The resolution is 1∘×1∘ in the horizontal, 102 depth levels (0–5500 m), and monthly (0–1500 m) to annual (1550–5500 m) temporal resolution, and it is centered around the year 1995. The uncertainty of the climatology is low when compared with climatological values derived from measured TCO2 in the largest time series stations. Furthermore, a computed climatology of partial pressure of CO2 (pCO2) from a previous climatology of total alkalinity and the present one of TCO2 supports the robustness of this product through the good correlation with a widely used pCO2 climatology (Landschützer et al., 2017). Our TCO2 climatology is distributed through the data repository of the Spanish National Research Council (CSIC; https://doi.org/10.20350/digitalCSIC/10551, Broullón et al., 2020).