1. Introduction
Climate variability on interdecadal scales has received considerable attention in the past decades. Policymakers and stakeholders would gain significantly from an extensive understanding of interdecadal climate variability because they would be able to formulate long-term plans for agriculture, water resources, environment, and infrastructure promptly. Observational evidence indicates a pronounced decadal–interdecadal variability in the East Asian summer rainfall (EASR) over the twentieth century (Nitta and Hu 1996; Wang et al. 2000; Chan and Zhou 2005; Zhou et al. 2006; Ding 2007). Several studies have revealed that the EASR underwent a striking transition around the late 1970s, featuring a southern flooding–northern drought pattern with more rainfall in the middle and lower reaches of the Yangtze River valley (YRV) and less rainfall in North China (Gong and Ho 2002; Yu et al. 2004; Wang et al. 2013). Concurrently, the East Asian summer monsoon (EASM) weakened in the late 1970s (Wang 2001; Jiang and Wang 2005; Wu et al. 2008). Since the late 1990s, the EASM has intensified (Zhu et al. 2014; Zhang 2015) and the onset has remarkably advanced (Xiang and Wang 2013); consequently, summer rainfall has decreased along the YRV’s middle and lower reaches but increased to the north (Si et al. 2009; Zhu et al. 2013; Xu et al. 2015).
There is considerable evidence suggesting that the interdecadal variability of EASR is closely associated with the Atlantic multidecadal variability (AMV) and atmospheric heat source over the Tibetan Plateau (TPHS). The AMV refers to a variability of the North Atlantic sea surface temperature (SST) with a period of about 60–80 years (Schlesinger and Ramankutty 1994; Kerr 2000), which is often thought to be driven by the Atlantic meridional overturning circulation (Delworth and Mann 2000; Hu et al. 2012). The AMV is also thought to be regulated by volcanic and anthropogenic aerosols during the twentieth century (Otterå et al. 2010; Booth et al. 2012; Qin et al. 2020). Many studies have presented evidence that the AMV can substantially modulate the climate variation in the Northern Hemisphere (Sutton and Hodson 2005; Zhang and Delworth 2006; Knight et al. 2006; Wang et al. 2009; Ting et al. 2011; Sun et al. 2015; Miao and Jiang 2021). In East Asia, a positive AMV was thought to enhance EASM through coupled ocean–atmosphere feedback in the western Pacific and Indian Oceans (Lu et al. 2006). A negative AMV (the cooling North Atlantic), such as that seen in the 1960s, significantly suppressed the EASM (Liu and Chiang 2012). The weakened monsoon circulation induced a long-lasting drought over the North China portion of the EASM region in the 1960s (Liu and Chiang 2012). Since the 1990s, the warming trend in the North Atlantic has caused increased hot temperature extremes over northern East Asia (Dong et al. 2016; Hong et al. 2017; Zhang et al. 2020). Moreover, beyond East Asia, the AMV affects many regions in the Northern Hemisphere through a circumglobal wave train pattern (Si and Ding 2016; Nicolì et al. 2020; Monerie et al. 2021). Recent studies have emphasized the synergistic impacts of the AMV and other oceans on the atmospheric circulation and climate anomalies over East Asia (Yang et al. 2017; Zhang et al. 2018; Huang et al. 2019).
The Tibetan Plateau (TP), with a mean elevation in excess of 4000 m, is located between the North Atlantic–European region and East Asia. The TP serves as an elevated heat source in boreal summer and plays a critical role in the downstream East Asian climate (Tao and Ding 1981; Luo and Yanai 1984; Ye and Wu 1998; Wang et al. 2008; Xu et al. 2010). Modeling study reveals the important role of the TP soil temperature in the lagged relationship between spring TPHS and EASR (Xu et al. 2022). The TPHS is correlated with the EASM and EASR, and a cooling TP since the late 1970s can cause above-normal summer rainfall in the middle and lower reaches of the YRV (Zhang et al. 2004; Wang et al. 2008; Ding et al. 2009). After the late 1990s, the TP has shown a notable warming trend, which intensified the EASM, thereby causing a northward migration of the summer rainfall belt over East Asia (Si and Ding 2013; Zhang et al. 2021). Moreover, the TPHS can cause significant climate fluctuations in other regions of the Northern Hemisphere, such as South Asia (Jiang and Ting 2017), west Asia, North Africa, and south Europe (Lu et al. 2018; Monerie et al. 2019; Sun et al. 2021; Wang et al. 2022), through the atmospheric teleconnection. Besides, the snow cover in the TP, through altering TPHS, exerts a modulating effect on the teleconnection between El Niño–Southern Oscillation and EASM (Wu et al. 2012). Recent studies additionally suggest that zonal teleconnections at midlatitude in the Northern Hemisphere also significantly affect the summer monsoon (Fan et al. 2022) and sensible heat flux (Zhu et al. 2019) over the TP. In particular, the positive phase of Pacific decadal oscillation since the 1980s weakens the sensible heat flux over the TP in the spring through a zonal wave train (Zhu et al. 2019).
The independent effects of the AMV and TPHS on the EASR have been studied extensively in the past; however, no attention has been given to their synergistic impact. Moreover, the relationship between the AMV and TPHS also remains unclear. Hence, we aim to address the relationship between the AMV and TPHS and the role of TP in the impact of AMV on the interdecadal variation of EASR. The rest of this article is organized as follows: Section 2 introduces the datasets, methods, the model, as well as the experiments. A distinctive relationship between the AMV and TPHS is presented in section 3. In section 4, we analyze the role of TP in the impact of AMV on the interdecadal variation of EASR. Finally, a summary is presented in section 5.
2. Data, methods, model, and experiments
a. Data
The observational and reanalysis data used in this study and their sources are as follows: 1) the monthly mean SST data from the Hadley Centre HadISST1.1 dataset on a 1° × 1° grid (Rayner et al. 2003); 2) the monthly mean land precipitation data from the Climatic Research Unit Time Series 3.26 (CRU-TS_3.26), with a horizontal resolution of 0.5° × 0.5° (Harris et al. 2014); and 3) the atmospheric pressure, sensible heat flux, and net radiation-flux data from Twentieth-Century Reanalysis (20CR), version 2, dataset of the National Oceanic and Atmospheric Administration for the period 1901–2012 (Compo et al. 2011) and the NCEP–NCAR reanalysis for the period 1949–2021 (Kalnay et al. 1996). We only focus on the summer [June–July–August (JJA)] season.
b. Method
The AMV index is defined as the linear detrended area-weighted average of the JJA mean North Atlantic (80°–0°W, 0°–60°N) SST anomalies.
To depict the TPHS intensity, we define the linear detrended vertically integrated Q1 averaged over the TP as the TPHS index.
In addition, the singular vector decomposition (SVD) method is used to acquire the most leading coupled mode between the North Atlantic SST and TPHS (Bretherton et al. 1992). The two-tailed Student’s t test is used to measure the significance of the composite analysis results. According to the test, the null hypothesis H0 assumes that the difference between two sample means is zero (supposing that the two samples with Gaussian distribution are independent); conversely, the alternative hypothesis H1 assumes that the difference is not zero. To obtain interdecadal variability of the AMV and TPHS indices, we filter out the variability shorter than 13 years by applying a low-pass Lanczos filter. The statistical significance of the regression and correlation analysis is evaluated by the Monte Carlo method (Livezey and Chen 1983).
c. Model and experiments
A 10-member ensemble of North Atlantic pacemaker simulations performed by the NCAR CESM1 model (Hurrell et al. 2013) is analyzed in this study. The pacemaker simulation’s fully restored time-evolving SST anomalies in the North Atlantic (5°–55°N, with two buffer zones along 0°–5°N and 55°–60°N) are nudged to observations. In this way, the observed variation of SSTs in the North Atlantic is maintained in each simulation, with the rest of the model’s coupled climate system free to evolve. The simulations are branched from 1920 to 2013. Several studies have suggested that the pacemaker simulations can capture the atmospheric response to the observed North Atlantic SST forcing (Yang et al. 2020; Meehl et al. 2021).
3. Distinctive relationship between the AMV and TPHS
a. Modulation effects of the AMV on the TPHS
The time series of the TPHS and AMV indices since the 1900s are shown in Fig. 1, both of which display significant interdecadal variabilities. The two indices are almost in phase during the entire period, with positive phases from the late 1990s to the 2000s and from the mid-1920s to the late 1960s, but negative ones from the late 1960s to the late 1990s and from the late 1910s to the mid-1920s. It appears that the TPHS has a close relationship with the AMV, and their correlation coefficient is 0.70, reaching the 95% confidence level. In addition, we extend the TPHS index to year 2021 by using the NCEP–NCAR reanalysis data. As seen in Fig. 2, the two indices are also in phase from 1949 to 2021, further confirming their in-phase relationship on the interdecadal scale.
We also examine the relationship between the TPHS and the interdecadal SST modes in the Pacific and Indian Oceans, respectively. Figure 3 shows the second leading mode of the empirical orthogonal function (EOF) analysis for the North Pacific and Indian Ocean SSTs. Their first EOF modes indicate a basinwide warming trend (figure not shown). The second EOF mode for the North Pacific and Indian Ocean SST variability corresponds to the Pacific decadal oscillation (PDO) and Indian Ocean dipole (IOD) mode, respectively (Fig. 3). The correlation coefficients between the TPHS index and the PDO and IOD indices are −0.28 and 0.35, respectively, and much lower than that between the TPHS and AMV indices (0.70). This result shows that the TPHS highly correlates with the AMV rather than the PDO and IOD in summer. Hence, this comparison further confirms that the interdecadal variation of the summer TPHS is significantly modulated by the AMV.
The SVD method is used to objectively capture the coupled mode between the North Atlantic (80°W–0°, 0°–60°N) SST and TPHS on the interdecadal scale. Before the analysis, the linear trend within each grid in the SST and TPHS fields is removed. Figure 4 displays the leading SVD mode for the SST and TPHS fields, accounting for 66.5% of the total squared covariance. Hence, this mode represents the coupled relationship between the AMV and TPHS. The results show that on the interdecadal scale, for the North Atlantic SST pattern, there are positive anomalies in most domains with two maximum centers in the low latitudes and subpolar regions but negative over the coasts of northeastern North America and west Europe (Fig. 4a). Meanwhile, the heat source is above normal in most domains of the TP, except for scattered areas over the southern and northern flanks (Fig. 4b). A high correlation coefficient of 0.86 between the two corresponding time coefficients of the two fields shows a close relation between the AMV and TPHS (Fig. 4c). In addition, the field significant test is also used to verify the relationship between the time series of the first SVD (SVD1) for the North Atlantic SST field and TPHS field, based on the method of Livezey and Chen (1983): 43.9% of the entire region exceeds the 95% significance level. The test results show that in less than 1% of the cases, the proportion of the correlation coefficient exceeding the 95% confidence level is greater than 43.9% of the entire region. Hence, it is highly unlikely that the positive relationship between the AMV and TPHS is a chance occurrence.
Next, we investigate the mechanism for the remote modulation of the AMV on the TPHS, which can be distinguished from the regression pattern of 500-hPa geopotential height and wave-activity flux upon the AMV index (Fig. 5a). The AMV, as noted in a previous study (Si and Ding 2016), is associated with the climate fluctuation over the globe via a circumglobal Atlantic–Northern Hemisphere teleconnection. The teleconnection exhibits a distinct circumglobal zonal wave train pattern with seven centers of action situated along the northern midlatitude (Fig. 5a), which originates from the midlatitude North Atlantic, propagates northeastward to western Europe, then extends eastward through the Eurasian continent to the western North Pacific along the zonal westerly (Fig. 5a). The evidence of AMV-associated atmospheric teleconnection has also been identified in several modeling studies. By analyzing idealized AMV experiments, Nicolì et al. (2020) suggested that this AMV-induced circumglobal stationary wave train could affect vast portions of northern Eurasia. The propagation of the Rossby wave follows the subtropical jet stream and regulates the climate from Europe to eastern Asia, which is triggered by the North Atlantic SST and is associated with the AMV (Monerie et al. 2021). Ruprich-Robert et al. (2017) distinguished different roles of the tropical and subpolar portions of the AMV in modulating the circumglobal stationary wave train. Since the nineteenth century, amplitude of the AMV has experienced an increasing trend under the anthropogenic forcings, which will intensify this zonal teleconnection (Moore et al. 2017).
In the midtroposphere, anomalous low pressure dominates the northern South Asia–TP–East Asia sector, while anomalous high pressure appears over the surrounding oceans, such as the western North Pacific and Indian Ocean (Fig. 5a). This “land low and ocean high” pressure pattern reinforces the Asian summer monsoon (Fig. 6). Moreover, a positive phase of AMV leads to a warming over the TP–Asia sector and enhanced land–sea thermal contrast between the Asian continent and the surrounding oceans through this wave train, which also strengthens the Asian summer monsoon. The southwesterly vapor transport associated with the reinforced South Asian summer monsoon (SASM) and the southeasterly vapor transport associated with the reinforced EASM strengthen the vapor transport to the TP, together with the anomalous low pressure over the TP, enhancing the summer precipitation and latent heating over the TP (Fig. 6).
To further explore the teleconnection between the AMV and TPHS, we examine the vertical structure of the atmospheric teleconnection pattern (Fig. 5b). The four centers of action from the North Atlantic to the TP reach their maximum strength in the middle and upper troposphere, which are accompanied by prominent eastward dispersion of the Rossby wave originating from the North Atlantic Ocean. Notably, the teleconnection pattern exhibits an equivalent barotropic structure in the vertical over the upstream regions of the TP, while it exhibits a baroclinic structure over the TP, with an anomalous low pressure in the lower levels accompanied by an anomalous high pressure in the middle and upper levels. This suggests that the warming of the North Atlantic (the positive phase of the AMV) can strengthen the baroclinicity over the TP, and vice versa.
Moreover, the AMV exerts a significant remote influence on the atmospheric heat source Q1 from the North Atlantic to the TP through the wave train pattern. In the Atlantic–west Europe sector, the negative center of action and North Atlantic warming enhance the convection and latent heating release, causing positive heat-source anomalies there (Fig. 5b). The positive center of action corresponds to an outflow of moisture toward the center and decreases precipitation and latent heating release, causing negative heat-source anomalies over east Europe and central Asia. In addition, a strengthened baroclinicity increases summer rainfall and latent heating (Fig. 6), yielding positive TPHS (Fig. 5b).
b. Booster effect of TP on the influence of AMV on the TPHS
Notably, the positive heat-source anomalies over the TP are stronger than those over other regions along the same latitudinal cross section (Fig. 5b). The boost in the TPHS is mainly due to the orographic effect of TP on the ASM latent release. On the one hand, the strong southwesterly SASM airflow driven by the positive phase of AMV, as described above, starts to split into two branches when reaching the South Asian subcontinent (Fig. 6). One branch flows through northwest India toward the southwestern flank of the TP. The other branch flows through the Bay of Bengal toward the southeastern flank of the TP. On the other hand, the strong southeasterly EASM airflow flows through Southwest China toward the eastern flank of the TP. When the three branches of warm and moist monsoon airflow approach the TP, they converge and climb along the slopes, resulting in an enormous ascending motion. Consequently, heavy precipitation occurs over the TP, particularly over the southern and eastern flanks (Fig. 6). The heavy precipitation eventually generates massive latent heating released in convection and condensation over the TP, with three maximum centers over the southern and eastern flanks (Fig. 7). Hence, the TPHS driven by the AMV is stronger than that over other regions along the same latitude, due to the booster effect of the TP (Fig. 5b). To quantify these results, the latent heat flux over the above three regions is investigated by analyzing four scenarios. The four scenarios are selected based on the AMV index values (strong scenarios: greater or less than one standard deviation; moderate scenarios: between 0.25 and one standard deviation). In general, the latent heating decreases in east Europe while it increases in the North Atlantic–west Europe sector and TP along with the increasing of the AMV index (Fig. 8). In the three regions, the latent heating anomalies are much larger in the case of strong scenarios compared to those in the case of moderate scenarios. It is also noteworthy that, in the case of strong scenarios, the latent heat anomalies in the TP are about 2.3 times and 1.25 times higher than those in North Atlantic–west Europe sector and east Europe, respectively.
From Eq. (3), the vertically integrated atmospheric heat source comprises three components: sensible heat flux, vertically integrated latent heat flux, and net radiation heat flux. We estimate the three components from the complete atmospheric heat-source equation [Eq. (3)] to identify the contribution to changes in the TPHS associated with the AMV. As shown in Fig. 7, both sensible heat flux and net radiation heat flux make small and negative contributions to the enhancement of TPHS. The magnitude and distribution of the latent heat flux resemble those of the actual heat source anomalies over the TP associated with the AMV. Moreover, the high coefficient is observed between the filtered time-series latent heat flux over the TP and the TPHS index (r = 0.77, significant at the 95% confidence level) for the period 1901–2012. Hence, the latent heat flux due to convection and precipitation plays a predominant role in the enhancement of TPHS modulated by the AMV. Previous study has found that the latent heating associated with precipitation is the primary factor in the TPHS in summer (Ye and Gao 1979). The results also indicate that the positive phase of the AMV can enhance the latent heating over the TP, resulting in an enhanced heat source there through the zonal teleconnection pattern, and vice versa.
c. Model results
To further validate the modulation effect of the AMV on the TPHS, the North Atlantic pacemaker ensemble simulations are analyzed in this section. Because the SSTs over the North Atlantic are nudged toward the observed anomalies, the simulated AMV index by the North Atlantic pacemaker simulations matches well with the observed one during 1920–2013 (figure not shown). The years in which the ensemble mean AMV index is below −0.25 and above +0.25 standard deviation are set as the negative (43 years) and positive (39 years) phases of the AMV, respectively. The composite simulated North Atlantic SST anomalies and the heat-source anomalies over the TP during the periods of the positive and negative phases and their difference are shown in Fig. 9. When the AMV is in the negative phase, the North Atlantic SST anomalies show a meridional “cold–mild–cold” tripole pattern with a significant maximum value of negative SST anomalies located in the subpolar gyre region (Fig. 9a). As for the TPHS, except for the weak positive anomalies over the northwest part of the TP, most parts of the TP have negative anomalies with significant maximum values over the southern and eastern flanks (Fig. 9b). In general, the SST and heat-source anomalies when the AMV is in the positive phase are mirror images of the negative phase of AMV with opposite signs, indicating that the AMV, indeed, is the primary cause of these anomalies (Figs. 9c,d). The differences between the positive and negative phases of the AMV (Figs. 9e,f) confirm the modulation effect of the AMV on the TPHS.
4. Synergistic influence of the AMV and TPHS on the EASR
The TPHS exerts a significant impact on the summer rainfall and atmospheric circulation over the TP and its downstream regions, including East Asia. Locally, the TP warming (the positive phase of the TPHS) significantly increases summer precipitation over the TP, particularly along its southern flank (Fig. 10a). In East Asia, the precipitation anomalies show a tripole pattern with two wet anomalies along the Yellow–Huaihe River valley to the Korean Peninsula and South China, while there is one dry anomaly along the YRV to southern Japan, which covers the prevailing mei-yu–baiu–changma frontal area (Fig. 10a). The mechanism of the TPHS influences on the EASM and rainfall anomalies is found to be through a distinct Rossby wave train. The wave train moves along the low-level SASM southwesterly that is excited in the southwest periphery of the TP (Figs. 10b,c). The development of the wave train involves a Rossby wave energy propagation in SASM regions (Fig. 10b). Because of the strong vertical shear in the southwesterly monsoon region, the Rossby wave energy propagation is trapped in the lower troposphere, which generates a stationary Rossby wave train (Wang and Xie 1996). The eastward propagation of this wave train induces an anticyclonic anomaly in the South China Sea (Fig. 10c) then propagates northward into East Asia. The anticyclonic anomaly over the South China Sea enhances the southwesterly along the south coast of China, increasing the warm and moist air float transport toward South China and summer rainfall there. In northeast Asia, the anomalous anticyclone converges the northward water vapor and southward cold air to the north of the YRV, enhancing the summer rainfall along the Yellow–Huaihe River valley to the Korean Peninsula and decreasing the summer rainfall along the lower reach of the YRV and southern Japan (Fig. 10a). Our results agree with a previous modeling study (Wang et al. 2008) that warming over the TP induces a distinct Rossby wave train propagated along the SASM monsoon southwesterly that is excited in the southwest periphery of the TP. This wave train propagates into East Asia and leads to an enhanced EASM. These analyses identify that the TPHS exerts a striking impact on the EASM on the interdecadal scale, with an enhanced TPHS causing a strong EASM, and vice versa.
The analyses in section 3 have identified that the AMV significantly modulates the variations of TPSH. Notably, through the zonal teleconnection pattern, the AMV can also affect the climate over regions downstream of the TP. Over the East Asia–western North Pacific regions, the meridionally arranged wave pattern is similar to that associated with the TPHS (Fig. 6). The anomalous anticyclone over the South China Sea and cyclone to the north of the YRV enhance the EASM, resulting in the northward shift of the summer rain belt in East Asia with dry anomaly along the YRV to southern Japan but wet anomaly along the Yellow–Huaihe River valley to the Korean Peninsula (Fig. 6). These results agree with previous observational (Si and Ding 2016) and numerical (Liu and Chiang 2012) studies, suggesting that the AMV can cause interdecadal changes in the EASM and EASR through the circumglobal teleconnection. For the EASR, although the impact of the AMV is similar to that of the TPHS, there are still differences between them. Compared with the AMV, the TPHS exerts a stronger impact on the rainfall over South China but a weaker impact on the rainfall over the YRV (Figs. 6 and 10a). Moreover, we also perform the field significance test (Livezey and Chen 1983) on the relationship among the AMV, TPHS index, and the summer rainfall in the East Asian monsoon region (10°–45°N, 100°–130°E), respectively. It shows that 54.3% and 49.4% of the entire region exceed the 90% significance level, respectively. The test results suggest that the correlation between the AMV, TPHS index, and the summer rainfall over the East Asian monsoon region exceeds the field significance test at the 99% confidence level.
The influences of the AMV and TPHS on the East Asian climate have been investigated separately by most previous studies, with little attention on their synergistic effect. To address this issue on the interdecadal scale, composite JJA precipitation anomalies are constructed for four combinations of the AMV and TPHS indices below and above a ±0.25 standard deviation (Fig. 1). The four combinations are classified as follows: 1) positive AMV (+AMV) and positive TPHS (+TPHS), 2) negative AMV (−AMV) and negative TPHS (−TPHS), 3) +AMV and −TPHS, and 4) −AMV and +TPHS.
For +AMV, the precipitation anomalies over East Asia generally show a tripole pattern with wet anomalies along the Yellow–Huaihe River valley to the northern Korean Peninsula and South China but dry anomalies along the YRV (Figs. 11a,c). The wet anomalies over South China are enhanced during +TPHS (Fig. 11a) but weakened during −TPHS (Fig. 11c). Meanwhile, for −AMV, the precipitation anomalies over East Asia also show tripole patterns but are of opposite signs (Figs. 11b,d). The dry anomalies over South China are enhanced during −TPHS (Fig. 11b) but weakened during +TPHS (Fig. 11d). Under −AMV and −TPHS, the rainfall anomalies in East Asia are more reminiscent of the southern flooding–northern drought in the 1970s–90s, with excessive rainfall along the middle and lower basins of the YRV but deficient to the north of the YRV (Fig. 11b). Under +AMV and +TPHS, the anomalous pattern is more reminiscent of the northward shift of the summer East Asian rain belt after the late 1990s (Si et al. 2009; Si and Ding 2013; Zhu et al. 2013; Xu et al. 2015), with deficient rainfall in the middle and lower basins of the YRV and excessive rainfall to the north of the YRV (Fig. 11a).
To explore the atmospheric process through which the AMV and TPHS synergistically impact the EASR, the composite atmospheric circulation anomalies under the four combinations of the AMV and TPHS phases are analyzed. For +AMV, the associated teleconnection wave train is strong during +TPHS (Fig. 12a) but weak during −TPHS (Fig. 12c). Under +AMV and +TPHS, the +TPHS corresponds to a warming over the TP, and +AMV also leads to a warming over the TP and Asia continent through the wave train (Si et al. 2023). Hence, the remote and local heating superimpose over the TP and the surrounding region, which reinforces the wave train over the TP and downstream over the Asia–Pacific region. The TPHS boosts the impacts of the AMV on the downstream climate, which jointly lead to a “high–low–high”-like pattern in the East Asia–western North Pacific regions (Fig. 12a). The anomalous high pressure over the western North Pacific intensifies the subtropical high there (WNPSH). In the lower troposphere, a meridional anomalous “anticyclone–cyclone–anticyclone”-like wave train extends from East Asia to the western North Pacific (Fig. 13a). The anomalous anticyclonic circulation over the western North Pacific and strong WNPSH significantly enhance the EASM, resulting in wet anomalies along the Yellow–Huaihe River valley to the northern Korean Peninsula and South China but dry anomalies along the YRV (Fig. 11a). Under +AMV and −TPHS, the heating over the TP is significantly suppressed, which weakens the wave train over the TP and downstream over the Asia–Pacific regions (Fig. 12c). The anomalous high pressures are located over the TP and western North Pacific, weakening the TPHS but enhancing the WNPSH. Although the EASM is slightly strong (Fig. 13c), the suppressed TPHS weakens the wet anomalies over South China (Fig. 11c).
For −AMV, the associated teleconnection pattern shifts to its opposite sign (Figs. 12b,d). Similarly, the teleconnection wave train is strong during −TPHS (Fig. 12b) but weak during +TPHS (Fig. 12d). Under −AMV and −TPHS, the wave train is also reinforced over the TP and downstream over the Asia–Pacific regions but becomes roughly opposite in sign to that under +AMO and +TPHS. The strong wave train reaches the Asia–Pacific regions, which leads to a “low–high–low”-like pattern in the East Asia–western North Pacific sector (Fig. 12b). Their synergistic impacts significantly weaken the WNPSH and EASM (Fig. 13b), causing more rainfall over the YRV, Korean Peninsula, and southern Japan but less over the Yellow–Huaihe River valley and Northeast China (Fig. 11b). Under −AMV and +TPHS, the wave train weakens, and the anomalous low pressure over the TP intensifies the latent heating and heat source there (Fig. 12d). The low pressure anomaly over the western North Pacific corresponds to a weakened WNPSH, which weakens the EASM (Fig. 13d). Moreover, the enhanced TPHS weakens the dry anomalies over South China (Fig. 11d).
5. Summary
Previous studies have shown that the EASR exhibits prominent interdecadal variability. Although the key roles of the AMV and TPHS in the interdecadal variability of the EASR have been identified in many previous studies, the studies on the relationship between the AMV and TPHS, and the role of TP in the impact of AMV on the EASR are inadequate.
In this study, we reveal a distinctive relationship between the AMV and TPHS, which is illustrated by the schematic diagram in Fig. 14. Under +AMV, the North Atlantic warming induces a zonal teleconnection wave train originating from the North Atlantic through Eurasia and extending to the North Pacific. This teleconnection pattern exhibits an equivalent barotropic structure in the vertical over the upstream regions of the TP and a baroclinic structure over the TP. The reinforced baroclinicity increases summer rainfall and latent heating over the TP, yielding an enhanced TPHS. Moreover, the +AMV enhances SASM and EASM. The enhanced warm and moist monsoon airflows approach the TP, converge, and climb up along its southern and eastern flanks and result in enormous latent heating release in convection and condensation over the TP. Eventually, heat-source anomalies over the TP are stronger than those over other regions along the same latitude, due to the orographic effect of the huge plateau. Hence, the AMV exerts a remote effect on the TPHS through the zonal teleconnection wave train, whereas the TP serves as a booster of this effect.
The +AMV can enhance the EASM via the zonal teleconnection pattern, yielding a tripole rainfall pattern with one dry anomaly along the YRV to southern Japan but two wet anomalies over the Yellow–Huaihe River valley–Korean Peninsula regions and South China, and vice versa. The enhanced TPHS modulated by +AMV intensifies the EASM via a meridional teleconnection wave pattern, yielding a tripole rainfall pattern, which is similar to that caused by the +AMV. Meanwhile, the TPHS exerts a stronger impact on the summer rainfall over South China but a weaker impact on that over the YRV compared with the AMV.
When the AMV and TPHS are in phase, the zonal wave train originating from the North Atlantic is strong and propagates eastward to the Asia–Pacific regions, significantly modulating not only the TPHS but also the downstream EASM. Notably, the TP boosts this remote impact from the AMV. Positive phases of the AMV and TPHS jointly result in a significantly strong EASM through the effect of superimposition. This situation is more reminiscent of the northward shift of the East Asian rain belt since the late 1990s. When the reversals in phases of the AMV and TPHS occur, that is, −AMV and −TPHS, there is a significantly weak EASM, which is more reminiscent of the southern flooding–northern drought in the 1970s–90s in East Asia. Meanwhile, when the AMV and TPHS are out of phase, the zonal teleconnection wave train is weak, particularly over Asia and its surrounding regions, contributing to a slightly weak or strong EASM through the effects of cancelation.
Acknowledgments.
The authors thank the three anonymous reviewers for their constructive comments and for significantly improving the manuscript. Our research is jointly supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant XDA20100304), the National Natural Science Foundation of China (Grant 41875104), and the Second Tibetan Plateau Scientific Expedition and Research Program of China (Grant 2019QZKK0208).
Data availability statement.
All data used in this study are publicly available. The HadISST1.1 SST data are openly available at Met Office Hadley Centre (https://www.metoffice.gov.uk/hadobs/hadisst/). The CRU-TS_3.26 precipitation data are openly available at Climatic Research Unit (CRU) at the University of East Anglia (https://research-portal.uea.ac.uk). The NOAA-20CR and NCEP/NCAR reanalysis data are openly available at National Oceanic and Atmospheric Administration (https://psl.noaa.gov/data/index.html). The CESM1 North Atlantic pacemaker simulations are accessible via the National Center for Atmospheric Research (NCAR) Climate Data Gateway (https://www.earthsystemgrid.org/).
REFERENCES
Booth, B. B. B., N. J. Dunstone, P. R. Halloran, T. Andrews, and N. Bellouin, 2012: Aerosols implicated as a prime driver of twentieth-century North Atlantic climate variability. Nature, 484, 228–232, https://doi.org/10.1038/nature10946.
Bretherton, C. S., C. Smith, and J. M. Wallace, 1992: An intercomparison of methods for finding coupled patterns in climate data. J. Climate, 5, 541–560, https://doi.org/10.1175/1520-0442(1992)005<0541:AIOMFF>2.0.CO;2.
Chan, J. C. L., and W. Zhou, 2005: PDO, ENSO, and the early summer monsoon rainfall over South China. Geophys. Res. Lett., 32, L08810, https://doi.org/10.1029/2004GL022015.
Compo, G. P., and Coauthors, 2011: The Twentieth Century Reanalysis Project. Quart. J. Roy. Meteor. Soc., 137 (654), 1–28, https://doi.org/10.1002/qj.776.
Delworth, T. L., and M. E. Mann, 2000: Observed and simulated multidecadal variability in the Northern Hemisphere. Climate Dyn., 16, 661–676, https://doi.org/10.1007/s003820000075.
Ding, Y., 2007: The variability of the Asian summer monsoon. J. Meteor. Soc. Japan, 85B, 21–54, https://doi.org/10.2151/jmsj.85B.21.
Ding, Y., Y. Sun, Z. Wang, Y. Zhu, and Y. Song, 2009: Inter-decadal variation of the summer precipitation in China and its association with decreasing Asian summer monsoon Part II: Possible causes. Int. J. Climatol., 29, 1926–1944, https://doi.org/10.1002/joc.1759.
Dong, B., R. T. Sutton, W. Chen, X. Liu, R. Lu, and Y. Sun, 2016: Abrupt summer warming and changes in temperature extremes over Northeast Asia since the mid-1990s: Drivers and physical processes. Adv. Atmos. Sci., 33, 1005–1023, https://doi.org/10.1007/s00376-016-5247-3.
Fan, W., Z. Hu, W. Ma, Y. Ma, Y. Yang, H. Yu, and X. Han, 2022: Impacts of mid-high latitude atmospheric teleconnection patterns on interannual variation of the Tibetan Plateau summer monsoon. Atmos. Res., 275, 106219, https://doi.org/10.1016/j.atmosres.2022.106219.
Gong, D.-Y., and C.-H. Ho, 2002: Shift in the summer rainfall over the Yangtze River valley in the late 1970s. Geophys. Res. Lett., 29, 1436, https://doi.org/10.1029/2001GL014523.
Harris, I., P. D. Jones, T. J. Osborn, and D. H. Lister, 2014: Updated high-resolution grids of monthly climatic observations—The CRU TS3.10 dataset. Int. J. Climatol., 34, 623–642, https://doi.org/10.1002/joc.3711.
Hong, X., R. Lu, and S. Li, 2017: Amplified summer warming in Europe–west Asia and northeast Asia after the mid-1990s. Environ. Res. Lett., 12, 094007, https://doi.org/10.1088/1748-9326/aa7909.
Hu, A., G. A. Meehl, W. Han, A. Abe-Ouchi, C. Morrill, Y. Okazaki, and M. O. Chikamoto, 2012: The Pacific-Atlantic seesaw and the Bering Strait. Geophys. Res. Lett., 39, L03702, https://doi.org/10.1029/2011GL050567.
Huang, D., A. Dai, B. Yang, P. Yan, J. Zhu, and Y. Zhang, 2019: Contributions of different combinations of the IPO and AMO to recent changes in winter East Asian jets. J. Climate, 32, 1607–1626, https://doi.org/10.1175/JCLI-D-18-0218.1.
Hurrell, J. W., and Coauthors, 2013: The Community Earth System Model: A framework for collaborative research. Bull. Amer. Meteor. Soc., 94, 1339–1360, https://doi.org/10.1175/BAMS-D-12-00121.1.
Jiang, D., and H. Wang, 2005: Natural interdecadal weakening of East Asian summer monsoon in the late 20th century. Chin. Sci. Bull., 50, 1923–1929, https://doi.org/10.1360/982005-36.
Jiang, X., and M. Ting, 2017: A dipole pattern of summertime rainfall across the Indian subcontinent and the Tibetan Plateau. J. Climate, 30, 9607–9620, https://doi.org/10.1175/JCLI-D-16-0914.1.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–472, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.
Kerr, R. A., 2000: A North Atlantic climate pacemaker for the centuries. Science, 288, 1984–1985, https://doi.org/10.1126/science.288.5473.1984.
Knight, J. R., C. K. Folland, and A. A. Scaife, 2006: Climate impacts of the Atlantic multidecadal oscillation. Geophys. Res. Lett., 33, L17706, https://doi.org/10.1029/2006GL026242.
Liu, Y., and J. C. H. Chiang, 2012: Coordinated abrupt weakening of the Eurasian and North African monsoons in the 1960s and links to extratropical North Atlantic cooling. J. Climate, 25, 3532–3548, https://doi.org/10.1175/JCLI-D-11-00219.1.
Livezey, R. E., and W. Y. Chen, 1983: Statistical field significance and its determination by Monte Carlo techniques. Mon. Wea. Rev., 111, 46–59, https://doi.org/10.1175/1520-0493(1983)111<0046:SFSAID>2.0.CO;2.
Lu, M., S. Yang, Z. Li, B. He, S. He, and Z. Wang, 2018: Possible effect of the Tibetan Plateau on the “upstream” climate over west Asia, North Africa, south Europe and the North Atlantic. Climate Dyn., 51, 1485–1498, https://doi.org/10.1007/s00382-017-3966-5.
Lu, R., B. Dong, and H. Ding, 2006: Impact of the Atlantic multidecadal oscillation on the Asian summer monsoon. Geophys. Res. Lett., 33, L24701, https://doi.org/10.1029/2006GL027655.
Luo, H., and M. Yanai, 1984: The large-scale circulation and heat sources over the Tibetan Plateau and surrounding areas during the early summer of 1979. Part II: Heat and moisture budgets. Mon. Wea. Rev., 112, 966–989, https://doi.org/10.1175/1520-0493(1984)112<0966:TLSCAH>2.0.CO;2.
Meehl, G. A., and Coauthors, 2021: Atlantic and Pacific tropics connected by mutually interactive decadal-timescale processes. Nat. Geosci., 14, 36–42, https://doi.org/10.1038/s41561-020-00669-x.
Miao, J., and D. Jiang, 2021: Multidecadal variations in the East Asian winter monsoon and their relationship with the Atlantic multidecadal oscillation since 1850. J. Climate, 34, 7525–7539, https://doi.org/10.1175/JCLI-D-21-0073.1.
Monerie, P.-A., J. Robson, B. Dong, D. L. R. Hodson, and N. P. Klingaman, 2019: Effect of the Atlantic multidecadal variability on the global monsoon. Geophys. Res. Lett., 46, 1765–1775, https://doi.org/10.1029/2018GL080903.
Monerie, P.-A., J. Robson, B. Dong, and D. L. R. Hodson, 2021: Role of the Atlantic multidecadal variability in modulating East Asian climate. Climate Dyn., 56, 381–398, https://doi.org/10.1007/s00382-020-05477-y.
Moore, G. W. K., J. Halfar, H. Majeed, W. Adey, and A. Kronz, 2017: Amplification of the Atlantic multidecadal oscillation associated with the onset of the industrial-ear warming. Sci. Rep., 7, 40861, https://doi.org/10.1038/srep40861.
Nicolì, D., A. Bellucci, D. Iovino, P. Ruggieri, and S. Gualdi, 2020: The impact of the AMV on Eurasian summer hydrological cycle. Sci. Rep., 10, 14444, https://doi.org/10.1038/s41598-020-71464-2.
Nitta, T., and Z.-Z. Hu, 1996: Summer climate variability in China and its association with 500 hPa height and tropical convection. J. Meteor. Soc. Japan, 74, 425–445, https://doi.org/10.2151/jmsj1965.74.4_425.
Otterå, O. H., M. Bentsen, H. Drange, and L. Suo, 2010: External forcing as a metronome for Atlantic multidecadal variability. Nat. Geosci., 3, 688–694, https://doi.org/10.1038/ngeo955.
Qin, M., A. Dai, and W. Hua, 2020: Quantifying contributions of internal variability and external forcing to Atlantic multidecadal variability since 1870. Geophys. Res. Lett., 47, e2020GL089504, https://doi.org/10.1029/2020GL089504.
Rayner, N. A., D. E. Parker, E. B. Horton, C. K. Folland, L. V. Alexander, D. P. Rowell, E. C. Kent, and A. Kaplan, 2003: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. J. Geophys. Res., 108, 4407, https://doi.org/10.1029/2002JD002670.
Ruprich-Robert, Y., R. Msadek, F. Castruccio, S. Yeager, T. Delworth, and G. Danabasoglu, 2017: Assessing the climate impacts of the observed Atlantic multidecadal variability using the GFDL CM2.1 and NCAR CESM1 global coupled models. J. Climate, 30, 2785–2810, https://doi.org/10.1175/JCLI-D-16-0127.1.
Schlesinger, M. E., and N. Ramankutty, 1994: An oscillation in the global climate system of period 65–70 years. Nature, 367, 723–726, https://doi.org/10.1038/367723a0.
Si, D., and Y. Ding, 2013: Decadal change in the correlation pattern between the Tibetan Plateau winter snow and the East Asian summer precipitation during 1979–2011. J. Climate, 26, 7622–7634, https://doi.org/10.1175/JCLI-D-12-00587.1.
Si, D., and Y. Ding, 2016: Oceanic forcings of the interdecadal variability in East Asian summer rainfall. J. Climate, 29, 7633–7649, https://doi.org/10.1175/JCLI-D-15-0792.1.
Si, D., Y. Ding, and Y. Liu, 2009: Decadal northward shift of the Meiyu belt and the possible cause. Chin. Sci. Bull., 54, 4742–4748, https://doi.org/10.1007/s11434-009-0385-y.
Si, D., A. Hu, D. Jiang, and X. Lang, 2023: Atmospheric teleconnection associated with the Atlantic multidecadal variability in summer: Assessment of the CESM1 model. Climate Dyn., 60, 1043–1060, https://doi.org/10.1007/s00382-022-06331-z.
Sun, C., J. Li, and S. Zhao, 2015: Remote influence of Atlantic multidecadal variability on Siberian warm season precipitation. Sci. Rep., 5, 16853, https://doi.org/10.1038/srep16853.
Sun, X., Y. Ding, and Q. Li, 2021: Interdecadal variation of the atmospheric heat source over the Tibetan Plateau and surrounding Asian monsoon region: Impact on the Northern Hemisphere summer circulation. J. Meteor. Res., 35, 238–257, https://doi.org/10.1007/s13351-021-0101-7.
Sutton, R. T., and D. L. R. Hodson, 2005: Atlantic Ocean forcing of North American and European summer climate. Science, 309, 115–118, https://doi.org/10.1126/science.1109496.
Takaya, K., and H. Nakamura, 2001: A formulation of a phase-independent wave-activity flux for stationary and migratory quasi-geostrophic eddies on a zonally varying basic flow. J. Atmos. Sci., 58, 608–627, https://doi.org/10.1175/1520-0469(2001)058<0608:AFOAPI>2.0.CO;2.
Tao, S.-Y., and Y.-H. Ding, 1981: Observational evidence of the influence of the Qinghai-Xizang (Tibet) Plateau on the occurrence of heavy rain and severe convective storms in China. Bull. Amer. Meteor. Soc., 62, 23–30, https://doi.org/10.1175/1520-0477(1981)062<0023:OEOTIO>2.0.CO;2.
Ting, M., Y. Kushnir, R. Seager, and C. Li, 2011: Robust features of Atlantic multi-decadal variability and its climate impacts. Geophys. Res. Lett., 38, L17705, https://doi.org/10.1029/2011GL048712.
Wang, B., and X. Xie, 1996: Low-frequency equatorial waves in vertically sheared zonal flow. Part I: Stable waves. J. Atmos. Sci., 53, 449–467, https://doi.org/10.1175/1520-0469(1996)053<0449:LFEWIV>2.0.CO;2.
Wang, B., Q. Bao, B. Hoskins, G. Wu, and Y. Liu, 2008: Tibetan Plateau warming and precipitation changes in East Asia. Geophys. Res. Lett., 35, L14702, https://doi.org/10.1029/2008GL034330.
Wang, H., 2001: The weakening of the Asian monsoon circulation after the end of 1970’s. Adv. Atmos. Sci., 18, 376–386, https://doi.org/10.1007/BF02919316.
Wang, S., D. Gong, J. Yi, and Z. Chen, 2000: Seasonal precipitation series of eastern China since 1880 and the variability (in Chinese). Acta Geogr. Sin., 55, 281–293, https://doi.org/10.11821/xb200003004.
Wang, T., H. J. Wang, O. H. Otterå, Y. Q. Gao, L. L. Suo, T. Furevik, and L. Yu, 2013: Anthropogenic agent implicated as a prime driver of shift in precipitation in eastern China in the late 1970s. Atmos. Chem. Phys., 13, 12 433–12 450, https://doi.org/10.5194/acp-13-12433-2013.
Wang, Y., S. Li, and D. Luo, 2009: Seasonal response of Asian monsoonal climate to the Atlantic multidecadal oscillation. J. Geophys. Res., 114, D02112, https://doi.org/10.1029/2008JD010929.
Wang, Z., S. Yang, H. Luo, and J. Li, 2022: Drying tendency over the southern slope of the Tibetan Plateau in recent decades: Role of a CGT-like atmospheric change. Climate Dyn., 59, 2801–2813, https://doi.org/10.1007/s00382-022-06262-9.
Wu, B., R. Zhang, Y. Ding, and R. D’Arrigo, 2008: Distinct modes of the East Asian summer monsoon. J. Climate, 21, 1122–1138, https://doi.org/10.1175/2007JCLI1592.1.
Wu, Z., J. Li, Z. Jiang, and T. Ma, 2012: Modulation of the Tibetan Plateau snow cover on the ENSO teleconnections: From the East Asian summer monsoon perspective. J. Climate, 25, 2481–2489, https://doi.org/10.1175/JCLI-D-11-00135.1.
Xiang, B., and B. Wang, 2013: Mechanisms for the advanced Asian summer monsoon onset since the mid-to-late 1990s. J. Climate, 26, 1993–2009, https://doi.org/10.1175/JCLI-D-12-00445.1.
Xu, H., X.-Z. Liang, and Y. Xue, 2022: Regional climate modeling to understand Tibetan heating remote impacts on East China precipitation. Climate Dyn., 58, 2807–2824, https://doi.org/10.1007/s00382-022-06266-5.
Xu, X., C. Lu, X. Shi, and Y. Ding, 2010: Large-scale topography of China: A factor for the seasonal progression of the Meiyu rainband? J. Geophys. Res., 115, D02110, https://doi.org/10.1029/2009JD012444.
Xu, Z., K. Fan, and H. Wang, 2015: Decadal variation of summer precipitation over China and associated atmospheric circulation after the late 1990s. J. Climate, 28, 4086–4106, https://doi.org/10.1175/JCLI-D-14-00464.1.
Yanai, M., and T. Tomita, 1998: Seasonal and interannual variability of atmospheric heat sources and moisture sinks as determined from NCEP–NCAR reanalysis. J. Climate, 11, 463–482, https://doi.org/10.1175/1520-0442(1998)011<0463:SAIVOA>2.0.CO;2.
Yang, D., J. M. Arblaster, G. A. Meehl, M. H. England, E.-P. Lim, S. Bates, and N. Rosenbloom, 2020: Role of tropical variability in driving decadal shifts in the Southern Hemisphere summertime eddy-driven jet. J. Climate, 33, 5445–5463, https://doi.org/10.1175/JCLI-D-19-0604.1.
Yang, Q., Z. Ma, X. Fan, Z.-L. Yang, Z. Xu, and P. Wu, 2017: Decadal modulation of precipitation patterns over eastern China by sea surface temperature anomalies. J. Climate, 30, 7017–7033, https://doi.org/10.1175/JCLI-D-16-0793.1.
Ye, D., and Y. Gao, 1979: Meteorology of the Qinghai-Xizang Plateau (in Chinese). Science Press, 278 pp.
Ye, D.-Z., and G.-X. Wu, 1998: The role of the heat source of the Tibetan Plateau in the general circulation. Meteor. Atmos. Phys., 67, 181–198, https://doi.org/10.1007/BF01277509.
Yu, R., B. Wang, and T. Zhou, 2004: Tropospheric cooling and summer monsoon weakening trend over East Asia. Geophys. Res. Lett., 31, L22212, https://doi.org/10.1029/2004GL021270.
Zhang, C., Y. Guo, and Z. Wen, 2021: Interdecadal change in the effect of Tibetan Plateau snow cover on spring precipitation over eastern China around the early 1990s. Climate Dyn., 58, 2807–2824, https://doi.org/10.1007/s00382-021-06035-w.
Zhang, J., Z. Chen, H. Chen, Q. Ma, and A. Teshome, 2020: North Atlantic multidecadal variability enhancing decadal extratropical extremes in boreal late summer in the early twenty-first century. J. Climate, 33, 6047–6064, https://doi.org/10.1175/JCLI-D-19-0536.1.
Zhang, R., 2015: Changes in East Asian summer monsoon and summer rainfall over eastern China during recent decades. Sci. Bull., 60, 1222–1224, https://doi.org/10.1007/s11434-015-0824-x.
Zhang, R., and T. L. Delworth, 2006: Impact of Atlantic multidecadal oscillations on India/Sahel rainfall and Atlantic hurricanes. Geophys. Res. Lett., 33, L17712, https://doi.org/10.1029/2006GL026267.
Zhang, Y., T. Li, and B. Wang, 2004: Decadal change of the spring snow depth over the Tibetan Plateau: The associated circulation and influence on the East Asian summer monsoon. J. Climate, 17, 2780–2793, https://doi.org/10.1175/1520-0442(2004)017<2780:DCOTSS>2.0.CO;2.
Zhang, Z., X. Sun, and X.-Q. Yang, 2018: Understanding the interdecadal variability of East Asian summer monsoon precipitation: Joint influence of three oceanic signals. J. Climate, 31, 5485–5506, https://doi.org/10.1175/JCLI-D-17-0657.1.
Zhou, W., C. Li, and J. C. L. Chan, 2006: The interdecadal variations of the summer monsoon rainfall over South China. Meteor. Atmos. Phys., 93, 165–175, https://doi.org/10.1007/s00703-006-0184-9.
Zhu, J., D.-Q. Huang, Y.-C. Zhang, A.-N. Huang, X.-Y. Kuang, and Y. Huang, 2013: Decadal changes of Meiyu rainfall around 1991 and its relationship with two types of ENSO. J. Geophys. Res. Atmos., 118, 9766–9777, https://doi.org/10.1002/jgrd.50779.
Zhu, L., G. Huang, G. Fan, X. Qu, Z. Wang, and W. Hua, 2019: Elevation-dependent sensible heat flux trend over the Tibetan Plateau and its possible causes. Climate Dyn., 52, 3997–4009, https://doi.org/10.1007/s00382-018-4360-7.
Zhu, Z., T. Li, and J. He, 2014: Out-of-phase relationship between boreal spring and summer decadal rainfall changes in southern China. J. Climate, 27, 1083–1099, https://doi.org/10.1175/JCLI-D-13-00180.1.