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

    Varying RossbyWave TrainsfromtheDevelopingtoDecayingPeriodof theUpper Atmospheric Heat Source over the Tibetan Plateau in Boreal Summer

    2018-06-29 08:24:30ChuandongZHURongcaiRENandGuoxiongWU
    Advances in Atmospheric Sciences 2018年9期

    Chuandong ZHU,Rongcai REN,and Guoxiong WU

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

    2Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters and KLME,Nanjing University of Information Science and Technology,Nanjing 210044,China

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

    1.Introduction

    As a huge heat source of the atmosphere in summer,the thermal in fl uence of the Tibetan Plateau(TP)is known to be crucial for the formation of the atmospheric circulation patterns and climate anomalies not only within the Asian region but also across the globe(Ye and Wu,1998;Duan et al.,2012).The highly elevated surface sensible heat over the TP can drive strong vertical motion in summer,which can act as a powerful“air pump”to pump the lower air upwards and result in divergence and anticyclonic circulation in the upper troposphere(Ye and Yang,1979;Wu et al.,1997).Other observational evidence(Flohn,1960;Reiter and Gao,1982;Yanai et al.,1992)and an early numerical experiment(Krishnamurti et al.,1973)have also con fi rmed that it is the diabatic heating over the summer TP that generates the South Asian high(SAH)—an important circulation system in the upper troposphere in summer(Mason and Anderson,1963).In particular,it is the inhomogeneous spatial distribution of the diabatic heating over the TP and eastern Asia that determines the locations of the SAH and the northwestern Paci fi c subtropical high(NWPSH)(Wu and Liu,2000;Liu et al.,2001a;Liu et al.,2004;Wu et al.,2015).

    Althoughthe land–sea thermalcontrast is knowntobethe main driver of monsoonal circulation(Murakami and Ding,1982;Luo and Yanai,1983,1984;Chen et al.,1991),the strong heat source over the TP has been recognized as one of the most important factors modulating the onset and intensity of the Asian summer monsoon(ASM)(e.g.Yanai and Wu,2006;Wu et al.,2007;Romatschke and Houze,2011).For example,the sensible heating over the TP originated in spring canleadtoanearlierreversalofthemeridionalgradientofthe atmospheric temperature over the east side of the TP,which leads to the fi rst onset of the ASM in the northeastern part of the Bay of Bengal(Wu and Zhang,1998).When the sensible heating over the TP is anomalously strong in spring,the seasonal reversal of the vertical shear of both the zonal and meridional wind in the East Asian monsoon region is anomalously early,which marks an earlier onset of the East Asian summer monsoon(EASM)(Duan et al.,2013).In particular,recent evidence based on data analysis and numerical experiments provided by Wu et al.(2012)further proved that the variability of the EASM,as well as the eastern part of the South Asian summer monsoon,can indeed be attributed to the thermal forcing anomalies over the TP.

    Owing to the high elevation and vast area of the TP heating,its in fl uence can be passed to remote regions outside Asia.Coupled with the SAH and strong upward motion over the TP,compensative downward motion always prevails over the vast surrounding regions(Ye and Wu,1998).The downward motion to the east of the TP can be as far away as the central North Paci fi c,and that to the west—over Afghanistan,Iran,and Saudi Arabia—is largely responsible for the dry climate and deserts in these regions(Yang et al.,1992;Rodwell and Hoskins,1996;Ye and Wu,1998;Wu et al.,2009).In the meridional direction,the divergent fl ow overtheTPmayextendsouthwardsacrosstheequator and induce downward motion in the southern subtropics,thus forming the well-known Hadley-monsoonal meridional circulation(Ye and Wu,1998).It may also extend northwards and then descend to the north of the TP,corresponding to the Taklimakan and Dzungaria deserts(Ye and Yang,1979).On the other hand,the in fl uence of the TP may also be passed to remote regions via a Rossby wave train(RWT).The modeling results of Hoskins and Karoly(1981)and Huang and Gambo(1981)suggested that the atmospheric responses to mechanical and thermal forcing of large-scale topographies may exhibit Rossby wave patterns.Two southeastward-propagating RWTscanbeidenti fi edinsummeracrosstheNorthernHemisphere,one of which is generated from the Rocky Mountains and the other from the TP(Grose and Hoskins,1979).Huang(1985)also indicated that the propagation of Rossby waves forced by a strong TP thermal forcing may result in anticyclonic circulation at 300 hPa over the TP,Sea of Okhotsk,and the northern part of North America,and cyclonic circulation over the west of Lake Baikal and the Aleutian region.In particular,modeling studies by Wu et al.(1997)and Ye and Wu(1998)emphasized that the elevated surface sensible heating over the TP can excite an RWT pattern across the globe.The RWT splits over the North Paci fi c,with its southern branch bending southward along the westerlies at the bottom of the tropical upper-tropospheric trough(TUTT);and its northern branch propagating further eastward to Alaska and Baffin Bay,and then bending southward over the North Atlantic to the north of Cape Verde,where the axis of the subtropical high lies.A similar RWT pattern emanating from the TP was also identi fi ed by Liu et al.(2002)and Wu et al.(2002).Theseresultsindicatedaclearpathwayfortheremote in fl uences of TP surface heating on the atmospheric circulation and climate from the subtropics to the high latitudes.

    However,the TP heat source in summer tends to be dominated by convective condensation heating in the upper troposphere,especially over the eastern TP region(Chen et al.,1985;Jiang et al.,2016),although the surface sensible heating is known to be the primary driver of the related vertical motion over the TP(Wu et al.,2002).This implies a further elevated atmospheric heat source over the TP,which may excite RWTs that are di ff erent from those excited by the surface sensible heat source.In addition,on top of the interannual variabilities of the TP heat source upon which most previous studies have focused(Krishnamurti et al.,1973;Reiter and Gao,1982;Wu et al.,1997;Ye and Wu,1998;Zhao and Chen,2001;Ueda et al.,2003;Duan et al.,2005,2013;Liuetal.,2007a;Zhangetal.,2009),signi fi cantintraseasonal timescale(10–30 days)variations have been identi fi ed for the column-integrated apparent heat source,especially over the eastern TP(Nitta,1983).Wang and Duan(2015)further indicated the dominance of a quasi-biweekly(10–20 days)oscillation mode of the precipitation,diabatic heating and tropospheric circulation over the TP.In terms of the origin of this quasi-biweekly oscillation over the TP,agreement is still lacking as to whether it originates locally or is in fl uenced by circulation in the tropics or in the middle and high latitudes(Wang and Duan,2015).Consistent with the quasi-biweekly oscillation of the heat source over the TP,the SAH,particularly its longitudinal position,seems also to exhibit quasibiweekly variation in boreal summer(Tao and Zhu,1964;Liu et al.,2007b;Ren et al.,2007).This longitudinal shifting of the SAH has been found to be closely related to an opposing longitudinal shift of the NWPSH in the lower and mid troposphere(Tao and Zhu,1964;Ren et al.,2007),which in turn is closely related with extreme changes in monsoonal rainfall over the East Asian region.The present study focuses on the daily variations of the upper-tropospheric heat source over the TP,and investigates the RWT related to this TP heat source in boreal summer.

    In particular,we examine the dynamic and thermodynamic processes related to the daily variations of the upper atmospheric heat source over the Tibetan Plateau(TPUHS)and the associated RWT,and demonstrate the possible climatic impacts of the RWT on the atmospheric circulation and weather anomalies in summer in the northern extratropics.The aim in carrying out this work is to help improve our understanding of the climatic e ff ects of TP heating,on multiple timescales,as well as the remote in fl uence of the TP on the northern extratropics in summer.

    Following this introduction,section 2 introduces the data and method;section 3 demonstrates the seasonal evolution of the TPUHS in summer and its dominant sub-monthly variability;section 4 illustrates the di ff erent RWTs related to the TPUHS from its developing to decaying period;section 5 explores the possible e ff ects of the TP-related RWTs on the climatological circulation systems and weather anomalies in the northern extratropics;and section 6 provides a summary and discussion.

    2.Data and method

    2.1.Data

    The daily atmospheric circulation fi elds are from the National Centers for Environmental Prediction and National Center for Atmospheric Research(NCEP–NCAR)reanalysis dataset(Kalnay et al.,1996),which covers the period from 1 January 1979 to 31 December 2014,and has a horizontal resolution of 2.5?latitude×2.5?longitude and 17 pressure levels spanning from 1000 to 10 hPa.The daily surface sensible heat fl ux fi elds from NCEP–NCAR covers the same period as the atmospheric circulation fi elds,and are interpolated from Gaussian grids to 2.5?latitude×2.5?longitude grids before analysis.The daily precipitation fi elds are from the Global Precipitation Climatology Project(GPCP),version 1.2,1?Daily Precipitation Dataset(Bolvin et al.,2009),covering the period from 1 January 1997 to 31 December 2014.We fi rst calculate the daily annual cycle of all variables by taking an average of them across the 36 years for each calendar day,and then we remove the daily annual cycle from the respective original fi elds to obtain the daily anomaly fi elds.The reliability of the NCEP–NCAR reanalysis in representing the diabatic heating over the TP and its variabilities has already been validated in many previous studies(e.g.,Nigam et al.,2000;Rodwell and Hoskins,2001;Duan,2003;Wang et al.,2011).We also con fi rmed the general consistency of the NCEP–NCAR reanalysis with other datasets(e.g.,JRA-55)in representing the TP heat source and the related RWTs.

    2.2.Method

    2.2.1. Heat source index

    Based on the large-scale heat budget equation(Yanai et al.,1992):

    where Tairdenotes air temperature;is the horizontal wind vector;ω is the vertical velocity in pressure(p)coordinates;θ is the potential temperature;p0equals 1000 hPa;and κ=R/cp,with R being the gas constant and cpthe speci fi c heat capacity of dry air at constant pressure,we fi rst obtain the apparent heat source(Q1) fi elds at 19 pressure levels(1000,925,850,700,600,550,500,450,400,300,250,200,150,100,70,50,30,20,and 10 hPa).Because the main focus of this study is the heat source in the upper troposphere,ratherthanthesurfacesensibleheatsource,wethenvertically integrate Eq.(1)from the tropopause pressure level(pt=100 hPa)to the lowest pressure level(pb)that is about 100 hPa away from the surface to obtain the column-integrated atmospheric heat source

    By averaging the(the angle brackets denote column integration)over the TP region(25?–40?N,70?–105?E)and above the altitude of 1500 m,we further de fi ne a daily TP apparent heat source index(THI).

    We also calculate the diabatic heating rate(K d?1)from dθ/dt(θ is the potential temperature)for each pressure level and the surface(Wu et al.,1999),to show the spatiotemporal evolution of the atmospheric heat source from the surface to the upper layer.

    2.2.2.Statistical signi fi cance test

    We apply the Student’s t-test to assess the statistical signi fi cance of our results,where the e ff ective sample size(N0)is estimated using the method discussed in Bretherton et al.(1999):

    where N is the original sample size,and r1and r2are the lag 1 auto-correlation coefficients of the two time series involved.As the lag 1 auto-correlation for the THI is 0.678,we have N0=0.37N,assuming that r1=r2,when we assess the statistical signi fi cance of the regressions against the THI.

    2.2.3. Wave activity fl ux

    To facilitate a subjective representation of Rossby waves,we adopt the three-dimensional wave activity fl ux in a spherical log-pressure coordinate,de fi ned by Plumb(1985):

    where

    represents the static stability(the caret denotes an area average over the region north of 20?N),Φ is the geopotential,u and v are the zonal and meridional components of horizontal wind,the primes denote the deviations from zonal means,Tairis the air temperature,H is the scale height,λ is longitude,? is latitude,a is Earth’s radius(6378.388 km),and ?is Earth’s rotational angular velocity(7.2685×10?5rad s?1).The primes denote the deviations from zonal means.Because the wave activity fl ux is parallel to the local group velocity of Rossby waves,it can e ff ectively indicate the propagating direction of the wave energy.

    3.Seasonal evolution of the TPUHS over the summer TP and its dominant quasibiweekly variability

    3.1.Seasonal evolution of the TPUHS

    Figure 1a shows the seasonal evolution of the area-mean column-integrated(from 100 hPa to the pressure level that is about 100 hPa away from the surface)apparent heat source()and the area-mean surface sensible heat fl ux(SH)over the TP region.It can be seen that the intensity of the atmospherichas already become lager than the surface SH since May.Thereafter,while the atmospherickeeps increasing after May and reaches 90–95 W m?2in June–July,the surface SH shows a weaker increase and remains at around 50 W m?2from May to October.This clearly indicates the dominance of strong latent condensation heating(LH)on the variation of the TP heat source in summer.After summer,both the SH anddecrease sharply from September,with the SH gradually becoming stronger than theand beginning to dominate the TP heat source again.

    Fig.1.Seasonal evolution(May–October)of(a)the area-mean surface sensible heat fl ux(blue;units:W m?2)and vertically accumulated apparent heat source(red;units:W m?2)over the TP(25?–40?N,70?–105?E)above 1500 m,and(b)the 200-hPa zonal wind(units:m s?1)averaged from 70?E to 105?E.The horizontal black lines in(b)roughly mark the south and north boundaries of the TP,and the vertical blue lines indicate 18 July and 23 August,respectively.

    The meridional location as well as the zonal structure of the upper-tropospheric westerly jet can a ff ect the propagation speed,amplitude and phase of the Rossby waves excited by a tropospheric heat source(Hoskins and Karoly,1981).By examining the successive seasonal shift of the East Asian upper-tropospheric westerly jet,Lin and Lu(2008)noted that it usually undergoes a signi fi cant northward shift around 20 July,and then stays stably to the north of 40?N,or to the north of the TP,for some time.To con fi rm this,we display the seasonal evolution of the 200-hPa westerly jet(>20 m s?1)over the TP region(70?–105?E)in Fig.1b.It is apparent that the westerly jet shifts gradually northward from May to mid June and begins to be centered at 40?N;then,it exhibits a further northward shift around 18 July and remains steadily at~ 42.5?N until 23 August,after which it begins to retreat southward.To consider the possible RWTs associated with the intensive deep convection and strong TPUHS covering most of the TP region,we focus mainly on the period from 20 July to 20 August in the following investigation.Actually,for the TPUHS in summer and its synoptic variations,the related RWTs during this period are the mostcompleteand strongest,compared with other time periods(data not shown).

    3.2.The dominant quasi-biweekly variability of the TPUHS

    Fig.2.Lead/lag regressions against the THI(units:W m?2)of the area-mean diabatic heating rate anomalies(units:K d?1)over the TP(a)at each level(contour interval:0.1;dashes for 0.05 and 0.15 are added)and(b)at the surface in boreal summer(20 July to 20 August).Stippling in(a)and dots in(b)indicate regressions above the 90%con fi dence based on the t-test.The abscissa is the lead time of the THI(units:d);the ordinate is the pressure level(units:hPa)in(a)and the amplitudes of the regression in(b).

    To demonstrate the dominant variability of the TPUHS in boreal summer,we display in Fig.2a the lead/lag regressions of the layer-by-layer area-mean diabatic heating rate anomalies over the TP against the THI(contours).It is apparent that,relative to the maximum values around the peak day of the THI,signi fi cant regression values for most layers decrease sharply with the increase in the lead/lag time,and consistently tend to approach zero until day±7.By referencing to the lead/lag auto-correlation of the THI time series and the corresponding e-folding time shown in Fig.3,we can estimate that the lifetime of an anomalous THI event is about seven days,or the dominant variability of the THI(i.e.,the upper-atmospheric heat source in the troposphere over the TP)is in a quasi-biweekly timescale.Spectral analysis confi rms that the THI does exhibit signi fi cant quasi-biweekly(or 10–20 days)oscillation and the variance explained by the quasi-biweekly variability is up to 35.8%of the total variance(~54.2)of the normalized THI timeseries(data not shown).As already mentioned in the introduction section,thequasi-biweeklyoscillationisthedominantmodeoftheintraseasonal variations of the column-integrated apparent heat source(Wang and Duan,2015).Consistent with this,Yang et al.(2016)identi fi ed a quasi-biweekly mode of the intraseasonal oscillations in the eastern TP summer rainfall.

    In addition,Fig.2a also shows that the maximum regression center of the diabatic heating rate lies in the upper troposphere,and a maximum center appears in the further higher layer up to 150 hPa around the peak day of the THI[day(0)].This clearly indicates the further vertical deepening of the convective LH with the maturing of TPUHS events.From the lead/lag regressions of the surface diabatic heating rate over the TP against the THI,it is easy to identify that the surface SH is always much weaker than the LH in the upper atmosphere throughout the entire period(Fig.2b).This again con fi rms that the sub-monthly variations of the TP heat source are indeed dominated by condensation LH in the upper atmosphere.

    Fig.3.Lead/lag auto-correlation of the THI(asterisks,left ordinate),and crosscorrelation(right ordinate)between the THI and WTI1(dots)and that between THI and WTI2(circles).The thickened curve indicates the THI auto-correlation that is above 1/e.The horizontal solid,long-dashed and short-dashed lines denote the 99%con fi dence level for the THI auto-correlation,the THI–WTI1 cross-correlation,and the THI–WTI2cross-correlation,respectively,based on t-test.

    Fig.4.Regression of(a)column-integrated apparent heat source(units:W m?2)and(b)surface sensible heat fl ux(units:W m?2)against the THI.Stippling indicates regressions above the 95%con fi dence based on the t-test.Blue contours delineate the TP’s topographic boundary of 1500 m.

    Figure 4 further shows the spatial regression patterns of the column-integrated heat source(Fig.4a)and that of the surface sensible heat fl ux(Fig.4b)in boreal summer.It can be seen that anomalous TPUHS events are characterized by signi fi cant anomalous hQ1i over most of the TP region,except its far western part(Fig.4a).In contrast,signi fi cant SH anomalies are only found over the western TP and along the southern boundary area of the TP.This further con fi rms the overwhelming dominance of the upper-atmospheric LH source over the SH heat source in boreal summer.

    3.3. Accompanying variations of the local circulation anomalies over the TP

    Figure 5 shows the lead/lag regressions of the local mean circulation anomalies over the TP region against the THI in boreal summer.As an anomalous TPUHS grows[before day(0)],anomalous convergence/divergence in the lower/upper troposphere(shading)and anomalous ascent in the troposphere intensi fi es(contours),indicating the development of anomalous deep convection over the TP(Fig.5a).And after the peak[day(0)]of the anomalous TPUHS,they begin to weaken and gradually diminish within about four days[day(+4)].Because of the greatly deepened convection and the elevated strong LH release,the anomalous warm center that was originally located around 300 hPa before day(?4)is further elevated to about 150 hPa and above until the peak day[day(0)](Fig.5b).Although gradually weakening after day(+2),the warm center stays in this layer throughout the decaying period of the TPUHS.Owing to the much-elevated anomalous warm center,the anomalous geopotential high center that was originally centered at 300 hPa in the upper troposphere before day(?2)is temporarily replaced by negative height anomalies from the surface to the upper layer(up to 100 hPa)around day(0)(Fig.5c).However,these negative height anomalies,resulting from the thermal e ff ect of greatly elevated LH release in the upper troposphere,are relatively weaker above 200 hPa and only exist within two days around the peak day of the heat source.This is due to the counteracting dynamic e ff ect of the persisting divergence/convergence in the upper/lower troposphere(Fig.5a),which is always supporting an anomalous high/low in the upper/lower troposphere.Also associated with the greatly elevated heat source,the horizontal vorticity,originally characterized by positive/negative anomalies in the lower/upper troposphere,begins to be dominated by positive vorticity anomalies throughout the troposphere up to 50 hPa from day(?1)(Fig.5d).These results clearly indicate that the strong thermal e ff ect of the greatly elevated LH around the peak day of the TPUHS may temporarily o ff set the anomalous circulation pattern(anticyclonic circulation above cyclonic circulation)that prevails during the earlier developing period of the TPUHS.

    Fig.5.Lead/lag regressions against the THI of the area-mean(a)horizontal divergence(shaded;units:10?6 s?1)and vertical velocity(contours;units:10?3Pa s?1;only those statistically signi fi cant above the 90%confi dence level are shown),(b)air temperature(units:?C),(c)geopotential height(units:gpm),and(d)relative vorticity(units:10?6s?1)in the atmospheric column over the TP.Stippling in all panels indicates regressions above the 90%con fi dence level based on the t-test.The abscissa is the lead time of the THI(units:d);the ordinate is the pressure level(units:hPa).

    4.Di ff erentRWT patterns during the developing/decaying period of anomalous TPUHS

    4.1.RWT patterns

    4.1.1.From vorticity anomalies

    Figure 6 shows the lead/lag regressions of the 200-hPa zonal wind(contours)and relative vorticity(shading)anomalies in the Northern Hemisphere against the THI during summer.Firstly,it is clear that,when the TPUHS or the convective LH anomaly is weak and begins to grow from day(?7)to day(?5)(Figs.6a and b),signi fi cant negative vorticity anomalies dominate the TP,indicating an anomalous SAH relatedtotherelativelyweakandshallowTPUHSaccompanied by minor surface SH over the summer TP.Meanwhile,positive vorticity anomalies exist in two areas from the Caspian region to Kazakhstan and from northeastern China to Japan.Correspondingly,positivezonalwind anomaliesappeartothe north of the TP and extend along the latitudinal band from the northwestern TP to southern Japan.Together with the negative zonal wind anomaly centers existing north of the positive vorticity anomaly centers mentioned above,the zonal wind anomalies favor a slightly weakening and southward-shifting of the subtropical westerly jet(SWJ)axis(thick dashed line)located just to the north of the TP(Fig.6b).Along the SWJ axis,a clear teleconnection pattern of vorticity anomaly centers can be identi fi ed downstream as well as upstream of the TPduringthistime.Speci fi cally,thedownstreamteleconnection pattern consists of fi ve centers(marked “P”,“Q”,“R”,“S”and “T”)from the negative vorticity center(“P”)over the TP to the negative vorticity center(“T”)over the northeastern Paci fi c and west coast of North America;and the upstream teleconnection pattern includes four centers(marked“M”,“N”,“O”and “P”)from “M”near the Baltic Sea to “P”over the TP(Fig.6b).As the TPUHS continues to grow,the negative vorticity center over the TP gradually shifts south-ward from day(?4)to day(?3)(Figs.6c and d).Accompanying this is the southward-shifting westerly anomaly center to its north,thus favoring a further southward shift of the SWJ axis.The negative vorticity center“R”tends to split and weakens and shifts westward,while the“Q”center becomes strengthened and shifts westward to North China.A new branch of the downstream teleconnection pattern seems to appear,from “Q”to“S”via the negative vorticity center“R1”over the central North Paci fi c(Fig.6d).Also associated with the weakening and southward shifting of the SWJ axis,the upstream teleconnection pattern is signi fi cantly weakened until day(?3)(Fig.6d).

    Fig.6.Lead/lag regressions against the THI of the relative vorticity anomalies(shaded;units:10?6s?1)and the zonal wind anomalies(contours;units:m s?1;interval:1;zero line omitted for clarity)at 200 hPa when the THI leads for(a)?7 days,(b)?5 days,(c)?4 days,(d)?3 days,(e)?1 day,(f)0 days,(g)2 days,(h)4 days,(i)6 days and(j)8 days.Hatched and stippled areas denote the 90%con fi dence level for relative vorticity anomalies and for zonal wind anomalies,respectively.The thick dashed line denotes the climatological westerly jet axis(>20 m s?1)during the period from 20 July to 20 August.The numbers in the upper-left corner of each panel indicate the lead time of the THI.The blue letters mark the RWT and the anomaly centers along it.

    After,when the TPUHS is further intensi fi ed,the negative vorticity center“P”shifts further southward,away from the TP region,until day(?1)(Fig.6e).A much stronger westerly anomaly center is then located just over the TP region,which corresponds to a maximum positive vorticity anomaly center to its north.Thereafter,a meridional tripole pattern of vorticity anomalies,as well as zonal wind anomalies,from the south to the north of the TP region,forms(Figs.6e and f).The tripole pattern exists from day(?1)to day(+1)(data not shown),while the original teleconnection pattern downstream of the TP no longer exists during this time(Figs.6e and f).This implies that the maturing of an anomalous TPUHS event can temporarily disrupt the teleconnection pattern that has developed during the developing period of the TPUHS,by o ff setting the anomalous atmospheric circulation structure over the TP,as already shown in Fig.5.

    With the decay of anomalous TPUHS events,and the weakening of the negative vorticity center south of the TP,both the strong westerly anomaly center over the TP and the associated meridional tripole pattern mentioned above quickly diminish from day(+2)to day(+4)(Figs.6g and h).However,the negative vorticity center north of the TP that originally extended eastward,begins to develop near the Baikal region(denoted as “B”).Later,we show that this anomalous anticyclone center is dynamically forced and exhibits a deep barotropic structure.With the intensi fi cation of this“B”center and the associated strong easterly anomaly center to its south afterwards,the SWJ tends to be mainly weakened in the upstream region from the eastern Caspian Sea to the northern TP,whereas it mainly shifts southward over its exit region(Figs.6g and h).Related to this,a new teleconnection pattern(B–C–D–E)seems to gradually form downstream along the SWJ axis.This new teleconnection pattern,with four vorticity anomaly centers(marked “B”,“C”,“D”and “E”)from region “B”near Baikal to “E”over the northeastern Paci fi c,becomes the most apparent until day(+8),when the TPUHS event has largely decayed(Fig.6j).

    4.1.2.From geopotential height and wave fl ux anomalies

    Fig.7.Lead/lag regressions against the THI of the geopotential height anomalies(contours;units:gpm;interval:4),the wave activity fl ux anomalies(vectors;units:10?2m2s?2;amplitudes larger than 0.6 are shown)calculated from the lead/lag regressions of stream function,and the wave activity fl ux divergence anomalies(shaded;units:10?5m s?2)at 200 hPa when the THI leads for(a)?7 days,(b)?5 days,(c)?4 days,(d)?3 days,(e)?1 day,(f)0 days,(g)2 days,(h)4 days,(i)6 days and(j)8 days.Stippling denotes the 90%con fi dence level for geopotential height anomalies,and the blue lines in(c),(d)and(j)delineate the Rossby wave rays.The numbers in the upper-left corner of each panel indicate the lead time of the THI.The blue letters and rays mark the RWT and the anomaly centers along it.

    To con fi rm that the TPUHS-related teleconnection patterns identi fi ed above are RWTs,we similarly display in Fig.7 the lead/lag regressions of geopotential height and horizontal Eliassen–Palm(EP) fl ux and its divergence anomalies.It is clear that the height anomaly centers exhibit teleconnection patterns similar to those of the vorticity anomaly centers in Fig.6;plus,the signi fi cant anomalous EP fl ux fl ow that links the teleconnection centers clearly indicates that the teleconnection patterns are indeed RWTs.Speci fi cally,for the teleconnection pattern downstream of the TP at day(?7)when the TPUHS is weak,the wave- fl ux fl ow is mainly in the midlatitudesalongthe zonaldirection.Subsequently,the teleconnection anomaly centers during this time are seen mainly to be linked by this midlatitude wave- fl ux fl ow from the Baikal region,while the wave- fl ux fl ow from the TP region is still relatively trivial(Fig.7a).With the growth of the TPUHS and the southward shifting of the SWJ later on,the midlatitude wave- fl ux fl ow gradually weakens,while that originating from the TP region increases and gradually becomes the main wave- fl ux fl ow linking the anomaly centers(P–Q–R1–S–T)(Figs.7b–d).This suggests that,during the developing period of the TPUHS,the TP-related wave- fl ux fl ow can strongly modulate the RWT pattern in the extratropics,and is mainly responsible for a signi fi cant RWT pattern that extends from the TP to the northeastern Paci fi c and western coastal region of North America.

    However,it should be noted that the wave- fl ux direction accompanying the upstream teleconnection pattern is southeastward rather than northwestward,indicating that this RWT is not originated from the TP region but from the North Atlantic region.Based on the climatological mean wave fl ux pattern(data not shown),we know that this RWT is a recurrent teleconnection pattern that plays an important role in linking the climate anomalies over the North Atlantic region to those over the Asian region(Chang et al.,2001;Buermann et al.,2005;Ren and Hu,2014).Nevertheless,this recurrent RWT pattern largely diminishes and breaks up around the TPUHS peak from day(?1)(Figs.7e and f),as also shown in Fig.6.

    Around the TPUHS peak from day(?1),the anomalous high over the TP during the developing period is gradually replaced by an anomalous low,as already shown in Fig.5.This displacement of the geopotential anomaly center over the TP also results in the break-up of the downstream RWT(P–Q–R1–S–T)due to the displacement of the wave fl ux source(Figs.7e and f).Furthermore,the anomalous low over the TP(especiallyitswesternpart)inturnseemstoactasanewwave source,which points to the Baikal region and corresponds to signi fi cant intensi fi cation of an anomalous high center over the Baikal region from day(+1)(data not shown)to day(+4)(Figs.7g and h).Although the anomalous TP low weakens quickly after day(+1),the wave- fl ux fl ow from the western TP to the Baikal region is persistently clear(Figs.7g–i).This implies a clear linkage between the anomalous TP low and the anomalous Baikal high(center “B”),despite the zonal wave- fl ux fl ow in the midlatitudes seeming also to contribute to the intensi fi cation of center“B”during this time.After,center“B”becomes the main wave- fl ux divergence center for a new RWT that extends eastward,and a continuous eastward wave- fl ux fl ow from “B”to“E”becomes clear until day(+8)(Figs.7j).

    4.2.Barotropic structure of the RWTs

    After showing the di ff erent RWT patterns from the developing to the decaying period of the TPUHS,next we diagnose the vertical structure of the anomaly centers along the RWTs.To do this,we choose day(?3)and day(+8),when the RWTs are the most pronounced,and show in Fig.8 the vertical cross sections of the regressed geopotential height anomalies along the RWTs.It can be seen that the geopotential height anomaly centers along both RWTs consistently exhibit an equivalent barotropic structure.The exception is the anomaly center“P”over the TP at day(?3),which is anomalously positive/negative above/below 500 hPa and shows a general baroclinic structure.This demonstrates that these anomaly centers are indeed forced by dynamical processes associated with the RWTs;meanwhile,center“P”over the TP,which acts as the main wave source of the RWT P–Q–R1–S–T from day(?5)to day(?3)during the developing period of the TPUHS,is primarily thermally forced.

    Fig.8.Vertical cross sections of the regressed geopotential height anomalies(contours;units:gpm;interval:2)(a)along the P–Q–R1–S–T RWT at day(?3),and(b)along the B–C–D–E RWT at day(+8).Shaded areas indicate the 90%con fi dence level of regression,based on the t-test.

    Furthermore,although the anomalously positive“P”center becomes relatively weaker with the growth of the TPUHS until day(?3),it is obvious that other anomaly centers along the RWT P–Q–R1–S–T are stronger and exhibit a much deeper vertical structure than those along the RWT B–C–D–E(Fig.8a vs Fig.8b).The signi fi cant anomalies along the former RWT can span vertically from the lower troposphere tothestratosphereupto50hPa(Fig.8a),whereasthosealong the RWT B–C–D–E are mostly con fi ned to below 100 hPa—except center“B”,which also spans beyond 30 hPa.The vertically deep structure of the anomaly centers along the RWTs implies that the RWTs may have signi fi cant e ff ects on the circulation systems in the northern extratropics.

    5.E ff ects of the RWTs on circulation and weather in the northern extratropics

    5.1.RWT indices

    To demonstrate the e ff ects of the TPUHS-related RWTs on circulation and weather in the extratropics,we de fi ne two RWT indices to capture the synoptic variabilities of the RWTs more e ff ectively than using the THI index.The indices are based mainly on the RWT patterns on day(?3)and day(+8),when the RWTs are the most pronounced and complete.For the RWT during the developing period of the TPUHS,we choose to de fi ne the index based on the P–Q–R1–S–T branch.Speci fi cally,with this RWT as an example,we fi rst calculate the weight coefficients(wi)for the fi ve anomaly centers based on their relative intensities on day(?3),

    where hiis the maximum value of the geopotential height anomaly center i(P,Q,R1,S,T).We then use the wito defi ne the daily RWT index for the developing period of the TPUHS(denoted as WTI1)through

    where Hiisthe dailynormalized geopotential height anomaly time series for the anomaly centers.Similarly,the daily RWT index for the decaying period of the TPUHS(denoted as WTI2)is de fi ned as

    Where wjis as same as wi,but for the geopotential height anomaly center j(B,C,D,E).Obviously,the values of the WTI1and WTI2daily indices fall within the range 0–1.As expected,Fig.3 shows that the two de fi ned wave indices are signi fi cantly correlated with the THI,and the maximum correlationappearsonthedaywhenWTI1leadstheTHIbythree days[day(?3)],and on the day when WTI2lags the THI by seven days[day(+7)].Below,based on these two wave indices,we diagnose the e ff ects of the TPUHS-related RWTs on circulation and weather in the extratropics.

    5.2.E ff ects of the TPUHS-related RWT during the TPUHS developing period

    Figure 9 shows the regressed anomaly patterns of the geopotential height at 200 hPa and 500 hPa(shading in Figs.9a and b),surface air temperature(Ts,Fig.9d),and 850 hPa wind and daily precipitation(Fig.9c)anomalies,against the fi rst RWT index(WTI1)de fi ned above.Firstly,it is apparent that the geopotential height anomaly patterns in Figs.9a and b are largely similar to that shown in Fig.7d,because of the signi fi cant relationship between the THI and WTI1shown in Fig.3.From the circulation patterns of the summer climatology(contours in Figs.9a and b),it is clear that the growth of the TPUHS can indeed exert a series of local and remote impacts on the circulation systems along the RWT in the troposphere.Speci fi cally,the positive “P”over the southern TP and the stronger negative“Q”together yield a slightly stronger but southward-shifted SAH along with a southward shifting of the SWJ(Fig.9a;also shown in Fig.7d).This corresponds to an anomalously stronger and southward-shifted NWPSH at 500 hPa over southeastern China(Fig.9b).On the other hand,the coexisting negative“Q”and downstream positive “R1”centers support an eastward extension of the SAH to the northwestern Paci fi c region(Fig.9a),and a strengthening of the NWPSH below in the middle troposphere over that region,which in turn favors a patternwithalowtothewestandahightotheeast(ora“l(fā)owwest–high-east”pattern)over the region from eastern China to the northwestern Paci fi c region(Fig.9b).Therefore,we can see that the wave-train-related anomalous precipitation zone in East Asia is southwest–northeast oriented—from the southeastern TP,the middle and lower reaches of the Yangtze and Huaihe river valleys,to Japan—where the anomalous lower-tropospheric fl ow converges(Fig.9c).Accompanying this is anomalously low Tsover the region from East China to Japan(Fig.9d).Meanwhile,negative precipitation and warm Tsanomalies exist to the south and north of the positive precipitation anomaly zone,respectively(Figs.9c and d).

    Similarly,by supporting the formation of a low-west–high-east pattern over the west coast of North America,the negative“S”anddownstreampositive“T”centersareresponsible for the positive precipitation anomaly center in that region,as well as the cold Tsanomalies over a wider region of the North Paci fi c,Alaska and northeast of Kamchatka(Figs.9a–d).Meanwhile,the positive geopotential height anomalies related to center“T”over the eastern North Paci fi c act to weaken the TUTT over the northern Paci fi c(Fig.9a),which is another important circulation system in the upper troposphere in northern summer(Kelly and Mock,1982).The weakening of the TUTT obviously weakens the southwest-erly fl ow from the eastern North Paci fi c to the western part of the United States,which should prevail in boreal summer(Fig.9a).This corresponds to the warm Tsanomalies over the North America region(Fig.9d),although no signi fi cant negative precipitation anomalies appear(Fig.9c).

    Fig.9.The regressed(a)200-hPa and(b)500-hPa geopotential height anomalies(shaded;units:gpm),(c)daily rainfall and 850-hPa wind(vectors;units:m s?1)anomalies(shaded;units:mm),and(d)surface air temperature anomalies(shaded;units:?C)against the RWT index WTI1.Contours denote the corresponding summertime(20 July to 20 August)climatology(units:gpm;interval:40).Stippling indicates the 90%con fi dence level of the regression,based on the t-test.

    5.3.E ff ects of the TPUHS-related RWT during the TPUHS decaying period

    Firstly,it can be seen from Fig.10 that the precipitation and Tsanomalies associated with the second TPUHS-related RWT are much weaker and largely scattered(Figs.10c and d),due to the relatively weaker anomaly centers along the RWT(Figs.10a and b).Since the SAH as well as the NWPSH show less signi fi cant changes during this time(Figs.10a and b),the precipitation anomalies are largely scattered over most of East and South China,and the Tsanomalies are only weakly negative over the region from East China to Japan(Figs.10c and d).Nevertheless,it is clear that the positive“B”and negative “C”centers support a high-west–low-east pattern from the Baikal region to the east of Japan,responsible for the strong and warm Tsanomalies over a wider regionofNortheastAsia(Fig.10d).Aroundtheregionbetween the negative“C”and downstream positive“D”center,significantly positive precipitation anomalies exist over the North Paci fi c to the east of the Japan region(Fig.10c).In addition,the geopotential height anomalies associated with the positive“D”,and downstream negative“E”centers,and the even weaker positive center southeast of center“E”,contribute to deepen the TUTT(Fig.10a);however,this only results in weak and cold Tsanomalies over the northwest region of North America and over the northeastern North Paci fi c(Fig.10d),due to the weak intensity of the centers along this RWT.

    Fig.10.As in Fig.9 but for the regression against the RWT index WTI2.

    6.Summary

    This study concerns the sub-monthly variability of the upper atmospheric heat source over the TP and the related RWTs in the northern extratropics in summer.To e ff ectively capture the variability of the TPUHS over the entire TP region,our study mainly focuses on the period from 20 July to 20 August,when the climatological westerly jet resides just to the north of the TP.We fi rst identify that the dominant variability of the TPUHS in summer(July–August)is quasibiweekly;and with the intensive development of an anomalous TPUHS event,the further elevated LH heating may temporarily o ff set the anomalous atmospheric circulation structure(anticyclonic circulation above cyclonic circulation)that prevails during the earlier developing period of the TPUHS,temporarily resulting in a low anomaly center over the TP in the upper troposphere where the SAH usually exists.However,this thermally forced anomalous low center diminishes rapidly due to the counteracting dynamic e ff ect of the persisting lower/upper convergence/divergence over the TP.

    During the developing period of a TPUHS event,an anomalous high corresponding to weak and shallow LH over the TP acts as the main wave source of a southernbranch RWT(P–Q–R1–S–T)along the SWJ from day(?5)to day(?3),which extends northeastward and eastward to the northeastern Paci fi c and west coast of North America.However,accompanying the southward shift of the anomalous TP high,and the replacement by an anomalous TP low around the TPUHS peak,this TPUHS-related RWT is broken up from day(?1).Nevertheless,the existing anomalous low during this time then temporarily acts as a wave source to connect to an anomalous high over the Baikal region.After the anomalous TP low quickly diminishes,this strong anomalous Baikal high becomes the main wave source of a new RWT during the decaying period of the TPUHS.This RWT extends eastward from the Baikal region to the northeastern Paci fi c around eight days after the TPUHS peak.

    Compared to that during the decaying period of the TPUHS,the RWT pattern during the developing period is stronger and presents a much deeper barotropic structure.Correspondingly,the signi fi cant extratropical in fl uences of the earlier RWT during the developing period of the TPUHS are much stronger and more systematic than those of the later RWT during the decaying period of the TPUHS.By modulating the circulation systems in the northern extratropics,the earlier RWT is associated with an anomalous southwest–northeast-oriented precipitation band extending from the southern TP and the middle and lower reaches of the Yangtze and Huaihe river valleys to Japan,whereas the later RWT only yields a signi fi cant precipitation anomaly center over the North Paci fi c—mainly to the east of Japan.The Tsanomalies associated with the earlier RWT are also generally stronger than those associated with the later one.

    This study demonstrates the varying RWT patterns from the developing to the decaying period of anomalous TPUHS events,and investigates the related thermodynamic processes and the remote climate e ff ects of the TPUHS-related RWTs.The dynamic processes responsible for the linkage between the TPUHS variation and the changes in the RWT pattern are also illustrated.However,as already indicated in section 4,the RWT patterns identi fi ed in this study might be implicitly contributed by the midlatitude RWT prevailing in boreal summer,including the upstream RWT originated from the Atlantic.However,the possible e ff ects of midlatitude circulation on the local circulation anomalies over the TP are subject to further investigations in the future.Nevertheless,the demonstrated RWT patterns are what we can always observe during the life cycle of anomalous TPUHS events.In our group’s future work,we will also use general circulation models to simulate the RWT patterns that are purely excited by the anomalous TPUHS,and provide further evidence of the mechanisms,processes and the variations of the TPUHS-related RWTs.

    Acknowledgements.Thanks to the two anonymous reviewers for their constructive comments and suggestions.This work was jointly supported by the National Science Foundation of China(Grant Nos.91437105,41575041 and 41430533)and the China Meteorological Administration Special Public Welfare Research Fund(Grant No.GYHY201406001).We are grateful for the availability of the NCEP–NCAR reanalysis dataset via the website http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html,and the GPCP dataset via http://precip.gsfc.nasa.nov/.

    REFERENCES

    Bolvin,D.T.,R.F.Adler,G.J.Hu ff man,E.J.Nelkin,and J.P.Poutiainen,2009:Comparison of GPCP monthly and daily precipitation estimates with high-latitude gauge observations.Journal of Applied Meteorology and Climatology,48,1843–1857,https://doi.org/10.1175/2009JAMC2147.1.

    Bretherton,C.S.,M.Widmann,V.P.Dymnikov,J.M.Wallace,and I.Blad′e,1999:The e ff ective number of spatial degrees of freedom of a time-varying fi eld.J.Climate,12,1990–2009,https://doi.org/10.1175/1520-0442(1999)012<1990:TENOSD>2.0.CO;2.

    Buermann,W.,B.Lintner,and C.Bon fi ls,2005:A wintertime Arctic Oscillation signature on early-season Indian Ocean monsoonintensity.J.Climate,18,2247–2269,https://doi.org/10.1175/JCLI3377.1.

    Chang,C.P.,P.Harr,and J.H.Ju,2001:Possible roles of Atlantic circulations on the weakening Indian monsoon rainfall-ENSO relationship.J.Climate,14,2376–2380,https://doi.org/10.1175/1520-0442(2001)014<2376:PROACO>2.0.CO;2.

    Chen,L.X.,E.R.Reiter,and Z.Q.Feng,1985:The atmospheric heat source over the Tibetan Plateau:May-August 1979.Mon.Wea.Rev.,113,1771–1790,https://doi.org/10.1175/1520-0493(1985)113<1771:TAHSOT>2.0.CO;2.

    Chen,L.X.,Q.G.Zhu,H.B.Luo,J.H.He,M.Dong,and Z.Q.Feng,1991:The East Asian Monsoon.Chinese Meteorological Press,362 pp.(in Chinese)

    Duan,A.M.,2003:The in fl uence of thermal and mechanical forcing of Tibetan Plateau upon the climate patterns in East Asia.PhD Dissertation,Institute of Atmospheric Physics,Chinese Academy of Sciences,23–31.(in Chinese)

    Duan,A.M.,Y.M.Liu,and G.X.Wu,2005:Heating status of the Tibetan Plateau from April to June and rainfall and atmospheric circulation anomaly over East Asia in midsummer.Science in China Series D:Earth Sciences,48,250–257,https://doi.org/10.1360/02yd0510.

    Duan,A.M.,W.R.Wang,Y.H.Lei,and Y.F.Cui,2013:Trends in summer rainfall over China associated with the Tibetan Plateau sensible heat source during 1980–2008.J.Climate,26,261–275,https://doi.org/10.1175/JCLI-D-11-00669.1.

    Duan,A.M.,G.X.Wu,Y.M.Liu,Y.M.Ma,and P.Zhao,2012:Weather and Climate E ff ects of the Tibetan Plateau.Adv.Atmos.Sci.,29,978–992,https://doi.org/10.1007/s00376-012-1220-y.

    Flohn,H.,1960:Recent investigation on the mechanism of the“summer monsoon”of southern and eastern Asia.Proc.Symp.Monsoon of the World,New Delhi,Hind Union Press,75–88.

    Grose,W.L.,and B.J.Hoskins,1979:On the in fl uence of orography on large-scale atmospheric fl ow.J.Atmos.Sci.,36,223–234,https://doi.org/10.1175/1520-0469(1979)036<0223:OTIOOO>2.0.CO;2.

    Hoskins,B.J.,and D.J.Karoly,1981:The steady linear response of a spherical atmosphere to thermal and orographic forcing.J.Atmos.Sci.,38,1179–1196,https://doi.org/10.1175/1520-0469(1981)038<1179:TSLROA>2.0.CO;2.

    Huang,R.H.,1985:The in fl uence of the heat source anomaly over Tibetan Plateau on the northern hemispheric circulation anomalies.Acta Meteorologica Sinica,43,208–220,https://doi.org/10.11676/qxxb1985.026.(in Chinese)

    Huang,R.H.,and K.Gambo,1981:The response of a model atmosphere in middle latitude to forcing by topography and stationary heat sources.J.Meteor.Soc.Japan,59,220–237,https://doi.org/10.2151/jmsj1965.59.2_220.

    Jiang,X.W.,Y.Q.Li,S.Yang,K.Yang,and J.W.Chen,2016:Interannual variation of summer atmospheric heat source over the Tibetan Plateau and the role of convection around the western Maritime Continent.J.Climate,29,121–138,https://doi.org/10.1175/JCLI-D-15-0181.1.

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

    Kelly,W.E.,Jr.,and D.R.Mock,1982:A diagnostic study of upper tropospheric cold lows over the western north paci fi c.Mon.Wea.Rev.,110,471–480,https://doi.org/10.1175/1520-0493(1982)110<0471:ADSOUT>2.0.CO;2.

    Krishnamurti,T.N.,S.M.Daggupaty,J.Fein,M.Kanamitsu,and J.D.Lee,1973:Tibetan high and upper tropospheric tropical circulations during northern summer.Bull.Amer.Meteor.Soc.,54,1234–1249.

    Lin,Z.D.,and R.Y.Lu,2008:Abrupt northward jump of the East Asian upper-tropospheric jet stream in mid-summer.J.Meteor.Soc.Japan,86,857–866,https://doi.org/10.2151/jmsj.86.857.

    Liu,X.,W.P.Li,and G.X.Wu,2002:Interannual variation of the diabatic heating over the Tibetan Plateau and the Northern Hemispheric circulation in summer.Acta Meteorologica Sinica,60,267–277,https://doi.org/10.3321/j.issn:0577-6619.2002.03.002.(in Chinese)

    Liu,X.,G.X.Wu,W.P.Li,and Y.M.Liu,2001a:Thermal adaptation of the large-scale circulation to the summer heating over the Tibetan Plateau.Progress in Natural Science,11,207–214.

    Liu,Y.M.,G.X.Wu,and R.C.Ren,2004:Relationship between the subtropical anticyclone and diabatic heating.J.Climate,17,682–698,https://doi.org/10.1175/1520-0442(2004)017<0682:RBTSAA>2.0.CO;2.

    Liu,Y.M.,B.J.Hoskins,and M.Blackburn,2007b:Impact of Tibetan orography and heating on the summer fl ow over Asia.J.Meteor.Soc.Japan,85B,1–19,https://doi.org/10.2151/jmsj.85B.1.

    Liu,Y.M.,Q.Bao,A.M.Duan,Z.A.Qian,and G.X.Wu,2007a:Recent progress in the impact of the Tibetan Plateau on climate in China.Adv.Atmos.Sci.,24,1060–1076,https://doi.org/10.1007/s00376-007-1060-3.

    Luo,H.B.,and M.Yanai,1983:The large-scale circulation and heat sources over the Tibetan Plateau and surrounding areas during the early summer of 1979.Part I:Precipitation and kinematic analyses.Mon.Wea.Rev.,111,922–944,https://doi.org/10.1175/1520-0493(1983)111<0922:TLSCAH>2.0.CO;2.

    Luo,H.B.,and M.Yanai,1984:The large-scale circulation and heat sources over the Tibetan Plateau and surrounding areas during the early summer of 1979.Part II:Heat and moisture budgets.Mon.Wea.Rev.,112,966–989,https://doi.org/10.1175/1520-0493(1984)112<0966:TLSCAH>2.0.CO;2.

    Mason,R.B.,and C.E.Anderson,1963:The development and decay of the 100-MB.Summertime anticyclone over southern Asia.Mon.Wea.Rev.,91,3–12,https://doi.org/10.1175/1520-0493(1963)091<0003:TDADOT>2.3.CO;2.

    Murakami,M.,and Y.H.Ding,1982:Wind and temperature changes over Eurasia during the early summer of 1979.J.Meteor.Soc.Japan,60,183–196,https://doi.org/10.2151/jmsj1965.60.1_183.

    Nigam,S.,C.Chung,and E.DeWeaver,2000:ENSO diabatic heating in ECMWF and NCEP-NCAR reanalyses,and NCAR CCM3 simulation.J.Climate,13,3152–3171,https://doi.org/10.1175/1520-0442(2000)013<3152:EDHIEA>2.0.CO;2.

    Nitta,T.,1983:Observational study of heat sources over the eastern Tibetan Plateau during the summer monsoon.J.Meteor.Soc.Japan,61,590–605,https://doi.org/10.2151/jmsj1965.61.4_590.

    Plumb,R.A.,1985:On the three-dimensional propagation of stationary waves.J.Atmos.Sci.,42,217–229,https://doi.org/10.1175/1520-0469(1985)042<0217:OTTDPO>2.0.CO;2.

    Reiter,E.R.,and D.Y.Gao,1982:Heating of the Tibet Plateau and movements of the South Asian high during spring.Mon.Wea.Rev.,110,1694–1711,https://doi.org/10.1175/1520-0493(1982)110<1694:HOTTPA>2.0.CO;2.

    Ren,R.C.,and J.G.Hu,2014:An emerging precursor signal in thestratosphereinrecentdecadesfortheIndiansummermonsoon onset.Geophys.Res.Lett.,41,7391–7396,https://doi.org/10.1002/2014GL061633.

    Ren,R.C.,Y.M.Liu,and G.X.Wu,2007:Impact of South Asia High on the short-term variation of the subtropical anticyclone over western Paci fi c in July 1998.Acta Meteorologica Sinica,65,183–197,https://doi.org/10.11676/qxxb2007.018.(in Chinese)

    Rodwell,M.J.,and B.J.Hoskins,1996:Monsoons and the dynamics of deserts.Quart.J.Roy.Meteor.Soc.,122,1385–1404,https://doi.org/10.1002/qj.49712253408.

    Rodwell,M.J.,and B.J.Hoskins,2001:Subtropical anticyclones and summer monsoons.J.Climate,14,3192–3211,https://doi.org/10.1175/1520-0442(2001)014<3192:SAASM>2.0.CO;2.

    Romatschke,U.,R.A.J.Houze Jr.,2011:Characteristics of precipitating convective systems in the South Asian monsoon.Journal of Hydrometeorology,12,3–26,https://doi.org/10.1175/2010JHM1289.1.

    Tao,S.Y.,and F.K.Zhu,1964:The 100-mb fl ow patterns in southern Asia in summer and its relation to the advance and retreat of the West-Paci fi c subtropical anticyclone over the far east.Acta Meteorologica Sinica,34,387–396,https://doi.org/10.11676/qxxb1964.039.(in Chinese)

    Ueda,H.,H.Kamahori,and N.Yamazaki,2003:Seasonal contrasting features of heat and moisture budgets between the eastern and western Tibetan Plateau during the GAME IOP.J.Climate,16,2309–2324,https://doi.org/10.1175/2757.1.

    Wang,M.R.,and A.M.Duan,2015:Quasi-biweekly oscillation over the tibetan plateau and its link with the Asian summer monsoon.J.Climate,28,4921–4940,https://doi.org/10.1175/JCLI-D-14-00658.1.

    Wang,T.M.,G.X.Wu,and M.Ying,2011:Comparison of diabatic heating data from NCEP/NCAR(I,II)and ERA40.Acta Scientiarum Naturalium Universitatis Sunyatseni,50,128–134.(in Chinese)

    Wu,G.X.,and Y.M.Liu,2000:Thermal adaptation,overshooting,dispersion,and subtropical anticyclone Part I:Thermal adaptation and overshooting.Chinese Journal of Atmospheric Sciences,24,433–446,https://doi.org/10.3878/j.issn.1006-9895.2000.04.01.(in Chinese)

    Wu,G.X.,Y.M.Liu,and P.Liu,1999:The e ff ect of spatially nonuniform heating on the formation and variation of subtropical high I:Scale analysis.Acta Meteorologica Sinica,57,257–263,https://doi.org/10.11676/qxxb1999.025.(in Chinese)

    Wu,G.X.,and Y.S.Zhang,1998:Tibetan Plateau Forcing and the Timing of the Monsoon Onset over South Asia and the South China Sea.Monthly Weather Review,126,913—927,https://doi.org/10.1175/1520-0493(1998)126<0913:TPFATT>2.0.CO;2.

    Wu,G.X.,W.P.Li,H.Guo,H.Li,J.Xue,and Z.Wang,1997:Sensible heat driven air-pump over the Tibetan Plateau and its impacts on the Asian summer monsoon.Collections on the Memory of Zhao Jiuzhang,D.Z.Ye,Ed.,Chinese Science Press,116–126.(in Chinese)

    Wu,G.X.,B.He,Y.M.Liu,Q.Bao,and R.C.Ren,2015:Location and variation of the summertime upper-troposphere temperature maximum over South Asia.Climate Dyn.,45,2757–2774,https://doi.org/10.1007/s00382-015-2506-4.

    Wu,G.X.,X.Liu,Q.Zhang,Y.F.Qian,J.Y.Mao,Y.M.Liu,and W.P.Li,2002:Progresses in the study of the climate impacts of the elevated heating over the Tibetan Plateau.Climatic and Environmental Research,7,184–201,https://doi.org/10.3878/j.issn.1006-9585.2002.02.06.(in Chinese)

    Wu,G.X.,Y.Liu,X.Zhu,W.Li,R.Ren,A.Duan,and X.Y.Liang,2009:Multi-scale forcing and the formation of subtropical desert and monsoon.Annales Geophysicae,27,3631–3644,https://doi.org/10.5194/angeo-27-3631-2009.

    Wu,G.X.,Y.M.Liu,B.He,Q.Bao,A.M.Duan,and F.F.Jin,2012:Thermal controls on the Asian summer monsoon.Scienti fi c Reports,2,404,https://doi.org/10.1038/srep00404.

    Wu,G.X.,and Coauthors,2007:The in fl uence of mechanical and thermal forcing by the Tibetan Plateau on Asian Climate.Journal of Hydrometeorology,8,770–789,https://doi.org/10.1175/JHM609.1.

    Yanai,M.,and G.X.Wu,2006:E ff ects of the Tibetan Plateau.Chap.13,Asian Monsoon,B.Wang,Ed.,Springer,Chichester,513–549.

    Yanai,M.,C.F.Li,and Z.S.Song,1992:Seasonal heating of the Tibetan Plateau and its e ff ects on the evolution of the Asian Summer monsoon.J.Meteor.Soc.Japan,70,319–351,https://doi.org/10.2151/jmsj1965.70.1B_319.

    Yang,J.,Q.Bao,B.Wang,H.Z.He,M.N.Gao,and D.Y.Gong,2016:Characterizing two types of transient intraseasonal oscillations in the Eastern Tibetan Plateau summer rainfall.Climate Dyn.,48,1749–1768,https://doi.org/10.1007/s00382-016-3170-z.

    Yang,W.Y.,D.Z.Ye,and G.X.Wu,1992:The in fl uence of the Tibetan Plateau on the summer thermal and circulation fi elds over East Asia I.The Humidity on the western Tibetan Plateau in the Height of Summer.Chinese Journal of Atmospheric Sciences,16,41–51,https://doi.org/10.3878/j.issn.1006-9895.1992.01.06.(in Chinese)

    Ye,D.Z.,and Y.X.Gao,1979:Meteorology of Tibetan Plateau.Science Press,278 pp.(in Chinese)

    Ye,D.Z.,and G.J.Yang,1979:The Average Vertical Circulation over the Qinghai-Xizang Plateau.Science Press.(in Chinese)

    Ye,D.Z.,and G.X.Wu,1998:The role of the heat source of the Tibetan Plateau in the general circulation.Meteor.Atmos.Phys.,67,181–198,https://doi.org/10.1007/BF01277509.

    Zhang,Y.,L.X.Chen,J.H.He,and W.Li,2009:A study of the characteristics of the low-frequency circulation over the Tibetan Plateau and its association with precipitation in the Yangtze River valley in 1998.Acta Meteorologica Sinica,23,175–190.

    Zhao,P.,and L.X.Chen,2001:Interannual variability of atmospheric heat source/sink over the Qinghai-Xizang(Tibetan)Plateau and its relation to circulation.Adv.Atmos.Sci.,18,106–116,https://doi.org/10.1007/s00376-001-0007-3.

    国产三级在线视频| 国产熟女午夜一区二区三区| av片东京热男人的天堂| 一卡2卡三卡四卡精品乱码亚洲| 久9热在线精品视频| 后天国语完整版免费观看| 久久精品国产亚洲av高清一级| 在线免费观看的www视频| 亚洲精品国产一区二区精华液| 久久天躁狠狠躁夜夜2o2o| 午夜福利高清视频| 久久香蕉精品热| 制服诱惑二区| 中文字幕av电影在线播放| 国产精品久久电影中文字幕| 欧美日韩亚洲综合一区二区三区_| 欧美av亚洲av综合av国产av| 亚洲一区中文字幕在线| 老司机福利观看| 两性夫妻黄色片| 亚洲一码二码三码区别大吗| 久久天堂一区二区三区四区| 此物有八面人人有两片| 国产成人精品久久二区二区免费| www国产在线视频色| 中文字幕精品亚洲无线码一区 | 国产一区在线观看成人免费| 热99re8久久精品国产| 黄色成人免费大全| tocl精华| 亚洲av熟女| 国产精品,欧美在线| 亚洲精品中文字幕在线视频| 欧美绝顶高潮抽搐喷水| 亚洲天堂国产精品一区在线| 久久这里只有精品19| 搡老妇女老女人老熟妇| 午夜福利18| 哪里可以看免费的av片| cao死你这个sao货| 免费看十八禁软件| 久久热在线av| 香蕉久久夜色| 美女免费视频网站| 久久精品影院6| 在线观看午夜福利视频| 国产97色在线日韩免费| 久久国产精品人妻蜜桃| 黄色成人免费大全| 18禁黄网站禁片免费观看直播| 免费电影在线观看免费观看| 欧美性猛交╳xxx乱大交人| 久久中文字幕一级| 一区二区三区精品91| 中文字幕精品亚洲无线码一区 | 中出人妻视频一区二区| 欧美精品啪啪一区二区三区| 亚洲专区中文字幕在线| 在线永久观看黄色视频| 老鸭窝网址在线观看| 午夜精品久久久久久毛片777| 国产视频内射| 免费观看精品视频网站| 妹子高潮喷水视频| avwww免费| 啦啦啦免费观看视频1| 久久人妻福利社区极品人妻图片| 欧美色视频一区免费| 中文在线观看免费www的网站 | 99久久综合精品五月天人人| 欧美精品啪啪一区二区三区| 亚洲精品一区av在线观看| 成人亚洲精品一区在线观看| 大香蕉久久成人网| 国产aⅴ精品一区二区三区波| 久久久久久久久中文| 欧美乱码精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品九九99| 久久午夜亚洲精品久久| 亚洲国产精品久久男人天堂| 国产精品1区2区在线观看.| 日本免费一区二区三区高清不卡| 亚洲第一av免费看| 精品久久久久久久末码| 狂野欧美激情性xxxx| 美国免费a级毛片| 色综合婷婷激情| 又黄又爽又免费观看的视频| 久久精品国产亚洲av香蕉五月| 国产精品日韩av在线免费观看| 好男人在线观看高清免费视频 | 久久国产乱子伦精品免费另类| 日韩欧美一区二区三区在线观看| 日韩视频一区二区在线观看| 久久久国产成人精品二区| 亚洲欧洲精品一区二区精品久久久| 亚洲一区中文字幕在线| 欧美av亚洲av综合av国产av| 色综合亚洲欧美另类图片| 人成视频在线观看免费观看| 最新在线观看一区二区三区| 欧美日韩瑟瑟在线播放| 国产男靠女视频免费网站| 淫秽高清视频在线观看| 中文字幕人妻熟女乱码| 国产v大片淫在线免费观看| 国产成人av激情在线播放| 欧美精品亚洲一区二区| 欧美国产日韩亚洲一区| 国产精品野战在线观看| 欧美又色又爽又黄视频| av视频在线观看入口| 99热这里只有精品一区 | 午夜精品在线福利| 成年女人毛片免费观看观看9| 亚洲一区二区三区不卡视频| 亚洲av成人不卡在线观看播放网| 亚洲狠狠婷婷综合久久图片| www.www免费av| 黄色片一级片一级黄色片| 高清在线国产一区| 少妇粗大呻吟视频| 免费高清视频大片| 在线观看免费午夜福利视频| 美女扒开内裤让男人捅视频| 久久国产精品影院| 88av欧美| 亚洲国产精品成人综合色| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品合色在线| 国产精品,欧美在线| 亚洲一码二码三码区别大吗| 丰满人妻熟妇乱又伦精品不卡| 91九色精品人成在线观看| 国产精品综合久久久久久久免费| 波多野结衣高清无吗| 精品久久久久久,| 亚洲精品在线观看二区| 亚洲全国av大片| 变态另类成人亚洲欧美熟女| 777久久人妻少妇嫩草av网站| 不卡av一区二区三区| 好看av亚洲va欧美ⅴa在| 黑人欧美特级aaaaaa片| 搡老妇女老女人老熟妇| 琪琪午夜伦伦电影理论片6080| 久久久精品欧美日韩精品| 国产91精品成人一区二区三区| 国产精品免费视频内射| 一级黄色大片毛片| 国产黄片美女视频| 757午夜福利合集在线观看| 国产亚洲精品久久久久久毛片| 在线观看66精品国产| 成人国产综合亚洲| 国产精品98久久久久久宅男小说| 午夜视频精品福利| 亚洲全国av大片| 女人高潮潮喷娇喘18禁视频| 亚洲国产欧美日韩在线播放| 一进一出好大好爽视频| 天堂动漫精品| 亚洲九九香蕉| 久久中文字幕人妻熟女| 老司机深夜福利视频在线观看| 一二三四在线观看免费中文在| 国产精品亚洲一级av第二区| 超碰成人久久| 欧美日韩精品网址| 成人三级黄色视频| 无遮挡黄片免费观看| 国产高清视频在线播放一区| 国产又色又爽无遮挡免费看| 男女做爰动态图高潮gif福利片| 久久99热这里只有精品18| 亚洲无线在线观看| 国产伦人伦偷精品视频| 香蕉丝袜av| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕人妻丝袜一区二区| 成人av一区二区三区在线看| 国产成人系列免费观看| 成人特级黄色片久久久久久久| 亚洲av成人不卡在线观看播放网| 日韩精品青青久久久久久| 美女大奶头视频| 精品人妻1区二区| 亚洲电影在线观看av| 亚洲成人精品中文字幕电影| 久久久国产欧美日韩av| 97碰自拍视频| 久久天堂一区二区三区四区| 男女床上黄色一级片免费看| 亚洲激情在线av| 亚洲精品美女久久av网站| 黄色视频不卡| 亚洲精华国产精华精| 久久久久九九精品影院| 欧美成人一区二区免费高清观看 | 日韩欧美国产在线观看| 黄网站色视频无遮挡免费观看| 精品国产美女av久久久久小说| 久久久国产成人免费| 国内精品久久久久精免费| 亚洲一区二区三区不卡视频| 精品国产美女av久久久久小说| 国产亚洲精品久久久久5区| 午夜成年电影在线免费观看| 精品少妇一区二区三区视频日本电影| 国产精品av久久久久免费| 国产成人av教育| 免费在线观看视频国产中文字幕亚洲| 国产精品香港三级国产av潘金莲| 久久香蕉国产精品| 一本一本综合久久| 日韩精品中文字幕看吧| 国产一区二区激情短视频| 99久久精品国产亚洲精品| 午夜激情av网站| 国内少妇人妻偷人精品xxx网站 | 老熟妇乱子伦视频在线观看| 中文字幕人妻丝袜一区二区| 精品无人区乱码1区二区| 免费观看精品视频网站| 国产乱人伦免费视频| 嫩草影视91久久| 国产精品久久久久久精品电影 | 亚洲人成网站在线播放欧美日韩| 国产私拍福利视频在线观看| 桃红色精品国产亚洲av| 99国产精品99久久久久| 91成人精品电影| 夜夜爽天天搞| 精品久久久久久成人av| 亚洲av成人一区二区三| 亚洲,欧美精品.| 精品欧美国产一区二区三| 中文字幕精品免费在线观看视频| 1024香蕉在线观看| 首页视频小说图片口味搜索| 日本熟妇午夜| 一区二区日韩欧美中文字幕| 99在线人妻在线中文字幕| 国产99白浆流出| 看免费av毛片| 国产又黄又爽又无遮挡在线| 国内毛片毛片毛片毛片毛片| 亚洲成av片中文字幕在线观看| av视频在线观看入口| a在线观看视频网站| 国产麻豆成人av免费视频| 欧美一级毛片孕妇| 亚洲av片天天在线观看| 激情在线观看视频在线高清| ponron亚洲| 99热只有精品国产| 国内毛片毛片毛片毛片毛片| 一区二区三区精品91| 国产精华一区二区三区| 欧美最黄视频在线播放免费| 极品教师在线免费播放| 欧美日韩亚洲综合一区二区三区_| 99国产精品一区二区蜜桃av| 男人的好看免费观看在线视频 | 亚洲一区二区三区不卡视频| 久久国产精品影院| 久久精品亚洲精品国产色婷小说| 国产精品亚洲一级av第二区| 精品国产美女av久久久久小说| 国产亚洲精品久久久久5区| 中文字幕人妻熟女乱码| 国产成人av激情在线播放| 成年版毛片免费区| 熟妇人妻久久中文字幕3abv| 国产1区2区3区精品| 天天一区二区日本电影三级| 人人妻人人澡欧美一区二区| 欧美日韩亚洲综合一区二区三区_| 久久久久久亚洲精品国产蜜桃av| 国产久久久一区二区三区| 亚洲av五月六月丁香网| 在线永久观看黄色视频| 18美女黄网站色大片免费观看| 首页视频小说图片口味搜索| 黄色毛片三级朝国网站| 亚洲色图av天堂| 欧美一区二区精品小视频在线| 妹子高潮喷水视频| 国产私拍福利视频在线观看| 俄罗斯特黄特色一大片| 十分钟在线观看高清视频www| 成人永久免费在线观看视频| 亚洲三区欧美一区| 国产真人三级小视频在线观看| 亚洲男人天堂网一区| 在线国产一区二区在线| 日日爽夜夜爽网站| 天堂影院成人在线观看| 淫妇啪啪啪对白视频| 在线观看舔阴道视频| 日本熟妇午夜| 黄片大片在线免费观看| 岛国视频午夜一区免费看| 国产成年人精品一区二区| 欧美激情高清一区二区三区| 老鸭窝网址在线观看| 亚洲色图 男人天堂 中文字幕| 国产不卡一卡二| 日本免费一区二区三区高清不卡| 欧美日韩精品网址| 日韩欧美一区视频在线观看| 精品少妇一区二区三区视频日本电影| 国产亚洲精品综合一区在线观看 | 国产精品久久电影中文字幕| 人人妻,人人澡人人爽秒播| 欧美日韩中文字幕国产精品一区二区三区| 黄色视频不卡| tocl精华| 久久精品国产99精品国产亚洲性色| 99久久精品国产亚洲精品| 久久人妻av系列| 一本精品99久久精品77| 国产片内射在线| 久久草成人影院| 黄片大片在线免费观看| 午夜福利视频1000在线观看| 久久精品成人免费网站| 变态另类丝袜制服| 亚洲精华国产精华精| 丝袜人妻中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 欧美黑人精品巨大| 天天躁狠狠躁夜夜躁狠狠躁| 搡老岳熟女国产| 成人午夜高清在线视频 | 免费av毛片视频| 精品欧美国产一区二区三| 麻豆成人av在线观看| 在线天堂中文资源库| 黑人欧美特级aaaaaa片| 韩国精品一区二区三区| 成人午夜高清在线视频 | 91字幕亚洲| 久久久久久国产a免费观看| 亚洲无线在线观看| 亚洲午夜精品一区,二区,三区| 国产三级在线视频| 成人国产一区最新在线观看| 久久婷婷成人综合色麻豆| 18禁观看日本| 免费电影在线观看免费观看| 久久亚洲真实| 久久香蕉激情| 在线播放国产精品三级| 国产激情欧美一区二区| 999久久久精品免费观看国产| 欧美激情 高清一区二区三区| 狠狠狠狠99中文字幕| 亚洲免费av在线视频| 亚洲国产精品久久男人天堂| 精品国产亚洲在线| 国产精品一区二区三区四区久久 | 国产伦人伦偷精品视频| 一本综合久久免费| 丰满的人妻完整版| 亚洲第一电影网av| 久久人人精品亚洲av| 国产国语露脸激情在线看| 757午夜福利合集在线观看| 黄频高清免费视频| 精品国产乱子伦一区二区三区| 在线观看一区二区三区| 国产精品美女特级片免费视频播放器 | 国产伦在线观看视频一区| 欧美黄色片欧美黄色片| 成人18禁高潮啪啪吃奶动态图| 国产av在哪里看| 国产午夜精品久久久久久| 国产99久久九九免费精品| 日本 av在线| 亚洲精品久久国产高清桃花| 真人做人爱边吃奶动态| 黄片小视频在线播放| 午夜久久久久精精品| 人人妻,人人澡人人爽秒播| tocl精华| 老汉色∧v一级毛片| 免费人成视频x8x8入口观看| 久久国产精品影院| 满18在线观看网站| 国产精品乱码一区二三区的特点| 国产精品日韩av在线免费观看| 国产成人精品久久二区二区免费| 中文字幕人妻丝袜一区二区| 国产高清videossex| 久久国产精品影院| 麻豆久久精品国产亚洲av| 欧美精品亚洲一区二区| 搡老岳熟女国产| 中文字幕另类日韩欧美亚洲嫩草| 俄罗斯特黄特色一大片| 露出奶头的视频| 亚洲第一欧美日韩一区二区三区| 女性被躁到高潮视频| 18禁国产床啪视频网站| 成人手机av| 高清在线国产一区| 好男人电影高清在线观看| 女生性感内裤真人,穿戴方法视频| 欧美中文日本在线观看视频| 成人精品一区二区免费| 国产高清激情床上av| 亚洲国产欧美一区二区综合| 一卡2卡三卡四卡精品乱码亚洲| 每晚都被弄得嗷嗷叫到高潮| 91麻豆精品激情在线观看国产| av福利片在线| 亚洲一区中文字幕在线| 亚洲狠狠婷婷综合久久图片| 国产v大片淫在线免费观看| 在线视频色国产色| 精品一区二区三区四区五区乱码| 一级a爱片免费观看的视频| 可以免费在线观看a视频的电影网站| 日韩精品中文字幕看吧| 男女之事视频高清在线观看| 曰老女人黄片| 黄频高清免费视频| 国产精品98久久久久久宅男小说| 成在线人永久免费视频| 色婷婷久久久亚洲欧美| 一个人免费在线观看的高清视频| 香蕉久久夜色| АⅤ资源中文在线天堂| 黑人操中国人逼视频| 午夜福利18| 可以免费在线观看a视频的电影网站| 亚洲精品美女久久av网站| 日韩欧美国产在线观看| 美女高潮喷水抽搐中文字幕| 正在播放国产对白刺激| 亚洲精品美女久久av网站| 国内精品久久久久精免费| 国产精品九九99| 日本一本二区三区精品| 日韩三级视频一区二区三区| 黄片播放在线免费| 国产一区二区在线av高清观看| 久久热在线av| 变态另类丝袜制服| 50天的宝宝边吃奶边哭怎么回事| 国产精品爽爽va在线观看网站 | 久久久久精品国产欧美久久久| 午夜免费观看网址| 欧美国产日韩亚洲一区| 亚洲熟妇中文字幕五十中出| 少妇的丰满在线观看| 久久亚洲精品不卡| 人妻久久中文字幕网| 免费av毛片视频| 午夜免费鲁丝| 精品国产乱码久久久久久男人| 欧美性长视频在线观看| 精品久久久久久久末码| 国产爱豆传媒在线观看 | 黄片播放在线免费| 亚洲av电影在线进入| 国产精品乱码一区二三区的特点| 男人的好看免费观看在线视频 | 午夜福利免费观看在线| 欧美精品啪啪一区二区三区| 国产亚洲av嫩草精品影院| 午夜福利高清视频| 久久精品国产99精品国产亚洲性色| av欧美777| 国产精品一区二区免费欧美| 99国产综合亚洲精品| 亚洲国产欧美网| 亚洲精品国产精品久久久不卡| 亚洲天堂国产精品一区在线| 美女午夜性视频免费| 在线观看免费日韩欧美大片| 午夜精品久久久久久毛片777| 久久精品国产综合久久久| 亚洲av美国av| 亚洲无线在线观看| 国产成人欧美在线观看| 日韩免费av在线播放| 国产精品久久久久久人妻精品电影| 国产精品野战在线观看| 一级a爱片免费观看的视频| 日韩大码丰满熟妇| 手机成人av网站| АⅤ资源中文在线天堂| 脱女人内裤的视频| 丁香欧美五月| www.www免费av| 一级黄色大片毛片| 91成年电影在线观看| 这个男人来自地球电影免费观看| svipshipincom国产片| 国产亚洲av高清不卡| 亚洲无线在线观看| 精品久久久久久成人av| 12—13女人毛片做爰片一| 1024手机看黄色片| 在线国产一区二区在线| 国产不卡一卡二| 精品午夜福利视频在线观看一区| 中文亚洲av片在线观看爽| 又大又爽又粗| 99久久99久久久精品蜜桃| 变态另类丝袜制服| 无限看片的www在线观看| 哪里可以看免费的av片| 99国产精品一区二区三区| 黄色视频不卡| 亚洲精品av麻豆狂野| 亚洲七黄色美女视频| 国产免费男女视频| 俺也久久电影网| 国产亚洲精品久久久久久毛片| 最近在线观看免费完整版| 老熟妇仑乱视频hdxx| 久久亚洲真实| 国产精品免费视频内射| 欧美成人一区二区免费高清观看 | 国产成人欧美在线观看| 亚洲av成人不卡在线观看播放网| 精品午夜福利视频在线观看一区| 亚洲国产看品久久| 欧美最黄视频在线播放免费| 久久香蕉国产精品| 婷婷亚洲欧美| 国产伦人伦偷精品视频| 久久久久久免费高清国产稀缺| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 香蕉丝袜av| 99国产精品99久久久久| 一级片免费观看大全| 白带黄色成豆腐渣| 久9热在线精品视频| 久热爱精品视频在线9| 在线观看午夜福利视频| 日韩欧美国产一区二区入口| 97超级碰碰碰精品色视频在线观看| 亚洲片人在线观看| 亚洲欧洲精品一区二区精品久久久| videosex国产| 国产久久久一区二区三区| 亚洲精品美女久久av网站| 日韩大尺度精品在线看网址| 免费观看人在逋| 国产1区2区3区精品| 琪琪午夜伦伦电影理论片6080| 白带黄色成豆腐渣| 国产成人影院久久av| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一小说| www.自偷自拍.com| 久久天躁狠狠躁夜夜2o2o| 亚洲专区中文字幕在线| 97超级碰碰碰精品色视频在线观看| 成人国产一区最新在线观看| 精品无人区乱码1区二区| 自线自在国产av| 成人三级做爰电影| 精品国产一区二区三区四区第35| 九色国产91popny在线| 天堂√8在线中文| 国产精品国产高清国产av| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产一区二区入口| 首页视频小说图片口味搜索| www.999成人在线观看| 可以在线观看毛片的网站| 满18在线观看网站| 757午夜福利合集在线观看| 99精品久久久久人妻精品| 亚洲自拍偷在线| 国产久久久一区二区三区| videosex国产| 无限看片的www在线观看| 一级作爱视频免费观看| av天堂在线播放| 男女那种视频在线观看| 岛国视频午夜一区免费看| 亚洲va日本ⅴa欧美va伊人久久| 女生性感内裤真人,穿戴方法视频| 国产黄片美女视频| 国产精品久久久久久人妻精品电影| 亚洲精品一卡2卡三卡4卡5卡| 老司机午夜十八禁免费视频| 国产又色又爽无遮挡免费看| 成人精品一区二区免费| 国产亚洲精品久久久久5区| 悠悠久久av| 一级作爱视频免费观看| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 美女高潮喷水抽搐中文字幕| 熟妇人妻久久中文字幕3abv| 中国美女看黄片| 波多野结衣av一区二区av| 久久草成人影院| 国产三级黄色录像| 国产亚洲精品综合一区在线观看 | 久久狼人影院| 看黄色毛片网站| 亚洲欧美一区二区三区黑人| 12—13女人毛片做爰片一| av片东京热男人的天堂| 一二三四社区在线视频社区8|