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 methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC) 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 global net uptake of CO2 by the ocean (SOCEAN, called the ocean sink) is estimated with global ocean biogeochemistry models and observation-based fCO2 products (fCO2 is the fugacity of CO2). The global net uptake of CO2 by the land (SLAND, called the land sink) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The sum of all sources and sinks results in the carbon budget imbalance (BIM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2023, EFOS increased by 1.3 % relative to 2022, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (10.3 ± 0.5 GtC yr−1 when the cement carbonation sink is not included), and ELUC was 1.0 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1 ± 0.9 GtC yr−1 (40.6 ± 3.2 GtCO2 yr−1). Also, for 2023, GATM was 5.9 ± 0.2 GtC yr−1 (2.79 ± 0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 2.3 ± 1.0 GtC yr−1, with a near-zero BIM (−0.02 GtC yr−1). The global atmospheric CO2 concentration averaged over 2023 reached 419.31 ± 0.1 ppm. Preliminary data for 2024 suggest an increase in EFOS relative to 2023 of +0.8 % (−0.2 % to 1.7 %) globally and an atmospheric CO2 concentration increase by 2.87 ppm, reaching 422.45 ppm, 52 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2023, with a near-zero overall budget imbalance, although discrepancies of up to around 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 the following: (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) 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 mean ocean sink. This living-data update documents changes in methods and datasets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2024 (Friedlingstein et al., 2024).

The Total Carbon Column Observing Network (TCCON) measures column-average mole fractions of several greenhouse gases (GHGs), beginning in 2004, from over 30 current or past measurement sites around the world using solar absorption spectroscopy in the near-infrared (near-IR) region. TCCON GHG data have been used extensively for multiple purposes, including in studies of the carbon cycle and anthropogenic emissions, as well as to validate and improve observations from space-based sensors. Here, we describe an update to the retrieval algorithm used to process the TCCON near-IR solar spectra and to generate the associated data products. This version, called GGG2020, was initially released in April 2022. It includes updates and improvements to all steps of the retrieval, including but not limited to the conversion of the original interferograms into spectra, the spectroscopic information used in the column retrieval, post hoc air mass dependence correction, and scaling to align with the calibration scales of in situ GHG measurements. All TCCON data are available through https://tccondata.org/ (last access: 22 April 2024) and are hosted on CaltechDATA (https://data.caltech.edu/, last access: 22 April 2024). Each TCCON site has a unique DOI for its data record. An archive of all the sites' data is also available with the DOI https://doi.org/10.14291/TCCON.GGG2020 (Total Carbon Column Observing Network (TCCON) Team, 2022). The hosted files are updated approximately monthly, and TCCON sites are required to deliver data to the archive no later than 1 year after acquisition. Full details of data locations are provided in the “Code and data availability” section.

The rate of increase of atmospheric carbon dioxide (CO2) in recent decades is only about half that expected if all anthropogenic emissions had stayed airborne. The “missing” CO2 is absorbed by the oceans and land vegetation1. Of these two sinks for CO2, the ocean sink is considered to be the better constrained2-4. Nevertheless, estimates of ocean uptake, using data products that integrate surface ocean observations5-9, suggest a sink that is larger, and increasing more rapidly since about 2000, than do global ocean models, so that they now disagree substantially1. Here we examine the recent history of the CO2 sinks using atmospheric observations of CO2 and oxygen/nitrogen ratios from globally distributed measurement locations10-12. These provide a robust and independent separation of the ocean and land contributions. This method indicates that the ocean sink increased rapidly through the years 2000-2015, to values consistent with, though generally larger than, the surface-ocean based data products, but inconsistent with models. Since 2015, our results suggest this increase has stalled, but it remains 50% higher than most models predict.

We present improved estimates of air–sea CO 2 exchange over three latitude bands of the Southern Ocean using atmospheric CO 2 measurements from global airborne campaigns and an atmospheric 4-box inverse model based on a mass-indexed isentropic coordinate (M θe ). These flux estimates show two features not clearly resolved in previous estimates based on inverting surface CO 2 measurements: a weak winter-time outgassing in the polar region and a sharp phase transition of the seasonal flux cycles between polar/subpolar and subtropical regions. The estimates suggest much stronger summer-time uptake in the polar/subpolar regions than estimates derived through neural-network interpolation of pCO 2 data obtained with profiling floats but somewhat weaker uptake than a recent study by Long et al. [ Science 374 , 1275–1280 (2021)], who used the same airborne data and multiple atmospheric transport models (ATMs) to constrain surface fluxes. Our study also uses moist static energy (MSE) budgets from reanalyses to show that most ATMs tend to have excessive diabatic mixing (transport across moist isentrope, θ e , or M θe surfaces) at high southern latitudes in the austral summer, which leads to biases in estimates of air–sea CO 2 exchange. Furthermore, we show that the MSE-based constraint is consistent with an independent constraint on atmospheric mixing based on combining airborne and surface CO 2 observations.

Global climate change is influencing the seasonal cycle amplitude of atmospheric CO2 (SCA), with the strongest increases at northern high latitudes (NHL; >45° N). In this Review, we explore the changes and underlying mechanisms influencing the NHL SCA, focusing on Arctic and boreal terrestrial ecosystems. Latitudinal gradients in the SCA are largely governed by seasonality in temperature and primary production, and their influence on ecosystem carbon dynamics. In the NHL, the SCA has increased by 50% since the 1960s, mostly due to enhanced seasonality in net carbon dioxide (CO2) exchange in NHL terrestrial ecosystems. Temperature most strongly influences this trend, owing to warming impacts on growing season length and plant productivity; CO2 fertilization effects have a secondary role. Eurasian boreal ecosystems exert the strongest influence on the SCA, and spring and summer are the most influential seasons. Enhanced ecosystem respiration during the non-growing season exhibits most uncertainty in the SCA response to global and landscape drivers. Observed changes in the seasonal amplitude are projected to continue. Key priorities include extending carbon flux and ecosystem observation networks, particularly in tundra ecosystems, and including drivers such as vegetation cover and permafrost in process models to better simulate seasonal dynamics of net CO2 exchange in the NHL. Changes in the seasonal cycle amplitude of atmospheric CO2 (SCA) reflect large-scale changes in the global carbon cycle. This Review summarizes the positive SCA trend in the northern high latitudes, where the signal is strongest, and explores the underlying mechanisms driving the trend and their relative importance.

ABSTRACT Observations of radiocarbon ( 14 C) in Earth’s atmosphere and other carbon reservoirs are important to quantify exchanges of CO 2 between reservoirs. The amount of 14 C is commonly reported in the so-called Delta notation, i.e., Δ 14 C, the decay- and fractionation-corrected departure of the ratio of 14 C to total C from that ratio in an absolute international standard; this Delta notation permits direct comparison of 14 C/C ratios in the several reservoirs. However, as Δ 14 C of atmospheric CO 2 , Δ 14 CO 2 is based on the ratio of 14 CO 2 to total atmospheric CO 2 , its value can and does change not just because of change in the amount of atmospheric 14 CO 2 but also because of change in the amount of total atmospheric CO 2 , complicating ascription of change in Δ 14 CO 2 to change in one or the other quantity. Here we suggest that presentation of atmospheric 14 CO 2 amount as mole fraction relative to dry air (moles of 14 CO 2 per moles of dry air in Earth’s atmosphere), or as moles or molecules of 14 CO 2 in Earth’s atmosphere, all readily calculated from Δ 14 CO 2 and the amount of atmospheric CO 2 (with slight dependence on δ 13 CO 2 ), complements presentation only as Δ 14 CO 2 , and can provide valuable insight into the evolving budget and distribution of atmospheric 14 CO 2 .

Seasonal change of atmospheric potential oxygen (APO ∼ O 2 + CO 2 ) is a tracer for air‐sea O 2 flux with little sensitivity to the terrestrial exchange of O 2 and CO 2 . In this study, we present the tropospheric distribution and inventory of APO in each hemisphere with seasonal resolution, using O 2 and CO 2 measurements from discrete airborne campaigns between 2009 and 2018. The airborne data are represented on a mass‐weighted isentropic coordinate ( M θe ) as an alternative to latitude, which reduces the noise from synoptic variability in the APO cycles. We find a larger seasonal amplitude of APO inventory in the Southern Hemisphere relative to the Northern Hemisphere, and a larger amplitude in high latitudes (low M θe ) relative to low latitudes (high M θe ) within each hemisphere. With a box model, we invert the seasonal changes in APO inventory to yield estimates of air‐sea flux cycles at the hemispheric scale. We found a larger seasonal net outgassing of APO in the Southern Hemisphere (518 ± 52.6 Tmol) than in the Northern Hemisphere (342 ± 52.1 Tmol). Differences in APO phasing and amplitude between the hemispheres suggest distinct physical and biogeochemical mechanisms driving the air‐sea O 2 fluxes, such as fall outgassing of photosynthetic O 2 in the Northern Hemisphere, possibly associated with the formation of the seasonal subsurface shallow oxygen maximum. We compare our estimates with four model‐ and observation‐based products, identifying key limitations in these products or in the tools used to create them.

The Southern Ocean is a major sink of atmospheric CO2, 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 CO2 sink from observation-based air-sea O2 fluxes. On interannual time scales, the variability in the air-sea fluxes of CO2 and O2 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 CO2 flux estimated from observations also tends to be supported by observation-based estimates of O2 flux variability. However, the large decadal variability in air-sea CO2 flux is absent from ocean models. Our analysis suggests that issues in representing the balance between the thermal and non-thermal components of the CO2 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'.

Tropospheric helium variations are tightly linked to CO2 due to the co-emission of He and CO2 from natural-gas burning. Recently, Birner et al. (2022a) showed that the global consumption of natural gas has measurably increased the He content of the atmosphere. Like CO2, He is also predicted to exhibit complex spatial and temporal variability on shorter timescales, but measurements of these short-term variations are lacking. Here, we present the development of an improved gas delivery and purification system for the semi-continuous mass spectrometric measurement of the atmospheric He-to-nitrogen ratio (He/N2). The method replaces the chemical getter used previously by Birner et al. (2021, 2022a) to preconcentrate He in an air stream with a cryogenic trap which can be more simply regenerated by heating and which improves the precision of the measurement to 22 per meg (i.e., 0.022 ‰) in 10 min (1σ). Using this “cryo-enrichment” method, we measured the He/N2 ratios in ambient air at La Jolla (California, USA) over 5 weeks in 2022. During this period, He/N2 was strongly correlated with atmospheric CO2 concentrations, as expected from anthropogenic emissions, with a diurnal cycle of 450–500 per meg (max–min) caused by the sea–land breeze pattern of local winds, which modulates the influence of local pollution sources.

Since 1999, Environment and Climate Change Canada (ECCC) has been coordinating a multi-laboratory comparison of measurements of long-lived greenhouse gases in whole air samples collected at the Global Atmosphere Watch (GAW) Alert Observatory located in the Canadian High Arctic (82∘28′ N, 62∘30′ W). In this paper, we evaluate the measurement agreement of atmospheric CO2, CH4, N2O, SF6, and stable isotopes of CO2 (δ13C, δ18O) between leading laboratories from seven independent international institutions. The measure of success is linked to target goals for network compatibility outlined by the World Meteorological Organization's (WMO) GAW greenhouse gas measurement community. Overall, based on ∼ 8000 discrete flask samples, we find that the co-located atmospheric CO2 and CH4 measurement records from Alert by CSIRO, MPI-BGC, SIO, UHEI-IUP, and ECCC versus NOAA (the designated reference laboratory) are generally consistent with the WMO compatibility goals of ± 0.1 ppm CO2 and ± 2 ppb CH4 over the 17-year period (1999–2016), although there are periods where differences exceed target levels and persist as systematic bias for months or years. Consistency with the WMO goals for N2O, SF6, and stable isotopes of CO2 (δ13C, δ18O) has not been demonstrated. Additional analysis of co-located comparison measurements between CSIRO and SIO versus NOAA or INSTAAR (for the isotopes of CO2) at other geographical sites suggests that the findings at Alert for CO2, CH4, N2O, and δ13C–CO2 could be extended across the CSIRO, SIO, and NOAA observing networks. The primary approach to estimate an overall measurement agreement level was carried out by pooling the differences of all individual laboratories versus the designated reference laboratory and determining the 95th percentile range of these data points. Using this approach over the entire data record, our best estimate of the measurement agreement range is −0.51 to +0.53 ppm for CO2, −0.09 ‰ to +0.07 ‰ for δ13C, −0.50 ‰ to +0.58 ‰ for δ18O, −4.86 to +6.16 ppb for CH4, −0.75 to +1.20 ppb for N2O, and −0.14 to +0.09 ppt for SF6. A secondary approach of using the average of 2 standard deviations of the means for all flask samples taken in each individual sampling episode provided similar results. These upper and lower limits represent our best estimate of the measurement agreement at the 95 % confidence level for these individual laboratories, providing more confidence for using these datasets in various scientific applications (e.g., long-term trend analysis).

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 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 CO2 concentration is measured directly, and its growth rate (G(ATM)) is computed from the annual changes in concentration. The ocean CO2 sink (S-OCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO(2) products. The terrestrial CO2 sink (S-LAND) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system 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 incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as +/- 1 sigma. For the year 2022, E-FOS increased by 0.9% relative to 2021, with fossil emissions at 9.9 +/- 0.5 GtC yr(-1) (10.2 +/- 0.5 GtC yr(-1) when the cement carbonation sink is not included), and E-LUC was 1.2 +/- 0.7 GtC yr(-1), for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1 +/- 0.8 GtC yr(-1) (40.7 +/- 3.2 GtCO(2) yr(-1)). Also, for 2022, G(ATM) was 4.6 +/- 0.2 GtC yr(-1) (2.18 +/- 0.1 ppm yr(-1); ppm denotes parts per million), S-OCEAN was 2.8 +/- 0.4 GtC yr(-1), and S-LAND was 3.8 +/- 0.8 GtC yr(-1), with a B-IM of 0.1 GtC yr(-1) (i.e. total estimated sources marginally too low or sinks marginally too high). The global atmospheric CO2 concentration averaged over 2022 reached 417.1 +/- 0.1 ppm. Preliminary data for 2023 suggest an increase in E-FOS relative to 2022 of +/- 1:1% (0.0% to 2.1 %) globally and atmospheric CO2 concentration reaching 419.3 ppm, 51% above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959-2022, with a near-zero overall budget imbalance, although discrepancies of up to around 1 Gt Cyr(-1) persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (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 methods and data sets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2023 (Friedlingstein et al., 2023).

The interannual to decadal variability in natural carbon sinks limits the explanation of recent changes in atmospheric CO 2 concentration. Here we account for interannual and decadal variability using a simple quasi-mechanistic model of the net land carbon exchange with terms scaling with atmospheric CO 2 and a weighted spatial average of temperature anomalies. This approach reduces the unexplained residual in Earth’s carbon cycle budget from ±0.76 GtC per year obtained using process models to ±0.50 GtC per year, with the largest improvements on decadal timescales despite assuming constant dynamics. Our findings reveal remarkable stability of the carbon cycle and allow verification of reported global emissions to within 4.4% (95% confidence level) over the five-year stocktake cycle of the Paris Agreement—half the uncertainty reported previously.

Understanding terrestrial ecosystems and their response to anthropogenic climate change requires quantification of land-atmosphere carbon exchange. However, top-down and bottom-up estimates of large-scale land-atmosphere fluxes, including the northern extratropical growing season net flux (GSNF), show significant discrepancies. We developed a data-driven metric for the GSNF using atmospheric carbon dioxide concentration observations collected during the High-Performance Instrumented Airborne Platform for Environmental Research Pole-to-Pole Observations and Atmospheric Tomography Mission flight campaigns. This aircraft-derived metric is bias-corrected using three independent atmospheric inversion systems. We estimate the northern extratropical GSNF to be 5.7 +/- 0.3 Pg C and use it to evaluate net biosphere productivity from the Coupled Model Intercomparison Project phase 5 and 6 (CMIP5 and CMIP6) models. While the model-to-model spread in the GSNF has decreased in the CMIP6 models relative to that of the CMIP5 models, there is still disagreement on the magnitude and timing of seasonal carbon uptake with most models underestimating the GSNF and overestimating the length of the growing season relative to the observations. We also use an emergent constraint approach to estimate annual northern extratropical gross primary productivity to be 56 +/- 17 Pg C, heterotrophic respiration to be 25 +/- 13 Pg C, and net primary productivity to be 28 +/- 12 Pg C. The flux inferred from these aircraft observations provides an additional constraint on large-scale gross fluxes in prognostic Earth system models that may ultimately improve our ability to accurately predict carbon-climate feedbacks. Plain Language Summary The exchange of carbon between the land and atmosphere is an important part of the Earth's climate, and this exchange might change due to human-caused climate change. However, estimates of land-atmosphere carbon fluxes made using different techniques do not agree with each other. We use atmospheric carbon dioxide observations collected during two flight campaigns to show that 5.7 Pg C is exchanged between the atmosphere and the land in the northern hemisphere during the summer growing season. This estimate is used to evaluate the performance of two generations of climate prediction models. The newer generation of models show less spread than the older generation, but there is still significant disagreement on the magnitude and timing of land-atmosphere carbon exchange among models. Most models underestimate the growing season net flux and overestimate the length of the growing season. We also use our observational estimate to reduce the spread on component fluxes of carbon exchange, namely uptake by photosynthesis and release by respiration.

The Los Angeles (LA) basin was responsible for approximately 20% of California’s methane emissions in 2016. Hence, curtailment of these emissions is required to meet California’s greenhouse gas emissions reduction targets. However, effective mitigation remains challenging in the presence of diverse methane sources like oil and gas production fields, refineries, landfills, wastewater treatment facilities, and natural gas infrastructure. In this study, we study the temporal variability in the surface concentrations from February 2015 to April 2022 to detect a declining trend in methane emissions. We quantify the reduction due to this declining trend through inverse modeling and show that methane emissions in the LA basin have declined by 15 Gg, or ∼7% over five years from January 2015 to May 2020.

A central challenge for sustaining international cooperation to cut global greenhouse gas emissions is confidence that national policy efforts are leading to a meaningful impact on the climate. Here, we apply a detection protocol to determine when the measurable signal of atmospheric CO _2 can be distinguished from the noise of the carbon cycle and uncertainties in emission trends. We test that protocol with a database of 226 emission mitigation scenarios—the universe of scenarios vetted by the Intergovernmental Panel on Climate Change. These scenarios are descriptive of ‘baseline’ trajectories of emissions trends in the absence of new policies along with trajectories that reflect substantial policy efforts to stop warming at 1.5 °C–2 °C above pre-industrial levels, as embodied in the Paris Agreement. The most aggressive mitigation scenarios (i.e. 1.5 °C) require 11–16 years to detect a signal of demonstrable progress from the noise; 2 °C scenarios lengthen detection by at least a decade. As more climate policy analysts face the reality that goals of 1.5 °C–2 °C seem infeasible, they have developed ‘overshoot’ scenarios with emissions that rise above the agreed goal and then, later on, fall aggressively to achieve it. These pathways come at the political cost of a 1–2 decade delay in detection, even for the 1.5 °C scenarios. The Paris Agreement requires a global ‘stocktake’ that interrogates national mitigation efforts; our results suggest that this effort must grapple with the question of when the world can gain confidence that the diplomacy on climate is demonstrably making an impact.

Atmospheric measurements of the O2/N2 ratio and the CO2 mole fraction (combined into the conceptual tracer "Atmospheric Potential Oxygen", APO) over continents have been proposed as a constraint on CO2 emissions from fossil-fuel burning. Here we assess the suitability of such APO data to constrain anthropogenic CO2 emissions in Western Europe, with particular focus on their decadal trends. We use an inversion of atmospheric transport to estimate spatially and temporally explicit scaling factors on a bottom-up fossil-fuel emissions inventory. Based on the small number of currently available observational records, our CO2 emissions estimates show relatively large apparent year-to-year variations, exceeding the expected uncertainty of the bottom-up inventory and precluding the calculation of statistically significant trends. We were not able to trace the apparent year-to-year variations back to particular properties of the APO data. Inversion of synthetic APO data, however, confirms that data information content and degrees of freedom are sufficient to successfully correct a counterfactual prior. Larger sets of measurement stations, such as the recently started APO observations from the Integrated Carbon Observation System (ICOS) European research infrastructure, improve the constraint and may ameliorate possible problems with local signals or with measurement or model errors at the stations. We further tested the impact of uncertainties in the O2:CO2 stoichiometries of fossil-fuel burning and land biospheric exchange and found they are not fundamental obstacles to estimating decadal trends in fossil-fuel CO2 emissions, though further work on fossil-fuel O2:CO2 stoichiometries seems necessary.

Concurrent measurements of atmospheric O-2 and CO2 amount fractions have been used for decades now to estimate fluxes of carbon to and from the oceans and the land biosphere. The equations used in these estimates explicitly include fossil fuel combustion but are built on the assumption that fossil fuels are oxidized solely by atmospheric O-2, ignoring the small fraction of fossil fuel oxidation by oxide ores that occurs during refining of metals, and thereby overestimating losses of atmospheric O-2. Here, we address this deficiency by quantifying the effective O-2 fluxes associated with the processing of iron, aluminium and copper. We also consider the potential impact of sulfur. We find that consideration of the oxygen mobilized during metal processing (equivalent to a net O-2 flux of 12.0(-0.4)(+0.2) Tmol a(-1)) leads to an increased estimate of ocean carbon uptake in the years 2000-2010 of (0.144(-0.005)(+0.002))Pg a(-1) of carbon with a corresponding decrease in estimated land uptake. A rough estimate of oxygen uptake due to sulfur chemistry during fossil fuel combustion (2.4 Tmol a(-1)) decreases ocean carbon uptake but by a much smaller amount. These corrections are small compared to existing estimates of the fluxes and their uncertainties ((2.27 +/- 0.60) and (1.05 +/- 0.84) Pg a(-1) of carbon for ocean and land respectively (Keeling and Manning, 2014)) but should be employed in future analyses.

Abstract. This study considers year-to-year and decadal variations as well as secular trends of the sea–air CO2 flux over the 1957–2020 period, as constrained by the pCO2 measurements from the SOCAT data base. In a first step, we relate interannual anomalies in ocean-internal carbon sources and sinks to local interannual anomalies in sea surface temperature (SST), the temporal changes of SST (dSST/dt), and squared wind speed (u2), employing a multi-linear regression. In the tropical Pacific, we find interannual variability to be dominated by dSST/dt, as arising from variations in the upwelling of colder and more carbon-rich waters into the mixed layer. In the eastern upwelling zones as well as in circumpolar bands in the high latitudes of both hemispheres, we find sensitivity to wind speed, compatible with the entrainment of carbon-rich water during wind-driven deepening of the mixed layer and wind-driven upwelling. In the Southern Ocean, the secular increase in wind speed leads to a secular increase in the carbon source into the mixed layer, with an estimated reduction of the sink trend in the range 17 to 42 %. In a second step, we combined the result of the multi-linear regression and an explicitly interannual pCO2-based additive correction into a “hybrid” estimate of the sea–air CO2 flux over the period 1957–2020. As a pCO2 mapping method, it combines (a) the ability of a regression to bridge data gaps and extrapolate into the early decades almost void of pCO2 data based on process-related observables and (b) the ability of an autoregressive interpolation to follow signals even if not represented in the chosen set of explanatory variables. The “hybrid” estimate can be applied as ocean flux prior for atmospheric CO2 inversions covering the whole period of atmospheric CO2 data since 1957.