HTML
--> --> -->A mounting body of observational and modeling evidence converges on a common conclusion that the occurrence of European cold winters becomes high when the solar activity reaches its minimum on the decadal time scale (Lockwood et al., 2010; Woollings et al., 2010; Ineson et al., 2011; Matthes, 2011). The surface pressure pattern in response to the 11-year solar cycle projects onto the North Atlantic Oscillation (NAO) (Kodera, 2003; Thiéblemont et al., 2015), a dominant mode of the Northern Hemisphere atmospheric variability influencing winter weather and climate in the Eurasia?Atlantic sector (Wallace, 2000; Bader et al., 2011; Matthes, 2011; Cohen et al., 2014; Screen, 2017a). The physical link between the solar ultraviolet (UV) forcing and the decadal NAO variability has been suggested to arise from the “top-down” mechanism (Kodera and Kuroda, 2002; Gray et al., 2010). A lagged and amplified NAO response to the solar cycle is further proposed involving the ocean?atmosphere interaction mechanism (Gray et al., 2013; Scaife et al., 2013). Thiéblemont et al. (2015) revealed a 1?2-year lagged solar?NAO linkage via a multidecadal experiment with 11-year solar forcing variability. Sunspot number, open solar flux, and F10.7 cm radio flux, representing the decadal solar variability or solar UV variations, have been widely used to investigate the solar contribution to Earth’s climate on the decadal time scale (Roy and Haigh, 2010; Woollings et al., 2010; Thiéblemont et al., 2015; Miao et al., 2018). However, the potential connections on the interannual time scale remain unclear.
In addition to radiative forcing, another important aspect of the solar activity influencing Earth’s climate is energetic particle forcing (Lilensten et al., 2015). Matthes et al. (2017) recently suggested that changes in solar spectral irradiance (most importantly in the UV) and energetic particle precipitation (EPP) are two major players for the climate from the aspect of solar variability. The EPP consists of particles originating from the sun, Earth’s magnetospheric field, and from beyond the solar system. Changes in EPP can lead to the production of NOx through ionization and dissociation and the subsequent intricate chemical processes in the middle atmosphere (Jackman et al., 2006; Rozanov et al., 2012). Then, the polar winter descent of EPP-generated NOx into the stratosphere affects the ozone abundance (Baumgaertner et al., 2011; Matthes et al., 2017). Sinnhuber et al. (2018) suggested that stratospheric EPP signals are more likely a result of dynamic processes than caused by downward descent of EPP-induced chemistry from the mesosphere and lower thermosphere. Rozanov et al. (2005) provided evidence that EPP-generated variations in stratospheric and tropospheric temperatures are comparable to those caused by solar UV forcing. The possible mechanisms linking the influence of EPP on the stratosphere to the surface are further suggested to be similar to those related to solar UV influence, including the accelerated stratospheric polar night jet and top-down coupling (Sepp?l? and Clilverd, 2014). Cnossen et al. (2016) also discussed the possible mechanism of downward connection from the thermosphere to the troposphere. However, a modeling study found that the EPP effect on tropospheric temperature is small and insignificant (Meraner and Schmidt, 2018).
As reviewed by Matthes et al. (2017), the influence of different EPP components on the surface climate is an emerging research topic, the challenge of which, however, is to quantify its climate impact. The solar wind, which cannot penetrate into the lower atmosphere directly but can violently interact with Earth’s magnetosphere, serves as a source of EPP (Mironova et al., 2015). The solar wind resembles the high-speed stream generated from the solar coronal regionand mainly encompass electrons, protons, and α-particles. Lu et al. (2008) found a significant relationship between the solar wind and the Northern Annular Mode. The energy input from the solar wind into Earth’s magnetosphere can cause space weather phenomena, such as magnetic storms, aurora, and so on (Akasofu, 1981). Magnetospheric particle precipitation associated with high-speed solar wind streams is more frequent during the decaying phase and near the minimum of the solar cycle (Sinnhuber et al., 2012).
Early studies suggest that the solar wind energy input is largely determined by the solar wind velocity (Crooker et al., 1977) or by the interplanetary magnetic field (IMF) (Dungey, 1961). Different opinions have been proposed since, but it remains a big challenge to quantitatively estimate the solar wind energy input into Earth’s magnetosphere (Ein) (Akasofu, 1981; Newell et al., 2008). Recently, Wang et al. (2014) quantified Ein as a function of interplanetary and solar wind conditions, using three-dimensional magnetohydrodynamics. The new Ein index shows both quasi-decadal variability and interannual variability (He et al., 2018, 2019b), which provides a direct motivation for us to explore the potential interannual connections between the solar activity?controlled magnetospheric energetic particle forcing and the surface climate. Furthermore, He et al. (2018) revealed for the first time that Ein has a tangible effect on the interannual variability of the subsequent winter ENSO. With this in mind, is it plausible that the year-to-year variability of Eurasian winter weather and climate is linked to the solar wind energy penetrating Earth’s magnetosphere. Thus, in this study, we attempt to understand the implications of the preceding Ein for winter extreme cold temperatures over Eurasia on the interannual time scale.
where nSW (units: cm?3) is the solar wind number density; VSW (units: km s?1) is the solar wind velocity; BT (units: nT) is the transverse magnetic field magnitude; θ is the IMF clock angle; and the solar wind data can be taken from the NASA OMNI project since 1963 (
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-1_mini.jpg)
Daily surface air temperature (SAT) and atmospheric circulation data are derived from the NCEP?NCAR Reanalysis 1 with a horizontal resolution of 2.5° × 2.5° (Kalnay et al., 1996). The daily NAO indices are obtained from the Climate Prediction Center of NOAA (
where N is the number of samples; r1 and r2 are the autocorrelation coefficients of the two variables at one time interval (Bretherton et al., 1999).
Cold events are defined when the daily SAT is below the climatological mean SAT by 1.5 or more standard deviation (Thompson and Wallace, 2001). This study focuses on the total number of cold events in winter, as well as the number of extreme circulation events. Changes in extreme circulation events may not be simply identical to changes in the mean climate and are of greater concern to the public (Screen, 2017a). The extreme sea level pressure (SLP) events, 850-hPa wind (UV850; cyclone) events, and 300-hPa zonal wind (U300) events, are defined as the daily SLP/UV850/U300 above its 90th percentile (Pfahl, 2014). The negative (positive) NAO events are defined as intervals in which the normalized daily NAO index is lesser (greater) than or equal to ?1.0 (+1.0) standard deviation for at least three consecutive days (Yao and Luo, 2015). Based on this definition, the frequency of negative (positive) NAO events in a winter is obtained by summing the duration of all negative (positive) NAO events in that winter. The 3?8-day bandpass-filtered 300-hPa transient eddies (v’2) are calculated according to Murakami (1979).
We use the Eliassen?Palm (EP) flux (F) and its divergence (DF), which are calculated according to Edmon et al. (1980), to diagnose the wave?mean flow interaction (Andrews et al., 1987):
where ρ is the air density;
![](http://www.iapjournals.ac.cn/aas/en/article/doi/10.1007/190153_M1.png)
![](http://www.iapjournals.ac.cn/aas/en/article/doi/10.1007/190153_M.png)
3.1. Interannual relationship between Ein and subsequent winter Eurasian cold extremes
The frequency of Eurasian cold events in the subsequent winter regressed onto the annual mean Ein index is firstly presented in Fig. 1a. Clearly, there is a significantly decreased frequency of cold extremes over northern Asia, Russia and northern Scandinavia, primarily north of 50°N, and an increased frequency over northern Africa and near the Black and Caspian seas (Fig. 1a). The relationship between Ein activity and winter cold events at lag 0 or +2 years is weaker and less significant (Fig. 2). To further investigate the linkage between Ein and the interannual variability of extreme cold weather in the subsequent winter, the dominant interannual modes of cold events over Eurasia (20°?80°N, 0°?150°E) during 1964?2018 are extracted by performing empirical orthogonal function (EOF) analysis. The first mode of EOF (EOF1) explains 23.8% of the total variance, which is 8.6% larger than the second mode. The corresponding principal component of the EOF1 mode (PC1) shows remarkable interdecadal and interannual variability (figure not shown). EOF1 exhibits consistent negative values across the northern Eurasian continent and positive values in northern Africa (Fig. 1b). This pattern of cold events is quite similar to that following higher Ein in the preceding year, with a high spatial correlation coefficient of 0.79. The regions that cover (45°?70°N, 40°?140°E) and (20°?45°N, 0°?50°E) (red and blue rectangles in Figs. 1a and b) are hereafter referred to as northern Eurasia (NE) and northern Africa (NA), respectively. The linear correlation between the annual mean Ein index and the frequency of cold events over NE and NA are ?0.38 and 0.30, respectively, significant at the 99% confidence level (Fig. 1c). Moreover, the probability for northern Eurasian winter to experience more than 30 extreme cold days following a higher (lower) Ein year is about 26.7% (48.2%), based on the probability density function (PDF) of the occurrence of cold extremes; the PDF for cold days in NA at lag +1 year shows more occurrences of more than 30 extreme cold days (from 33.0% to 49.4%), in response to higher Ein relative to lower Ein (figure not shown). It is thus hypothesized that the magnetospheric particle precipitation from solar wind is associated with the interannual variations of Eurasian cold events in the subsequent winter.![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-2_mini.jpg)
2
3.2. Extreme circulation events related to Ein variations
Extreme cold temperatures over Eurasia are affected by anomalous circulations that are related to various forcing (Liu et al., 2012; Cohen et al., 2014). Figure 3a shows a decrease in the frequency of SLP extremes associated with PC1 in the Arctic region and a pronounced increase in the frequency at midlatitudes of the Europe?Atlantic sector (Fig. 3a). More frequent days with the lower-level anomalous anticyclone appear at midlatitudes (Fig. 3b). In the upper troposphere, more (less) frequent extreme strong westerly wind is shown to the north (south) of the climatological westerly jet (Fig. 3c), consistent with the significant positive (negative) U300 anomalies in the north of 50°N (south of 40°N) (figure not shown). This suggests the poleward shift of the North Atlantic westerly jet. The more (less) frequent strong westerly wind near the Norwegian (Mediterranean) Sea would drive milder (colder) conditions in NE (NA). Clearly, the anomalous circulation extremes related to the higher-than-normal Ein in the preceding year (Figs. 3d-f) are qualitatively in good agreement with those related to PC1 (Figs. 3a-c). This provides further support for the interannual connections between the solar wind?magnetosphere energy flux, the subsequent winter circulation extremes in the Europe?Atlantic sector, and extreme temperatures over Eurasia.![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-3_mini.jpg)
The North Atlantic synoptic-scale eddy variations play an important role in the highly variable midlatitude weather patterns (Cohen et al., 2014). Figure 4 presents the anomalous storm-track activity response to PC1 and the preceding annual mean Ein, separately. As expected, significantly enhanced storm-track activity emerges primarily in the high-latitude North Atlantic (Fig. 4). Previous studies have suggested that positive storm-track activity anomalies are intimately linked to westerly wind anomalies in situ, cyclonic eddy forcing to the north, and anticyclonic eddy forcing to the south (Lau, 1988; Gong et al., 2011; He et al., 2019a). The northward shift of the synoptic-scale eddies favors the northward shift of the North Atlantic westerly jet and more frequent strong westerly wind at high latitudes (Figs. 3c and f). The increases in the strength of the synoptic-scale eddy forcing at high latitudes are also tied to the cyclone extremes to the north and anticyclone extremes to the south (Figs. 3b and e). Considering the wider literature, changes in the North Atlantic storm tracks are consistent with the shift of the NAO phase (Bader et al., 2011; Cohen et al., 2014). A poleward shift of the storm track occurs when the NAO is in its positive phase (NAO|+) and an equatorward shift is observed in negative NAO (NAO|?) winters. The influence of geomagnetic activity on the NAO has been examined in previous studies (Baumgaertner et al., 2011; Li et al., 2011). A dynamical change of the positive shift of the NAO is thus expected in association with the higher solar wind?magnetosphere energy flux input.
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-4_mini.jpg)
Next, we turn our attention to the frequency of strong NAO events following higher and lower Ein years. As shown in Fig. 5a, the composite occurrence of 1-year-lagged NAO|+ events is higher-than-normal following higher Ein, and lower-than-normal following lower Ein. Similarly, the NAO|? events have a higher frequency following lower Ein compared to higher Ein (Fig. 5a). These differences support the notion that the subsequent winter NAO|+ (NAO|?) events are intensified (weakened) following higher Ein activity, favoring less frequent cold events in NE and more frequent cold events in NA, and vice versa (Fig. 5b).
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-5_mini.jpg)
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-6_mini.jpg)
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-7_mini.jpg)
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2020/6/PIC/190153-8_mini.jpg)
The solar wind serves as a source of magnetospheric energetic particle forcing, in the main form of electrons (Mironova et al., 2015). Recently, Wang et al. (2014) quantified the solar wind energy input into Earth’s magnetosphere (Ein) using three-dimensional magnetohydrodynamics. It is worth noting that this new Ein index exhibits both interdecadal and interannual variability (He et al., 2018, 2019b). The interannual relationship between the solar wind energy penetrating Earth’s magnetosphere and Eurasian cold extremes in the subsequent winter is thus explored in this study. The results show that the frequency of cold winter events following years with increased solar wind?magnetosphere energy flux shows a significant decrease over NE and increase in NA, and vice versa. This is broadly similar to the dominant EOF mode of interannual variability in Eurasian cold extremes, with a high spatial correlation coefficient of 0.79. The high variance of winter cold events over NE and NA explained by the preceding annual mean Ein can reach up to 0.2 (figure not shown). Moreover, the frequent positive NAO (NAO|+) events and poleward shift of the North Atlantic storm track activities following the preceding higher-than-normal Ein activity, favors the anomalous circulation extremes and less frequent extreme cold temperatures over northern Eurasia. It is further hypothesized that the wave?mean flow interaction in the stratosphere and troposphere is favorable for the connection of Ein signals to tropospheric circulation and climate in the following winter. In general, our study provides further insight into the potential interannual linkage between EPP variations and the surface extreme cold weather.
Acknowledgements. This work was supported by the National Key R&D Program of China (Grant No. 2016YFA0600703), the National Natural Science Foundation of China (Grant Nos. 41875118, 41605059 and 41505073), the Postgraduate Research & Practice Innovation Program of Jiangsu Province (Grant No. KYCX18_0997), the funding of the Jiangsu Innovation & Entrepreneurship Team, and the Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions. This work was partially supported by the International Space Science Institute (ISSI) in Beijing, through the working team “Dynamical signatures of energetic particle precipitation in atmospheric re-analyses” (ID:32, 2019).