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

    Simple Metrics for Representing East Asian Winter Monsoon Variability: Urals Blocking and Western Pacific Teleconnection Patterns

    2016-11-14 05:55:59HoffmanCHEUNGandWenZHOU
    Advances in Atmospheric Sciences 2016年6期

    Hoffman H.N.CHEUNGand Wen ZHOU*

    1Guy Carpenter Asia-Pacific Climate Impact Centre,School of Energy and Environment, City University of Hong Kong,Hong Kong SAR

    2City University of Hong Kong Shenzhen Research Institute,Shenzhen 518057

    Simple Metrics for Representing East Asian Winter Monsoon Variability: Urals Blocking and Western Pacific Teleconnection Patterns

    Hoffman H.N.CHEUNG1,2and Wen ZHOU*1,2

    1Guy Carpenter Asia-Pacific Climate Impact Centre,School of Energy and Environment, City University of Hong Kong,Hong Kong SAR

    2City University of Hong Kong Shenzhen Research Institute,Shenzhen 518057

    Instead of conventional East Asian winter monsoon indices(EAWMIs),we simply use two large-scale teleconnection patterns to represent long-term variations in the EAWM.First,the Urals blocking pattern index(UBI)is closely related to cold air advection from the high latitudes towards western Siberia,such that it shows an implicit linkage with the Siberian high intensity and the surface air temperature(SAT)variations north of 40°N in the EAWM region.Second,the well-known western Pacific teleconnection index(WPI)is connected with the meridional displacement of the East Asian jet stream and

    the East Asian trough.This is strongly related to the SAT variations in the coastal area south of 40°N in the EAWM region. The temperature variation in the EAWM region is also represented by the two dominant temperature modes,which are called the northern temperature mode(NTM)and the southern temperature mode(STM).Compared to 19 existing EAWMIs and other well-known teleconnection patterns,the UBI shows the strongest correlation with the NTM,while the WPI shows an equally strong correlation with the STM as four EAWMIs.The UBI–NTM and WPI–STM relationships are robust when the correlation analysis is repeated by(1)the 31-year running correlation and(2)the 8-year high-pass and low-pass filter.Hence, these results are useful for analyzing the large-scale teleconnections of the EAWM and for evaluating this issue in climate models.In particular,more studies should focus on the teleconnection patterns over extratropical Eurasia.

    East Asian winter monsoon,Urals blocking high,western Pacific teleconnection pattern,teleconnection

    1.Introduction

    The intensity of the East Asian winter monsoon(EAWM) is closely related to severe conditions in densely populated regions in East Asia,including China,Japan,and Korea.A simple and predictable EAWM index(EAWMI)is important for seasonal forecasting and for studying the long-term variations of the EAWM.Based on large-scale circulation features over the EAWM region(Fig.1),many previous works have constructed different EAWMIs using:(1)the mean SLP (MSLP)gradient between East Asia and the Pacific Ocean; (2)the lower tropospheric wind along the coastal region of East Asia;(3)the East Asian trough;and(4)the upper tropospheric zonal wind shear(Wang and Chen,2010).Yet,the EAWM is a complicated system involving the air–sea interaction between the Asian continent and the Pacific Ocean, and the topographic forcing exerted by the Tibetan Plateau. The interannual variations of surface air temperature(SAT)in the northern and southern parts of the EAWM region are uncorrelated(Wang et al.,2010a).Therefore,it is challenging to construct a unified index that can successfully capture the SAT variations in the entire EAWM region.

    The evolution of the EAWM system is related to middle and upper tropospheric remote signals from both upstream and downstream(e.g.,Chang and Lau,1980;Joung and Hitchman,1982;Lau and Lau,1984;Hsu and Wallace, 1985;Takaya and Nakamura,2005,2013;Song et al.,2014). A typical severe cold-air outbreak in the EAWM region is characterized by an inverted-omega geopotential height pattern.This is composedoftwoblockingridgesovertheUrals–Siberia region and the North Pacific(Yeh et al.,1962;Ding, 1994).A recent study by Cheung et al.(2015)showed that the long-term variation in the number of cold days(not the year-to-year SAT variability)in Hong Kong,a coastal city in Southeast China,is closely related to the frequency of the Urals blocking pattern index(UBI)and the western Pacific teleconnection pattern index(WPI).In this paper,we will further show how both the UBI and WPI constitute one of the most robust signals accounting for the year-to-yearSAT variability in the EAWM region.Unlike most previous works,which tendedto focuson the relationshipbetween blocking/teleconnectionsandEAWMvariability(Gonget al., 2001;Lee and Jhun,2006;Liu et al.,2014;Lim and Kim, 2015),our analysis begins with the SAT in the EAWM region.

    ?Institute of Atmospheric Physics/Chinese Academy of Sciences,and Science Press and Springer-Verlag Berlin Heidelberg 2016

    Fig.1.Large-scale circulation features of the EAWM:(a)MSLP(shading;units:hPa)and the 1000-hPa wind(vectors;units:m s-1);and(b)U250(shading;units:m s-1)and Z500(contour interval:60 gpm).In(a),the Tibetan Plateau is colored black,while A,B and C denote the Siberian high,Aleutian low and near-surface northeasterly flow over the East Asian coast,respectively.In(b),D and E denote the East Asian jet stream and East Asian trough,respectively.

    Following the definition of Wang et al.(2010a),the SAT variation in the EAWM region will also be described by the two major temperature modes.The northern temperature mode(NTM)captures the temperature variation north of 40°N well,whereas the southern temperature mode(STM) successfully depicts the temperature variation south of 40°N. Compared to 19 EAWMIs and other well-known teleconnection indices,we will demonstrate that the UBI shows the strongest correlationwith the NTM,which is poorlycaptured by most of the EAWMIs.The WPI shows an equally strong correlation with the STM with four EAWMIs.Because the UBI–NTM and WPI–STM correlations exceed 0.7 in magnitude on both interannual and interdecadal timescales,we suggest that Urals blocking(UB)and western Pacific(WP) teleconnection are crucial for EAWM studies.Based on the findings of this study and other recent studies,we will also discuss the key to defining a representative EAWMI.

    The present study investigates 66 winters,from 1948/49 to 2013/14,using the NCEP–NCAR reanalysis datasets. These include the daily field of geopotential height(Z),the monthly field of SAT,MSLP,and the zonal and meridional components of wind(U and V).Following this introduction,the teleconnectionindices(UBI and WPI)are defined in section 2.Section 3 presents the relationship between these teleconnection indices and the temperature variations in the EAWMregion.Section4highlightsthestrengthoftheseteleconnection indices to represent the EAWM variability.The results are summarized and discussed in section 5.

    2.Definition of teleconnection indices

    The prominent large-scale circulation features related to theSAT intheEAWMregionaredepictedbyperformingcorrelation analysis for the area-averagedSAT in the EAWM region(20°–50°N and 100°–140°E;green box in Fig.2a)with the 500-hPa geopotential height(Z500;shading in Fig.2a). We use the following three steps to highlight the key regions responsible for the year-to-year variability of the SAT in the EAWM region(white contours in Fig.2a):

    ·First,we repeat the correlation analysis for each of the grids in the EAWM region(221 grids in total).The output in the analysis is r1(λ,φ,i),whereλ∈[0,357.5]°E, φ∈[0,90.0]°N and i∈[1,221],and r1is the linear correlation coefficient between Z500and the SAT grid in the EAWM region.

    Fig.2.Correlation between the area-averaged SAT in the EAWM region(green box)and(a)the Z500grids(shading),(b)the blocking frequency in the NH in the December–January–February period,where the percentage of SATgrids inside the box showing a linear correlation coefficient larger than 0.574 in magnitude with each of the Z500grids is denoted by the white contours in (a),and that with blocking frequency is shown in(c).(see text for description)

    ·Second,to identify the regions accounting for a significant fraction of the SAT variability in the EAWM region, we set a threshold of 0.574 for r1(in magnitude;corresponding to 33.3%or one-third of the total variance).The thresh-old(0.574)is much greater than commonly used confidence levels(95%or 99%)with 65 degrees of freedom(where the length of study period is 66 years).The output in the analysis is N1(λ,φ,i)(1 if exceeding the threshold and 0 otherwise).

    ·Third,at each Z500grid point(λ,φ)we count the percentage of SAT grids in the EAWM region with r1exceeding the threshold,where N1(λ,φ)=(λ,φ,i)]/(221× 100)(%;white contours in Fig.2a).

    As shown in Fig.2a,the year-to-year variability of the SAT over the EAWM region is teleconnected with both upstream and downstream signals over the extratropical region. Upstream of the East Asian continent is a tripole pattern, where a low region is centered over the European continent and the Mediterranean Sea(~15°E),while a high region is centered near the Ural Mountains(~60°E)and another low region can be found over western Siberia(~90°–100°E). This is analogous to the dominant UB pattern(Wang et al., 2010b;Cheung et al.,2012).A significant negative correlation is also found over the Arabian Sea(~20°N,60°E),but further analysis reveals that this signal is significantly correlated with UB(figure not shown).On the other hand,downstream of the East Asian continent is a dipole pattern overthe western North Pacific.This resembles the WP teleconnection pattern(Wallace and Gutzler,1981).The above results suggest that UB and WP account for a significant fraction of the SAT variability in the EAWM region.

    The Z500anomaly pattern resembles an inverted omega pattern(Yeh et al.,1962;Ding,1994),suggesting that lower SAT in the EAWM region is related to the more frequent occurrence of blocking.The area-averaged SAT of the EAWM region is significantly correlated with the blocking frequency over 45°–90°E(the Urals sector)and 105°–130°E(the WPsector;Fig.2b).Note that the blocking frequency is defined in the same way as the algorithm listed in Cheung et al.(2015),which identifies the regions showing a reversal of Z500gradient over the midlatitudes for at least four consecutive days.

    To identify the key blocking region accounting for a significant fraction of the SAT variability in the EAWM region, we follow the three steps that constructed the white contours in Fig.2a.First,we correlate the blocking frequency at each longitude with each of the 221 SAT grids inside the EAWM region.The output of this correlation analysis is r2(λ,i),whereλ∈[0,357.5]°E and i∈[1,221].Second, we use N2(λ,i)to record whether the correlation exceeds the threshold(0.574 in magnitude),where N2(λ,i)=1 if exceeding this threshold and N2(λ,i)=0 otherwise.Third, for each longitude,we count the percentage of SAT grids in the EAWM region with r2exceeding the threshold,where

    As shown in Fig.2c,more than 30%of the SAT grids are found to show a linear correlation coefficient larger than 0.574 in magnitude with the UB sector,but much less than 10%of the SAT grids have such a strong correlation with the WP blocking.Therefore,the Urals sector is the major blocking sector strongly linked to SAT variability in the EAWM region.

    Accordingly,we deduce that a blocking pattern centered near the Urals and a north–south-oriented dipole over the western NorthPacific are two large-scaleatmosphericsignals related to SAT variability in the EAWM region:

    ·Signal(1)can be represented by the UBI.This is defined as the first leading EOF pattern obtained from the covariance matrix of Z500enclosing the region(30°–80°N, 0°–120°E)(Fig.3a).Each grid of the covariance matrix is weighed by cosφ.The definition is similar to that in Cheung et al.(2012),but a larger domain is chosen here in order to capture the tripole pattern in Fig.2a.Unlike in Cheung et al. (2015),the area-averaged blocking frequency is not adopted here because the EOF pattern can be easily reproduced,and this EOF pattern is closely related to the UB frequency.

    ·Signal(2)is captured well by the western Pacific index(WPI).This is defined as the first leading EOF obtained from the correlation matrix of the Z500anomalies enclosing the regionof(20°–80°N,90°E–120°W),whichresemblesthe WP pattern(Wallace and Gutzler,1981;Cheunget al.,2015). Because the positive sign of the WPI represents a lower Z500over the high-latitude region of the North Pacific,its sign is reversed in Fig.3b.

    3.Distinct relationships with EAWM temperatures

    Fig.3.Eigenvector of the(a)UB and(b)WP teleconnection pattern.Linear correlation coefficients between SAT and the principal component of(c)UB and(d)WP.In(c,d),only the regions significant at the 95%confidence level are shaded,and the values larger than 0.574 in magnitude(equivalent to one-third of the explained variance)are enclosed by thick contours.

    After defining the two teleconnection patterns in Figs.3a and b,their standardized principal component is correlated with the SAT field in order to illustrate their relationship with the EAWM temperatures.As shown in Figs.3c and d,thetwo patterns are associated with two distinct temperaturepatterns in Asia.The UBI is significantly negatively correlated with the SAT over the midlatitude region in Asia(including Siberia,central and northern parts of China,Korea,and Japan),with the strongest negative correlation southeast of Lake Baikal(Fig.3c).The WPI,on the other hand,is accompanied by a temperature dipole anomaly pattern over the Asia–Pacific region.In particular,the negative phase of the WPI is significantly negatively correlated with the SAT over the coastal region in East Asia and Southeast Asia,including southeastern China,South Korea,southern Japan,and Vietnam(Fig.3d).In short,the UBI forms a very strong linkage with the SAT north of 40°N,whereas the WPI is strongly linked to the SAT south of 40°N in the EAWM region.

    During the observational period,the interannual and interdecadal variations of the SAT in the EAWM region can be described mainly by two dominant modes(Wang et al., 2010a).Applying an EOF analysis to the covariance matrix of SAT in the EAWM region(0°–60°N,100°–140°E),Wang et al.(2010a)called the first two EOFs the NTM and STM. The spatial patterns of the two temperature modes are given in Figs.4a and b.Following their definitions,the two temperature modes were obtained using the ERA-40 datasets. Then,the principal component of the two SAT modes was obtained by projecting the two eigenvectors onto the NCEP data,where the climatological mean difference between theNCEP and ERA-40 data was subtracted before the projection.The two modes account for 72%of the total variance, where the NTM(STM)explains well the SAT variability north(south)of 40°N.The region explained by the UBI (WPI)seems to coincide with that of the NTM(STM).

    Fig.4.(a,b)Spatial patterns of the first two leading EOF modes(EOF1 and EOF2)of SAT over the EAWM region(0°–60°N,100°–140°E)using the ERA-40 dataset for the period 1957/58–2000/01:(a)NTM;(b)STM (units:°C)[redrawn from Wang et al.(2010a,Figs.3a and c)].(c,d)Standardized principal component of the NTM/STM(black line)and the time series of the UBI(blue line)/-WPI(red line).The linear correlation coefficient between the two time series is shown in the top-right.

    The UBI–NTM and WPI–STM relationshipsare revealed in Figs.4c and d,where their year-to-year correlations are larger than 0.7 in magnitude and the explained variances are greater than 50%.Because the intensity of the EAWM has been found to have undergone a strong decadal variation in the late 1980s(e.g.,Jhun and Lee,2004;Wang et al.,2009b; Wang et al.,2010a),one might wonder if these relationships are still strong on interannual timescales.In other words, the relationship in Figs.4c and d might be mainly due to the decadal variation.In this regard,we applied an 8-year highpass(low-pass)Lanczos filter with 21 weights to extract the interannual(interdecadal)component of each time series in Figs.4c and d(Fig.5).

    As can be seen in Figs.5a and b,both the UBI–NTM and WPI–STM relationships still have a linear correlation coefficient of larger than 0.7 in magnitude on interannual timescales.On the other hand,NTM and STM undergo a decadal variation around 1986/87,where they change from a predominantly positive to negative phase(Fig.5c–d).Correspondingly,the SAT in the EAWM region tends to be below normal before this period and above normal afterward(Lee et al.,2013).In the late 2000s,the two temperature modes become less negative due to some severe winters,which is consistent with the strengthening tendency of the EAWM after 2004/05(Wang and Chen,2014a).In comparison,both the UBI andWPI can capturethe decadalchangein the1980s (dashed colored lines in Figs.4c and d).The correlation between the interdecadal time series of the UBI and NTM is 0.887,and that between the WPI and STM is 0.883.Apparently,the UB and WP are two large-scale teleconnection patterns that can account for a significant fraction of SAT variability in different parts of the EAWM region on interannual and interdecadal timescales.

    4.Representativeness of the UBI and WPI

    Fig.5.As in Figs.4c and d but for(a,b)the interannual component(<8 years)and(c,d)the interdecadal component(>8 years)of the time series.

    To measure how well the UBI and WPI can capture SAT variability in the EAWM region compared to other EAWMIsand teleconnection patterns[Arctic Oscillation(AO),ENSO Ni?o3 index;Eurasian pattern(EU),Scandinavian pattern (SCAN)],all indices are correlated with the NTM and STM in Figs.6a–d.The definition of all EAWMIs is given in Table 1,whereas the definitions of the teleconnection patterns are the same as those of the CPC/NOAA(http://www.cpc. ncep.noaa.gov/data/teledoc/telecontents.shtml).Because all climate indices on the CPC/NOAA website are available for 1950 onwards,the correlation analysis is confined to the period 1950/51–2013/14.

    Among the 19 EAWMIs listed in Table 1,only five show a correlation of 0.50 or above in magnitude with the NTM in the entire study period(Fig.6a).Two of them(WC14 and H15)are new indices,following the review of Wang and Chen(2010).They consider the meridional MSLP or Z500gradient over East Asia instead of solely the zonal gradient. The remaining three EAWMIs describe either the intensity of the Siberian high(G01)or the East Asian trough(S97 and CS99b).In contrast,only three EAWMIs have a correlation smaller than 0.50 in magnitude with the STM(Fig.6b).These results are consistent with those of previous studies(e.g.,Wang et al.,2010a;Wang and Chen,2014a),where most of the EAWMIswere foundto showa strongcorrelation with the STM but not the NTM.

    Fig.6.Correlation between EAWMIs and the(a)NTM and(b)STM,where the first column of each plot shows the correlation throughout the entire study period,and the second column shows the 31-year running correlation.(c,d)As in(a,b)but for the correlation with teleconnection indices.

    Table 1.List of EAWMIs.The EAWMIs considered in Wang and Chen’s(2010)study are classified into four groups(EAWMI1 to EAWMI4)based on the parameter defining the index.The two remaining EAWMIs(WC14 and H15)have been defined recently.

    As mentionedin Wanget al.(2010a),a coldNTMis characterized by a northwestward shift in the East Asian trough toward Lake Baikal(~50°N,110°E)and a northwestward shift in the Siberian high[which is located at(40°–65°N, 80°–120°E)in the climatology].A cold STM,on the other hand,is related to a deepening of the East Asian trough and a stronger surface MSLP gradient between the Siberian high and the subtropical western North Pacific.Most of the EAWMIs listed in Table 1 belongto groups1 and 2(13 outof 19),wheregroup1 is definedbythe zonalMSLP gradientbetween the Asian continentand the Pacific Ocean,and group2 is defined by the lower troposphericV over the coastal region in East Asia.These EAWMIs describe mainly the large-scale circulation features south of 40°N.Accordingly,they tend to show a strong correlation with the STM that represents the SAT variability south of 40°N(Fig.4b).

    ComparedtotheEAWMIsandotherwell-knownteleconnectionindices,the UBI has thestrongestcorrelationwiththe NTM,whereas the WPI has a correlation equally as strong (~0.7)as that of four of the EAWMIs with the STM.Such strong year-to-year relationships are stable throughout the study period,as evidenced by slight changes in the 31-year running correlation.The above conclusion is still valid if we repeat the correlation analysis for the interannual and interdecadal component of the indices(same as Fig.5),except that more EAWMIs show a correlation with the STM greaterthan0.7ontheinterannualtimescales,andH15shows a stronger correlation with the NTM(0.913)than the UBI (0.888)on the interdecadal timescales(figures not shown). Therefore,both the UBI and WPI are representative for measuring EAWM temperatures.Because many EAWMIs can represent the STM well,the remaining part of this section focuses on the NTM.

    The UBI–NTM relationship is related to the linkage be-tween UB and both the Siberian high and the East Asian trough.As can be inferred from Fig.3a,the positive phase of the UBI corresponds to stronger northerly geostrophic wind from the polar region towards western Siberia.This enhances the cold air intensity(i.e.,negative SAT anomaly) over the climatological Siberian high region.Moreover,the strongermeridional-typecirculationduringthepositivephase of the UBI tends to intensify the cyclonic flow over the East Asian continent,which can be deduced by the negative Z500anomalyoverEastAsia(Fig.3a).ThisnegativeZ500anomaly pattern is analogous to the northwestward shift of the East Asian trough from the south of Japan towards the East Asian continent.

    Compared to the UBI,SCAN shows a weaker correlation with the NTM throughout the study period.However,the 31-year moving correlation suggests that the SCAN–NTM relationship has become stronger in recent decades,which is comparable to the UBI–NTM relationship.Previous studies have also identified its impact on the EAWM(Bueh and Nakamura,2007;Bueh et al.,2011;Sohn et al.,2011;Liu et al.,2014).The difference between SCAN and the UBI arises from their centers of action over Eurasia;those of SCAN are locatedat a higherlatitude.The aboveresults suggest that the large-scale circulations associated with the NTM are characterizedbya quasi-stationaryRossby wave-trainpatternrecurring over Eurasia.This is also related to the location of UB. In future work,we will look to further explain the dynamic mechanisms responsible for the recurring position of blocking and the quasi-stationary Rossby wave train over Eurasia.

    5.Summary and discussion

    Based on recent advances in EAWM studies(e.g.,Lee and Jhun,2006;Wang et al.,2010a;Takaya and Nakamura, 2013;Kim et al.,2014;Liu et al.,2014;Hu et al.,2015;Leung and Zhou,2015)and the results of this study,an EAWMI should be able to represent the following anomalous largescale circulation features in order to capture the SAT variability in the majority of the EAWM region:

    ·ThestrengthoftheSiberianhigh,whichis thecold-core system related to cold air activity over the EAWM region;

    ·The meridional-type circulation(or zonal index)over East Asia,which is related to the meridional displacement of the East Asian troughand the cold-air pathwayin the EAWM region.

    The former feature is captured well by any index enclosing the climatological Siberian high region(e.g.,Hu et al., 2015).This is closely related to the teleconnection pattern centered near the Ural Mountains,which can be described by the UBI[Urals blocking pattern(Cheung et al.,2012)].The latter feature is closely related to the deepening of the East Asian trough(Wang et al.,2009a;Leung and Zhou,2015). This resembles the large-scale circulation features associated with the WPI(Lim and Kim,2013;Takaya and Nakamura, 2013;Wang and Chen,2014b).Therefore,variability in the EAWM can be explained well by the large-scale circulation features upstream and downstream of the EAWM region.

    However,the SAT variability north and south of 40°N in the EAWM region(NTM and STM)is quite different on interannual timescales.As a result,very few EAWMIs show an equally strong correlation with both the NTM and STM. Among the EAWMIs considered in this study,only CS99b(a trough index)shows a correlation above 0.60 with both the NTM and STM.We have also tried to use a linear combination of the UBI and WPI to represent the NTM and STM,and the optimization is 3×UBI-2×WPI.The correlation between this linear combination and the NTM(STM)is 0.614 (0.626),which is close to that of CS99b.In other words,a unified and representative EAWMI(say,with an explained variance larger than 50%in the entire EAWM region)is not easy to establish.One possible way is to consider the latitude and longitude dependence of SAT on the two important EAWM circulation features mentioned above(i.e.,EAWMI is a function of latitude and longitude).

    To improve the seasonal forecasting of temperature anomalies in the EAWM and to project future changes in the EAWM,more in-depth analyses should be carried out to explore the underlying dynamics regarding the large-scale teleconnections related to the EAWM,such as the UB and WP. The WP is related to ENSO,which has been discussed extensively in previous studies(e.g.,Li,1990;Wang et al.,2000; Zhou et al.,2007;Wang et al.,2008;Takaya and Nakamura, 2013;Chen et al.,2014).On the other hand,UB has been studied less extensively.

    On synoptic timescales,UB is an important dynamic precursor of severe cold-air outbreaks in the EAWM.The occurrence of UB strengthens cold-air advection from the polar region to western Siberia.This enhances the Siberian high and subsequent cold-air outbreaks in the EAWM(Tao,1957; Takaya and Nakamura,2005;Lu and Chang,2009;Zhou et al.,2009;Cheunget al.,2013).On seasonaltimescales,when UB occurs frequentlyin a particular winter,the Siberian high intensity tends to be higherthan normal.Because of the close linkage between the Siberian high intensity and the NTM, the temperature tends to be lower in the northern part of the EAWM(Wang et al.,2010b;Chang and Lu,2012;Cheung et al.,2012).Moreover,the occurrence of UB is also characterized by a stronger meridional-type flow over East Asia. The more frequent occurrence of UB partly enhances the southward intrusion of cold air into the southern part of the EAWM.This potentially causes more frequent cold extremes in theEAWM,as in HongKong(Cheunget al.,2015).Therefore,the dynamics of UB should be investigated in depth in order to improve the understanding of the EAWM on different timescales.In addition to the major dynamic factors that lead to the occurrence of UB,we should understand the factors related to the recurring location of UB.

    Acknowledgements.The authors greatly appreciate the valuable comments provided by the two anonymous reviewers,which helped improve the clarity of the manuscript.The work described in this paper was supported by Shenzhen Research Project(Grant No.GJHS20120820144245169),and the French/Hong Kong Joint Research Project(No.F-HK002/12T).

    REFERENCES

    Bueh,C.,and H.Nakamura,2007:Scandinavian pattern and its climatic impact.Quart.J.Roy.Meteor.Soc.,133,2117–2131, doi:10.1002/qj.173.

    Bueh,C.,N.Shi,and Z.W.Xie,2011:Large-scale circulation anomalies associated with persistent low temperature over southern China in January 2008.Atmos.Sci.Lett.,12,273–280,doi:10.1002/asl.333.

    Chan,J.C.L.,and C.Y.Li,2004:The EastAsian winter monsoon. East Asian Monsoon,C.-P.Chang,Ed.,World Scientific Publishing Co.Pet.Ltd.,54–106.

    Chang,C.-P.,and K.M.W.Lau,1980:Northeasterly cold surges and near-equatorial disturbances over the winter MONEX area during December 1974.Part II:Planetary-scale aspects. Mon.Wea.Rev.,108,298–312.

    Chang,C.-P.,and M.-M.Lu,2012:Intraseasonal predictability of Siberian high and East Asian winter monsoon and its interdecadal variability.J.Climate,25,1773–1778,doi:10.1175/ JCLI-D-11-00500.1.

    Chen,J.,and S.Q.Sun,1999:Eastern Asian winter monsoon anomaly and variation of global circulation Part I:A comparison study on strong and weak winter monsoon.Chinese J.Atmos.Sci.,23,101–111.(in Chinese).

    Chen,W.,H.-F.Graf,and R.H.Huang,2000:The interannual variability of East Asian winter monsoon and its relation to the summer monsoon.Adv.Atmos.Sci.,17,48–60,doi: 10.1007/s00376-000-0042-5.

    Chen,Z.,R.G.Wu,and W.Chen,2014:Distinguishing interannual variations of the northern and southern modes of the east Asian winter monsoon.J.Climate,27,835–851,doi: 10.1175/JCLI-D-13-00314.1.

    Cheung,H.H.N.,and W.Zhou,2015:Implications of Ural blocking for East Asian winter climate in CMIP5 GCMs.Part I: Biases in the historical scenario.J.Climate,28,2203–2216, doi:10.1175/JCLI-D-14-00308.1.

    Cheung,H.N.,W.Zhou,H.Y.Mok,and M.C.Wu,2012:Relationship between Ural–Siberian blocking and the East Asian winter monsoon in relation to the Arctic oscillation and the El Ni?o–southern oscillation.J.Climate,25,4242–4257,doi: 10.1175/JCLI-D-11-00225.1.

    Cheung,H.N.,W.Zhou,Y.P.Shao,W.Chen,H.Y.Mok,and M.C.Wu,2013:Observational climatology and characteristics of wintertime atmospheric blocking over Ural–Siberia. Climate Dyn.,41,63–79.

    Cheung,H.H.N.,W.Zhou,S.M.Lee,and H.W.Tong,2015:Interannual and interdecadal variability of the number of cold days in Hong Kong and their relationship with large-scale circulation.Mon.Wea.Rev.,143,1438–1454,doi:10.1175/ MWR-D-14-00335.1.

    Cui,X.P.,and Z.B.Sun,1999:East Asian winter monsoon index and its variation analysis.Journal of Nanjing Institute of Meteorology,22,321–325.(in Chinese).

    Ding,Y.H.,1994:Monsoon over China.Kluwer Academic Publishers,420 pp.

    Gong,D.Y.,S.W.Wang,and J.H.Zhu,2001:East Asian winter monsoon and Arctic oscillation.Geophys.Res.Lett.,28, 2073–2076.

    Guo,Q.W.,1994:Relationship between the variations of East Asian winter monsoon and temperature anomalies in China. Quarterly Journal of Applied Meteorology,5,218–225.(in Chinese).

    Hsu,H.-H.,and J.M.Wallace,1985:Vertical structure of wintertime teleconnection patterns.J.Atmos.Sci.,42,1693–1710.

    Hu,C.D.,S.Yang,and Q.G.Wu,2015:An optimal index formeasuring the effect of East Asian winter monsoon on China winter temperature.Climate Dyn.,45,2751–2589,doi:10.1007/ s00382-015-2493-5.

    Hu,Z.-Z.,L.Bengtsson,and K.Arpe,2000:Impact of global warming on the Asian winter monsoon in a coupled GCM. J.Geophys.Res.,105,4607–4624,doi:10.1029/1999JD 901031.

    Jhun,J.G.,and E.J.Lee,2004:A new East Asian winter monsoon index and associated characteristics of the winter monsoon. J.Climate,17,711–726,doi:10.1175/1520-0442(2004)017<0711:ANEAWM>2.0.CO;2.

    Ji,L.R.,S.Q.Sun,K.Arpe,and L.Bengtsson,1997:Model study on the interannual variability of Asian winter monsoon and its influence.Adv.Atmos.Sci.,14,1–22,doi:10.1007/s00376-997-0039-4.

    Joung,C.H.,and M.H.Hitchman,1982:On the role of successive downstream development in East Asian polar air outbreaks. Mon.Wea.Rev.,110,1224–1237.

    Kim,Y.,K.Y.Kim,and S.Park,2014:Seasonal scale variability of the East Asian winter monsoon and the development of a two-dimensional monsoon index.Climate Dyn.,42,2159–2172,doi:10.1007/s00382-013-1724-x.

    Lau,N.-C.,and K.-M.Lau,1984:The structure and energetics of midlatitude disturbances accompanying cold-air outbreaks over East Asia.Mon.Wea.Rev.,112,1309–1327,doi: 10.1175/1520-0493(1984)112<1309:TSAEOM>2.0.CO;2.

    Lee,H.S.,and J.G.Jhun,2006:Two types of the Asian continental blocking and their relation to the east Asian monsoon during the boreal winter.Geophys.Res.Lett.,33,L22707,doi: 10.1029/2006GL027948.

    Lee,S.-S.,S.-H.Kim,J.-G.Jhun,K.-J.Ha,and Y.-W.Seo,2013: Robust warming over East Asia during the boreal winter monsoon and its possible causes.Environ.Res.Lett.,8,034001, doi:10.1088/1748-9326/8/3/034001.

    Leung,M.Y.-T.,and W.Zhou,2015:Variation of circulation and East Asian climate associated with anomalous strength and displacement of the East Asian trough.Climate Dyn.,45, 2731–2732,doi:10.1007/s00382-015-2504-6.

    Li,C.Y.,1990:Interaction between anomalous winter monsoon in East Asia and El Ni?o events.Adv.Atmos.Sci.,7,36–46,doi: 10.1007/BF02919166.

    Li,Y.Q.,and S.Yang,2010:A dynamical index for the East Asian winter monsoon.J.Climate,23,4255–4262,doi:10.1175/ 2010JCLI3375.1.

    Lim,Y.K.,and H.D.Kim,2013:Impact of the dominant largescale teleconnections on winter temperature variability over East Asia.J.Geophys.Res.,118,7835–7848,doi:10.1002/ jgrd.50462.

    Lim,Y.-K.,and H.-D.Kim,2015:Comparison of the impact of the Arctic Oscillation and Eurasian teleconnection on interannual variation in East Asian winter temperatures and monsoon.Theor.Appl.Climatol.,1–13,doi:10.1007/s00704-015-1418-x.

    Liu,Y.Y.,L.Wang,W.Zhou,and W.Chen,2014:Three Eurasian teleconnection patterns:Spatial structures,temporal variability,and associated winter climate anomalies.Climate Dyn., 42,2817–2839,doi:10.1007/s00382-014-2163-z.

    Lu,E.,and J.C.L.Chan,1999:A unified monsoon index for South China.J.Climate,12,2375–2385,doi:10.1175/1520-0442(1999)012<2375:AUMIFS>2.0.CO;2.

    Lu,M.M.,and C.-P.Chang,2009:Unusual late-season cold surges during the 2005 Asian winter monsoon:Roles of Atlantic blocking and the central Asian anticyclone.J.Climate, 22,5205–5217,doi:10.1175/2009JCLI2935.1.

    Shi,N.,1996:Features of the East Asian winter monsoon intensity on multiple time scale in recent 40 years and their relation to climate.Journal of Applied Meteorological Science,7,175–182.(in Chinese).

    Sohn,S.-J.,C.-Y.Tam,and C.-K.Park,2011:Leading modes of East Asian winter climate variability and their predictability: An assessment of the APCC multi-model ensemble.J.Meteor.Soc.Japan,89,455–474,doi:10.2151/jmsj.2011-504.

    Song,J.,C.Y.Li,and W.Zhou,2014:High and low latitude types of the downstream influences of the North Atlantic Oscillation.Climate Dyn.,42,1097–1111,doi:10.1007/s00382-013-1844-3.

    Sun,B.-M.,and C.-Y.Li,1997:Relationship between the disturbances of East Asian trough and tropical convective activity in boreal winter.Chinese Science Bulletin,42,500–504.(in Chinese).

    Takaya,K.,and H.Nakamura,2005:Mechanisms of intraseasonal amplification of the cold Siberian high.J.Atmos.Sci.,62, 4423–4440,doi:10.1175/JAS3629.1.

    Takaya,K.,and H.Nakamura,2013:Interannual variability of the East Asian winter monsoon and related modulations of the planetary waves.J.Climate,26,9445–9461,doi:10.1175/ JCLI-D-12-00842.1.

    Tao,S.Y.,1957:A study of activities of cold airs in East Asian winter.Handbook of Short-Term Forecast,China Meteorological Administration,Ed.,Meteorology Press,60–92.

    Wallace,J.M.,and D.S.Gutzler,1981:Teleconnections in the geopotential height field during the northern hemisphere winter.Mon.Wea.Rev.,109,784–812.

    Wang,B.,R.G.Wu,and X.H.Fu,2000:Pacific–East Asian teleconnection:how does ENSO affect East Asian climate?J. Climate,13,1517–1536,doi:10.1175/1520-0442(2000)013<1517:PEATHD>2.0.CO;2.

    Wang,B.,Z.W.Wu,C.-P.Chang,J.Liu,J.P.Li,and T.J.Zhou, 2010a:Another look at interannual-to-interdecadal variations of the east Asian winter monsoon:The northern and southern temperature modes.J.Climate,23,1495–1512,doi: 10.1175/2009JCLI3243.1.

    Wang,L.,and W.Chen,2010:How well do existing indices measure the strength of the East Asian winter monsoon?Adv. Atmos.Sci.,27,855–870,doi:10.1007/s00376-009-9094-3.

    Wang,L.,and W.Chen,2014a:The East Asian winter monsoon: re-amplification in the mid-2000s.Chinese Science Bulletin, 59,430–436,doi:10.1007/s11434-013-0029-0.

    Wang,L.,and W.Chen,2014b:An intensity index for the East Asian winter monsoon.J.Climate,27,2361–2374,doi: 10.1175/JCLI-D-13-00086.1.

    Wang,L.,W.Chen,and R.H.Huang,2008:Interdecadal modulation of PDO on the impact of ENSO on the east Asian winter monsoon.Geophys.Res.Lett.,35,doi:10.1029/2008GL 035287.

    Wang,L.,W.Chen,W.Zhou,and R.H.Huang,2009a:Interannual variations of East Asian trough axis at 500 hPa and its association with the East Asian winter monsoon pathway.J. Climate,22,600–614,doi:10.1175/2008JCLI2295.1.

    Wang,L.,R.H.Huang,L.Gu,W.Chen,and L.H.Kang,2009b: Interdecadal variations of the East Asian winter monsoon and their association with quasi-stationary planetary wave activity.J.Climate,22,4860–4872,doi:10.1175/2009JCLI 2973.1.

    Wang,L.,W.Chen,W.Zhou,J.C.L.Chan,D.Barriopedro,and R.H.Huang,2010b:Effect of the climate shift around mid 1970s on the relationship between wintertime Ural blocking circulation and East Asian climate.Int.J.Climatol.,30,153–158,doi:10.1002/joc.1876.

    Wu,B.Y.,and J.Wang,2002:Winter Arctic oscillation,Siberian high and East Asian winter monsoon.Geophys.Res.Lett.,29, 1897,doi:10.1029/2002GL015373.

    Xu,S.Y.,and J.J.Ji,1965:The climate and weather features during the outbreak period of China’s winter monsoon.Geographical Symposium,9,85–101.(in Chinese).

    Yang,S.,K.-M.Lau,and K.-M.Kim,2002:Variations of the East Asian jet stream and Asian-Pacific-American winter climate anomalies.J.Climate,15,306–325,doi:10.1175/1520-0442 (2002)015<0306:VOTEAJ>2.0.CO;2.

    Yeh,T.-C.,S.Y.Tao,and B.Z.Zhu,1962:Studies on the Blocking Situation in the Northern Hemisphere in Winter.Science Press,130 pp.(in Chinese).

    Zhou,W.,X.Wang,T.J.Zhou,C.Li,and J.C.L.Chan,2007: Interdecadal variability of the relationship between the East Asian winter monsoon and ENSO.Meteor.Atmos.Phys.,98, 283–293,doi:10.1007/s00703-007-0263-6.

    Zhou,W.,J.C.L.Chan,W.Chen,J.Ling,J.G.Pinto,and Y.P. Shao,2009:Synoptic-scale controls of persistent low temperature and icy weather over southern China in January 2008. Mon.Wea.Rev.,137,3978–3991,doi:10.1175/2009MWR 2952.1.

    Cheung,H.H.N.,and W.Zhou,2016:Simple metrics for representing East Asian winter monsoon variability: Urals blocking and western Pacific teleconnection patterns.Adv.Atmos.Sci.,33(6),695–705,

    10.1007/s00376-015-5204-6.

    18 September 2015;revised 11 November 2015;accepted 30 November 2015)

    Wen ZHOU

    Email:wenzhou@cityu.edu.hk

    日韩电影二区| 黑人猛操日本美女一级片| 男的添女的下面高潮视频| 人妻人人澡人人爽人人| 女性被躁到高潮视频| 日本一区二区免费在线视频| bbb黄色大片| 99精品久久久久人妻精品| 亚洲欧洲国产日韩| 精品午夜福利在线看| 七月丁香在线播放| 精品第一国产精品| 国产亚洲最大av| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 欧美日韩一级在线毛片| 黄色视频不卡| 伊人久久国产一区二区| 热re99久久国产66热| 大码成人一级视频| 亚洲激情五月婷婷啪啪| 99re6热这里在线精品视频| 999久久久国产精品视频| 日韩欧美精品免费久久| √禁漫天堂资源中文www| 久久精品亚洲av国产电影网| 美女主播在线视频| 一区在线观看完整版| 国产av精品麻豆| 欧美在线黄色| 99热全是精品| 高清不卡的av网站| 少妇的丰满在线观看| 亚洲国产欧美一区二区综合| 日韩熟女老妇一区二区性免费视频| 黑人巨大精品欧美一区二区蜜桃| 精品卡一卡二卡四卡免费| 精品午夜福利在线看| 免费人妻精品一区二区三区视频| www.av在线官网国产| 丝袜喷水一区| 侵犯人妻中文字幕一二三四区| 欧美少妇被猛烈插入视频| 久久精品亚洲av国产电影网| 激情视频va一区二区三区| 久久综合国产亚洲精品| 精品国产一区二区久久| 国产伦人伦偷精品视频| 赤兔流量卡办理| 国产欧美日韩一区二区三区在线| 欧美日韩国产mv在线观看视频| 一级毛片 在线播放| 欧美在线一区亚洲| 黄片小视频在线播放| 亚洲伊人色综图| 男男h啪啪无遮挡| 18禁裸乳无遮挡动漫免费视频| 久久精品亚洲av国产电影网| 亚洲精品乱久久久久久| 国产亚洲午夜精品一区二区久久| 熟女少妇亚洲综合色aaa.| av网站免费在线观看视频| 国产又爽黄色视频| 欧美成人午夜精品| 久久99热这里只频精品6学生| 亚洲熟女精品中文字幕| 美女中出高潮动态图| 国产av国产精品国产| 日韩伦理黄色片| 色94色欧美一区二区| 丝瓜视频免费看黄片| 91精品伊人久久大香线蕉| 爱豆传媒免费全集在线观看| 亚洲精品aⅴ在线观看| 男女无遮挡免费网站观看| 亚洲精品成人av观看孕妇| 色播在线永久视频| 一边亲一边摸免费视频| 欧美日韩视频精品一区| 国产有黄有色有爽视频| 国产不卡av网站在线观看| 韩国精品一区二区三区| 少妇人妻精品综合一区二区| 美国免费a级毛片| 亚洲国产最新在线播放| 国产亚洲最大av| 国产麻豆69| 交换朋友夫妻互换小说| 亚洲成av片中文字幕在线观看| 亚洲人成网站在线观看播放| 天天躁夜夜躁狠狠躁躁| 国产精品熟女久久久久浪| 久久久久人妻精品一区果冻| 久久久久精品久久久久真实原创| 精品国产露脸久久av麻豆| 秋霞伦理黄片| 伊人久久大香线蕉亚洲五| 国产黄色视频一区二区在线观看| 啦啦啦视频在线资源免费观看| 免费av中文字幕在线| 日韩一卡2卡3卡4卡2021年| 哪个播放器可以免费观看大片| 欧美在线一区亚洲| 亚洲人成77777在线视频| 国产又色又爽无遮挡免| 中文字幕亚洲精品专区| 热99国产精品久久久久久7| 黄片播放在线免费| 日韩av免费高清视频| 亚洲欧美成人综合另类久久久| 黄频高清免费视频| 91精品三级在线观看| 秋霞伦理黄片| 国产成人系列免费观看| 国产日韩欧美亚洲二区| 九九爱精品视频在线观看| 日韩一区二区三区影片| av不卡在线播放| 一区二区三区四区激情视频| 亚洲欧洲日产国产| 日韩一卡2卡3卡4卡2021年| 黄片无遮挡物在线观看| 午夜福利网站1000一区二区三区| 亚洲av日韩精品久久久久久密 | 日韩人妻精品一区2区三区| 亚洲人成77777在线视频| 国产av一区二区精品久久| 欧美日韩亚洲高清精品| 99久国产av精品国产电影| 久久久国产一区二区| 精品一区二区三区四区五区乱码 | 久久久久久久国产电影| 一二三四在线观看免费中文在| netflix在线观看网站| 国产一区二区 视频在线| 亚洲综合色网址| 精品一区二区免费观看| av免费观看日本| 超色免费av| 亚洲精品aⅴ在线观看| 免费女性裸体啪啪无遮挡网站| 国产熟女欧美一区二区| 国产av精品麻豆| 成人亚洲精品一区在线观看| 亚洲在久久综合| 国产又色又爽无遮挡免| tube8黄色片| 欧美成人精品欧美一级黄| 天堂俺去俺来也www色官网| av不卡在线播放| 大片免费播放器 马上看| 啦啦啦中文免费视频观看日本| 少妇人妻 视频| 99精品久久久久人妻精品| 国产国语露脸激情在线看| 国产黄频视频在线观看| 免费人妻精品一区二区三区视频| 视频在线观看一区二区三区| 国产女主播在线喷水免费视频网站| 51午夜福利影视在线观看| 免费女性裸体啪啪无遮挡网站| 51午夜福利影视在线观看| 男女之事视频高清在线观看 | 一级黄片播放器| 狂野欧美激情性bbbbbb| 少妇人妻精品综合一区二区| 麻豆精品久久久久久蜜桃| 韩国高清视频一区二区三区| 久久久久久久大尺度免费视频| 一区二区日韩欧美中文字幕| 国产精品.久久久| 一级片免费观看大全| 午夜日本视频在线| 国产黄频视频在线观看| 成人手机av| 久久99精品国语久久久| 久久久亚洲精品成人影院| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品美女久久av网站| 亚洲人成网站在线观看播放| 另类亚洲欧美激情| 另类亚洲欧美激情| 午夜福利一区二区在线看| 久久毛片免费看一区二区三区| 一区二区日韩欧美中文字幕| kizo精华| 不卡视频在线观看欧美| 岛国毛片在线播放| 久久久久久久久久久免费av| 国产成人精品福利久久| 亚洲精品美女久久av网站| 亚洲男人天堂网一区| 老汉色av国产亚洲站长工具| 亚洲,欧美精品.| 亚洲在久久综合| 五月开心婷婷网| 欧美日韩亚洲综合一区二区三区_| 国产福利在线免费观看视频| 母亲3免费完整高清在线观看| 欧美在线一区亚洲| 老鸭窝网址在线观看| 亚洲美女黄色视频免费看| 男人操女人黄网站| 国产成人91sexporn| 一本色道久久久久久精品综合| 国产亚洲欧美精品永久| 欧美国产精品一级二级三级| 国产高清国产精品国产三级| 男女国产视频网站| 在线观看国产h片| 国产视频首页在线观看| 国产视频首页在线观看| 国精品久久久久久国模美| av卡一久久| 亚洲欧美清纯卡通| 亚洲欧美一区二区三区国产| 免费看av在线观看网站| 蜜桃国产av成人99| 亚洲精品国产av蜜桃| 乱人伦中国视频| 欧美精品高潮呻吟av久久| 成人午夜精彩视频在线观看| 18禁国产床啪视频网站| 国产成人精品福利久久| 精品免费久久久久久久清纯 | 夫妻午夜视频| www.熟女人妻精品国产| 久久久久人妻精品一区果冻| 精品一区二区三区av网在线观看 | 久久久国产精品麻豆| 国产免费福利视频在线观看| 这个男人来自地球电影免费观看 | 我要看黄色一级片免费的| 老司机深夜福利视频在线观看 | 久久国产精品大桥未久av| 午夜福利网站1000一区二区三区| 国产精品无大码| 精品久久久久久电影网| 国产亚洲欧美精品永久| 欧美精品av麻豆av| 校园人妻丝袜中文字幕| 看免费av毛片| 免费在线观看视频国产中文字幕亚洲 | 精品一区在线观看国产| 熟妇人妻不卡中文字幕| 亚洲欧美成人精品一区二区| 婷婷色麻豆天堂久久| 国产精品免费大片| www日本在线高清视频| 天天影视国产精品| 国产男女超爽视频在线观看| 最新在线观看一区二区三区 | 大片电影免费在线观看免费| 狠狠精品人妻久久久久久综合| 国产精品.久久久| 国产麻豆69| 亚洲婷婷狠狠爱综合网| 涩涩av久久男人的天堂| 看免费成人av毛片| 香蕉国产在线看| 精品一区二区三区四区五区乱码 | 亚洲少妇的诱惑av| 欧美黄色片欧美黄色片| 色94色欧美一区二区| 欧美最新免费一区二区三区| 久久99一区二区三区| 中文欧美无线码| 精品一区二区免费观看| 伦理电影免费视频| 天堂8中文在线网| 亚洲精品国产av成人精品| 日韩中文字幕欧美一区二区 | 中文乱码字字幕精品一区二区三区| 老汉色∧v一级毛片| 国产精品麻豆人妻色哟哟久久| 亚洲成人国产一区在线观看 | h视频一区二区三区| 一本一本久久a久久精品综合妖精| 大陆偷拍与自拍| 国产极品天堂在线| 成人国产av品久久久| 精品国产乱码久久久久久男人| 欧美精品一区二区大全| 男人操女人黄网站| 亚洲av福利一区| 人人澡人人妻人| 亚洲一区中文字幕在线| 熟妇人妻不卡中文字幕| 日韩 欧美 亚洲 中文字幕| 国产精品二区激情视频| 亚洲精品国产区一区二| 精品人妻一区二区三区麻豆| 91精品三级在线观看| 男女免费视频国产| 操出白浆在线播放| 国产成人啪精品午夜网站| 一区二区三区四区激情视频| 国产精品久久久久成人av| 天堂俺去俺来也www色官网| tube8黄色片| 国产探花极品一区二区| 男的添女的下面高潮视频| 少妇人妻久久综合中文| 欧美久久黑人一区二区| 亚洲欧美清纯卡通| 亚洲三区欧美一区| 新久久久久国产一级毛片| 人人妻人人爽人人添夜夜欢视频| 中文字幕人妻丝袜制服| 欧美日韩福利视频一区二区| 搡老乐熟女国产| 日本欧美国产在线视频| 欧美日韩亚洲高清精品| 美女中出高潮动态图| 欧美日本中文国产一区发布| 亚洲国产精品一区二区三区在线| bbb黄色大片| 国产视频首页在线观看| av不卡在线播放| 久久久久久人人人人人| 午夜免费男女啪啪视频观看| 香蕉丝袜av| av免费观看日本| 免费久久久久久久精品成人欧美视频| 黄色毛片三级朝国网站| 天堂俺去俺来也www色官网| 欧美av亚洲av综合av国产av | 日本猛色少妇xxxxx猛交久久| 99久久综合免费| 久久久精品区二区三区| 久久久国产欧美日韩av| 老汉色av国产亚洲站长工具| 亚洲av国产av综合av卡| 各种免费的搞黄视频| 国产1区2区3区精品| 性色av一级| 午夜免费男女啪啪视频观看| 亚洲欧美一区二区三区黑人| 日韩伦理黄色片| 激情五月婷婷亚洲| 2021少妇久久久久久久久久久| 中文字幕亚洲精品专区| 亚洲国产精品999| 一级,二级,三级黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 午夜91福利影院| 悠悠久久av| 欧美变态另类bdsm刘玥| 亚洲七黄色美女视频| 久久这里只有精品19| 最近中文字幕高清免费大全6| 欧美成人午夜精品| 80岁老熟妇乱子伦牲交| 亚洲av日韩在线播放| av片东京热男人的天堂| xxxhd国产人妻xxx| 少妇人妻精品综合一区二区| 久久婷婷青草| 免费观看av网站的网址| 少妇人妻 视频| 啦啦啦视频在线资源免费观看| 妹子高潮喷水视频| 另类亚洲欧美激情| 男女免费视频国产| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 悠悠久久av| 亚洲精品美女久久久久99蜜臀 | av.在线天堂| 一本—道久久a久久精品蜜桃钙片| 一本大道久久a久久精品| tube8黄色片| av.在线天堂| 国产精品亚洲av一区麻豆 | 80岁老熟妇乱子伦牲交| 免费观看人在逋| 中国三级夫妇交换| 国语对白做爰xxxⅹ性视频网站| 伊人久久国产一区二区| 19禁男女啪啪无遮挡网站| 亚洲一码二码三码区别大吗| 女人精品久久久久毛片| 精品卡一卡二卡四卡免费| 老熟女久久久| 精品国产露脸久久av麻豆| 亚洲av男天堂| 99九九在线精品视频| 午夜精品国产一区二区电影| 亚洲久久久国产精品| 亚洲天堂av无毛| 成人漫画全彩无遮挡| 久久久精品区二区三区| 一本一本久久a久久精品综合妖精| 亚洲精品视频女| 欧美国产精品一级二级三级| 免费日韩欧美在线观看| 一区福利在线观看| 欧美精品一区二区大全| 日韩av免费高清视频| 亚洲欧美精品综合一区二区三区| 少妇被粗大的猛进出69影院| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| 女性被躁到高潮视频| 精品一品国产午夜福利视频| 亚洲国产欧美一区二区综合| 久久影院123| 超碰97精品在线观看| 亚洲成av片中文字幕在线观看| 免费高清在线观看视频在线观看| 国产片内射在线| 亚洲精品第二区| av一本久久久久| 超碰97精品在线观看| www.av在线官网国产| 成人国语在线视频| 久久久久网色| 考比视频在线观看| 欧美亚洲 丝袜 人妻 在线| 国产高清国产精品国产三级| 啦啦啦啦在线视频资源| 在线天堂中文资源库| 国产成人免费观看mmmm| 欧美 亚洲 国产 日韩一| 精品国产一区二区三区四区第35| 国产福利在线免费观看视频| 一边亲一边摸免费视频| 麻豆乱淫一区二区| 久久鲁丝午夜福利片| 人人妻人人澡人人看| 亚洲综合精品二区| 免费少妇av软件| 伊人亚洲综合成人网| 久久综合国产亚洲精品| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| 精品午夜福利在线看| 亚洲欧美一区二区三区国产| 成人午夜精彩视频在线观看| 麻豆av在线久日| 精品少妇久久久久久888优播| 99久久人妻综合| 亚洲一码二码三码区别大吗| 亚洲精品国产色婷婷电影| 妹子高潮喷水视频| 看免费成人av毛片| 可以免费在线观看a视频的电影网站 | 无限看片的www在线观看| 国产男人的电影天堂91| 老司机影院成人| 久久天堂一区二区三区四区| 男人操女人黄网站| 啦啦啦在线免费观看视频4| 亚洲美女搞黄在线观看| 男女下面插进去视频免费观看| 中文字幕最新亚洲高清| 免费黄网站久久成人精品| 国产成人精品久久久久久| 丰满乱子伦码专区| 国产精品一二三区在线看| 国产欧美亚洲国产| av又黄又爽大尺度在线免费看| 美女国产高潮福利片在线看| 亚洲,欧美精品.| 精品国产一区二区久久| 女性被躁到高潮视频| 亚洲精华国产精华液的使用体验| av免费观看日本| 五月天丁香电影| 一本一本久久a久久精品综合妖精| 亚洲美女视频黄频| 一边摸一边抽搐一进一出视频| 国产不卡av网站在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 热re99久久国产66热| 秋霞伦理黄片| 日本午夜av视频| 免费日韩欧美在线观看| 欧美黄色片欧美黄色片| 日本wwww免费看| 男人爽女人下面视频在线观看| 国产毛片在线视频| 一本色道久久久久久精品综合| 麻豆精品久久久久久蜜桃| 亚洲人成77777在线视频| 日本猛色少妇xxxxx猛交久久| 精品国产超薄肉色丝袜足j| 欧美日韩综合久久久久久| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区蜜桃| 欧美精品人与动牲交sv欧美| 人人澡人人妻人| videos熟女内射| 国产极品天堂在线| 无限看片的www在线观看| 一本一本久久a久久精品综合妖精| 日本vs欧美在线观看视频| 97精品久久久久久久久久精品| 女人久久www免费人成看片| 国产亚洲欧美精品永久| 亚洲精品国产av蜜桃| 纵有疾风起免费观看全集完整版| 一区二区三区四区激情视频| 日韩 欧美 亚洲 中文字幕| 最黄视频免费看| 色婷婷久久久亚洲欧美| 久久精品国产亚洲av涩爱| 美女扒开内裤让男人捅视频| 免费黄网站久久成人精品| 日韩一本色道免费dvd| 男男h啪啪无遮挡| 久久精品熟女亚洲av麻豆精品| 日本色播在线视频| 国产精品一区二区精品视频观看| 看十八女毛片水多多多| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 国产极品天堂在线| 亚洲综合精品二区| 在线天堂中文资源库| 男女国产视频网站| 国产野战对白在线观看| 久久性视频一级片| 免费在线观看视频国产中文字幕亚洲 | 国产精品.久久久| 丰满乱子伦码专区| 女人高潮潮喷娇喘18禁视频| 只有这里有精品99| 我要看黄色一级片免费的| 伦理电影大哥的女人| 久久久精品94久久精品| 国产亚洲av片在线观看秒播厂| 男人操女人黄网站| 国产精品久久久久久精品电影小说| 大香蕉久久网| 人人澡人人妻人| 精品人妻在线不人妻| 操美女的视频在线观看| 妹子高潮喷水视频| 成人毛片60女人毛片免费| 人人澡人人妻人| 国产成人午夜福利电影在线观看| 成人国语在线视频| 晚上一个人看的免费电影| 久久久亚洲精品成人影院| 一级毛片我不卡| 在线天堂中文资源库| 亚洲国产精品一区三区| 精品国产一区二区三区四区第35| 少妇精品久久久久久久| 成人手机av| 一本久久精品| 亚洲国产日韩一区二区| 综合色丁香网| 国产成人午夜福利电影在线观看| e午夜精品久久久久久久| 悠悠久久av| 十八禁人妻一区二区| 亚洲欧美一区二区三区国产| 欧美激情极品国产一区二区三区| 国产亚洲最大av| 国产免费又黄又爽又色| 亚洲美女视频黄频| 亚洲av综合色区一区| 国产精品国产三级专区第一集| 老司机深夜福利视频在线观看 | 久久久国产一区二区| 久久人人97超碰香蕉20202| 国产精品国产三级专区第一集| 中文字幕精品免费在线观看视频| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 亚洲人成77777在线视频| 亚洲国产最新在线播放| 99国产综合亚洲精品| 国语对白做爰xxxⅹ性视频网站| 大陆偷拍与自拍| 婷婷色综合www| 精品免费久久久久久久清纯 | 国产亚洲午夜精品一区二区久久| 国产又爽黄色视频| 国产成人91sexporn| 熟妇人妻不卡中文字幕| 亚洲熟女精品中文字幕| 欧美成人午夜精品| av网站免费在线观看视频| 又黄又粗又硬又大视频| 精品国产露脸久久av麻豆| 色播在线永久视频| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 国产又色又爽无遮挡免| 成人影院久久| 热99久久久久精品小说推荐| 夫妻午夜视频| 久久狼人影院| 国产激情久久老熟女| 日韩视频在线欧美| 久久久久视频综合| 国产日韩欧美亚洲二区| 精品一区二区免费观看| 新久久久久国产一级毛片| 少妇人妻久久综合中文| av线在线观看网站| 国产国语露脸激情在线看| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧洲国产日韩| 九九爱精品视频在线观看| 在线观看免费午夜福利视频| 国产欧美日韩一区二区三区在线| 国产 一区精品| 亚洲精品国产av蜜桃| 久久久久精品久久久久真实原创| 日韩制服丝袜自拍偷拍| 在线观看免费午夜福利视频|