the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Interannual global carbon cycle variations linked to atmospheric circulation variability
Na Li
Sebastian Sippel
Alexander J. Winkler
Miguel D. Mahecha
Markus Reichstein
Ana Bastos
One of the least understood temporal scales of global carbon cycle (C-cycle) dynamics is its interannual variability (IAV). This variability is mainly driven by variations in the local climatic drivers of terrestrial ecosystem activity, which in turn are controlled by large-scale modes of atmospheric variability. Here, we quantify the fraction of global C-cycle IAV that is explained by large-scale atmospheric circulation variability, which is quantified by spatiotemporal sea level pressure (SLP) fields. C-cycle variability is diagnosed from the global detrended atmospheric CO2 growth rate and the land CO2 sink from 16 dynamic global vegetation models and two atmospheric inversions in the Global Carbon Budget 2018. We use a regularized linear regression model, which represents a statistical learning technique apt to deal with the large number of atmospheric circulation predictors (p≥800, each representing one pixel-based time series of SLP anomalies) in a relatively short observed record (n<60 years). We show that boreal winter and spring SLP anomalies allow predicting IAV in the atmospheric CO2 growth rate and the global land sink, with Pearson correlations between reference and predicted values between 0.70 and 0.84 for boreal winter SLP anomalies. This is comparable to or higher than that of a similar model using 15 traditional teleconnection indices as predictors. The spatial patterns of regression coefficients of the model based on SLP fields show a predominant role of the tropical Pacific and over Southeast Asia extending to Australia, corresponding to the regions associated with the El Niño–Southern Oscillation variability. We also identify another important region in the western Pacific, roughly corresponding to the West Pacific pattern.
We further evaluate the influence of the time series length on the predictability of IAV and find that reliable estimates of global C-cycle IAV can be obtained from records of 30–54 years. For shorter time series (n<30 years), however, our results show that conclusions about CO2 IAV patterns and drivers need to be evaluated with caution. Overall, our study illustrates a new data-driven and flexible approach to model the relationship between large-scale atmospheric circulation variations and C-cycle variability at global and regional scales, complementing the traditional use of teleconnection indices.
- Article
(14692 KB) - Full-text XML
- BibTeX
- EndNote
The global carbon cycle (C-cycle) varies at multiple timescales ranging from minutes to a millennium (Ciais et al., 2013). Quantifying and understanding the patterns of variability in the C-cycle and their drivers are crucial to better understanding C-cycle dynamics and better constraining future climate projections (Cox et al., 2013; Friedlingstein et al., 2014). Primarily driven by the land sink dynamics (Piao et al., 2020), interannual variability (IAV) is one of the most uncertain and poorly understood terms of the global C-cycle in the observational period (Friedlingstein et al., 2019).
A fundamental challenge is that variability in land–atmosphere carbon exchange is affected in complex ways by large-scale atmosphere circulation modes but also land use change, forced climate changes, and direct physiological CO2 effects on ecosystems, among others (IPCC, 2013). Separating these effects is difficult because of the large uncertainties associated with some processes, for example land use change (Friedlingstein et al., 2022) or processes not directly observable at global scale (e.g., photosynthesis or heterotrophic respiration) (Schimel et al., 2015; Basile et al., 2020). The second challenge is that the land sink, as a balance of carbon uptake and release, responds differently to variations in the climatic forcing (Jung et al., 2017; Piao et al., 2020). This makes it hard to attribute induced land sink IAV to specific drivers, which is crucial for process understanding (Jung et al., 2017; Humphrey et al., 2018, 2021; Wang et al., 2022). Last, the limited length of observational records may hamper robust statistical analysis (IPCC, 2013): the longest continuous observations of atmospheric CO2 at the South Pole and Mauna Loa observatory only exist from 1958 onwards (Dlugokencky and Tans, 2019).
Land biospheric CO2 uptake results from the net balance of carbon uptake from gross primary productivity, release from multiple respiration terms, and disturbance-induced fluxes such as fires, amongst other smaller terms (IPCC, 2013). Most of these processes are primarily driven by temperature, water, and radiation availability (Jung et al., 2017). These meteorological drivers are, in turn, modulated by large-scale modes of atmospheric circulation on multiple timescales, such as the El Niño–Southern Oscillation (ENSO; Gu and Adler, 2011) and the Pacific Decadal Oscillation (PDO; Newman et al., 2016). These climate variability modes are generated within coupled atmosphere–ocean systems (Ghil, 2002) and considered irreducible noise in climate projections (Madden, 1976; Schneider and Kinter, 1994; Deser et al., 2012). Because these modes typically interact and affect weather dynamics in regions beyond those where they emerge, such modes are collectively referred to as teleconnections (IPCC, 2013). Bacastow (1976) showed that ENSO is highly correlated with annual variations in observed atmospheric CO2 measured at the South Pole and Mauna Loa, Hawaii. Keeling et al. (1995) attributed these correlations to the ENSO impact on the biospheric sink. In addition to ENSO, Zhu et al. (2017) showed that the PDO and the Atlantic Multidecadal Oscillation (AMO; Enfield et al., 2001; Rayner et al., 2003) may also influence global terrestrial ecosystem carbon fluxes and that other modes of variability in the Northern Hemisphere also have local impacts on carbon cycling (Zhu et al., 2017).
A common approach to diagnose the impacts of natural climate variability is to use ensembles of Earth system model simulations with perturbed initial conditions to quantify the impacts of natural climate variability at decadal to millennium scales (Frölicher et al., 2013). However, the inherently chaotic atmosphere, in combination with model structural uncertainty, implies large uncertainties for future projections (Deser et al., 2020). In addition, Earth system model projections can be compromised by limited representation of the full complexity of physical processes involved, lack of observational constraints, and high computational demands when aiming to resolve high resolutions (Randall et al., 2007; IPCC, 2013).
Statistical approaches are a simplified but effective way to reveal physical processes in observations (von Storch, 1995). A traditional approach consists of evaluating relationships between the variables of interest (e.g., CO2 time series) and teleconnection indices (Bacastow, 1976; Bastos et al., 2013; Zhu et al., 2017). As a simple representation of the large-scale atmospheric circulation modes, teleconnection indices are extracted mainly from sea surface temperature or atmospheric anomalies (Kumar and Hoerling, 1997; IPCC, 2013). Such indices are an effective way to reduce the complexity of the spatiotemporal variability in multiple variables (Stenseth et al., 2003; Wills et al., 2017), but they may not be able to capture spatial variations in the large-scale atmospheric circulation modes themselves.
Recently, Sippel et al. (2019) applied ridge regression, a regularized linear regression method (Hastie et al., 2009; Friedman et al., 2010), to quantify the component of precipitation and temperature variability driven by atmospheric variations based on sea level pressure (SLP) fields rather than teleconnection indices. Their approach allowed them to robustly infer the main spatiotemporal patterns of atmospheric variability influencing these two climate variables. On the one hand, including a field of circulation-based predictors, avoids considering predefined assumptions about their spatial configurations as they are common to teleconnection indices, while compensating for relatively short historical records. The regularization approach, on the other hand, allows overcoming overfitting and multicollinearity issues due to short time series and a very large number of spatial predictors.
In this study, we adopt the ridge regression approach in Sippel et al. (2019), aiming to quantify the fraction of global C-cycle IAV influenced by large-scale atmospheric circulation variability. We use observation-based time series of global atmospheric CO2 growth rate (AGR) and land CO2 surface fluxes from atmospheric inversions and dynamic global vegetation models (DGVMs), as well as the land sink estimated as a residual of other terms in the Global Carbon Budget 2018 (Le Quéré et al., 2018). We first evaluate and compare the predictive skill of the ridge regression model when using SLP as predictors versus commonly used traditional teleconnection indices (Sect. 3.1). Next, we analyze and discuss the global C-cycle sensitivity to atmospheric circulation variability from various latitudinal domains of SLP anomaly fields (Sect. 3.2). Finally, we evaluate the sensitivity of the results to the length of the time series (Sect. 3.3) by comparing the fraction of C-cycle IAV that can be explained by large-scale atmospheric circulation variability based on these datasets with that of a very long time series (4000 years) of land CO2 fluxes simulated by the Community Earth System Model (CESM).
2.1 CO2 datasets for the recent past
We select the CO2 time series datasets from the Global Carbon Budget 2018 version 1.0 (Le Quéré et al., 2018): the atmospheric CO2 growth rate (AGR), the land sink from models (SLDGVMs), the residual land sink (SLResid), and the land sink from two atmospheric inversions.
In the Global Carbon Budget 2018 (Le Quéré et al., 2018), the global CO2 balance is calculated based on the carbon emissions from fossil fuel (FF) (Boden et al., 2017; UNFCCC, 2018; Peters et al., 2011) and land use change (FLUC) (Houghton and Nassikas, 2017; Hansis et al., 2015), the AGR (Dlugokencky and Tans, 2018), and the carbon uptake by the ocean sink (SO) and the land sink (SL) (references for individual models of SO and SL can be found in Table 4 of Le Quéré et al., 2018).
The difference of annual atmospheric CO2 in a given year and the previous year (Ballantyne et al., 2012; Dlugokencky and Tans, 2018; Le Quéré et al., 2018) corresponds to the AGR, which is based on direct observations. The AGR is based on the average of well-mixed CO2 measurements at multiple global stations from the US National Oceanic and Atmospheric Administration Earth System Research Laboratory (NOAA ESRL) (Dlugokencky and Tans, 2018).
FF emissions are based on inventories, while FLUC, SL, and SO are estimated by models (the latter two are denoted by SLDGVMs and SO in Eq. 1), all of which contain uncertainties (Le Quéré et al., 2018). The total emissions from FF and FLUC minus AGR should equal the SO and SLDGVMs (Eq. 1). Due to uncertainties in modeled land and/or ocean sinks or in land use estimations (Bastos et al., 2020; Hauck et al., 2020), the budget cannot be balanced, and thus an imbalance term (IMB) is introduced to the budget.
The annual land sink of CO2 (SLDGVMs) is the average net biome production (NBP) simulated by 16 dynamic global vegetation models (DGVMs) forced with historical CO2 and changing climate (Le Quéré et al., 2018). The residual land CO2 (SLResid) is calculated from emissions, AGR, and ocean sinks, as described in Eq. (1). SLResid corresponds to the balance of the fossil fuel and land use change emissions as well as the sinks in the atmosphere and ocean, and it provides an alternative estimate of the global land sink.
The time series of AGRR, SLDGVMs, and SLResid in the Global Carbon Budget 2018 (Le Quéré et al., 2018) are provided at annual time steps over the period 1959–2017. In the following analysis, we invert the AGR time series (AGRR for reversed AGR, i.e., ) for sign consistency with the land sink datasets used (defined as a positive flux from the atmosphere to the land).
Additionally, we use the globally aggregated net atmosphere to land CO2 flux (positive sign as a sink in the biosphere) estimated from two atmospheric CO2 inversions in the Global Carbon Budget 2018 (Le Quéré et al., 2018): the Jena CarboScope SLCarboScope (Rödenbeck, 2005; Rödenbeck et al., 2018) and the Copernicus Atmosphere Monitoring Service inversion SLCAMS (Chevallier et al., 2005), which cover the periods 1976–2017 and 1979–2017, respectively. Here we use the global annual CO2 fluxes of these two inversions adjusted for fossil fuel emissions and lateral fluxes from Bastos et al. (2020). The period common to the CO2 time series (1980–2017) is selected.
2.1.1 Sea level pressure
We use global monthly mean SLP fields from ERA5 reanalysis with the spatial resolution of (Bell et al., 2020) at monthly time steps covering the period 1950–1978 (Bell et al., 2020) and 1979–present (Hersbach et al., 2019). The period common to other datasets of 1958–2017 is selected here.
2.1.2 Teleconnection indices
In addition to SLP fields, we select 15 teleconnection indices from the atmosphere–ocean variability, Northern Hemisphere, and Southern Hemisphere.
Three important atmosphere–ocean coupled variability modes influence global climate and the C-cycle: the El Niño-Southern–Oscillation (SOI), the Pacific Decadal Oscillation (PDO), and the Atlantic Multidecadal Oscillation (AMO) (Zhu et al., 2017).
In the Northern Hemisphere, the most relevant indices are the Arctic Oscillation (AO), the North Atlantic Oscillation (NAO), the Pacific North American pattern (PNA), the East Atlantic (EA), East Atlantic/West Russia (EAWR), the Scandinavian pattern (SCAND), polar–Eurasia (polarEA), and the West Pacific (WP). These indices are calculated and provided by the Climate Prediction Centre (CPC) of the National Oceanic and Atmospheric Administration (NOAA) (CPC, 2008). Detailed information on calculation procedures is described in NOAA CPC (2008) and Barnston and Livezey (1987).
In the Southern Hemisphere, important indices are the Antarctic Oscillation (AAO), the tropical Atlantic dipole (TAD), the Dipole Mode Index (DMI) of the Indian Ocean Dipole, and the Trans Polar Index (TPI).
The teleconnection indices used here have been summarized in Table 1. All the indices are provided as monthly means and selected for the period of 1958–2017 except the AAO, which is available for 1979–2017 only.
(McBride and Nicholls, 1983; Ropelewski and Jones, 1987)(Mantua et al., 1997; Mantua and Hare, 2002)(Mantua et al., 1997)(Rayner et al., 2003; Enfield et al., 2001)(Enfield et al., 2001)(Higgins et al., 2000, 2002)(Barnston and Livezey, 1987)Barnston and Livezey, 1987Chen and Van den Dool, 2003(Barnston and Livezey, 1987)(Barnston and Livezey, 1987)(Barnston and Livezey, 1987)(Barnston and Livezey, 1987)(Barnston and Livezey, 1987)(Mo, 2000)(CPC, 2008)Enfield et al., 1999Reboita et al., 2021(Saji and Yamagata, 2003; Reboita et al., 2021)(Pittock, 1980, 1984; Jones et al., 1999)Henley et al. (2015)2.1.3 Long-term pre-industrial control simulations for statistical benchmarking
Here we select the SLP fields and global net biome production fields (NBP) from simulations by the Community Earth System Model (CESM) version 1.2.2 (in the B1850C5CN configuration), which has been used by Stolpe et al. (2019). This experiment corresponds to a 4000-year control run. The simulation was run at an atmospheric resolution of using the Community Atmosphere Model version 5 (CAM5.3; Neale et al., 2012) with 30 vertical levels. The model consists of fully coupled atmosphere, ocean, sea ice, and land surface components (Hurrell et al., 2013; Meehl et al., 2013b) and did not include dynamic vegetation. This simulation includes no external forcing, so it is ideal to analyze patterns driven by internal variability.
2.2 Data pre-treatment
For all historical datasets (CO2 time series, SLP fields, and teleconnection indices), we first remove years corresponding to volcanic eruptions (1963, 1982, 1983, 1991, 1992). We then pre-treat the datasets as follows.
2.2.1 Trend removal
The long-term trend of CO2 time series is removed by locally weighted scatterplot smoothing (LOWESS) (Cleveland et al., 1991) of the annual time series with a fixed window size that is 25 % of the longer interval (1959–2017) and 45 % for the shorter period (1980–2017). For monthly teleconnection indices, we first calculate the seasonal mean values of DJF, MAM, JJA, and SON; we then remove the seasonal long-term trends by applying the LOWESS as for the CO2 time series and further include DJF and MAM combined (DJF+MAM) as treated in SLP (as described below).
2.2.2 Spatial and temporal aggregation
The monthly mean SLP fields are area-weighted and aggregated to , , and spatial resolution, and the seasonal cycle is removed by subtracting the monthly mean values for each pixel. We then aggregate SLP values in seasonal means for December of the previous year to February (DJF), March–May (MAM), June–August (JJA), and September–November (SON) of each given year and further consider DJF and MAM combined (DJF+MAM), so the number of pixel-based time series (predictors) in DJF+MAM is double DJF. Note that a large fraction of the pixel-based time series of seasonal SLP anomalies show no long-term trend, and the predicted differences between LOWESS detrended and non-detrended SLP are small. Here we keep the analysis of SLP anomalies with no LOWESS detrending. We refer to DJF and MAM as boreal winter and boreal spring.
For the CESM simulations, the SLP fields are originally provided at spatial resolution at monthly mean time steps, which we then resample to spatial resolution. Annual mean NBP is calculated from the monthly fields. NBP and SLP fields are selected for the simulation period of 1000–5000 years.
2.3 Statistical analysis
The overall goal is to characterize annual variations in the global C-cycle that can be explained by large-scale atmospheric circulation variability. Here, the pixel-based time series of SLP anomalies are used as predictors (p≥800) of CO2 time series (n≤ 54 years) in a linear regression model. However, the small sample size relative to the large number of predictors (n<p) can cause severe overfitting problems and result in unstable predictions (Hastie et al., 2009). Moreover, the existing spatial correlations among the neighboring pixels of SLP anomalies might cause multicollinearity among the predictors (von Storch and Zwiers, 1999). The potential multicollinearity problem results in unstable ridge regression coefficients in least square estimation, making it difficult to diagnose the most sensitive spatial patterns of predictors (von Storch and Zwiers, 1999).
Sippel et al. (2019) applied ridge regression to avoid these overfitting and multicollinearity problems. Ridge regression is a regularized linear regression, the fundamental principle of which is to introduce a constraint (hyper-parameter λ) to regularize the varying regression coefficients in least squares estimation (Hastie et al., 2009; Friedman et al., 2010). The regularized variance comes with a compromise of biased predictions and is addressed as the bias–variance trade-off (Hastie et al., 2009). When selecting the best hyper-parameter λ, this trade-off is considered to achieve stable (low variance) while slightly biased predictions (Hastie et al., 2009).
Model performance is evaluated by the R2, Pearson's correlation R, and mean squared error of the original CO2 time series against predicted values. Pearson's correlation R is selected as the main measure of predictability, and the significance P<0.05 is selected. Given the relatively short period (n<60), here we use leave-one-out cross-validation to achieve optimal model training and testing. For each training and testing group splitting, we select the training group as all years excluding 3 consecutive years and the middle year of those 3 is then selected as the test sample. We exclude the preceding and following years to reduce the potential influence of temporal autocorrelation.
A schematic description of the workflow from model training, validation, and selection with the selected interval 1959–2017 is shown in Fig. 1. In Fig. 1a, the global maps represent an example of the spatial distribution of SLP with a resolution of , where m varies in [2, 5, 9]. Each pixel corresponds to a time series of SLP, so that p predictors (xik for i in 1∼n years and k in 1∼p, where p depends on the spatial resolution selected) with n=54 time steps are defined. Each predictor is assigned a coefficient ωk to collectively predict CO2 time series, with length n=54. The cost function (Harrington, 2012) is the sum of all the squared errors of original yi minus estimated (, ω represents the vector of all ωk). At the same time, the constraint function (Harrington, 2012) suppresses the coefficient variations under a regularized range defined by the hyper-parameter λ. Figure 1b shows the model training and validation processes. In this study, SLP is aggregated to , so p=800. Note that to reduce the heavy computational load, when conducting the spatial and temporal sensitivity study as described in Sect. 2.4, the range of λ is lower and with smaller steps than the range and step shown in Fig. 1b. For example, when only selecting the tropical domain of SLP, the number of predictors is much less than 800. So λ is selected in the range [10, 1000] with a step of 50. When using teleconnection indices instead of SLP anomalies, the predictors are equal to or less than 15, and the range of λ is selected from [1, 200] with a step of 2.
The first step is to divide datasets into training and testing groups, as shown in Fig. 1c. The grouped training datasets are then used for model training to tune the best λ through 5-fold cross-validation. The best λ that achieves the optimal prediction (highest R2) is then selected by the model to predict test datasets. The model then starts another iteration with training and testing grouping. Figure 1c describes the leave-one-out training and testing grouping. The years before and after the selected test year are removed for each grouping to reduce the impact of temporal autocorrelations in CO2 time series.
Ridge regression leave-one-out cross-validation is performed using the Python package Scikit-learn “Ridge”, and λ is tuned by Scikit-learn “RidgeCV” (Pedregosa et al., 2011). The global maps (Figs. 3, A4, A11, and A12) are plotted by Cartopy (Met Office, 2010–2015).
2.4 Experimental design
-
Preliminary dependency tests. To evaluate the robustness of the results for different characteristics of the datasets and methodological choices, we perform several preliminary tests. (1) A resolution dependency test is performed to evaluate the sensitivity of results to the SLP spatial resolution under , , and . (2) A seasonality dependency test is performed to evaluate the dependence of results on the definition of particular seasons, with each season being the combination of 3 consecutive months (from November last year to July in the given year). (3) Temporal autocorrelation of the CO2 time series is performed to ensure no significant trend remains in the detrended CO2 time series.
Here we directly use the results following the preliminary dependency test (Appendix A). The spatial resolution does not influence the results considerably (see Fig. A1), and therefore we select SLP spatial resolution given its smaller number of grid points. The seasonal dependency test shows that DJF and MAM are seasonal combinations more representative of boreal winter and spring (see Fig. A2). JJA and SON are found to have lower or no predictability to CO2 time series; therefore, we limit our results to DJF and MAM. As shown in Fig. A3, the temporal autocorrelation of all CO2 time series is mostly less than 0.4 with lag ranging from 1 to 35 years. With a lag of 1 year, absolute values of autocorrelation are below 0.2 so that we can exclude strong temporal autocorrelation effects.
-
Model training and evaluation. We evaluate the predictability of annual CO2 time series using SLP anomalies, teleconnection indices, and SOI independently in the DJF, MAM, and DJF+MAM seasons with the approach described above. We compare Pearson's correlation of observations and predicted values for different periods (1959–2017 and 1980–2017, Fig. 2) and show the corresponding ridge regression coefficient distribution maps (Fig. 3a).
-
Spatial sensitivity study. We evaluate the predictability of historical annual CO2 time series using DJF SLP anomalies under different spatial domains in the periods of 1959–2017 and 1980–2017. Then we use a 30-year sliding window with annual AGRR to depict how the predictability under various SLP domains evolves in the period 1959–2017.
-
Temporal sensitivity study. We evaluate the predictability of annual CO2 time series AGRR, SLDGVMs, and CESM using DJF and MAM SLP anomalies under different time intervals. Sliding windows are employed at time intervals of 15, 20, 30, and 40 years for historical datasets and CESM, as well as 100, 500, and 2000 years for CESM only. For the interval of 100 and 500 years, we use the sliding window of a 50-year step and a 500-year step for the 2000-year interval. The intervals shorter than 100 years are all in the 1-year step. We also evaluate the error rate of the model in each sliding window of 15-, 20-, 30-, and 40-year lengths. The error rate is calculated by the number of invalid predictions that have significance P>0.05 divided by the number of total predictions within a given window.
-
Comparison with empirical orthogonal function analysis. We compare the predictability of AGRR by applying ridge regression and a generalized linear model that uses the principal components of SLP fields estimated by empirical orthogonal function (see, e.g., von Storch and Zwiers, 1999) decomposition as predictors. We first select three different DJF SLP spatial domains (global, 18∘ N–18∘ S, and 18∘ N–72∘ S). For each spatial domain, we select the first 10 components of SLP anomalies. When selecting the global SLP field, the first 10 components explain up to 75 % of the SLP variance (see Fig. A14). We then reconstruct SLP fields based on an increasing number of components from 1 (based on component 1) to 10 components. These components are used as predictors of AGRR in a simple linear regression (also with leave-one-out cross-validation). The periods 1959–2017 and 1980–2017 are included (see Figs. A13 and A14).
3.1 Global IAV patterns
In this section, we test the predictability of the global C-cycle IAV by using different predictors: global SLP fields by ridge regression, teleconnection indices by ridge regression, and the single SOI index by simple linear regression. In this context, each detrended annual CO2 time series (including AGRR, SLResid, SLDGVMs, SLCAMS, and SLCarboScope) is predicted by the above predictors over DJF, MAM, and DJF + MAM separately (Fig. 2). The predictability is evaluated by the Pearson correlation (ρ) between the original and predicted detrended annual CO2 time series. ρSLP, ρTele, and ρSOI represent the predictability by using different predictors of global SLP fields, teleconnection indices, and SOI, respectively. Accordingly, the relevant ridge regression coefficients with the above different predictors are represented as ωSLP, ωTele, and ωSOI.
First, we find that the detrended annual CO2 time series are generally consistent with each other except that the SLDGVMs shows slight deviation (Fig. 2a). We find two anomalous years (1987 and 1998), which show deviations larger than 2 standard deviations in most CO2 time series, both signifying apparent AGR increases and subsequent lower land sink (Fig. 2a). These two years correspond to strong El Niño events, which are usually associated with below-average land CO2 uptake (Keeling et al., 1995; van der Werf et al., 2004; Bonan, 2016).
Global SLP and teleconnection indices show comparable predictive skill of global C-cycle IAV in winter, while teleconnection indices have higher predictive skill in spring (Fig. 2b). In both periods, the value of ρSLP (except SLDGVMs) is higher in DJF (0.51–0.70) than in MAM (0.29–0.60). On the other hand, the values of ρTele are higher in MAM (0.57–0.79) than in DJF (0.53–0.70). The relatively low predictive skill of global SLP anomalies compared to teleconnection indices might result from (1) limited sample size (less than 60 years) and a large number of predictors (p=800) for ridge regression training with global SLP anomalies. But for teleconnection indices and SOI, the predictive skills are much less influenced by the limited sample size due to their limited predictors (p≤15 for teleconnection indices and p=1 for SOI). As we increase the sample size to over 100 years, the predictive skill of SLP anomalies increases considerably, as is shown in temporal sensitivity study (Fig. 6), and (2) the predictive skill of SLP anomalies in explaining global C-cycle IAV can be reduced in domains with large local rather than global impacts of atmospheric variations to land carbon sinks (Jung et al., 2017). In such domains, the SLP anomalies might show a strong relationship to local C-cycle variations but a weaker link to global C-cycle variations. Selecting the SLP domains with a higher contribution to the global C-cycle variability could improve the predictability, as is shown by the analyzes of the sensitivity of the results to the SLP spatial domains (Fig. 4).
We find that the ridge regression using the full SLP fields as predictors yields generally higher predictability than that of linear regression based on principal components of SLP (see Figs. A13 and A14). In winter and in the SLP domain of 18∘ N to 72∘ S, the linear regression based on the leading principal components of SLP shows lower predictability that the ridge regression for different numbers of components retained (see Fig. A13, ρ<0.7 in 1959–2017 and ρ<0.75 in 1980–2017). The lower predictive skill of the linear regression based on the leading components of SLP might be due to (1) principal component analysis capturing the main variances of the SLP field but not necessarily that of the main patterns that influence CO2 IAV and the fact that, (2) mathematically, ridge regression and principal component analysis are deeply connected. Principal component analysis cuts off all components with small variance beyond a certain threshold, while ridge regression shrinks them, which allows for information in low-variance components to be used for prediction (van Wieringen, 2021). Thus, ridge regression might reveal hidden components of SLP variability that are nevertheless important for CO2 IAV.
Compared to the predictive skill of teleconnection indices, which includes a set of 14 teleconnection indices for the period 1959–2017 and 15 for the period 1980–2017 as predictors, the predictive skill of SOI is slightly lower or similar in both seasons, with 0.53–0.67 in DJF and 0.60–0.74 in MAM (Fig. 2). This is consistent with the dominant role of ENSO in driving global C-cycle IAV, with other modes showing fewer contributions. Such interpretation requires caution as the indices cannot fully represent the complex atmospheric dynamics.
The predictive skill of the combined winter and spring global SLP anomalies reveals the different seasonal responses of global C-cycle IAV to large-scale atmospheric circulation variability (Fig. 2). The predictive skill of SLP and teleconnection indices in DJF+MAM is within the values for DJF and MAM for most datasets and slightly higher than the best-performing season for the predictive skill of SOI.
The predictive skill of SLP to SLResid is similar to SLDGVMs in MAM and slightly higher than SLDGVMs in DJF. The difference in the predictive skill of SLP to SLResid and SLDGVMs in DJF may due to the fact that (1) compared to SLResid, land sink IAV simulated by DGVMs is less sensitive to DJF climate forcing (Bastos et al., 2018), and (2) SLResid implicitly includes the variability from land use change as well as ocean sink variations (Dufour et al., 2013; DeVries et al., 2017; Friedlingstein et al., 2019).
We next compare the spatial patterns of the ridge regression coefficients of SLP and teleconnection indices in the period of 1959–2017; the results of the period 1980–2017 can be found in Figs. A4 and A5. The spatial patterns of the ωSLP are similar for the three CO2 time series: positive coefficients over the eastern tropical Pacific Ocean, negative coefficients from Southeast Asia extending to Australia that together are roughly consistent with ENSO, and negative coefficients from the western Pacific (Fig. 3a). In winter the positive coefficients over the eastern tropical Pacific are higher than in other regions, which are influenced by El Niño and La Niña, respectively (Monahan, 2001; Hsieh, 2004; Rodgers et al., 2004; Schopf and Burgman, 2006; Sun and Yu, 2009; Yu and Kim, 2011): El Niño induces negative SLP anomalies over the eastern Pacific and positive SLP anomalies over the western Pacific (see King et al., 2020, Fig. 5). The results are consistent with the land sink being negatively driven by ENSO in winter: strong El Niño, decreased land sink; strong La Niña, increased land sink. In spring, the area over the central and western tropical Pacific also shows stronger coefficients and likely corresponds to a mix of different modes, such as the ENSO, the western Pacific teleconnection, and the Interdecadal Pacific Oscillation, all showing strong coefficients in Fig. 3b (SOI, WP, and TPI indices). In Fig. A11 we show the anomalies in temperature and precipitation associated with these patterns, as well as those in NBP from the two atmospheric inversions (see Fig. A12). Generally, the temperature anomalies over the tropics show negative correlations with annual land sink (SLP-driven AGRR) in both winter (as high as −0.85) and spring (as high as −0.73), while weaker but positive correlations are found in Eurasia. Tropical precipitation anomalies show roughly positive correlations in winter (as high as 0.73) and in spring (as high as 0.67). This pattern indicates that AGRR is generally higher for cooler and wetter conditions over the tropics and Southern Hemisphere semi-arid regions in both seasons, which result in increased NBP (see Fig. A12), as well as cooler but predominantly drier conditions over Eurasia, which result in a complex pattern of NBP anomalies (see Fig. A12). These results are consistent with the strong ENSO fingerprint on the IAV of global CO2 atmospheric growth rate and global land sink, e.g., as pointed out by Piao et al. (2020) and with the importance of southern semi-arid ecosystems (Ahlström et al., 2015), for IAV in the global land sink.
The higher ridge regression coefficients of teleconnection indices are consistent with the high-sensitivity domains of the ridge regression coefficients of SLP corresponding to the patterns of ENSO and WP (Fig. 3b). Our results show that C-cycle IAV reveals a high positive sensitivity to SOI (0.18 to 0.21) and negative sensitivity to TPI (−0.24 to −0.28) and WP (−0.09 to −0.15) in winter. High sensitivities are also found for DMI in winter (negative) and AMO in spring (negative) (Fig. 3b). We find that the global C-cycle IAV is very sensitive to TPI as well as SOI, which is not so obvious for the spatial patterns of SLP ridge regression coefficients. TPI is a hemispheric-scale index defined as the pressure anomaly differences between the locations Hobart (43∘ S, 147∘ E) and Stanley (52∘ S, 58∘ W) (Pittock, 1980, 1984). We find the TPI to be strongly anticorrelated with SOI in winter and spring (−0.89 and −0.85, respectively). This might indicate an amplification of ENSO impacts on C-cycle IAV due to large-scale atmospheric circulation variability in the Southern Hemisphere.
However, the observed patterns of the ridge regression coefficients in teleconnection indices and SLP need to be interpreted with caution, since these patterns are not necessarily independent from each other. For example, the area from Southeast Asia extending to Australia corresponds to a region influenced by several large-scale atmospheric circulation modes: ENSO, the Indian Ocean Dipole (IOD), and the Southern Annular Mode (SAM) (Cleverly et al., 2016). Interactions between these modes have been shown to modulate the occurrence of drought and extreme precipitation in semi-arid areas of Australia and thus induce large interannual variability in gross primary productivity in the region (Cleverly et al., 2016).
Compared to other CO2 datasets, the predictability of SLDGVMs is higher when using SLP fields in MAM as predictors rather than DJF (Fig. 2b). Moreover, ωSLP of SLDGVMs exhibits distinct spatial patterns, especially in winter; ωSLP for SLDGVMs shows higher values in the southern Pacific than over the tropical Pacific region (Fig. 3a). Compared to historical results, predictability of SLDGVMs is lower by using winter SLP as predictors in ridge regression rather than by using spring SLP, and the spatial patterns of the ridge regression coefficients for SLDGVMs are slightly different. These differences might be an indication of shortcomings of DGVMs in simulating the sensitivity of land sink to climatic drivers.
The general match of spatial patterns of the ridge regression coefficients using SLP and the teleconnection indices as predictors of global C-cycle IAV indicates that SLP can capture the spatial distribution of the atmospheric patterns that influence IAV, with the advantage of being more flexible than teleconnection indices, since it does not require predefined definitions. However, the short sample size and the large number of predictors for ridge regression training hinder the performance of SLP anomalies, especially the lower predictability when using SLP anomalies rather than teleconnection indices in spring. Reducing the number of predictors (smaller domains of SLP anomalies) or increasing the sample size (longer time interval) for ridge regression training could improve the predictive skill of SLP anomalies. Therefore, in the next subsections, we conduct a spatial and temporal sensitivity study of the global C-cycle to SLP anomalies.
3.2 Sensitivity to the SLP domains
Here, we test the sensitivity of the global C-cycle IAV predictability by using different spatial domains of SLP as predictors in ridge regression. The SLP domains are selected over different latitudinal bands in DJF and MAM separately (Fig. 4). We find improved predictability in both seasons when selecting smaller spatial domains (particularly including the tropics to high latitudes of the Southern Hemisphere) rather than global SLP anomalies. MAM fields show lower predictability in general and are less sensitive to the spatial domain considered. In the following, we show the results for DJF; the results for MAM can be found in Fig. A7. Here we only show the results of the period 1980–2017. The results of the period 1959–2017 show a similar trend (see Fig. A6).
Consistent with previous studies (Zeng et al., 2005; Piao et al., 2020), the tropical domain corresponds to higher predictability for all datasets, but stronger predictability is found for regions extending from the tropics to the Southern Hemisphere (Fig. 4). Including the Northern Hemisphere regions results in lower predictability. The domain 18∘ N–72∘ S shows the highest predictability, with ρSLP of 0.81 for AGRR and 0.85 for SLResid in 1980–2017.
The results for net atmosphere–land fluxes estimated by atmospheric inversions are consistent with those of AGRR, with ρSLP of 0.70 for SLCAMS and 0.76 for SLCarboScope in the same domain of 18∘ N–72∘ S. The values of ρSLP of SLDGVMs are systematically lower than the other datasets, independently of the domain.
The weaker values of predictability when extending SLP domains from the tropics to the Northern Hemisphere (Fig. 4) might be due to the local rather than global impacts of large-scale atmospheric circulation variability in the Northern Hemisphere to land sink IAV. Additional explanations include the fact that carbon fluxes are weaker in the winter Northern Hemisphere so that large-scale atmospheric circulation variability exerts a weaker influence in the global land sink and that there are strong compensatory effects of gross primary productivity versus terrestrial ecosystem respiration in the Northern Hemisphere in response to water and temperature variations (Jung et al., 2017; Wang et al., 2022).
The increasing predictability of global C-cycle IAV when using SLP fields extending from the tropics to the Southern Hemisphere (Fig. 4) is likely due to the strong contribution of semi-arid regions in the Southern Hemisphere extratropics to the global sink through their drought–wet anomalies (Poulter et al., 2014; Ahlström et al., 2015). The drought–wet anomalies in these regions are controlled by large-scale atmospheric circulation variability in the Southern Hemisphere (ENSO and ENSO related modes) and to the interactions between ENSO and other large-scale atmospheric circulation modes in the Southern Hemisphere, such as the synergistic effects from ENSO, IOD, and the SAM on Australia C-cycle variability (Cleverly et al., 2016).
3.3 Sensitivity to the temporal domains
Because of multidecadal variability in the climate system, it is possible that the relationships found for short intervals are not stable. In order to investigate whether these results depend on the temporal domain considered, we additionally analyze the influence of different temporal domains (30-year interval sliding window) on the predictability of global C-cycle IAV by selecting different temporal domains and performing the ridge regression separately for the global and the tropical domains (Fig. 5).
Results show stronger predictability of AGRR confined to the tropics of SLP domains in earlier periods and an intensification of this predictability for the SLP domains extending to the Southern Hemisphere over the study period. In some periods, the tropics and southern extratropics show the highest values of ρSLP, for example in 1978–2011 and 1984–2015 (Fig. 5). There is, however, high temporal variability in the predictability and the most relevant spatial domain, with other periods showing higher global coherence (e.g., 1968–2001). It is unclear whether these temporal variations occur randomly due to internal variability in the climate system or are influenced by external forcing. Potential explanations for this pattern include trends found in SLP variability over the Pacific and southern Atlantic (Schneider et al., 2012; IPCC, 2013; Roxy et al., 2019) or enhanced sensitivity of C-cycle variability to climatic drivers, particularly in semi-arid areas, under progressive climate change (Wang et al., 2014; Poulter et al., 2014).
Understanding and attributing these changes to given processes are beyond the scope of this study, but these results highlight the importance of the temporal domain when analyzing IAV in the global C-cycle. Since the observed CO2 time series are short and cover only limited temporal domains, results are likely to be affected by multidecadal internal climate variability, in addition to external forcing. Moreover, the data-driven ridge regression method to quantify circulation-induced global C-cycle variability uses a large number of predictors, while only relatively short time series are available for training, which may negatively affect the model's performance. Therefore, we further test the sensitivity of the predictability to the length of the time series (Fig. 6). We test the predictability of global C-cycle IAV for different lengths of the temporal domain: 15, 20, 30, and 40 years for the datasets in the Global Carbon Budget 2018 as well as CESM simulations and 100, 500, and 2000 years for CESM only.
The box plots in Fig. 6 show the distribution of predictability calculated for multiple time intervals, with each time interval using a sliding window over the whole period of the respective time series of AGRR, SLDGVMs, and CESM. The spread of predictability provides an indication of internal variability in the predictability of global C-cycle IAV due to the choice of temporal domain and the uncertainty in the ridge regression fit for a large number of predictors and comparatively small number of training samples.
We find that the longer the time interval the higher the mean predictability and the smaller the variation, i.e., the less dependent the results are on the temporal domain considered (Fig. 6). However, the mean value tends to be lower than the median for intervals shorter than 30 years and similar to the median for longer intervals. The lower mean is influenced by some domains with very low or even negative predictability from invalid predictions in shorter time intervals.
The mean predictability of AGRR in winter for the global domain increases from 0.45 to 0.61 from 15 to 40 years, respectively, while the spread (maximum ρSLP − minimum ρSLP) decreases from 1.04 to 0.18. Predictability for SLDGVMs is consistent with those of AGRR, with systematically lower mean predictability for winter and higher for spring, but similar spread in both. The mean value of predictability for CESM in SLP DJF over the global domain increases from 0.42 (15 years) to 0.57 (40 years) and to 0.84 (2000 years), and the spread decreases from 1.72 to 0.72 and to 0.008, respectively.
At global scale, the predictive skill of SLP anomalies with AGRR and with models from SLDGVMs and CESM is different in winter and spring (Fig. 6). Predictability of AGRR is higher with winter SLP (ρSLP is 0.61 in 40-year DJF), but predictability of SLDGVMs and CESM is higher with spring SLP (ρSLP in 40-year MAM is 0.70 and 0.65, respectively).
When limiting the SLP domain to the tropics, results follow the same patterns as those at the global scale, but with better predictive skill: AGRR shows the highest mean predictability of 0.74 for intervals of 40 years in winter (Fig. 6). SLDGVMs shows the highest mean predictability, with ρSLP of 0.73 for 40 years in spring (Fig. 6), a result that is very similar to those of CESM for the same temporal length (ρSLP=0.70). We find that with different time intervals, tropical SLP in winter leads to higher predictability of AGRR than global SLP fields. The spring tropical SLP only shows slightly higher predictive skill than global SLP for SLDGVMs and CESM. AGRR is highly influenced by winter tropical SLP, while SLDGVMs and CESM are more sensitive to spring tropical SLP (Fig. 6). This is consistent with the results of Fig. 4 even when different time intervals are considered. This might be due to the fact that Earth system models, as CESM used here, have been found not to reliably simulate the seasonal timing of ENSO occurrence (Sheffield et al., 2013). The predictive skill of SLP in seasonal and domain differences between observation-based data and models (SLDGVMs and CESM) could be used as an indicator to reveal their different driving mechanisms.
We evaluate the error rate of valid predictions within each time interval sliding window in Fig. 6 (i.e., the fraction of predictions in one time interval sliding window with ρSLP of significance P>0.05) for AGRR, SLDGVMs, and CESM (see Fig. A10). We find that at least 30-year-long intervals are needed for robust prediction of global C-cycle IAV. For periods shorter than 30 years, the rate of invalid predictions (in the sense given above) can be higher than 40 % for most datasets and SLP domains and seasons. It is worth noting that even when predicting AGRR by using winter SLP in the tropical domain, the error rate can still be as low as 13 % in a 15-year interval. From the 15- to 30-year interval, the error rate of AGRR is reduced from 0.4 to 0 in winter global and 0.13 to 0 in winter tropical. The error rate decreases to less than 0.16 in a 30-year interval except that AGRR decreases to 0.24 in spring global, which also matches the relative low predictability of spring SLP anomalies in the period 1980–2017 with different spatial domains (see Fig. A7b). All error rates are reduced to almost zero in a 40-year interval.
The major objective of this study is to explore the relationship between SLP anomalies (as a proxy for large-scale atmospheric circulation variability) and the global C-cycle IAV. Specifically, our goals are (1) to investigate the skill of SLP to predict global C-cycle IAV using ridge regression as well as to compare with traditional teleconnection indices and (2) to establish statistical links at different spatiotemporal scales between large-scale atmospheric circulation variability and global C-cycle IAV.
First, we find that boreal winter and spring SLP anomalies allow predicting IAV of atmospheric CO2 growth rate and IAV of the global land sink, with correlations between predicted and reference values between 0.70 and 0.84 when using winter SLP fields. This is comparable to or higher than the predictive skill of a similar model using 15 teleconnection indices as predictors. The spatial patterns of the ridge regression coefficients reveal a strong influence of large-scale atmospheric circulation variability on global C-cycle IAV, particularly in the El Niño–Southern Oscillation and western Pacific domains. Second, the comprehensive spatiotemporal sensitivity analysis indicates an increasing sensitivity of global C-cycle IAV to large-scale atmospheric circulation variability during boreal winter in the Southern Hemisphere extratropics in the recent decades. This increased sensitivity may be influenced by internal climate variability or by enhanced sensitivity of global C-cycle variability to externally forced changes but requires further research. Finally, we find that time series of at least 30 years are needed for robust predictability of global C-cycle IAV. For shorter time series, predictability is highly dependent on the particular period considered and thus largely due to statistical artifacts of random internal climate variability in the fitting process.
Overall, ridge regression using seasonal SLP fields as predictors of global C-cycle variability provides a novel and efficient data-driven approach for detecting the relationship of large-scale atmospheric circulation variability to global C-cycle variability. Compared to teleconnection indices, this approach requires no predefined spatial configurations, is more flexible to the particular domain considered, and shows equal or higher predictability of global C-cycle IAV. This method allows quantifying the contribution of atmospheric dynamical processes in driving variability in the C-cycle at global and regional scales and may further be useful for attributing observed changes to internal climate variability versus anthropogenic climate change.
The Python scripts are available at https://doi.org/10.17617/3.DMUQZY (Li, 2022).
The Global Carbon Budget 2018 dataset (Le Quéré et al., 2018) is available at https://www.icos-cp.eu/science-and-impact/global-carbon-budget/2018. The two atmospheric inversion datasets are available in Bastos et al. (2019). The monthly sea level pressure from ERA5 reanalysis available at the Climate Data Store for the period 1959–1978 is from Bell et al. (2020) and the period 1979–2017 from Hersbach et al. (2019). The download links for the 15 teleconnection indices are available in Table 1 of Sect. 2.1.2. The sea level pressure and NBP from CESM1.2 are available in Stolpe et al. (2019).
NL processed the data, performed the analysis, and wrote the first draft of the paper. SS, AB, and MR designed the study. SS pre-processed the CESM1.2 data. All authors contributed to designing the methodology and writing the paper.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We are grateful for the support from the International Max Planck Research School for Global Biogeochemical Cycles (IMPRS) for providing relevant scientific seminars and courses.
Markus Reichstein and Alexander J. Winkler acknowledge support by the European Research Council (ERC) Synergy Grant “Understanding and Modelling the Earth System with Machine Learning (USMILE)” under the Horizon 2020 research and innovation program (grant agreement no. 855187).
The article processing charges for this open-access publication were covered by the Max Planck Society.
This paper was edited by Kira Rehfeld and reviewed by Ingeborg Levin and Pradeebane Vaittinada Ayar.
Ahlström, A., Raupach, M. R., Schurgers, G., Smith, B., Arneth, A., Jung, M., Reichstein, M., Canadell, J. G., Friedlingstein, P., Jain, A. K., Kato, E., Poulter, B., Sitch, S., Stocker, B. D., Viovy, N., Wang, Y. P., Wiltshire, A., Zaehle, S., and Zeng, N.: The dominant role of semi-arid ecosystems in the trend and variability of the land CO2 sink, Science, 348, 895–899, https://doi.org/10.1126/science.aaa1668, 2015. a, b
Bacastow, R. B.: Modulation of atmospheric carbon dioxide by the Southern Oscillation, Nature, 261, 116–118, https://doi.org/10.1038/261116a0, 1976. a, b
Ballantyne, A. P., Alden, C. B., Miller, J. B., Tans, P. ., and White, J. W. C.: Increase in observed net carbon dioxide uptake by land and oceans during the past 50 years, Nature, 488, 70–72, https://doi.org/10.1038/nature11299, 2012. a
Barnston, A. G. and Livezey, R. E.: Classification, seasonality and persistence of low-frequency atmospheric circulation patterns, Mon. Weather Rev., 115, 1083–1126, https://doi.org/10.1175/1520-0493(1987)115<1083:CSAPOL>2.0.CO;2, 1987. a, b, c, d, e, f, g, h
Basile, S. J., Lin, X., Wieder, W. R., Hartman, M. D., and Keppel-Aleks, G.: Leveraging the signature of heterotrophic respiration on atmospheric CO2 for model benchmarking, Biogeosciences, 17, 1293–1308, https://doi.org/10.5194/bg-17-1293-2020, 2020. a
Bastos, A., Running, S. W., Gouveia, C., and Trigo, R. M.: The global NPP dependence on ENSO: La Niña and the extraordinary year of 2011, J. Geophys. Res.-Biogeo., 118, 1247–1255, https://doi.org/10.1002/jgrg.20100, 2013. a
Bastos, A., Friedlingstein, P., Sitch, S., Chen, C., Mialon, A., Wigneron, J.-P., Arora, V. K., Briggs, P. R., Canadell, J. G., Ciais, P., Chevallier, F., Cheng, L., Delire, C., Haverd, V., Jain, A. K., Joos, F., Kato, E., Lienert, S., Lombardozzi, D., Melton, J. R., Myneni, R., Nabel, J. E. M. S., Pongratz, J., Poulter, B., Rödenbeck, C., Séférian, R., Tian, H., van Eck, C., Viovy, N., Vuichard, N., Walker, A. P., Wiltshire, A., Yang, J., Zaehle, S., Zeng, N., and Zhu, D.: Impact of the 2015/2016 El Niño on the terrestrial carbon cycle constrained by bottom-up and top-down approaches, Philos. T. Roy. Soc. B, 373, 20170304, https://doi.org/10.1098/rstb.2017.0304, 2018. a
Bastos, A., Sullivan, M., Ciais, P., Makowski, D., Sitch, S., Friedlingstein, P., Chevalier, F., Rödenbeck, C., Pongratz, J., Luijkx, I., Patra, P., Peylin, P., Canadell, J., Lauerwald, R., Li, W., Smith, N., Peters, W., Goll, D., Jain, A., Kato, E., Lienert, S., Lombardozzi, D., Haverd, V., Nabel, J., Tian, H., Walker, A., and Zaehle, S.: Aggregated regional estimates of net atmosphere-land CO2 fluxes from the five atmospheric inversions and 16 Dynamic Global Vegetation Models, supplemental data to Bastos et al. (2019), ICOS ERIC – Carbon Portal [data set], https://doi.org/10.18160/1SVH-3DNB, 2019. a
Bastos, A., O'Sullivan, M., Ciais, P., Makowski, D., Sitch, S., Friedlingstein, P., Chevallier, F., Rödenbeck, C., Pongratz, J., Luijkx, I. T., Patra, P. K., Peylin, P., Canadell, J. G., Lauerwald, R., Li, W., Smith, N. E., Peters, W., Goll, D. S., Jain, A., Kato, E., Lienert, S., Lombardozzi, D. L., Haverd, V., Nabel, J. E. M. S., Poulter, B., Tian, H., Walker, A. P., and Zaehle, S.: Sources of uncertainty in regional and global terrestrial CO2 exchange estimates, Global Biogeochem. Cy., 34, e2019GB006393, https://doi.org/10.1029/2019GB006393, 2020. a, b, c, d
Bell, B., Hersbach, H., Berrisford, P., Dahlgren, P., Horányi, A., Muñoz Sabater, J., Nicolas, J., Radu, R., Schepers, D., Simmons, A., Soci, C., and Thépaut, J.-N.: ERA5 monthly averaged data on pressure levels from 1950 to 1978 (preliminary version), Copernicus Climate Change Service (C3S) Climate Data Store (CDS), https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels-monthly-means-preliminary, last access: 17 December 2020. a, b, c
Boden, T. A., Marland, G., and Andres, R. J.: Global, Regional, and National Fossil-Fuel CO2 Emissions, Tech. rep., Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tenn., USA, http://cdiac.ornl.gov/trends/emis/overview_2014.html, last access: July 2017. a
Bonan, G. B.: Ecological climatology, concepts and applications, Cambridge University Press, ISBN 978-1-107-04377-0, 2016. a
Chen, W. Y. and Van den Dool, H.: Sensitivity of Teleconnection Patterns to the Sign of Their Primary Action Center, Mon.Weather Rev., 131, 2885–2899, https://doi.org/10.1175/1520-0493(2003)131<2885:SOTPTT>2.0.CO;2, 2003. a
Chevallier, F., Fisher, M., Peylin, P., Serrar, S., Bousquet, P., Bréon, F.-M., Chédin, A., and Ciais, P.: Inferring CO2 sources and sinks from satellite observations: Method and application to TOVS data, J. Geophys. Res., 110, D24309, https://doi.org/10.1029/2005JD006390, 2005. a, b, c
Ciais, P., Sabine, C., Bala, G., Bopp, L., Brovkin, V., Canadell, J., Chhabra, A., DeFries, R., Galloway, J., Heimann, M., Jones, C., Le Quéré, C., Myneni, B. R., Piao, S., and Thornton, P.: Carbon and Other Biogeochemical Cycles, in: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK and New York, NY, USA, https://www.ipcc.ch/site/assets/uploads/2018/02/WG1AR5_all_final.pdf (last access: 20 March 2022), 2013. a
Cleveland, W., Grosse, E., and Shyu, W.: Local regression models, in: Statistical Models in S, Chapman and Hall, 309–376, ISBN 9780412830402, 1991. a
Cleverly, J., Eamus, D., Luo, Q., Coupe, N. R., Kljun, N., Ma, X., Ewenz, C., Li, L., Yu, Q., and Huete, A.: The importance of interacting climate modes on Australia's contribution to global carbon cycle extremes, Scient. Rep., 6, 23113, https://doi.org/10.1038/srep23113, 2016. a, b, c
Cox, P. M., Pearson, D., Booth, B. B., Friedlingstein, P., Huntingford, C., Jones, C. D., and Luke, C. M.: Sensitivity of tropical carbon to climate change constrained by carbon dioxide variability, Nature, 494, 341–344, https://doi.org/10.1038/nature11882, 2013. a
CPC: NOAA/National Weather Service NOAA Center for Weather and Climate Prediction, Climate Prediction Center, https://www.cpc.ncep.noaa.gov/data/teledoc/teleindcalc.shtml (last access: 11 August 2021), 2008. a, b, c
Deser, C., Phillips, A., Bourdette, V., and Teng, H.: Uncertainty in climate change projections: the role of internal variability, Clim. Dynam., 38, 527–546, https://doi.org/10.1007/s00382-010-0977-x, 2012. a
Deser, C., Lehner, F., Rodgers, K. B., Ault, T., Delworth, T. L., DiNezio, P. N., Fiore, A., Frankignoul, C., Fyfe, J. C., Horton, D. E., Kay, J. E., Knutti, R., Lovenduski, N. S., Marotzke, J., McKinnon, K. A., Minobe, S., Randerson, J., Screen, J. A., Simpson, I. R., and Ting, M.: Insights from Earth system model initial-condition large ensembles and future prospects, Nat. Clim. Change, 10, 277–286, https://doi.org/10.1038/s41558-020-0731-2, 2020. a
DeVries, T., Holzer, M., and Primeau, F.: Recent increase in oceanic carbon uptake driven by weaker upper-ocean overturning, Nature, 542, 215–218, https://doi.org/10.1038/nature21068, 2017. a
Dlugokencky, E. and Tans, P.: Trends in atmospheric carbon dioxide, NOAA/ESRL – National Oceanic and Atmospheric Administration, Earth System Research Laboratory, http://www.esrl.noaa.gov/gmd/ccgg/trends/global.html, last access: 4 September 2018. a, b, c
Dlugokencky, E. and Tans, P.: Trends in atmospheric carbon dioxide, NOAA/ESRL – National Oceanic and Atmospheric Administration, Earth Syem Research Laboratory, http://www.esrl.noaa.gov/gmd/ccgg/trends/global.html, last access: 3 November 2019. a
Dufour, C. O., Le Sommer, J., Gehlen, M., Orr, J. C., Molines, J. M., Simeon, J., and Barnier, B.: Eddy compensation and controls of the enhanced sea-to-air CO2 flux during positive phases of the Southern Annular Mode, Global Biogeochem. Cy., 27, 950–961, https://doi.org/10.1002/gbc.20090, 2013. a
Enfield, D. B., Mestas-Nuñez, A. M., Mayer, D. A., and Cid-Serrano, L.: How ubiquitous is the dipole relationship in tropical Atlantic sea surface temperatures?, J. Geophys. Res., 104, 7841–7848, https://doi.org/10.1029/1998JC900109, 1999. a
Enfield, D. B., Mestas-Nuñez, A. M., and Trimble, P. J.: The Atlantic Multidecadal Oscillation and its relation to rainfall and river flows in the continental U.S., Geophys. Res. Lett., 28, 2077–2080, https://doi.org/10.1029/2000GL012745, 2001. a, b, c
Friedlingstein, P., Meinshausen, M., Arora, V. K., Jones, C. D., Anav, A., Liddicoat, S. K., and Knutti, R.: Uncertainties in CMIP5 climate projections due to carbon cycle feedbacks, J. Climate, 27, 511–526, https://doi.org/10.1175/JCLI-D-12-00579.1, 2014. a
Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Hauck, J., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Le Quéré, C., Bakker, D. C. E., Canadell, J. G., Ciais, P., Jackson, R. B., Anthoni, P., Barbero, L., Bastos, A., Bastrikov, V., Becker, M., Bopp, L., Buitenhuis, E., Chandra, N., Chevallier, F., Chini, L. P., Currie, K. I., Feely, R. A., Gehlen, M., Gilfillan, D., Gkritzalis, T., Goll, D. S., Gruber, N., Gutekunst, S., Harris, I., Haverd, V., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Joetzjer, E., Kaplan, J. O., Kato, E., Klein Goldewijk, K., Korsbakken, J. I., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Marland, G., McGuire, P. C., Melton, J. R., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Neill, C., Omar, A. M., Ono, T., Peregon, A., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Séférian, R., Schwinger, J., Smith, N., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Werf, G. R., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2019, Earth Syst. Sci. Data, 11, 1783–1838, https://doi.org/10.5194/essd-11-1783-2019, 2019. a, b
Friedlingstein, P., Jones, M. W., O'Sullivan, M., Andrew, R. M., Bakker, D. C. E., Hauck, J., Le Quéré, C., Peters, G. P., Peters, W., Pongratz, J., Sitch, S., Canadell, J. G., Ciais, P., Jackson, R. B., Alin, S. R., Anthoni, P., Bates, N. R., Becker, M., Bellouin, N., Bopp, L., Chau, T. T. T., Chevallier, F., Chini, L. P., Cronin, M., Currie, K. I., Decharme, B., Djeutchouang, L. M., Dou, X., Evans, W., Feely, R. A., Feng, L., Gasser, T., Gilfillan, D., Gkritzalis, T., Grassi, G., Gregor, L., Gruber, N., Gürses, O., Harris, I., Houghton, R. A., Hurtt, G. C., Iida, Y., Ilyina, T., Luijkx, I. T., Jain, A., Jones, S. D., Kato, E., Kennedy, D., Klein Goldewijk, K., Knauer, J., Korsbakken, J. I., Körtzinger, A., Landschützer, P., Lauvset, S. K., Lefèvre, N., Lienert, S., Liu, J., Marland, G., McGuire, P. C., Melton, J. R., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S.-I., Niwa, Y., Ono, T., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rödenbeck, C., Rosan, T. M., Schwinger, J., Schwingshackl, C., Séférian, R., Sutton, A. J., Sweeney, C., Tanhua, T., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F., van der Werf, G. R., Vuichard, N., Wada, C., Wanninkhof, R., Watson, A. J., Willis, D., Wiltshire, A. J., Yuan, W., Yue, C., Yue, X., Zaehle, S., and Zeng, J.: Global Carbon Budget 2021, Earth Syst. Sci. Data, 14, 1917–2005, https://doi.org/10.5194/essd-14-1917-2022, 2022. a
Friedman, J., Hastie, T., and Tibshirani, R.: Regularization paths for generalized linear models via coordinate descent, J. Stat. Softw., 33, 1–22, https://doi.org/10.18637/jss.v033.i01, 2010. a, b
Frölicher, T. L., Joos, F., Raible, C. C., and Sarmiento, J. L.: Atmospheric CO2 response to volcanic eruptions: The role of ENSO, season, and variability, Global Biogeochem. Cy., 27, 239–251, https://doi.org/10.1002/gbc.20028, 2013. a
Ghil, M.: Natural climate variability, in: Volume 1, The Earth system: physical and chemical dimensions of global environmental change, from Encyclopedia of Global Environmental Change, John Wiley and Sons, Ltd, Chichester, 544–549, http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.15.3765&rep=rep1&type=pdf (last access: 20 February 2022), 2002. a
Gu, G. and Adler, R. F.: Precipitation and temperature variations on the interannual time scale: Assessing the impact of ENSO and volcanic eruptions, J. Climate, 24, 2258–2270, https://doi.org/10.1175/2010JCLI3727.1, 2011. a
Hansis, E., Davis, S. J., and Pongratz, J.: Relevance of methodological choices for accounting of land use change carbon fluxes, Global Biogeochem. Cy., 29, 1230–1246, https://doi.org/10.1002/2014GB004997, 2015. a
Harrington, P.: Machine learning in action, Manning Publications, p. 155, 164, 167, ISBN 9781617290183, 2012. a, b
Harris, I., Osborn, T. J., Jones, P., and Lister, D.: Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset, Sci. Data, 7, 109, https://doi.org/10.1038/s41597-020-0453-3, 2020. a
Hastie, T., Tibshirani, R., and Friedman, J.: The Elements of Statistical Learning: Data Mining, Inferencee, and Prediction, in: Springer Series in Statistics, Springer, 61–67, https://doi.org/10.1007/b94608, 2009. a, b, c, d, e
Hauck, J., Zeising, M., Le Quéré, C., Gruber, N., Bakker, D. C. E., Bopp, L., Chau, T. T. T., Gürses, O., Ilyina, T., Landschützer, P., Lenton, A., Resplandy, L., Rödenbeck, C., Schwinger, J., and Séférian, R.: Consistency and challenges in the ocean carbon sink estimate for the Global Carbon Budget, Front. Mar. Sci., 7, 852, https://doi.org/10.3389/fmars.2020.571720, 2020. a
Henley, B. J., Gergis, J., Karoly, D. J., Power, S. B., Kennedy, J., and Folland, C. K.: A Tripole Index for the Interdecadal Pacific Oscillation, Clim. Dynam., 45, 3077–3090, https://doi.org/10.1007/s00382-015-2525-1, 2015. a
Hersbach, H., Bell, B., Berrisford, P., Biavati, G., Horányi, A., Muñoz Sabater, J., Nicolas, J., Peubey, C., Radu, R., Rozum, I., Schepers, D., Simmons, A., Soci, C., Dee, D., and Thépaut, J.-N.: ERA5 monthly averaged data on pressure levels from 1979 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS), https://doi.org/10.24381/cds.6860a573, 2019. a, b
Higgins, R. W., Leetmaa, A., Xue, Y., and Barnston, A.: Dominant factors influencing the seasonal predictability of U.S. precipitation and surface air temperature, J. Climate, 13, 3994–4017, https://doi.org/10.1175/1520-0442(2000)013<3994:DFITSP>2.0.CO;2, 2000. a
Higgins, R. W., Leetmaa, A., and Kousky, V. E.: Relationships between climate variability and winter temperature extremes in the United States, J. Climate, 15, 1555–1572, https://doi.org/10.1175/1520-0442(2002)015<1555:RBCVAW>2.0.CO;2, 2002. a
Houghton, R. A. and Nassikas, A. A.: Global and re- gional fluxes of carbon from land use and land cover change 1850–2015, Global Biogeochem. Cy., 31, 456–472, https://doi.org/10.1002/2016GB005546, 2017. a
Hsieh, W. W.: Nonlinear multivariate and time series analysis by neural network methods, Rev. Geophys., 42, RG1003, https://doi.org/10.1029/2002RG000112, 2004. a
Humphrey, V., Zscheischler, J., Ciais, P., Gudmundsson, L., Sitch, S., and Seneviratne, S. I.: Sensitivity of atmospheric CO2 growth rate to observed changes in terrestrial water storage, Nature, 560, 628–631, https://doi.org/10.1038/s41586-018-0424-4, 2018. a
Humphrey, V., Berg, A., Ciais, P., Gentine, P., Jung, M., Reichstein, M., Seneviratne, S. I., and Frankenberg, C.: Soil moisture–atmosphere feedback dominates land carbon uptake variability, Nature, 592, 65–69, https://doi.org/10.1038/s41586-021-03325-5, 2021. a
Hurrell, J. W., Holland, M. M., Gent, P. R., Ghan, S., Kay, J. E., Kushner, P. J., Lamarque, J.-F., Large, W. G., Lawrence, D., Lindsay, K., Lipscomb, W. H., Long, M. C., Mahowald, N., Marsh, D. R., Neale, R. B., Rasch, P., Vavrus, S., Vertenstein, M., Bader, D., Collins, W. D., Hack, J. J., Kiehl, J., and Marshall, S.: The community earth system model: a framework for collaborative research, B. Am. Meteorol. Soc., 94, 1339–1360, https://doi.org/10.1175/Bams-D-12-00121.1, 2013. a
IPCC: Climate Change 2013: The Physical Science Basis, in: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, UK and New York, NY, USA, p. 223, 232, 233, 470, 473, 489, 502, 504, 745, 749, 1535, https://www.ipcc.ch/site/assets/uploads/2018/02/WG1AR5_all_final.pdf (last access: 20 March 2022), 2013. a, b, c, d, e, f, g
Jones, P. D., Salinger, M. J., and Mullan, A. B.: Extratropical circulation indices in the Southern Hemisphere based on station data, Int. J. Climatol., 19, 1301–1317, https://doi.org/10.1002/(SICI)1097-0088(199910)19:12<1301::AID-JOC425>3.0.CO;2-P, 1999. a
Jung, M., Reichstein, M., Schwalm, C. R., Huntingford, C., Sitch, S., Ahlström, A., Arneth, A., Camps-Valls, G., Ciais, P., Friedlingstein, P., Gans, F., Ichii, K., Jain, A. K., Kato, E., Papale, D., Poulter, B., Raduly, B., Rödenbeck, C., Tramontana, G., Viovy, N., Wang, Y.-P., Weber, U., Zaehle, S., and Zeng, N.: Compensatory water effects link yearly global land CO2 sink changes to temperature, Nature, 541, 516–520, https://doi.org/10.1038/nature20780, 2017. a, b, c, d, e
Keeling, C. D., Whorf, T. P., Wahlen, M., and van der Plichtt, J.: Interannual extremes in the rate of rise of atmospheric carbon dioxide since 1980, Nature, 375, 666–670, https://doi.org/10.1038/375666a0, 1995. a, b
King, M. P., Yu, E., and Sillmann, J.: Impact of strong and extreme El Niños on European hydroclimate, Tellus A, 72, 1–10, https://doi.org/10.1080/16000870.2019.1704342, 2020. a
Kumar, A. and Hoerling, M. P.: Interpretation and implications of the observed inter-El Niño variability, J. Climate, 10, 83–91, https://doi.org/10.1175/1520-0442(1997)010<0083:IAIOTO>2.0.CO;2, 1997. a
Le Quéré, C., Andrew, R. M., Friedlingstein, P., Sitch, S., Hauck, J., Pongratz, J., Pickers, P. A., Korsbakken, J. I., Peters, G. P., Canadell, J. G., Arneth, A., Arora, V. K., Barbero, L., Bastos, A., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Doney, S. C., Gkritzalis, T., Goll, D. S., Harris, I., Haverd, V., Hoffman, F. M., Hoppema, M., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Johannessen, T., Jones, C. D., Kato, E., Keeling, R. F., Goldewijk, K. K., Landschützer, P., Lefèvre, N., Lienert, S., Liu, Z., Lombardozzi, D., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., Neill, C., Olsen, A., Ono, T., Patra, P., Peregon, A., Peters, W., Peylin, P., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rocher, M., Rödenbeck, C., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Steinhoff, T., Sutton, A., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., Wright, R., Zaehle, S., and Zheng, B.: Global Carbon Budget 2018, Earth Syst. Sci. Data, 10, 2141–2194, https://doi.org/10.5194/essd-10-2141-2018, 2018. a, b, c, d, e, f, g, h, i, j
Li, N.: Inter-annual global carbon cycle variations linked to atmospheric circulation variability_script_NaLi, Edmond [code], https://doi.org/10.17617/3.DMUQZY, 2022. a
Madden, R. A.: Estimates of the natural variability of time-averaged sea-level pressure, Mon. Weather Rev., 104, 942–952, https://doi.org/10.1175/1520-0493(1976)104<0942:EOTNVO>2.0.CO;2, 1976. a
Mantua, N. J. and Hare, S. R.: The Pacific Decadal Oscillation, J. Oceanogr., 58, 35–44, https://doi.org/10.1023/A:1015820616384, 2002. a
Mantua, N. J., Hare, S. R., Zhang, Y., Wallace, J. M., and Francis, R. C.: A Pacific interdecadal climate oscillation with impacts on salmon production, B. Am. Meteorol. Soc., 78, 1069–1080, https://doi.org/10.1175/1520-0477(1997)078<1069:APICOW>2.0.CO;2, 1997. a, b
McBride, J. L. and Nicholls, N.: Seasonal relationships between Australian rainfall and the Southern Oscillation, Mon. Weather Rev., 111, 1998–2004, https://doi.org/10.1175/1520-0493(1983)111<1998:SRBARA>2.0.CO;2, 1983. a
Meehl, G. A., Washington, W. M., Arblaster, J. M., Hu, A., Teng, H., Kay, J. E., Gettelman, A., Lawrence, D. M., Sanderson, B. M., and Strand, W. G.: Climate change projections in CESM1 (CAM5) compared to CCSM4, J. Climate, 26, 6287–6308, https://doi.org/10.1175/jcli-d-12-00572.1, 2013b. a
Met Office: Cartopy: a cartographic python library with a Matplotlib interface, Met Office, Exeter, Devon, https://scitools.org.uk/cartopy (last access: 22 August 2022), 2010–2015. a
Mo, K. C.: Relationships between low-frequency variability in the Southern Hemisphere and sea surface temperature anomalies, J. Climate, 13, 3599–3610, https://doi.org/10.1175/1520-0442(2000)013<3599:RBLFVI>2.0.CO;2, 2000. a
Monahan, A. H.: Nonlinear principal component analysis: Tropical Indo-Pacific sea surface temperature and sea level pressure, J. Climate, 14, 219–233, https://doi.org/10.1175/1520-0442(2001)013<0219:NPCATI>2.0.CO;2, 2001. a
Neale, R. B., Gettelman, A., Park, S., Chen, C. C., Lauritzen, P. H., Williamson, D. L., Conley, A. J., Kinnison, D., Marsh, D., Smith, A. K., Vitt, F. M., Garcia, R., Lamarque, J.-F., Mills, M. J., Tilmes, S., Morrison, H., Cameron, P., Collins, W. D., Lacono, M. J., Easter, R. C., Liu, X., Ghan, S. J., Rasch, P. J., and Taylor, M. A.: Description of the NCAR community atmosphere model (CAM 5.0), NCAR Technical Note No. NCAR/TN-486+STR, NCAR, https://doi.org/10.5065/wgtk-4g06, 2012. a
Newman, M., Alexander, M. A., Ault, T. R., Cobb, K. M., Deser, C., Lorenzo, E. D., Mantua, N. J., Miller, A. J., Minobe, S., Nakamura, H., Schneider, N., Vimont, D. J., Phillips, A. S., Scott, J. D., and Smith, C. A.: The Pacific Decadal Oscillation, Revisited, J. Climate, 29, 4399–4427, https://doi.org/10.1175/JCLI-D-15-0508.1, 2016. a
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., Blondel, M., Prettenhofer, P., Weiss, R., Dubourg, V., Vanderplas, J., Passos, A., Cournapeau, D., Brucher, M., Perrot, M., and Duchesnay, E.: Scikit-learn: Machine Learning in Python, J. Mach. Learn. Res., 12, 2825–2830, 2011. a
Peters, G. P., Minx, J. C., Weber, C. L., and Edenhofer, O.: Growth in emission transfers via international trade from 1990 to 2008, P. Natl. Acad. Sci. USA, 108, 8903–8908, https://doi.org/10.1073/pnas.1006388108, 2011. a
Piao, S., Wang, X., Wang, K., Li, X., Bastos, A., Canadell, J. G., Ciais, P., Friedlingstein, P., and Sitch, S.: Interannual variation of terrestrial carbon cycle: Issues and perspectives, Global Change Biol., 26, 300–318, https://doi.org/10.1111/gcb.14884, 2020. a, b, c, d
Pittock, A. B.: Patterns of climatic variation in Argentina and Chile, I: Precipitation, 1931–60, Mon. Weather Rev., 108, 1347–1361, https://doi.org/10.1175/1520-0493(1980)108<1347:POCVIA>2.0.CO;2, 1980. a, b
Pittock, A. B.: On the reality, stability and usefulness of Southern Hemisphere teleconnections, Aust. Meteorol. Mag., 32, 75–82, 1984. a, b
Poulter, B., Frank, D., Ciais, P., Myneni, R. B., Andela, N., Bi, J., Broquet, G., Canadell, J. G., Chevallier, F., Liu, Y. Y., Running, S. W., Sitch, S., and van der Werf, G. R.: Contribution of semi-arid ecosystems to interannual variability of the global carbon cycle, Nature, 509, 600–603, https://doi.org/10.1038/nature13376, 2014. a, b
Randall, D. A., Wood, R. A., Bony, S., Colman, T., Fichefet, T., Fyfe, J., Kattsov, V., Pitman, A., Shukla, J., Srinivasan, J., Stouffer, R. J., Sumiand, A., and Taylor, K. E.: Climate models and their evaluation, in: Climate Change 2007: The Physical Science Basis, Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Solomon, S., Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tignor, M., and Miller, H. L., Cambridge University Press, Cambridge, UK and New York, NY, USA, 589–662, https://www.ipcc.ch/site/assets/uploads/2018/02/ar4-wg1-chapter8-1.pdf (last access: 23 October 2021), 2007. a
Rayner, N. A., Parker, D. E., Horton, E. B., Folland, C. K., Alexander, L. V., Rowell, D. P., Kent, E. C., and Kaplan, A.: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res.-Atmos., 108, 4407, https://doi.org/10.1029/2002JD002670, 2003. a, b
Reboita, M. S., Ambrizzi, T., Crespo, N. M., Dutra, L. M. M., de S. Ferreira, G. W., Rehbein, A., Drumond, A., da Rocha, R. P., and de Souza, C. A.: Impacts of teleconnection patterns on South America climate, Ann. NY Acad. Sci., 504, 116–153, https://doi.org/10.1111/nyas.14592, 2021. a, b
Rödenbeck, C.: Estimating CO2 sources and sinks from atmospheric mixing ratio measurements using a global inversion of atmospheric transport, Technical Report 6, Max Planck Institute for Biogeochemistry, Jena, http://www.bgc-jena.mpg.de/CarboScope/s/tech_report6.pdf (last access: 1 November 2019), 2005. a
Rödenbeck, C., Houweling, S., Gloor, M., and Heimann, M.: CO2 flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport, Atmos. Chem. Phys., 3, 1919–1964, https://doi.org/10.5194/acp-3-1919-2003, 2003. a, b
Rödenbeck, C., Zaehle, S., Keeling, R., and Heimann, M.: How does the terrestrial carbon exchange respond to inter-annual climatic variations? A quantification based on atmospheric CO2 data, Biogeosciences, 15, 2481–2498, https://doi.org/10.5194/bg-15-2481-2018, 2018. a
Rodgers, K. B., Friederichs, P., and Latif, M.: Tropical Pacific decadal variability and tts relation to decadal modulations of ENSO, J. Climate, 17, 3761–3774, https://doi.org/10.1175/1520-0442(2004)017<3761:TPDVAI>2.0.CO;2, 2004. a
Ropelewski, C. F. and Jones, P. D.: An extension of the Tahiti–Darwin Southern Oscillation Index, Mon. Weather Rev., 115, 2161–2165, https://doi.org/10.1175/1520-0493(1987)115<2161:AEOTTS>2.0.CO;2, 1987. a
Roxy, M. K., Dasgupta, P., McPhaden, M. J., Suematsu, T., Zhang, C., and Kim, D.: Twofold expansion of the Indo-Pacific warm pool warps the MJO life cycle, Nature, 575, 647–651, https://doi.org/10.1038/s41586-019-1764-4, 2019. a
Saji, N. H. and Yamagata, T.: Possible impacts of Indian Ocean Dipole mode events on global climate, Clim. Res., 25, 151–169, 2003. a
Schimel, D., Stephens, B. B., and Fisher, J. B.: Effect of increasing CO2 on the terrestrial carbon cycle, P. Natl. Acad. Sci. USA, 112, 436–441, https://doi.org/10.1073/pnas.1407302112, 2015. a
Schneider, D. P., Okumura, Y., and Deser, C.: Observed Antarctic interannual climate variability and tropical linkages, J. Climate, 25, 4048–4066, https://doi.org/10.1175/JCLI-D-11-00273.1, 2012. a
Schneider, E. K. and Kinter, J. L.: An examination of internally generated variability in long climate simulations, Clim. Dynam., 10, 181–204, https://doi.org/10.1007/BF00208987, 1994. a
Schopf, P. S. and Burgman, R. J.: A simple mechanism for ENSO residuals and asymmetry, J. Climate, 19, 3167–3179, https://doi.org/10.1175/JCLI3765.1, 2006. a
Sheffield, J., Camargo, S. J., Fu, R., Hu, Q., Jiang, X., Johnson, N., Karnauskas, K. B., Kim, S. T., Kinter, J., Kumar, S., Langenbrunner, B., Maloney, E., Mariotti, A., Meyerson, J. E., Neelin, J. D., Nigam, S., Pan, Z., Ruiz-Barradas, A., Seager, R., Serra, Y. L., Sun, D.-Z., Wang, C., Xie, S.-P., Yu, J.-Y., Zhang, T., and Zhao, M.: North American climate in CMIP5 experiments. Part II: evaluation of historical simulations of intraseasonal to decadal variability, J. Climate, 26, 9247–9290, https://doi.org/10.1175/JCLI-D-12-00593.1, 2013. a
Sippel, S., Meinshausen, N., Merrifield, A., Lehner, F., Pendergrass, A. G., Fischer, E., and Knutti, R.: Uncovering the Forced Climate Response from a Single Ensemble Member Using Statistical Learning, J. Climate, 32, 5677–5699, https://doi.org/10.1175/JCLI-D-18-0882.1, 2019. a, b, c
Stenseth, N. C., Ottersen, G., Hurrell, J. W., Mysterud, A., Lima, M., Chan, K.-S., Yoccoz, N. G., and Ådlandsvik, B.: Review article. Studying climate effects on ecology through the use of climate indices: the North Atlantic Oscillation, El Niño Southern Oscillation and beyond, P. Roy. Soc. Lond. B, 270, 2087–2096, https://doi.org/10.1098/rspb.2003.2415, 2003. a
Stolpe, M. B., Medhaug, I., Beyerle, U., and Knutti, R.: Weak dependence of future global mean warming on the background climate state, Clim. Dynam., 53, 5079–5099, https://doi.org/10.1007/s00382-019-04849-3, 2019. a, b
Sun, F. and Yu, J.-Y.: A 10–15-yr modulation cycle of ENSO intensity, J. Climate, 22, 1718–1735, https://doi.org/10.1175/2008JCLI2285.1, 2009. a
UNFCCC: National Inventory Submissions, Tech. rep., UNFCCC, https://unfccc.int/process/transparency-and-reporting/reporting-and-review-under-the-convention/greenhouse-gas-inventories-annex-i-parties/national-inventory-submissions-2018, last access: June 2018. a
van der Werf, G. R., Randerson, J. T., Collatz, G. J., Giglio, L., Kasibhatla, P. S., Arellano Jr., A. F., Olsen, S. C., and Kasischke, E. S.: Continental-scale partitioning of fire emissions during the 1997 to 2001 El Niño/La Niña period, Science, 303, 73–76, https://doi.org/10.1126/science.1090753, 2004. a
van Wieringen, W. N.: Lecture notes on ridge regression, arXiv [preprint], arXiv:1509.09169v7, https://arxiv.org/pdf/1509.09169.pdf (last access: 23 June 2022), 2021. a
von Storch, H.: Analysis of climate variability applications of statistical techniques, Springer-Verlag, Berlin, Heidelberg, p. 10, https://doi.org/10.1007/978-3-662-03167-4, 1995. a
von Storch, H. and Zwiers, W. F.: Statistical analysis in climate research, Cambridge University Press, p. 136, 139, 293, ISBN 0521450713, 1999. a, b, c
Wang, K., Bastos, A., Ciais, P., Wang, X., Rödenbeck, C., Gentine, P., Chevallier, F., Humphrey, V. W., Huntingford, C., O'Sullivan, M., Seneviratne, S. I., Sitch, S., and Piao, S.: Regional and seasonal partitioning of water and temperature controls on global land carbon uptake variability, Nat. Commun., 13, 3469, https://doi.org/10.1038/s41467-022-31175-w, 2022. a, b
Wang, X., Piao, S., Ciais, P., Friedlingstein, P., Myneni, R. B., Cox, P., Heimann, M., Miller, J., Peng, S., Wang, T., Yang, H., and Chen, A.: A two-fold increase of carbon cycle sensitivity to tropical temperature variations, Nature, 506, 212–215, https://doi.org/10.1038/nature12915, 2014. a
Wills, R. C., Battisti, D. S., Hartmann, D. L., and Schneider, T.: Extracting modes of variability and change from climate model ensembles, in: Proceedings of the 7th International Workshop on Climate Informatics: CI 2017, NCAR Technical Note NCAR/TN-536+PROC, NCAR, Boulder, USA, https://climate-dynamics.org/wp-content/uploads/2017/12/ci2017_Wills_et_al.pdf (last access: 30 September 2021), 2017. a
Yu, J.-Y. and Kim, S. T.: Reversed spatial asymmetries between El Niño and La Niña and their linkage to decadal ENSO modulation in CMIP3 models, J. Climate, 24, 5423–5434, https://doi.org/10.1175/JCLI-D-11-00024.1, 2011. a
Zeng, N., Mariotti, A., and Wetzel, P.: Terrestrial mechanisms of interannual CO2 variability, Global Biogeochem. Cy., 19, GB1016, https://doi.org/10.1029/2004GB002273, 2005. a
Zhu, Z., Piao, S., Xu, Y., Bastos, A., Ciais, P., and Peng, S.: The effects of teleconnections on carbon fluxes of global terrestrial ecosystems, Geophys. Res. Lett., 44, 3209–3218, https://doi.org/10.1002/2016gl071743, 2017. a, b, c, d