{'fraction_non_alphanumeric': 0.04506996420436056, 'fraction_numerical': 0.045639440286365116, 'mean_word_length': 4.450339934968962, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 0, 'https://': 2, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Impact of Changing Winds on the Mauna Loa CO2 Seasonal Cycle in Relation to the Pacific Decadal Oscillation

Abstract 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.

The SCA increase is an important indicator of enhanced growing-season ecosystem productivity over the northern temperate and boreal forest (Graven et al., 2013;Lin et al., 2020). Increasing SCA has been used as an emergent constraint on CO 2 fertilization of primary production (Wenzel et al., 2016), which is possible because simulations from land ecosystem models that are coupled with atmospheric transport models attribute much of the simulated SCA increase to CO 2 fertilization (Bastos et al., 2019;Ito et al., 2016;Piao et al., 2018;Zhao et al., 2016). However, most of these coupled models systematically underestimate the SCA increase at high latitudes and in the midtroposphere (500 mbar) from 1960 to 2010, suggesting that the ecosystem models are missing key mechanisms (Graven et al., 2013;Thomas et al., 2016). Possibilities include underestimating CO 2 fertilization on light-use efficiency (Haverd et al., 2020;Thomas et al., 2016), or underestimating response of changing land carbon cycle to climate change. Early studies which examined SCA changes up until the mid-1990s attributed the SCA increase to enhanced photosynthetic CO 2 uptake due to warming-induced lengthening of the growing season, especially in spring (Keeling et al., 1996;Randerson et al., 1999). In contrast, recent studies focusing on trends after 1990 suggest that the CO 2 uptake has become increasingly limited by water stress tied to continued warming (Jiao et al., 2021;Lian et al., 2020;Peñuelas et al., 2017;Wang et al., 2018;Zhang et al., 2020). The SCA trend has likely also been strongly influenced by warming-driven changes in vegetation cover (Forkel et al., 2016;Liu et al., 2020), and changes in agriculture and other land-use changes (Gray et al., 2014;Zeng et al., 2014). Interannual to decadal climate change, including soil temperature and water supply availability, has been suggested as a major source of SCA variability (Buermann et al., 2007).
Our ability to infer changes in growing-season productivity or constrain CO 2 fertilization from changes in SCA may be complicated, however, by variations in atmospheric circulation. Taguchi et al. (2003) and Lintner et al. (2006) found that the MLO SCA is sensitive to the interannual variations in the airflow from Siberia, where the biospheric CO 2 signal is especially strong. Murayama et al. (2007) found significant year-to-year variabilities of the seasonal cycle at two high-latitude stations (BRW and Alert) caused by changing atmospheric circulation alone. Higuchi et al. (2002) showed that atmospheric transport alone produced significant variations in SCA on interannual to quasi-decadal time scales at both MLO and at high-latitude station BRW. Wang et al. (2020) showed that, since 1959, MLO has become more influenced by CO 2 exchanges from Amazonia, which could be connected with an expansion of the Hadley cell in the Southern Hemisphere. These studies, which explore the connection between SCA changes and atmospheric circulation, however, did not discuss connections between decadal climate modes and wind influences on SCA, and did not discuss possible impacts of long-term wind shifts within the Northern Hemisphere.
Here we study the contribution of varying winds and varying ecological CO 2 fluxes to the trend and variability of SCA, using surface-to-atmosphere CO 2 fluxes transported by the Atmospheric Tracer Model 3 (TM3, Heimann & Körner, 2003), and the MIROC4.0-based Atmospheric Chemistry-Transport Model (ACTM, Patra et al., 2018). We focus primarily on impacts at MLO, but also briefly examine impacts on other low-latitude stations as well as high-latitude station BRW. For MLO, we examine the correlation between MLO SCA and several climate indices, including the Pacific Decadal Oscillation (PDO; Mantua et al., 1997), a closely related wind index U-PDO (based on the first EOF of subtropical 700 mbar winds), land temperature and Palmer Drought Severity Index (PDSI). We also examine MLO 10-day back trajectories under different phases of the PDO to better understand transport pathways. We discuss these results in connection with changes in MLO SCA variability and circulation changes on the interannual to decadal scale and the long-term trend. The results highlight the significant role of changing circulation on the variability and SCA trend at MLO as well as other low-latitude stations, which complicates using SCA at these stations to quantify large-scale ecological changes (e.g., Wenzel et al., 2016).

CO 2 Observation and Its SCA
We use monthly averaged atmospheric in situ CO 2 measurements at MLO, from the Scripps CO 2 Program (Keeling et al., 2001). All the missing values from February to April in 1964 are filled by linear interpolation. We compute the SCA by using the Thoning function following the steps below (Thoning et al., 1989): 1. Fit the observed MLO SCA with a cubic polynomial for the trend and 4-harmonics for the climatological seasonal cycle (Equation 1). Each time series is detrended by removing its long-term linear trend. We also show the correlation (r) between each detrended simulation and detrended observations, and the corresponding p-value based on 1,000 iterations of two-tailed random phase test of each correlation (Ebisuzaki, 1997 where t is the decimal year; are the coefficients for the long-term trend; and are the coefficients for the climatological seasonal cycle; and ( ) is the residual term.
2. Apply a Butterworth band-pass filter to the residual ( (t)) with 90-day and 540-day cutoffs.
3. Compute the seasonal amplitude of the sum of the harmonic and band-passed residual, using the peak and trough values for each calendar year.
We note that our choice of parameters is slightly different compared to the default setup of the Thoning function, but the impact on calculated SCA is negligibly small with RMSE of 0.023 ppm from 1959 to 2019. We then compute the long-term growth rate (ppm decade −1 ) of CO 2 SCA using generalized least-square (GLS) regression (Fox & Weisberg, 2019;Pinheiro et al., 2021). We choose this method in order to address serial correlation, as we find significant first-order autocorrelation of the SCA time series suggested by partial autocorrelation function.

Modeled MLO CO 2 Seasonal Cycle
We simulate CO 2 concentrations at MLO using two models. The first model is the global Atmospheric Tracer Model 3 (TM3; Heimann & Körner, 2003), which was run at 5° in longitude and 4° in latitude with 19 vertical levels, using 6-hourly NCEP/NCAR reanalysis meteorological fields (Kalnay et al., 1996) from 1951 to 2019, and surface fluxes from the Jena CarboScope CO 2 Inversion (version ID: s57Xoc_v2020, Rödenbeck et al., 2003Rödenbeck et al., , 2018. The first six years (1951)(1952)(1953)(1954)(1955)(1956)  This correction only has a small impact on our analysis, since we primarily focus on the SCA variability relative to each other, and the long-term SCA trend.
In order to account for uncertainties of the transport model and reanalysis wind fields, we also simulate CO 2 with MIRCO-ACTM (Patra et al., 2005(Patra et al., , 2018). The ACTM model was run at ∼2.8° in longitude and ∼2.8° in latitude with 67 vertical levels driven by JRA-55 winds (Kobayashi et al., 2015) from 1981 to 2019, with the first four years (1981)(1982)(1983)(1984) as model spin-up period. The simulation (ACTM-W) uses cyclostationary fluxes for the ocean from Takahashi et al. (2009), and for the land from CASA (Randerson et al., 1997). It also uses fossil-fuel fluxes from Jones et al. (2021) which vary both by month and by year. The impact of changing fossil-fuel emissions on MLO SCA is very small (Nevison et al., 2008;Wang et al., 2020). We also found the impact of biomass burning to the MLO SCA is negligibly small in the ACTM simulations.
For the modeled MLO CO 2 time series, we compute SCAs and the growth rate of the amplitude using the same approach as for observations (Section 2.1).

Climate Data
As in Buermann et al. (2007), we compute 11-year moving-window correlations between MLO SCA and several climate indices to understand which time periods show the strongest correlation. These climate indices include PDO, land temperature, and PDSI, with each time series detrended by removing the long-term linear trend and smoothed by a five-point binomial filter to focus on interannual to decadal variability. The p-value of correlations is computed using 1,000 iterations of 2-tail random phase test (Ebisuzaki, 1997 We use the monthly land temperature data at 1º × 1º spatial resolution provided by Berkeley Earth (Muller et al., 2013). We use the monthly PDSI at 0.5º × 0.5º spatial resolution provided by the Climatic Research Unit (Van Der Schrier et al., 2013). We compute the growing-season average land temperature and PDSI as the average from May to October of each year, and an area-weighted average for Eurasia (LandT EU , PDSI EU ) and North America (LandT NA , PDSI NA ). Eurasia is defined as all land within 30 • -70 • N and 120 • -160 • E, while North America is defined as all land between 20 • -70 • N over the North America land.
We also generate a wind-based climate index, which we call the U-PDO, based on the first EOF of the 700 mbar monthly wind field from NCEP/NCAR within the domain of 10 • -40 • N and 120 • E-100 • W (method see Legler, 1983 and Appendix Text S1 in Supporting Information S1). The U-PDO is the U wind component (West-East) of this first EOF, as discussed below in Figure 4 and Section 3.3. The choice of 700 mbar is to approximate the pressure level of MLO at ∼680 mbar.

Backward Trajectories
We calculate 10-day back trajectories at MLO using the HYSPLIT program (Draxler & Rolph, 2010;Stein et al., 2015) driven by NCEP/NCAR winds every 6 hr from 1985 to 2019 (Kalnay et al., 1996). The choice of 10-day conforms with prior studies of air-mass originates at MLO (Harris & Kahl, 1990;Harris et al., 1992;Hess, 1996;Miller, 1981). Particles are released at 19.53°N, 155.58°W, and 3,437.00 m.a.s.l. The vertical motion is forced to be along isentropes. All trajectories that reach above 10,000 m altitudes are excluded in our analysis. We compute the fraction of backward trajectories that end above Eurasia or North America on day 10 as a function of the month and the monthly PDO index value.

Observed and Modeled SCA at MLO
In Figure  In Figure 1b, we show the two components of MLO SCA, one driven by varying winds alone (Jena-W), and one driven by varying ecological CO 2 fluxes alone (Jena-F). Jena-F shows a persistent increase of SCA over time with a growth rate at 0.31 ± 0.024 ppm decade −1 over the whole record and 0.29 ± 0.086 ppm decade −1 since 1986. In contrast, the Jena-W shows an SCA decrease of −0.10 ± 0.022 ppm decade −1 over the whole record, accelerating to −0.15 ± 0.057 ppm decade −1 after 1986. In comparison, ACTM-W simulates a relatively smaller SCA decrease of −0.09 ± 0.068 ppm decade −1 from 1986 to 2019 (see Table 1, not shown in figure). The growth rate of the SCA trend based on observation, Jena-WF, Jena-W, and ACTM-W is also summarized in Table 1.
In Figure 1c, we show detrended time series for the observations, Jena-WF, Jena-W, and Jena-F, where the detrending involves subtracting a long-term linear trend. Projecting the detrended Jena-W and detrended Jena-F onto the detrended Jena-WF (method see Text S2 in Supporting Information S1), we find the varying winds alone and varying ecological CO 2 fluxes alone play equal roles in driving the SCA variability from 1959 to 2019, accounting for 51% and 49% of the variability, respectively. The SCA variability explained by winds alone increases slightly to 60% during the period from 1986 to 2019.

Correlation Between MLO CO 2 SCA and Climate Indices
To understand the sources of variability of MLO SCA, we examine the correlation between the observed MLO SCA and a suite of climate indices in Figure 2. This analysis duplicates the method of Buermann et al. (2007), but updates to 2019 and adds PDO and U-PDO to the suite of indices. Over the full record from 1959 to 2019, the correlation with PDO (r = 0.71, p-value < 0.005) and U-PDO (r = 0.52, p-value < 0.005) are both relatively strong. A marginally significant correlation is also found between MLO SCA and PDSI NA (r = 0.42, p-value = 0.021). There is no significant correlation between MLO SCA and other climate indices (PDSI EU , LandT EU , and LandT NA ).
The moving-window correlation with the PDO remains uniformly high (r ≳ 0.5) throughout the record, except for the period from 1984 to 1995. The correlation with U-PDO shows similar time variations. There is no persistent high correlation between MLO SCA and other climate indices. As noted in Buermann et al. (2007), the correlation between the MLO SCA and PDSI NA increased after the late 1970s, with a high correlation found from 1980 to 1990, from 1995 to 2004. Our analysis also shows that correlation with the PDSI NA is similarly high over the recent period from 2011 to 2019.

Correlation Between Components of Amplitude Variability and PDSI or PDO
To address whether the SCA-PDSI NA and SCA-PDO correlations are flux-driven or wind-driven, we examine the separate correlations with ACTM-W, Jena-W, and Jena-F, as shown in Figure 3. Each time series is detrended, normalized, and smoothed, using the same method as described in Section 3.2. As shown in Figure 3a, we find a high correlation between PDO and ACTM-W (r = 0.68, p-value = 0.007), and Jena-W (r = 0.73, p-value = 0.008), showing the importance of wind shifts in driving variability in SCA associated with the PDO. The correlation between PDO and Jena-W is stronger since 1980 (r = 0.89, p-value < 0.005). The 11-year moving-window correlation in the lower panel of Figure 3a suggests that the connection between PDO and Jena-W is relatively low before 1977 and around 1990, while the connection between PDO and ACTM-W is relatively low before 1993 and from 2003 to 2009. The PDO index shows a relatively weak correlation between Jena-F at 0.31 (p-value = 0.142), as in Figure 3a. The 11-year moving-window correlation (lower panel of Figure 3a) switches between positive and negative, with a positive high correlation only found around 1975 and since 2008, and a negative high correlation found from 1979 to 1994. Replacing PDO with U-PDO, the moving-window correlations between ACTM-W, Jena-W, and Jena-F are very similar ( Figure S1 in Supporting Information S1).
As shown in Figure 3b, we do not find a high correlation between PDSI NA with either Jena-W (r = 0.29, p-value = 0.302) or Jena-F (r = 0.06, p-value = 0.958) over the full record from 1959 to 2018. However, the 11-year moving-window correlation in the bottom panel suggests a relatively tight connection between Jena-W and PDSI NA around 1980NA around , around 1990NA around , from 1995NA around to 2005, and from 2009 onward. We also find a relatively high correlation between PDSI NA and ACTM-W (r = 0.56, p-value = 0.041) from 1986 to 2018, with a relatively low connection before 1995.  (Schuldt et al., 2021). Trends are computed using GLS method (see Section 2.1).

Large-Scale Wind-Patterns Associated With PDO Variability
To explore the mechanism by which wind shifts impact the MLO SCA variability in association with the PDO, we examine the related U-PDO index, defined based on the first EOF of the Pacific wind field at 700 mbar (Section 2.3). The spatial pattern of the wind EOF is shown in Figure 4a and the temporal variations are shown in Figure 4b. The temporal variation demonstrates a very close association between the U-PDO and the PDO (r = 0.83, p-value < 0.005). The high temporal variability of the east-west (U) component of the EOF (i.e., the U-PDO) shows the importance of shifting zonal wind strength in connection with the PDO. When the PDO and U-PDO are positive, a stronger westerly anomaly occurs in the domain from 20 • N to 40 • N, together with a stronger northerly anomaly in the spatial domain west of 160 • E but a stronger southerly anomaly in the spatial domain from 130 • W to 160 • W. This pattern makes MLO more exposed to air transported from Eurasia in a positive PDO year. We note that, around MLO, the flow pattern shows both western and eastern air-mass origins at all times of a year, but generally progresses from more westerly in winter to easterly in summer (Harris & Kahl, 1990;Harris et al., 1992).
Similar insight is found using HYSPLIT back trajectories at Mauna Loa ( Figure 5 and Figure S2 in Supporting Information S1), which show a larger Eurasian footprint and smaller North American footprint associated with positive PDO (Figure 5 and Figure S2 in Supporting Information S1). The back trajectories also resolve seasonality in this pattern. From September to May, the back trajectories favor Eurasian air-mass origins, and from July to August, they favor North American origins. This is true regardless of the phase of the PDO, but the North American influence is strengthened at all seasons during low PDO years, while the Eurasian influence is strengthened from September to May during high PDO years. Analysis of individual back trajectories ( Figure S2 in Supporting Information S1) shows that during the winter (DJF) of positive PDO years, most trajectories originate to the west of MLO centered between 20 • N and 40 • N, suggesting strong zonal flow. In contrast, during the winter of negative PDO years, the trajectories originate over a much wider range of latitudes, both north and south. These winters still have a cluster of trajectories from the west, typically following a curved pathway involving eastward flow at higher latitudes (e.g., 30-50 • N) and bending south or southwest on the approach to MLO.
The PDO-related wind shifts are closely associated with the strength of the Aleutian low (Newman et al., 2016), as shown in Figure 6. In a positive PDO year (Figure 6a), the Aleutian low strengthens and moves equatorward, exposing MLO to the influence of westerly winds. In a negative PDO year (Figure 6b), the Aleutian low shrinks and moves poleward, leaving MLO more exposed to air from North America and the eastern Pacific.

The Impact of Varying Winds on the SCA Trend
Our analysis allows dividing the observed MLO SCA increase of 0.22 ± 0.034 ppm decade −1 (3.3 ± 0.51% decade −1 ) from 1959 to 2019 into components driven by variable winds and variable ecological CO 2 fluxes. From winds alone, we find a decreasing trend of −0.10 ± 0.022 ppm decade −1 (−1.5 ± 0.33% decade −1 ) since 1959 based on Jena-W simulations, offsetting ∼30% of the increase of 0.22 + 0.10 = 0.3 Each time series (points) is first detrended by removing the long-term linear trend and normalized by its standard deviation. We also smooth each time series using a five-point binomial filter, shown as curves. We compute the correlation between each smoothed time series and MLO SCA (observation) in the right of the upper panel, while the p-value is computed by 1,000 iterations of two-tailed random phase test (Ebisuzaki, 1997). In the bottom, we show the 11-year moving-window correlation between each smoothed time series and MLO SCA (observation).
The moving-window correlations are plotted by using the same color assignments as in the upper panel. The correlation and moving-window correlation between MLO SCA (Observation) and land temperature are computed with land temperature leading by 1 year, following Keeling et al. (1996). Black dashed lines in the lower panel denote the ± 0.5 value.

10.1029/2021JD035892
8 of 15 2 ppm decade −1 (4.8% decade −1 ) expected from ecological flux changes alone. This decreasing trend accelerates to −0.15 ± 0.057 ppm decade −1 from 1986 to 2019. Simulations with ACTM-W suggests a slightly smaller decreasing trend (−0.09 ± 0.068 ppm decade −1 ) from 1986 to 2019. Over the whole record from 1959 to 2019, our results based on Jena-W show a larger wind effect than found by Wang et al. (2020), who suggested a negligibly small overall trend due to changing winds alone. Over the period from 1986 to 2019, the impacts of wind in our study are consistent with Wang et al. (2020), who reported a trend of −0.1 to −0.2 ppm decade −1 . This wind impact appears to be an important factor contributing to the relatively smaller SCA increase at MLO (low latitude at 19.5°N) compared to BRW (high latitude at 71.2°N), where Jena-W suggests an insignificant wind-induced trend of 0.01 ± 0.054 ppm decade −1 (p-value = 0.97 based on Cox-Stuart trend test) from 1961 to 2019. Similarly, ACTM-W suggests an insignificant trend at BRW of −0.07 ± 0.127 ppm decade −1 (p-value = 0.83 based on Cox-Stuart trend test) from 1986 to 2019, as shown in Table 1. The long-term winds impact at MLO is mostly driven by reduced Eurasia influence, with additional small contributions from reduced North America influence ( Figure S3 in Supporting Information S1).
Wind shifts also lead to a smaller SCA at other low-latitude stations, offsetting the observed increases. As shown in Table 1 Wind shifts thereby offset 16% (at MID), 30% (at KUM), and 50% (at CHR) of the SCA trend expected from ecological flux changes alone. ACTM-W also suggests decreasing SCA from 1986 to 2019 at these three stations, but the magnitude is smaller than that of Jena-W (Table 1). We note that, although Jena-WF underestimated the SCA trend at MID by 50% compared to observation (Table 1), this discrepancy is not necessarily due to error in winds, as it might alternately result from errors in the CO 2 flux fields.
The impact of winds at KUM is of importance in relation to the study of Wenzel et al. (2016), who applied the SCA trend at KUM as an emergent constraint on the magnitude of global CO 2 fertilization effect on gross primary production. However, their study did not evaluate the impact of wind shifts on this constraint.  Figure 2. We compute the correlation between each smoothed time series and the corresponding target (thick black curve) in the right of the upper panel. The target is annual average PDO in (a), and PDSI NA in (b). The p-value of the correlation is computed by 1,000 iterations of two-tailed random phase test (Ebisuzaki, 1997). In the bottom panel, we show the 11-year moving-window correlation between each smoothed time series and corresponding target (thick black curve in the upper panel). The moving-window correlations are plotted by using the same color assignments as in the upper. Black dashed lines in the lower panel denote the ± 0.5 value.
The long-term wind-driven influence on MLO SCA trend and at other stations may have a connection with a poleward shift and weakening of the subtropical jet stream in the Northern Hemisphere (Archer & Caldeira, 2008). This change leads to a weakened influence of fluxes from northern high latitudes of both Eurasia and North America at MLO ( Figure S3c in Supporting Information S1). These northern landmasses are the main contributors to the MLO SCA variability and trend ( Figure S3b in Supporting Information S1), explaining 58.6% of the total SCA variabilities (method see Text S2 in Supporting Information S1).
The changes of subtropical jet stream are also reflected in the U-PDO index (Figure 4). A notable feature of the U-PDO, which is distinct from the PDO, is a stronger downward trend since 1980, which indicates a progressive weakening of westerly anomalies at 700 mbar within the domain from 20 • N to 40 • N. To assess the contribution of the U-PDO index to the long-term wind-driven SCA trend since 1980, we use a simple statistical model SCA(year) = ⋅ year + + ⋅ U − PDO(year) , where U − PDO is the annual average U-PDO. The parameters are computed using GLS method (similar to Section 2.1). The parameter (0.21 ± 0.105 ppm U-PDO −1 ) is largely independent of , indicated by a negligibly small correlation between and (| | < 0.05), and the same model with detrended SCA only leads to a slightly smaller (0.19 ± 0.132 ppm U-PDO −1 ). Using β = 0.21, the trend  in U-PDO from 1980 to 2019 yields an SCA trend of −0.07 ± 0.032 ppm decade −1 which explains ∼44% of the full wind-driven effect estimated by Jena-W over this period. Over the full record from 1959 to 2019, the U-PDO term yields an SCA trend of −0.01 ± 0.041 ppm decade −1 , which is not significant. The poleward shift of the subtropical jet stream is possibly tied to a long-term expansion of the Hadley cell. Wang et al. (2020) argued that a southward expansion of the Hadley cell in the Southern Hemisphere made MLO more influenced by CO 2 fluxes from Amazonia, which reduce the CO 2 SCA at MLO because the phase is approximately opposite of the signal from the Northern Hemisphere. Our analysis partly agrees with Wang et al. (2020), but we only find small increasing tropical influence since 2010, as shown in Figure S3 in Supporting Information S1. Here, we suggest that the northward displacement of the northern boundary of the Hadley cell may also have been important ( Figure S5 in Supporting Information S1, method see Text S3 in Supporting Information S1 and Adam et al., 2018), associated with northward displacement of Aleutian low under greenhouse warming, corresponding to a wider region of easterly trade winds (Gan et al., 2017;Grise et al., 2019;Staten et al., 2019).
One important caveat in our study is the reliability of winds fields and transport models. The Jena-W (based on TM3 transport model and NCEP reanalysis), the ACTM-W (based on MIROC-ACTM model and JRA-55 reanalysis) simulations and the Wang et al. (2020) study all yield slightly different estimates of the wind impact on the MLO SCA trend from 1986 to 2019 (Table 1). Even larger uncertainties may apply to our estimated trends from 1959 to 2019, which entirely depends on TM3 and NCEP reanalysis. NCEP might have larger errors before the satellite era starting in 1979 (Bromwich & Fogt, 2004). Although climate models and reanalysis products generally show an expanding Hadley cell in the Southern Hemisphere, the similar expansion in the Northern Hemisphere is highly uncertain, with a larger trend from reanalyzes than from model ensemble means (Allen & Kovilakam, 2017;Allen et al., 2014;Grise et al., 2019;Johanson & Fu, 2009;Lu et al., 2007;Staten et al., 2018;Tao et al., 2016).
Will shifting winds continue to contribute to decreasing SCA at low latitudes in the future? The answer may depend on whether the Hadley cell expansion over the past decades has resulted from an internal climate mode or from anthropogenic forcing. This topic remains controversial, with some studies arguing for connection with increasing greenhouse gases (Lu et al., 2007(Lu et al., , 2009, and anthropogenic aerosols (Allen & Ajoku, 2016;Allen et al., 2012Allen et al., , 2014, and others pointing to connections with PDO and ENSO indices (Allen & Kovilakam, 2017;Allen et al., 2014;Amaya et al., 2018;Grassi et al., 2012).

The Impact of Varying Winds on the Variability of MLO SCA
We find that varying winds alone contribute significantly (∼51%) to the interannual/decadal variability of MLO SCA. Much of the wind-related variability is tied to the PDO, which modulates the relative impact of Eurasian versus North American influences on the MLO SCA variability. Years with high PDO tend to have large Eurasian influence and larger seasonal cycles. This result is consistent with Lintner et al. (2006), who reported that during the months from April to June (near the seasonal CO 2 peak at MLO), synoptic patterns favoring strong Eurasia influence are associated with high springtime CO 2 anomalies, and patterns favoring North American influence are associated with low springtime CO 2 anomalies. Buermann et al. (2007) similarly attributed part of the long-term MLO SCA decreasing in the early 1990s to the progressively declining Eurasia influence in the northern spring. However, neither of these previous studies noted the close association with the PDO. As shown in Figure S6 in Supporting Information S1, we also find high correlation between MLO SCA and ENSO index (NOAA, 2019) in some periods, but not consistently over the whole period from 1959 to 2019. Lintner et al. (2006) suggested that the connection between SCA and ENSO index can be explained by the flow patterns associated with ENSO phases. This connection can also be tied to the PDO phases, as PDO resembles ENSO on the decadal scale.
We find that the contribution from variations in ecological CO 2 fluxes to the PDO-SCA correlation is small, as indicated by the weak correlation between PDO and Jena-F of 0.33 (Figure 3b). A relatively weak correlation is expected considering that the changes in land temperature or precipitation/evaporation associated with the PDO are quite heterogeneous in space and time ( Figure S4 in Supporting Information S1).
Our analysis suggests that the controls of winds on the MLO SCA variability were unusually weak during the period from the mid-1980s to the early 1990s. The SCA was low in the late 1980s and higher in the early 1990s. During the 1985-1995 period, the SCA-PDO correlation (Figure 3a, r = −0.18, not smoothed) and the SCA-Jena-W correlation (Figure 1c, r = −0.03, not smoothed) was unusually weak, while the SCA-Jena-F correlation was relatively stronger (Figure 1c, r = 0.45, not smoothed), suggesting that the variability within this period was dominated by fluxes that were not associated with the PDO. A severe drought in North America in 1988 may have reduced the North American contribution to the SCA variability (Schwalm et al., 2012), while the eruption of the Mt. Pinatubo in 1991 may have strengthened summer CO 2 uptake throughout the Northern Hemisphere (Angert et al., 2004;Gu et al., 2003;Jones & Cox, 2005).
The long-term wind shifts might complicate the attribution of the MLO SCA trend to ecological changes. For example, Buermann et al. (2007) suggested that the decrease or stable SCA trend from 1990 to 2003 was mostly due to more frequent droughts over North America, with additional contributions from wind shifts that led to an overall weakening of springtime Eurasia influence. Our results partially agree with Buermann et al. (2007), as Jena-F driven by midlatitude North America fluxes alone ( Figure S3c in Supporting Information S1) shows decreasing SCA anomalies from 1990 to 2000. Over the long-term, however, Jena-F shows North American fluxes contributing to the overall increasing trend in SCA, with this trend persisting after 2000. Also, we do not find a persistent high correlation between PDSI NA and Jena-F, as shown in Figure 3b, but we do find a generally high correlation of PDSI NA with Jena-W since 1980, which is also when PDSI NA show a significant correlation with observed MLO SCA, as in Figure 2. This suggests that the higher PDSI-SCA correlation since the mid-1980s (Figure 2 and in Buermann et al. (2007)) was not driven primarily by CO 2 fluxes tied to drought but rather by wind variability that happens to be correlated with the PDSI NA . On the other hand, the PDSI NA -Jena-F correlation is high after 2008, possibly implying an increased role of hydrological variability after 2008. More generally, Jena-WF shows that the SCA variability is dominated by the Eurasia influence due to both varying winds and varying ecological CO 2 fluxes, with relatively small contributions from North America ( Figure S3b in Supporting Information S1).

Summary and Conclusion
We find that varying winds alone contribute significantly to the trend and variability of MLO SCA.
Since 1959, the Jena-W simulations suggest that varying winds by themselves caused a downward SCA trend of −0.10 ± 0.022 ppm decade −1 , offsetting ∼30% of the long-term SCA increase (0.22 + 0.1 = 0.32 ppm dec ade −1 ) due to ecological CO 2 flux changes alone. The wind-driven decrease simulated by Jena-W accelerated after 1986 to −0.15 ± 0.057 ppm decade −1 , while this decrease simulated by ACTM-W since 1986 is relatively smaller at −0.09 ± 0.068 ppm decade −1 . Both of these two simulations are consistent within uncertainties to Wang et al. (2020), who reported SCA decreasing between −0.1 and −0.2 ppm decade −1 driven by changing circulation alone after 1986. Roughly 44% of the post-1980 downward trend is explained by a secular trend of U-PDO (Figure 4b). We also find a wind-driven decreasing tendency in SCA at other low-latitude stations (MID, KUM, and CHR), with Jena-W suggests a larger decrease than ACTM-W (Table 1). Whereas Wang et al. (2020) speculated that the wind effect may be related to Hadley cell expansion in the Southern Hemisphere, we suggest a possibly important role of poleward expansion in the Northern Hemisphere, which acts to isolate the air at low latitudes from the influence of northern high-latitude ecosystems.
On interannual/decadal time scales, we find a high correlation between detrended MLO SCA and the annual average PDO index, as well as a closely related wind index, U-PDO. This high correlation is closely tied to PDO-related wind patterns: (a) In a positive PDO year (Figure 6a), the Aleutian low is stronger and moves equatorward, resulting in a stronger and equatorward-displaced westerly jet over the Pacific which transports air from Eurasia to MLO. (b) In a negative PDO year (Figure 6b), the Aleutian low shrinks and moves poleward and westward, leading to a weaker and poleward-displaced westerly jet, and leaving MLO more exposed to CO 2 transported from North America and the Eastern Pacific.
Our analyses suggest that long-term circulation changes can partly explain the smaller relative SCA growth rate at low latitudes compared to high latitudes, observed at both surface stations and midtroposphere airborne data, as shown in Graven et al. (2013). Especially at low latitudes, the changing circulation makes it hard to simply relate SCA changes to ecological changes, such as the enhanced magnitude of CO 2 fertilization effect (Wenzel et al., 2016) and North America drought (Buermann et al., 2007).
Wind data that are used to conduct EOF analysis are from NCEP/NCAR reanalysis and downloaded from https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.pressure.html (last access: 10 December 2020; Kalnay et al., 1996). Simulated CO 2 time series data based on Jena CO 2 inversion and TM3 transport model, and MIROC-ACTM are available from https://doi.org/10.5281/zenodo.5514400. Acknowledgments Y. Jin was supported under a grant from NSF (AGS-1623748) and under a grant from Earth Networks. The historical measurements of the Scripps CO 2 program at MLO, KUM, CHR, BRW, have been supported by numerous grants from NSF and DOE and by staff at the Mauna Loa Observatory, which is supported by NOAA. Recent measurements at these stations have been supported by NASA (NNX17AE74G), NSF (OPP-1304270, OPP-1922922) and by Eric and Wendy Schmidt via recommendation of the Schmidt Futures program, with support for the instrumentation at MLO provided by Earth Networks. We also thank two anonymous reviewers and Dr. Benjamin Lintner for their valuable comments and efforts. Any opinions, findings, conclusions, and recommendations expressed in this material are those of the authors and do not necessarily reflect the views of NSF, NOAA, or DOE.