• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Seasonal Cumulative Effect of Ural Blocking Episodes on the Frequent Cold events in China during the Early Winter of 2020/21※

    2022-04-06 08:49:08YaoYAOWenqiZHANGDehaiLUOLinhaoZHONGandLinPEI
    Advances in Atmospheric Sciences 2022年4期

    Yao YAO, Wenqi ZHANG, Dehai LUO, Linhao ZHONG, and Lin PEI

    1Key Laboratory of Regional Climate-Environment for Temperate East Asia, Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China

    2University of Chinese Academy of Sciences, Beijing 100029, China

    3Institute of Urban Meteorology, China Meteorological Administration, Beijing 100029, China

    ABSTRACT Starting in mid-November, China was hit by several cold events during the early winter of 2020/21. The lowest temperature observed at Beijing station on 7 January reached -19.6°C. In this paper, we show that the outbreak of the record-breaking extreme cold event can be attributed to a huge merging Ural blocking (UB) ridge over the Eurasian region.The sea-ice cover in the Kara and East Siberia Seas (KESS) in autumn was at its lowest value since 1979, which could have served as a precursor signal. Further analysis shows that several successive UB episodes occurred from 1 September 2020 to 10 January 2021. The persistent UB that occurred in late September/early October 2020 may have made an important contribution to the October historical minimum of sea ice in the KESS region. Our results also show that, after each UB episode in winter, significant upward propagation of wave activity occurred around 60°E, which resulted in weakening the stratospheric vortex. Meanwhile, each UB episode also caused a significant reduction in sea-ice extent in KESS and a significant weakening of the westerly jet in mid-high-latitude Eurasia. Results suggest that the Arctic vortex, which is supposed to enhance seasonally, became weaker and more unstable than the climatic mean under the seasonal cumulative effects of UB episodes, KESS warming, and long-lasting negative-phase North Atlantic Oscillation (NAO-). Those seasonal cumulative effects, combined with the impact of La Ni?a winter, led to the frequent occurrence of extreme cold events.

    Key words: extreme cold events, Ural blocking, Arctic sea ice, Arctic vortex, cumulative effect

    1. Introduction

    On 6 January 2021, China was impacted by the strongest cold air outbreak of that winter, accompanied by widespread strong cooling and windy weather. Several meteorological stations observed their lowest temperatures on record for the same period during 6-8 January. For example, the 2-m temperature recorded at Beijing station on 7 January reached -19.6°C, which is the second-lowest temperature since 1951 (http://www.cma.gov.cn/2011wmhd/2011wzbft/2011wftzb/202101/t20210107_569795.html). In addition, China also experienced several frequent cold events in the early winter of 2020/21.

    In recent years, along with continued global warming and regional climate anomalies, extreme weather-especially extreme cold weather in winter-has been frequently observed over continental Eurasia (Ding et al., 2008; Hui,2009; Wu et al., 2011, 2017; Cohen et al., 2014; Yao et al.,2016, 2017) and North America (Whan et al., 2016; Cohen et al., 2018a). Extreme cold surges are an important research topic because they have a significant impact on many aspects of people’s daily lives-for example, their clothing, food, housing, and economic activities. Many studies have examined the physical causes of such extreme weather processes (Ding et al., 2008; Hui, 2009; Wu et al., 2011,2017; Luo et al., 2014; Yao et al., 2017; Zheng et al., 2021).For the invasion of cold air from the Arctic to mid-low-latitude land areas, anomalous meridional circulation is necessary to provide a suitable background such as blocking and a negative-phase of the North Atlantic Oscillation (NAO-)(Herring et al., 2016; Whan et al., 2016; Yao et al., 2016; Martineau et al., 2017). The meridional circulation can guide the southward movement of cold air in the polar regions and is the conceptual equivalent to opening a gap through the strong westerly jet. In addition, a climatic context is necessary for such extreme anomalous meridional circulation to occur-for example, ENSO, Arctic amplification, etc.(Ding et al., 2008; Hui, 2009; Wu et al., 2011; Yao et al.,2017; Zheng et al., 2021). In other words, the occurrence of an extreme weather process is often the result of a combination of several factors (Yao et al., 2016, 2017; Li et al.,2020; Zheng et al., 2021).

    Many studies (Luo et al., 2016; Gong and Luo, 2017;Chen et al., 2018; Tyrlis et al., 2019) have examined the important role of Ural blocking (UB) in driving Arctic seaice loss and Eurasian cold extremes and their link with the stratospheric vortex in autumn and early winter 2016-17 (Tyrlis et al., 2019). Large-scale atmospheric circulations such as UB and the Siberian high may increase the frequency of cold waves over Eurasia by intensifying the amplitude of the westerly jet meandering under the forcing of global warming, as revealed by Ma and Zhu (2019). Lü et al. (2020) suggested that the accumulation of snow in Siberia can contribute to sudden stratospheric warming and its downward transmission, which in turn facilitates the formation and development of a negative phase of the NAO, ultimately causing the onset of extreme cold weather in Eurasia. Zhang et al.(2019) suggested that the dipole mode of Eurasian winterspring snow cover can induce large-scale teleconnection patterns by influencing sea-ice concentrations (SICs), thus further prolonging the cold conditions in Eurasia. The link between Arctic amplification and mid-latitude extreme cold weather in the context of global warming has become a hot topic of research (Cohen et al., 2014; Wu et al., 2015, 2017;Overland et al., 2016; Francis et al., 2017; Wu, 2017).Many studies have examined how anomalies in pre-existing sea ice can affect the state of atmospheric circulation (Gong and Luo, 2017; Wu, 2017; Chen et al., 2018) or precipitation in later periods, and on time scales that are likely to be cross-seasonal (Li and Wang, 2012; Wu et al., 2017; Han and Li, 2018; Li et al., 2018; Lü et al., 2019; Ding et al.,2021). Han and Li (2018) indicated that the sea-ice anomaly in the Labrador Sea in winter seems to be a better precursor for predicting the spring precipitation anomalies over southeastern North America and Western Europe. Several authors have published detailed reviews on Arctic sea ice and its possible link with Eurasian climate, from which we can see that the possible link between sea ice in the Arctic and atmospheric circulation is a consistent assertion, but that controversies also still exist (Gao et al., 2015; Cohen et al., 2018b). Indeed, due to the uncertainties and chaotic nature of the atmosphere and other systems, research in this area is still very challenging, and the resultant academic debate is ongoing (Wu et al., 2017; Dai and Song, 2020).

    In any case, the extreme cold surges affecting the Eurasian midlatitudes are very closely related to meridional circulation anomalies (such as blocking). Although most cold surges have certain commonalities, each individual case is unique given the nonlinear and chaotic nature of the atmosphere. This uniqueness is what makes related in-depth investigations a continually worthwhile pursuit. In fact, in many cases, extreme weather events are dominated by one kind of weather process but at the same time require the involvement of other weather systems. For instance, Yao et al.(2016) proposed a mechanism for the combined effect of the NAO and European blocking on extreme snowstorms in Europe. In addition, their latest study highlights the unique combined effect of the NAO and European blocking on European heat waves (Li et al., 2020). As mentioned in these studies (Luo et al., 2015; Yao et al., 2016; Li et al.,2020), the combined effect of the NAO and blocking circulation mainly manifests as the NAO changing the position and shape of the downstream blocking by altering the strength,orientation, and position of the westerly jet. In addition, the ways in which the processes of the two systems can combine to impose a collective effect can be both spatial and temporal in form, and then possibly involve atmospheric teleconnection processes. In the case of a temporally combined process, the time scale may extend from a contemporaneous one to an interseasonal, interannual, or even an interdecadal one. This study examines the specific physical processes involved in this cold event and the possible combined effect that existed within these processes.

    The remainder of this paper is structured as follows:The datasets and methods are described in section 2. In section 3 we describe the evolution of the atmospheric circulation during the January 2021 cold event and its detailed features from a daily perspective. The historical minimum value of the sea ice in the Kara and East Siberia Seas(KESS) in autumn (October) 2020 is investigated in section 4, which is considered as a precursor to this cold event.Also in section 4, we show how the UB was organized into five successive events from 1 September 2020 to 6 January 2021. Meanwhile, the upward flux of wave activity induced by the UB can weaken the stratospheric vortex. The anomalies of the West Pacific sea surface temperature (SST) associated with La Ni?a are presented in section 5. Finally, conclu-

    sions and some further discussion are provided in section 6.

    2. Data and methods

    2.1. Datasets

    The main datasets used in this study are the latest daily reanalysis (ERA5) dataset (Hersbach and Dee, 2016) from the European Centre for Medium-Range Weather Forecasts,covering the period from 1979 to 2021 and including geopotential height, temperature, horizontal wind components (u and v), surface (2-m) air temperature (SAT), surface pressure, and SIC, at a 2.5° × 2.5° spatial resolution, and the monthly SIC and SST data derived from the Met Office Hadley Centre at a resolution of 1° × 1°. To corroborate the robustness of the sea-ice variability, the monthly and daily(1 September 2020 to 1 January 2021) datasets of the National Snow & Ice Data Center (NSIDC) in the northern polar stereographic projection, on a grid of 25 km × 25 km,are also used in this paper.

    To measure the SAT variability over the region of China under the impact of the cold event, the time-mean(1-8 January 2021) temperature data of 2414 stations nationwide are used, as well as the mean temperature data of nationwide stations from 1 September 2020 to 1 January 2021. Furthermore, the monthly tropical SST series, including the Ni?o-1+2, Ni?o-3.4, and Ni?o-4 indexes, and the daily NAO index, from the Climate Prediction Center(https://psl.noaa.gov/data/climateindices), are used to examine the impact of multi-factor synergy on this cold-air outbreak over East Asia.

    2.2. Methods

    The occurrence of this cold event was closely related to the UB, and so it is necessary to examine the evolution of the UB during this period. To prevent omitted UB cases,two kinds of identification methods are used in this research, the geopotential height anomaly, and the one-dimensional blocking index (Tibaldi and Molteni, 1990), abbreviated as TM index. A UB case is defined when at least one of these two conditions is met. Specifically, the first method works based on the maximum value of the 500-hPa geopotential height anomaly over the region (40°-85°N, 30°-150°E)being larger than 250 gpm for at least five consecutive days.The TM index method, meanwhile, refers to the reverse of the 500-hPa geopotential height gradient, which is defined by

    where Z500 indicates the 500-hPa geopotential height,φN=80°N+Δ, φ0=60°N+Δ, φS=40°N+Δ, and Δ=-5°,0°or 5°. A UB case can be recognized if two criteria,GHGS>0 gpm (°lat)-1and GHGN<-10 gpm (°lat)-1,are both met for at least 10 consecutive longitudes ( λ) and five continuous days within the zonal range of 30°-150°E.Lag-0 is the strongest day of the UB, which denotes when the maximum value of the geopotential height anomaly is at its largest during the lifespan of the UB.

    3. Extreme cold weather event in January 2021

    3.1. Evolution of the atmospheric circulation

    Fig. 1. Spatial patterns of time-averaged (6-8 January 2021) SAT anomalies(units: K) relative to corresponding daily means for 1981-2010 of 2414 stations across China.

    On 6 January 2021, China was hit by a widespread cold wave accompanied by strong cooling and windy weather.Figure 1 shows the SAT anomaly distribution averaged from 6 to 8 January based on the station-observed data(2414 stations). It can be seen that nearly all of China experienced a sharp SAT decline during 6 to 8 January, except for a small part in southwestern China. The average SAT in northern China from 6 to 8 January was more than -10°C colder than climatology. Even over the most southern region of China, the SAT still exhibited a clear decline during this cold wave. To visualize the physical processes during this cold event, it is necessary to examine the evolution of the atmospheric circulation, especially from the daily perspective. Figure 2 shows the daily evolution of the geopotential height at 500 hPa and the SAT anomalies from 29 December 2020 to 9 January 2021. On 29 December, a huge anticyclone already existed over mid-high-latitude Eurasian land areas, which can be identified as a UB process. It can be seen that the anticyclonic center covered a large area from around 0° to 120°E, with the maximum center located at approximately 60°E. The cyclonic center located to the southeast of the anticyclone was much weaker in intensity compared to the anticyclonic center. However, the anticyclonic center had a southward spread around 30°E to 60°E and reached 30°N on 4 January and even to 30°N on 9 January(Fig. 2). This may have been due to the mid-low-latitude wave train that developed from 29 December to 9 January.Starting from 29 December, we can see that there was a latitudinal wave train starting to develop from the North Atlantic to the Mediterranean Sea, which then extended to near the Caspian Sea. The wave train had a structure, and the last positive center of the wave train was the part of the UB anticyclonic center that extended southward during this period. In other words, the wave train and the UB both developed in the same period, and some parts of the two merged and coupled together during their development.

    Therefore, there were quite a few differences between this UB (Fig. 2) and the classical UB structure described in previous observational studies (see Fig. 8b in Yao et al.(2017)). In the North Atlantic region, the latitudinal wave train mentioned above had an unstable structure during its development from 29 December to 9 January (Fig. 2). From 2 to 9 January, the positive and negative anomaly centers over the North Atlantic region gradually shifted from a latitudinal to a longitudinal orientation and developed into a northwest-southeast-tilted dipole mode (Fig. 2). This dipole mode had a classic NAO- structure as shown in many previous studies (e.g., Benedict et al., 2004; Yao and Luo, 2014), only at a lower overall latitude than the mean state. An NAO in a more shifted mean state will have a different impact on the local weather, leading to an increased probability of extreme weather processes (Yao and Luo, 2014; Yao et al.,2016; Li et al., 2020). In fact, the defining criterion for an NAO or blocking event is the unified description standard for their most common characteristics. However, each individual case has its own deviations from this criterion. At the same time, each individual case is likely to be accompanied by other systems when it occurs, and such deviations and combined impacts between systems are probably the most important nonlinear and uncertain mechanisms leading to the occurrence of extreme weather, as revealed by previous research(Luo et al., 2015; Yao et al., 2016; Li et al., 2020). The collective processes during this cold event are further described in section 5.

    Meanwhile, the positive and negative SAT anomalies that accompanied the UB geopotential height anomalies can also be examined (Fig. 2, shading). It can be seen that positive (negative) geopotential height anomalies were associated with positive (negative) SAT anomalies, which is the classical characteristic of temperature distribution due to blocking circulation as presented in previous studies (Luo et al., 2016; Yao and Luo, 2018). The positive SAT anomaly at high latitudes around the KESS region would have favored a decline in sea ice during the UB life cycle on the synoptic time scale (Luo et al., 2017), which is demonstrated in the following section. As the blocking developed,the negative SAT anomaly associated with the cyclonic center, shown in Fig. 2, moved slowly from the Siberian region to the middle and lower latitudes, eventually causing significant cooling to occur in areas of east-central and southern China on 6 January. In fact, from the end of December to the 8th of January, China experienced two extreme cold events, which has been indicated by Zheng et al. (2021).This will be mentioned later in the further analysis.

    Fig. 2. Evolutions of daily spatial distributions of 500-hPa geopotential height (contours, CI = 80 gpm) and SAT (shading, units: K) anomalies from 29 December 2020 to 9 January 2021.

    3.2. Source and pathways of the cold air

    Fig. 3. Spatial patterns of daily cold air mass (shading; in hPa) and its flux vector (arrows; in hPa m s-1) from 29 December 2020 to 9 January 2021.Coldair mass,DP=PS-P(θT),is defined as thepressure difference between the 280-K potential temperature surface P(θT)and theearth’ssurface PS,andits fluxvectorcan be expressedas

    To provide a more visual representation of the origin of this cold air and its pathways, we calculated the cold air flux to show the movement of the cold air mass in Fig. 3. It can be seen that at the end of December the source of the cold air (shaded in red) was over the Arctic, at around 120°E. Then, beginning on 1 January, the cold air started to move gradually towards the lower latitudes. The direction of movement of the cold air can be observed through the wind vectors (arrows) in Fig. 3. By 3 January, the cold air had moved pointedly into the middle and high latitudes and was distributed between 90° and 120°E. Then, from 5 January onwards, the cold air mass split from the higher latitudes and moved rapidly to the middle and lower latitudes.From 6 to 9 January, it can be seen that a large, relatively independent cold air mass had emerged at 120°E in the middle and low latitudes. This cold air mass was entrenched over the central-eastern and southern parts of China, causing an outbreak of extremely cold weather as shown in Fig. 1. It should be pointed out that on 29 December 2020, a strong mass of cold air had settled over northeast China,which caused a prominent cooling process to take place over a wide area of the country. In other words, about a week before this extreme cold event (6 Jan), an appreciable cooling process had already been experienced by most of China, therefore, this extreme cold wave was preceded by the combined impacts of at least one or even a few cold waves. Several cooling processes were closely associated with the UB, and these coactive processes are examined in detail together in section 5.

    4. Signals of pre-autumn KESS sea ice

    To investigate the potential role of the Arctic sea-ice anomaly on the cold event, we show the sea-ice cover distributions for each month in autumn 2020 in Fig. 4. To avoid uncertainty between different datasets, two widely used seaice datasets (Hadley sea-ice data and the NSIDC dataset)were used simultaneously. The data suggest that, for the entire Arctic, the negative sea-ice cover anomaly was mainly concentrated in the KESS (75°-85°N, 60°E-180°)for each month in autumn. The negative sea-ice anomaly was most evident in October during autumn 2020. The negative sea-ice anomaly in autumn may have caused a weakening of the meridional temperature gradient in the middle and high latitudes of the Eurasian region, causing the westerly jet to weaken and become unstable in the Eurasian region,providing a favorable background condition for the excitation of future meridional circulation anomalies, as revealed in Tyrlis et al. (2019). It can be seen that the negative anomaly of the sea ice moved gradually from the KESS (Figs. 4b and e) towards the Barents-Kara Seas as the season became colder (Figs. 4c and f).

    Fig. 4. Spatial patterns of SIC anomalies (units: %) for (a, d) September, (b, e) October, and (c, f) November 2020, based on monthly SIC data from (a-c) the Met Office Hadley Centre and (d-f) NSIDC. The region marked by the black frame is KESS(75°-85°N, 60°E-180°).

    To determine whether the sea-ice cover in KESS in autumn 2020 was extremely low compared to historical values, we present the Arctic sea-ice variabilities from 1979 to the present day in Fig. 5. Figures 5a and c show the variabilities of sea-ice cover in the entire Arctic for autumn (and each month of autumn). It can be seen that the autumn seaice cover in the entire Arctic has been decreasing rapidly since 1979, due to the effects of Arctic amplification(Figs. 5a and c). The Arctic sea-ice cover in September reached its lowest value on record in 2012 as shown in Figs.5a and c (red dashed lines), which is consistent with the associated official report (https://www.nasa.gov/topics/earth/features/2012-seaicemin.html). In addition, it can be observed that the Arctic sea-ice cover in October 2020 (blue dashed lines) reached its lowest value on record since 1979.Although the Arctic sea-ice cover in autumn 2020 was very close to that of 2012 (Hadley) and 2016 (NSIDC), it was also numerically at its lowest level on record since 1979.Thus, the extremely low sea-ice anomaly in the Arctic may have provided a pre-existing signal for the atmospheric circulation to act upon in later months. However, for the Eurasian region, the sea ice in the KESS may have a closer link with the atmospheric circulation and extreme weather, according to previous studies (Luo et al., 2016, 2017; Tyrlis et al.,2019). Therefore, we show the KESS sea-ice variabilities since 1979 in Figs. 5b and d. The results show that the KESS sea-ice cover in both October and autumn 2020 exhibited their lowest recorded values since 1979. In particular,the KESS sea-ice cover in October 2020 (blue dashed lines in Figs. 5b and d) was well below the second-lowest value on record. In summary, the conclusion that the KESS seaice cover in October and autumn 2020 were at their lowest values since 1979 is solid and reliable, both for the Hadley and NSIDC sea-ice datasets. In the following section, we discuss several crucial processes that preceded the cold event on a cross-seasonal time scale and then identify their combined effects.

    5. Combined effects of UB and other climatic factors

    5.1. Cumulative effects of several consecutive and persistent UBs

    Fig. 5. Time series during 1980-2020 of the spatially averaged SIC anomaly (units: %) over (a, c) the whole Arctic and (b, d) KESS (75°-85°N, 60°E-180°). The red, blue, and purple dashed lines represent the SIC variations in September, October, and November, respectively, and the black solid line represents the autumn mean. The monthly SIC dataset used here is derived from (a, b) the Met Office Hadley Centre and (c, d) NSIDC.

    Fig. 6. Time series of (a) the 300-hPa zonal wind (units: m s-1) averaged over the Ural region(45°-65°N, 30°-120°E), where the red solid line represents autumn/winter 2020, the red dashed line represents the corresponding daily mean during 1980-2020, and the blue line is the normalized NAO index, (b) the KESS (75°-85°N, 60°E-180°) SIC derived from the monthly data of NSIDC (red line) and ERA5 (blue line) and mean SAT of 2414 stations nationwide (black line). The black (blue) dashed boxes indicate the sharp SAT (SIC) decline associated with UB. The gray shading represents the UB episodes and the vertical gray line represents the Lag-0 day for each UB episode.

    It has already been mentioned above that UBs exert an important influence on Eurasian temperatures and KESS sea-ice changes. Therefore, it is necessary to examine the frequency of UBs from autumn 2020 to January 2021. We identified five UB episodes using the identification method mentioned earlier from 1 September 2020 to 16 January 2021,as shown in Fig. 6a (gray shading). From mid-November to January 2021, three consecutive UB episodes occurred with only a few days between each episode. The Ural regional mean 300-hPa zonal wind exhibited frequent fluctuations from 1 September to 16 January (red solid line in Fig. 6a).This suggests that each UB episode contributed significantly to the weakening of the westerly zonal wind on synoptic time scales. The onset of a UB, accompanied by the establishment and strengthening of the meridional circulation,would inevitably have led to a weakening of the mid- and high-latitude regional zonal wind (Luo, 2005; Luo and Cha,2012; Yao et al., 2017). The regional mean climatic zonal wind, as shown in Fig. 6a (dashed line), slowly strengthened from the warm to the cold season. However,due to the seasonal cumulative effect of consecutive UB episodes, the zonal wind during mid-November to January became unstable and weakened dramatically, to a level much lower than the climatic mean. The weakening of the high-altitude zonal jet would have affected the stability of the polar vortex and changed the normal seasonal strengthening process of the polar vortex from autumn to winter (Tyrlis et al., 2019). This may have provided the background conditions for the early jet stream dynamics to change and eventually allow for the cold wave outbreak in early January. In addition, the NAO daily index is also plotted in Fig. 6a(blue line), from which a clear positive NAO episode(NAO+) can be observed before or at the beginning of each UB episode (gray shading). This is consistent with previous findings that outlined the physical mechanism whereby a NAO+ episode favors the development of a downstream UB (Luo et al., 2016; Yao et al., 2017; Yao and Luo, 2018).At the same time, the development of NAO+ episode usually occurs a few days ahead of blocking over Europe (Yao and Luo, 2018). We further composite the NAO and UB index based on the peak day of the four UB episodes and calculate the lead-lag correlation between them (not shown). Results show that the transition of the NAO to its positive phase typically occurs 10 days ahead of a UB episode with a correlation coefficient of 0.78, which also confirms the conclusion mentioned above that each UB is preceded by a NAO+ episode. Here, the NAO index had been negative in value since the start of December, and, although there was a brief positive phase at the end of December, it turned negative again immediately thereafter. This may be attributed to the frequent UB episodes that began in mid-November and the associated weakening of the westerly jet mentioned above.Three consecutive and long-lasting UBs provided favorable westerly jet conditions for the NAO to change from positive to negative phase and to be maintained over a long period. At the same time, the positive feedback (Luo and Cha, 2012; Yao and Luo, 2018) of the NAO- contributes to a notable weakening of the mid-to-high-latitude westerly jet in the North Atlantic and western Europe. We also composite the UB index and regional (Eurasian region) mean 300-hPa zonal wind based on the peak day of the four UB episodes (not shown). Results show that a UB occurs three days ahead of the Ural regional mean 300-hPa zonal wind with a correlation coefficient of -0.69. Thus, the seasonal cumulative and combined positive feedbacks of the NAO- and a UB will contribute to a weakening of the mid-to-high-latitude westerly jet across the North Atlantic and Eurasian region, providing favorable conditions for the merger and development of the large-scale anomalous meridional circulation illustrated in Fig. 2.

    Figure 6b shows the variations of the KESS sea ice and SAT over China. The station-observed SAT showed a normal seasonal decline from autumn to winter (black line in Fig. 6b); and within this seasonal change, four instances of dramatic SAT decline can be observed during each UB episode, which were closely linked to the frequent occurrence of cold events during early winter of 2020/21. Three significant SAT declines before 6 January set the stage for the record low temperatures on that day. The physical mechanisms related to the UB, which plays an important role in SAT declines in East Asia, especially in China, have been examined in previous studies (Yao et al., 2017; Tyrlis et al.,2019). In addition, when a UB occurs, it causes significant warming to occur at high latitudes, particularly in the Barents-Kara Seas, which causes rapid melting of sea ice on synoptic time scales (Luo et al., 2016, 2017; Tyrlis et al.,2019). In this study, however, we are concerned with the KESS region, which partially overlaps with the Barents-Kara Seas, but is a little further east and contains the East Siberian Sea. This is because the sea-ice anomalies in Figs. 4 and 5 were in the KESS region. In addition, the location of this UB episode was higher in latitude relative to the climatic norm and more extensive, so the extent of impact was different. The variation of KESS sea ice from 1 September to 16 January is shown in Fig. 6b (red and blue lines). Sea-ice coverage gradually increased as the seasonal decline in temperature in KESS progressed from September to January. Then, amid the gradual increase in KESS sea ice, there were six relatively significant short periods of decline (blue boxes), as shown in Fig. 6b. Each of these short-term reductions in KESS sea ice was accompanied by a UB episode. A UB had a more pronounced effect on the reduction in sea ice in the Barents-Kara Seas during this period (not shown). However, because the abnormal signal of sea ice in the early autumn was concentrated in the KESS region, we focus on the KESS region in this study. The data suggest that, from 1 to 16 September, accompanied by the seasonal reduction in KESS sea ice, the UB episode also made an important positive feedback to the sea-ice decline. In addition, the notable short-term sea-ice reduction during 20-26 September, due to the UB episode, may have been the main reason for the historical sea-ice minimum in KESS. In any event, several UB episodes made an important contribution(positive feedback) to the reduction in sea ice in the KESS region.

    5.2. Upward signals caused by UBs

    Fig. 7. Horizontal patterns of vaTa (units: K m s-1) at the 100-hPa pressure level for (a) 16 November 2020, (b) 10 December 2020, and (c) 1 January 2021.

    Several UB episodes preceding the onset of the recordbreaking cold event not only caused a reduction in KESS sea ice and weakened the regional westerly jet, but may also have affected the stratosphere through modulating the upward progress of planetary-scale waves (Kodera et al.,2013; Shen et al., 2020), which in turn affected the stability of the polar vortex. Given that the meridional EHTv′T′is directly proportional to the vertical component of the EP flux,we, therefore, show the horizontal distribution of EHTv′T′at 100 hPa in Fig. 7 to explore the vertical propagation of wave energy associated with several crucial UB episodes.The results show that, during or after each UB episode, a significant positive (upward) EHT could be observed over the Ural region (near 60°E to 90°E), especially the UB episode that occurred in January 2021. This suggests that UB events can play an important role in the upward propagation of EP fluxes from the troposphere to the stratosphere over the Ural region. Meanwhile, anomalous signals of EP fluxes in the stratosphere are often followed by a UB by a few days to a week, as indicated in Fig. 8a (contours). We composite the associated blocking index and Fpbased on the peak day of the four UB episodes (not shown). The lead-lag correlation indicated that UB is ahead of Fpby five days with a correlation coefficient of 0.77. We also calculated the lead-lag correlation between them for the entire period from 1 September to 16 January. Results show that the UB is ahead of Fpby 5 days with a correlation coefficient of 0.6 (which is above the 90% confidence level based on a two-tailed Student’s t-test). These results are consistent with related research (Tyrlis et al., 2019). Figure 8a shows the detailed vertical component of the EP flux, Fpat every vertical level,and the zonal wind tendency, Δuˉ, caused by the convergence or divergence of EP flux. The data suggest that the upward propagation of the EP flux caused by UB in October was not obvious (only up to 150 hPa). The upward EP flux due to the UB on 16 November was stronger than the former UB. Subsequently, two UB episodes in early December and January resulted in highly significant uploading of EP flux which penetrated up to the 10-hPa level, which then led to the strong aggregation of planetary-wave energy in the middle stratosphere, as shown in Fig. 8a. Meanwhile,the shading in Fig. 8a indicates the daily evolution of zonal wind tendency caused by the EP flux variance within one day. From mid-December 2020 to early January 2021, the stratosphere featured significant negative anomalies in zonal wind tendency associated with the upward propagation and stratospheric convergence of EP flux. Meanwhile, Fig. 8b shows the zonal mean u-wind anomaly at 55°N from autumn 2020 to January 2021. The data suggest that the notable upward propagation of EP flux (Fig. 8a) on 16 December and 1 January contributed greatly to the weakening of the upper-level (stratosphere) westerly jet and further weakened the zonal wind downward to the troposphere (as indicated by the arrow in Fig. 8b). The signal, which characterized the attenuated polar vortex in the stratosphere, penetrated down into the troposphere after 1 January to modulate the local weather over East Asia (Fig. 8b, arrows). This means that the UB episodes-especially the UBs that occurred during early December and January-markedly weakened the strength and stability of the polar vortex. The weakened polar vortex further stripped the cold air mass guided by the tropospheric UB, breaking up and peeling away from the main body, which in turn affected the middle and lower latitudes.

    Fig. 8. Time-pressure evolutions of the (a) domain-mean (50°-80°N, 30°-120°E) zonal wind tendency (shading, units: m s-1 d-1) and vertical component of EP flux (contours, × 106 m2 s-2 Pa),and (b) zonal mean (90°-150°E) westerly wind anomaly at 55°N. Similar to Taguchi and Hartmann (2006), the variables in (a) are scaled by a factor inversely proportional to pressure 1000/P. The vertical dashed line represents the Lag-0 days of several UB events.

    Finally, it is worth noting that the Pacific SST had been developing anomalously since the summer of 2020, which was a typical La Ni?a year (Zheng et al., 2021). The distribution of SST anomalies in the Pacific Ocean demonstrates a typical La Ni?a SST distribution (Fig. 9). The SST in the central-eastern equatorial Pacific shows a widespread negative anomaly, which contrasts sharply with the KESS warming(sea-ice decline) at high latitudes (Figs. 4 and 5). This contrast would have caused a weakening of the meridional temperature gradient in middle and high latitudes of Eurasia,favoring the development of a meridional circulation (i.e., a UB) (Yao et al., 2017; Tyrlis et al., 2019; Zheng et al.,2021). In addition, there may be anomalous oceanic or atmospheric signals that indirectly influenced the cooling processes of the 2020/21 winter. However, they will not be discussed in this study.

    6. Conclusions and discussion

    This paper focuses on the frequent occurrence of cold events during the early winter of 2020/21, in particular the record-breaking extreme cold event in early January in China. We have examined, in detail, the physical mechanisms of the processes involved (Fig. 10), mainly from the perspective of the large-scale atmospheric circulation (i.e., UB episodes in Fig. 10). At the same time, given the complex climatic background of this cold event, we have also explored other climatic factors that may have provided an abnormal background for these processes. In short, the historical minimum value of sea-ice cover in the KESS region that occurred in autumn was a pre-existing background condition, and several UB episodes from autumn to winter were specific seasonal cumulative physical processes (as indicated in Fig. 10, dashed lines and arrows). In addition, anomalously cold SSTs in the central-eastern Pacific due to an established La Ni?a, and a long-lasting NAO- event, were also important general background conditions. The seasonal cumulative and combined impacts (Fig. 10) of these processes led to the occurrence of these extreme cold events.

    The sea-ice cover across the entire Arctic was at a historic minimum value (since 1979) in autumn 2020 (especially in October). The spatial distribution of sea ice in the Arctic further suggests that the historical minimum value of sea ice was mainly attributed to the KESS region (KESS shading in Fig. 10), where sea ice-both in autumn and in October-was at its lowest since 1979. The abnormally low KESS sea-ice conditions in autumn may have provided a key precursor to the jet dynamics that occurred later that winter. This differs from many previous studies (Luo et al.,2016, 2017) that focused only on sea ice in the Barents-Kara Seas.

    At the same time, from September 2020 to January 2021, several consecutive persistent UB episodes occurred(gray flow lines in Fig. 10), especially during the cold season, which, on the one hand, significantly slowed down the seasonal intensification of the westerly jet and the polar vortex, and on the other hand created a strong tropospheric-tostratospheric EP flux upward from the Ural region (arrows a in Fig. 10), which further weakened the polar vortex (polar vortex in Fig. 10) through seasonal cumulative effects. The prolonged NAO- that persisted from December 2020 to January 2021 also led to a weakening of the mid- and high-latitude westerly jet over the North Atlantic and western Europe(NAO- in Fig. 10). In addition, on synoptic time scales,each UB caused a significant SAT decline (cold events) in China and a notable loss of sea ice in the KESS region (red and blue shading in Fig. 10). The physical mechanisms underlying the effects of the UB on sea ice and SAT have received much attention among researchers (Luo et al.,2016, 2017; Tyrlis et al., 2019). Here, we have further demonstrated the seasonal cumulative and combined effect of several UB processes across seasons. The seasonal cumulative impact of several UB episodes from early to late December 2020 led to several consecutive cooling processes in China as indicated in Fig. 6b. Several episodes of upward EP flux(Fig. 8a) and sea ice decline (Fig. 6b) due to several consecutive UBs (Fig. 6, gray vertical stripes) exert a seasonal cumulative impact upon the occurrence of several extreme cold episodes (Fig. 6b, black boxes). In the schematic diagram(Fig. 10), we use the dashed lines to represent the several UBs and associated processes that occur across seasons from September 2020 to January 2021), which is referred to as the seasonal cumulative effect in this study. In fact, several UBs also caused three significant cooling episodes between December and early January (Fig. 6b), although their specific physical processes would have differed. In addition, the three UBs and their corresponding processes provided a favorable meridional anomalous circulation(weak westerly jet) for the record-breaking extreme cold event outbreak in early January, through a seasonal cumulative effect.

    Fig. 9. Spatial distributions of December SST (units: K) for the 1979-2020 mean climatology (contours) and anomaly pattern in 2020 (shading).

    Fig. 10. Schematic diagram of the physical processes leading to the cold events in early winter 2020/21. The thick gray solid line over the Eurasian continent represents the flow line of the UB. The letter H (L) represents the anticyclonic (cyclonic) center of the UB. Red(blue) shading indicates the warming (cooling) caused by UB. KESS marks the Kara and East Siberia Seas with abnormal sea-ice cover. The upward-pointing arrows (a) indicate the propagation of EP fluxes from the troposphere to the stratosphere (weakening of the polar vortex) associated with several UB episodes. The downward-pointing arrows (b) highlight the break-up of the polar vortex and its associated several cold air outbreaks in East Asia. The light blue shading in the central-eastern Pacific represents the cold SST anomaly due to La Ni?a. NAO- indicates the negative NAO phase starting from December 2020. The dashed lines represent the several UBs and associated processes that occur across seasons (from Sep 2020 to Jan 2021), which is referred to as the seasonal cumulative effect.

    It is worth noting that there was also a long-lasting UB episode in early September and one in early October. Significant KESS sea-ice decline was observed during these two UB episodes (Fig. 6b). In general, the Arctic sea ice reaches its minimum coverage by mid-September. Here, the seasonal reduction in KESS sea ice from 1 September to 16 September included an additional reduction due to warming caused by the UB. This reduction, combined with the persistent UB episode in late September/early October,weakened the rapid seasonal increase in sea ice that would have been normally been expected in October. This was likely one of the important factors contributing to the historical sea-ice minimum in the KESS region in October(autumn), at least from an atmospheric circulation perspective.

    Finally, the combination of abnormally cold SSTs in the central-eastern Pacific (a La Ni?a year), starting in August, and an abnormal KESS warming, starting in autumn, led to a further weakening of the meridional temperature gradient across the Eurasian region. This unstable temperature structure and weakened westerly jet provided for a favorable climatic background for the establishment and intensification of an anomalous meridional circulation (i.e.,UB) in later winter. Although the winter SATs in China were warmer in 2020/21 compared to climatology because of the extreme warm condition observed during the later winter (not shown), the frequency and amplitude of SAT fluctuations (decline) during pre-winter (Dec 2020 to Jan 2021)have become greater due to several UB episodes and associated processes. However, the numerous cold events of the pre-winter and the warm anomalies experienced later that winter were not contradictory. This is likely related to the shift of the AO (NAO) from a predominantly negative phase to a positive phase in late February (February), with the positive phase AO characterized by a stable polar vortex structure that is also not conducive to the appearance of rapid disturbances. It is also highly likely that some other,stronger atmospheric or oceanic signal (e.g., the winter monsoon) intervened to influence the continuation of the cumulative seasonal effect, which needs to be further explored in the future.

    In fact, each extreme cold spell has certain commonalities, but there are also unusual aspects. These unusual aspects are similar to the nonlinear problem in atmospheric dynamics, which has a high degree of uncertainty and chaos. Nonetheless, this study provides a reference for understanding the physical mechanisms of such extreme cold weather. Also, the pre-existing anomalous climate signal in early autumn can be regarded as a potential predictor to be used in forecasting. In addition, the issue concerning the sea-ice cover minimum in autumn 2020 will need to be further examined in future work.

    Acknowledgements. The authors acknowledge the financial support from the National Natural Science Foundation of China(Grants Nos. 41975068, 41790473, and 41975099) and the Chinese Academy of Sciences Strategic Priority Research Program (Grant No. XDA19070403).

    Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

    亚洲专区字幕在线| 日本免费一区二区三区高清不卡| 美女免费视频网站| 99国产综合亚洲精品| 青草久久国产| av黄色大香蕉| 亚洲国产高清在线一区二区三| 男女下面进入的视频免费午夜| 亚洲精品粉嫩美女一区| 欧美色视频一区免费| 国产精品永久免费网站| 欧美乱色亚洲激情| 成人三级黄色视频| 天天添夜夜摸| 亚洲成人久久爱视频| 久久久久久久久免费视频了| 在线观看免费视频日本深夜| 免费电影在线观看免费观看| 51午夜福利影视在线观看| 久久久色成人| 亚洲欧美日韩高清在线视频| 全区人妻精品视频| 久久婷婷人人爽人人干人人爱| 精品国产乱码久久久久久男人| 亚洲成人久久爱视频| 免费在线观看日本一区| 亚洲专区中文字幕在线| or卡值多少钱| 在线看三级毛片| 成人欧美大片| 日韩人妻高清精品专区| 波多野结衣高清作品| 精品久久久久久久末码| 久久草成人影院| 搡老熟女国产l中国老女人| 成人国产一区最新在线观看| 精品久久久久久久毛片微露脸| 亚洲人与动物交配视频| 99热精品在线国产| 男人舔奶头视频| 长腿黑丝高跟| 亚洲色图 男人天堂 中文字幕| 日韩欧美在线二视频| 超碰成人久久| 免费看a级黄色片| 成在线人永久免费视频| 亚洲成a人片在线一区二区| 成人精品一区二区免费| 国产精品影院久久| 91av网站免费观看| 成年免费大片在线观看| 欧美乱色亚洲激情| 国产蜜桃级精品一区二区三区| 天天躁日日操中文字幕| 丰满人妻熟妇乱又伦精品不卡| 日韩精品中文字幕看吧| 免费av毛片视频| a级毛片在线看网站| 午夜福利成人在线免费观看| 国产精品99久久99久久久不卡| 久久午夜综合久久蜜桃| 久久久久国产精品人妻aⅴ院| 亚洲乱码一区二区免费版| or卡值多少钱| 精品久久蜜臀av无| 在线观看免费视频日本深夜| 九色成人免费人妻av| 香蕉国产在线看| 黄色丝袜av网址大全| 三级男女做爰猛烈吃奶摸视频| 色噜噜av男人的天堂激情| 国产成人精品无人区| 国产精品久久电影中文字幕| 动漫黄色视频在线观看| 久久久久久久午夜电影| 97超级碰碰碰精品色视频在线观看| 9191精品国产免费久久| 少妇丰满av| 亚洲aⅴ乱码一区二区在线播放| 欧美一级毛片孕妇| 久久国产精品影院| 老汉色av国产亚洲站长工具| 国产成人aa在线观看| 看片在线看免费视频| 亚洲专区中文字幕在线| 欧美日韩乱码在线| 又爽又黄无遮挡网站| 免费看光身美女| 日韩 欧美 亚洲 中文字幕| 99国产精品一区二区蜜桃av| 欧美极品一区二区三区四区| 一级黄色大片毛片| 亚洲国产欧洲综合997久久,| 国内精品久久久久久久电影| 欧美日韩瑟瑟在线播放| 88av欧美| 日本黄色片子视频| 在线十欧美十亚洲十日本专区| av国产免费在线观看| 51午夜福利影视在线观看| 熟女少妇亚洲综合色aaa.| 少妇熟女aⅴ在线视频| 成人三级黄色视频| 亚洲欧美日韩卡通动漫| 亚洲欧美日韩高清在线视频| 午夜视频精品福利| 国产精品女同一区二区软件 | 国产精品乱码一区二三区的特点| 色综合亚洲欧美另类图片| www.www免费av| xxx96com| 91字幕亚洲| 综合色av麻豆| 人妻久久中文字幕网| 国内精品一区二区在线观看| 亚洲美女视频黄频| 99精品久久久久人妻精品| 国产成人精品久久二区二区免费| 男女午夜视频在线观看| 中亚洲国语对白在线视频| 99久久成人亚洲精品观看| 又紧又爽又黄一区二区| 国产黄a三级三级三级人| 中文字幕精品亚洲无线码一区| 看黄色毛片网站| 午夜影院日韩av| 91久久精品国产一区二区成人 | 亚洲美女视频黄频| 真人一进一出gif抽搐免费| 脱女人内裤的视频| 久久久国产欧美日韩av| 国产成人影院久久av| 午夜精品一区二区三区免费看| 老汉色∧v一级毛片| 国产精品九九99| 中文字幕精品亚洲无线码一区| 好看av亚洲va欧美ⅴa在| 中文字幕人成人乱码亚洲影| 精品一区二区三区视频在线 | 国产精品综合久久久久久久免费| 噜噜噜噜噜久久久久久91| 欧美中文日本在线观看视频| 老司机在亚洲福利影院| 国产亚洲欧美在线一区二区| 人妻久久中文字幕网| 又紧又爽又黄一区二区| a级毛片a级免费在线| 国产亚洲欧美98| 天天一区二区日本电影三级| 日韩欧美三级三区| 18禁黄网站禁片午夜丰满| 免费在线观看视频国产中文字幕亚洲| 日本 av在线| 一本久久中文字幕| 99久久精品一区二区三区| 午夜免费成人在线视频| 亚洲五月婷婷丁香| 一本综合久久免费| 天堂动漫精品| 日本撒尿小便嘘嘘汇集6| 又黄又爽又免费观看的视频| 波多野结衣高清作品| 两人在一起打扑克的视频| 国产精品一区二区三区四区免费观看 | 变态另类成人亚洲欧美熟女| 亚洲中文av在线| xxxwww97欧美| 全区人妻精品视频| 久久久久久国产a免费观看| 成人鲁丝片一二三区免费| 亚洲欧美精品综合一区二区三区| 淫秽高清视频在线观看| 国产高清videossex| 国内少妇人妻偷人精品xxx网站 | 久久精品91蜜桃| 免费高清视频大片| 欧美一级毛片孕妇| 动漫黄色视频在线观看| 欧美zozozo另类| 男女那种视频在线观看| 日韩免费av在线播放| 叶爱在线成人免费视频播放| 亚洲国产精品久久男人天堂| 神马国产精品三级电影在线观看| 色播亚洲综合网| 国产探花在线观看一区二区| 亚洲国产精品sss在线观看| 天天躁日日操中文字幕| 亚洲一区二区三区不卡视频| 亚洲精品美女久久久久99蜜臀| 老司机午夜十八禁免费视频| 国产亚洲精品久久久久久毛片| 午夜福利免费观看在线| 最近最新中文字幕大全电影3| 大型黄色视频在线免费观看| 亚洲精品久久国产高清桃花| www日本在线高清视频| 久久久国产欧美日韩av| 国产真实乱freesex| 久久精品亚洲精品国产色婷小说| 日本精品一区二区三区蜜桃| 亚洲第一欧美日韩一区二区三区| 国产爱豆传媒在线观看| 午夜福利在线观看吧| 国产精品电影一区二区三区| 色噜噜av男人的天堂激情| 久9热在线精品视频| 亚洲专区字幕在线| 久久香蕉国产精品| 在线播放国产精品三级| 久久久国产欧美日韩av| 香蕉国产在线看| 欧美绝顶高潮抽搐喷水| 国产午夜精品久久久久久| 国产v大片淫在线免费观看| 国产精品久久久久久久电影 | 亚洲欧美激情综合另类| 成人特级黄色片久久久久久久| 欧美性猛交╳xxx乱大交人| 亚洲avbb在线观看| 亚洲 国产 在线| 亚洲精华国产精华精| 国产成人精品久久二区二区免费| 亚洲激情在线av| 国产视频内射| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久末码| 两性午夜刺激爽爽歪歪视频在线观看| 久久香蕉国产精品| 黑人欧美特级aaaaaa片| 啦啦啦免费观看视频1| 一级作爱视频免费观看| 亚洲中文字幕日韩| 亚洲黑人精品在线| 亚洲一区二区三区不卡视频| 亚洲精品美女久久久久99蜜臀| 婷婷丁香在线五月| 欧美日本亚洲视频在线播放| 99在线视频只有这里精品首页| 欧美激情久久久久久爽电影| 成在线人永久免费视频| 日日夜夜操网爽| 久久精品国产亚洲av香蕉五月| 亚洲成人久久爱视频| av福利片在线观看| 成人午夜高清在线视频| 久久精品国产综合久久久| 香蕉av资源在线| 国产精品亚洲av一区麻豆| 国产99白浆流出| 亚洲欧美日韩卡通动漫| 波多野结衣巨乳人妻| 久久99热这里只有精品18| 99精品久久久久人妻精品| 国产精品98久久久久久宅男小说| 亚洲九九香蕉| 国产激情欧美一区二区| 国产成年人精品一区二区| 成熟少妇高潮喷水视频| 无遮挡黄片免费观看| 久久九九热精品免费| 老司机深夜福利视频在线观看| 亚洲熟妇中文字幕五十中出| 久久午夜综合久久蜜桃| 黄色女人牲交| 老司机福利观看| 日韩精品青青久久久久久| 黄片大片在线免费观看| 久久国产精品影院| 日韩av在线大香蕉| 成人永久免费在线观看视频| 亚洲熟妇熟女久久| 国产成人av激情在线播放| 亚洲片人在线观看| 免费在线观看亚洲国产| 一进一出抽搐gif免费好疼| 欧美成人一区二区免费高清观看 | 性色avwww在线观看| 特级一级黄色大片| 叶爱在线成人免费视频播放| 高清在线国产一区| 一区二区三区高清视频在线| 老鸭窝网址在线观看| 久久国产精品人妻蜜桃| 亚洲真实伦在线观看| 毛片女人毛片| 亚洲熟女毛片儿| 亚洲aⅴ乱码一区二区在线播放| 99久久99久久久精品蜜桃| 欧美黄色片欧美黄色片| 国产人伦9x9x在线观看| av国产免费在线观看| 很黄的视频免费| 三级国产精品欧美在线观看 | 久久久久久九九精品二区国产| 中亚洲国语对白在线视频| 搞女人的毛片| 少妇裸体淫交视频免费看高清| 看片在线看免费视频| 美女高潮喷水抽搐中文字幕| 最新美女视频免费是黄的| 一级a爱片免费观看的视频| 热99re8久久精品国产| 久久精品91蜜桃| 亚洲国产欧美一区二区综合| 老司机在亚洲福利影院| 国产精品久久久久久人妻精品电影| 9191精品国产免费久久| 噜噜噜噜噜久久久久久91| 脱女人内裤的视频| 国产av不卡久久| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 国产久久久一区二区三区| 美女大奶头视频| 日本撒尿小便嘘嘘汇集6| 久久久久久九九精品二区国产| 午夜福利视频1000在线观看| 夜夜看夜夜爽夜夜摸| 中亚洲国语对白在线视频| 国产精品乱码一区二三区的特点| 国产黄色小视频在线观看| 少妇的逼水好多| 床上黄色一级片| 久久久久久久久免费视频了| 国产成年人精品一区二区| 黄色丝袜av网址大全| 国产精品香港三级国产av潘金莲| 欧美日韩亚洲国产一区二区在线观看| 午夜福利欧美成人| 成人亚洲精品av一区二区| 亚洲自拍偷在线| 午夜激情福利司机影院| 国内少妇人妻偷人精品xxx网站 | 国产真人三级小视频在线观看| 两人在一起打扑克的视频| www.自偷自拍.com| 三级毛片av免费| 久久久久久九九精品二区国产| 久久午夜亚洲精品久久| 两个人的视频大全免费| 在线观看免费午夜福利视频| 日本免费一区二区三区高清不卡| 女警被强在线播放| 国产精品影院久久| 两人在一起打扑克的视频| 久久精品人妻少妇| 很黄的视频免费| 嫩草影视91久久| 久久久国产成人免费| 夜夜看夜夜爽夜夜摸| 亚洲精品久久国产高清桃花| av福利片在线观看| 俺也久久电影网| 黑人巨大精品欧美一区二区mp4| 国产精品 国内视频| 白带黄色成豆腐渣| 亚洲精华国产精华精| 曰老女人黄片| 欧美xxxx黑人xx丫x性爽| 欧美不卡视频在线免费观看| 亚洲精华国产精华精| 国产乱人伦免费视频| 99国产精品一区二区蜜桃av| 一个人看的www免费观看视频| 69av精品久久久久久| 日本熟妇午夜| 亚洲欧美日韩卡通动漫| 国产精品一区二区精品视频观看| 欧美xxxx黑人xx丫x性爽| 国产精品99久久99久久久不卡| 免费搜索国产男女视频| 久久国产精品影院| 国产亚洲av嫩草精品影院| 精品国产亚洲在线| 国产精品久久电影中文字幕| 久久99热这里只有精品18| 日本在线视频免费播放| 18禁国产床啪视频网站| 国内久久婷婷六月综合欲色啪| 亚洲片人在线观看| 欧美大码av| 在线永久观看黄色视频| 午夜成年电影在线免费观看| 丰满的人妻完整版| 精品久久久久久久久久免费视频| 长腿黑丝高跟| 国产激情偷乱视频一区二区| 成在线人永久免费视频| 日本熟妇午夜| 中亚洲国语对白在线视频| 一级黄色大片毛片| 日韩av在线大香蕉| 精品不卡国产一区二区三区| 午夜免费激情av| 在线a可以看的网站| 十八禁人妻一区二区| 一区二区三区激情视频| 久久久色成人| 此物有八面人人有两片| 色视频www国产| 午夜福利18| 女人被狂操c到高潮| 成人三级黄色视频| 欧美国产日韩亚洲一区| 久久久久久久久久黄片| 窝窝影院91人妻| 亚洲无线在线观看| 国产伦一二天堂av在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩精品网址| 免费观看的影片在线观看| 成年女人看的毛片在线观看| 可以在线观看的亚洲视频| 97碰自拍视频| 国产v大片淫在线免费观看| 在线观看一区二区三区| 国产 一区 欧美 日韩| 亚洲熟妇中文字幕五十中出| 色综合婷婷激情| 我要搜黄色片| 国产乱人伦免费视频| 99riav亚洲国产免费| 999久久久精品免费观看国产| 国产真实乱freesex| 日韩精品青青久久久久久| 欧美丝袜亚洲另类 | 美女扒开内裤让男人捅视频| 香蕉av资源在线| 欧美一级a爱片免费观看看| 哪里可以看免费的av片| 少妇的逼水好多| 日韩欧美免费精品| 国产精品1区2区在线观看.| 一本久久中文字幕| 亚洲专区中文字幕在线| 国产单亲对白刺激| 丰满的人妻完整版| 成人一区二区视频在线观看| 欧美性猛交╳xxx乱大交人| 亚洲最大成人中文| 一本久久中文字幕| 欧美+亚洲+日韩+国产| 亚洲一区高清亚洲精品| 久久亚洲精品不卡| 久久国产精品人妻蜜桃| 欧美性猛交黑人性爽| 一本综合久久免费| 午夜福利在线观看吧| 欧美日韩黄片免| 午夜激情欧美在线| 午夜福利欧美成人| 精品久久久久久久久久免费视频| 90打野战视频偷拍视频| 欧美激情久久久久久爽电影| 欧美日本亚洲视频在线播放| 国产视频内射| 巨乳人妻的诱惑在线观看| 99久久成人亚洲精品观看| 亚洲国产欧洲综合997久久,| 91九色精品人成在线观看| 最近最新中文字幕大全免费视频| 国产亚洲精品久久久久久毛片| 麻豆国产av国片精品| 精品人妻1区二区| 午夜福利成人在线免费观看| 88av欧美| 韩国av一区二区三区四区| 99国产精品一区二区蜜桃av| 国产一级毛片七仙女欲春2| 狂野欧美白嫩少妇大欣赏| 美女大奶头视频| 国产高清videossex| 国产黄a三级三级三级人| 一本综合久久免费| 久久久久国内视频| 国产乱人视频| 精品日产1卡2卡| 国产精品,欧美在线| 窝窝影院91人妻| 少妇裸体淫交视频免费看高清| 黄色日韩在线| 操出白浆在线播放| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全电影3| 国产三级黄色录像| 欧美日韩瑟瑟在线播放| 欧美成狂野欧美在线观看| 97碰自拍视频| 法律面前人人平等表现在哪些方面| 亚洲18禁久久av| 韩国av一区二区三区四区| 母亲3免费完整高清在线观看| 国产99白浆流出| 两人在一起打扑克的视频| 九九久久精品国产亚洲av麻豆 | 热99在线观看视频| 国产精品电影一区二区三区| 亚洲欧美一区二区三区黑人| 欧美三级亚洲精品| 日本黄色片子视频| 国产91精品成人一区二区三区| 波多野结衣巨乳人妻| 亚洲欧美日韩卡通动漫| 黄色 视频免费看| 欧美中文日本在线观看视频| 美女黄网站色视频| 中文资源天堂在线| 日韩高清综合在线| 午夜福利在线在线| av黄色大香蕉| 天天添夜夜摸| 99视频精品全部免费 在线 | 男人舔女人下体高潮全视频| 脱女人内裤的视频| 最近最新中文字幕大全免费视频| 色哟哟哟哟哟哟| 又爽又黄无遮挡网站| 美女黄网站色视频| 好看av亚洲va欧美ⅴa在| 又大又爽又粗| 婷婷精品国产亚洲av在线| svipshipincom国产片| 亚洲精品国产精品久久久不卡| 一个人看视频在线观看www免费 | 欧美一级a爱片免费观看看| 国内毛片毛片毛片毛片毛片| 亚洲国产高清在线一区二区三| xxxwww97欧美| 美女黄网站色视频| 国产精品亚洲美女久久久| 午夜免费激情av| 国产三级在线视频| 亚洲天堂国产精品一区在线| 99国产综合亚洲精品| 国产亚洲精品久久久久久毛片| aaaaa片日本免费| 香蕉av资源在线| netflix在线观看网站| 免费看十八禁软件| 免费看a级黄色片| 国产aⅴ精品一区二区三区波| 综合色av麻豆| 黑人巨大精品欧美一区二区mp4| 免费在线观看亚洲国产| 国产91精品成人一区二区三区| 国产熟女xx| 一本久久中文字幕| 日韩欧美精品v在线| 国产又黄又爽又无遮挡在线| 成在线人永久免费视频| 日本黄色视频三级网站网址| 狂野欧美激情性xxxx| 亚洲精品美女久久久久99蜜臀| 少妇熟女aⅴ在线视频| 国产欧美日韩一区二区精品| 国产不卡一卡二| 丁香六月欧美| 色播亚洲综合网| 久久精品人妻少妇| 国产人伦9x9x在线观看| 国产精品女同一区二区软件 | 在线免费观看不下载黄p国产 | 亚洲 欧美一区二区三区| cao死你这个sao货| 久久久久久九九精品二区国产| 国产精品久久电影中文字幕| 露出奶头的视频| 熟女电影av网| 亚洲精品456在线播放app | 久久精品国产综合久久久| 99riav亚洲国产免费| 亚洲成人精品中文字幕电影| 久久久久久久久免费视频了| av视频在线观看入口| 日韩人妻高清精品专区| 久久久久精品国产欧美久久久| 欧美性猛交╳xxx乱大交人| 精品一区二区三区视频在线观看免费| 久久精品国产99精品国产亚洲性色| 天堂动漫精品| 亚洲色图av天堂| 久久这里只有精品19| 美女黄网站色视频| 亚洲最大成人中文| 黄色片一级片一级黄色片| 午夜免费观看网址| 欧美绝顶高潮抽搐喷水| 国产真实乱freesex| 国产蜜桃级精品一区二区三区| 国产熟女xx| 1024手机看黄色片| 亚洲美女视频黄频| 久久久久免费精品人妻一区二区| 色在线成人网| 中文字幕久久专区| 老鸭窝网址在线观看| 国产男靠女视频免费网站| 真实男女啪啪啪动态图| 国产高清有码在线观看视频| 成人国产一区最新在线观看| 最新美女视频免费是黄的| 91在线精品国自产拍蜜月 | 天堂动漫精品| 欧美日韩精品网址| 国产97色在线日韩免费| 99在线人妻在线中文字幕| 变态另类成人亚洲欧美熟女| 首页视频小说图片口味搜索| 九色成人免费人妻av| 女同久久另类99精品国产91| 中出人妻视频一区二区| 九九热线精品视视频播放| www日本黄色视频网| 国产成人精品久久二区二区免费|