• <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视频精品全部免费 在线| 成人毛片a级毛片在线播放| 国产黄频视频在线观看| 久久女婷五月综合色啪小说| 久久99一区二区三区| 婷婷色综合大香蕉| 偷拍熟女少妇极品色| 热99国产精品久久久久久7| 国产精品久久久久久久久免| 国产日韩欧美在线精品| 国产伦理片在线播放av一区| 亚洲精品aⅴ在线观看| 日本av免费视频播放| 午夜免费男女啪啪视频观看| 国产精品久久久久久av不卡| 亚洲精品日韩av片在线观看| 免费看不卡的av| 夜夜骑夜夜射夜夜干| 欧美精品人与动牲交sv欧美| 80岁老熟妇乱子伦牲交| 热re99久久精品国产66热6| 夫妻午夜视频| 国产av精品麻豆| 成年美女黄网站色视频大全免费 | 亚洲精品成人av观看孕妇| 国产极品天堂在线| 一级av片app| 国产日韩一区二区三区精品不卡 | 91精品伊人久久大香线蕉| 国产av精品麻豆| 欧美成人午夜免费资源| 人人妻人人澡人人看| 另类亚洲欧美激情| 大香蕉97超碰在线| 亚洲伊人久久精品综合| 国产精品人妻久久久久久| 久久韩国三级中文字幕| 99热全是精品| 夫妻性生交免费视频一级片| 一本一本综合久久| 亚洲怡红院男人天堂| 亚洲人成网站在线观看播放| 亚洲国产最新在线播放| 亚洲真实伦在线观看| 久久久久久久亚洲中文字幕| 晚上一个人看的免费电影| 麻豆精品久久久久久蜜桃| 成年人午夜在线观看视频| 一本久久精品| 亚洲av男天堂| 亚洲成色77777| 国产精品蜜桃在线观看| 日日爽夜夜爽网站| 中文精品一卡2卡3卡4更新| 日韩欧美 国产精品| 久久久久国产网址| 99re6热这里在线精品视频| 两个人免费观看高清视频 | 亚洲成人手机| 国产精品人妻久久久影院| 国产极品天堂在线| 精品人妻偷拍中文字幕| 在线天堂最新版资源| 精品少妇黑人巨大在线播放| 中文字幕久久专区| 亚洲精品456在线播放app| 特大巨黑吊av在线直播| 极品教师在线视频| 晚上一个人看的免费电影| 亚洲va在线va天堂va国产| 六月丁香七月| 国精品久久久久久国模美| 激情五月婷婷亚洲| 国产精品欧美亚洲77777| 日本vs欧美在线观看视频 | 中文乱码字字幕精品一区二区三区| 国产黄片视频在线免费观看| 日韩三级伦理在线观看| 欧美日韩精品成人综合77777| 寂寞人妻少妇视频99o| 人妻一区二区av| 久久久a久久爽久久v久久| 欧美成人精品欧美一级黄| 亚洲国产色片| 高清欧美精品videossex| 91精品伊人久久大香线蕉| 久久人人爽人人爽人人片va| 一区二区av电影网| 久久久国产一区二区| 99热国产这里只有精品6| 国产欧美亚洲国产| 国产在线视频一区二区| 欧美老熟妇乱子伦牲交| 日日摸夜夜添夜夜添av毛片| 97在线人人人人妻| 丰满乱子伦码专区| 亚洲自偷自拍三级| 久久久久久伊人网av| 国产伦精品一区二区三区视频9| av在线观看视频网站免费| 久久鲁丝午夜福利片| 丝瓜视频免费看黄片| 中文字幕免费在线视频6| 久久精品国产亚洲网站| 亚洲色图综合在线观看| 久久国产亚洲av麻豆专区| 国国产精品蜜臀av免费| 极品人妻少妇av视频| 日日啪夜夜撸| 午夜日本视频在线| 免费高清在线观看视频在线观看| 熟女电影av网| 成人亚洲精品一区在线观看| 精品久久久噜噜| 亚洲av在线观看美女高潮| 晚上一个人看的免费电影| 亚洲综合色惰| 国产高清不卡午夜福利| 亚洲国产精品999| 99国产精品免费福利视频| 久久av网站| 久久人妻熟女aⅴ| 亚洲精品乱码久久久久久按摩| 日本色播在线视频| 18禁裸乳无遮挡动漫免费视频| 国国产精品蜜臀av免费| 欧美国产精品一级二级三级 | 亚洲成人手机| 午夜福利影视在线免费观看| 久久国产亚洲av麻豆专区| 国精品久久久久久国模美| 91在线精品国自产拍蜜月| 久久午夜福利片| 麻豆成人午夜福利视频| 中文字幕制服av| 精品视频人人做人人爽| 人妻制服诱惑在线中文字幕| 午夜福利在线观看免费完整高清在| 日本与韩国留学比较| 搡女人真爽免费视频火全软件| 日本av免费视频播放| 最新的欧美精品一区二区| 欧美日韩一区二区视频在线观看视频在线| h日本视频在线播放| 国产男女超爽视频在线观看| 另类亚洲欧美激情| av卡一久久| 日韩欧美精品免费久久| 国产精品一区二区在线不卡| 色视频在线一区二区三区| 久久久亚洲精品成人影院| 免费看光身美女| 国产毛片在线视频| av有码第一页| 制服丝袜香蕉在线| 男人舔奶头视频| 国产黄频视频在线观看| 高清黄色对白视频在线免费看 | 国产欧美日韩综合在线一区二区 | 午夜福利网站1000一区二区三区| 国产av精品麻豆| 日韩av不卡免费在线播放| 国产极品天堂在线| 国产熟女午夜一区二区三区 | 日韩一本色道免费dvd| 久久免费观看电影| 国产色爽女视频免费观看| 国产日韩欧美亚洲二区| 亚洲三级黄色毛片| av在线老鸭窝| 一级毛片aaaaaa免费看小| 秋霞伦理黄片| 日韩电影二区| 免费人成在线观看视频色| 在线观看一区二区三区激情| 亚洲欧美成人精品一区二区| 久久这里有精品视频免费| 高清毛片免费看| 精品亚洲乱码少妇综合久久| 亚洲美女搞黄在线观看| 国产欧美日韩一区二区三区在线 | av卡一久久| 毛片一级片免费看久久久久| 偷拍熟女少妇极品色| 桃花免费在线播放| 日韩三级伦理在线观看| 黑人高潮一二区| 日韩av免费高清视频| 亚洲自偷自拍三级| 国产亚洲最大av| 久久精品久久精品一区二区三区| 欧美日韩视频高清一区二区三区二| 亚洲国产精品一区三区| 少妇裸体淫交视频免费看高清| 男女国产视频网站| 观看美女的网站| 欧美日韩国产mv在线观看视频| 六月丁香七月| 麻豆成人午夜福利视频| 简卡轻食公司| 国产精品伦人一区二区| 极品教师在线视频| 精品人妻一区二区三区麻豆| 三上悠亚av全集在线观看 | 一级,二级,三级黄色视频| 91午夜精品亚洲一区二区三区| 边亲边吃奶的免费视频| av网站免费在线观看视频| 成人毛片a级毛片在线播放| 美女大奶头黄色视频| 美女大奶头黄色视频| 成人免费观看视频高清| h日本视频在线播放| 日韩,欧美,国产一区二区三区| 九草在线视频观看| av福利片在线观看| 国产免费福利视频在线观看| 国产精品久久久久久久电影| 夜夜看夜夜爽夜夜摸| 99re6热这里在线精品视频| 国产成人freesex在线| 人妻制服诱惑在线中文字幕| 精品酒店卫生间| 在线天堂最新版资源| 寂寞人妻少妇视频99o| 22中文网久久字幕| av免费观看日本| 最近最新中文字幕免费大全7| 免费久久久久久久精品成人欧美视频 | 久久久国产欧美日韩av| 97在线人人人人妻| 久久久欧美国产精品| 内射极品少妇av片p| 看十八女毛片水多多多| 久久久久久久久久久久大奶| 99热国产这里只有精品6| 免费av中文字幕在线| 午夜福利影视在线免费观看| 国产淫片久久久久久久久| 超碰97精品在线观看| 老熟女久久久| 不卡视频在线观看欧美| 黄色视频在线播放观看不卡| 老司机亚洲免费影院| 成年人午夜在线观看视频| 中文字幕制服av| 国产伦在线观看视频一区| 夜夜看夜夜爽夜夜摸| 在线天堂最新版资源| 亚洲精品自拍成人| 草草在线视频免费看| 91久久精品国产一区二区三区| 观看美女的网站| 精品人妻熟女av久视频| 亚洲国产欧美日韩在线播放 | 色婷婷av一区二区三区视频| 五月玫瑰六月丁香| 日日摸夜夜添夜夜爱| 一本一本综合久久| 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| 久久婷婷青草| 看免费成人av毛片| 午夜日本视频在线| 日日撸夜夜添| av天堂久久9| 中文字幕精品免费在线观看视频 | av天堂中文字幕网| 日产精品乱码卡一卡2卡三| 国产 一区精品| 一本—道久久a久久精品蜜桃钙片| 久久综合国产亚洲精品| 久久国内精品自在自线图片| 亚洲精品一区蜜桃| h日本视频在线播放| 国产午夜精品久久久久久一区二区三区| 偷拍熟女少妇极品色| 亚洲va在线va天堂va国产| 亚洲欧美日韩另类电影网站| 熟妇人妻不卡中文字幕| 国产精品国产av在线观看| 日韩熟女老妇一区二区性免费视频| 2022亚洲国产成人精品| 视频区图区小说| 日本猛色少妇xxxxx猛交久久| 综合色丁香网| 我要看黄色一级片免费的| 国产欧美日韩一区二区三区在线 | 亚洲国产欧美日韩在线播放 | 91aial.com中文字幕在线观看| av又黄又爽大尺度在线免费看| 免费大片18禁| 免费观看无遮挡的男女| 久久99热6这里只有精品| 少妇精品久久久久久久| 亚洲国产日韩一区二区| 男人舔奶头视频| 2018国产大陆天天弄谢| 一区二区三区乱码不卡18| 美女视频免费永久观看网站| 人妻 亚洲 视频| 亚洲精品aⅴ在线观看| 日韩伦理黄色片| 一区二区三区四区激情视频| 国产乱来视频区| 日日摸夜夜添夜夜添av毛片| 亚洲精品日韩av片在线观看| a级毛片免费高清观看在线播放| 王馨瑶露胸无遮挡在线观看| 久久人人爽人人片av| 中文字幕人妻熟人妻熟丝袜美| a级片在线免费高清观看视频| 人人妻人人澡人人爽人人夜夜| 亚洲精品乱码久久久久久按摩| 亚洲av成人精品一二三区| 久久久精品94久久精品| 精品卡一卡二卡四卡免费| 水蜜桃什么品种好| tube8黄色片| 伦理电影大哥的女人| 久久久亚洲精品成人影院| 中国三级夫妇交换| 91久久精品国产一区二区成人| 国产美女午夜福利| 日韩视频在线欧美| 大陆偷拍与自拍| 国产熟女午夜一区二区三区 | 2022亚洲国产成人精品| 美女国产视频在线观看| av一本久久久久| 老司机亚洲免费影院| 纯流量卡能插随身wifi吗| 99久久精品一区二区三区| 国产伦在线观看视频一区| 天堂中文最新版在线下载| 久久国产乱子免费精品| 一级a做视频免费观看| 美女福利国产在线| 亚洲成人一二三区av| 一级毛片aaaaaa免费看小| 黄色日韩在线| 成人国产麻豆网| 亚洲色图综合在线观看| 亚洲人与动物交配视频| 亚洲精品一区蜜桃| 亚洲av二区三区四区| av在线观看视频网站免费| 大片电影免费在线观看免费| 女性被躁到高潮视频| 99热全是精品| 国产精品人妻久久久久久| 欧美日韩一区二区视频在线观看视频在线| 美女国产视频在线观看| 黑丝袜美女国产一区| 亚洲欧美精品自产自拍| 18禁在线播放成人免费| 欧美日韩亚洲高清精品| 亚洲欧美成人精品一区二区| 男女免费视频国产| xxx大片免费视频| 亚洲欧美一区二区三区国产| 久久婷婷青草| a级毛片在线看网站| 国产国拍精品亚洲av在线观看| 久久国内精品自在自线图片| 涩涩av久久男人的天堂| 在线精品无人区一区二区三| 欧美精品一区二区大全| 亚洲av成人精品一二三区| 免费大片18禁| 纯流量卡能插随身wifi吗| 十八禁网站网址无遮挡 | 9色porny在线观看| 亚洲精品456在线播放app| 国产精品99久久久久久久久| 国产黄色免费在线视频| 一级片'在线观看视频| 亚洲va在线va天堂va国产| 国产日韩欧美亚洲二区| 成年人免费黄色播放视频 | 偷拍熟女少妇极品色| 日本黄色片子视频| 亚洲欧美成人精品一区二区| 交换朋友夫妻互换小说| 中国三级夫妇交换| 亚洲欧美中文字幕日韩二区| 亚洲av在线观看美女高潮| 女人久久www免费人成看片| 亚洲欧美成人综合另类久久久| 国产av国产精品国产| 亚洲人成网站在线播| av天堂中文字幕网| 一级二级三级毛片免费看| 丰满饥渴人妻一区二区三| 精品视频人人做人人爽| 午夜福利网站1000一区二区三区| 亚洲欧洲日产国产| 婷婷色av中文字幕| 全区人妻精品视频| 日韩av免费高清视频| 狂野欧美白嫩少妇大欣赏| 欧美高清成人免费视频www| 免费人妻精品一区二区三区视频| 亚洲av电影在线观看一区二区三区| 男人和女人高潮做爰伦理| 蜜桃在线观看..| 五月伊人婷婷丁香| av国产久精品久网站免费入址| 亚洲美女搞黄在线观看| 成人免费观看视频高清| 国产欧美日韩综合在线一区二区 | 精品少妇内射三级| 熟妇人妻不卡中文字幕| 欧美日本中文国产一区发布| 国产熟女午夜一区二区三区 | a级毛片免费高清观看在线播放| 亚洲国产精品成人久久小说| 一级爰片在线观看| 成人毛片a级毛片在线播放| 日韩av免费高清视频| 国产 精品1| 日韩在线高清观看一区二区三区| 久久久亚洲精品成人影院| 久久久精品94久久精品| 免费黄色在线免费观看| 国产在线男女| av免费观看日本| 精品久久久精品久久久| 少妇精品久久久久久久| 插逼视频在线观看| 国产精品三级大全| 欧美3d第一页| 成人亚洲欧美一区二区av| 亚洲欧美一区二区三区国产| videos熟女内射| 日韩av在线免费看完整版不卡| 校园人妻丝袜中文字幕| 九九爱精品视频在线观看| 最近中文字幕2019免费版| 精品久久久久久久久av| 亚洲精品第二区| 少妇高潮的动态图| 日韩av免费高清视频| 日韩一区二区三区影片| 国产伦精品一区二区三区四那| 久久99一区二区三区| 成人特级av手机在线观看| 亚洲av不卡在线观看| 你懂的网址亚洲精品在线观看| 五月玫瑰六月丁香| 人妻系列 视频| www.av在线官网国产| 色婷婷av一区二区三区视频| 高清欧美精品videossex| 十八禁高潮呻吟视频 | 9色porny在线观看| 蜜桃久久精品国产亚洲av| 全区人妻精品视频| xxx大片免费视频| 国产精品久久久久成人av| 亚洲丝袜综合中文字幕| 亚洲精品久久久久久婷婷小说| 人人澡人人妻人| 欧美精品高潮呻吟av久久| 汤姆久久久久久久影院中文字幕| 91在线精品国自产拍蜜月| 国产精品麻豆人妻色哟哟久久| 亚洲av中文av极速乱| 日韩成人av中文字幕在线观看| 国产精品蜜桃在线观看| 日本免费在线观看一区| 精品一区在线观看国产| h日本视频在线播放| 啦啦啦在线观看免费高清www| 青春草国产在线视频| 色94色欧美一区二区| 有码 亚洲区| 嘟嘟电影网在线观看| 欧美+日韩+精品| 国产成人91sexporn| 日韩 亚洲 欧美在线| 成年人免费黄色播放视频 | 插阴视频在线观看视频| 精品一品国产午夜福利视频| 男女免费视频国产| 亚洲欧美日韩卡通动漫| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美 国产精品| 在线观看国产h片| 国产在线一区二区三区精| 国内少妇人妻偷人精品xxx网站| 久久ye,这里只有精品| 国产精品久久久久久精品古装| 日本黄色日本黄色录像| 色吧在线观看| 国产视频内射| 免费观看在线日韩| 一区二区三区四区激情视频| 少妇精品久久久久久久| 免费人妻精品一区二区三区视频| 毛片一级片免费看久久久久| 欧美97在线视频| 最新中文字幕久久久久| 国产精品嫩草影院av在线观看| 久久av网站| 少妇的逼好多水| 伊人久久国产一区二区| 日韩人妻高清精品专区| 中文字幕av电影在线播放| 国产高清三级在线| 最新的欧美精品一区二区| 丰满少妇做爰视频| 赤兔流量卡办理| 我的女老师完整版在线观看| 欧美xxxx性猛交bbbb| av天堂中文字幕网| av卡一久久| 国产成人精品福利久久| 久久久久久久久大av| 日本黄色片子视频| 三级国产精品欧美在线观看| 在线观看国产h片| 大又大粗又爽又黄少妇毛片口| 成人毛片a级毛片在线播放| 热re99久久国产66热| 熟女av电影| 黑人猛操日本美女一级片| 精品人妻一区二区三区麻豆| 自拍偷自拍亚洲精品老妇| 亚洲精品国产av成人精品| 亚洲欧美日韩东京热| 免费久久久久久久精品成人欧美视频 | 人妻人人澡人人爽人人| 亚洲,欧美,日韩| 成人影院久久| 美女中出高潮动态图| 国产欧美日韩综合在线一区二区 | av在线app专区| 午夜日本视频在线| 下体分泌物呈黄色| 寂寞人妻少妇视频99o| 久久国产亚洲av麻豆专区| 成人午夜精彩视频在线观看| 国产成人免费观看mmmm| 亚洲精品国产av蜜桃| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 看十八女毛片水多多多| 久久影院123| 久久精品熟女亚洲av麻豆精品| 久久精品国产自在天天线| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 中文字幕精品免费在线观看视频 | 久久久久久久精品精品| 日韩精品免费视频一区二区三区 | 一二三四中文在线观看免费高清| 久久久精品免费免费高清| 国产伦精品一区二区三区视频9| 麻豆乱淫一区二区| 男女边摸边吃奶| 熟女人妻精品中文字幕| 成人亚洲欧美一区二区av| 久久人人爽av亚洲精品天堂| 免费人成在线观看视频色| 菩萨蛮人人尽说江南好唐韦庄| 久久99热6这里只有精品| 寂寞人妻少妇视频99o| 国产精品久久久久久av不卡| 精品亚洲成a人片在线观看| 亚洲国产精品999| 亚洲国产精品一区三区| 大又大粗又爽又黄少妇毛片口| 搡女人真爽免费视频火全软件| 少妇熟女欧美另类| 日韩av不卡免费在线播放| 2021少妇久久久久久久久久久| 午夜免费观看性视频| 日日撸夜夜添| 又爽又黄a免费视频| 亚洲天堂av无毛| 女的被弄到高潮叫床怎么办| 如日韩欧美国产精品一区二区三区 | 久久久欧美国产精品| 国产一区亚洲一区在线观看| 久久精品国产a三级三级三级| 亚洲精品日韩av片在线观看| 在线免费观看不下载黄p国产| 亚洲欧美一区二区三区黑人 | 成人无遮挡网站| 免费av不卡在线播放| 亚洲色图综合在线观看| 在线精品无人区一区二区三| 97在线视频观看| 日韩欧美一区视频在线观看 | 在线观看人妻少妇| 十分钟在线观看高清视频www | 午夜福利,免费看| 久久综合国产亚洲精品| 亚洲人成网站在线观看播放| 一级片'在线观看视频| 观看免费一级毛片| 桃花免费在线播放| 一级爰片在线观看| 精品久久国产蜜桃| 久久毛片免费看一区二区三区| 国产亚洲欧美精品永久| 少妇精品久久久久久久| 我的女老师完整版在线观看|