Accurate assessment of anthropogenic carbon dioxide (CO₂) 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 CO₂ emissions (E_FOS) are based on energy statistics and cement production data, while emissions from land-use change (E_LUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO₂ concentration is measured directly, and its growth rate (G_ATM) is computed from the annual changes in concentration. The ocean CO₂ sink (S_OCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO₂ sink (S_LAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (B_IM), 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 E_LUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, E_FOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 ± 0.5 GtC yr⁻¹ (9.3 ± 0.5 GtC yr⁻¹ when the cement carbonation sink is included), and E_LUC was 0.9 ± 0.7 GtC yr⁻¹, for a total anthropogenic CO₂ emission of 10.2 ± 0.8 GtC yr⁻¹ (37.4 ± 2.9 GtCO₂). Also, for 2020, G_ATM was 5.0 ± 0.2 GtC yr⁻¹ (2.4 ± 0.1 ppm yr⁻¹), S_OCEAN was 3.0 ± 0.4 GtC yr⁻¹, and S_LAND was 2.9 ± 1 GtC yr⁻¹, with a B_IM of −0.8 GtC yr⁻¹. The global atmospheric CO₂ concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in E_FOS 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⁻¹ persist for the representation of annual to semi-decadal variability in CO₂ 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 CO₂ 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 land and ocean absorb on average just over half of the anthropogenic emissions of carbon dioxide (CO₂) every year. These CO₂ "sinks" are modulated by climate change and variability. Here we use a suite of nine dynamic global vegetation models (DGVMs) and four ocean biogeochemical general circulation models (OBGCMs) to estimate trends driven by global and regional climate and atmospheric CO₂ in land and oceanic CO₂ exchanges with the atmosphere over the period 1990–2009, to attribute these trends to underlying processes in the models, and to quantify the uncertainty and level of inter-model agreement. The models were forced with reconstructed climate fields and observed global atmospheric CO₂; land use and land cover changes are not included for the DGVMs. Over the period 1990–2009, the DGVMs simulate a mean global land carbon sink of −2.4 ± 0.7 Pg C yr⁻¹ with a small significant trend of −0.06 ± 0.03 Pg C yr⁻² (increasing sink). Over the more limited period 1990–2004, the ocean models simulate a mean ocean sink of −2.2 ± 0.2 Pg C yr⁻¹ with a trend in the net C uptake that is indistinguishable from zero (−0.01 ± 0.02 Pg C yr⁻²). The two ocean models that extended the simulations until 2009 suggest a slightly stronger, but still small, trend of −0.02 ± 0.01 Pg C yr⁻². Trends from land and ocean models compare favourably to the land greenness trends from remote sensing, atmospheric inversion results, and the residual land sink required to close the global carbon budget. Trends in the land sink are driven by increasing net primary production (NPP), whose statistically significant trend of 0.22 ± 0.08 Pg C yr⁻² exceeds a significant trend in heterotrophic respiration of 0.16 ± 0.05 Pg C yr⁻² – primarily as a consequence of widespread CO₂ fertilisation of plant production. Most of the land-based trend in simulated net carbon uptake originates from natural ecosystems in the tropics (−0.04 ± 0.01 Pg C yr⁻²), with almost no trend over the northern land region, where recent warming and reduced rainfall offsets the positive impact of elevated atmospheric CO₂ and changes in growing season length on carbon storage. The small uptake trend in the ocean models emerges because climate variability and change, and in particular increasing sea surface temperatures, tend to counter\-act the trend in ocean uptake driven by the increase in atmospheric CO₂. Large uncertainty remains in the magnitude and sign of modelled carbon trends in several regions, as well as regarding the influence of land use and land cover changes on regional trends.

The Surface Ocean CO₂ Atlas (SOCAT) is a synthesis of quality-controlled fCO₂ (fugacity of carbon dioxide) values for the global surface oceans and coastal seas with regular updates. Version 3 of SOCAT has 14.7 million fCO₂ values from 3646 data sets covering the years 1957 to 2014. This latest version has an additional 4.6 million fCO₂ 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 fCO₂ 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 fCO₂ 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.

Accurate assessment of anthropogenic carbon dioxide (CO₂) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere 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 a 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, consistency within and among components, alongside methodology and data limitations. CO₂ emissions from fossil fuel combustion and cement production (E_FF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (E_LUC), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO₂ concentration is measured directly and its rate of growth (G_ATM) is computed from the annual changes in concentration. The mean ocean CO₂ sink (S_OCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in S_OCEAN is evaluated with data products based on surveys of ocean CO₂ measurements. The global residual terrestrial CO₂ sink (S_LAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO₂, and land-cover-change (some including nitrogen–carbon interactions). 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 (2004–2013), E_FF was 8.9 ± 0.4 GtC yr⁻¹, E_LUC 0.9 ± 0.5 GtC yr⁻¹, G_ATM 4.3 ± 0.1 GtC yr⁻¹, S_OCEAN 2.6 ± 0.5 GtC yr⁻¹, and S_LAND 2.9 ± 0.8 GtC yr⁻¹. For year 2013 alone, E_FF grew to 9.9 ± 0.5 GtC yr⁻¹, 2.3% above 2012, continuing the growth trend in these emissions, E_LUC was 0.9 ± 0.5 GtC yr⁻¹, G_ATM was 5.4 ± 0.2 GtC yr⁻¹, S_OCEAN was 2.9 ± 0.5 GtC yr⁻¹, and S_LAND was 2.5 ± 0.9 GtC yr⁻¹. G_ATM was high in 2013, reflecting a steady increase in E_FF and smaller and opposite changes between S_OCEAN and S_LAND compared to the past decade (2004–2013). The global atmospheric CO₂ concentration reached 395.31 ± 0.10 ppm averaged over 2013. We estimate that E_FF will increase by 2.5% (1.3–3.5%) to 10.1 ± 0.6 GtC in 2014 (37.0 ± 2.2 GtCO₂ yr⁻¹), 65% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the global economy. From this projection of E_FF and assumed constant E_LUC for 2014, cumulative emissions of CO₂ will reach about 545 ± 55 GtC (2000 ± 200 GtCO₂) for 1870–2014, about 75% from E_FF and 25% from E_LUC. This paper documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this living data set (Le Quéré et al., 2013, 2014). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2014).

Accurate assessment of anthropogenic carbon dioxide (CO₂) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere 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 a 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, consistency within and among components, alongside methodology and data limitations. CO₂ emissions from fossil-fuel combustion and cement production (E_FF) are based on energy statistics, while emissions from land-use change (E_LUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO₂ concentration is measured directly and its rate of growth (G_ATM) is computed from the annual changes in concentration. The mean ocean CO₂ sink (S_OCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in S_OCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO₂ measurements. The global residual terrestrial CO₂ sink (S_LAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO₂ and land cover change (some including nitrogen–carbon interactions). 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 (2003–2012), E_FF was 8.6 ± 0.4 GtC yr⁻¹, E_LUC 0.9 ± 0.5 GtC yr⁻¹, G_ATM 4.3 ± 0.1 GtC yr⁻¹, S_OCEAN 2.5 ± 0.5 GtC yr⁻¹, and S_LAND 2.8 ± 0.8 GtC yr⁻¹. For year 2012 alone, E_FF grew to 9.7 ± 0.5 GtC yr⁻¹, 2.2% above 2011, reflecting a continued growing trend in these emissions, G_ATM was 5.1 ± 0.2 GtC yr⁻¹, S_OCEAN was 2.9 ± 0.5 GtC yr⁻¹, and assuming an E_LUC of 1.0 ± 0.5 GtC yr⁻¹ (based on the 2001–2010 average), S_LAND was 2.7 ± 0.9 GtC yr⁻¹. G_ATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high E_FF. The global atmospheric CO₂ concentration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that E_FF will increase by 2.1% (1.1–3.1%) to 9.9 ± 0.5 GtC in 2013, 61% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO₂ will reach about 535 ± 55 GtC for 1870–2013, about 70% from E_FF (390 ± 20 GtC) and 30% from E_LUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2013_V2.3).

Using measurements of the surface-ocean CO₂ partial pressure (pCO₂) and 14 different pCO₂ mapping methods recently collated by the Surface Ocean pCO₂ Mapping intercomparison (SOCOM) initiative, variations in regional and global sea–air CO₂ fluxes are investigated. Though the available mapping methods use widely different approaches, we find relatively consistent estimates of regional pCO₂ seasonality, in line with previous estimates. In terms of interannual variability (IAV), all mapping methods estimate the largest variations to occur in the eastern equatorial Pacific. Despite considerable spread in the detailed variations, mapping methods that fit the data more closely also tend to agree more closely with each other in regional averages. Encouragingly, this includes mapping methods belonging to complementary types – taking variability either directly from the pCO₂ data or indirectly from driver data via regression. From a weighted ensemble average, we find an IAV amplitude of the global sea–air CO₂ flux of 0.31 PgC yr⁻¹ (standard deviation over 1992–2009), which is larger than simulated by biogeochemical process models. From a decadal perspective, the global ocean CO₂ uptake is estimated to have gradually increased since about 2000, with little decadal change prior to that. The weighted mean net global ocean CO₂ sink estimated by the SOCOM ensemble is −1.75 PgC yr⁻¹ (1992–2009), consistent within uncertainties with estimates from ocean-interior carbon data or atmospheric oxygen trends.