Improved atmospheric constraints on Southern Ocean CO2 exchange.

We present improved estimates of air-sea CO2 exchange over three latitude bands of the Southern Ocean using atmospheric CO2 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 CO2 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 pCO2 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 CO2 exchange. Furthermore, we show that the MSE-based constraint is consistent with an independent constraint on atmospheric mixing based on combining airborne and surface CO2 observations.

Estimating regional fossil fuel CO2 concentrations from 14CO2 observations: challenges and uncertainties.

The direct way to estimate the regional fossil fuel CO2 surplus (ΔffCO2) at a station is by measuring the Δ14CO2 depletion compared with a respective background. However, this approach has several challenges, which are (i) the choice of an appropriate Δ14CO2 background, (ii) potential contaminations through nuclear 14CO2 emissions and (iii) masking of ΔffCO2 by 14C-enriched biosphere respiration. Here we evaluate these challenges and estimate potential biases and typical uncertainties of 14C-based ΔffCO2 estimates in Europe. We show that Mace Head (MHD), Ireland, is a representative background station for the Integrated Carbon Observation System (ICOS) atmosphere station network. The mean ΔffCO2 representativeness bias when using the MHD Δ14CO2 background for the whole observation network is of order 0.1 ± 0.3 ppm. At ICOS sites, the median nuclear contamination leads to 25% low-biased ΔffCO2 estimates if not corrected for. The ΔffCO2 masking due to 14C-enriched heterotrophic CO2 respiration can lead to similar ΔffCO2 biases as the nuclear contaminations, especially in summer. Our evaluation of all components contributing to the uncertainty of ΔffCO2 estimates reveals that, due to the small ffCO2 signals at ICOS stations, almost half of the 14C-based ΔffCO2 estimates from integrated samples have an uncertainty that is larger than 50%. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Sparse observations induce large biases in estimates of the global ocean CO2 sink: an ocean model subsampling experiment.

Estimates of ocean [Formula: see text] uptake from global ocean biogeochemistry models and [Formula: see text]-based data products differ substantially, especially in high latitudes and in the trend of the [Formula: see text] uptake since 2000. Here, we assess the effect of data sparsity on two [Formula: see text]-based estimates by subsampling output from a global ocean biogeochemistry model. The estimates of the ocean [Formula: see text] uptake are improved from a sampling scheme that mimics present-day sampling to an ideal sampling scheme with 1000 evenly distributed sites. In particular, insufficient sampling has given rise to strong biases in the trend of the ocean carbon sink in the [Formula: see text] products. The overestimation of the [Formula: see text] flux trend by 20-35% globally and 50-130% in the Southern Ocean with the present-day sampling is reduced to less than [Formula: see text] with the ideal sampling scheme. A substantial overestimation of the decadal variability of the Southern Ocean carbon sink occurs in one product and appears related to a skewed data distribution in [Formula: see text] space. With the ideal sampling, the bias in the mean [Formula: see text] flux is reduced from 9-12% to 2-9% globally and from 14-26% to 5-17% in the Southern Ocean. On top of that, discrepancies of about [Formula: see text] (15%) persist due to uncertainties in the gas-exchange calculation. 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'.

Diagnosing destabilization risk in global land carbon sinks.

Global net land carbon uptake or net biome production (NBP) has increased during recent decades1. Whether its temporal variability and autocorrelation have changed during this period, however, remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink2,3. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO2 concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.

Impact of Changing Winds on the Mauna Loa CO2 Seasonal Cycle in Relation to the Pacific Decadal Oscillation.

Long-term measurements at the Mauna Loa Observatory (MLO) show that the CO2 seasonal cycle amplitude (SCA) increased from 1959 to 2019 at an overall rate of 0.22  ±  0.034 ppm decade-1 while also varying on interannual to decadal time scales. These SCA changes are a signature of changes in land ecological CO2 fluxes as well as shifting winds. Simulations with the TM3 tracer transport model and CO2 fluxes from the Jena CarboScope CO2 Inversion suggest that shifting winds alone have contributed to a decrease in SCA of -0.10  ±  0.022 ppm decade-1 from 1959 to 2019, partly offsetting the observed long-term SCA increase associated with enhanced ecosystem net primary production. According to these simulations and MIROC-ACTM simulations, the shorter-term variability of MLO SCA is nearly equally driven by varying ecological CO2 fluxes (49%) and varying winds (51%). We also show that the MLO SCA is strongly correlated with the Pacific Decadal Oscillation (PDO) due to varying winds, as well as with a closely related wind index (U-PDO). Since 1980, 44% of the wind-driven SCA decrease has been tied to a secular trend in the U-PDO, which is associated with a progressive weakening of westerly winds at 700 mbar over the central Pacific from 20°N to 40°N. Similar impacts of varying winds on the SCA are seen in simulations at other low-latitude Pacific stations, illustrating the difficulty of constraining trend and variability of land CO2 fluxes using observations from low latitudes due to the complexity of circulation changes.

Regional and seasonal partitioning of water and temperature controls on global land carbon uptake variability.

Global fluctuations in annual land carbon uptake (NEEIAV) depend on water and temperature variability, yet debate remains about local and seasonal controls of the global dependences. Here, we quantify regional and seasonal contributions to the correlations of globally-averaged NEEIAV against terrestrial water storage (TWS) and temperature, and respective uncertainties, using three approaches: atmospheric inversions, process-based vegetation models, and data-driven models. The three approaches agree that the tropics contribute over 63% of the global correlations, but differ on the dominant driver of the global NEEIAV, because they disagree on seasonal temperature effects in the Northern Hemisphere (NH, >25°N). In the NH, inversions and process-based models show inter-seasonal compensation of temperature effects, inducing a global TWS dominance supported by observations. Data-driven models show weaker seasonal compensation, thereby estimating a global temperature dominance. We provide a roadmap to fully understand drivers of global NEEIAV and discuss their implications for future carbon-climate feedbacks.

Strong Southern Ocean carbon uptake evident in airborne observations.

The Southern Ocean plays an important role in determining atmospheric carbon dioxide (CO2), yet estimates of air-sea CO2 flux for the region diverge widely. In this study, we constrained Southern Ocean air-sea CO2 exchange by relating fluxes to horizontal and vertical CO2 gradients in atmospheric transport models and applying atmospheric observations of these gradients to estimate fluxes. Aircraft-based measurements of the vertical atmospheric CO2 gradient provide robust flux constraints. We found an annual mean flux of ­0.53 ± 0.23 petagrams of carbon per year (net uptake) south of 45°S during the period 2009­2018. This is consistent with the mean of atmospheric inversion estimates and surface-ocean partial pressure of CO2 (Pco2)­based products, but our data indicate stronger annual mean uptake than suggested by recent interpretations of profiling float observations.

Decadal variability in land carbon sink efficiency.

BACKGROUND: The climate mitigation target of limiting the temperature increase below 2 °C above the pre-industrial levels requires the efforts from all countries. Tracking the trajectory of the land carbon sink efficiency is thus crucial to evaluate the nationally determined contributions (NDCs). Here, we define the instantaneous land sink efficiency as the ratio of natural land carbon sinks to emissions from fossil fuel and land-use and land-cover change with a value of 1 indicating carbon neutrality to track its temporal dynamics in the past decades. RESULTS: Land sink efficiency has been decreasing during 1957-1990 because of the increased emissions from fossil fuel. After the effect of the Mt. Pinatubo eruption diminished (after 1994), the land sink efficiency firstly increased before 2009 and then began to decrease again after 2009. This reversal around 2009 is mostly attributed to changes in land sinks in tropical regions in response to climate variations. CONCLUSIONS: The decreasing trend of land sink efficiency in recent years reveals greater challenges in climate change mitigation, and that climate impacts on land carbon sinks must be accurately quantified to assess the effectiveness of regional scale climate mitigation policies.

State of the science in reconciling top-down and bottom-up approaches for terrestrial CO2 budget.

Robust estimates of CO2 budget, CO2 exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic CO2 emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: "top-down" atmospheric inversions and "bottom-up" biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere-biosphere CO2 exchange. Here, we conducted a thorough comparison of CO2 budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO2 budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO2 budgets in the 2000s. Regionally, improved agreement in CO2 budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO2 budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink-source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top-down and bottom-up approaches and yield more robust knowledge of regional CO2 budgets.

Global atmospheric CO2 inverse models converging on neutral tropical land exchange, but disagreeing on fossil fuel and atmospheric growth rate.

We have compared a suite of recent global CO2 atmospheric inversion results to independent airborne observations and to each other, to assess their dependence on differences in northern extratropical (NET) vertical transport and to identify some of the drivers of model spread. We evaluate posterior CO2 concentration profiles against observations from the High-Performance Instrumented Airborne Platform for Environmental Research (HIAPER) Pole-to-Pole Observations (HIPPO) aircraft campaigns over the mid-Pacific in 2009-2011. Although the models differ in inverse approaches, assimilated observations, prior fluxes, and transport models, their broad latitudinal separation of land fluxes has converged significantly since the Atmospheric Carbon Cycle Inversion Intercomparison (TransCom 3) and the REgional Carbon Cycle Assessment and Processes (RECCAP) projects, with model spread reduced by 80% since TransCom 3 and 70% since RECCAP. Most modeled CO2 fields agree reasonably well with the HIPPO observations, specifically for the annual mean vertical gradients in the Northern Hemisphere. Northern Hemisphere vertical mixing no longer appears to be a dominant driver of northern versus tropical (T) annual flux differences. Our newer suite of models still gives northern extratropical land uptake that is modest relative to previous estimates (Gurney et al., 2002; Peylin et al., 2013) and near-neutral tropical land uptake for 2009-2011. Given estimates of emissions from deforestation, this implies a continued uptake in intact tropical forests that is strong relative to historical estimates (Gurney et al., 2002; Peylin et al., 2013). The results from these models for other time periods (2004-2014, 2001-2004, 1992-1996) and reevaluation of the TransCom 3 Level 2 and RECCAP results confirm that tropical land carbon fluxes including deforestation have been near neutral for several decades. However, models still have large disagreements on ocean-land partitioning. The fossil fuel (FF) and the atmospheric growth rate terms have been thought to be the best-known terms in the global carbon budget, but we show that they currently limit our ability to assess regional-scale terrestrial fluxes and ocean-land partitioning from the model ensemble.

Reassessing Southern Ocean Air-Sea CO2 Flux Estimates With the Addition of Biogeochemical Float Observations.

New estimates of pCO2 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 CO2 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 pCO2 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.

Impact of the 2015/2016 El Niño on the terrestrial carbon cycle constrained by bottom-up and top-down approaches.

Evaluating the response of the land carbon sink to the anomalies in temperature and drought imposed by El Niño events provides insights into the present-day carbon cycle and its climate-driven variability. It is also a necessary step to build confidence in terrestrial ecosystems models' response to the warming and drying stresses expected in the future over many continents, and particularly in the tropics. Here we present an in-depth analysis of the response of the terrestrial carbon cycle to the 2015/2016 El Niño that imposed extreme warming and dry conditions in the tropics and other sensitive regions. First, we provide a synthesis of the spatio-temporal evolution of anomalies in net land-atmosphere CO2 fluxes estimated by two in situ measurements based on atmospheric inversions and 16 land-surface models (LSMs) from TRENDYv6. Simulated changes in ecosystem productivity, decomposition rates and fire emissions are also investigated. Inversions and LSMs generally agree on the decrease and subsequent recovery of the land sink in response to the onset, peak and demise of El Niño conditions and point to the decreased strength of the land carbon sink: by 0.4-0.7 PgC yr-1 (inversions) and by 1.0 PgC yr-1 (LSMs) during 2015/2016. LSM simulations indicate that a decrease in productivity, rather than increase in respiration, dominated the net biome productivity anomalies in response to ENSO throughout the tropics, mainly associated with prolonged drought conditions.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'.

Spring photosynthetic onset and net CO2 uptake in Alaska triggered by landscape thawing.

The springtime transition to regional-scale onset of photosynthesis and net ecosystem carbon uptake in boreal and tundra ecosystems are linked to the soil freeze-thaw state. We present evidence from diagnostic and inversion models constrained by satellite fluorescence and airborne CO2 from 2012 to 2014 indicating the timing and magnitude of spring carbon uptake in Alaska correlates with landscape thaw and ecoregion. Landscape thaw in boreal forests typically occurs in late April (DOY 111 ± 7) with a 29 ± 6 day lag until photosynthetic onset. North Slope tundra thaws 3 weeks later (DOY 133 ± 5) but experiences only a 20 ± 5 day lag until photosynthetic onset. These time lag differences reflect efficient cold season adaptation in tundra shrub and the longer dehardening period for boreal evergreens. Despite the short transition from thaw to photosynthetic onset in tundra, synchrony of tundra respiration with snow melt and landscape thaw delays the transition from net carbon loss (at photosynthetic onset) to net uptake by 13 ± 7 days, thus reducing the tundra net carbon uptake period. Two global CO2 inversions using a CASA-GFED model prior estimate earlier northern high latitude net carbon uptake compared to our regional inversion, which we attribute to (i) early photosynthetic-onset model prior bias, (ii) inverse method (scaling factor + optimization window), and (iii) sparsity of available Alaskan CO2 observations. Another global inversion with zero prior estimates the same timing for net carbon uptake as the regional model but smaller seasonal amplitude. The analysis of Alaskan eddy covariance observations confirms regional scale findings for tundra, but indicates that photosynthesis and net carbon uptake occur up to 1 month earlier in evergreens than captured by models or CO2 inversions, with better correlation to above-freezing air temperature than date of primary thaw. Further collection and analysis of boreal evergreen species over multiple years and at additional subarctic flux towers are critically needed.

Land use change and El Niño-Southern Oscillation drive decadal carbon balance shifts in Southeast Asia.

An integrated understanding of the biogeochemical consequences of climate extremes and land use changes is needed to constrain land-surface feedbacks to atmospheric CO2 from associated climate change. Past assessments of the global carbon balance have shown particularly high uncertainty in Southeast Asia. Here, we use a combination of model ensembles to show that intensified land use change made Southeast Asia a strong source of CO2 from the 1980s to 1990s, whereas the region was close to carbon neutral in the 2000s due to an enhanced CO2 fertilization effect and absence of moderate-to-strong El Niño events. Our findings suggest that despite ongoing deforestation, CO2 emissions were substantially decreased during the 2000s, largely owing to milder climate that restores photosynthetic capacity and suppresses peat and deforestation fire emissions. The occurrence of strong El Niño events after 2009 suggests that the region has returned to conditions of increased vulnerability of carbon stocks.

Compensatory water effects link yearly global land CO2 sink changes to temperature.

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO2) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

European land CO2 sink influenced by NAO and East-Atlantic Pattern coupling.

Large-scale climate patterns control variability in the global carbon sink. In Europe, the North-Atlantic Oscillation (NAO) influences vegetation activity, however the East-Atlantic (EA) pattern is known to modulate NAO strength and location. Using observation-driven and modelled data sets, we show that multi-annual variability patterns of European Net Biome Productivity (NBP) are linked to anomalies in heat and water transport controlled by the NAO-EA interplay. Enhanced NBP occurs when NAO and EA are both in negative phase, associated with cool summers with wet soils which enhance photosynthesis. During anti-phase periods, NBP is reduced through distinct impacts of climate anomalies in photosynthesis and respiration. The predominance of anti-phase years in the early 2000s may explain the European-wide reduction of carbon uptake during this period, reported in previous studies. Results show that improving the capability of simulating atmospheric circulation patterns may better constrain regional carbon sink variability in coupled carbon-climate models.

Enhanced seasonal CO2 exchange caused by amplified plant productivity in northern ecosystems.

Atmospheric monitoring of high northern latitudes (above 40°N) has shown an enhanced seasonal cycle of carbon dioxide (CO2) since the 1960s, but the underlying mechanisms are not yet fully understood. The much stronger increase in high latitudes relative to low ones suggests that northern ecosystems are experiencing large changes in vegetation and carbon cycle dynamics. We found that the latitudinal gradient of the increasing CO2 amplitude is mainly driven by positive trends in photosynthetic carbon uptake caused by recent climate change and mediated by changing vegetation cover in northern ecosystems. Our results underscore the importance of climate-vegetation-carbon cycle feedbacks at high latitudes; moreover, they indicate that in recent decades, photosynthetic carbon uptake has reacted much more strongly to warming than have carbon release processes.

The reinvigoration of the Southern Ocean carbon sink.

Several studies have suggested that the carbon sink in the Southern Ocean-the ocean's strongest region for the uptake of anthropogenic CO2 -has weakened in recent decades. We demonstrated, on the basis of multidecadal analyses of surface ocean CO2 observations, that this weakening trend stopped around 2002, and by 2012, the Southern Ocean had regained its expected strength based on the growth of atmospheric CO2. All three Southern Ocean sectors have contributed to this reinvigoration of the carbon sink, yet differences in the processes between sectors exist, related to a tendency toward a zonally more asymmetric atmospheric circulation. The large decadal variations in the Southern Ocean carbon sink suggest a rather dynamic ocean carbon cycle that varies more in time than previously recognized.

Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate.

Terrestrial gross primary production (GPP) is the largest global CO(2) flux driving several ecosystem functions. We provide an observation-based estimate of this flux at 123 +/- 8 petagrams of carbon per year (Pg C year(-1)) using eddy covariance flux data and various diagnostic models. Tropical forests and savannahs account for 60%. GPP over 40% of the vegetated land is associated with precipitation. State-of-the-art process-oriented biosphere models used for climate predictions exhibit a large between-model variation of GPP's latitudinal patterns and show higher spatial correlations between GPP and precipitation, suggesting the existence of missing processes or feedback mechanisms which attenuate the vegetation response to climate. Our estimates of spatially distributed GPP and its covariation with climate can help improve coupled climate-carbon cycle process models.