HTML
--> --> --> -->3.1. Reversal of the SH between November and DJ
To capture the spatial phase relationships and the coherent variations between November and DJ SLP, firstly, the multivariate empirical orthogonal function (MV-EOF) analysis method (Wang, 1992) is employed. In this paper, a covariance matrix is constructed for the combined normalized November and DJ monthly averaged SLP to carry out the MV-EOF. The domain (10°-60°N, 60°-150°E), which includes most of Asia, is chosen for the analyses. The leading mode accounts for 26.42% of the total covariance for the November and DJ SLP together. The first mode exhibits an in-phase spatial pattern with consistent anomalies over most of central and eastern China and the entire tropics from November to DJ (figure omitted). (Chang and Lu, 2012) ascribed this mode to ENSO because of the significant correlations between this mode and the SST in the Ni?o3.4 region.This study primarily focuses on the second empirical mode, which presents different spatial patterns between November and DJ and accounts for 14.04% of the total covariance. The SLP exhibits inverse signs in November and DJ in the mid-high latitudes (Figs. 2a and b). In November, the maximum amplitude lies to the northwest of Lake Balkhash, and the anomalies extend eastward to Northeast Asia. By contrast, the pattern in DJ shows nearly opposite anomalies over the analyzed domain. The maximum anomalies appear to the southeast of Lake Baikal and extend to southeastern China, which is in accordance with the aforementioned cold-air outbreak over East Asia associated with the SH. The second empirical mode displays large interannual variability, which can be seen through the second principal component (PC2) (Fig. 2c).
Figure2. Spatial patterns of the second MV-EOF modes of (a) November and (b) DJ SLP. (c) Normalized time series of the second principal component (PC2). The region where the SHI is calculated is denoted by the black rectangular frame in panels (a, b).
The SHI is applied to quantitatively describe the intensity of the reversion zone in Figs. 2a and b. Accordingly, the SHI in November and DJ are respectively referred to as SHI_N and SHI_DJ in the following. To further illustrate the atmospheric differences between November and DJ over the Eurasian continent, the regressions of SLP/SAT anomalies on the SHI from 1979 to 2015 are shown in Fig. 3. When the SHI_N is positive, significant positive SLP anomalies prevail in November north of 40°N and extend southeast to eastern China, Japan and the Northwest Pacific (Fig. 3a). Maximum amplitude appears in northwestern Siberia and adjacent areas, indicating prominent resemblance to the MV_EOF2 pattern in November. Correspondingly, anomalously low SAT appears in November to the west of Lake Baikal and stretches westward to the Caspian Sea, northeastward to the Far East, and southeastward to southern China (Fig. 3b). Interestingly, in the following DJ, noticeable positive anomalies replace those positive SLP (negative SAT) anomalies in November, although the anomalous area shrinks slightly and moves a little southward (Figs. 3c and d). These results further confirm the reversal of signals associated with the SH between November and DJ. Besides the simultaneous reflection, the impacts of November SH on the lagged atmospheric circulation anomalies are explored. The regressions of uv850, Z500 and u300 in DJ on the SHI_N are presented in Fig. 4. With the occurrence of a positive November SH, a significant anomalous cyclone (anticyclone) appears over Siberia (the West Pacific) during the following DJ, when anomalous southerly wind prevails along coastal East Asia (Fig. 4a). At 500 hPa, a meridional dipole of geopotential height anomalies (Fig. 4b), with significant negative anomalies throughout Siberia and positive anomalies expanding from the West Pacific to the hinterland of Asia, indicates a weakening of both the Ural high and East Asian trough (He and Wang, 2013a). Meanwhile, the upper-tropospheric westerly jet stream decelerates over 30°-35°N, and a wide range of zonal wind accelerates to the north of the jet core (Fig. 4c). Thus, the meridional shear of zonal winds on the north side of the jet stream decreases, which results in a weaker EAWM and SH (Li and Yang, 2010). By contrast, when a weaker SH occurs in November, the EAWM-related circulation is more likely to be stronger in DJ. Overall, it is proved that the extratropical circulation processes in DJ over Siberia are negatively correlated with the preceding November SH.
Figure3. Regressions of November (a) SLP (units: hPa) and (b) SAT (units: K) anomalies on the SHI_N during 1979-2015. Panels (c, d) are similar to (a, b), respectively, except that they show the SLP and SAT in DJ. Dotted regions indicate the anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
Figure4. Normalized and detrended time series of the SIAI in September, along with the SHI_DJ and SHI_N during 1979-2015. The black (grey) dots in the bottom panel indicate the years when the SHI_N and SHI_DJ (SHI_N and SIAI) are in reverse phase, and the brown dots indicate the years when the SHI_DJ and SIAI are in the same phase.
In addition, Fig. 5 displays the time series of SHI_N and SHI_DJ during 1979-2015. The two indices barely keep the same phase with each other on the interannual timescale. During the 37-year period, more than half of the time (23 years) witnesses opposite relationships. To be specific, the correlation coefficient between SHI_N and SHI_DJ is -0.50 (above the 99% confidence level).
Figure5. Regressions of (a) wind at 850 hPa (units: m s-1), (b) geopotential height at 500 hPa (units: gpm), and (c) zonal wind at 300 hPa (units: m s-1) in DJ on the SHI_N. Bold arrows in (a) indicate either component of the wind anomalies is significant at the 90% confidence level. Dotted regions in (b, c) indicate anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
2
3.2. Possible physical mechanisms
The above analyses reveal the reversal of the SH between November and DJ. However, the underlying physical mechanism is not clear. As suggested by many previous studies (Thompson and Wallace, 2000; Wang and Ikeda, 2000; Gong et al., 2001), the SH is closely related to the AO. To examine the possible impact of the AO on the reversal of the SH, the November/DJ AO signal (obtained from the National Oceanic and Atmospheric Administration's website: https://www.esrl.noaa.gov/psd/data/correlation/ao.data) is removed from the November/DJ SLP via linear regression. The correlation coefficient between the November and DJ non-AO SHI turns out to be -0.56 (above the 99% confidence level), indicating that the sign change between the SHI_N and SHI_DJ is unlikely due to the effect of the AO.(Jung et al., 2014) suggested that the involvement of subpolar processes may lead to more skillful forecasts than the midlatitude and tropical dynamics over eastern Europe and northern Asia, at least for medium-range and extended-range forecasts. Thus, it is important to connect the sub-seasonal variations of the SH to the Arctic change. Given the potential connection, regression maps of September and October sea-ice concentration on the SHI_N and SHI_DJ are presented in Fig. 6. The sea-ice concentration anomalies show high coherence from September to October with regard to SHI_N. Significant positive sea-ice concentration anomalies appear in the Kara Sea and the Laptev Sea in September (Fig. 6a) and persist into October (Fig. 6c). Besides, zonally elongated negative sea-ice concentration anomalies arise along the Siberian coast in September and October with regard to the SHI_DJ (Figs. 6b and d). The indication is that, when the SIA is reduced in September and October, the SH tends to be weaker in the subsequent November and stronger in DJ. Therefore, the following analyses mainly focus on the Arctic sea ice to explore the possible mechanism of the reversal of the SH between November and DJ. Following most earlier studies (Francis et al., 2009; Honda et al., 2009; Overland and Wang, 2010; Wu et al., 2011), we choose the September SIAI as representative of the autumn sea-ice variability. The normalized and detrended September SIAI is plotted in Fig. 5 along with SHI_N and SHI_DJ. There is a significant negative correlation (-0.42) between the SIAI and SHI_N and a significant positive correlation (0.49) between the SIAI and SHI_DJ (both above the 99% confidence level). Moreover, during the 37 years of 1979-2015, 25 cases have opposite signs for the SIAI and SHI_N, while 23 cases have the same sign for the SIAI and SHI_DJ (the dots below in Fig. 5).
Figure6. Regressions of September sea-ice concentration (SIC) anomalies (units: %) on the (a) SHI_N and (b) SHI_DJ. Panels (c, d) are similar to (a, b), respectively, except they show the October SIC anomalies. Dotted regions indicate anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test. The region for the definition of the SIAI is denoted in panel (d).
To better understand the influence of the autumn sea ice on the atmospheric circulation in the following winter, a composite analysis of the November/DJ SLP and Z500 anomalies between light and heavy sea-ice years is performed, based on a standard deviation of 0.5. In November, significant negative SLP anomalies can be seen throughout Siberia, with maximum amplitude of more than -5 hPa (Fig. 7a). The positive SLP anomalies during DJ show a significant positive anomaly center to the west of Siberia (Fig. 7b). The Z500 anomalies in November depict a hemispheric wave pattern resembling the wave-3 structure (Fig. 7c). By contrast, the 500 hPa anomalies display a dipole pattern over the Arctic-Eurasian sector in DJ (Fig. 7d). Considering the light sea ice in autumn, the DJ Z500 shows significant positive responses over northern Siberia and the adjacent Arctic sea, and negative responses over the midlatitude Asian continent, which bears some resemblance to the negative phase of the AO over the Arctic-Eurasian area (Tang et al., 2013; Mori et al., 2014).
Figure7. Composite differences of SLP (units: hPa) in (a) November and (b) DJ between the years of positive and negative September SIAI. Panels (c, d) are similar to (a, b), respectively, except they show the Z500 anomalies (units: gpm). Dotted regions indicate anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
The coupled interactions between the Arctic sea ice and the SHF, as well as the associated circulation anomalies in November and DJ are presented in Figs. 8a and b. Positive (negative) values of SHF are indicative of upward (downward) heat fluxes. The regressions of SHF are characterized by a sharp contrast over the Barents Sea. It is familiar that the reduction of sea ice (which means more open water) corresponds to increased ocean-to-atmosphere SHF (Honda et al., 2009). As open water has a lower albedo than ice, the surface ocean absorbs more radiative energy during summertime (Perovich et al., 2007), and the growth of sea ice in November is hindered (Francis et al., 2009). Consequently, during November, the absorbed heat is transferred upward via turbulent fluxes, and the SHF exhibits significant positive anomalies in the Barents Sea and western Kara Sea (Fig. 8a). However, reversal appears in DJ with significant negative SHF anomalies (Fig. 8b). The downward SHF in absolute terms is larger in magnitude (~40 W m-2) than the November upward SHF (~30 W m-2). Results derived from ERA-Interim data (figure omitted) are consistent with those displayed in Fig. 8. (Sorokina et al., 2016) demonstrated that the leading pattern of winter (here, December-February) interannual variability of turbulent heat flux over the Barents Sea is weakly correlated with sea ice. The ocean does no longer plays a dominant role in the ocean-atmosphere interaction; that is, the SHF anomalies are mainly driven by atmospheric changes in DJ. The warming of the overlying atmosphere due to Arctic amplification (Serreze and Barry, 2011) reduces the sea-air temperature differences, and thus the SHF is reduced (Sorokina et al., 2016).
Figure8. Regressions of November (a) SHF (shading; units: W m-2) and 850 hPa winds (vectors; units: m s-1), and (c) snow depth (SD) anomalies (units: cm), on the SIAI. Positive SHF anomalies indicate anomalous upward fluxes. Panels (b, d) are similar to (a, c), respectively, except for DJ. Winds where either component of the wind anomaly exceeds 0.3 m s-1 are displayed. Dotted regions indicate anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
Based on the above observational analyses, a physical mechanism that relates the autumn Arctic sea-ice loss to the reversal of the winter SH is investigated. Many previous studies have pointed out that the interaction between the storm track activity (also called the synoptic-scale eddy) and the low-frequency flow plays a key role in the variability of the atmospheric circulation (Lau and Nath, 1991; Hartmann and Lo, 1998; Hoskins, 2001; Grise et al., 2013; Lehmann and Coumou, 2015). (Lau, 1988) demonstrated that enhanced storm track activity is accompanied by westerly wind accelerating locally, as well as cyclonic vorticity forcing to the north and anticyclonic vorticity forcing to the south, and vice versa. The modification of storm tracks can influence not only nearby areas but also downstream atmospheric circulation through synoptic transient eddy forcing (Tyrlis and Hoskins, 2008). In short-term climate change, the extratropical storm track plays an important role through transporting heat, moisture and momentum (Seierstad and Bader, 2008). Moreover, the increased temperature, moisture and the latent heating can enhance the eddy energy and storm track activity in a warmer climate (Hall et al., 1994). Previous analyses suggest that the change in surface heat flux over the Barents Sea is huge, so it is reasonable to associate sea-ice anomalies with storm tracks, as well as downstream intraseasonal circulation anomalies, which will be discussed next.
Figure 9 presents the regressions of storm track intensity and blocking frequency on the SIAI. The storm track activity is calculated as the variance of the 2.5-6-day bandpass-filtered Z250 [$\overline{(z')^2}$, where z' denotes the synoptic-scale geopotential height anomalies and the overbar indicates the time mean in November/DJ] (Lau and Nath, 1991). A blocking high event is defined if the daily Z500 exceeds one standard deviation above the monthly mean for each grid cell and persists for at least five consecutive days (Liu et al., 2012). The local blocking frequency for November (DJ) is measured by the ratio between the number of blocking days and the total number of days in November (DJ). Associated with reduced September SIA, increasing storm tracks appear over the Barents Sea coast and northeastern Europe in November (Fig. 9a) and significantly decreased storm tracks extend from eastern Europe southeastward to western Asia in DJ (Fig. 9b). Decreased blocking frequency is seen over the Urals in November (Fig. 9c) and, conversely, more extended increased blockings are seen over the Urals and the Barents Sea in DJ (Fig. 9d), which is consistent with the circulation anomalies shown in Figs. 7 and 8. The zonal mean anomalies over the Eurasian region (20°-140°E) are further depicted in Fig. 10. In November, a smaller SIA favors the release of heat and moisture from the Arctic Ocean. Therefore, the static stability in the lower troposphere will decline and the baroclinicity aloft will be enhanced, which then intensify the storm track activity over the ocean and adjacent European continent (Figs. 9a and 10a). The wave activity propagates upward and equatorward at the latitude of the Arctic Ocean (65°-80°N), and then it propagates downward and diverges at 50°-65°N, which could result in the acceleration of westerly wind (Andrews et al., 1987). According to previous studies, positive storm track anomalies are accompanied by westerly wind anomalies and cyclonic vorticity forcing to their north (Lau, 1988; Gong et al., 2011; Chen et al., 2017). Correspondingly, large-scale cyclonic circulations prevail over high-latitude Eurasia (Fig. 8a). In DJ, in contrast, the upward SHF tightly related to the sea-ice loss disappears, with decreased heat and moisture supplies for the adjacent European region. Moreover, the meridional temperature gradient decreases because of the effect of the Arctic warming. The weakened meridional temperature gradient restrains the growth of baroclinic eddies and decelerates the westerlies (Fig. 10b). In this situation, the input of the westerly angular momentum is weakened over high-latitude Eurasia and less eddy energy is required to maintain the momentum balance (Held, 1993). Unsurprisingly, the transportation of wave activity from the Arctic to the region between 50°N and 65°N, which exists in November, disappears in DJ. In this latitudinal zone, in turn, a net loss of energy can be seen as the WAF diverges to lower latitudes. Finally, the intensity of storm tracks decreases in DJ. To the north of the negative storm track anomalies, anticyclonic circulation anomalies replace the cyclonic winds prevailing in November (Fig. 8b). Notably, the distinct differences between November and DJ exist mostly in the troposphere.
Figure9. Regressions of November (a) mean storm track intensity (units: dgpm2) and (c) blocking frequency (units: %) on the SIAI. Panels (b, d) are similar to (a, c), respectively, except for DJ. Dotted regions indicate anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
Figure10. Regressions of zonal wind (shading; units: m s-1), WAF (vectors; units: m2 s-2) and storm track intensity (black contours; units: dgpm2) averaged between 20°E and 140°E in (a) November and (b) DJ on the SIAI. The vertical component of the WAF is multiplied by 125 (Castanheira and Graf, 2003). Red and blue contours indicate zonal wind anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
Finally, the question as to the way that the SH responds to the variability of upstream storm tracks is resolved. Many studies have focused on this issue. For instance, (Joung and Hitchman, 1982) revealed a clear sequence whereby the surface synoptic situation over East Asia during wintertime is characterized by an intense eastward-moving perturbation beginning six or seven days in advance from the western North Atlantic across Eurasia. The synoptic-scale perturbations over Siberia share similar characteristics to the North Atlantic storm tracks (Hoskins and Hodges, 2002). By carrying out model experiments, (Magnusdottir et al., 2004) confirmed that the December-March mean cold anomaly response appears in eastern Siberia when a weaker, southward-shifted, and more zonal storm track is found over the North Atlantic sector. More detailed mechanisms are documented as follows. When the North Atlantic-Europe storm tracks are active, the synoptic transient waves extending from the Atlantic to Europe intensify the anticyclonic wave breaking (Rivière and Orlanski, 2007; Zeng et al., 2015). Through the feedback of the transient waves on the monthly mean flow, the circumpolar eddy-driven jet is strengthened, and thus the southward movement of Arctic cold air is obstructed. The background field is disadvantageous for the maintenance of the Ural blocking high. As a result, the air over Siberia and East Asia becomes warmer and the SH becomes weaker, which shrinks and moves southward. The situation for DJ is basically reversed, and a stronger SH occurs. All these studies shed light on the impacts of upstream storm tracks on the SH and are capable of proving that, under the condition of light autumn Arctic sea ice, the increased (decreased) storm track activity favors a weaker (stronger) SH in November (DJ).
In addition, the possible role of Siberian snow in the evolution of the SH during November and DJ is also investigated. As the Arctic sea ice is reduced, higher snow depth anomalies can be seen over most of the Eurasian land in DJ, especially in Siberia and northern China (Fig. 8d). The presence of an anomalously thick snow cover leads to near-surface cooling due to the high albedo of snow, which favors a strengthened and more expansive SH (Cohen and Entekhabi, 1999). However, there are barely any significant snow depth anomalies associated with low SIA in November (Fig. 8c). In particular, the northerly cold air advection along the eastern and southern flanks of the enhanced Ural high favors the thicker snowpack in DJ. Conversely, the warmer air advection from lower latitudes in November barely favor effective snow accumulation locally. Therefore, the large snowfall in DJ could further intensify the SH in DJ and enlarge the discrepancy of the SH in November and DJ.
2
3.3. Validation by model simulations
To test whether coupled models can represent the reversed response of the SH to autumn Arctic sea ice, we use the simulations (historical experiments) of five models [ACCESS1.0, ACCESS1.3, CESM1(CAM5), HadGEM2-CC, and GFDL CM2.1] from phase 5 of the Coupled Model Intercomparison Project (CMIP5) for the period 1979-2005. The selection of the CMIP5 models is primarily driven by their capability to reproduce the variability of observed sea-ice extent, SLP and uv850 fields since 1979 (e.g., Gleckler et al., 2008; Sperber et al., 2012; Liu et al., 2013). The available variables used include the monthly SLP, uv850, Z500, sea-ice concentration, and surface snow amount. Before performing the analyses, we interpolate all the model simulation data to a horizontal resolution of 1°. Here, the ensemble mean results of the five model simulations are provided.First, we calculate the correlation coefficient between SHI_N and SHI_DJ, which is -0.44 (above the 95% confidence level), based on the CMIP5 simulations during 1979-2004. Figure 11 illustrates the regression maps of the September and October ensemble-mean sea-ice concentration on the SHI_N and SHI_DJ. As expected, significant sea-ice anomalies persist from September to October along the Siberian coast with regard to either SHI_N or SHI_DJ, whereas an apparent reversal can be seen between the regressions upon SHI_N and SHI_DJ. The responses of Z500 and uv850 to the reduced Arctic sea ice are also shown in Fig. 12. The negative Z500 anomalies to the east of the Urals and the coherent cyclonic winds can be seen in November (Figs. 12a and c). Although the CMIP5-simulated anomalous amplitudes are smaller than observed, it is evident that the models simulate the dipole pattern of Z500 over the Arctic-Eurasian sector and anticyclonic winds over northern Europe in DJ (Figs. 12b and d), consistent with Figs. 7d and 8b. Moreover, no robust linkage between autumn sea ice and November snowfall in Siberia is found, but a significant increase in snow cover can be seen in DJ (Figs. 12c and d). To some extent, the physical linkages between autumn Arctic sea ice and the reversed SH in November and DJ can be reproduced by the historical simulations.
Figure11. As in Fig. 6, except for CMIP5 historical simulations (1979-2004).
Figure12. Composite differences of Z500 anomalies (units: gpm) in (a) November and (b) DJ between the years of positive and negative SIAI; regressions of (c) November and (d) DJ surface snow amount (SA) (shading; units: kg m-2) and 850 hPa winds (vectors; units: m s-1) anomalies on the SIAI, during 1979-2004, derived from the CMIP5 historical simulations. Winds where either component of the wind anomaly exceeds 0.15 m s-1 are displayed. Dotted regions indicate anomalies significant at the 90% confidence level, based on a two-tailed Student's t-test.
Figure13. Schematic diagram illustrating the mechanisms for the low autumn Arctic sea ice leading to the reversal of the SH in November and DJ.