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

    Impact of Planetary Wave Reflection on Tropospheric Blocking over the Urals-Siberia Region in January 2008

    2016-11-25 03:10:19DebashisNATHandWenCHEN
    Advances in Atmospheric Sciences 2016年3期

    Debashis NATH and Wen CHEN

    Center for Monsoon System Research,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100190

    Impact of Planetary Wave Reflection on Tropospheric Blocking over the Urals-Siberia Region in January 2008

    Debashis NATH and Wen CHEN?

    Center for Monsoon System Research,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100190

    Planetary wave reflection from the stratosphere played a significant role in changing the tropospheric circulation pattern over Eurasia in mid-January 2008.We studied the 2008 event and compared with composite analysis(winters of 2002/2003, 2004/2005,2006/2007,2007/2008,2010/2011 and 2011/2012),when the downward coupling was stronger,by employing time-lagged singular value decomposition analysis on the geopotential height field.In the Northern Hemisphere,the geopotential fields were decomposed into zonal mean and wave components to compare the relative covariance patterns.It was found that the wavenumber 1(WN1)component was dominant compared with the wavenumber 2(WN2)component and zonal mean process.For the WN1 field,the covariance was much higher(lower)for the negative(positive)lag,with a prominent peak around+15 days when the leading stratosphere coupled strongly with the troposphere.It contributed to the downward coupling due to reflection,when the stratosphere exhibited a partially reflective background state.We also analyzed the evolution of the WN1 anomaly and heat flux anomaly,both in the troposphere and stratosphere,during January–March 2008.The amplitude of the tropospheric WN1 pattern reached a maximum and was consistent with a downward wave coupling event influenced by the stratospheric WN1 anomaly at 10 hPa.This was consistent with the reflection of the WN1 component over Eurasia,which triggered an anomalous blocking high in the Urals–Siberia region.We further clarified the impact of reflection on the tropospheric WN1 field and hence the tropospheric circulation pattern by changing the propagation direction during and after the event.

    planetary wave reflection,singular value decomposition(SVD),Eliassen–Palm(EP)flux,blocking,circulation

    1.Introduction

    In recent decades,several studies have related tropospheric variability with the downward propagation of stratospheric anomalies through planetary wave reflection(Kodera et al.,2008).Due to the gradual increase in atmospheric pressure downward,the reflected component gets attenuated faster and its impact is considered to be minimal on the tropospheric regime.The theory of planetary wave reflection on tropospheric fields was initially proposed by Hines (1974)and Geller and Alpert(1980).Later,several authors discussed the coupling processes and downward propagation in the light of Northern Annular modes(Baldwin and Dunkerton,1999,2001),wave–mean flow interaction(Baldwin and Dunkerton,2001;Christiansen,2001)and anomalous propagation of the mean zonal wind field(Kodera et al., 1990;Kuroda and Kodera,1999;Christiansen,2000)down to the troposphere.Perlwitz and Graf(2001)and Perlwitz and Harnik(2003,2004)statistically described the verticalcoupling of the wavenumber 1(WN1)and wavenumber 2 (WN2)components downward with the tropospheric height fields.In winter the stratosphere is either reflective or nonreflective based on the strength of the polar vortex(Perlwitz and Harnik,2004).The WN1 and zonal mean component play a major role when the background is reflective.Furthermore,employing time-lagged correlation analysis,Perlwitz and Harnik(2004)explained the gradual tilt in phase of the WN1 regression pattern(leading mode)vertically upward.With negative time lag the tilt is westward,but with positive lag it is eastward,and this feature is consistent with the downward propagation of planetary waves due to reflection.

    Shaw et al.(2010)described the characteristics of downward wave coupling between the stratosphere and troposphere,using the 40-year ECMWF(European Centre for Medium-Range Weather Forecasts)reanalysis dataset.Employing cross-spectral correlation analysis and wave geometry diagnostics,they found that the downward WN1 coupling occurs both in the presence of a vertical reflecting surface in the mid-to-upper stratosphere and a high-latitude meridional waveguide in the lower stratosphere.They also discussed theimportanceof theseasonal cycleof thewavegeometry for the proper representation of downward wave coupling between the stratosphere and troposphere,both in the Northern and Southern Hemisphere.In a separate study,Shaw and Perlwitz(2013)statistically investigated the life cycle of Northern Hemisphere wave coupling events and found that it occurs over a period of 28 days.Furthermore,they showed that during the downward coupling process,there is a transition in stratospheric WN1 heat flux,from positive to negative,and the WN1 phase tilts from westward to eastward. Shaw et al.(2014)established a new dynamical metric of troposphere–stratosphere coupling,based on extreme stratospheric planetary-scale wave heat flux events.

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

    Coughlin and Tung(2005)demonstrated the possibility of wave reflection in the context of major sudden stratospheric warming(SSW)events and discussed its impact on the tropospheric weather regime(Nath et al.,2013).They illustrated the changes in the tropospheric WN1 field in response to the reflected component from the stratosphere. Separately,Kodera et al.(2008)related the occurrence of an extreme cold event in March 2007 over the northeast coast of the North American continent,to the upward and reflected component of planetary waves over Eurasia and the North American sector,respectively.During an SSW event in 1984–85,Kodera and Chiba(1995)investigated the changes in circulation pattern due to downward and equatorwardpropagationofmidlatitudeplanetarywavestothetroposphere.Geopotential anomalies that propagate downward to the troposphere have a significant impact on the tropospheric weather regime,particularly in non-reflective years(Perlwitz and Harnik,2004).Perlwitz and Harnik(2004)categorized the reflective and non-reflective basic states for planetary WN1 reflection based on the zonal-mean zonal wind difference between 2 and 10 hPa,averaged over 58?–74?N and over time.The reflective basic state corresponds to a negative index with the polar night jet peaking in the mid-stratosphere; whereas,for the non-reflective state,the zonal wind increases with increasing height.The planetary waves propagate upwardalongthestratosphericwesterlyjet,weakeningthepolar night jet in the upper stratosphere.This inhibits further propagation of the planetary waves high up in the stratosphere and it reflects back to the troposphere.During reflective winters the stratospheric signals are weak and get attenuated above the tropopause(Perlwitz and Harnik,2004).Moreover,in strong polar vortex winters,the WN1 reflection pattern is more prominent(Perlwitz and Graf,2001);whereas,in weak vortex years,stratosphere–troposphere coupling is relatively strong(Baldwin and Dunkerton,1999,2001).

    Here,we considered a specific case in the pre-warming phase of a major SSW event in January 2008.As reported previously(Hui,2009;Zhouetal.,2009;Nathetal.,2014),in January and early February 2008,parts of Eurasia and China experienced extreme cold events,snowfall and freezing rain, particularly in the southern part of China.These phenomena caused excessive damage,disruption and major infrastructure loss,resulting in broken power transmission lines and chaotic traffic conditions(Zhou et al.,2009).China experienced substantial economic losses of 53.8 billion RMB due to freezing rain alone.Between 10 January and 2 February 2008 there were four episodes of severe and persistent snow over the Yangtze River basin,South China,and Southwest China.The 2008 event was the coldest event since at least 1979,bringing about 107 casualties,according to the Ministry of Civil Affairs.

    Zhou et al.(2009)indicated the key factors as the occurrence of a persistent blocking high over Siberia,as well as strong and persistent southwesterly flow,which triggered moisture advection from the Bay of Bengal to southerncentral China,and the formation of a deep inversion layer in the lower troposphere.Furthermore,Hui(2009)attributed these adverse meteorological conditions with abnormal circulation anomalies at high latitudes.In a separate study, Nath et al.(2014)demonstrated that when the stratospheric basic state is partially reflective,a wave packet emanating from Baffin Island/the coast of Labrador propagates eastward,equatorward and reflects back over central Eurasia and parts of China,which in turn triggers the advection of cold wind from the northern part of the boreal forest region and Siberia to the subtropics.The extraordinary persistence of this particular cold event has been linked with anomalous blocking high over the Urals–Siberia region.

    Despite many previous studies having investigated the key factors(like tropospheric blocking)that triggered the extremecoldeventinJanuary2008,noneexploredthecausative mechanism underpinning the occurrence of the anomalous blocking high over the Urals–Siberia region.In the present study,we performed lagged correlation analysis to understand the respective contribution of the wave(WN1 and WN2)and zonal mean flow in conjunction with stratosphere–troposphere coupling processes.By analyzing the squared covariance between the stratospheric and tropospheric height field,we identified the key dates in January 2008 when the impact of the stratospheric basic state on the troposphere was at a maximum(or vice versa).We also compared the singular value decomposition(SVD)pattern of the 2008 event with the composite mean pattern,which included several winters for which downward wave coupling has been reported. Based on Shaw and Perlwitz(2013),Kodera et al.(2013), and Dunn-Sigouin and Shaw(2015),we chose the winters of 2003,2005,2007,2008,2011 and 2012 for the composite analysis.We also estimated the changes in the tropospheric WN1 field and hence the tropospheric circulation,during and after the event.Furthermore,we clarified the role of planetary wave reflection on the tropospheric circulation pattern and formation of strong Urals–Siberia blocks.

    2.Data and methodology

    2.1.Data

    Daily mean ECMWF Interim Reanalysis(ERA-Interim) data(Uppala et al.,2008;Dee et al.,2011)were used for potential vorticity,geopotential height,zonal wind,meridional wind,and temperature.The individual parameters werearchived from December to April from 2002/2003 to 2011/ 2012.The ERA-Interim data are available at 37 pressure levels from 1000 hPa to 1 hPa,with a horizontal resolution of 1.5?×1.5?.

    2.2.Methodology

    2.2.1.Time-lagged SVD

    Time-lagged SVD analysis was used to establish the dynamical connection between the geopotential height fields in the stratosphere and troposphere(Perlwitz and Harnik,2003, 2004).The leading coupled modes were extracted from the spatiotemporal structures of the geopotential height perturbations.We performed the analysis and estimated the covariancewiththetemporalseriesoftwoheightfieldsatindividual time lags separately.The height perturbations were arranged such that each column was the time series for a given location.Based on Perlwitz and Harnik(2003,2004),the geopotential height fields could be expanded orthogonally,which can be expressed as

    where H1and H2are the geopotential height fields at time t and t+τ,unand vnare the singular matrices,N is the number of modes,and s2nis the square of nth singular value of the covariance matrix between H1and H2–constructed by taking the covariance between the two expansion coefficients a(t) and b(t+τ).The coupled modes are arranged with increasing n and decreasing covariance.S is the total squared covariance between the two coupled fields:

    In our analysis,the reference height was fixed at 10 hPa and the SVD analysis was performed with the levels descending downward from 10 hPa to 1000 hPa at different lags(τ)in time lags from-30 to+30 days,i.e.,61 time lags.At the reference height,the time span from 1 January to 30 March, i.e.,90 days,remained fixed;whereas,the other levels were shifted temporally by-30 to 30 days with an interval of 1 day.Thus,positive lags indicated that the stratosphere was leading and the troposphere was lagging,and vice versa for negative lags.In order to detect the time lag at which the dynamical relation between H1and H2was maximal,correlationcoefficientsbetweentheleadingcoupledmodes,a1and b1,were computed for each of the 61 SVD analyses.

    To understand the contribution of wave processes and the zonalmeanfield,SVDanalysiswasperformedbetweenzonal mean fields,the eddy field(i.e.,deviation from the zonal mean),and the WN1 and WN2 height fields,separately.Prior to the SVD analysis,we removed the mean seasonal cycle and multiplied the data by the square root of the density and the cosine of latitude along the altitude and latitude,respectively.In order to concentrate on the intra-annual variability and exclude the influence of a trend in the covariance, the annual mean averages of the geopotential height fields were removed.Although we did not use any temporal filtering,strong spatial filtering was applied in the wavenumber domain,both for the H1and H2zonal mean fields.To extract the WN1 and WN2 components from the geopotential height perturbations,we applied the least squares fitting (LSF)method for spectral analysis.This method was used to fit a set of zero mean observations,yi,at times i=1,2,···,N, to the equation given by

    where w is the wavenumber,λirepresents the longitudes,and A and B are the coefficients to be fitted.The individual wave components were then computed using the empirical relation

    where A is the amplitude andφis the phase,estimated by means of LSF analysis.

    In order to estimate a grid size independent measure of S, the mean squared covariance,C,between two grid points of the H1and H2fields can be defined as

    where m1and m2are the grid points of H1and H2fields. Here,we interpolated the height fields to 4.5?in longitude and 3?in latitude and performed the SVD analysis between 30?and 85?N.Hence,m1=m2=1600,corresponding to 80 and 20 longitudinal and latitudinal grid points,respectively.

    2.2.2.Blocking index

    Midlatitude blocking is characterized by local formation ofanomalouseasterlyflowduetotheblockingofthewesterly jet and mass transfer from high-to midlatitudes(Namias and Clapp,1951;Treidl et al.,1981;Barriopedro et al.,2006). In general,the blockings are quasi-stationary patterns that persist for several weeks and have a significant impact on rainfall redistribution and the occurrence of extreme weather events at regional scales.A persistent blocking pattern also induces strong advection of polar air,southward,leading to extreme cold weather in boreal winter months(Nath et al., 2014).For January 2008,Zhou et al.(2009)and Nath et al. (2014)reported an anomalous and persistent blocking pattern in the Urals–Siberia region.The frequency of blocking exceeded the climatological high over 55?–70?E.We computed the blocking index from Tibaldi and Molteni(1990),with the additional criteria proposed by Barriopedro et al.(2006).The 500 hPa geopotential height gradients in the north and south (GHGN and GHGS)(units:gpm/latitude)were simultaneously computed using the following expressions:

    where H(λ,θ)is the 500 hPa geopotential height,δis the shift in latitude,GHGS is the measure of the zonal geostrophic wind component,and GHGN is imposed to exclude the non-blocked flows(Barriopedro et al.,2006).An arbitrary longitude was considered to be blocked if the following conditions were satisfied:

    To identify the potential blocks,a three-day running mean filter was applied at each longitude.

    3.Results and discussion

    3.1.Covariability of the stratospheric and tropospheric height fields

    First,we compared the relative dominance of the zonal mean and the height wave fields for January–February–March(JFM),as obtained from the SVD analysis.The squared covariance between the 10 hPa and various pressure levels(10 to 1000 hPa)at different time lags(-30 to 30 days),for the zonal mean,deviation from the zonal mean, WN1 height,and WN2 height are shown in Fig.1.The composite mean patterns for six winters(2002/2003,2004/2005, 2006/2007,2007/2008,2010/2011 and 2011/2012)when the downward coupling was stronger is shown in the upper panels,Figs.1a–d.The lower panels,Figs.1e–h,exhibit the SVD patterns for the 2008 event,in order to compare the consistencywiththe compositemeanpattern.Forthezonalmean field(Figs.1a and e),in the positive time lag(stratosphere leads),the covariance is stronger and extended(longer time scale)in the lower stratospheric heights.Meanwhile,in the negative time lag(troposphere leads),the covariance is relatively weaker and less persistent below 20 km.In the present analysis,the covariability is maximum and dominant for the leading coupled mode(first),because it explains around 80% of the squared covariance in all height regions.

    Figures 1b and f depict the covariability of the deviation from the zonal mean field for the composite and 2008 case, respectively.Unlike the zonal mean field,the covariance is relatively stronger for the negative and positive time lag in themid-tropospheric(3–11km)andlowerstratospheric(>16 km)heights,respectively.The mid-tropospheric covariance is well extended over all time lags(-30 to 30 days),with a maximum around±15 days;whereas,in the lower stratospheric heights,the covariance is biased towards the positive side,with peaks around 0–5 day(>20 km)both for the composite and 2008 case.The zonal deviation field includes the contribution of various wave processes,and at given latitudes the WN1 and WN2 components were separated out using the LSF method,described in section 2.2.1.

    Fig.1.The covariance(units:gpm2)between the geopotential height fields at 10 hPa and all pressure levels between 1000 and 10 hPa,for time lags ranging from-30 to 30 days:(a–d)the covariance of the zonal mean,deviations from the zonal mean,WN1 height,and WN2 height,respectively,for the composite winter of 2002/2003,2004/2005, 2006/2007,2007/2008,2010/2011 and 2011/2012;(e–h)the same,but for 2008 case.A positive time lag indicates that the stratospheric field is leading.The left and right axes represent the height(units:km)and pressure(units:hPa), respectively.

    As is clear from Fig.1c,the WN1 covariance(~300gpm2)for the composite case compares well with the deviation from the zonal mean field(~400 gpm2).We can see the humps with larger covariance at the lags of-3 days(troposphere leads)and+15 days(stratosphere leads)in the WN1 field.The features are quite consistent with the 2008 case (Fig.1g),at least on the positive side,with stronger downward coupling at+15 days’lag.Another noticeable feature is the intense downward coupling due to the WN1 field,with a persistent covariance pattern down to the surface at+15 days’lag.But,for the WN2 field,both the composite and 2008 event(Figs.1d and h)exhibit much weaker covariance (~100 gpm2)throughout the height range.The WN2 covariance,meanwhile,although weaker than WN1,exhibits a dominant peak at around+5 days’lag(stratosphere leads) and 400 hPa,both for the composite and 2008 event.Perlwitz and Harnik(2004)linked the planetary wave reflection with the humps in the positive lags when the basic state of the stratosphere was reflective.In Nath et al.(2014),it was shown that,apart from the zonal mean reflective index,the longitudinal variation too has a severe impact on regional weather extremes,and the concept of a partially reflective stratospheric background state was introduced.All six winters chosen for the composite analysis–based on Shaw and Perlwitz(2013),Kodera et al.(2013)and Dunn-Sigouin and Shaw(2015)–except 2006/2007,exhibit a partially reflective stratospheric background state.Therefore,during these years,when the downward coupling was stronger,the relative dominance of the WN1 covariance from the stratosphere should have had a significant impact on the tropospheric circulation pattern.

    To further elucidate the relative contribution of the zonal mean and the WN1 field,we plotted the covariance for the zonal mean(850 hPa)and WN1(400 hPa)field,both for the composite(Fig.2a)and 2008 event(Fig.2b).As can be seen,for the WN1 field,the covariance is much higher (lower)for the negative(positive)lag,with a prominent peak at around+15 days when the leading stratosphere coupled strongly with the troposphere.Meanwhile,in the negative lag(troposphere leads),the peaks are prominent at lags of-4 days and-10 days for the composite and 2008 case,respectively.This difference is obvious and can be attributed to the difference in upward wave propagation during the six winters included in the composite analysis.The winters were chosen based on Shaw and Perlwitz(2013),Kodera et al.(2013)and Dunn-Sigouin and Shaw(2015),when downward coupling was prominent,irrespective of any precursory upward wave propagation events.For example,according to Dunn-Sigouin and Shaw(2015,Table 1),there were three downward propagating events,on 6 February,25 February and 31 March 2003,but the upward propagation occurred long before the downward coupling events on 14 January 2003.Similarly, for the 2008 case,there is no upward propagation prior to the wave reflection event.In the zonal mean field,the covariance is much lower compared with the WN1 field,and the absence of any significant maxima,either in the positive or negative time lag,is prominent.

    Fig.2.The covariance[gpm2]between the tropospheric and stratospheric height fields at 30?–85?N for time lags between -30 and 30 days for the composite winters(upper panel)and 2008 case(bottom panel)between the leading coupled mode. Blueline:covariabilitybetween10hPaand850hPazonalmean fields;Black line:covariability between 10 and 400 hPa height WN1 fields.A positive time lag indicates that the stratospheric field is leading.The maxima are significant at least to the 99% confidence level.

    3.2.Evolution of WN1 height and heat flux for the 2008 event Based on the SVD timed lagged analysis,we identified that—both in the composite and 2008 case—downward coupling due to reflection was strongest at+15 days’lag.This corresponds to 15 of January as the key date for the 2008 case study.The evolution of the high-latitude WN1 pattern can be illustrated using a Hovm¨oller plot.Figure 3a shows the total WN1 pattern averaged between 60?and 80?N at 400 hPa (black contours)and 10 hPa(coloring)as a function of longitude and time from-20 to+20 days(15 January as the start date).Downward WN1 coupling events clearly coincide with changes in the tropospheric wave pattern.In the first stage (-20 to-10 days),the 400 hPa WN1 pattern over the Urals–Siberia region is very weak.During the second stage(-10 to 0 days),the amplitude of the 10 hPa high-latitude WN1 pattern reaches a maximum and precedes the maximum amplitude at 400 hPa,which occurs during stage three(0 to10 days).The features are consistent with Shaw and Perlwitz (2013,Fig.5).The amplitude of the 400 hPa WN1 pattern reaches a maximum during the third stage and at the same time the pattern continues to move westward.Finally,in the fourth stage(10 to 20 days),the amplitude of the WN1 pattern decreases significantly and,overall,the wave pattern evolution is very consistent with a downward wave couplingevent:a stratospheric WN1 anomaly at 10 hPa precedes a tropospheric WN1 anomaly at 400 hPa.All the features are highly consistent with Shaw and Perlwitz(2013,Fig.5).

    Fig.3.(a)Evolution of the total 400 hPa(black contours)and 10 hPa(color shading)WN1 pattern averaged from 60?to 80?N for the 2008 event as a function of time from-20 to 20 days and longitude.The contour interval is 10 m,and the ranges are-100 to 100 m and-75 to 75 m for the black contours and color shading,respectively.(b)As in(a)but for the WN1 heat flux anomaly.The contour interval is 0.01 m2s-2and 0.003 m2s-2,and the ranges are from-0.1 to 0.1 m2s-2and-0.018 to 0.018 m2s-2for the black contour and color shading,respectively.The bold dotted and light normal black contours in(a,b)represent the negative and positive anomalies,respectively.The longitude–time section of the blocking index(gray contour lines)is overplotted in both(a)and(b).

    A downward coupling event due to reflection is linked withthetransitionoftheheatfluxanomaly(productofmeridional wind and temperature anomaly)from positive to negative in the stratosphere(Shaw and Perlwitz,2013).Therefore,we also plotted the evolution of the WN1 heat flux anomaly for the 2008 case(Fig.3b).Like the WN1 anomaly, the heat flux pattern averaged between 60?N to 80?N at 400 hPa(black contours)and 10 hPa(coloring)as a function of longitude and time from-20 to+20 days is shown.In the first stage(-20 to-10 days),the heat flux anomaly in the stratosphere is strongly positive,particularly over the Urals–Siberia region,which is indicative of an upward wave coupling precursor.Meanwhile,in the troposphere,there is no significant heat flux anomaly during this stage.In the second stage(-10 to 0 days),the heat flux anomaly in the stratosphere changes sign from positive to negative.In addition,a positive heat flux anomaly starts to develop in the troposphere.In the third stage(0 to 10 days),the negative anomaly in the stratosphere attains its maximum,with subsequent development of a strong positive heat flux anomaly in the troposphere.The tropospheric maxima clearly lag the minima in the stratosphere.Finally,in the fourth stage(10–20 days),the heat flux anomaly weakens,but remains negative in the stratosphere;while in the troposphere,it fades out completely.The features are highly consistent with Shaw and Perlwitz(2013),completely describing the evolution of the WN1 anomaly during the downward wave coupling events.

    3.3.Urals-Siberia blocking index

    To illustrate the zonal propagation of planetary waves, Nath et al.(2014)computed the eddy component of Plumb fluxes(Plumb,1985)in 3D space.From the vertical component of wave activity flux at 200 hPa,they showed that the reflection phenomena were more prominent from 10 to 19 January 2008.Furthermore,the upward propagation was stronger over the Labrador coast and Baffin Island.Whereas, the reflected components were prominent over the Eurasian continent and eastern parts of China.Upward and reflected fluxes were prominent in the upstream and downstream regions of the reflecting surfaces(Nath et al.,2014),respectively,indicating the impact of the polar jet stream,which preferentially guided the planetary waves(WN1)eastward and downward of the source region(Nath et al.,2014).

    Fig.4.Heterogeneous regression pattern(units:gpm)of the leading coupled mode of the(a,b)10 hPa and(c, d)400 hPa WN1 fields at time lags of(a,c)-10 days and(b,d)+15 days for the 2008 case.The color shading varies from-400 to 400 gpm in(a,b)and from-40 to 40 gpm in(c,d),with an interval of 40 gpm and 4 gpm, respectively.These maps were constructed by regressing the time series of the 10 hPa(400 hPa)WN1 fields onto the temporal expansion coefficients of the leading mode of 400 hPa(10 hPa).The percentage in the title of the individual subplots indicates the variance accounting for the leading coupled mode.

    Fig.5.Phase difference(degrees)between the associated WN1 regression patterns at 10 hPa and 400 hPa,averaged over the latitudinal band of 30?–85?N,as a function of time lags.Negative and positive values indicate westward and eastward phase shifts with heights,respectively. The red line indicates zero phase difference.

    Several authors have delineated the role of blocking as a precursor to SSW events(Martius et al.,2009).Recently, Kodera et al.(2013)elucidated the relationship between stratospheric planetary wave reflection and blocking formation in the troposphere during SSW events.The upward propagation of the planetary waves in the pre-warming stage involves a Euro-Atlantic block;whereas,the downward propagation promotes the formation of Pacific blocking during the warming event.The longitude–time section of the blocking index were overplotted(gray contour lines)in both Figs.3a and b to address the coincidence of the blocking event with tropospheric WN1 evolution due to reflection.The blocking index(dimensionless;contour lines)represents the zonal and temporal spread over which the mid-tropospheric flow is blocked,i.e.,for which all three criteria[Eqs.(7a–c)]are satisfied simultaneously.A strong Urals–Siberia blocking is predominant in the third stage,when the amplitude of the 400 hPa WN1 pattern reaches its maximum value and the pattern continues to move westward.Consistently,the negative heat flux anomaly in the stratosphere attains its maximum withsubsequent development of a strong positive anomaly in the troposphere.This feature is analogous to the formation of an anomalous Urals–Siberia blocking high,downstream of the reflected fluxes(Nath et al.,2014)from 15 to 25 January. The blocking index was calculated using Eqs.(7a–d).Furthermore,it is clear from the Eliassen–Palm(EP)flux vector (Nathetal.,2014)that,aroundmid-January,theupwardcomponent of the high latitude wave guide was very weak;it was only the downward component that could have contributed to the development of the blocking high over the Urals–Siberia region.

    Fig.6.Daily 500 hPa geopotential height maps from(a–l)10–21 January 2008,respectively.The contours are from 5100 to 5900 gpm,with an interval of 60 gpm.Heights between 5100 to 5160 gpm are marked with red contours.

    3.4.Heterogeneous regression patterns and tropospheric circulation

    Perlwitz and Graf(2001)showed that the stratospheric WN1 field(50 hPa)leads the tropospheric WN1 field(500 hPa)by 6 days;whereas,using time-lagged SVD analysis, Perlwitz and Harnik(2004)investigated the close relationship between the 500 hPa(tropospheric field)and 50 hPa,30 hPa and 10 hPa(stratospheric field)levels,individually,in composites of all winter seasons.Again,Kodera and Chiba (1995)showed that,in the troposphere,the circulation pattern changed significantly in relation to SSW events during 1984–85.They also linked the generation of anomalous cold surges and synoptic-scale eddies due to changes in the meridional propagation of tropospheric waves around the 500 hPa level.They further suggested that the changes in planetary wave structure could trigger enhanced baroclinic waves in the troposphere.

    In the present study,we expected the heterogeneous regression pattern at dominant positive and negative time lags to exhibit gradual eastward and westward shift,in phase, along the vertical direction,respectively.As discussed in section3.1,strongdownwardcouplingduetotheWN1fieldwith a persistent covariance(between 10 hPa and 400 hPa)pattern down to the surface at lags of-10 days(troposphere leads) and+15 days(stratosphere leads)is noticeable in Fig.2a. The two maxima exceed at least the 99%confidence level. The strength of coupling between the modes at different time lags is well illustrated by the heterogeneous regression patterns of the WN1 field at 400 hPa and 10 hPa,constructed when the 400 hPa(10 hPa)field leads the 10 hPa(400 hPa)at a time lag of-10 days(+15 days).The regression patterns in the lag of-10 days(+15 days),both for 10 hPa and 400 hPa,are plotted in Figs.4a and b(4c and d),respectively,to understand the evolution of the entire process.The regression patterns associated in the negative and positive time lags exhibit completely different structure.At-10 days’lag,the ridge of the WN1 field at 10 hPa shifts westward with respect to the 400 hPa level;whereas,at+15 days’lag,the regression pattern at 400 hPa level shifts eastward relative to the 10 hPa level.This eastward shift in phase of the WN1 fieldat 400 hPa is consistent with the downward reflected wave. Since the square of the correlations represents the variance explained locally,the leading coupled mode accounts for up to~70%and more of the variance in the region of large amplitude.The variance,in terms of percentage,is noted in the titles of the individual subplots.In addition,at 400 hPa and for positive time lags,the phase tilts largely eastward with increasing latitude.We further compared the phase shift(degrees)between the WN1 field at the 10 hPa level and 400 hPa level;the differences at different time lags are plotted in Fig.5.A continuous westward phase shift(negative)of the WN1 ridge,averaged over the latitude band of 30?–85?N,is clearly visible for all negative time lags;whereas,after+10 days’lag,the phase shift tends towards the zero mark,becomes eastward(positive)after+15 days’lag,and westward again from+22 days’lag.

    To further establish the link between anomalous blocking patterns in the Urals–Siberia region and the tropospheric WN1 field due to changes in circulation pattern,the daily march of the 500 hPa geopotential height fields from 10 to 21 January are plotted in Fig.6.The contours from 5100 to 5900 gpm,with an interval of 60 gpm,are shown,and the geopotential height between 5100 and 5160 gpm is marked in red.Up until 12 January,the WN1 pattern in the troposphere is not very clear and is mainly concentrated at high latitudes.But,from 13 January onward,the WN1 pattern starts to develop slowly and exhibit a meridionally elongated pattern stretching from the Bering Sea to the North European plains.The pattern becomes clearer and is fully developed by 15 January.From 16 January,the trough starts to propagate equatorward and eastward over the Urals–Siberia region;and by 20 January,it reaches the Asian landmass,close to the Tibetan plateau.Moreover,it is clear from the EP flux(Nath et al.,2014)that,aroundmid-January,theupwardcomponentof the high latitude wave guide was very weak;it was only the downward component that could contribute to the development of the blocking high over the Urals–Siberia region.Furthermore,this development of the tropospheric WN1 pattern is consistent with the wave reflection from the stratosphere since 11 January.Moreover,the southeastward propagation and the intensification of the trough is consistent with the development of the strong blocking pattern in the Urals–Siberia region from 16 January onward.Although the reflection of the planetary waves ceases by 19 January,the blocking event persists until the last week of January with gradual eastward propagation.This extraordinary persistence had a significant impact on the extreme cold event over Eurasia and parts of China(Nath et al.,2014).

    4.Summary and discussion

    Time-laggedSVDanalysiswasperformedtocomparethe covariance and correlation coefficients for the zonal mean and wave processes in the latitude band of 30?–85?N.We comparedthe2008winterwiththecompositemeanpatternof six winters(2002/2003,2004/2005,2006/2007,2007/2008, 2010/2011and2011/2012)whenthedownwardcouplingwas stronger.Thefeatureswereconsistentinthezonalmeanfield, with stronger and extended(longer time scale)covariance in the lower stratospheric heights(positive lag side).For the WN1 field,the covariance was much higher(lower)for the negative(positive)lag,with a prominent peak around+15 days when the leading stratosphere coupled strongly with the troposphere.In all six winters except 2007,the basic state was partially reflective of the WN1 field;and during 2008, the reflective index was strongly negative over the Atlantic Ocean and Eurasian continent(Nath et al.,2014),favorable for the propagation of planetary waves down to the troposphere.We also analyzed the evolution of the WN1 anomaly and heat flux,both in the troposphere and stratosphere,during JFM 2008.The amplitude of the tropospheric WN1 pattern reached a maximum and was consistent with a downward wave coupling event.A stratospheric WN1 anomaly at 10 hPa preceded a tropospheric WN1 anomaly at 400 hPa. Similarly,the negative heat flux anomaly in the stratosphere attained its maximum,with subsequent development of a strong positive heat flux anomaly in the troposphere.The tropospheric maxima clearly lagged the minima in the stratosphere.

    To interpret the occurrence of the anomalous blocking high over the Urals–Siberia region due to wave reflection, we focused on the period 10–21 January 2008.The blocking anomaly developed strongly when the amplitude of the 400 hPa WN1 pattern reached its maximum value,and the pattern continued to move westward in response to the reflection of planetary waves down to the troposphere.From the zonal-mean EP flux vectors(Nath et al.,2014),on 10 January,the high-latitude wave guide pointed vertically upward. On 13 January,the upward fluxes exhibited a gradual weakening trend,while a strengthening in the downward component was evident.This reflection/overturning was stronger on 16 January and continued until 19 January.In 3D space, wave fluxes propagated upward from the Labrador coast and reflected back to the Eurasian continent.

    To illustrate the strength of the coupling between 10 hPa and the 400 hPa WN1 fields,we plotted the heterogeneous regression pattern at negative and positive time lags.The associated regression patterns for the 10 hPa and 400 hPa levels at negative and positive time lags shifted westward and eastward relative to the 400 hPa and 10 hPa levels,respectively.This eastward phase shift of the WN1 ridge at 400 hPa was consistent with the reflection of the WN1 field from the stratosphere.TherelationshipbetweenthetroposphericWN1 field and the anomalous blocking high in the Urals–Siberia region was further established based on the daily march of the 500 hPa geopotential height fields.From the EP flux it was quite clear that,around mid-January,the upward component of the high-latitude wave guide was very weak;it was only the downward component that could have contributed to the development of the blocking high over the Urals–Siberia region.Moreover,we found that,from 13 January onward,the WN1 pattern started to develop slowly and exhibit a meridionally elongated pattern stretching from the Bering Sea to the North European plains.The pattern developed fully by 15January and,from 16 January onward,the trough started to propagate equatorward and eastward over the Urals–Siberia region.And by 20 January,it reached the Asian landmass, close to the Tibetan Plateau.This development of the tropospheric WN1 pattern was due to the reflection of the WN1 field from the stratosphere from 11 January onwards.

    Acknowledgements.This work was supported jointly by the National Natural Science Foundation of China(Grant Nos. 41350110331 and 41450110431)and the China Postdoctoral Science Foundation(Grant No.2013M541010).

    REFERENCES

    Baldwin,M.P.,and T.J.Dunkerton,1999:Propagation of the Arctic Oscillation from the stratosphere to the troposphere.J. Geophys.Res.,104,30 937–30 946.

    Baldwin,M.P.,and T.J.Dunkerton,2001:Stratospheric harbingers of anomalous weather regimes.Science,294,581–584.

    Barriopedro,D.,R.Garc′?a-Herrera,A.R.Lupo,andE.Hern′andez, 2006:A climatology of Northern Hemisphere Blocking.J. Climate,19,1042–1063.

    Christiansen,B.,2000:A model study of the dynamical connection between the Arctic Oscillation and stratospheric vacillations.J.Geophys.Res.,105,29 461–29 474.

    Christiansen,B.,2001:Downward propagation of zonal mean zonal wind anomalies from the stratosphere to the troposphere:Model and reanalysis.J.Geophys.Res.,106,27 307–27 322.

    Coughlin,K.,and K.K.Tung,2005:Tropospheric wave response to decelerated stratosphere seen as downward propagation in northern annular mode.J.Geophys.Res.,110,D01103,doi: 10.1029/2004JD004661.

    Dee,D.P.,and Coauthors,2011:The ERA-Interim reanalysis: configuration and performance of the data assimilation system.Quart.J.Roy.Meteor.Soc.,137,553–597.

    Dunn-Sigouin,E.,and T.A.Shaw,2015:Comparing and contrasting extreme stratospheric events,including their coupling to the tropospheric circulation.J.Geophys.Res.:Atmos.,120, 1374–1390,doi:10.1002/2014JD022116.

    Geller,M.A.,and J.C.Alpert,1980:Planetary wave coupling between the troposphere and the middle atmosphere as a possible sun-weather mechanism.J.Atmos.Sci.,37,1197–1214.

    Hines,C.O.,1974:A possible mechanism for the production of sun-weather correlations.J.Atmos.Sci.,31,589–591.

    Hui,G.,2009:China’s snow disaster in 2008,who is the principal player?International Journal of Climatology,29,2191–2196,doi:10.1002/joc.1859.

    Kodera,K.,and M.Chiba,1995:Tropospheric circulation changes associated with stratospheric sudden warmings:A case study. J.Geophys.Res.,100,11 055–11 068.

    Kuroda,Y.,and K.Kodera,1999:Role of planetary waves in the stratosphere-troposphere coupled variability in the Northern Hemisphere winter.Geophys.Res.Lett.,26,2375–2378.

    Kodera,K.,H.Mukougawa,and S.Itoh,2008:Tropospheric impact of reflected planetary waves from the stratosphere.Geophys.Res.Lett.,35,L16806,doi:10.1029/2008GL034575.

    Kodera,K.,H.Mukougawa,and A.Fujji,2013:Influence of the vertical and zonal propagation of stratospheric planetary waves on tropospheric blockings.J.Geophys.Res.:Atmos., 118,8333–8345.

    Kodera,K.,K.Yamazaki,M.Chiba,and K.Shibata,1990:Downwardpropagationofupperstratosphericmeanzonalwindperturbation to the troposphere.Geophys.Res.Lett.,17,1263–1266.

    Martius,O.,L.M.Polvani,and H.C.Davies,2009:Blocking precursors to stratospheric sudden warming events.Geophys. Res.Lett.,36,L14806,doi:10.1029/2009GL038776.

    Namias,J.,and P.F.Clapp,1951:Observational studies of general circulation patterns.Compendium of Meteorology,T.F. Malone,Ed.,Amer.Meteor.Soc.,551–568.

    Nath,D.,S.Sridharan,S.Sathishkumar,S.Gurubaran,and W. Chen,2013:Lower stratospheric gravity wave activity over Gadanki(13.5?N,79.2?E)during the stratospheric sudden warming of 2009:Link with potential vorticity intrusion near Indian sector.Journal of Atmospheric and Solar-Terrestrial Physics,94,54–64.

    Nath,D.,W.Chen,L.Wang,and Y.Ma,2014:Planetary wave reflection and its impact on tropospheric cold weather over Asia during January 2008.Adv.Atmos.Sci.,31,851–862,doi: 10.1007/s00376-013-3195-8.

    Perlwitz,J.,and H.F.Graf,2001:Troposphere-stratosphere dynamic coupling under strong and weak polar vortex conditions.Geophys.Res.Lett.,28,271–274.

    Perlwitz,J.,and N.Harnik,2003:Observational evidence of a stratospheric influence on the troposphere by planetary wave reflection.J.Climate,16,3011–3026.

    Perlwitz,J.,and N.Harnik,2004:Downward coupling between the stratosphere and troposphere:The relative roles of wave and zonal mean processes.J.Climate,17,4902–4909.

    Plumb,R.A.,1985:On the three-dimensional propagation of stationary waves.J.Atmos.Sci.,42,217–229.

    Shaw,T.A.,and J.Perlwitz,2013:The life cycle of Northern Hemisphere downward wave coupling between the stratosphere and troposphere,J.Climate,26,1745–1763.

    Shaw,T.A.,J.Perlwitz,and N.Harnik,2010:Downward wave coupling between the stratosphere and troposphere:The importance of meridional wave guiding and comparison with zonal-mean coupling.J.Climate,23,6365–6381.

    Shaw,T.A.,J.Perlwitz,and O.Weiner,2014:Tropospherestratosphere coupling:Links to North Atlantic weather and climate,including their representation in CMIP5 models.J. Geophys.Res.:Atmos.,119,5864–5880.

    Tibaldi,S.,and F.Molteni,1990:On the operational predictability of blocking.Tellus A,42,343–365.

    Treidl,R.A.,E.C.Birch,and P.Sajecki,1981:Blocking action in the Northern Hemisphere:A climatological study.Atmos.-Ocean,19,1–23.

    Uppala,S.M.,D.Dee,S.Kobayashi,P.Berrisford,and A.Simmons,2008:Towards a climate data assimilation system:Status update of ERA-Interim.ECMWF Newsletter,115,12–18.

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

    Nath,D.,and W.Chen,2016:Impact of planetary wave reflection on tropospheric blocking over the Urals–Siberia region in January 2008.Adv.Atmos.Sci.,33(3),309–318,

    10.1007/s00376-015-5052-4.

    24 March 2015;revised 28 July 2015;accepted 20 August 2015)

    ?Wen CHEN

    Email:cw@post.iap.ac.cn

    亚洲av欧美aⅴ国产| 亚洲精品av麻豆狂野| 亚洲人成77777在线视频| av电影中文网址| 精品熟女少妇八av免费久了| 亚洲av日韩精品久久久久久密 | 亚洲av男天堂| 50天的宝宝边吃奶边哭怎么回事| 亚洲专区国产一区二区| 男女无遮挡免费网站观看| 岛国毛片在线播放| 亚洲色图 男人天堂 中文字幕| 一本色道久久久久久精品综合| 久久国产精品人妻蜜桃| 一级a爱视频在线免费观看| 涩涩av久久男人的天堂| 看十八女毛片水多多多| 成人黄色视频免费在线看| 男女国产视频网站| 亚洲欧美一区二区三区国产| 精品一品国产午夜福利视频| 99国产精品一区二区蜜桃av | 日韩制服丝袜自拍偷拍| 黄色一级大片看看| 免费在线观看黄色视频的| 91老司机精品| 久久久久久免费高清国产稀缺| 午夜福利乱码中文字幕| 一级,二级,三级黄色视频| 一本—道久久a久久精品蜜桃钙片| 免费一级毛片在线播放高清视频 | 亚洲一区中文字幕在线| 国产xxxxx性猛交| 精品国产一区二区三区久久久樱花| 国产成人影院久久av| 国产精品三级大全| 9热在线视频观看99| av在线老鸭窝| 大陆偷拍与自拍| 亚洲人成电影免费在线| 中文字幕高清在线视频| 午夜免费成人在线视频| 国产熟女欧美一区二区| 成年女人毛片免费观看观看9 | 久久久久精品人妻al黑| 久9热在线精品视频| 久久久久久久久免费视频了| 热re99久久精品国产66热6| 国产91精品成人一区二区三区 | 色94色欧美一区二区| 别揉我奶头~嗯~啊~动态视频 | 精品国产超薄肉色丝袜足j| 人妻 亚洲 视频| 男女免费视频国产| 夫妻午夜视频| 国产视频一区二区在线看| avwww免费| 久久人妻福利社区极品人妻图片 | avwww免费| 亚洲国产精品一区二区三区在线| 高清av免费在线| 免费在线观看视频国产中文字幕亚洲 | 99久久综合免费| 亚洲国产欧美日韩在线播放| 美女高潮到喷水免费观看| 国产成人影院久久av| 久久久久精品国产欧美久久久 | 丝袜在线中文字幕| 19禁男女啪啪无遮挡网站| 免费观看a级毛片全部| 久久热在线av| 在线亚洲精品国产二区图片欧美| 亚洲av日韩精品久久久久久密 | 久热爱精品视频在线9| 国产免费又黄又爽又色| 久久久久久久久久久久大奶| 久久久久精品国产欧美久久久 | 悠悠久久av| 亚洲精品美女久久av网站| 国产精品一国产av| 免费一级毛片在线播放高清视频 | 超碰成人久久| 美女中出高潮动态图| 永久免费av网站大全| 日韩av免费高清视频| 下体分泌物呈黄色| 久久av网站| 国产精品一国产av| 王馨瑶露胸无遮挡在线观看| 少妇精品久久久久久久| 1024视频免费在线观看| 久久久久久久久免费视频了| 成年av动漫网址| 久久女婷五月综合色啪小说| 各种免费的搞黄视频| 夜夜骑夜夜射夜夜干| 亚洲欧美色中文字幕在线| 亚洲欧美精品自产自拍| 成人影院久久| 纯流量卡能插随身wifi吗| 国产精品.久久久| 日本wwww免费看| 精品少妇内射三级| 免费少妇av软件| 亚洲国产精品一区二区三区在线| 欧美精品av麻豆av| 女人被躁到高潮嗷嗷叫费观| 久久中文字幕一级| 欧美日韩亚洲国产一区二区在线观看 | 新久久久久国产一级毛片| 最近中文字幕2019免费版| 99热网站在线观看| 美女扒开内裤让男人捅视频| 国产黄色视频一区二区在线观看| 国产在线视频一区二区| 国产精品亚洲av一区麻豆| 2021少妇久久久久久久久久久| 亚洲第一av免费看| 亚洲一区二区三区欧美精品| 看免费av毛片| 女人高潮潮喷娇喘18禁视频| 久久久亚洲精品成人影院| 国产高清国产精品国产三级| 久久久久久久国产电影| 精品免费久久久久久久清纯 | 宅男免费午夜| 黄色片一级片一级黄色片| 中文字幕高清在线视频| 男女无遮挡免费网站观看| 日韩一本色道免费dvd| 成人18禁高潮啪啪吃奶动态图| 婷婷色av中文字幕| 亚洲国产精品999| 乱人伦中国视频| 亚洲欧洲国产日韩| 久久九九热精品免费| 精品久久久久久电影网| www.熟女人妻精品国产| av一本久久久久| 国产成人欧美在线观看 | 亚洲成人手机| 久热爱精品视频在线9| 免费黄频网站在线观看国产| 国产午夜精品一二区理论片| 亚洲欧美一区二区三区黑人| a级片在线免费高清观看视频| www.av在线官网国产| 久久九九热精品免费| 国产高清videossex| 欧美大码av| 亚洲精品在线美女| 波多野结衣av一区二区av| 看免费成人av毛片| 别揉我奶头~嗯~啊~动态视频 | 国产人伦9x9x在线观看| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 岛国毛片在线播放| 美女大奶头黄色视频| 日本猛色少妇xxxxx猛交久久| 亚洲欧美精品综合一区二区三区| 另类亚洲欧美激情| 免费看十八禁软件| 久久国产精品人妻蜜桃| 色视频在线一区二区三区| 国产视频一区二区在线看| 国产亚洲av片在线观看秒播厂| 日本a在线网址| 久久久国产一区二区| 亚洲专区国产一区二区| 少妇精品久久久久久久| 最近中文字幕2019免费版| 一区在线观看完整版| 一区二区三区精品91| 一级毛片电影观看| 男人爽女人下面视频在线观看| 国产日韩欧美在线精品| 欧美日本中文国产一区发布| 国产欧美日韩精品亚洲av| 国产熟女欧美一区二区| 啦啦啦在线观看免费高清www| 亚洲自偷自拍图片 自拍| 中文字幕制服av| 国产成人系列免费观看| 在线精品无人区一区二区三| 亚洲一码二码三码区别大吗| 69精品国产乱码久久久| 免费不卡黄色视频| 国产精品三级大全| 中文字幕精品免费在线观看视频| 免费在线观看影片大全网站 | 国产女主播在线喷水免费视频网站| 亚洲国产av影院在线观看| 国产成人av激情在线播放| 亚洲欧美日韩另类电影网站| 久热这里只有精品99| 纵有疾风起免费观看全集完整版| 十分钟在线观看高清视频www| 日日爽夜夜爽网站| 在线亚洲精品国产二区图片欧美| 99国产精品免费福利视频| 久久99精品国语久久久| av天堂久久9| 国产视频一区二区在线看| 亚洲精品久久成人aⅴ小说| 男人操女人黄网站| 美国免费a级毛片| 午夜福利一区二区在线看| 视频区图区小说| 777久久人妻少妇嫩草av网站| www.av在线官网国产| 大片电影免费在线观看免费| 亚洲伊人色综图| 1024视频免费在线观看| 精品国产国语对白av| 少妇被粗大的猛进出69影院| 免费不卡黄色视频| 久久久久久久大尺度免费视频| 亚洲av日韩在线播放| 午夜免费鲁丝| 久久久欧美国产精品| 观看av在线不卡| 女性被躁到高潮视频| 日韩,欧美,国产一区二区三区| av又黄又爽大尺度在线免费看| 亚洲精品国产区一区二| 亚洲av欧美aⅴ国产| 校园人妻丝袜中文字幕| 亚洲人成网站在线观看播放| 亚洲国产欧美在线一区| 纯流量卡能插随身wifi吗| 母亲3免费完整高清在线观看| 真人做人爱边吃奶动态| 国产精品二区激情视频| 免费av中文字幕在线| 亚洲精品一区蜜桃| 99香蕉大伊视频| 亚洲午夜精品一区,二区,三区| 欧美激情极品国产一区二区三区| 黄色毛片三级朝国网站| 国产深夜福利视频在线观看| 亚洲人成电影观看| 99久久99久久久精品蜜桃| 国产午夜精品一二区理论片| 亚洲国产欧美在线一区| 亚洲精品国产av蜜桃| 日韩人妻精品一区2区三区| 91精品伊人久久大香线蕉| 久久天堂一区二区三区四区| 亚洲综合色网址| 91字幕亚洲| 欧美日韩成人在线一区二区| 97人妻天天添夜夜摸| 老司机影院毛片| 久久综合国产亚洲精品| 亚洲国产看品久久| 18禁裸乳无遮挡动漫免费视频| 久久久久久免费高清国产稀缺| 男女下面插进去视频免费观看| 亚洲中文日韩欧美视频| 成年人免费黄色播放视频| 一级毛片电影观看| 日韩视频在线欧美| 国产成人免费观看mmmm| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区大全| 欧美性长视频在线观看| 国产高清不卡午夜福利| 免费在线观看日本一区| 欧美黑人欧美精品刺激| 青春草亚洲视频在线观看| 亚洲欧美色中文字幕在线| bbb黄色大片| 一区二区三区激情视频| 日韩一本色道免费dvd| 人妻人人澡人人爽人人| 啦啦啦 在线观看视频| 久久久久精品国产欧美久久久 | 亚洲av美国av| 一级毛片电影观看| av欧美777| 欧美黄色淫秽网站| 午夜免费成人在线视频| 一区二区av电影网| 国产欧美日韩一区二区三区在线| 亚洲色图 男人天堂 中文字幕| 亚洲欧美一区二区三区国产| 欧美在线一区亚洲| 国产老妇伦熟女老妇高清| 夫妻性生交免费视频一级片| 在线观看免费高清a一片| 欧美老熟妇乱子伦牲交| 国产视频首页在线观看| 午夜福利一区二区在线看| 可以免费在线观看a视频的电影网站| 成年av动漫网址| 中文欧美无线码| 97人妻天天添夜夜摸| 狠狠婷婷综合久久久久久88av| 久久久久久久精品精品| 中文字幕人妻丝袜一区二区| 在线观看国产h片| 久久狼人影院| 欧美黄色淫秽网站| 国产亚洲精品第一综合不卡| 久久久久国产精品人妻一区二区| 一本一本久久a久久精品综合妖精| 午夜福利乱码中文字幕| 国产爽快片一区二区三区| 国产亚洲精品久久久久5区| 大码成人一级视频| 久久久久国产精品人妻一区二区| 国产亚洲欧美在线一区二区| 亚洲欧洲日产国产| 精品国产超薄肉色丝袜足j| 国产欧美亚洲国产| 桃花免费在线播放| 国产精品久久久久久精品电影小说| 巨乳人妻的诱惑在线观看| 亚洲一区二区三区欧美精品| 丰满迷人的少妇在线观看| 纯流量卡能插随身wifi吗| 欧美日韩亚洲综合一区二区三区_| 两人在一起打扑克的视频| 咕卡用的链子| av在线app专区| 波野结衣二区三区在线| 国产精品人妻久久久影院| 色婷婷久久久亚洲欧美| 精品一区二区三卡| 国产色视频综合| 大片免费播放器 马上看| 亚洲精品久久成人aⅴ小说| 各种免费的搞黄视频| 女人高潮潮喷娇喘18禁视频| 国产成人91sexporn| 中文字幕亚洲精品专区| bbb黄色大片| 欧美日韩亚洲高清精品| 日本一区二区免费在线视频| 久久久久久久国产电影| 最黄视频免费看| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 秋霞在线观看毛片| 欧美精品一区二区大全| 欧美精品啪啪一区二区三区 | 成人免费观看视频高清| 一区二区三区四区激情视频| 亚洲中文av在线| 国精品久久久久久国模美| 亚洲精品美女久久av网站| 50天的宝宝边吃奶边哭怎么回事| 老汉色av国产亚洲站长工具| 国产成人免费无遮挡视频| 国产成人一区二区在线| 亚洲欧美日韩高清在线视频 | 免费黄频网站在线观看国产| 狂野欧美激情性xxxx| 国产在线观看jvid| 美女扒开内裤让男人捅视频| av视频免费观看在线观看| 一边亲一边摸免费视频| 日本欧美视频一区| 日本猛色少妇xxxxx猛交久久| 日韩 欧美 亚洲 中文字幕| 各种免费的搞黄视频| 丝袜喷水一区| 精品久久久久久电影网| 免费黄频网站在线观看国产| 久久久久视频综合| 国产精品香港三级国产av潘金莲 | 亚洲av美国av| 日本av手机在线免费观看| 深夜精品福利| 亚洲精品av麻豆狂野| 性高湖久久久久久久久免费观看| 久久国产亚洲av麻豆专区| 欧美激情极品国产一区二区三区| 日日夜夜操网爽| 久久鲁丝午夜福利片| 高清不卡的av网站| 看十八女毛片水多多多| 精品久久久精品久久久| 国产高清不卡午夜福利| 日日摸夜夜添夜夜爱| 日韩熟女老妇一区二区性免费视频| 亚洲中文日韩欧美视频| 肉色欧美久久久久久久蜜桃| 美女国产高潮福利片在线看| 免费看十八禁软件| 亚洲精品国产区一区二| 午夜老司机福利片| 涩涩av久久男人的天堂| 在线av久久热| 在线观看人妻少妇| 日韩av免费高清视频| 久热这里只有精品99| 欧美国产精品一级二级三级| 免费在线观看完整版高清| 人人澡人人妻人| 成年人免费黄色播放视频| 国产成人免费无遮挡视频| 久久久国产精品麻豆| 欧美日韩国产mv在线观看视频| 中文字幕另类日韩欧美亚洲嫩草| 热99国产精品久久久久久7| 国产女主播在线喷水免费视频网站| 下体分泌物呈黄色| 丝袜人妻中文字幕| 丰满饥渴人妻一区二区三| 久久性视频一级片| 999久久久国产精品视频| 老司机深夜福利视频在线观看 | 在线观看免费午夜福利视频| 日韩制服丝袜自拍偷拍| 久久久久国产一级毛片高清牌| 男女无遮挡免费网站观看| 欧美另类一区| 精品少妇黑人巨大在线播放| 91精品国产国语对白视频| 久久国产精品大桥未久av| 中文字幕精品免费在线观看视频| 久久免费观看电影| 国产国语露脸激情在线看| 香蕉国产在线看| av电影中文网址| 亚洲精品美女久久av网站| 少妇的丰满在线观看| 婷婷成人精品国产| 纵有疾风起免费观看全集完整版| 亚洲中文日韩欧美视频| 久久av网站| 2018国产大陆天天弄谢| av在线app专区| 91国产中文字幕| 欧美日韩精品网址| 欧美激情高清一区二区三区| 丁香六月欧美| 尾随美女入室| √禁漫天堂资源中文www| 日本91视频免费播放| 又大又黄又爽视频免费| 性色av一级| 一级片免费观看大全| av不卡在线播放| 久久女婷五月综合色啪小说| 午夜日韩欧美国产| 国产女主播在线喷水免费视频网站| 亚洲精品第二区| 成人黄色视频免费在线看| 乱人伦中国视频| 成人国产av品久久久| 两个人看的免费小视频| 色94色欧美一区二区| 超碰成人久久| 国产精品 国内视频| 蜜桃国产av成人99| 桃花免费在线播放| 另类亚洲欧美激情| 久久青草综合色| 99精国产麻豆久久婷婷| 亚洲国产精品成人久久小说| 两性夫妻黄色片| 黑人巨大精品欧美一区二区蜜桃| 国产成人欧美在线观看 | av在线播放精品| 精品国产一区二区三区久久久樱花| 中文乱码字字幕精品一区二区三区| 丝袜人妻中文字幕| 成人国产一区最新在线观看 | 欧美老熟妇乱子伦牲交| 精品高清国产在线一区| 纵有疾风起免费观看全集完整版| 热re99久久精品国产66热6| 午夜av观看不卡| 久久久久久免费高清国产稀缺| 大香蕉久久成人网| 亚洲,一卡二卡三卡| 成人18禁高潮啪啪吃奶动态图| 老熟女久久久| 天天躁夜夜躁狠狠久久av| 精品第一国产精品| 在线精品无人区一区二区三| 久久影院123| 国产欧美日韩一区二区三 | 免费在线观看日本一区| 国产精品久久久人人做人人爽| svipshipincom国产片| 考比视频在线观看| 成人18禁高潮啪啪吃奶动态图| 午夜福利,免费看| 一级毛片女人18水好多 | 久久99热这里只频精品6学生| 国产在线免费精品| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久人人做人人爽| 亚洲精品国产色婷婷电影| 午夜免费男女啪啪视频观看| 亚洲欧洲日产国产| av在线老鸭窝| 成人国产av品久久久| 免费观看a级毛片全部| 久久久久久免费高清国产稀缺| 在线观看免费日韩欧美大片| 亚洲成人免费电影在线观看 | 侵犯人妻中文字幕一二三四区| 男女免费视频国产| 亚洲一区中文字幕在线| 天天操日日干夜夜撸| 97人妻天天添夜夜摸| 亚洲激情五月婷婷啪啪| 黄网站色视频无遮挡免费观看| 久久99精品国语久久久| 精品国产超薄肉色丝袜足j| 亚洲五月婷婷丁香| 亚洲欧美激情在线| 只有这里有精品99| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| 国产精品久久久久久精品古装| 国产精品偷伦视频观看了| 丰满饥渴人妻一区二区三| 久久精品亚洲av国产电影网| 在线观看www视频免费| 国产色视频综合| 美女视频免费永久观看网站| 如日韩欧美国产精品一区二区三区| 婷婷成人精品国产| 中文字幕人妻熟女乱码| 91字幕亚洲| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 最新在线观看一区二区三区 | 一级毛片黄色毛片免费观看视频| 夜夜骑夜夜射夜夜干| 最近最新中文字幕大全免费视频 | 黄色视频在线播放观看不卡| 亚洲av日韩精品久久久久久密 | 乱人伦中国视频| 国产成人免费观看mmmm| 国产真人三级小视频在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲精品久久成人aⅴ小说| 大片电影免费在线观看免费| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品久久久久久噜噜老黄| 韩国精品一区二区三区| 极品人妻少妇av视频| 国产在线一区二区三区精| 成人黄色视频免费在线看| 亚洲精品中文字幕在线视频| 亚洲国产欧美一区二区综合| 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 晚上一个人看的免费电影| 国产成人av教育| av线在线观看网站| 少妇的丰满在线观看| 一区福利在线观看| 丰满迷人的少妇在线观看| 19禁男女啪啪无遮挡网站| 精品卡一卡二卡四卡免费| 黄频高清免费视频| 日韩 亚洲 欧美在线| 欧美日韩成人在线一区二区| 飞空精品影院首页| 日韩制服骚丝袜av| 欧美在线黄色| 久久ye,这里只有精品| 成年人免费黄色播放视频| 美女脱内裤让男人舔精品视频| 欧美日韩精品网址| 菩萨蛮人人尽说江南好唐韦庄| 黄片小视频在线播放| 2021少妇久久久久久久久久久| 日本色播在线视频| 欧美黑人精品巨大| netflix在线观看网站| 91国产中文字幕| 欧美在线黄色| 国产成人a∨麻豆精品| 1024视频免费在线观看| 精品亚洲成国产av| 一级毛片女人18水好多 | 丝袜在线中文字幕| 99久久精品国产亚洲精品| 美女午夜性视频免费| 制服人妻中文乱码| 啦啦啦 在线观看视频| 一本大道久久a久久精品| 亚洲精品日韩在线中文字幕| 18禁观看日本| 高清av免费在线| 宅男免费午夜| 久久 成人 亚洲| 免费观看av网站的网址| 老司机亚洲免费影院| 十八禁人妻一区二区| 乱人伦中国视频| 欧美国产精品一级二级三级| 国产精品九九99| 亚洲一区中文字幕在线| 免费观看a级毛片全部| 自线自在国产av| 欧美国产精品一级二级三级| www.av在线官网国产| 热99国产精品久久久久久7| 久久99精品国语久久久| 午夜免费成人在线视频| 国产成人啪精品午夜网站| 久久久久精品人妻al黑| 国产在视频线精品| 男女边摸边吃奶| 欧美 日韩 精品 国产|