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

    Asymmetric Relationship between the Meridional Displacement of the Asian Westerly Jet and the Silk Road Pattern

    2018-03-06 03:35:53XiaoweiHONGRiyuLUandShuanglinLIClimateChangeResearchCenterChineseAcademyofSciencesBeijing00029China
    Advances in Atmospheric Sciences 2018年4期

    Xiaowei HONG,Riyu LU,and Shuanglin LIClimate Change Research Center,Chinese Academy of Sciences,Beijing 00029,China

    2State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

    3University of the Chinese Academy of Sciences,Beijing 100049,China

    4Nansen-Zhu International Research Centre,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

    5Department of Atmospheric Science,School of Environmental Studies,China University of Geosciences,Wuhan 430074,China

    1.Introduction

    The meridional displacement of the Asian jet(JMD)in summer,which manifests as the leading mode of zonal wind anomalies in the upper troposphere over the midlatitude Eurasian continent,has been revealed by recent studies(Du et al.,2016;Hong and Lu,2016;Wei et al.,2017).The JMD is characterized by a north–south seesaw pattern,with out-of-phase anomalies in zonal wind to the north of the climatological mean Asian jet axis and to the south.In other words,a northward JMD corresponds to westerly anomalies to the north of the axis and easterly anomalies to the south,and vice versa.Hong and Lu(2016)emphasized that the circulation anomalies associated with the JMD appear inhomogeneous along the zonal direction,with much stronger amplitudes over West Asia and East Asia but weak amplitudes over Central Asia.

    Hong and Lu(2016)found that the JMD is significantly related to the Silk Road Pattern(SRP),which is a well-known teleconnection pattern in summer along the Asian jet in the form of alternate southerly and northerly anomalies,and behaves as the leading mode of upper-tropospheric meridional wind anomalies(e.g.,Lu et al.,2002;Sato and Takahashi,2006;Kosaka et al.,2009;Yasui and Watanabe,2010;Chen and Huang,2012).Corresponding to a northward(southward)JMD,the SRP tends to appear as a Specific phase,that is,there are southerly(northerly)anomalies over the Caspian Sea.Hong and Lu(2016)suggested that the JMD may trigger the meridional wind anomalies at the entrance of the Asian jet,and further induce the downstream SRP.This JMD–SRP connection therefore results in a combination of the JMD-and SRP-related circulation anomalies manifested by two significant anomalous anticyclones(cyclones)in the midlatitudes,with one located over West Asia and the other over East Asia,but a weak cyclone(anticyclone)over Central Asia(Hong and Lu,2016).These circulation anomalies—particularly the two anticyclonic(or cyclonic)anomalies over West Asia and East Asia—have been identified as the dominant pattern in the leading mode of upper-tropospheric hori-zontal winds(Hong and Lu,2016).

    The SRP propagates along the Asian jet,and its wave source has been demonstrated in various previous studies to be around the entrance of the Asian jet(e.g.,Enomoto et al.,2003;Sato and Takahashi,2006;Kosaka et al.,2009;Yasui and Watanabe,2010;Chen and Huang,2012;Lin and Lu,2016).For instance,through model experiments,Yasui and Watanabe(2010)indicated that the diabatic forcing over the Black Sea and the Caspian Sea contributes to more than 60%of the SRP response.Chen and Huang(2012)suggested that the SRP is triggered through the advection of vorticity by the divergent flow over this domain.Therefore,changes in the wave source here might be crucial for the variability of the SRP and the JMD–SRP relationship.

    The present study is an extension of Hong and Lu(2016).We show that a positive JMD–SRP relationship exists during northward JMD years but not during southward years(section 3).Furthermore,the possible mechanisms underlying this asymmetry in the JMD–SRP relationship are explored(section 4).Conclusions are given in section 5.

    2.Data and methods

    The main data source for this study is the monthly NCEP–NCAR reanalysis dataset(Kalnay et al.,1996),with a 2.5°×2.5°horizontal resolution.The analysis period for all data is 1958–2014,and the average calculations for June–July–August(JJA)are used to represent summer.

    We define separate indexes for the JMD and SRP,as follows.In line with Yasui and Watanabe(2010),we define the SRP index(SRPI)as the principal component(PC)of the leading mode for the 200-hPa meridional wind(V200)anomalies,which is obtained through empirical orthogonal function(EOF)analysis in the domain(20°–60°N,0°–150°E).A positive(negative)SRPI indicates an SRP phase with southerly(northerly)anomalies over the Caspian Sea.In addition,we define the JMD index (JMDI) as the difference in 200-hPa zonal wind(U200)anomalies between the domains(40°–55°N,40°–150°E)and(25°–40°N,40°–150°E),which are located to the north and south of the climatological jet axis,respectively,following Hong and Lu(2016).A positive(negative)JMDI indicates a northward(southward)JMD.To obtain robust results,we also use another index for the JMD,i.e.,the PC of the first EOF mode of U200 anomalies over the domain(20°–60°N,0°–150°E)(hereafter,Upc1).

    The main method used in the present study is composite analysis.Positive and negative cases are defined as the maximum 14 years and minimum 14 years,respectively,based on the JMDI and Upc1(Table 1).The same number for positive and negative cases is used deliberately to guarantee a reliable description of the asymmetry in the JMD–SRP relationship.This total of 28 years constitutes about half of the entire analysis period,and the residual 29 years are considered as“normal”cases.These normal cases are used as a reference,i.e.,the anomalies for positive/negative cases are the differences between these cases and normal cases.Besides using 14 selected years,we also checked other thresholds andobtained similar results(not shown).The Student’s t-test is used to test the statistical significance of the analyzed results.

    Table 1.Selected cases for the positive and negative JMDI and Upc1.Bold and italic font delineate the discrepancies for the positive and negative phases,respectively,when selecting the cases.

    3.Asymmetry in the relationship between the JMD and SRP

    Figure 1 shows the composite V200 and U200 anomalies for the positive and negative JMDI cases(hereafter referred to as“JMDI+”and “JMDI?”cases),respectively.The V200 anomalies in the JMDI+cases(Fig.1a)are characterized by a well-organized wave-like pattern along the Asian jet,from West Asia to East Asia.The amplitudes of these alternate anomalies are roughly the same for each cell,which are all above 3 m s?1.The wave-like pattern along the Asian jet agrees with the SRP-related anomalies, which are represented by the red contours,and coincide well with the anomalies related to the SRP in previous studies(e.g.,Lu et al.,2002;Yasui and Watanabe,2010).By contrast,there is no clear wavelike pattern along the Asian jet for the JMDI?cases(Fig.1b).The V200 anomalies tend to be scattered to the north and south of the Asian jet,but are weak along the jet.These anomalies differ remarkably from the SRP-related anomalies(also shown as red contours in Fig.1b),even showing opposite signs over the Caspian Sea and Central Asia.Therefore,the relationship between the JMD and SRP is close for positive JMDI years but vague for negative years.

    On the other hand,the composite U200 anomalies for the JMDI+and JMDI?cases also show evident distinctions.For the JMDI+cases(Fig.1c),the most prominent feature of the U200 anomalies is a north–south seesaw pattern,with westerly and easterly anomalies on the two sides of the climatological mean jet axis,respectively.Consistent with previous studies(Du et al.,2016;Hong and Lu,2016),this seesaw pattern appears much stronger over West Asia and East Asia,with the amplitudes of the centers around these regions being above 7 m s?1.Meanwhile,the anomalies over Central Asia are quite weak,as indicated by the strength being only about 2 m s?1.Hong and Lu(2016)suggested that the JMD can be further promoted by the SRP,manifested by intensi fi ed wind anomalies over West Asia and East Asia but weakened anomalies over Central Asia.Therefore,the inhomogeneous U200 anomalies in the zonal direction for the JMDI+cases(Fig.1c)can also be viewed as a manifestation of a close JMD–SRP relationship.The U200 anomalies in the JMDI?cases(Fig.1d)also exhibit a seesaw pattern,similar to that for the JMDI+cases.However,these anomalies are much weaker and tend to shift westward.

    Fig.1.Composite JJA mean 200-hPa anomalies of(a,b)meridional wind and(c,d)zonal wind(contours;units:m s?1)for(a,c)positive and(b,d)negative JMDI cases.Contour intervals are 1 m s?1and zero contours are omitted.Solid and dashed contours indicate positive and negative anomalies,respectively.Thick lines delineate the jet axes for the(a,c)positive and(b,d)negative JMDI cases,which are determined by the maximum zonal wind being greater than 25 m s?1.Red solid(dashed)contours denote regions where the composite V200 anomalies exceed 2 m s?1(fall below?2 m s?1)for the(a)maximum 14 years and(b)minimum 14 years,based on the SRPI.Shading indicates regions of significance at the 0.05 level,according to the Student’s t-test.

    To check that the above mentioned asymmetry in the JMD–SRP relationship is not merely the result of a particular index used,we adopt Upc1 as another index for the JMD.This index is similar to the JMDI,with the correlation coefficient between them being as high as 0.96,but the two indexes differ appreciably,as explained later.Figure 2 shows the composite V200 and U200 anomalies for the positive and negative Upc1 cases(hereafter referred to as“Upc1+”and“Upc1?”cases),respectively.It is clear that for the Upc1+cases(Figs.2a and c),there is a well-organized wave-like pattern in V200 anomalies along the Asian jet,which nearly overlaps with the anomalies related to the SRP(red contours in Fig.2a),and the north–south seesaw pattern of U200 anomalies is stronger over West Asia and East Asia.These anomalies show great resemblance to those for the JMDI+cases(Figs.1a and c),verifying the robustness of the JMD–SRP relationship in northward JMD years.By comparison,the anomalies for the Upc1?cases(Figs.2b and d)show much weaker intensity compared to those for the Upc1+cases.Hardly any significant V200 anomalies exist along the Asian jet over the western Eurasian continent,deviating greatly from those related to the SRP(red contours in Fig.2b),and the intensities of the U200 anomalies are roughly only around half those for the Upc1+cases.

    Figures 1 and 2 indicate that when different indexes are used for the JMD,the related anomalies are similar for the northward JMD but quite different for the southward JMD cases.The reason for this similarity and difference is that the cases selected by the two indexes are similar for the northward JMD cases but quite different for the southward JMD cases,as shown in Fig.3 and Table 1.Almost all cases(13 of 14)are the same for the northward JMD cases,but 10 cases in total are distinct for the southward JMD cases.specifically,1960,1968,1983,1988 and 2003 are included in the JMDI?cases but are absent in the Upc1?cases;while 1965,1974,1991,2002 and 2012 are selected for the Upc1?cases but are not chosen for the JMDI?cases.

    In addition, the anomalies tend to be more spatially coherent for the northward JMD cases than the southward cases.The pattern correlation coefficients in the U200 anomalies within the domain(20°–60°N,0°–150°E)between individual cases and the corresponding composites(e.g.,positive cases with the positive composite and negative cases with the negative composite)are generally greater for the positive cases than the negative cases(Table 2).

    Fig.2.As in Fig.1 but for composites based on the Upc1.

    Fig.3.Time series of standardized JMDI(bars)and Upc1(marked line).The thin horizontal lines indicate the thresholds selecting the JMDI+and JMDI?cases.

    Table 2.Pattern correlations coefficients(c.c.)in U200 anomalies within the domain(20°–60°N,0°–150°E)between individual selected cases(see Table 1)and the corresponding composites.

    Scatterplots of standardized JMDI and Upc1 with the SRPI,as shown in Fig.4,better illustrate the asymmetry in the JMD–SRP relationship.In general,there is a positive relationship between the JMD and SRP:the SRPI is positive(negative)in most positive(negative)JMDI/Upc1 years.This positive relationship is confirmed by the correlation coefficient between the JMDI(Upc1)and SRPI,which is 0.40(0.54)and significant at the 0.01 level.However,Fig.4 also indicates an interesting detail of this positive relationship:with increment of the JMDI,the SRPI increases obviously when the JMDI is positive but tends to decrease when the JMDI is negative(Fig.4a).This can be verified by the correlation coefficients between the two indexes,which drop dramatically from 0.60 for the positive JMDI years to?0.16 for the negative JMDI years.The correlation coefficient is negative,although weak, for the negative JMDI years. The scatterplot of Upc1 and SRPI(Fig.4b)shows quite similar features to Fig.4a,with contrasting correlation coefficients of 0.62 and?0.04 for positive and negative Upc1 years,respectively.

    All the above results on the asymmetry in the JMD–SRP relationship are obtained from the viewpoint of JMD indexes.If this relationship is viewed from the perspective of positive and negative SRPI years,the asymmetry can also be found.The correlation coefficient between JMDI(Upc1)and SRPI is 0.55(0.59)and significant at the 0.01 level in positive SRPI years,but is only 0.04/0.20 in negative SRPI years.For brevity,we only show the results based mainly on the JMD indexes in this paper.

    4.Role of the Rossby wave source in the asymmetry of the JMD–SRP relationship

    The Rossby wave source(RWS),which plays a primary role in inducing the SRP,has been demonstrated to be located around the entrance of the Asian westerly jet(Enomoto et al.,2003;Sato and Takahashi,2006;Kosaka et al.,2009;Yasui and Watanabe,2010;Chen and Huang,2012;Lin and Lu,2016).Consequently,changes in the jet stream—the JMD,for instance—may initiate corresponding changes of the RWS, and eventually influence the SRP and its connection with the JMD.Therefore,we next analyze the role played by the RWS in the asymmetry of the JMD–SRP relationship.

    Fig.4.Scatterplots of standardized JMDI and Upc1 with the SRPI.“cc”indicates the correlation coefficient between the JMDI/Upc1 and SRPI,and the subscripts“+”and“?”denote the positive JMD years and negative JMD years,respectively.

    Fig.5.As in Figs.1c–d and Figs.2c–d but for the JJA mean RWS anomalies at 200 hPa(contours;units:10?10s?2).Contour intervals are 0.3×10?10s?2,and zero contours are omitted.Thick lines delineate the jet axes for the(a,c)northward and(b,d)southward JMD cases,respectively,which are determined by the maximum zonal wind being greater than 25 m s?1.

    Figure 5 shows the composite 200-hPa RWS(Sardeshmukh and Hoskins,1988)anomalies,based on the JMD indexes.In the above expression,f and ζ are the planetary vorticity and relative vorticity,respectively;denotes the divergent component of horizontal winds.

    The most prominent feature of the RWS for the JMD+cases is the significant anomalies over Central Asia:negative anomalies to the east of the Caspian Sea and positive ones to the south of the Lake Balkhash.There are also positive anomalies over the Mediterranean Sea and eastern Europe.These RWS anomalies are very similar to the SRP-related anomalies(Fig.6),and thus it can be deduced that the SRP is effectively triggered in northward JMD years and results in the close JMD–SRP connection.It should be mentioned that the SRP-related RWS anomalies are basically opposite in sign between the positive and negative SRPI years,so we show in Fig.6 the differences between the positive and negative years.In comparison with the JMD+cases,the RWS anomalies for the JMD?cases(Fig.5b)are much weaker and tend to be of the same sign,which may be associated with the negative correlation coefficient between the JMDI and SRPI when the JMDI is negative.On the other hand,the weak RWS anomalies for these cases suggest much less efficiency in the triggering of the SRP and the resultant JMD–SRP relationship.Results for the Upc1+composite(Fig.5c)present a high degree of similarity with those for the JMDI+composite(Fig.5a),in terms of both distribution and intensity.However,the distribution for the Upc1?cases(Fig.5d)appears much more chaotic.The noisy RWS anomalies for the southward JMD cases(Figs.5b and d)are possibly related to the smaller spatial consistency between the southward JMD cases than the northward cases.

    The RWS can be rewritten as

    Here,we ignore the small terms of vorticity tendency,vertical advection and twisting.Dχis the divergence and β =df/dy.

    Fig.6.Composite differences in the JJA mean 200-hPa RWS anomalies(contours;units:10?10s?2)between the maximum 14 and minimum 14 SRPI years.Contour intervals are 0.6×10?10s?2,and zero contours are omitted.Shading indicates regions of statistical significance at the 0.05 level,according to the Student’s t-test.The three boxes denote domains where the SRP-related RWS anomalies are most significant.The thick line delineates the climatological jet axis,which is determined by the maximum zonal wind being greater than 25 m s?1.

    Although we analyzed the horizontal distribution of each RWS term,for brevity we only show the values averaged over three domains(boxes in Fig.6),i.e.,(35°–42.5°N,10°–35°E),(35°–45°N,55°–65°E)and(35°–45°N,70°–85°E),to estimate the relative contribution of the four RWS terms.These values are shown in Fig.7.It is evident that,for both the JMDI+and Upc1+cases,term1,i.e.,the planetary vortex stretching term ?fDχ,plays the primary role in determining the RWS,for all three domains(Figs.7a and c).Term2 and term4,i.e.,?ζDχand ?Vχ·?ζ,also contribute to the RWS,but they are much weaker than term1.On the other hand,the values for the southward JMD cases are much weaker(Figs.7b and d),particularly for the Upc1?cases.The RWS anomalies for the JMDI?cases tend to be in phase with those for the northward JMD cases(Fig.7b),and those for the Upc1? cases are much weaker—consistent with the meridional wind anomalies shown in Figs.1 and 2.

    Fig.7.Domain-averaged values of the four terms of the JJA mean RWS anomalies at 200 hPa for the(a)positive JMDI cases,(b)negative JMDI cases,(c)positive Upc1 cases,and(d)negative Upc1 cases.

    Therefore,we focus on the anomalies of the planetary vortex stretching term(?fDχ)and show their horizontal distribution in Fig.8.significant and positive anomalies appear to the south of the Lake Balkhash,the Mediterranean Sea and eastern Europe,and negative anomalies appear around the Caspian Sea,for the northward JMD cases(Figs.8a and c).These anomalies resemble the RWS anomalies well in both distribution and intensity,con fi rming?fDχis the term that contributes most to the RWS,as shown in Fig.7.The?fDχanomalies are much weaker for the southward JMD cases(Figs.8b and d).

    Fig.8.As in Fig.5 but for anomalies in term1 of the RWS.

    Fig.9.As in Fig.5 but for anomalies in the JJA mean 200-hPa divergence(contours;units:10?6s?1).

    The planetary vortex stretching term is determined by the divergence,although can be modified by the planetary vorticity.Therefore,200-hPa divergence anomalies play a crucial role in inducing the asymmetry in the JMD–SRP relationship.The upper-tropospheric divergence anomalies are generally associated with precipitation anomalies through ascent or de-scent.Therefore,we analyzed the precipitation anomalies for the northward JMD cases,and found that there are negative anomalies in Central Asia(not shown).The region for these negative precipitation anomalies includes both anomalous convergence and divergence,shown in Figs.9a and c.Thus,the relationship between divergence and precipitation anomalies may be complicated in Central Asia,where the amount of climatological rainfall is generally light.The reason for 200-hPa divergence anomalies requires further investigation.

    5.Conclusions and discussion

    In this study,we have carried out further investigation into the relationship between the JMD and SRP in summer,which had been earlier identified by Hong and Lu(2016).The present results indicate that the JMD–SRP relationship is asymmetric:the SRP becomes stronger with a more northward JMD,resulting in a significant JMD–SRP relationship.However,the relationship is weak in southward JMD years.This asymmetric relationship is confirmed by the correlation coefficients between the indexes of JMD and SRP in northward and southward JMD years,which are 0.60 and?0.16 with respect to JMDI–SRPI,and 0.62 and ?0.04 with respect to Upc1–SRPI(Upc1 being another indicator quantifying the JMD),during the analysis period of this study(1958–2014).Due to this asymmetry of the relationship,the uppertropospheric zonal wind anomalies are much stronger in West Asia and East Asia and weaker in Central Asia for northward JMD years,but tend to be not well-organized for southward JMD years.

    Further results suggest that northward JMD years correspond to much stronger RWS anomalies,which are primarily determined by the planetary vortex stretching term?fDχ,around the entrance of the Asian jet.specifically,there are significant positive anomalies to the south of Lake Balkhash and over the Mediterranean Sea,and negative anomalies to the east of the Caspian Sea—closely coherent with those associated with the SRP.However,the RWS anomalies are much weaker and not well-organized for southward JMD years. Therefore,we conclude that the RWS anomalies around the entrance of the Asian jet play a crucial role in inducing the asymmetry of the JMD–SRP relationship.

    However,the mechanisms responsible for the asymmetry of the RWS or vortex stretching term ?fDχanomalies between northward and southward JMD years,remain unknown.This may be related to the internal processes between the JMD and RWS,but possibly also to the greater spatial coherence between northward JMDs than the southward ones.In addition,the asymmetry in the relationship between the JMD and SRP implies that their impacts on climate anomalies may also be distinct between positive and negative JMDI/SRPI years.These issues require further investigation in future studies.

    Acknowledgements.We thank the two anonymous reviewers for their comments,which triggered the results shown in section 4 and helped considerably in improving the expression of our findings.This work was supported by the National Natural Science Foundation of China(Grant Nos.41320104007,41421004,and 41731177).

    Chen,G.S.,and R.H.Huang,2012:Excitation mechanisms of the teleconnection patterns affecting the July precipitation in Northwest China.J.Climate,25,7834–7851,https://doi.org/10.1175/JCLI-D-11-00684.1.

    Du,Y.,T.Li,Z.Q.Xie,and Z.W.Zhu,2016:Interannual variability of the Asian subtropical westerly jet in boreal summer and associated with circulation and SST anomalies.Climate Dyn.,46,2673–2688,https://doi.org/10.1007/s00382-015-2723-x.

    Enomoto,T.,B.J.Hoskins,and Y.Matsuda,2003:The formation mechanism of the Bonin high in August.Quart.J.Roy.Meteor.Soc.,129,157–178,https://doi.org/10.1256/qj.01.211.

    Hong,X.W.,and R.Y.Lu,2016: The meridional displacement of the summer Asian jet,Silk Road Pattern,and tropical SST anomalies.J.Climate,29,3753–3766,https://doi.org/10.1175/JCLI-D-15-0541.1.

    Kalnay,E.,and Coauthors,1996:The NCEP/NCAR 40-year reanalysis project.Bull.Amer.Meteor.Soc.,77,437–471,https://doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

    Kosaka,Y.,H.Nakamura,M.Watanabe,and M.Kimoto,2009:Analysis on the dynamics of a wave-like teleconnection pattern along the summertime Asian jet based on a reanalysis dataset and climate model simulations.J.Meteor.Soc.Japan,87,561–580,https://doi.org/10.2151/jmsj.87.561.

    Lin,Z.D.,and R.Y.Lu,2016:Impact of summer rainfall over southern-central Europe on the circumglobal teleconnection.Atmospheric Science Letters,17,258–262,https://doi.org/10.1002/asl.652.

    Lu,R.-Y.,J.-H.Oh,and B.-J.Kim,2002:A teleconnection pattern in upper-level meridional wind over the North African and Eurasian continent in summer.Tellus A,54,44–55,https://doi.org/10.1034/j.1600-0870.2002.00248.x10.3402/tellusa.v54i1.12122.

    Sardeshmukh,P.D.,and B.J.Hoskins,1988:The generation of global rotational fl ow by steady idealized tropical divergence.J.Atmos.Sci.,45,1228–1251,https://doi.org/10.1175/1520-0469(1988)045<1228:TGOGRF>2.0.CO;2.

    Sato,N.,and M.Takahashi,2006:Dynamical processes related to the appearance of quasi-stationary waves on the subtropical jet in the midsummer northern hemisphere.J.Climate,19,1531–1544,https://doi.org/10.1175/JCLI3697.1.

    Wei,W.,R.H.Zhang,M.Wen,and S.Yang,2017:Relationship between the Asian westerly jet stream and summer rainfall over Central Asia and North China:Roles of the Indian monsoon and the South Asian High.J.Climate,30,537–552,https://doi.org/10.1175/JCLI-D-15-0814.1.

    Yasui,S.,and M.Watanabe,2010:Forcing processes of the summertime circumglobal teleconnection pattern in a dry AGCM.J.Climate,23,2093–2114,https://doi.org/10.1175/2009JCLI3323.1.

    国产精品福利在线免费观看| 久久综合国产亚洲精品| 色94色欧美一区二区| 男人狂女人下面高潮的视频| 久久精品久久精品一区二区三区| 在线观看国产h片| 亚洲av成人精品一二三区| 免费高清在线观看视频在线观看| 亚洲成人av在线免费| 99re6热这里在线精品视频| 另类精品久久| 色婷婷av一区二区三区视频| 亚洲一区二区三区欧美精品| 国产永久视频网站| 亚洲第一区二区三区不卡| 日本黄色日本黄色录像| 国产日韩欧美亚洲二区| 一本大道久久a久久精品| 亚洲av电影在线观看一区二区三区| 美女主播在线视频| 欧美日韩综合久久久久久| 国产欧美日韩精品一区二区| 欧美 日韩 精品 国产| videossex国产| 久久久久久伊人网av| 如日韩欧美国产精品一区二区三区 | 97在线人人人人妻| 少妇人妻精品综合一区二区| 中文字幕av电影在线播放| 五月玫瑰六月丁香| 欧美日韩在线观看h| av.在线天堂| 蜜桃在线观看..| 亚洲精品一区蜜桃| 少妇人妻精品综合一区二区| 欧美日韩国产mv在线观看视频| a级毛片免费高清观看在线播放| 看非洲黑人一级黄片| 纵有疾风起免费观看全集完整版| 99re6热这里在线精品视频| 欧美一级a爱片免费观看看| 婷婷色av中文字幕| 色婷婷久久久亚洲欧美| 精品少妇久久久久久888优播| 又粗又硬又长又爽又黄的视频| 国产真实伦视频高清在线观看| 亚洲第一区二区三区不卡| 欧美日本中文国产一区发布| 国产色婷婷99| 少妇被粗大的猛进出69影院 | 新久久久久国产一级毛片| 久久人人爽人人片av| 插阴视频在线观看视频| 我要看日韩黄色一级片| 一级毛片 在线播放| 汤姆久久久久久久影院中文字幕| 欧美最新免费一区二区三区| 蜜桃在线观看..| 91成人精品电影| 久热久热在线精品观看| 中文乱码字字幕精品一区二区三区| 国产综合精华液| 国产男人的电影天堂91| 又大又黄又爽视频免费| 亚洲国产精品一区二区三区在线| 欧美日本中文国产一区发布| 久久婷婷青草| 国产成人精品一,二区| 欧美成人精品欧美一级黄| 99热6这里只有精品| 婷婷色综合www| 91精品国产国语对白视频| 99九九线精品视频在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 免费观看的影片在线观看| 日本av手机在线免费观看| 精品一区在线观看国产| 最近最新中文字幕免费大全7| 亚洲欧洲国产日韩| 自拍偷自拍亚洲精品老妇| 精品人妻偷拍中文字幕| 人体艺术视频欧美日本| 国产黄片美女视频| 人妻制服诱惑在线中文字幕| 亚洲精华国产精华液的使用体验| 少妇精品久久久久久久| 纯流量卡能插随身wifi吗| 丰满饥渴人妻一区二区三| 日韩亚洲欧美综合| 好男人视频免费观看在线| 观看美女的网站| 精品一区在线观看国产| 在线观看三级黄色| 久久久久久久久大av| 看非洲黑人一级黄片| 日本与韩国留学比较| 午夜日本视频在线| 日韩欧美 国产精品| 美女脱内裤让男人舔精品视频| 日韩一本色道免费dvd| 精品久久久久久久久亚洲| 少妇人妻 视频| 亚洲美女黄色视频免费看| 男的添女的下面高潮视频| 精品熟女少妇av免费看| 成人亚洲精品一区在线观看| 免费大片黄手机在线观看| 色5月婷婷丁香| 人妻 亚洲 视频| 免费少妇av软件| 九九爱精品视频在线观看| 国产无遮挡羞羞视频在线观看| 久久久久国产精品人妻一区二区| 亚洲欧美清纯卡通| av不卡在线播放| 91久久精品国产一区二区成人| 亚洲av电影在线观看一区二区三区| a级毛片免费高清观看在线播放| 五月玫瑰六月丁香| 亚洲国产精品专区欧美| 日韩精品有码人妻一区| 亚洲天堂av无毛| 久久国产乱子免费精品| 欧美亚洲 丝袜 人妻 在线| 婷婷色麻豆天堂久久| 亚洲精品aⅴ在线观看| 在线看a的网站| 男女免费视频国产| 国产成人精品无人区| 三级国产精品欧美在线观看| 肉色欧美久久久久久久蜜桃| 免费大片黄手机在线观看| 国产毛片在线视频| a级片在线免费高清观看视频| 99国产精品免费福利视频| 欧美+日韩+精品| 久久av网站| 欧美+日韩+精品| 久久精品国产亚洲av涩爱| 亚洲国产精品999| 国产又色又爽无遮挡免| a级一级毛片免费在线观看| 国产精品一二三区在线看| 久久久久国产精品人妻一区二区| 在线免费观看不下载黄p国产| 最近的中文字幕免费完整| 人妻系列 视频| 一级黄片播放器| 久久精品国产亚洲网站| 在线观看www视频免费| 亚洲成人一二三区av| 91精品国产九色| 亚洲国产欧美日韩在线播放 | 一级毛片黄色毛片免费观看视频| 一级毛片黄色毛片免费观看视频| 亚洲第一区二区三区不卡| 22中文网久久字幕| 男人舔奶头视频| 久久99精品国语久久久| 中文欧美无线码| 男女啪啪激烈高潮av片| 久久热精品热| 亚洲精品视频女| 国产极品天堂在线| 在线观看免费日韩欧美大片 | 天天躁夜夜躁狠狠久久av| 国产成人精品久久久久久| 国产免费福利视频在线观看| 国产女主播在线喷水免费视频网站| 三级国产精品欧美在线观看| 在线观看av片永久免费下载| 日韩人妻高清精品专区| 午夜日本视频在线| 日本av免费视频播放| 看非洲黑人一级黄片| 亚洲怡红院男人天堂| 黑人高潮一二区| 大片电影免费在线观看免费| 丰满乱子伦码专区| 欧美性感艳星| 亚洲,一卡二卡三卡| 亚洲色图综合在线观看| 中文字幕av电影在线播放| 欧美日韩国产mv在线观看视频| 99久国产av精品国产电影| 国产成人aa在线观看| 国产精品一二三区在线看| 嘟嘟电影网在线观看| 久久久久久人妻| 少妇精品久久久久久久| 日韩一区二区视频免费看| 在线播放无遮挡| 亚洲国产日韩一区二区| 欧美日韩综合久久久久久| 永久网站在线| 男女啪啪激烈高潮av片| 如日韩欧美国产精品一区二区三区 | 日日摸夜夜添夜夜爱| 99热全是精品| 免费在线观看成人毛片| 嘟嘟电影网在线观看| 自线自在国产av| 国产乱人偷精品视频| 欧美精品高潮呻吟av久久| 日韩电影二区| 欧美老熟妇乱子伦牲交| 人妻制服诱惑在线中文字幕| 国产av精品麻豆| 午夜福利网站1000一区二区三区| 亚洲精品乱码久久久久久按摩| 国产一区二区三区av在线| 亚洲欧洲国产日韩| 亚洲一级一片aⅴ在线观看| 亚洲av福利一区| 2018国产大陆天天弄谢| 国国产精品蜜臀av免费| 久久ye,这里只有精品| 日韩成人av中文字幕在线观看| 中文字幕av电影在线播放| 婷婷色综合www| 色婷婷久久久亚洲欧美| 亚洲精品第二区| 嘟嘟电影网在线观看| 久久精品国产a三级三级三级| videos熟女内射| 国产一区有黄有色的免费视频| 男人添女人高潮全过程视频| 亚洲欧美成人精品一区二区| 成年人免费黄色播放视频 | 国产男女超爽视频在线观看| 天堂中文最新版在线下载| av网站免费在线观看视频| 亚洲真实伦在线观看| 亚洲伊人久久精品综合| 一级a做视频免费观看| 观看美女的网站| 精品国产国语对白av| 国产精品熟女久久久久浪| 99热全是精品| 纯流量卡能插随身wifi吗| 国产日韩欧美在线精品| 赤兔流量卡办理| 一区二区三区乱码不卡18| 人妻少妇偷人精品九色| 免费看日本二区| 涩涩av久久男人的天堂| 久热这里只有精品99| 国产一区二区三区综合在线观看 | 精品国产乱码久久久久久小说| 日韩制服骚丝袜av| 黑人猛操日本美女一级片| av国产久精品久网站免费入址| 最近的中文字幕免费完整| 亚洲人成网站在线观看播放| 欧美xxⅹ黑人| 亚洲欧洲国产日韩| 久久人人爽av亚洲精品天堂| 一区二区三区乱码不卡18| 亚洲国产精品国产精品| 亚洲久久久国产精品| 日韩人妻高清精品专区| 亚洲综合色惰| 午夜免费鲁丝| 欧美xxxx性猛交bbbb| 久久久久久久亚洲中文字幕| 中文字幕人妻丝袜制服| 美女主播在线视频| 亚洲国产精品一区三区| 久久精品夜色国产| h视频一区二区三区| 成人毛片a级毛片在线播放| 色婷婷av一区二区三区视频| 色视频在线一区二区三区| 简卡轻食公司| 熟妇人妻不卡中文字幕| 国产毛片在线视频| 男女啪啪激烈高潮av片| 国产精品一区二区在线不卡| 日韩亚洲欧美综合| 欧美成人午夜免费资源| 国产一区二区三区综合在线观看 | 日本欧美国产在线视频| 自线自在国产av| 十八禁网站网址无遮挡 | 2021少妇久久久久久久久久久| 欧美成人午夜免费资源| 伊人久久精品亚洲午夜| 日产精品乱码卡一卡2卡三| 噜噜噜噜噜久久久久久91| 99热6这里只有精品| 日日摸夜夜添夜夜添av毛片| 日本av免费视频播放| 青青草视频在线视频观看| 丰满饥渴人妻一区二区三| 国产又色又爽无遮挡免| 男人舔奶头视频| 九色成人免费人妻av| 午夜福利影视在线免费观看| 亚洲美女搞黄在线观看| 欧美日本中文国产一区发布| 欧美日韩在线观看h| 少妇裸体淫交视频免费看高清| 黄色视频在线播放观看不卡| 一级毛片aaaaaa免费看小| 国产免费福利视频在线观看| 国产亚洲91精品色在线| 精品国产国语对白av| 少妇精品久久久久久久| av天堂中文字幕网| xxx大片免费视频| 日韩免费高清中文字幕av| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产精品一区三区| 欧美xxxx性猛交bbbb| 日日啪夜夜爽| 免费av不卡在线播放| 26uuu在线亚洲综合色| 美女xxoo啪啪120秒动态图| 特大巨黑吊av在线直播| 99国产精品免费福利视频| 国产黄色免费在线视频| 国产亚洲最大av| 插阴视频在线观看视频| 精品99又大又爽又粗少妇毛片| 日本黄色日本黄色录像| 只有这里有精品99| av天堂久久9| 最近的中文字幕免费完整| 青春草视频在线免费观看| 九九在线视频观看精品| 大码成人一级视频| 亚洲av免费高清在线观看| 亚洲av在线观看美女高潮| 五月开心婷婷网| 观看av在线不卡| 久久久久久久精品精品| 亚洲国产成人一精品久久久| 国产成人freesex在线| 国产精品偷伦视频观看了| 久久午夜综合久久蜜桃| 国产美女午夜福利| 一级毛片aaaaaa免费看小| 久久久午夜欧美精品| 在线观看国产h片| 少妇熟女欧美另类| 另类亚洲欧美激情| 久久热精品热| 中国国产av一级| 日韩一区二区三区影片| 欧美日本中文国产一区发布| 久久精品久久久久久噜噜老黄| 99九九在线精品视频 | 美女国产视频在线观看| 高清不卡的av网站| 午夜精品国产一区二区电影| 国产成人精品福利久久| 亚洲欧美成人综合另类久久久| 内地一区二区视频在线| 看十八女毛片水多多多| 日日撸夜夜添| 亚洲精品日本国产第一区| 青春草视频在线免费观看| 日韩av在线免费看完整版不卡| 国产高清国产精品国产三级| 亚洲国产成人一精品久久久| 久久鲁丝午夜福利片| 一区二区av电影网| 亚洲va在线va天堂va国产| 国产高清三级在线| 在线观看免费日韩欧美大片 | 久久久久人妻精品一区果冻| 激情五月婷婷亚洲| 大香蕉97超碰在线| 丝袜在线中文字幕| 久久av网站| 亚洲av国产av综合av卡| 亚洲欧洲精品一区二区精品久久久 | 亚洲高清免费不卡视频| 乱人伦中国视频| 插逼视频在线观看| 国产精品三级大全| 国产日韩欧美视频二区| 看十八女毛片水多多多| 菩萨蛮人人尽说江南好唐韦庄| 夜夜骑夜夜射夜夜干| 久久久精品免费免费高清| 最新的欧美精品一区二区| 亚洲精品久久午夜乱码| 久久久久久人妻| 大陆偷拍与自拍| 亚洲国产欧美在线一区| 国产乱来视频区| www.av在线官网国产| 最近最新中文字幕免费大全7| 曰老女人黄片| 国产成人一区二区在线| 国产极品粉嫩免费观看在线 | 中文欧美无线码| 91精品伊人久久大香线蕉| 简卡轻食公司| 国产91av在线免费观看| av.在线天堂| 9色porny在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品乱码久久久久久按摩| 一区二区三区乱码不卡18| 91精品国产九色| 亚洲欧美精品自产自拍| av一本久久久久| 在线观看人妻少妇| 男人舔奶头视频| 日日爽夜夜爽网站| 久久久国产欧美日韩av| 日韩一本色道免费dvd| 成人18禁高潮啪啪吃奶动态图 | 精品久久国产蜜桃| 亚洲国产欧美日韩在线播放 | 亚洲欧美清纯卡通| 91精品一卡2卡3卡4卡| 精品一区二区免费观看| 亚洲国产精品成人久久小说| 伦精品一区二区三区| 国产精品国产三级国产av玫瑰| 狂野欧美白嫩少妇大欣赏| 人妻制服诱惑在线中文字幕| 亚洲自偷自拍三级| 嘟嘟电影网在线观看| 中文字幕人妻丝袜制服| 亚洲欧洲精品一区二区精品久久久 | 又爽又黄a免费视频| 欧美亚洲 丝袜 人妻 在线| 亚洲成人手机| 国产日韩欧美在线精品| 只有这里有精品99| 精品卡一卡二卡四卡免费| 国产伦理片在线播放av一区| 精华霜和精华液先用哪个| 国语对白做爰xxxⅹ性视频网站| h视频一区二区三区| 日本色播在线视频| 久久久久网色| 五月开心婷婷网| 大片电影免费在线观看免费| 国产日韩欧美视频二区| 老司机影院毛片| 亚洲精品国产成人久久av| 极品教师在线视频| 久久久久久久久久久久大奶| 91精品国产国语对白视频| 大又大粗又爽又黄少妇毛片口| 久久ye,这里只有精品| 日本爱情动作片www.在线观看| 极品少妇高潮喷水抽搐| 亚洲av国产av综合av卡| 91精品伊人久久大香线蕉| 男女无遮挡免费网站观看| 亚洲美女黄色视频免费看| 欧美 亚洲 国产 日韩一| 中文资源天堂在线| 美女国产视频在线观看| 国产精品久久久久久精品电影小说| 中文资源天堂在线| 中国三级夫妇交换| 成人二区视频| 欧美xxⅹ黑人| 色视频在线一区二区三区| 一级爰片在线观看| 日韩中字成人| 水蜜桃什么品种好| 久久久久久久久久久免费av| 黑人猛操日本美女一级片| 国产日韩欧美亚洲二区| freevideosex欧美| 性色avwww在线观看| 欧美高清成人免费视频www| 日本vs欧美在线观看视频 | 99精国产麻豆久久婷婷| 精品一区在线观看国产| 99久久人妻综合| av播播在线观看一区| 性色avwww在线观看| 在线播放无遮挡| 国产亚洲91精品色在线| 亚洲av不卡在线观看| 亚洲精品aⅴ在线观看| 18禁在线无遮挡免费观看视频| 国产在线男女| 99久久中文字幕三级久久日本| 欧美日韩亚洲高清精品| 91午夜精品亚洲一区二区三区| 一区二区三区免费毛片| 欧美人与善性xxx| 少妇的逼好多水| 久久久国产欧美日韩av| 又大又黄又爽视频免费| 久久精品国产亚洲网站| 永久网站在线| 亚洲精品国产av成人精品| 精品人妻熟女毛片av久久网站| 久久久久视频综合| 欧美 日韩 精品 国产| 91aial.com中文字幕在线观看| 男人爽女人下面视频在线观看| 国产成人精品福利久久| 中国国产av一级| 日韩视频在线欧美| 99久久精品热视频| 国产亚洲一区二区精品| 天天躁夜夜躁狠狠久久av| 欧美日韩国产mv在线观看视频| 五月开心婷婷网| 国产高清有码在线观看视频| 91精品一卡2卡3卡4卡| 亚洲av福利一区| 国产色爽女视频免费观看| 免费播放大片免费观看视频在线观看| 免费看av在线观看网站| av网站免费在线观看视频| 插阴视频在线观看视频| 国产成人精品福利久久| 国产高清国产精品国产三级| 午夜免费观看性视频| 久久久久久久久大av| 欧美 日韩 精品 国产| 亚洲精品乱码久久久久久按摩| 高清欧美精品videossex| 大陆偷拍与自拍| 国精品久久久久久国模美| 国产午夜精品一二区理论片| 亚洲精品乱久久久久久| 九九爱精品视频在线观看| 深夜a级毛片| 欧美丝袜亚洲另类| 免费高清在线观看视频在线观看| 乱系列少妇在线播放| 午夜福利在线观看免费完整高清在| 黄色毛片三级朝国网站 | av免费在线看不卡| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 日韩 亚洲 欧美在线| 中文天堂在线官网| 日韩 亚洲 欧美在线| 久久免费观看电影| 老司机影院成人| 久久人妻熟女aⅴ| 插阴视频在线观看视频| 亚洲成色77777| 高清不卡的av网站| 久久久久久久久久成人| 狂野欧美激情性xxxx在线观看| 免费少妇av软件| 成人国产av品久久久| 熟女人妻精品中文字幕| 国产精品久久久久久久电影| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| 少妇裸体淫交视频免费看高清| 欧美老熟妇乱子伦牲交| 一级a做视频免费观看| 两个人免费观看高清视频 | 乱码一卡2卡4卡精品| 草草在线视频免费看| 18禁动态无遮挡网站| 美女主播在线视频| 国产精品蜜桃在线观看| 3wmmmm亚洲av在线观看| 亚洲美女搞黄在线观看| 99热国产这里只有精品6| 在线亚洲精品国产二区图片欧美 | 如日韩欧美国产精品一区二区三区 | 成人亚洲精品一区在线观看| 国产成人免费无遮挡视频| 国产精品一区二区三区四区免费观看| 日韩成人av中文字幕在线观看| 久久亚洲国产成人精品v| 在线观看国产h片| 国产成人午夜福利电影在线观看| 欧美少妇被猛烈插入视频| 国产午夜精品一二区理论片| 一边亲一边摸免费视频| 久久婷婷青草| a级毛片免费高清观看在线播放| 赤兔流量卡办理| 亚洲欧洲精品一区二区精品久久久 | www.av在线官网国产| 啦啦啦视频在线资源免费观看| 黄色配什么色好看| 一二三四中文在线观看免费高清| 黄色一级大片看看| 免费观看无遮挡的男女| 久久久久久久久久久免费av| 色婷婷久久久亚洲欧美| 美女大奶头黄色视频| 国产无遮挡羞羞视频在线观看| 看非洲黑人一级黄片| 又爽又黄a免费视频| 中文精品一卡2卡3卡4更新| 80岁老熟妇乱子伦牲交| 18禁在线播放成人免费| 国产一区二区三区综合在线观看 | 久久久久久久久久久久大奶| 少妇 在线观看| 综合色丁香网| 国产精品不卡视频一区二区| 精品少妇久久久久久888优播| 欧美日本中文国产一区发布| 少妇猛男粗大的猛烈进出视频| 亚洲成色77777| 亚洲综合精品二区| 日韩欧美 国产精品| 人妻系列 视频| 欧美日韩在线观看h| 高清毛片免费看| 黑人高潮一二区|