• <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片在线观看爽| 又黄又爽又免费观看的视频| а√天堂www在线а√下载| 国产精品一区二区性色av| 欧美另类亚洲清纯唯美| 亚洲av.av天堂| 欧美日本视频| or卡值多少钱| 三级男女做爰猛烈吃奶摸视频| 日韩 亚洲 欧美在线| 成人国产综合亚洲| 成人av在线播放网站| 亚洲精品色激情综合| 老熟妇仑乱视频hdxx| 麻豆一二三区av精品| 一二三四社区在线视频社区8| 欧美性感艳星| 男人舔奶头视频| av在线老鸭窝| 一区二区三区免费毛片| 亚洲精品成人久久久久久| 桃色一区二区三区在线观看| 一进一出好大好爽视频| 亚洲三级黄色毛片| 亚洲av日韩精品久久久久久密| 日韩大尺度精品在线看网址| 97热精品久久久久久| 美女免费视频网站| 国产精品不卡视频一区二区 | 久久久久免费精品人妻一区二区| 久久久精品大字幕| 男女之事视频高清在线观看| 免费大片18禁| 欧美性感艳星| 久久天躁狠狠躁夜夜2o2o| 真人一进一出gif抽搐免费| 精品一区二区三区视频在线| 在线免费观看不下载黄p国产 | 老司机午夜十八禁免费视频| 99久国产av精品| 天堂影院成人在线观看| 精品国内亚洲2022精品成人| 亚洲成av人片在线播放无| 国产69精品久久久久777片| 观看美女的网站| 色播亚洲综合网| 99riav亚洲国产免费| av女优亚洲男人天堂| 99热这里只有是精品在线观看 | 1000部很黄的大片| av在线观看视频网站免费| 欧美在线黄色| 深夜精品福利| av在线天堂中文字幕| 色播亚洲综合网| 一级作爱视频免费观看| 蜜桃亚洲精品一区二区三区| 乱码一卡2卡4卡精品| 少妇的逼水好多| 又黄又爽又免费观看的视频| 国产极品精品免费视频能看的| 大型黄色视频在线免费观看| 亚洲国产精品999在线| 国产精品久久久久久精品电影| 天堂动漫精品| 国产精品精品国产色婷婷| 精品久久久久久久末码| 日韩欧美免费精品| 国产爱豆传媒在线观看| 成人毛片a级毛片在线播放| 亚洲精品色激情综合| 中文字幕熟女人妻在线| 日韩av在线大香蕉| 给我免费播放毛片高清在线观看| 一进一出抽搐动态| 亚洲av不卡在线观看| 老女人水多毛片| 国内精品久久久久精免费| 国产成+人综合+亚洲专区| 亚洲专区国产一区二区| av在线蜜桃| 日韩成人在线观看一区二区三区| 宅男免费午夜| 亚洲性夜色夜夜综合| 偷拍熟女少妇极品色| 午夜福利在线观看吧| 国产精品美女特级片免费视频播放器| 亚洲第一欧美日韩一区二区三区| 美女高潮的动态| 男人狂女人下面高潮的视频| 99久久成人亚洲精品观看| a级一级毛片免费在线观看| 欧美不卡视频在线免费观看| 亚洲成人久久性| 亚洲精品成人久久久久久| АⅤ资源中文在线天堂| av视频在线观看入口| 久久久国产成人免费| 亚洲精品456在线播放app | 国产av一区在线观看免费| 成人精品一区二区免费| 亚洲av五月六月丁香网| 好男人电影高清在线观看| 热99在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 免费搜索国产男女视频| 免费人成视频x8x8入口观看| 国内少妇人妻偷人精品xxx网站| 我的老师免费观看完整版| 亚洲国产高清在线一区二区三| 午夜影院日韩av| 国产精品,欧美在线| 免费大片18禁| 亚洲av免费高清在线观看| 一级av片app| 人妻丰满熟妇av一区二区三区| 欧美日本视频| 九九在线视频观看精品| 精品国产亚洲在线| 动漫黄色视频在线观看| 色在线成人网| 99热只有精品国产| 一夜夜www| 国产真实乱freesex| 色尼玛亚洲综合影院| netflix在线观看网站| 精品国产亚洲在线| 最新中文字幕久久久久| 中文字幕熟女人妻在线| 一本一本综合久久| 欧洲精品卡2卡3卡4卡5卡区| netflix在线观看网站| 三级国产精品欧美在线观看| 热99在线观看视频| 欧美性猛交黑人性爽| 一区二区三区四区激情视频 | 亚洲欧美精品综合久久99| 成人精品一区二区免费| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| a级毛片a级免费在线| 久久午夜亚洲精品久久| 欧洲精品卡2卡3卡4卡5卡区| 老司机深夜福利视频在线观看| 一进一出抽搐gif免费好疼| 免费av毛片视频| 51午夜福利影视在线观看| 国产成人av教育| 午夜精品一区二区三区免费看| .国产精品久久| av国产免费在线观看| 欧美性猛交黑人性爽| 午夜精品在线福利| www.熟女人妻精品国产| 国产精品乱码一区二三区的特点| 91午夜精品亚洲一区二区三区 | 琪琪午夜伦伦电影理论片6080| 亚洲最大成人手机在线| 久久伊人香网站| 伊人久久精品亚洲午夜| 亚洲精品一区av在线观看| 久久久久久久午夜电影| 国产亚洲精品久久久com| 日韩欧美在线二视频| 亚洲欧美日韩东京热| 色在线成人网| 一级毛片久久久久久久久女| 国内毛片毛片毛片毛片毛片| 少妇丰满av| 嫩草影院精品99| 国产精品美女特级片免费视频播放器| 十八禁国产超污无遮挡网站| 欧美最新免费一区二区三区 | 淫秽高清视频在线观看| 别揉我奶头~嗯~啊~动态视频| 九色国产91popny在线| 国产精品亚洲一级av第二区| 日韩中字成人| 深爱激情五月婷婷| 国产大屁股一区二区在线视频| 国产激情偷乱视频一区二区| 国产亚洲精品久久久久久毛片| 精品久久久久久久久亚洲 | 中亚洲国语对白在线视频| 亚洲欧美精品综合久久99| 国产伦在线观看视频一区| 久久精品人妻少妇| 国产av不卡久久| 日韩亚洲欧美综合| 国产高清三级在线| av福利片在线观看| 国产一区二区在线观看日韩| 一级毛片久久久久久久久女| 免费在线观看成人毛片| 午夜视频国产福利| 国产精品综合久久久久久久免费| 亚洲av美国av| 18禁裸乳无遮挡免费网站照片| 欧美高清性xxxxhd video| 黄色一级大片看看| 欧美+亚洲+日韩+国产| 国产精品一区二区三区四区久久| 最后的刺客免费高清国语| 国产成人啪精品午夜网站| 亚洲最大成人中文| 午夜福利18| 国内毛片毛片毛片毛片毛片| 两人在一起打扑克的视频| 久久久久久久久中文| 又粗又爽又猛毛片免费看| 亚洲内射少妇av| 欧美在线一区亚洲| 色av中文字幕| 成人高潮视频无遮挡免费网站| 国产精品影院久久| 亚洲最大成人av| 老司机午夜十八禁免费视频| 国产亚洲精品av在线| 国产男靠女视频免费网站| 亚洲欧美日韩无卡精品| 国产精品1区2区在线观看.| 午夜影院日韩av| 久久久久久国产a免费观看| 亚洲人成电影免费在线| 欧美最黄视频在线播放免费| 999久久久精品免费观看国产| 丰满人妻一区二区三区视频av| 久久午夜福利片| 丰满人妻熟妇乱又伦精品不卡| 激情在线观看视频在线高清| av国产免费在线观看| 亚洲欧美日韩卡通动漫| 亚洲av电影不卡..在线观看| 少妇熟女aⅴ在线视频| 成人三级黄色视频| 男女之事视频高清在线观看| 91麻豆av在线| 婷婷亚洲欧美| 国产免费av片在线观看野外av| 午夜免费激情av| 精品一区二区三区av网在线观看| 久久精品91蜜桃| 日本精品一区二区三区蜜桃| 午夜精品在线福利| 99热精品在线国产| 日本免费一区二区三区高清不卡| 色综合站精品国产| 在现免费观看毛片| 免费人成在线观看视频色| 成人毛片a级毛片在线播放| 成年女人永久免费观看视频| 非洲黑人性xxxx精品又粗又长| 91九色精品人成在线观看| 人人妻人人澡欧美一区二区| 在线观看美女被高潮喷水网站 | 日日摸夜夜添夜夜添小说| 亚洲专区中文字幕在线| 国产黄a三级三级三级人| 亚洲人成电影免费在线| 亚洲av中文字字幕乱码综合| 婷婷精品国产亚洲av在线| 观看美女的网站| 亚洲久久久久久中文字幕| 波野结衣二区三区在线| 97超视频在线观看视频| 久久久色成人| 赤兔流量卡办理| 美女大奶头视频| 久久亚洲真实| 亚洲在线观看片| 国产黄片美女视频| 午夜视频国产福利| 欧美又色又爽又黄视频| 久久久精品欧美日韩精品| 老司机午夜十八禁免费视频| 老司机深夜福利视频在线观看| 久久精品久久久久久噜噜老黄 | 亚洲激情在线av| 国产欧美日韩一区二区精品| 日日夜夜操网爽| 99国产极品粉嫩在线观看| 精品免费久久久久久久清纯| 99热这里只有精品一区| 国产精品亚洲一级av第二区| 国产国拍精品亚洲av在线观看| x7x7x7水蜜桃| 一个人看视频在线观看www免费| 日韩中字成人| 男人的好看免费观看在线视频| 可以在线观看毛片的网站| 自拍偷自拍亚洲精品老妇| 天美传媒精品一区二区| 亚洲不卡免费看| 日韩中文字幕欧美一区二区| 色哟哟·www| 如何舔出高潮| 国产精品久久久久久久电影| 国产欧美日韩精品亚洲av| 在线看三级毛片| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件 | 欧美三级亚洲精品| 国产精品日韩av在线免费观看| 国产日本99.免费观看| 看十八女毛片水多多多| 哪里可以看免费的av片| 99riav亚洲国产免费| 亚洲av中文字字幕乱码综合| 无遮挡黄片免费观看| 亚洲精品影视一区二区三区av| 别揉我奶头 嗯啊视频| 999久久久精品免费观看国产| 综合色av麻豆| 色综合欧美亚洲国产小说| 在线观看66精品国产| 国产高清视频在线观看网站| 精品无人区乱码1区二区| 国内精品一区二区在线观看| 国产精品免费一区二区三区在线| bbb黄色大片| 在线十欧美十亚洲十日本专区| 丰满的人妻完整版| 一本一本综合久久| 国产伦在线观看视频一区| 九九久久精品国产亚洲av麻豆| 1024手机看黄色片| 国产伦一二天堂av在线观看| 国产伦精品一区二区三区视频9| 给我免费播放毛片高清在线观看| 国产精品人妻久久久久久| av在线观看视频网站免费| 51午夜福利影视在线观看| 中文在线观看免费www的网站| 大型黄色视频在线免费观看| 女人十人毛片免费观看3o分钟| 精品一区二区三区av网在线观看| 国产男靠女视频免费网站| 日本 欧美在线| 国产精品电影一区二区三区| 美女免费视频网站| 午夜精品在线福利| 国产精品99久久久久久久久| 国产伦精品一区二区三区视频9| 国产精品嫩草影院av在线观看 | 国产久久久一区二区三区| 九九在线视频观看精品| 舔av片在线| 久久久久精品国产欧美久久久| 久久热精品热| 老鸭窝网址在线观看| 久久人妻av系列| 日韩欧美精品v在线| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区三区四区久久| 国产欧美日韩一区二区精品| 小蜜桃在线观看免费完整版高清| 最新中文字幕久久久久| 精品熟女少妇八av免费久了| 久久久久久久久久成人| 免费电影在线观看免费观看| 在线看三级毛片| 午夜两性在线视频| 国内毛片毛片毛片毛片毛片| 国产aⅴ精品一区二区三区波| 国产爱豆传媒在线观看| 亚洲黑人精品在线| 国产在视频线在精品| 精品人妻1区二区| 俺也久久电影网| 亚洲,欧美,日韩| 精品久久国产蜜桃| 免费黄网站久久成人精品 | 色精品久久人妻99蜜桃| 精品久久久久久久久久免费视频| av中文乱码字幕在线| 日本一二三区视频观看| 99国产精品一区二区三区| 国产精品三级大全| 国产久久久一区二区三区| 日韩欧美国产一区二区入口| 午夜老司机福利剧场| 亚洲午夜理论影院| 97超视频在线观看视频| 蜜桃亚洲精品一区二区三区| 免费在线观看影片大全网站| 久久6这里有精品| 国产欧美日韩精品一区二区| 久久这里只有精品中国| 狂野欧美白嫩少妇大欣赏| 伦理电影大哥的女人| 欧美在线黄色| 国产精品久久视频播放| 日日摸夜夜添夜夜添av毛片 | www日本黄色视频网| 毛片女人毛片| 波多野结衣巨乳人妻| 亚洲电影在线观看av| 亚洲av一区综合| 啦啦啦观看免费观看视频高清| 男人舔女人下体高潮全视频| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 我要搜黄色片| 美女被艹到高潮喷水动态| 一进一出好大好爽视频| 一级av片app| 日韩人妻高清精品专区| 欧美黑人巨大hd| 午夜免费男女啪啪视频观看 | 特级一级黄色大片| 亚洲av成人av| 亚洲男人的天堂狠狠| 久久欧美精品欧美久久欧美| 欧美激情在线99| 亚洲午夜理论影院| 天堂动漫精品| 色尼玛亚洲综合影院| 中文字幕人成人乱码亚洲影| 夜夜夜夜夜久久久久| 一进一出抽搐gif免费好疼| 精品一区二区免费观看| 长腿黑丝高跟| 国产精品久久久久久久久免 | 香蕉av资源在线| 国产亚洲欧美98| 一区二区三区高清视频在线| 91久久精品国产一区二区成人| 人人妻人人看人人澡| 日本精品一区二区三区蜜桃| 亚洲精品乱码久久久v下载方式| 99国产极品粉嫩在线观看| 欧美色视频一区免费| 国产黄片美女视频| 亚洲av一区综合| 赤兔流量卡办理| 亚洲国产欧洲综合997久久,| 欧美国产日韩亚洲一区| 丁香欧美五月| 亚洲人与动物交配视频| 国产极品精品免费视频能看的| 久久亚洲真实| 伊人久久精品亚洲午夜| 一个人免费在线观看的高清视频| 欧美一级a爱片免费观看看| www.色视频.com| 国产单亲对白刺激| 自拍偷自拍亚洲精品老妇| 亚洲精华国产精华精| 99精品久久久久人妻精品| 国产黄片美女视频| 女人被狂操c到高潮| 麻豆成人av在线观看| 最近视频中文字幕2019在线8| 人妻丰满熟妇av一区二区三区| 免费av毛片视频| 国产精品久久电影中文字幕| 国产野战对白在线观看| 欧美日韩福利视频一区二区| 欧美最新免费一区二区三区 | 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2| 一区二区三区免费毛片| 免费无遮挡裸体视频| 精品免费久久久久久久清纯| 人人妻人人看人人澡| 波多野结衣高清作品| 亚洲第一电影网av| 少妇被粗大猛烈的视频| 两人在一起打扑克的视频| 好男人电影高清在线观看| 日韩亚洲欧美综合| 亚洲第一区二区三区不卡| 18禁裸乳无遮挡免费网站照片| 国产精品98久久久久久宅男小说| 女人十人毛片免费观看3o分钟| 欧美乱色亚洲激情| 俄罗斯特黄特色一大片| 在线国产一区二区在线| АⅤ资源中文在线天堂| 99热这里只有是精品50| 国产精品一及| 日本三级黄在线观看| 国内揄拍国产精品人妻在线| 色综合欧美亚洲国产小说| 少妇人妻一区二区三区视频| 很黄的视频免费| 精品人妻偷拍中文字幕| 一级作爱视频免费观看| 国产精品一区二区性色av| 一本一本综合久久| 国产一区二区三区视频了| 久久性视频一级片| 国产一级毛片七仙女欲春2| 给我免费播放毛片高清在线观看| 麻豆国产av国片精品| 成年版毛片免费区| 亚洲av一区综合| 18+在线观看网站| 国产精品日韩av在线免费观看| 啦啦啦观看免费观看视频高清| 99热只有精品国产| 亚洲在线自拍视频| 天堂动漫精品| 观看免费一级毛片| 深夜精品福利| 怎么达到女性高潮| 美女高潮喷水抽搐中文字幕| 中文字幕人妻熟人妻熟丝袜美| 男女做爰动态图高潮gif福利片| 色综合欧美亚洲国产小说| 精华霜和精华液先用哪个| 女人被狂操c到高潮| 亚洲欧美日韩无卡精品| 欧美中文日本在线观看视频| 久久国产乱子伦精品免费另类| 亚洲欧美激情综合另类| 久久久久亚洲av毛片大全| 51国产日韩欧美| 欧美日韩综合久久久久久 | 精品人妻偷拍中文字幕| 中文字幕精品亚洲无线码一区| 在线观看av片永久免费下载| 精品熟女少妇八av免费久了| 国产爱豆传媒在线观看| 无人区码免费观看不卡| 嫩草影视91久久| 久久久久性生活片| 亚洲成人精品中文字幕电影| 91麻豆av在线| 亚洲片人在线观看| 丰满人妻熟妇乱又伦精品不卡| 香蕉av资源在线| 国产成人a区在线观看| 亚洲内射少妇av| 日日摸夜夜添夜夜添小说| 90打野战视频偷拍视频| 日日摸夜夜添夜夜添av毛片 | 国产成人影院久久av| 日本黄色视频三级网站网址| 真人一进一出gif抽搐免费| 国产爱豆传媒在线观看| 国内精品久久久久久久电影| 亚洲成人精品中文字幕电影| 成熟少妇高潮喷水视频| 丰满的人妻完整版| 国产成人a区在线观看| 国产午夜福利久久久久久| 亚洲熟妇中文字幕五十中出| 日日摸夜夜添夜夜添小说| 亚洲狠狠婷婷综合久久图片| 少妇的逼水好多| a级毛片免费高清观看在线播放| 韩国av一区二区三区四区| 欧美一区二区亚洲| 少妇高潮的动态图| a在线观看视频网站| 色哟哟哟哟哟哟| 无人区码免费观看不卡| 精品久久久久久,| 日韩欧美在线二视频| 国产私拍福利视频在线观看| 国产av麻豆久久久久久久| 国产精品爽爽va在线观看网站| 亚洲午夜理论影院| 91在线精品国自产拍蜜月| 国产毛片a区久久久久| 日本 欧美在线| 在线观看免费视频日本深夜| 国产亚洲av嫩草精品影院| av专区在线播放| 久久久久久九九精品二区国产| 18禁黄网站禁片午夜丰满| 国产三级在线视频| 亚洲欧美精品综合久久99| 黄色女人牲交| 真实男女啪啪啪动态图| 久久精品综合一区二区三区| 韩国av一区二区三区四区| 成人午夜高清在线视频| av天堂中文字幕网| 精品欧美国产一区二区三| 国产精品综合久久久久久久免费| 老司机午夜十八禁免费视频| 婷婷色综合大香蕉| 精品久久久久久久久久久久久| 国产欧美日韩精品亚洲av| 午夜日韩欧美国产| 观看美女的网站| av天堂中文字幕网| 精品一区二区三区视频在线| 国产av在哪里看| 欧美极品一区二区三区四区| 美女黄网站色视频| 麻豆一二三区av精品| 搡老妇女老女人老熟妇| 亚洲色图av天堂| 亚洲av免费高清在线观看| 婷婷精品国产亚洲av在线| 亚洲av电影不卡..在线观看| 在线免费观看的www视频| 国产欧美日韩精品一区二区| 一级毛片久久久久久久久女| 亚洲av二区三区四区| 日韩亚洲欧美综合| 欧美日本亚洲视频在线播放| 亚洲精品乱码久久久v下载方式| 五月玫瑰六月丁香| 在线播放无遮挡| 亚洲国产精品成人综合色| 日韩中文字幕欧美一区二区| 国产美女午夜福利| 精华霜和精华液先用哪个| 久久国产精品人妻蜜桃| 欧美一区二区国产精品久久精品| 日韩精品中文字幕看吧|