• <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电影| 久久久精品94久久精品| videosex国产| 日韩视频在线欧美| av不卡在线播放| 一级片'在线观看视频| 波野结衣二区三区在线| 亚洲,欧美,日韩| 一区二区三区激情视频| 亚洲精品国产av成人精品| 色播在线永久视频| 最近中文字幕高清免费大全6| 国产精品女同一区二区软件| 人人澡人人妻人| 一区二区三区四区激情视频| 最近中文字幕高清免费大全6| 一区二区三区四区激情视频| 卡戴珊不雅视频在线播放| 啦啦啦在线免费观看视频4| 日本爱情动作片www.在线观看| 国产精品熟女久久久久浪| 亚洲国产欧美日韩在线播放| 在线观看人妻少妇| 19禁男女啪啪无遮挡网站| 男女边摸边吃奶| 十分钟在线观看高清视频www| 国产视频首页在线观看| 亚洲欧美清纯卡通| 午夜精品国产一区二区电影| 久热这里只有精品99| 国产欧美亚洲国产| 中文字幕人妻熟女乱码| 亚洲精品国产色婷婷电影| 人妻人人澡人人爽人人| 极品少妇高潮喷水抽搐| 亚洲成av片中文字幕在线观看| 亚洲av电影在线进入| 国产一区二区三区av在线| e午夜精品久久久久久久| 久久人人爽av亚洲精品天堂| 性色av一级| 国产色婷婷99| 国产伦人伦偷精品视频| 中国三级夫妇交换| 一级毛片我不卡| 亚洲专区中文字幕在线 | 久久免费观看电影| 51午夜福利影视在线观看| 国产精品香港三级国产av潘金莲 | 国产成人啪精品午夜网站| 一边亲一边摸免费视频| 日日撸夜夜添| 国产成人免费观看mmmm| 别揉我奶头~嗯~啊~动态视频 | 美女视频免费永久观看网站| 性色av一级| 久热爱精品视频在线9| 9191精品国产免费久久| 久久精品亚洲熟妇少妇任你| 熟女少妇亚洲综合色aaa.| 亚洲精品一二三| 久久久国产欧美日韩av| 日本91视频免费播放| 建设人人有责人人尽责人人享有的| 亚洲成人手机| 丰满乱子伦码专区| 街头女战士在线观看网站| 久久国产精品大桥未久av| 黑人巨大精品欧美一区二区蜜桃| 中文字幕精品免费在线观看视频| 国产成人啪精品午夜网站| 91国产中文字幕| 精品酒店卫生间| 美女高潮到喷水免费观看| 欧美精品高潮呻吟av久久| 国产亚洲精品第一综合不卡| 欧美成人精品欧美一级黄| www日本在线高清视频| 中文字幕高清在线视频| 日日撸夜夜添| 一级,二级,三级黄色视频| 欧美精品人与动牲交sv欧美| 久久99一区二区三区| 狠狠婷婷综合久久久久久88av| 久久毛片免费看一区二区三区| 热re99久久精品国产66热6| 看免费av毛片| 一本大道久久a久久精品| 在线观看www视频免费| 纯流量卡能插随身wifi吗| 久久久久久人妻| 精品一区在线观看国产| 丰满乱子伦码专区| 一二三四中文在线观看免费高清| av在线播放精品| 国产人伦9x9x在线观看| 国产视频首页在线观看| 亚洲国产av影院在线观看| 国产伦人伦偷精品视频| 成年美女黄网站色视频大全免费| 欧美精品一区二区大全| 男女床上黄色一级片免费看| 看免费成人av毛片| 久久婷婷青草| 久久av网站| 日韩熟女老妇一区二区性免费视频| 国产精品av久久久久免费| 天天操日日干夜夜撸| 99精品久久久久人妻精品| 七月丁香在线播放| 日本欧美国产在线视频| 国产不卡av网站在线观看| 亚洲一区二区三区欧美精品| 人人妻人人添人人爽欧美一区卜| 新久久久久国产一级毛片| 午夜福利免费观看在线| 少妇人妻 视频| 91精品伊人久久大香线蕉| 伊人久久大香线蕉亚洲五| 亚洲av电影在线进入| 99热网站在线观看| 老司机深夜福利视频在线观看 | 国产精品一区二区精品视频观看| 日韩电影二区| 亚洲人成电影观看| 色视频在线一区二区三区| 在线亚洲精品国产二区图片欧美| 久久久久久久久久久免费av| av片东京热男人的天堂| 亚洲欧美中文字幕日韩二区| 国产在线视频一区二区| 各种免费的搞黄视频| 一级毛片我不卡| 国产又爽黄色视频| 91国产中文字幕| 久久精品久久精品一区二区三区| 丁香六月天网| 欧美精品人与动牲交sv欧美| 日韩成人av中文字幕在线观看| 午夜精品国产一区二区电影| 免费观看av网站的网址| 多毛熟女@视频| 日韩中文字幕欧美一区二区 | 日韩精品有码人妻一区| 好男人视频免费观看在线| 中文字幕亚洲精品专区| 亚洲人成77777在线视频| 亚洲av成人精品一二三区| 亚洲国产最新在线播放| 少妇被粗大的猛进出69影院| 免费观看a级毛片全部| 国产av国产精品国产| 18禁观看日本| 亚洲欧美一区二区三区黑人| 国产亚洲最大av| 少妇 在线观看| 18禁国产床啪视频网站| 看非洲黑人一级黄片| netflix在线观看网站| 人人妻,人人澡人人爽秒播 | 亚洲美女黄色视频免费看| 国产成人av激情在线播放| 丝袜在线中文字幕| 久久人妻熟女aⅴ| 一级片'在线观看视频| 在线免费观看不下载黄p国产| 欧美人与善性xxx| 天天躁夜夜躁狠狠久久av| 欧美精品人与动牲交sv欧美| 女人高潮潮喷娇喘18禁视频| 亚洲国产毛片av蜜桃av| 一区二区日韩欧美中文字幕| 国产野战对白在线观看| 久久久久久人妻| 午夜福利一区二区在线看| 99国产精品免费福利视频| 欧美精品一区二区免费开放| 久久久久久久精品精品| 少妇 在线观看| 国产精品久久久av美女十八| 中文字幕人妻丝袜制服| 精品午夜福利在线看| xxxhd国产人妻xxx| 中文字幕高清在线视频| 色婷婷久久久亚洲欧美| 另类精品久久| 赤兔流量卡办理| 欧美激情高清一区二区三区 | 婷婷色综合www| 女的被弄到高潮叫床怎么办| 宅男免费午夜| 大话2 男鬼变身卡| 操出白浆在线播放| 男女边吃奶边做爰视频| av国产久精品久网站免费入址| 啦啦啦啦在线视频资源| 十八禁网站网址无遮挡| 日韩大片免费观看网站| 在线观看www视频免费| 黄色一级大片看看| 高清av免费在线| 免费看av在线观看网站| 精品国产一区二区三区久久久樱花| 亚洲成国产人片在线观看| 最近手机中文字幕大全| av视频免费观看在线观看| 男女下面插进去视频免费观看| 美女扒开内裤让男人捅视频| 性色av一级| 国产成人精品福利久久| 丰满乱子伦码专区| 午夜精品国产一区二区电影| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久电影网| 女人精品久久久久毛片| 99九九在线精品视频| 飞空精品影院首页| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲国产日韩| 九色亚洲精品在线播放| 国产免费现黄频在线看| 在线观看国产h片| av线在线观看网站| av在线观看视频网站免费| av国产久精品久网站免费入址| 巨乳人妻的诱惑在线观看| 亚洲精品,欧美精品| 久久99热这里只频精品6学生| 老熟女久久久| 在线看a的网站| 亚洲欧洲精品一区二区精品久久久 | av一本久久久久| 老司机靠b影院| 欧美日韩综合久久久久久| 国产亚洲欧美精品永久| 另类精品久久| 国产成人精品久久二区二区91 | 国产亚洲av片在线观看秒播厂| 亚洲一区二区三区欧美精品| 高清欧美精品videossex| 久久久欧美国产精品| 操出白浆在线播放| 国产精品99久久99久久久不卡 | 精品卡一卡二卡四卡免费| 两个人免费观看高清视频| 国产欧美日韩综合在线一区二区| 久热爱精品视频在线9| 国产成人精品在线电影| 亚洲国产精品一区三区| 国产成人系列免费观看| 国产欧美日韩综合在线一区二区| 国产日韩欧美视频二区| 男女之事视频高清在线观看 | 丰满迷人的少妇在线观看| 涩涩av久久男人的天堂| 国产精品一国产av| 九色亚洲精品在线播放| 精品国产露脸久久av麻豆| 久久午夜综合久久蜜桃| 久久狼人影院| 国产成人91sexporn| 最近的中文字幕免费完整| 老汉色av国产亚洲站长工具| 大片电影免费在线观看免费| 日韩制服丝袜自拍偷拍| 日韩大码丰满熟妇| 欧美国产精品va在线观看不卡| 91精品三级在线观看| 日韩大片免费观看网站| 中文字幕制服av| 永久免费av网站大全| 久久 成人 亚洲| 午夜免费观看性视频| 七月丁香在线播放| 自线自在国产av| 亚洲av国产av综合av卡| 亚洲欧美清纯卡通| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 国产xxxxx性猛交| 久久久久精品久久久久真实原创| 超碰97精品在线观看| 亚洲一码二码三码区别大吗| 欧美人与性动交α欧美软件| 老司机影院毛片| 青青草视频在线视频观看| 天天影视国产精品| 日韩伦理黄色片| 九色亚洲精品在线播放| 成年美女黄网站色视频大全免费| 亚洲国产欧美网| 亚洲男人天堂网一区| 最近最新中文字幕免费大全7| 黄色 视频免费看| 母亲3免费完整高清在线观看| 亚洲国产欧美网| 18禁国产床啪视频网站| 高清av免费在线| 欧美亚洲日本最大视频资源| 国产不卡av网站在线观看| 国产老妇伦熟女老妇高清| 啦啦啦中文免费视频观看日本| 久久久久网色| 高清视频免费观看一区二区| 一级毛片电影观看| 欧美在线黄色| 七月丁香在线播放| 乱人伦中国视频| av在线观看视频网站免费| 日本欧美国产在线视频| 国产av码专区亚洲av| 国产免费视频播放在线视频| 亚洲成人国产一区在线观看 | 亚洲熟女精品中文字幕| 青春草国产在线视频| 亚洲国产精品一区二区三区在线| 精品久久久精品久久久| 成人免费观看视频高清| 黄色怎么调成土黄色| 一区福利在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美日韩在线播放| 欧美精品一区二区免费开放| 超碰97精品在线观看| 老鸭窝网址在线观看| 咕卡用的链子| 精品国产一区二区三区久久久樱花| 国产99久久九九免费精品| 亚洲欧美成人精品一区二区| 成人18禁高潮啪啪吃奶动态图| 久久精品熟女亚洲av麻豆精品| 亚洲免费av在线视频| 国产欧美日韩一区二区三区在线| 日本午夜av视频| 国产亚洲av高清不卡| 国产成人精品无人区| 精品国产国语对白av| 久久久久久久久久久免费av| 久久韩国三级中文字幕| 国产精品.久久久| 精品少妇久久久久久888优播| 宅男免费午夜| 久久影院123| 午夜福利乱码中文字幕| 日本猛色少妇xxxxx猛交久久| 国产精品麻豆人妻色哟哟久久| 日本wwww免费看| 亚洲精品国产区一区二| 亚洲精品,欧美精品| 秋霞在线观看毛片| 丝袜人妻中文字幕| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看| 亚洲自偷自拍图片 自拍| 亚洲av日韩精品久久久久久密 | 精品国产国语对白av| 岛国毛片在线播放| 中文字幕av电影在线播放| 久久av网站| 国产精品三级大全| 久久久精品免费免费高清| 又粗又硬又长又爽又黄的视频| 男女午夜视频在线观看| 极品人妻少妇av视频| 男人舔女人的私密视频| 亚洲欧美成人综合另类久久久| 国产精品蜜桃在线观看| 看非洲黑人一级黄片| 女人被躁到高潮嗷嗷叫费观| 欧美人与性动交α欧美精品济南到| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 性色av一级| 亚洲熟女精品中文字幕| 少妇 在线观看| 伊人亚洲综合成人网| 欧美日韩精品网址| 亚洲激情五月婷婷啪啪| 99精国产麻豆久久婷婷| 咕卡用的链子| 欧美日韩综合久久久久久| 99久久综合免费| 精品酒店卫生间| 卡戴珊不雅视频在线播放| 九草在线视频观看| 热re99久久国产66热| 国产成人精品在线电影| 亚洲精品国产av成人精品| 最近中文字幕2019免费版| 在线观看免费高清a一片| www.av在线官网国产| 午夜福利视频精品| 黑丝袜美女国产一区| 美国免费a级毛片| 叶爱在线成人免费视频播放| 欧美 亚洲 国产 日韩一| 亚洲七黄色美女视频| 99久国产av精品国产电影| 精品免费久久久久久久清纯 | 日韩一区二区视频免费看| 多毛熟女@视频| 午夜免费男女啪啪视频观看| 青青草视频在线视频观看| 亚洲国产日韩一区二区| netflix在线观看网站| 国产有黄有色有爽视频| 中文乱码字字幕精品一区二区三区| 国产福利在线免费观看视频| 日韩 亚洲 欧美在线| 超色免费av| 丰满饥渴人妻一区二区三| 亚洲色图 男人天堂 中文字幕| 亚洲精品中文字幕在线视频| 在线观看免费午夜福利视频| 乱人伦中国视频| 蜜桃国产av成人99| 美女主播在线视频| 亚洲少妇的诱惑av| 大码成人一级视频| 午夜免费男女啪啪视频观看| 亚洲精品久久午夜乱码| 国产成人精品在线电影| 99久国产av精品国产电影| 国产女主播在线喷水免费视频网站| 国产黄色视频一区二区在线观看| 黄片无遮挡物在线观看| 亚洲专区中文字幕在线 | 国产成人午夜福利电影在线观看| 天堂8中文在线网| 亚洲国产日韩一区二区| 一区二区三区乱码不卡18| 制服丝袜香蕉在线| 国产黄色视频一区二区在线观看| 日本av免费视频播放| 国产精品蜜桃在线观看| e午夜精品久久久久久久| 国产精品麻豆人妻色哟哟久久| 又粗又硬又长又爽又黄的视频| 国产成人午夜福利电影在线观看| 亚洲 欧美一区二区三区| 人人妻人人澡人人爽人人夜夜| 午夜福利影视在线免费观看| 永久免费av网站大全| 在线免费观看不下载黄p国产| 成人手机av| 中文字幕制服av| 操美女的视频在线观看| 国产又爽黄色视频| 各种免费的搞黄视频| 亚洲国产精品国产精品| 狠狠婷婷综合久久久久久88av| 亚洲av日韩精品久久久久久密 | 午夜日韩欧美国产| av视频免费观看在线观看| 国产精品人妻久久久影院| 另类精品久久| 国产精品国产三级专区第一集| 一级黄片播放器| 成人毛片60女人毛片免费| 色播在线永久视频| 亚洲成人国产一区在线观看 | 男女国产视频网站| 国产精品免费视频内射| 亚洲欧美成人精品一区二区| 亚洲国产精品一区三区| 亚洲成av片中文字幕在线观看| 我的亚洲天堂| 久久久久精品人妻al黑| 久久综合国产亚洲精品| 久久久欧美国产精品| 亚洲国产精品一区三区| a级片在线免费高清观看视频| www.熟女人妻精品国产| 国产免费又黄又爽又色| e午夜精品久久久久久久| 晚上一个人看的免费电影| 国产精品无大码| 超碰97精品在线观看| 国产黄色免费在线视频| 免费黄色在线免费观看| 久久久久久久久久久免费av| 亚洲精品国产色婷婷电影| 在线 av 中文字幕| 丰满饥渴人妻一区二区三| 国产亚洲精品第一综合不卡| 亚洲久久久国产精品| 十分钟在线观看高清视频www| 国产视频首页在线观看| 91精品国产国语对白视频| 欧美精品高潮呻吟av久久| 自拍欧美九色日韩亚洲蝌蚪91| 国产无遮挡羞羞视频在线观看| 在线天堂最新版资源| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品第二区| 亚洲欧美一区二区三区黑人| 美女国产高潮福利片在线看| 飞空精品影院首页| 国产精品秋霞免费鲁丝片| 高清在线视频一区二区三区| 丝袜人妻中文字幕| 男人爽女人下面视频在线观看| www.av在线官网国产| 黄色怎么调成土黄色| 久久久久久久大尺度免费视频| 国产免费现黄频在线看| 女人精品久久久久毛片| 成人午夜精彩视频在线观看| 久久久国产精品麻豆| 日日爽夜夜爽网站| 搡老岳熟女国产| 久久天堂一区二区三区四区| 国产人伦9x9x在线观看| 亚洲美女搞黄在线观看| 免费人妻精品一区二区三区视频| bbb黄色大片| 青春草国产在线视频| 日韩精品免费视频一区二区三区| 女人久久www免费人成看片| 丰满乱子伦码专区| 80岁老熟妇乱子伦牲交| 亚洲欧美成人综合另类久久久| 国精品久久久久久国模美| 国产av国产精品国产| 老司机亚洲免费影院| 亚洲成av片中文字幕在线观看| 一级毛片电影观看| 久久天堂一区二区三区四区| 亚洲精品乱久久久久久| 在线亚洲精品国产二区图片欧美| 欧美亚洲 丝袜 人妻 在线| 国产高清国产精品国产三级| 亚洲成人一二三区av| 国产精品久久久久久人妻精品电影 | 久久久久国产精品人妻一区二区| 亚洲一码二码三码区别大吗| av卡一久久| 岛国毛片在线播放| 熟女少妇亚洲综合色aaa.| 精品国产一区二区三区四区第35| 欧美日韩一级在线毛片| 18禁动态无遮挡网站| 99国产综合亚洲精品| 咕卡用的链子| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| 日本91视频免费播放| 日韩伦理黄色片| 免费黄色在线免费观看| 国产熟女欧美一区二区| 午夜91福利影院| 亚洲在久久综合| 久久久久精品国产欧美久久久 | 久久人人97超碰香蕉20202| 男女免费视频国产| 国产老妇伦熟女老妇高清| 亚洲av日韩精品久久久久久密 | 亚洲欧美中文字幕日韩二区| 电影成人av| 精品国产国语对白av| 亚洲av电影在线进入| 日韩免费高清中文字幕av| www.av在线官网国产| 91精品国产国语对白视频| 久久av网站| 久久久久久久精品精品| 国产免费又黄又爽又色| 国产精品人妻久久久影院| 如何舔出高潮| 精品一区二区三区四区五区乱码 | 美女扒开内裤让男人捅视频| 另类亚洲欧美激情| 99香蕉大伊视频| 亚洲专区中文字幕在线 | 国产福利在线免费观看视频| 免费在线观看视频国产中文字幕亚洲 | 纵有疾风起免费观看全集完整版| 97人妻天天添夜夜摸| 久久国产精品大桥未久av| av卡一久久| 亚洲视频免费观看视频| 国产日韩一区二区三区精品不卡| 久久久欧美国产精品| 各种免费的搞黄视频| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 日韩欧美一区视频在线观看| 各种免费的搞黄视频| 国产精品麻豆人妻色哟哟久久| 亚洲伊人久久精品综合| 国产淫语在线视频| 中国国产av一级| 国产av一区二区精品久久| 久久影院123| 国产一卡二卡三卡精品 | 在现免费观看毛片| 亚洲av国产av综合av卡| 精品久久久精品久久久| 亚洲精品aⅴ在线观看| a级毛片黄视频| 丝袜人妻中文字幕| 热99国产精品久久久久久7| av.在线天堂| 最近的中文字幕免费完整| 十八禁人妻一区二区| 在线观看www视频免费| 久久久亚洲精品成人影院| 黄网站色视频无遮挡免费观看| 69精品国产乱码久久久| 日本欧美国产在线视频| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久久久免| 十八禁人妻一区二区| 男人舔女人的私密视频|