1. Introduction
The East Asian winter monsoon (EAWM) is one of the most active climate systems in the boreal winter. It has great impacts on the weather and climate of China and surrounding regions (Lau and Li 1984; Ding 1994; He et al. 2007). Based on the observational studies, the EAWM exhibits significant interannual variabilities that are usually accompanied by anomalous cold surges (Gu et al. 2008; Yang et al. 2018), snowfall (Sun et al. 2010; Wang and He 2013), rainfall (Jia and Ge 2017), and haze pollution (Li et al. 2016; Wang and Chen 2016) over East Asia. On this time scale, the EAWM is strongly affected by El Niño–Southern Oscillation (ENSO) (Zhang et al. 1996; Wang et al. 2000; Zhou et al. 2007; Wang and He 2012). Usually, the EAWM is weaker (stronger) when El Niño (La Niña) occurs (Wang et al. 2000). In addition, the North Atlantic Oscillation (Wu and Huang 1999), the Arctic Oscillation (AO) (Gong et al. 2001; Wu and Wang 2002; He et al. 2017; Li et al. 2018), and autumn Arctic sea ice (Liu et al. 2012; C. Sun et al. 2016; Wang and Liu 2016) can also influence the EAWM through both thermodynamic and dynamic processes.
In recent decades, research on the interdecadal changes of the EAWM has raised more concerns. They are the background conditions of the interannual variations and modulate long-term changes in the East Asian climate (Wang and Fan 2013). Some studies indicate that the EAWM experienced significant weakening in the mid-1980s (Xu et al. 1999; Wang and He 2012; Miao et al. 2018b). Usually, reasons for the interdecadal changes of the East Asian climate are very complex (Ding et al. 2014). Influences from external forcings on the East Asian climate (Wang et al. 2013; Miao et al. 2016; Jiang et al. 2017; Zhou 2017; Chen and Dong 2019; Miao et al. 2018a; Tian et al. 2018) and even on the large-scale Pacific climate (Wang et al. 2012; Xiao et al. 2017; Miao et al. 2018c) are significant. A recent modeling study suggested that both greenhouse gases and natural external forcings played important roles in weakening the EAWM and related changes over the North Pacific in the mid-1980s (Miao et al. 2018b). In this process, influences from the internal variabilities cannot be neglected. For example, Ding et al. (2014) noted that the Pacific decadal oscillation (PDO) and the EAWM are significantly negatively correlated on the interdecadal time scale in the observations. In addition, some previous studies indicated that on multidecadal time scales, a weaker EAWM is linked to a warm phase of the Atlantic multidecadal oscillation (AMO), and vice versa (Li and Bates 2007; Wang et al. 2009). However, in the observations, the time of phase transient for the PDO and AMO is not consistent with the interdecadal weakening of the EAWM around the mid-1980s. Therefore, how these internal modes affect decadal variations of the EAWM needs further investigation.
In fact, it is difficult to use only an observation-based analysis to investigate the internal variabilities of the climate system at the interdecadal and longer time scales due to the limitation of observational data length and the influence of mixed external forcing signals. Thus, coupled climate models become a useful tool to explore such processes and their associated mechanisms (Lei et al. 2014). As the external forcings remain constant, the internal dynamics of the coupled atmosphere–ocean–land–cryosphere system induce internal variability on various time scales (Hasselmann 1976). Thus, many previous studies compared unperturbed control simulations with observations to detect internal variabilities and externally forced climate changes. Stouffer et al. (2000) analyzed the internal variability of surface air temperature (SAT) in three 1000-yr coupled model integrations and suggested that the observed warming of ~0.5°C over the twentieth century is not due to internally generated variability. Lei et al. (2014) explored the interdecadal natural variability of summer rainfall over China using multicentury simulations of the HadCM3 model and concluded that both anthropogenic and natural factors could play a role in recent changes in Chinese summer rainfall. Furthermore, Polvani and Smith (2013) suggested that the recent observed positive trend in Antarctic sea ice lies well within the natural variability shown in CMIP5 coupled models. Control simulations can also be used to investigate the mechanisms of observed phenomena (e.g., Kug et al. 2006; Cui et al. 2013; Yan et al. 2018). For example, a recent study analyzed the dynamic processes accounting for the observed linkage between the Atlantic meridional overturning circulation and Atlantic multidecadal variability using control simulations from the CMIP3 and CMIP5 (Yan et al. 2018). Focusing on the EAWM, can the internal variability of the climate system cause the interdecadal variation of EAWM as observed around the mid-1980s? And what is the reason for the process? These issues remain unresolved.
To address the abovementioned issues, in this study, we investigate the interdecadal variations of the EAWM using multiple long-term preindustrial simulations within the CMIP5. In section 2, we describe the data and methods used here. In section 3, we evaluate the multimodels’ performance in simulating the EAWM. In section 4, we investigate the interdecadal variations of the EAWM in the preindustrial simulations and the associated physical mechanisms. The conclusions and discussion are given in section 5.
2. Data and method
Model results, including historical and 500-yr preindustrial control (piControl) simulations from 28 coupled models within the CMIP5 (Taylor et al. 2012), are used in this study. The details for these models are illustrated in Table 1. The historical simulations are integrated with both anthropogenic and natural external forcings. For the piControl simulations, however, the models impose nonevolving preindustrial conditions. The last 500 years in each piControl simulation are analyzed in this study.
Descriptions of the CMIP5 models used in this study.
The monthly sea level pressure (SLP) from the Hadley Centre’s HadSLP2 dataset (Allan and Ansell 2006) and the geopotential height and the wind fields from the NCEP–NCAR reanalysis data (Kalnay et al. 1996) are used here to evaluate the models’ performances in simulating the East Asian winter (December–February) climate and to analyze interdecadal climate changes. To obtain the EAWM index for a long period, the Twentieth Century Reanalysis, version 2 (20CR V2), data covering the period of 1871–2012 are adopted (Compo et al. 2011). In addition, the AO index calculated using these data was downloaded from the website (https://www.esrl.noaa.gov/psd/data/20thC_Rean/timeseries/monthly/AO/) to examine the AO–EAWM relationship in the observation. Furthermore, the PDO index is obtained from the website (http://research.jisao.washington.edu/pdo/PDO.latest.txt) to illustrate the observed PDO–EAWM relationship.
Additionally, some numerical experiments are performed to investigate influences from the PDO-like SST anomalies on the East Asian winter climate. The model used is the Community Atmosphere Model, version 4 (CAM4), which is the atmospheric component of the Community Earth System Model, version 1.0.5 (CESM1.0.5) (Gent et al. 2011). The horizontal resolution is 1.9° × 2.5°, and the vertical direction has 26 layers.
3. Evaluation of the CMIP5 models and CAM4
First, we evaluate the CMIP5 models’ performances in simulating winter climatology over East Asia. Focusing on variations of the EAWM, especially on the circulation, we mainly assess the simulated SLP, Z500, and U300 over the East Asian region (20°–50°N, 100°–145°E) in the CMIP5 historical simulations. The Taylor diagram is used to show the quantitative evaluation (Fig. 1), which includes the spatial correlation coefficient (SCC) and the ratio of standard deviation (RSD) between the models and the observations. In the following sections, we analyze the model results, which exhibit better performances in simulating the East Asia winter climate based on the Taylor diagram. That is, the SCC must be at least 0.9, and the RSD must be between 0.75 and 1.25, which are employed in previous studies (e.g., Jiang and Tian 2013).
Taylor diagram of climatological winter SLP (hPa), Z500 (m), and U300 (m s−1) during 1967–2004 over East Asia (20°–50°N, 100°–145°E).
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
For the SLP, the SCC range is from 0.74 (MIROC-ESM) to 0.98 (ACCESS1.0; HadGESM2-ES). We exclude six models with an SCC less than 0.9. They are IPSL-CM5A-LR, FGOALS-g2, MIROC-ESM, MIROC4h, CCSM4, and CESM1(BGC). Most RSD values are between 0.75 and 1.25, except for GISS-E2-R (0.65) and IPSL-CM5A-LR (1.28). For Z500, the SCCs are larger than 0.99, and the RSDs are near 1.0 in most models. For U300, the SCCs are larger than 0.95 in all models. However, the RSDs are larger than 1.25 in CSIRO Mk3.6.0, IPSL-CM5A-LR and MRI-CGCM3, and these three models are thus excluded. Overall, nine models [CSIRO Mk3.6.0, IPSL-CM5A-LR, FGOALS-g2, MIROC-ESM, MIROC4h, MRI-CGCM3, GISS-E2-R, CCSM4, and CESM1(BGC)] are excluded, and the remaining 19 models, which show better performances in simulating the climatology of the EAWM, are used to investigate the internal variabilities of the EAWM on the interdecadal time scale.
In addition, some previous studies have indicated that the PDO patterns (e.g., Yim et al. 2015; Wang and Miao, 2018), the AO pattern (e.g., Zuo et al. 2013), and the AMO pattern (e.g., Han et al. 2016) in most CMIP5 coupled models are in good agreement with the observation. According to detailed information from these previous evaluations, the models analyzed in the present study have reliable simulations of these modes. Thus, we will not repeat the evaluation of these internal modes in this study.
For the CAM4, We carried out a “F_2000” simulation forced by the observed climatological (1986–2005) SST and sea ice concentration (SIC). The results show that CAM4 can simulate well the present East Asian winter climate (figure not shown). It is a good starting point to address the following issues.
4. Internal variability of the EAWM
a. The observed and simulated interdecadal weakening of the EAWM
Figure 2a illustrates the observed EAWM index and its variation on an interdecadal time scale. The larger (smaller) values of this index denote a stronger (weaker) observed EAWM. As shown in the 9-yr running mean of the EAWM index, the EAWM undergoes significant interdecadal weakening during 1967–2006. The EAWM is relatively strong during 1967–86 but becomes weaker during 1987–2006. The t value of the moving t test reaches a minimum in 1986 (Fig. 2b), which suggests that there is an abrupt change in the EAWM during the mid-1980s. We thus choose two 20-yr subperiods (i.e., 1967–86 and 1987–2006) to investigate this observed interdecadal change in the EAWM. As shown in Fig. 3a, negative SLP anomalies can be found over the northern Asian continent, whereas positive SLP anomalies are evident over southern China, indicating that the Siberian high is weakened during the latter period. Correspondingly, significant southwesterly wind anomalies are observed over northern East Asia in the lower troposphere (Fig. 3b). It is not conducive to the southward invasion of the cold air from high latitudes to East Asia. In addition, during the latter period, the EAT becomes shallower, and the EAJS shifts poleward (Figs. 3c,d). Therefore, the results indicate that the EAWM underwent a significant interdecadal weakening after the mid-1980s. As noted in a previous study, external forcings are likely to play an important role in this weakening (Miao et al. 2018b). However, how the internal variabilities contribute to interdecadal changes of the EAWM remains unclear. In the following section, this observed interdecadal weakening of the EAWM is used to make the model–data comparison.
(a) The EAWM index (gray bars) and its decadal component (9-yr running mean; line) for the period of 1948–2017 in the observation. (b) The moving t test values (gray bars) with two subseries (each duration is 20 years), and the 0.05 significance level (dashed lines). The maximum value corresponds to two subseries from 1967 to 2006 [black bar in (a)].
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
The differences in winter (a) SLP (hPa), (b) 850-hPa wind fields (UV850) (m s−1), (c) Z500 (m), and (d) U300 (m s−1) between two subperiods of 1987–2006 and 1967–86. They are calculated based on (a) the HadSLP2 and (b)–(d) the NCEP–NCAR datasets, respectively. Areas exceeding the 95% confidence level are denoted with dots in (a), (c), and (d) or shaded gray in (b) based on the Student’s t test.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
In fact, similar interdecadal weakening of the EAWM can be found in the CMIP5 piControl simulations. Based on the detection by using the t-test method, there are in total 53 significant interdecadal weakening cases of EAWM in the selected 19 piControl simulations (Table 2). As shown in Fig. 4a, the 9-yr running mean of the EAWM indices all show interdecadal weakening of the EAWM, especially their MCE. We further examine the large-scale common characteristics of these interdecadal weakening of the EAWM in the different coupled models.
Number of EAWM weakening and strengthening cases selected in each model. The gaps between two cases are at least 40 years.
The EAWM indices (9-yr running mean; gray) of the selected (a) weakening and (b) strengthening cases of the EAWM in the 19 CMIP5 models and their ensemble mean (black). The red line in (a) stands for the observed EAWM index for the period of 1967–2006.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
Figure 5a shows the MCE of the differences in SLP between two subperiods in all weakening cases. Negative SLP anomalies can be found over the Arctic and Siberian regions, which show a positive AO pattern. These negative SLP anomalies are surrounded by the positive ones. Two positive SLP anomaly centers are located at the North Pacific and the Mediterranean, respectively. Particularly over the North Pacific, the positive SLP anomalies are much stronger, suggesting a positive NPO and a weakened Aleutian low. The multicase consistency indicates that more than 90% of the weakening cases show positive AO and NPO patterns (Fig. 5b). Correspondingly, anticyclonic circulation anomalies are evident over the North Pacific, and southerly wind anomalies can be found over the Asian continent in the lower troposphere (Fig. 6a). The multicase consistency for wind anomalies is larger than 90% over East Asia and the North Pacific (Figs. 6b,c). These results suggest that the common features of the interdecadal weakening cases of the EAWM are a positive AO/NPO and an associated anomalously anticyclonic circulation. In the middle troposphere, more than 90% of the weakening cases show positive Z500 anomalies over the midlatitude Asian–Pacific region (Figs. 7a,b). In contrast, negative Z500 anomalies can be found over high-latitude regions; the multicase consistency is also larger than 90% in those regions. This meridional dipole pattern indicates that the EAT is weakened during the latter period. In the upper troposphere, negative (positive) U300 anomalies are evident south (north) of 40°N over the Asian–Pacific regions (Fig. 8a). The multicase consistency is larger than 90% over these regions (Fig. 8b), suggesting that the EAJS shifts poleward in most EAWM weakening cases. Overall, all the members of the EAWM system, including the Siberian high, the low-level circulation, the EAT, and the EAJS, match well the anomalous AO/NPO patterns and show an interdecadal weakening of the EAWM in these cases.
(a) The simulated MCE of the differences in winter SLP (hPa) between the two subperiods in the interdecadal weakening cases of the EAWM and (b) their multicase consistency (%). (c),(d) As in (a) and (b), but for the interdecadal strengthening cases of the EAWM.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
(a) The simulated MCE of the differences in winter UV850 (m s−1) between the two subperiods in the interdecadal weakening cases of the EAWM, and (b),(c) their multicase consistency (%). (d)–(f) As in (a)–(c), but for the interdecadal strengthening cases of the EAWM.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
(a) The simulated MCE of the differences in winter Z500 (m) between the two subperiods in the interdecadal weakening cases of the EAWM and (b) their multicase consistency (%). (c),(d) As in (a) and (b), but for the interdecadal strengthening cases of the EAWM.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
(a) The simulated MCE of the differences in winter U300 (m s−1) between the two subperiods in the interdecadal weakening cases of the EAWM and (b) their multicase consistency (%). (c),(d) As in (a) and (b), but for the interdecadal strengthening cases of the EAWM.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
b. Simulated interdecadal strengthening of the EAWM
In addition, we also find 49 cases of significant interdecadal strengthening of the EAWM in the CMIP5 simulations (Table 2). In these cases, the 9-yr running means of the EAWM indices increase during the latter subperiod (Fig. 4b). As indicated by anomalous spatial patterns, the SLPs increase over the mid- and high-latitude Asian regions (Fig. 5c), and the Z500s decrease over northeastern China and southern Japan (Fig. 7c), implying a strengthened Siberian high and a deepened East Asian trough. Northeasterly wind anomalies are evident over East China (Fig. 6d), and the EAJS shifts equatorward (Fig. 8c). In these cases, the EAWM experiences significant interdecadal strengthening.
At the same time, positive SLP anomalies are evident over the Arctic regions, which are accompanied by negative SLP anomalies over the North Pacific and North Europe (Fig. 5c). It suggests a negative AO phase and a negative NPO phase in most interdecadal strengthening cases of the EAWM (Fig. 5d). As a result, cyclonic circulation anomalies are seen over the North Pacific (Fig. 6d). The multicase consistencies for the wind fields are larger than 90% over the Asian–Pacific regions (Figs. 6e,f). Overall, significant interdecadal strengthening cases of the EAWM can also be found in the CMIP5 piControl simulations. The EAWM members and the associated large-scale circulation show changes that are opposite to the weakening cases.
c. Mechanisms for the EAWM interdecadal variations
Based on observational studies, the decadal variations of the East Asian winter climate are closely related to the natural variabilities of the ocean (e.g., Li et al. 2009; J. Sun et al. 2016), and phase shifting of the PDO is likely one of the most important reasons for that (e.g., Zhu and Yang 2003; Ding et al. 2014). Indeed, in most cases of interdecadal weakening of the EAWM, we can find similar, anomalously negative PDO-like SST patterns in the later subperiods relative to the former ones. Positive SST anomalies are evident over the western North Pacific (Fig. 9a). The maximum positive anomalies are observed over the Kuroshio–Oyashio Extension region. Additionally, these positive SST anomalies are surrounded by negative ones in the North Pacific. Particularly in the high-latitude North Pacific, the negative SST anomalies are relatively stronger. It suggests that the PDO enters its negative phase from the former subperiod to the latter one. Changes in the MCE of the PDO indices further confirm this phase shifting (Fig. 10a). In the interdecadal strengthening cases of the EAWM, the changes in the PDO are reversed (Figs. 9c and 10b). Thus, the weakening (strengthening) of the EAWM generally correlates with the negative (positive) phase of the PDO. However, the PDO phase transition occurs later than the interdecadal variation of the EAWM, implying that the PDO may not be the cause of the interdecadal variations of the EAWM. We calculate the correlation coefficients of the low-pass-filtered PDO and EAWM indices in each 500-yr piControl simulation (Fig. 11a). The results show that the relationships between the PDO and the EAWM are close in most models (14 of 19 models). The correlation coefficients are positive in 18 models and significant (p < 0.05) in 11 models. Nevertheless, changes in the EAWM lead the PDO in these long-term simulations (Fig. 12). In most models, the simulated EAWM indices lead the PDO indices by 2–4 years, suggesting that the PDO cannot be the reason for the interdecadal weakening (strengthening) of the EAWM in the piControl simulations.
(a) The simulated MCE of the differences in winter SST (°C) between the two subperiods in the interdecadal weakening cases of the EAWM and (b) their multicase consistency (%). (c),(d) As in (a) and (b), but for the interdecadal strengthening cases of the EAWM.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
The winter (a) PDO, (c) AO, and (e) NPO indices (9-yr running mean; gray) of the EAWM weakening cases in the models and their MCE (black). The blue lines are the MCE of the EAWM index in the EAWM weakening cases in the models that presented in Fig. 4. (b),(d),(f) As in (a), (c), and (e), but for the EAWM strengthening cases.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
The correlation coefficient between the winter values of (a) PDO, (b) AO, and (c) NPO and the EAWM index for the 500-yr piControl simulations (the 9-yr low-pass filter is employed on the all indices). The symbols indicate the significance level.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
The cross-correlation coefficient between the winter PDO and EAWM index for the 500-yr piControl simulations (9-yr low-pass filtered). The positive value of x axis means that PDO leads EAWM index, and vice versa. The dashed lines indicate the 0.05 significance level.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
We carry out some additional CAM4 experiments to examine the influence of the PDO on the EAWM. The simulation is set up using the “F_1850” configuration. The atmospheric composition constant is set in the year 1850. A total of five simulations are carried out (Table S1 in the online supplemental material). All the simulations are integrated for 40 years, and the averages of winters in the last 30 years are analyzed. The first simulation is a control simulation, which is forced by the CAM4’s preindustrial climatological SST and sea ice boundary conditions. This control simulation is herein abbreviated as EXP1. In the second simulation (EXP2), an anomalously negative PDO-like SST pattern is added to the climatological SST boundary condition (shown in Fig. 9a). Specifically, only the SST anomalies where the multicase consistency is larger than 70% are added into the North Pacific region (20°–60°N, 110°E–100°W). Similarly, an anomalously positive PDO-like SST pattern (Fig. 9c) is added into the climatological SST in the third simulation (EXP3). The differences between EXP2 (EXP3) and EXP1 reflect the impacts of the PDO-like SST anomalies in the piControl simulations on the climate. Furthermore, we perform sensitivity simulations named EXP4 and EXP5. They are similar to EXP2 and EXP3, but the added SST anomalies are doubled in these two simulations. The CAM4 simulations show that the PDO-like SST pattern could not cause significant changes in the EAWM in the AGCM experiments, even when the SST anomalies are doubled (s. S1 and S2). These results confirm that the PDO is unlikely to be the reason for the interdecadal variations of the EAWM, as suggested by the CMIP5 coupled model results.
In addition to the natural variabilities in the ocean, the interdecadal variations of the EAWM are also related to major modes of natural variability in atmospheric circulation (summarized in Ding et al. 2014). As mentioned above, the SLP anomalies in the EAWM weakening (strengthening) cases resemble the positive (negative) phase of the AO/NPO, respectively (Fig. 5). In each coupled model, we calculate the winter AO and NPO indices in the 500-yr piControl simulations and the correlation coefficients between the 9-yr low-pass-filtered AO/NPO and EAWM indices (Figs. 11b,c). We find that the correlation coefficients between AO and EAWM indices are significantly negative in all the models at the 0.05 significance level. For the NPO and EAWM indices, the correlation coefficients are all negative and significant in the most models. The weakening (strengthening) of the EAWM is thus highly correlated with the positive (negative) phase of the AO/NPO. We further examine the AO and NPO indices in the models’ EAWM weakening (strengthening) cases (Figs. 10c–f) and find that the AO and NPO also experience an interdecadal phase shift. The phase shifting of the MCEs of the AO and NPO indices is consistent with those of the EAWM, and the interdecadal changes of the MCEs of the AO indices are larger. This means that the AO could play an important role in causing interdecadal variations of the EAWM. Both the NPO phase shifting and PDO-like SST anomalies are possibly caused by the interdecadal variations of the AO, which needs further investigation.
We further show the cross-correlation coefficients between the low-pass-filtered AO and EAWM indices for the 500-yr simulation in each model. As shown in Fig. 13, the correlation coefficients are minimized and become significant at the zero point in all the simulations, which implies that the AO is in phase with the EAWM at a decadal time scale. To confirm this inferred linkage, we also examine the ensemble mean of the regression map of the EAWM-related variables on the AO index (9-yr low-pass filtered) in each model’s 500-yr piControl simulation. As shown in Fig. 14a, during the positive phase of the AO, negative SLP anomalies can be seen over the Arctic and Siberian regions. In contrast, positive SLP anomalies are evident over the North Pacific and the North Atlantic. Both the Siberian high and the Aleutian low are thus weakened. Correspondingly, anomalous anticyclonic circulation can be found over the North Pacific, and southerly winds are obvious over the northern Asian continent and eastern China (Fig. 14b). In the middle troposphere, positive Z500 anomalies over the midlatitude Asian–Pacific region reveal that the EAT is weakened (Fig. 14c). In the upper troposphere, the EAJS is weakened and shifted poleward (Fig. 14d). The patterns for the anomalous atmospheric circulations resemble those in the EAWM weakening cases shown in Figs. 5–8, and vice versa (not shown). Therefore, the phasing shifting of the AO should be the main reason for the interdecadal variations of the EAWM in the piControl simulations.
The cross-correlation coefficient between the winter AO and EAWM index for the 500-yr piControl simulations (9-yr low-pass filtered). The positive value of x axis means that AO leads the EAWM index, and vice versa. The dashed lines indicate the 0.05 significance level.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
The ensemble mean of the regression map of the winter (a) SLP (hPa), (b) UV850 (m s−1), (c) Z500 (m), and (d) U300 (m s−1) on the AO index in each model’s 500-yr piControl simulation. The atmospheric variables and AO index are both 9-yr low-pass filtered.
Citation: Journal of Climate 33, 2; 10.1175/JCLI-D-19-0148.1
5. Discussion and conclusions
In this study, we investigate the internal variability of the EAWM on the interdecadal time scale using the piControl simulations output from CMIP5 coupled models. Nineteen simulations from the selected coupled models are analyzed here. A total of 53 cases of significant interdecadal weakening of the EAWM can be found in these 500-yr piControl simulations. We further illustrate the ensemble differences in the EAWM-related atmospheric circulations between two subperiods of these cases to find the commonalities among them. The result indicates that when the EAWM weakens in the piControl simulations, positive AO and NPO patterns can be found in more than 90% of cases. Correspondingly, anticyclonic circulation anomalies are obvious over the North Pacific, and southerly wind anomalies can be found over the Asian continent in the lower troposphere. In addition, the EAT in the midtroposphere is weakened, and the EAJS in the upper troposphere is weakened and shifted poleward. These simulated changes resemble the observed interdecadal weakening of EAWM around the mid-1980s. In instances of interdecadal strengthening of the EAWM, the anomalous patterns are opposite. We further investigate the mechanisms for the interdecadal variations of the EAWM in the piControl simulations. The result suggests that the interdecadal phase shifting of the AO is likely the main reason for the EAWM interdecadal variations in the piControl simulations.
As shown in Table 2, there are few cases of significant interdecadal weakening and strengthening of the EAWM in each 500-yr simulation. The average number of weakening and strengthening cases is approximately 2.8 (500 yr)−1 and 2.6 (500 yr)−1, respectively. This means that the probability of interdecadal variations of the EAWM is very low in an unforced climate system.
Though the internal variability of the climate system can cause interdecadal weakening of the EAWM, the associated large-scale circulation anomalies exhibit some differences from those observed in the mid-1980s mainly over the North Pacific. It implies that their root causes could be different. As noted by our previous study (Miao et al. 2018b), the multimodel ensemble of the CMIP5 historical simulations, which excludes internal variability of the climate system, can reproduce well the weakening shift of the EAWM in the observation. It means that the external forcing agencies are likely to play a dominate role in causing the observed interdecadal weakening of EAWM in the mid-1980s. Compared with the piControl simulations, in addition, the simulated weakening in the CMIP5 historical simulations indeed resembles the observed one better. The differences between the piControl simulations and observation/historical simulations also confirm the important influences from external forcings on the interdecadal weakening of the EAWM around the mid-1980s. In fact, it is difficult to quantify the relative contributions of the external forcing and the internal variability to the observed interdecadal weakening of the EAWM due to nonlinearity of the climate system. This present result suggests that the internal variability of the climate system may also contribute to the observed interdecadal weakening of EAWM in the mid-1980s. If it is true, the role of the internal variability could be secondary based on the CMIP5 multimodel results.
Our study emphasizes the AO’s role in causing the interdecadal variations of the EAWM in the no transient external agents forced climate system. The correlation coefficients between the 9-yr low-pass-filtered AO and EAWM indices are significantly negative in all the models at the 0.05 significance level (Fig. 11b). In the observation, during the period of 1951–2012, the correlation coefficient between the EAWM and the AO indices (9-yr running mean) is −0.84 (p = 0.08; Fig. S3), similar to the relationship between them in the piControl simulations. And the observed AO phase shifting could contribute to the interdecadal weakening of the EAWM around the mid-1980s (He and Wang 2012). However, during a longer period (1872–2012), the correlation coefficient between them is −0.43 (nonsignificant). It means that the AO–EAWM relationship is unstable in the observation. We further examined the relationship between the AO and EAWM in the CMIP5 historical simulations. More than half the models cannot reproduce the significant correlation between the AO and EAWM (Fig. S4), which is very different from the piControl simulations. The external forcings, which can regulate both the AO (e.g., Shindell et al. 2001; Wang et al. 2012) and the EAWM (e.g., Miao et al. 2016, 2018a), could be one of the reasons for the observed unstable AO–EAWM relationships. Some further studies are needed on this point.
In the observation, the correlation coefficient between EAWM and PDO is only 0.32 (nonsignificant) for the period of 1901–2012 (Fig. S5). In addition, the PDO phase shifting occurred around the late 1970s, which is earlier by 5–6 years relative to the interdecadal weakening of the EAWM around the mid-1980s. This also confirms our previous conclusion that the PDO phase shifting could not be the reason for the interdecadal variations of the EAWM.
As noted in the introduction, ENSO has significant impacts on the EAWM at the interannual time scale. The positive (negative) phase of ENSO favors a weaker (stronger) EAWM. On the interdecadal time scale, however, the influence from ENSO is not significant. As shown in Fig. 9, there are no significant SST anomalies over the tropical eastern Pacific. This suggests that ENSO has no influence on the interdecadal variations of the EAWM in the piControl simulations.
Acknowledgments
We thank two anonymous reviewers and the editor for their valuable comments and suggestions, which helped to improve the quality of this paper significantly. This research was supported by the National Natural Science Foundation of China (41421004, 41575086, and 41661144005), the CAS–PKU Joint Research Program, and the Scientific Research Foundation of Joint Laboratory of Climate and Environment Change from Chengdu University of Information Technology (JLCEC201802). We thank the IPCC for providing the CMIP5 datasets (http://www.ipcc-data.org/sim/gcm_monthly/AR5/Reference-Archive.html).
REFERENCES
Allan, R., and T. Ansell, 2006: A new globally complete monthly historical gridded mean sea level pressure dataset (HadSLP2): 1850–2004. J. Climate, 19, 5816–5842, https://doi.org/10.1175/JCLI3937.1.
Bretherton, C. S., M. Widmann, V. P. Dymnikov, J. M. Wallace, and I. Bladé, 1999: The effective number of spatial degrees of freedom of a time-varying field. J. Climate, 12, 1990–2009, https://doi.org/10.1175/1520-0442(1999)012<1990:TENOSD>2.0.CO;2.
Chen, W., and B. Dong, 2019: Anthropogenic impacts on recent decadal change in temperature extremes over China: Relative roles of greenhouse gases and anthropogenic aerosols. Climate Dyn., 52, 3643–3660, https://doi.org/10.1007/S00382-018-4342-9.
Compo, G. P., and Coauthors, 2011: The Twentieth Century Reanalysis Project. Quart. J. Roy. Meteor. Soc., 137, 1–28, https://doi.org/10.1002/qj.776.
Cui, X., Y. Gao, D. Gong, D. Guo, and F. Tore, 2013: Teleconnection between winter Arctic Oscillation and Southeast Asian summer monsoon in the pre-industry simulation of a coupled climate model. Atmos. Ocean. Sci. Lett., 6, 349–354, https://doi.org/10.3878/j.issn.1674-2834.12.0114.
Ding, Y., 1994: Monsoons over China. Kluwer Academic Publishers, 419 pp.
Ding, Y., and Coauthors, 2014: Interdecadal variability of the East Asian winter monsoon and its possible links to global climate change. J. Meteor. Res., 28, 693–713, https://doi.org/10.1007/s13351-014-4046-y.
Gent, P. R., and Coauthors, 2011: The Community Climate System Model version 4. J. Climate, 24, 4973–4991, http0073://doi.org/10.1175/2011JCLI4083.1.
Gong, D. Y., S. W. Wang, and J. H. Zhu, 2001: East Asian winter monsoon and Arctic Oscillation. Geophys. Res. Lett., 28, 2073–2076, https://doi.org/10.1029/2000GL012311.
Gu, L., K. Wei, and R. H. Huang, 2008: Severe disaster of blizzard, freezing rain and low temperature in January 2008 in China and its association with the anomalies of East Asian monsoon system (in Chinese). Climatic Environ. Res., 13, 405–418.
Guo, D., and Z. B. Sun, 2004: Relationship of winter North Pacific oscillation anomalies with the East Asian winter monsoon and the weather and climate in China (in Chinese). J. Nanjing Inst. Meteor., 27, 461–470.
Han, Z., F. Luo, S. Li, Y. Gao, T. Furevik, and L. Svendsen, 2016: Simulation by CMIP5 models of the Atlantic multidecadal oscillation and its climate impacts. Adv. Atmos. Sci., 33, 1329–1342, https://doi.org/10.1007/s00376-016-5270-4.
Hasselmann, K., 1976: Stochastic climate models. Part I. Theory. Tellus, 28, 473–485, https://doi.org/10.3402/tellusa.v28i6.11316.
He, J. H., J. H. Ju, Z. P. Wen, J. M. Lu, and Q. H. Jin, 2007: A review of recent advances in research on Asian monsoon in China. Adv. Atmos. Sci., 24, 972–992, https://doi.org/10.1007/s00376-007-0972-2.
He, S., and H. Wang, 2012: An integrated East Asian winter monsoon index and its interannual variability (in Chinese). Chin. J. Atmos. Sci., 36, 523–538.
He, S., Y. Gao, F. Li, H. Wang, and Y. He, 2017: Impact of Arctic Oscillation on the East Asian climate: A review. Earth Sci. Rev., 164, 48–62, https://doi.org/10.1016/j.earscirev.2016.10.014.
Jia, X., and J. Ge, 2017: Interdecadal changes in the relationship between ENSO, EAWM, and the wintertime precipitation over China at the end of the twentieth century. J. Climate, 30, 1923–1937, https://doi.org/10.1175/JCLI-D-16-0422.1.
Jiang, D. B., and Z. P. Tian, 2013: East Asian monsoon change for the 21st century: Results of CMIP3 and CMIP5 models. Chin. Sci. Bull., 58, 1427–1435, https://doi.org/10.1007/s11434-012-5533-0.
Jiang, Y., and Coauthors, 2017: Anthropogenic aerosol effects on East Asian winter monsoon: The role of black carbon-induced Tibetan Plateau warming. J. Geophys. Res. Atmos., 122, 5883–5902, https://doi.org/10.1002/2016JD026237.
Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437–471, https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.
Kug, J.-S., T. Li, S.-I. An, I.-S. Kang, J.-J. Luo, S. Masson, and T. Yamagata, 2006: Role of the ENSO–Indian Ocean coupling on ENSO variability in a coupled GCM. Geophys. Res. Lett., 33, L09710, https://doi.org/10.1029/2005GL024916.
Lau, K.-M., and M.-T. Li, 1984: The monsoon of East Asia and its global associations—A survey. Bull. Amer. Meteor. Soc., 65, 114–125, https://doi.org/10.1175/1520-0477(1984)065<0114:TMOEAA>2.0.CO;2.
Lei, Y., B. Hoskins, and J. Slingo, 2014: Natural variability of summer rainfall over China in HadCM3. Climate Dyn., 42, 417–432, https://doi.org/10.1007/s00382-013-1726-8.
Li, Q., R. Zhang, and Y. Wang, 2016: Interannual variation of the wintertime fog–haze days across central and eastern China and its relation with East Asian winter monsoon. Int. J. Climatol., 36, 346–354, https://doi.org/10.1002/joc.4350.
Li, S., and G. T. Bates, 2007: Influence of the Atlantic Multidecadal Oscillation on the winter climate of East China. Adv. Atmos. Sci., 24, 126–135, https://doi.org/10.1007/s00376-007-0126-6.
Li, S., Y. Wang, and Y. Gao, 2009: A review of the researches on the Atlantic multidecadal oscillation (AMO) and its climate influence (in Chinese). Trans. Atmos. Sci., 32, 458–465.
Li, S., S.-P. He, F. Li, and H.-J. Wang, 2018: Simulated and projected relationship between the East Asian winter monsoon and winter Arctic Oscillation in CMIP5 models. Atmos. Ocean. Sci. Lett., 11, 417–424, https://doi.org/10.1080/16742834.2018.1512356.
Liu, J., J. A. Curry, H. Wang, M. Song, and R. M. Horton, 2012: Impact of declining Arctic sea ice on winter snowfall. Proc. Natl. Acad. Sci. United States, 109, 4074–4079, https://doi.org/10.1073/pnas.1114910109.
Mantua, N. J., S. R. Hare, Y. Zhang, J. M. Wallace, and R. C. Francis, 1997: A Pacific interdecadal climate oscillation with impacts on salmon production. Bull. Amer. Meteor. Soc., 78, 1069–1080, https://doi.org/10.1175/1520-0477(1997)078<1069:APICOW>2.0.CO;2.
Miao, J. P., T. Wang, Y. L. Zhu, J. Z. Min, H. J. Wang, and D. Guo, 2016: Response of the East Asian winter monsoon to strong tropical volcanic eruptions. J. Climate, 29, 5041–5057, https://doi.org/10.1175/JCLI-D-15-0600.1.
Miao, J. P., T. Wang, H. J. Wang, and Y. Q. Gao, 2018a: Influence of low-frequency solar forcing on the East Asian winter monsoon based on HadCM3 and observations. Adv. Atmos. Sci., 35, 1205–1215, https://doi.org/10.1007/s00376-018-7229-0.
Miao, J. P., T. Wang, H. J. Wang, Y. L. Zhu, and J. Q. Sun, 2018b: Interdecadal weakening of the East Asian winter monsoon in the mid-1980s: The roles of external forcings. J. Climate, 31, 8985–9000, https://doi.org/10.1175/JCLI-D-17-0868.1.
Miao, J. P., T. Wang, H. J. Wang, and J. Q. Sun, 2018c: Interannual weakening of the tropical Pacific Walker circulation due to strong tropical volcanism. Adv. Atmos. Sci., 35, 645–658, https://doi.org/10.1007/s00376-017-7134-y.
Polvani, L. M., and K. L. Smith, 2013: Can natural variability explain observed Antarctic sea ice trends? New modeling evidence from CMIP5. Geophys. Res. Lett., 40, 3195–3199, https://doi.org/10.1002/grl.50578.
Shindell, D. T., G. A. Schmidt, M. E. Mann, D. Rind, and A. Waple, 2001: Solar forcing of regional climate change during the Maunder Minimum. Science, 294, 2149–2152, https://doi.org/10.1126/science.1064363.
Stouffer, R. J., G. Hegerl, and S. Tett, 2000: A comparison of surface air temperature variability in three 1000-yr coupled ocean–atmosphere model integrations. J. Climate, 13, 513–537, https://doi.org/10.1175/1520-0442(2000)013<0513:ACOSAT>2.0.CO;2.
Sun, C., S. Yang, W. Li, R. Zhang, and R. Wu, 2016: Interannual variations of the dominant modes of East Asian winter monsoon and possible links to Arctic sea ice. Climate Dyn., 47, 481–496, https://doi.org/10.1007/s00382-015-2851-3.
Sun, J., H. Wang, W. Yuan, and H. Chen, 2010: Spatial-temporal features of intense snowfall events in China and their possible change. J. Geophys. Res., 115, D16110, https://doi.org/10.1029/2009JD013541.
Sun, J., S. Wu, and J. Ao, 2016: Role of the North Pacific sea surface temperature in the East Asian winter monsoon decadal variability. Climate Dyn., 46, 3793–3805, https://doi.org/10.1007/s00382-015-2805-9.
Taylor, K. E., R. J. Stouffer, and G. A. Meehl, 2012: An overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc., 93, 485–498, https://doi.org/10.1175/BAMS-D-11-00094.1.
Thompson, D. W. J., and J. M. Wallace, 1998: The Arctic oscillation signature in the wintertime geopotential height and temperature fields. Geophys. Res. Lett., 25, 1297–1300, https://doi.org/10.1029/98GL00950.
Tian, F., B. Dong, J. Robson, and R. J. C. D. Sutton, 2018: Forced decadal changes in the East Asian summer monsoon: The roles of greenhouse gases and anthropogenic aerosols. Climate Dyn., 51, 3699–3715, https://doi.org/10.1007/S00382-018-4105-7.
Wang, B., R. Wu, and X. Fu, 2000: Pacific–East Asian teleconnection: How does ENSO affect East Asian climate? J. Climate, 13, 1517–1536, https://doi.org/10.1175/1520-0442(2000)013<1517:PEATHD>2.0.CO;2.
Wang, H. J., and S. P. He, 2012: Weakening relationship between East Asian winter monsoon and ENSO after mid-1970s. Chin. Sci. Bull., 57, 3535–3540, https://doi.org/10.1007/s11434-012-5285-x.
Wang, H. J., and K. Fan, 2013: Recent changes in the East Asian monsoon (in Chinese). Chin. J. Atmos. Sci., 37, 313–318.
Wang, H. J., and S. P. He, 2013: The increase of snowfall in Northeast China after mid-1980s (in Chinese). Chin. Sci. Bull., 58, 629–633, https://doi.org/10.1360/csb2013-58-8-629.
Wang, H. J., and H. P. Chen, 2016: Understanding the recent trend of haze pollution in eastern China: Roles of climate change. Atmos. Chem. Phys., 16, 4205–4211, https://doi.org/10.5194/acp-16-4205-2016.
Wang, S., and J. Liu, 2016: Delving into the relationship between autumn Arctic sea ice and central–eastern Eurasian winter climate. Atmos. Ocean. Sci. Lett., 9, 366–374, https://doi.org/10.1080/16742834.2016.1207482.
Wang, T., and J. Miao, 2018: Twentieth-century Pacific decadal oscillation simulated by CMIP5 coupled models. Atmos. Ocean. Sci. Lett., 11, 94–101, https://doi.org/10.1080/16742834.2017.1381548.
Wang, T., O. H. Ottera, Y. Q. Gao, and H. J. Wang, 2012: The response of the North Pacific decadal variability to strong tropical volcanic eruptions. Climate Dyn., 39, 2917–2936, https://doi.org/10.1007/s00382-012-1373-5.
Wang, T., H. J. Wang, O. H. Ottera, 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.
Wu, B., and R. Huang, 1999: Effects of the extremes in the North Atlantic Oscillation on East Asia winter monsoon (in Chinese). Chin. J. Atmos. Sci., 23, 641–651.
Wu, B., and J. Wang, 2002: Winter Arctic oscillation, Siberian high and East Asian winter monsoon. Geophys. Res. Lett., 29, 1897, https://doi.org/10.1029/2002GL015373.
Xiao, Z.-N., D.-L. Li, L.-M. Zhou, L. Zhao, and W.-J. Huo, 2017: Interdisciplinary studies of solar activity and climate change. Atmos. Oceanic Sci. Lett., 10, 325–328, https://doi.org/10.1080/16742834.2017.1321951.
Xu, J., Q. Zhu, and T. Zhou, 1999: Sudden and periodic changes of East Asian winter monsoon in the past century (in Chinese). Quart. J. Appl. Meteor., 10, 1–8.
Yan, X., R. Zhang, and T. R. Knutson, 2018: Underestimated AMOC variability and implications for AMV and predictability in CMIP models. Geophys. Res. Lett., 45, 4319–4328, https://doi.org/10.1029/2018GL077378.
Yang, Z., W. Huang, B. Wang, R. Chen, J. S. Wright, and W. Ma, 2018: Possible mechanisms for four regimes associated with cold events over East Asia. Climate Dyn., 51, 35–56, https://doi.org/10.1007/s00382-017-3905-5.
Yim, B. Y., M. Kwon, H. S. Min, and J.-S. Kug, 2015: Pacific decadal oscillation and its relation to the extratropical atmospheric variation in CMIP5. Climate Dyn., 44, 1521–1540, https://doi.org/10.1007/s00382-014-2349-4.
Zhang, R., A. Sumi, and M. Kimoto, 1996: Impact of El Niño on the East Asian monsoon: A diagnostic study of the ‘86/87 and ‘91/92 events. J. Meteor. Soc. Japan, 74, 49–62, https://doi.org/10.2151/JMSJ1965.74.1_49.
Zhou, W., 2017: Impact of Arctic amplification on East Asian winter climate. Atmos. Ocean. Sci. Lett., 10, 385–388, https://doi.org/10.1080/16742834.2017.1350093.
Zhou, W., X. Wang, T. J. Zhou, and J. C. L. Chan, 2007: Interdecadal variability of the relationship between the East Asian winter monsoon and ENSO. Meteor. Atmos. Phys., 98, 283–293, https://doi.org/10.1007/s00703-007-0263-6.
Zhu, Y., and X. Yang, 2003: Relationships between Pacific Decadal Oscillation (PDO) and climate variabilities in China (in Chinese). Acta Meteor. Sin., 61, 641–654.
Zuo, J., W. Li, and H. Ren, 2013: Representation of the Arctic Oscillation in the CMIP5 Models. Adv. Climate Change Res., 4, 242–249, https://doi.org/10.3724/SP.J.1248.2013.242.