Abstract

The frequency distribution of winter extreme cold events (ECEs) in North China and the influences of mid-latitude sea surface temperature anomalies (SSTAs) in the Northern Hemisphere are studied. The results show that (1) the frequency of single station ECEs (SSECEs) in winter increases from southeast to northwest, with a decrease before 2008 and then a significant increase. This trend abrupt change occurs in late winter. (2) When the SST in the North Pacific shows an “El-Niño-like” anomaly in winter, it triggers the negative Arctic Oscillation (−AO), positive Pacific North America (+PNA), and positive Eurasia Pacific (+EUP) atmospheric teleconnection patterns in the mid-lower troposphere. As a result, the ridge to south of Lake Baikal becomes stronger. Meanwhile, SST in the North Atlantic shows a “reversed C” negative anomaly with North Atlantic Oscillation (+NAO), (+PNA)-like and (+EUP)-like patterns, and the ridge to southwest of Lake Baikal becomes stronger. Furthermore, both cause the Siberian High to become weaker in the north and stronger in the south. With the weaker East Asia subtropical jet and stronger East Asia winter monsoon, these factors lead to a significant increase of SSECE frequency in North China. (3) When the SSTA shows an “El Niño-like” developing pattern from summer to autumn in the North Pacific, the winter SSECE frequency will be higher. (4) The purported mechanism between the mid-latitude SSTA and the winter SSECE frequency in North China is the following: the SSTA in the North Pacific in summer and autumn excites atmospheric teleconnection wave trains, and the Atlantic stores these anomaly signals. In winter, the interaction between the SSTAs in the North Pacific and the North Atlantic enhances the Eurasian teleconnection wave train. With the upstream fluctuation energy dispersing downstream, the wave train centers move eastward with the season, resulting in an increase in the frequency of the SSECEs.

1. Introduction

According to the sixth special report of the Intergovernmental Panel on Climate Change (IPCC), the increase in global mean surface temperature (GMST) reached 0.87°C in 2006–2015 [1]. Under the background of global warming, extreme cold events (ECEs), such as freezing rain and snow [2], extreme low temperatures [35], and persistent low temperatures [6], are more frequent and threaten people’s lives and property. It is important to study extreme weather events in North China as a major area of agricultural economic development and an important transportation hub in China [7, 8].

The frequency and intensity of winter ECEs in Northern China are vulnerable to the winter monsoon circulations systems [912]. Ding [9] studied the formation and spread of the Siberian High and its relationship with cold waves in East Asia. In winter, the intensification and northward extension of the Urals blocking high and the positive longitude anomaly of the North Atlantic storm track can strengthen the Siberian High and deepen the East Asian trough, which cause persistent low-temperature events in East Asia [5, 10, 11]. However, the frequency of cold waves does not increase monotonically with the strengthening of the Siberian High [12]. ECEs occur when Siberian High and the Aleutian Low are both strong [13] and are related to the location of these two systems [14]. Furthermore, Tao [15] divided the paths of cold waves invading China into a northwest path, a superpole cold wave path, and a southward invasion path of cold air from the east of Lake Baikal. Different cold air paths correspond to different configurations of atmospheric systems. Besides, the frequency before 2008 has obviously decreased [16, 17], while it is closely related to the weakening of the East Asian winter monsoon and the warming of the low-level cold air stack, especially in the Inner Mongolia, North China, and Jianghuai areas [4, 12, 18]. At the same time, the interaction between the East Asian subtropical jet and temperate jet is an important factor affecting the winter ECEs in China [19].

ECEs are also related to atmospheric teleconnection, e.g., the Arctic Oscillation (AO) and the North Atlantic Oscillation (NAO) [2026]. The AO affects the Siberian High by modulating quasi-stationary planetary waves. When the AO is in a positive phase, the East Asian trough, Siberian High, and East Asian winter monsoon are weaker, and the temperature in China is higher [2024]. In winter, when the NAO is in a positive phase, the temperature in China is higher [25, 26]. In addition, external forcing anomalies such as SSTAs also have important effects on the winter temperature in China. For example, the El Niño-Southern Oscillation (ENSO) [27, 28], SSTAs in the tropical Indian Ocean [29], Kuroshio Sea areas [30], North Pacific [31], and North Atlantic [31] all have strong correlations with the winter temperature in China.

In summary, there have been many studies on atmospheric circulation and external forcing anomalies affecting the winter temperature in China or North China. In recent years, cold waves in the mid-latitudes of Eurasia have been frequent and persistent [32]. North China is one of the regions most affected by cold waves, and it is necessary to further research the new characteristics and possible causes of ECEs, such as cold waves in North China. In this paper, the temporal and spatial characteristics of the winter ECE frequency in North China since 1980 are analyzed, and the influence of the mid-latitude SSTA and its atmospheric teleconnection patterns in the Northern Hemisphere on the interannual anomalies of ECE is emphasized, to identify interannual prediction signals and provide a reference for short-term climate prediction in North China.

2. Data and Methods

2.1. Data and Preprocessing

The data include (1) daily minimum temperatures at 2,480 stations provided by the National Climate Center in China; (2) the NCEP/NCAR monthly sea-level pressure field, the 850 hPa temperature and wind fields, the 200 hPa zonal wind field, and the 500 hPa geopotential height field, with a resolution of 2.5 × 2.5 [33]; and (3) the NOAA monthly SST, with a resolution of 1 × 1 [34]. The data are for the 1980–2016 period.

North China is defined as the study region (35°–42.5°N, 110°–120°E). To ensure the continuity and integrity of the data, the stations from which the missing data of the daily minimum temperature does not exceed 1% are retained. Excluding the high mountain stations which are higher than 2500 m, 421 stations are finally selected. The distribution is shown in Figure 1.

2.2. Methods

Based on the regulations of “Cold Air Level” (GB/T 20484-2017) and “Cold Wave Level” (GB/T 21987-2017) published by the National Meteorological Center of China Meteorological Administration in 2017, a single station extreme cold event (i.e., a single station cold wave) is defined by a daily minimum air temperature lower than and equal to 4°C and either one of three conditions: a decrease by no less than 8°C within 24 hours, or a decrease by no less than 10°C within 48 hours, or a decrease by no less than 12°C within 72 hours. For convenience, an extremely cold event can be abbreviated as an ECE, with a single station ECE as a SSECE.

The frequency of winter ECEs in a given year is the total number of ECEs occurring from December to February of the following year.

In this paper, statistical analysis methods including correlation, composite analysis, and linear regression are mainly used. Student’s t-test method is used to determine the statistical significance of the correlation coefficients and linear trend.

3. The Spatial and Temporal Characteristics of the ECE Frequency in Winter in North China

3.1. Single Station ECEs
3.1.1. Spatial Distribution

The multi-year average of the SSECE frequency during the period of 1980 to 2016 is calculated (Figure 2). The frequency of SSECEs increases from southeast to northwest in winter. More than 4 SSECEs occurred in western Inner Mongolia, northern Shanxi, northeastern Shaanxi, northwestern Hebei, and northern Liaoning, with the maximum frequency 8 in western Inner Mongolia (Figure 2(a)). The frequency distribution of SSECEs in December (Figure 2(b)), January (Figure 2(c)), and February (Figure 2(d)) is similar to that in winter.

3.1.2. Temporal Variations

The regional averages of the SSECE frequency in North China are calculated for the winters of 1980–2016. According to the statistical results, the 1.92 SSECEs occurred per winter on average, and the maximum (minimum) value is 3.4 (0.8) in 2008 (2007), with larger interannual differences. The average frequency is 0.64, 0.66, and 0.62 per year in December, January, and February, with the maximum (minimum) value being 1.94 (0.16), 1.31 (0.09), and 1.10 (0.17), respectively.

The frequency of SSECEs in winter decreases from 1980 to 2007 (Phase I) and increases significantly from 2009 to 2016 (Phase II) (Figure 3(a)). The frequency of SSECEs in January (Figure 3(c)) and February (Figure 3(d)) shows an abrupt trend similar to that in winter. However, December (Figure 3(b)) is an exception in that the frequency still decreases in Phase II. Therefore, the abrupt change in the SSECE frequency in winter is mainly caused by the SSECEs in later winter (January and February). Why is the frequency trend in December different from that in winter? It is worth further study.

Figure 4 further shows the regional distribution of the linear trend of the SSECE frequency. Except for the junction of southwest Shanxi and Shaanxi, the frequency of the SSECEs decreases significantly in most areas of North China in Phase I in winter (Figure 4(a)). In Phase II (Figure 4(e)), except for a significant increase at the junction of the southwest of Inner Mongolia, and central and western Shanxi and Shaanxi, most areas show a decrease. The trends of the SSECE frequency in December (Figure 4(b)), January (Figure 4(c)), and February (Figure 4(d)) in Phase I are similar to those in winter, but the significant areas of the trend are slightly different. In Phase II, the frequency decreases over the entire region in December (Figure 4(f)) but increases in most areas in January (Figure 4(g)) and February (Figure 4(h)). The high value areas are located at the junction of central and western Shanxi, Shaanxi, and Inner Mongolia, which indicates once again that the increase in the SSECE frequency in winter in Phase II is caused by those in later winter (January and February).

4. The Influence of Mid-Latitude SSTA on the Winter SSECE Frequency in North China

Previous research suggested that the mid-latitude SSTA in the Northern Hemisphere has a significant influence on the winter temperature in China [30, 31]. How are the anomalies of the mid-latitude SST in the Northern Hemisphere related to the frequency of the SSECEs in North China in winter? Moreover, what kind of anomalous patterns of atmospheric circulation is caused by the SSTA? Through the study of these questions, this paper attempts to reveal the possible causes for the frequent ECEs in North China. It also looks for its possible interannual prediction signals, to provide a reference for short-term climate prediction in North China. The first question will be studied in Section 4.1, and the second will be discussed in Sections 4.2 and 4.3.

4.1. The Relationship between SSTA and SSECE Frequency

Figure 5 shows the lead-lag and contemporaneous correlation coefficients between the regional average frequency of the winter SSECEs (Figure 3(a)) and SSTA over the northern hemisphere in spring, summer, autumn, and winter, respectively. According to Figure 5, the significant correlation area in the North Pacific is too small in spring, so we select the key areas in summer, autumn, and winter. For the North Atlantic, the correlation coefficient between the regional average SSTA in the significant correlation areas in spring (Figure 5(a)) and the regional average frequency of the winter SSECEs in North China is not statistically significant, and the significant correlation areas are too small in summer and autumn, so we select the key areas in winter. For convenience, the significant correlation areas of each season in Figure 5 are indicated by the boxes labeled A, B, and C. A is the region (8–15°N, 145–175°E), B is (25–45°N, 180°–200°W), and C consists of three subareas of (55–60°N, 15–50°W), (30–50°N, 11–20°W), and (20–25°N, 20–50°W).

The correlation coefficient between the regional average SSTA over the above significant correlation area of the North Pacific and the Niño3 SST index in winter is −0.52, which represents statistical significance at the 99% confidence level using Student’s t-test. The significant correlation areas are positive in the east and negative in the west of the North Pacific in Figure 5(d), which is similar to an El Niño pattern over the tropical Pacific. Because the two patterns appear in different parts of the ocean, the SSTA pattern over the North Pacific in Figure 6(d) is called the “El Niño-like” pattern for convenience, while the SSTA in summer and autumn is collectively called the “developing El Niño-like” pattern.

To further understand the correlation relationships between SSTA and the frequency of the SSECEs, the SST indices of North Pacific and North Atlantic for each season are defined as INP and IATL, respectively. That is, the regional average series of the SSTA over the significant correlation area from 1980 to 2016 is calculated and is then standardized as SST index series. If the significant correlation area contains two or more subareas, the subareas will be combined and considered as one. The correlation coefficients between these indices for each season and the regional average frequency of the winter SSECEs are calculated, and those with statistical significance at the 95% or 99% confidence level are listed in Table 1. When the SST in these key sea areas is anomalously cold (warm) in the corresponding seasons, the SSECEs in North China tend to be more (less) frequent in winter. The SSTA in the key areas of the North Pacific in summer and autumn has predictive significance for the winter SSECEs in North China.

The SSTAs over the North Pacific and North Atlantic key areas in winter are both negatively correlated with the winter SSECE frequency. Furthermore, the partial correlations among them are calculated, and whether the SSTAs of the two key areas influence each other is analyzed. The partial correlation coefficient between INP, IATL and the regional average frequency of the winter SSECEs in North China are −0.29 and −0.41, respectively. The corresponding the correlation coefficient in Table 1 is −0.43 and −0.48, respectively. Therefore, the partial correlation coefficient for the North Atlantic SST changes very little, and their absolute value obviously decreases for the North Pacific SST but is still statistically significant at the 95% confidence level. Therefore, it can be concluded that the impacts of the SST over the two key areas on the winter SSECE frequency in North China are relatively independent but have a certain combined effect.

In order to further understand the regional correlation distribution between the SST indices of the North Pacific and North Atlantic and the winter SSECE frequency in North China, Figure 6 shows the coefficient distribution among INP in summer, autumn, winter, IATL in winter and the winter SSECE frequency in North China. According to Figure 6, the significant negative correlation areas between INP and the frequency are mainly located in the south-central part of North China, namely, the junction areas of Hebei, Shanxi, Henan, and Shandong (Figures 6(a)6(c)). The IATL has a significant negative correlation with the frequency not only in the regions mentioned above but also in the north-central part of North China (Figure 6(d)).

4.2. The Response of Atmospheric Circulation to the North Pacific SSTA

To understand how the SSTA in the North Pacific and North Atlantic affects the atmospheric circulation anomalies and thus the frequency of the SSECEs in North China in winter, the following section will analyze the linear regression coefficients of INP in summer, fall, and winter (Section 4.2) and IATL in winter (Section 4.3) on the atmospheric circulation factors in winter in the Northern Hemisphere, respectively.

4.2.1. Summer

Figure 7 shows the linear regression coefficients between INP in summer and the atmospheric circulation variables in winter. To facilitate understanding of the circulation anomalies when the SST in the key areas is anomalously cold and the SSECE frequency in winter in North China is high, all regression coefficients in Figure 7 are multiplied by −1, and similar processing is done in Figures 810.

As seen in Figure 7, when the SST in the tropical northwest Pacific in summer is anomalously cold (Figure 5(b)), the sea level pressure (SLP) field in winter (Figure 7(a)) is likely to show an (-AO) anomalous pattern. Eurasia has an anomalous pressure gradient force, with its direction from polar region to the mid-high latitudes, which strengthens the meridional circulation and is conducive to the southward movement of cold air in the polar region. The PNA wave train shows a pattern “−” “+” “−” “+” starting in the North Pacific to North America, while in the North Atlantic to Eurasia, the wave train shows a pattern “+” “−” “+” “−”. The two wave trains converge in the subtropical northwest Atlantic and cause a north-negative and south-positive anomalous distribution from the Eurasian mid-high latitudes (50–80°N) to the southeast side of Lake Balkhash, which is similar to the empirical orthogonal function (EOF) second mode of the interannual anomaly of the Siberian High developed by Zhu et al. [35]. Its interaction with the Aleutian Low located southwest of the normal position produces a strong East Asian winter monsoon (Figure 7(e)).

The circulation patterns of the 500 hPa geopotential height field (H500, Figure 7(b)) and SLP field (Figure 7(a)) basically show a quasi-barotropic structure in the mid-high latitudes. The Urals ridge is stronger, and the East Asian trough is toward the west and south. A cold air stack in the 850 hPa temperature anomalous field (T850) lies near Greenland, which provides conditions for the southward eruption of cold air (Figure 7(c)). The East Asia temperate jets and the subtropical jets in the 200 hPa zonal wind field (U200) over the Middle East to the Qinghai-Tibet Plateau are weaker (Figure 7(d)) and conducive to moving southward for the northwest air flow in front of the Urals ridge and behind the East Asian trough, which leads the 850 hPa low-level cold air (UV850) to North China through Greenland, Novaya Zemlya, and the west side of Lake Baikal (Figure 7(e)) and causes frequent SSECEs in North China.

4.2.2. Fall

When the negative SSTA expands from the tropical northwest Pacific to the central North Pacific in fall (Figure 5(c)), the (−AO), (+PNA)-like, and (+EUP)-like quasi-barotropic teleconnection patterns or wave trains in the middle and lower troposphere (Figures 8(a) and 8(b)) are more prominent. The north-south seesaw anomaly of the Siberian High is clearer (Figure 8(a)), the ridge previously located in the Urals (Figure 7(b)) moves to between Lake Balkhash and Lake Baikal (Figure 8(b)), and the cold air stack between Greenland and Severnaya Zemlya is further strengthened (Figure 8(c)), as is the stronger temperate jet and weaker subtropical jet (Figure 8(d)). All of these are favorable for the cold air to arrive in North China along Western Europe, the East European Plain, West Siberia, and the west of Lake Baikal (Figure 8(e)).

4.2.3. Winter

In winter, along with the formation of an “El Niño-like” pattern in the North Pacific (Figure 5(d)), the atmospheric circulation anomalies (Figure 9) continue developing and are very similar to those in autumn (Figure 8). However, the anomalies of the system centers are stronger, e.g., the (−AO), (+PNA), and (+EUP) teleconnection patterns or wave trains (Figure 9(a) and 9(b)), leading to stronger meridional circulation. Meanwhile, there are stronger temperate and weaker East Asian subtropical jets (Figure 9(d)), stronger low-level cold air stacks in the vicinity of Novaya Zemlya (Figure 9(c)), and stronger East Asian winter monsoons (Figure 9(e)). The above circulation anomalies are all favorable for cold air to invade North China from Novaya Zemlya, the East European Plain, the Caspian Sea, West Siberia, and Lake Baikal (Figure 9(e)), resulting in frequent ECEs in North China.

4.3. The Response of Atmospheric Circulation to the North Atlantic SSTA

Except for (+NAO) rather than (−AO), the circulation anomalies in the Pacific, North America, and Eurasia that are excited by the “reversed C” negative SSTA in the North Atlantic in winter (Figure 10) are very similar to those excited by the “El Niño-like” SSTA in the North Pacific (Figure 9). The atmospheric circulation displays “(+PNA)-like” and “(+EUP)-like” quasi-barotropic teleconnection wave trains in the mid-lower troposphere (Figures 10(a) and 10(b)). The Siberian High presents the north-south seesaw anomaly (Figure 10(a)), and a ridge lies between Lake Balkhash and Lake Baikal (Figure 10(b)). The East Asian temperate is stronger and subtropical jet is weaker (Figure 10(d)). There is a cold air stack between Victoria Island and the Greenland Sea (Figure 10(c)). The configuration of upper and lower circulation systems is conducive to the arrival of cold air to in North China from the Greenland Sea, East European Plain, West Siberia, and west of Lake Baikal (Figure 10(e)), resulting in frequent ECEs.

5. Conclusion and Discussion

This paper analyzes the characteristics of the SSECE frequency in North China in winter. Furthermore, a possible way of SSTA in the northern mid-latitudes affecting the SSECE frequency is described. The major conclusions are given and the possible mechanism are discussed as follows.(1)The winter SSECE frequency in North China increases gradually from southeast to northwest, decreasing before 2008 and then significantly increasing, which is caused by the SSECEs in late winter.(2)The SSTA in the North Pacific and North Atlantic affects the winter SSECE frequency in North China by exciting atmospheric teleconnection patterns or wave trains. When the SST in the North Pacific shows an “El Niño-like” anomaly in winter, it triggers the −AO, +PNA, and +EUP teleconnection patterns or wave trains in the mid-lower troposphere. Meanwhile, when the SST in the North Atlantic shows a “reversed C” negative anomaly in winter, it triggers + NAO, +PNA-like, and +EUP-like patterns. Together, the two deepen the Urals trough and strengthen the ridge to the south or southwest of Lake Baikal. Consequently, the Siberian High is weaker in the north and stronger in the south. Along with the weaker subtropical jet in the upper troposphere and the stronger East Asia winter monsoon, these factors cause cold air to invade North China, resulting in frequent SSECEs in North China.Zhu et al. [35] indicated that when the Siberian High is abnormally weak in the north and strong in the south, the northeast part of East Asia is colder and vice versa, which is consistent with our conclusion. In this paper, it is further pointed out that the interannual anomalies of the Siberian high seesaw are related not only to AO, PNA, and EUP teleconnection waves induced by ENSO-type SSTA in the North Pacific Ocean, but also to NAO, PNA-type, and EUP teleconnection waves induced by inverted C-type SSTA in the North Atlantic Ocean.(3)When a negative SSTA develops from the tropical northwest Pacific to the central North Pacific from summer to fall and the SSTA shows an “El Niño-like” development trend, the SSECE frequency will be higher in winter.(4)A possible mechanism for the impacts of the northern mid-latitude SSTA on the winter SSECE frequency in North China is discussed and speculated below.

The previous research showed that the interdecadal variation of the mid-latitude air-sea system in the northern hemisphere may have a coupling mechanism of “two oceans and one atmosphere” [36, 37]. Liu and Song [31] pointed out that the positive (negative) SSTA in the eastern (Midwest) part of the North Pacific triggered an anomalous circulation similar to Figure 8(b) in the North Pacific-North America-North Atlantic in autumn, resulting in the circulation anomaly over the North Atlantic Ocean, and then the atmosphere forced the ocean, which caused the North Atlantic Ocean to store this anomalous signal in autumn. Then, the Atlantic-Eurasian teleconnection pattern anomaly affects the circulation in the East Asian region through the energy dispersion of upstream and downstream fluctuation in winter [38, 39], thereby affecting the temperature in winter in China [40].

Sections 4.2 and 4.3 indicate that the SSTAs in the North Pacific and North Atlantic have combined effects on the atmospheric circulation anomalies and jointly affect the SSECE frequency in North China. Combined with the above existing research results, the possible mechanism of northern mid-latitude SSTA affecting the winter SSECE frequency in North China is speculated: the negative SST anomalies in the Midwest part of North Pacific excite the teleconnection wave train in the North Pacific-North America-North Atlantic-Eurasia in summer and autumn, and then the atmosphere forces the ocean, resulting in the North Atlantic storing these atmospheric signals. In winter, the persistent SST anomalies in the North Pacific and North Atlantic interact with each other, and their superposition enhances the Eurasian teleconnection wave train. The upstream wave energy disperses downstream, and the wave train centers move eastward from summer to winter, leading to an increase of the SSECEs frequency in North China. Of course, further time scale separation [41], dynamic diagnosis, and simulation verification are needed to support this speculation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This research was supported jointly by the National Key R&D Program of China (2018YFC1505602) and NSFC (41790471). The authors acknowledge the United States National Centers for Environmental Prediction/National Center for Atmospheric Research for providing the reanalysis data.