HTML
--> --> -->Previous studies have shown that sea surface temperature (SST) variability in the tropical Pacific (Wang et al., 2000; He and Wang, 2013a; Jin et al., 2016), North Pacific (He and Wang, 2013b; Ye and Chen, 2016), and Indian Ocean (Yang et al., 2010) is closely associated with the East Asian winter monsoon (EAWM) system. (Wang et al., 2000) once revealed that El Ni?o events cause a weaker than normal EAWM and predominantly warm conditions over East Asia through a Pacific-East Asian teleconnection pattern. However, recent studies have demonstrated that the interannual relationship between the El Ni?o-Southern Oscillation (ENSO) and EAWM clearly undergoes a low-frequency oscillation and weakened around the mid-1970s (Wang and He, 2012; He and Wang, 2013a; He and Wang, 2013a). Furthermore, it has been noted that the Pacific Decadal Oscillation (PDO) exerts an important modulation effect on the ENSO-EAWM relationship on low-frequency timescales; the negative correlation between ENSO and EAWM is remarkably enhanced when they are in phase (Kim et al., 2014b). For the East Asian summer monsoon (EASM) system, a decadal variation in the EASM-ENSO connection occurred in the mid-1990s (Yim et al., 2008), and the decadal shift of the summer rainfall variability over southern China after the mid-1990s is also related to the tropical Pacific SST distribution through wave-like trains from both the western Pacific and Eurasia to East Asia (Chang et al., 2014).
In addition to the SST contribution, East Asian winter climate is greatly influenced by changes in Arctic sea-ice cover (Honda et al., 2009; Wu et al., 2011; Liu et al., 2012; Li and Wang, 2014; Li et al., 2015c; Wang et al., 2015b; Wang and Liu, 2016) and Eurasian snow (Gong et al., 2003; Cohen et al., 2007; Xin et al., 2014; Xu et al., 2017) during autumn and winter. (Yeo et al., 2014) indicated that the recent climate variability in the Bering and Chukchi Seas induced by recent diminishing sea-ice cover is strongly connected with large-scale climate patterns in the Pacific. Sea-ice reduction in the Barents-Kara Seas is responsible for recent severe winters through eliciting vertical propagation of planetary waves and the subsequent weakening of the stratospheric polar vortex (Kim et al., 2014a) and downstream response of cyclonic activity (Tang et al., 2013). (Liu et al., 2012) revealed that declining autumn Arctic sea ice induces much broader meridional meanders of atmospheric circulations at midlatitudes and more frequent blocking episodes that generate heavy snowfall in the European winter. Eurasian snow cover has been widely considered as an important predictor of East Asian winter climate due to its significant impacts on atmospheric circulations through regional or large-scale radiative budgets (Jeong et al., 2011; Furtado et al., 2015). High October Eurasian snow cover extent is generally accompanied by a weakened polar vortex, a negative Arctic Oscillation (AO), and cold spells in winter through stationary planetary wave activities (Cohen et al., 2007, 2012). From the perspective of vertical-propagating planetary waves with zonal wavenumber-1, however, (Xu et al., 2017) argued the influence of reduced October northern Eurasian snow cover on the "warm Arctic-cold Eurasia" pattern in the following January.
Early studies mostly focused on seasonal mean climate anomalies, paying little attention to subseasonal variability (Wang et al., 2000; Cohen et al., 2012; Kug et al., 2015). However, observational evidence shows that a reverse variability in winter climate on the subseasonal timescale sometimes occurs over East Asia, resulting in huge economic losses and severe social impacts. For example, in January 2015, the surface air temperature (SAT) in China reached another historical high record (Wang et al., 2015a), while cold conditions prevailed in December 2014, with the amplitude of temperature anomalies reaching approximately 5°C. Furthermore, as part of the present work we noticed that in winter 2014/15 the Pacific SST anomalies resembled El Ni?o conditions in the central and eastern tropical Pacific, along with a positive PDO phase in the extratropical North Pacific. Meanwhile, a dramatic reduction in sea-ice cover over the Laptev-East Siberian Sea occurred in autumn 2014. The opposite sign in December and January-February SAT anomalies in the same context of winter Pacific SST and autumn Arctic sea-ice anomalies motivated us to consider that the effect of winter tropical Pacific SST variability (autumn Arctic sea-ice changes) might be different on the subseasonal timescale under the positive phase of PDO modulation. Therefore, in this study we attempt to explore how the winter warm Ni?o4 phase (autumn Arctic sea-ice reduction) exerted an influence on the East Asian SAT in the early and late winter of 2014/15 under the modulation of a positive PDO phase, an understanding of which might be helpful for climate prediction on the subseasonal timescale.
The December, January and February PDO indices are obtained from http://research.jisao.washington.edu/pdo/. The December, January and February Ni?o4 indices are defined as the area-averaged SST anomalies in (5°S-5°N, 160°E-30°W). The September-October Laptev-East Siberian sea-ice cover (LE-SIC) index is defined as the area-averaged sea-ice cover anomalies for (72°-82°N, 90°-180°E) (black frame in Fig. 4d). All the defined indices are standardized. To emphasize the interannual variability, the linear trend is removed from all data and indices prior to the correlation and composite analyses. The correlation coefficients between the Ni?o4 and September-October LE-SIC indices are 0.07 in December, 0.14 in January, and 0.16 in February, suggesting that the impacts of winter Ni?o4 SST anomalies and autumn LE-SIC reduction on the East Asian winter climate are independent of each other.
To investigate the interdecadal modulation of the positive PDO phase on winter warm Ni?o4 SST anomalies (autumn Arctic sea-ice reduction), the PDO index is applied by a 5-yr running mean to extract a decadal variability of PDO. High and low Ni?o4 (LE-SIC) cases are determined when the Ni?o4 (LE-SIC) index is above and below standard deviations of 0.5 and -0.5, respectively. Positive and negative PDO (PDO\(\vert\)+ and PDO\(\vert\)-) phases correspond to cases in which the 5-yr running mean PDO index is above and below zero, respectively. The classification of high and low Ni?o4 (LE-SIC) years according to the different phases of PDO is shown in Table 2a (2b). We categorize high PDO conditions when the 5-yr running mean PDO index is above a standard deviation of 0.5 (Table 2c).

Figure 2 is the same as Fig. 1, but displays SLP, UV850, SH indices, and EAWM indices (defined in the following analysis). In December 2014, positive SLP anomalies are located in (80°-120°N, 40°-60°E), where the SH domain is (Wu and Wang, 2002), suggesting a strengthened SH; negative values are observed over Japan and the eastern North Pacific (Fig. 2a; shaded). In the wind field, anomalous northerly surface winds prevail along the eastern flank of the SH (Fig. 2a; vectors), indicating that the prevailing northerlies that convey cold air from high latitudes to East Asia are strengthened (He and Wang, 2012). In contrast, the SH is suppressed in January 2015 because of the negative SLP anomalies over the Asian continent and the eastern North Pacific (Fig. 2b; shaded), accompanied by anomalous southerly winds along the eastern flank (i.e., the weakened EAWM; Fig. 2b; vectors). The SLP pattern in February is similar to that in January, and mid- to high-latitude Asia is occupied by anomalous southwesterly and southeasterly winds (Fig. 2c; vectors). We employ an SH index defined by the area-averaged SLP of (40°-60°N, 80°-120°E) (Wu and Wang, 2002) and an EAWM index defined by the area-averaged 850-hPa wind speed of (25°-50°N, 115°-140°E) (Wang and Jiang, 2004). The SH indices in December, January and February are 0.7, -1.3 and -0.6, respectively, and the EAWM indices are 2.3, -0.3 and -0.7, respectively (Fig. 2d). That is, the low-level atmospheric circulation anomalies in the Asian-western North Pacific region reverse on the subseasonal timescale, consistent with the East Asian SAT variability. Furthermore, the East Asian SAT indices and SH indices are highly correlated, with coefficients of -0.77 in December, -0.85 in January, and -0.81 in February (above the 99% confidence level) for the period 1979/80-2014/15.
The subseasonal U200 climatology during 1979/80-2014/15, U200 anomalies in winter 2014/15 relative to the average of winters 1979/80-2014/15, and the corresponding EAJS indices [defined by (Li and Yang, 2010)] are provided in Fig. 3. Positive U200 anomalies extend from the North Pacific to northern China in December (Fig. 3a; shaded), suggesting an acceleration and westward extension of the EAJS that is usually centered near southern Japan in the wintertime (Figs. 3a-c; contours). In January-February 2015, the EAJS decelerates as negative U200 anomalies are observed in China, Japan, and the western North Pacific (Figs. 3b and c; shaded). The corresponding EAJS indices in December, January and February are 2.4, -1.1 and -1.4, respectively (Fig. 3d). Concomitant with low-level circulation variations, the upper-level circulation anomalies are also of opposite sign in December 2014 and January-February 2015. Taken together, in winter 2014/15, December East Asian cold conditions are associated with the strengthened SH, enhanced surface northerlies, and an accelerated and westward-extended EAJS; warm anomalies in January-February are also closely related to the suppressed SH, weakened low-level northerlies, and a decelerated EAJS. (Wang et al., 2010) suggested that the EAWM includes southern and northern modes. Interestingly, the time series of the southern mode also indicated that the EAWM experienced subseasonal diversity in winter 2014/15, with a negative temperature anomaly in December and positive temperature anomaly in January and February (figures not shown). This suggests that the subseasonal reversal of East Asian SAT variability might be related to the southern mode.



2
4.1. How do the PDO\(\vert \)+ phase and Ni?o4 SST anomalies influence January-February SAT variability?
The left-hand panel of Fig. 5 illustrates the composite subseasonal zonal-mean mass stream function anomalies between high and low Ni?o4 cases (based on a standard deviation of 0.5; Table 2a) under a PDO\(\vert\)+ phase. As we can see, the Hadley cell in the Northern Hemisphere is relatively weak in December, with an insignificant positive anomaly center in the low troposphere in the tropical region (0°-15°N) (Fig. 5a). In January-February, the Hadley cell develops along with a significant Ferrel cell: the positive anomaly center in the tropical region extends vertically into the upper troposphere, with a negative anomaly center in the low- to mid-troposphere in the subtropical region (15°-35°N) (Figs. 5b and c). It has been documented that El Ni?o is linked to an intensified Hadley cell through surface heating at the equator (Li et al., 2015b); in this sense, the enhanced Hadley cell and Ferrel cell in January-February may represent a subseasonal delay of tropical surface heating, i.e., a stronger coupling between the warm Ni?o4 phase and extratropical circulations in January-February relative to the weaker coupling in December. The right-hand panel of Fig. 5 displays the evolution of composite SST anomalies (150°E-160°W mean) and U300 anomalies (60°E-120°W mean) between December high and low Ni?o4 cases during the PDO\(\vert\)+ phase. Statistically significant negative SST anomalies develop in the subtropics (20°N-35°N) in December, which lessen the north-south temperature gradient, thereby suppressing the development of the subtropical westerly jet (Li et al., 2015b), which shifts slightly southward (15°N-25°N) in January-February (Fig. 5d). As a result, under the regulation of the PDO\(\vert \)+ phase, there is less acceleration of the Northern Hemisphere westerly jet in December relative to January-February (Fig. 5e).
We further present the subseasonal U200 anomalies between high and low Ni?o4 cases under a PDO\(\vert\)+ phase in Fig. 6. The common characteristics exhibited in January-February are the significant "positive-negative-positive" U200 anomaly bands from lower (south to 10°N) to higher (north to 70°N) latitudes over East Asia, indicative of a decelerated EAJS and stronger zonal wind speed to the south and north of the EAJS. What distinguishes December are less significant negative values in northern China-southern Japan, within the EAJS domain (Yang et al., 2002). As (Zhang et al., 1997) once suggested an abnormally weak EAJS during El Ni?o events, the difference of the EAJS in December and January-February further supports the idea of a delayed linkage between warm Ni?o4 SST anomalies and extratropical atmospheric circulations under the modulation of a PDO\(\vert\)+ phase, which disappears under a PDO\(\vert\)- phase (figure not shown). According to earlier studies (Jhun and Lee, 2004; He and Wang, 2013b), such a significant upper-tropospheric "positive-negative-positive" anomaly pattern in January-February may lead to a weaker meridional shear of the EAJS and weaker surface northerlies in January-February.

The next diagnostic is generated by presenting the corresponding 1000-hPa stream function and SAT (Fig. 7). A notable discrepancy over Siberia is the anomalous cyclonic circulation observed only in January-February, which manifests the subseasonal variability of the SH and a weakened one in January-February (Figs. 7a-c). The anticyclone in the western North Pacific also shows a subseasonal variability (Figs. 7a-c). In the SAT field, it is apparent that statistically significant warm anomalies dominate over most of East Asia in January-February, while fewer significant values occur in December (Figs. 7d-f). To conclude, under the modulation of a PDO\(\vert\)+ phase, the subseasonal delay of tropical surface heating associated with warm Ni?o4 SST anomalies may result in strong Hadley and Ferrel cells in January-February, which are responsible for the coupling between the tropical and extratropical regions. Consequently, the decelerated EAJS coincides with the weakened SH and warm conditions over East Asia in January-February.

2
4.2. How do the PDO\(\vert\)+ phase and Laptev-East Siberian sea-ice reduction influence December SAT variability?
Figure 8 illustrates the evolution of composite SST anomalies (120°E-180° mean) and U300 anomalies (60°E-120°W mean) under high PDO conditions in December (based on a standard deviation of 0.5; Table 2c) (top panel), and between December low and high LE-SIC cases under the PDO\(\vert\)+ phase (based on a standard deviation of 0.5; Table 2b) (bottom panel), respectively. Under high PDO conditions, cold SST anomalies in the western North Pacific persist from November to the following March and become the most significant in December (Fig. 8a). The response of atmospheric circulation to such strong SST variability is characterized by an accelerated westerly jet in December through an enlarged meridional temperature gradient (Fig. 8b). In addition, under the coincidence of a PDO\(\vert\)+ phase and September-October LE-SIC declines, cold SST anomalies remain in the western North Pacific in winter, accompanied by a statistically enhanced westerly jet in December (Figs. 8c and d). That is, the relationship between sea-ice loss and the westerly jet might be modulated by a PDO\(\vert\)+ phase on the subseasonal timescale.To clarify the regulation from the PDO\(\vert\)+ phase, we display the composite subseasonal U200 anomalies (left-hand panel) and the GPH500 and horizontal wave activity flux anomalies [computed according to the equation from (Takaya and Nakamura, 2001)] (right-hand panel) between low and high LE-SIC cases under the PDO\(\vert\)+ phase in Fig. 9. The "negative-positive-negative" U200 anomaly structure from the lower (south to 10°N) to higher latitudes (north to 60°N) over East Asia occurs only in December (Fig. 9a), which is consistent with Fig. 8d, and indicates an accelerated and westward-shifted EAJS and weaker zonal wind speed to the south and north of it. As depicted in (Li et al., 2014), the reduction in autumn sea-ice cover is responsible for the westward penetration of the EAJS, exciting the rearrangement of eastward-propagating Rossby waves with a much wider horizontal structure. In this sense, the acceleration and westward shift of the December EAJS may be conducive to the propagation of Rossby waves induced by the decreased sea-ice cover through wave-mean flow interaction (Wallace, 2000; Honda et al., 2009). As shown in Fig. 9d, the GPH500 response displays a remarkable wave pattern, with "negative-positive-negative" anomaly centers in Europe, central Asia, and East Asia (contours); a stationary Rossby wave train that propagates eastward from Europe to East Asia can also be estimated from the divergence of wave activity flux (vectors). The configurations of the EAJS and wave activities suggest robust polar-extratropical coupling and a strong influence of reduced autumn sea-ice cover on the East Asian climate in December during the PDO\(\vert\)+ phase. Under the PDO\(\vert\)- phase, the connection between changes in sea-ice cover and EAJS/wave activities becomes much weaker (figure not shown).
For the corresponding low-level circulation and SAT variability, in December, the positive SLP anomalies in the SH domain and negative SLP anomalies over the Sea of Japan increase the pressure gradient in the East Asia-western North Pacific region, and favor anomalous northerly winds and significant low temperatures over East Asia (Figs. 10a and d). In January-February, the anomalous SLP and SAT are quantitatively smaller and become insignificant (Figs. 10b and c, 10 and f). Taken together, under a PDO\(\vert\)+ phase, autumn LE-SIC reduction is associated with the westward shift of the EAJS and subsequently the eastward propagation of Rossby waves in December, which is responsible for the coupling between the polar and extratropical regions. As a result, the accelerated and westward-shifted EAJS concurs with a deepened East Asian trough (Fig. 9d), strengthened SH and significant cold conditions over East Asia in December.



In this section, we use the CMIP5 historical simulations to further test our hypothesis. We select six (five) models from CMIP5 (see Table 1), in which the response of winter SAT variations over East Asia to warm Ni?o4 phase (autumn LE-SIC reduction) under the PDO\(\vert\)+ phase for the period of 1979/80-2004/05 are successfully captured. When the winter warm Ni?o4 phase (autumn LE-SIC reduction) occurs with the PDO\(\vert\)+ phase, the spatial correlation coefficient of SAT in the region of (20°-80°N, 60°-180°E) between the reanalysis data and the six (five) model simulations is 0.43 (0.24).
Figure 11 illustrates the ensemble-mean SAT anomalies composited with the LE-SIC and Ni?o4 indices, respectively, during the PDO\(\vert\)+ phase. As expected, for the LE-SIC reduction, in December only, the simulated SAT anomalies show negative values from the east of Lake Baikal to coastal East Asia (Fig. 11a). For the warm Ni?o4 phase, positive SAT anomalies are located over mid- to high-latitude Eurasia in January, although the values in February are less significant (Figs. 11e and f). To some extent, the historical simulations are able to reproduce the subseasonal reversal in winter SAT variability over East Asia (cold December and warm January-February) based on the winter warm Ni?o-4 SST anomalies and autumn LE-SIC reduction under a PDO\(\vert\)+ phase.

Earlier studies indicated that the influence of the AO on the East Asian SAT can be regulated by the phase of the WP (Park and Ahn, 2016). Specifically, the AO-SAT relationship is strengthened (weakened) when the AO and WP are in-phase (out-of-phase), as a result of the significant (insignificant) zonal wavenumber-2 pattern of the EAWM-related atmospheric circulation. In addition to the decadal modulation of PDO (Kim et al., 2014b), the ENSO-EAWM relationship also varies depending on the AMO phase. La Ni?a events coincide with a strengthened (weakened) EAWM for a positive (negative) AMO, but the El Ni?o-EAWM relationship is irrelevant to the AMO phase (Geng et al., 2017). (He and Wang, 2013a) also attributed the low-frequency oscillation of the ENSO-EAWM relationship to a combination of PDO and AMO. In response to solar forcing, the ENSO-related winter climate anomaly varies with the 11-yr solar cycle and becomes significant during low solar activity winters (Zhou et al., 2013; Huo and Xiao, 2016). For January 2016, significantly colder conditions occurred over East Asia under the combined effect of the super El Ni?o and extreme Arctic warming, an effect that has been further verified using statistical prediction models (He et al., 2016). Therefore, the co-variability of impacts from different factors on the East Asian winter climate should be investigated in more detail in further work.
For other seasons and regions, a combined effect of different factors on the climate also exists. For example, the EASM exhibits strong interannual variation in response to ENSO (Wang et al., 2008); however, EASM variability can also be influenced by Tibetan Plateau diabatic heating through the western Pacific subtropical high (Zhang et al., 2004). Considering the forcing from the Tibetan Plateau, the EASM-ENSO relationship is highly correlated only when the Tibetan Plateau snow cover is reduced in summer (Wu et al., 2012). (Lee et al., 2015) recently attributed the extreme 2013/14 winter circulation over North America to the compounding effects of the warm SST anomalies in the tropical Pacific and extratropical North Pacific and the diminished sea-ice cover in the Arctic.
In terms of PDO modulation, it is suggested that the decadal change in the intensity of the interannual variability of the South China sea summer monsoon in the 20th century (Fan and Fan, 2017) and the predictability of the interannual variation in rainfall during early summer in southern China (Duan et al., 2013) can be modulated by PDO, which are higher during a PDO\(\vert\)+ phase than during a PDO\(\vert\)- phase. (Yu et al., 2015) also proposed a modulation of the "southern flood and northern drought" pattern by PDO in eastern China summers using 600-yr control simulations. The observed recent shift to a PDO\(\vert\)+ phase (Screen and Francis, 2016) motivated us to examine how the PDO\(\vert\)+ phase modulates the influences of warm Ni?o4 SST anomalies and Arctic sea-ice reduction on the winter SAT variability over East Asia on the subseasonal time scale. As depicted in (Li et al., 2015a), however, extratropical ocean warming plays an important role in changes in the winter circulation in the Arctic. Much more research is needed on the interdecadal modulation by the PDO\(\vert\)- phase.
(1) Under the modulation of a PDO\(\vert\)+ phase, warm Ni?o4 SST anomalies can induce a subseasonal delay of tropical surface heating, and hence a strong January-February Hadley cell and Ferrel cell, which are responsible for the coupling between tropical SST anomalies and extratropical atmospheric circulation. The resulting "positive-negative-positive" anomaly pattern in U200 suggests weakened meridional shear of the EAJS, and favors a situation in which the northerlies and SH are weakened. Therefore, significant warm conditions occur over East Asia in January-February.
(2) Under a PDO\(\vert\)+ phase, September-October LE-SIC reduction is related to significant cold SST anomalies in the western North Pacific and the subsequently larger meridional temperature gradient between the tropical and midlatitude regions. The EAJS is thereby accelerated and westward-shifted in December, and favors the eastward propagation of Rossby waves induced by sea-ice declines that strengthen the coupling between the polar and extratropical regions. Therefore, the SH and East Asian trough are intensified and East Asia becomes colder in December.
In summary, the subseasonal reversal of East Asian SAT variability in winter of 2014/15 may be attributable to winter warm Ni?o4 SST anomalies and autumn Arctic sea-ice reduction under the PDO\(\vert\)+ phase modulation. Moreover, such a "reversal" phenomenon has sometimes occurred under similar SST anomaly patterns over the past century, such as the winters of 1913/14 and 1965/66 (figure not shown).