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

    Upper-Tropospheric Environment–Tropical Cyclone Interactions over the Western North Pacific:A Statistical Study

    2016-12-07 07:41:02YuKunQIANChangXiaLIANGZhuojianYUANShiqiuPENGJunjieWUandSihuaWANG5CenterforMonsoonandEnvironmentalResearchDepartmentofAtmosphericScienceSunYatsenUniversityGuangzhou5075
    Advances in Atmospheric Sciences 2016年5期

    Yu-Kun QIAN,Chang-Xia LIANG,Zhuojian YUAN,Shiqiu PENG,Junjie WU,and Sihua WANG5Center for Monsoon and Environmental Research/Department of Atmospheric Science,Sun Yat-sen University,Guangzhou 5075

    2State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology, Chinese Academy of Sciences,Guangzhou 510301

    3South China Sea Marine Prediction Center,South Oceanic Administration,Guangzhou 510300

    4Civil Aviation Flight University of China,Guanghan 618307

    5Luogang District Meteorological Bureau,Guangzhou 510530

    Upper-Tropospheric Environment–Tropical Cyclone Interactions over the Western North Pacific:A Statistical Study

    Yu-Kun QIAN?1,2,Chang-Xia LIANG3,Zhuojian YUAN1,Shiqiu PENG2,Junjie WU4,and Sihua WANG51Center for Monsoon and Environmental Research/Department of Atmospheric Science,Sun Yat-sen University,Guangzhou 510275

    2State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology, Chinese Academy of Sciences,Guangzhou 510301

    3South China Sea Marine Prediction Center,South Oceanic Administration,Guangzhou 510300

    4Civil Aviation Flight University of China,Guanghan 618307

    5Luogang District Meteorological Bureau,Guangzhou 510530

    Based on 25-year(1987–2011)tropical cyclone(TC)best track data,a statistical study was carried out to investigate the basic features of upper-tropospheric TC–environment interactions over the western North Pacific.Interaction was defined

    as the absolute value of eddy momentum flux convergence(EFC)exceeding 10 m s-1d-1.Based on this definition,it was found that 18%of all six-hourly TC samples experienced interaction.Extreme interaction cases showed that EFC can reach

    ~120 m s-1d-1during the extratropical-cyclone(EC)stage,an order of magnitude larger than reported in previous studies. Composite analysis showed that positive interactions are characterized by a double-jet flow pattern,rather than the traditional trough pattern,because it is the jets that bring in large EFC from the upper-level environment to the TC center.The role of the outflow jet is also enhanced by relatively low inertial stability,as compared to the inflow jet.Among several environmental factors,it was found that extremely large EFC is usually accompanied by high inertial stability,low SST and strong vertical wind shear(VWS).Thus,the positive effect of EFC is cancelled by their negative effects.Only those samples during the EC stage,whose intensities were less dependent on VWS and the underlying SST,could survive in extremely large EFC environments,or even re-intensify.For classical TCs(not in the EC stage),it was found that environments with a moderate EFC value generally below~25 m s-1d-1are more favorable for a TC’s intensification than those with extremely large EFC.

    Tropical cyclone–trough interaction,upper troposphere,eddy momentum flux convergence

    1.Introduction

    Track and intensity forecasts are two major concerns in tropical cyclone(TC)research.Particular emphasis has been placed on the former of these two aspects,leading to a substantial increase in the skill of TC track prediction. Meanwhile,TC intensity forecasts have received virtually no improvement(e.g.,Emanuel,2000).This lagging behind of intensity forecast skill,which is a widely accepted fact (e.g.,Elsberry et al.,1992;Fitzpatrick,1997;DeMaria et al.,2007;Montgomery and Smith,2014;Smith and Montgomery,2015),means that greater effort should be put into researching TC intensity.

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

    It has long been known that upper-tropospheric environmental flow is an important factor that could modulate the intensity change of a TC.Two reasons are proposed,based on Sawyer–Eliassen balance(SEB)vortex theory(Eliassen, 1952;Sawyer,1956).The first is that,in the upper troposphere,a TC does not retain its axisymmetric structure and thus environmental asymmetric forcing is generally stronger than that of the middle or lower troposphere(e.g.,Pfeffer and Challa,1981).The second is that the inertial stability in the upper layer,representing the resistance of the vortex axisymmetric response to environmental forcings,is much weaker thanthatinthemiddleorlowerlayers(e.g.,HollandandMerrill,1984).Therefore,upper-level asymmetric forcings will excite larger vortex responses,which could easily penetrate into the core region and thus result in TC’s intensity change. This phenomenonis also termedas upper-levelenvironment–TC interaction.

    A large number of interaction cases have been reported in the literature,such as Hurricanes(or Tropical Storms) Dorothy(1966)(Erickson,1967),Elena(1985)(Molinariand Vollaro,1989),Danny(1985)(Molinari et al.,1998),Florence(1988)(Rodgers et al.,1991;Shi et al.,1997),Opal (1995)(Rodgers et al.,1998;Bosart et al.,2000),Bertha (1996)(Hanley,2002)and Gabrielle(2001)(Molinari et al., 2006),overtheAtlanticbasin;Typhoons(orTropicalStorms) Flo(1990)(Wu and Cheng,1999),Gene(1990)(Wu and Cheng,1999),Winnie(1997)(Li et al.,2006),Olga(1999) (Yu and Kwon,2005),Prapiroon(2000)(Yu and Kwon, 2005)and Haima(2004)(Qian et al.,2011),over the western North Pacific basin;as well as Tropical Cyclone Dora(2007) (Lerouxet al.,2013)overthe southwest IndianOcean.These cases also show that intensification induced by upper-level environments can be found in all stages of TCs,such as tropical depression formation(Bracken and Bosart,2000),tropical depression(Bosart and Bartlo,1991)or tropical storm (Shi et al.,1997;Molinari et al.,1998)to hurricane transition,rapid intensification to a category 5 hurricane(Bosart et al.,2000)or supertyphoon(Titley and Elsberry,2000),and tropical storm(Qian et al.,2011)or typhoon(Li et al.,2006) to extratropical cyclone transition.

    Early case studies(e.g.,Erickson,1967)identified interaction visually from weather maps,by observing an upperlevel westerly trough approaching a TC.Thus,such interaction is also termed TC–trough interaction(e.g.,Kimball and Evans,2002;Yu andKwon,2005).Later,a typical case,Hurricane Elena(1985),was studied at length by Molinari and colleagues(Molinari and Vollaro,1989,1990;Molinari et al.,1995).In this case,they adopted the SEB vortex model to diagnose the relationship between eddy fluxes and Elena’s intensity change and found that 200-hPa eddy angular momentumflux convergence(EFC)significantly increasedas an upper-level trough swept over Elena’s northern part.From then on,the EFC value,which can be computed quantitatively and objectively,has been commonlyused as a diagnostic for identifying TC–trough interaction,although the role of EFC in TC intensification was also explored in pioneering works(e.g.,Sundqvist,1970;Challa and Pfeffer,1980; Pfeffer and Challa,1981)that used idealized models.

    Environmental EFC larger than a certain value is usually adopted to identify whether interaction has occurred(e.g., Molinari et al.,2006;Chen et al.,2015).Although many of the above cases underwent intensification when they interacted with upper-level flow using the EFC criterion,there were also weakened or decayed cases(e.g.,Lewis and Jorgensen,1978).Therefore,interaction with upper-level flows does not guarantee TC intensification.In order to verify the effect of upper-level flow on TC intensity,statistical studies are also carried out.Merrill(1988a,b)systematically investigated a large sample of observed outflow-layer wind data near Atlantic hurricanes.The common features of the upper-tropospheric flow pattern for intensifying and nonintensifying hurricanes were summarized,using a rotatedcoordinatecompositetechnique rather than the EFC diagnostic.DeMaria et al.(1993)investigated the relation between intensity change and 200-hPa EFC with 3-yr records of Atlantic hurricanes.They foundthat about 1/3 of the TCs intensified just after the enhancedEFC.The reason for TCs not intensifying was mostly increased vertical wind shear(VWS). Hanley et al.(2001)also examined TC–trough interaction, but using 12-yr records of Atlantic hurricanes and the composite technique.After excluding records that were over cold water or close to land,they separated the interaction events into four composites.Their results showed that 78%of superposition and 61%of distant interaction cases deepened. Although these results confirm that interaction favors TC intensification in a statistical sense,to identify whether or not a particular trough interaction is conducive to TC intensification is still not straightforward(Hanley et al.,2001;Leroux et al.,2013).Besides,these two studies ignored the role of upper-levelinertial stability,emphasizedin a numberof other studies(e.g.,Holland and Merrill,1984;Rappin et al.,2011).

    Since the comprehensive work of Hoskins et al.(1985), the perspectiveof isentropic potential vorticity(PV)has been introduced for viewing dynamical processes in the atmosphere.The PV perspective is not only applicable to TC motion(e.g.,Wu and Wang,2000;Chan et al.,2002;Choi et al.,2013),but also to TC intensity change(e.g.,Molinari et al.,1998).Molinari et al.(1998,1995)also applied the PV perspectiveto interactioncases,in which interactioncouldbe viewed as an upper-level synoptic PV anomaly(the environment)interacting with a lower-level mesoscale PV anomaly (the TC).From such a perspective,although interaction is clearly and effectively shown as a superposition of two PV anomalies,thereisnoquantitativewayofidentifyingwhether interaction has occurred.Therefore,the EFC diagnostic is more frequently used than the PV perspective.

    To date,statistical studies in this field have been relatively less common—especially for the western North Pacific,where TCs are most active compared with other basins. Therefore,the climatology of TC–upper-level-flow interactions and its possible impact on TC intensity over this basin remain unclear.Before addressing the“good trough/bad trough”issue,we should first carry out a preliminary but thorough study to reveal more basic features of upperlevelenvironment–TCinteractions,especiallyforthe western North Pacific basin.For example:What are the characteristics of upper-level flow patterns during interactions,besides the well-known trough signature?What is the relationship betweenthe upper-levelenvironmentalfactorandotherlargescale factors(e.g.,SST and VWS)in controllingTC intensity over this basin?What is the nature of inertial stability during interactions?The present paper aims to address these aspects by examining25-yrTC records and revealing statistical characteristics of upper-troposphericenvironment–TCinteraction over the western North Pacific.

    Section 2 describes the data and method.Section 3 presents the climatology of interactions.Section 4 describes the characteristics of upper-level flow patterns.The relationship between environmental factors and TC intensity change follows in section 5,and conclusions are given in section 6.

    2.Data and method

    2.1.Data description

    The TC“best track”data over the western North Pacific were obtained from the Regional Specialized Meteorological Centre(RSMC),Tokyo.Given that aircraft reconnaissance of TCs in this basin terminated in 1987,making observations of TCs thereafter completely dependent on satellite retrieval(Martin and Gray,1993),for consistency,only data after 1987 were used.Specifically,25 years of data,from 1987 to 2011,including six-hourly TC positions,near-center maximum surface wind speed,and minimum sea level pressure,were chosen for the present study.

    ERA-Interim(Dee et al.,2011)data were employed to describe the TC environment.The benefits of using these reanalysisdataforstudyingtheupper-levelenvironmentofTCs have been demonstratedin a numberof studies(Molinari and Vollaro,1990;Molinari et al.,1992).In the present study, six-hourly pressure-level wind fields,as well as SST data,on a 1.5?×1.5?grid,were used.

    2.2.Method

    In order to identify upper-level environment–TC interactions,EFC was computed in storm-relative cylindrical coordinates.The cylindrical coordinates contained 36 grids in the azimuthal direction,with 10?intervals and 28 grids in the radial direction with 0.3?intervals(about 33.3 km).Latitude/longitude gridded data were interpolated to the cylindrical gridded data using a 16-neighbouring-pointbicubic polynomial.EFC was calculated following Molinari and Vollaro (1990):

    where v and u are the tangential and radial componentsof the storm-relative wind vector(i.e.,after subtracting the translation velocity of a TC from the full wind vector),r is the distance to the TC center,the overbar is the azimuthal mean, andthe primeis the deviationfromthe mean.Consistent with previous studies(e.g.,Hanley et al.,2001),the units of measurement for EFC is m s-1d-1.Molinari and Vollaro(1989) showed that the errors of EFC calculated in the inner radii (within 300 km)may exceed 40%.EFC evaluated at large radii is more reliable;however,it may not have an immediate impact on TC intensity.Therefore,the present study used a radial averagefrom 300 to 600 km to identifyTC–upper-flow interactions,following Hanley et al.(2001).

    According to the SEB theory,within a slowly evolving TC,the TC’s axisymmetric secondary circulation is largely balanced by several forcings,including EFC.This implies that if other forcings are not important,larger EFC will result in a stronger response of the TC’s secondary circulation. If the strength of the TC’s secondary circulation is positively correlated with its intensity,which is generally true,then one may expect a simple rule that larger EFC is more likely to enhance the outflow and thus intensify the TC(e.g.,Challa and Pfeffer,1980;PfefferandChalla,1981;MolinariandVollaro, 1990;Qian et al.,2011).However,when EFC becomes large, other factors may come into play.

    One of the important factors concerning a TC’s interaction with upper-level environmentalflow,besides EFC,is the inertial stability,usuallydefinedin the SEB theory(e.g.,Rappin et al.,2011)as

    in which the azimuthal-averaged absolute vorticityζais defined as

    whereζis relative vorticity and f is planetary vorticity.Notice that the factor of(f+2v/r)also appears in the momentum forcing functions[e.g.,Molinari and Vollaro,1990;see the first term on the rhs of their Eq.(2)].When simultaneously taking into account the forcing effect of momentum source and the resistance effect of inertial stability,the factor of(f+2v/r)can be eliminated.Therefore,the axisymmetric outflow is proportional to a forcing function in terms of EFC[rather than EFC multiplied by(f+2v/r)]and inversely proportional to the inertial stability,proxied byζa(rather than I).In the SEB theory,it is generally required that the azimuthal-averaged absolute vorticityζa>0(inertial stable).Smallerζaindicates less resistance to the environmental forcing(e.g.,EFC)and a larger radial extent of the vortex response(e.g.,Holland and Merrill,1984). In the present study,the absolute vorticityζabefore azimuthal averaging was used to represent localized inertial stability:

    whereλis the azimuth.

    Besides EFC and inertial stability,two other environmental factors,SST and VWS,were also computed in cylindrical coordinates for each six-hourly TC record.According to Molinari and Vollaro(1993),VWS is defined as the difference between the 200-and 850-hPa wind vector averaged within 500 km from the TC’s center:

    In this definition,the angled brackets indicate area-weighted averaging:

    whereβj=0.3j is the radial angle,in degrees(interval of~33.3 km),and j is the radial grid index,starting from J0to J.J0=0 and J=15 are used for averaging within 500 km. Similarly,SST is also averaged within 500 km.

    3.Climatology of upper-tropospheric TC–

    environment interaction

    3.1.Characteristics of the dataset

    The“best track”dataset from the RSMC used in the resent study contains 628 TCs.Records at 0000,0600,1200 and 1800 UTC were chosen,yielding 21 685 six-hourlysamples,including tropical depressions(TD),tropical storms (TS),typhoons(TY),and extratropical cyclones(EC).Since theinteractioncouldoccurat anystageof a TC(Hanleyet al., 2001),keeping all these records yields more samples for reliable statistics.TC intensity change is defined as the forward difference of sea-level minimum central pressure(SLP)i.e.,?p(t)=p(t+?t)-p(t)where?t=6 h.Negative(positive) pressure change means intensification(weakening),and zero means no change.

    Firstly,the 200-hPa EFC averaged within the 300–600 km radial band for all 21 685 samples was calculated.Figure 1 shows the distribution of the results.Following previous studies(e.g.,DeMaria et al.,1993;Hanley et al.,2001),an EFC value of 10 m s-1d-1is usually defined for identifying interactions.In the present study,interaction is defined similarly,but an EFC value above 10 m s-1d-1means positive interaction,while EFC below-10 m s-1d-1means negative interaction.According to this definition,about 17.7%of the samples experienceinteraction:3.7%negativeand14.0% positive.The proportion of positive interaction(14.0%)is smaller than that(~23%)over the Atlantic basin,as shown by Hanley et al.(2001).More positive than negative samples over the western North Pacific basin,which is quite similar to the situation over the Atlantic basin,indicates that a positive EFC environment is more favorable for the maintenance of TCs than a negative one.However,a larger EFC value does not guarantee a greater number(or larger proportion)of intensifying samples,as can be seen from Fig.1.

    Figure 2 shows the percentage of different TC types,including TDs,TSs,TYs and ECs,as identified by the RSMC. The TS and TY samples account for 37%and 25%of the total,respectively,resulting in 62%of samples being canonical tropical cyclones(i.e.,TSs and TYs).Besides,26%of samples are TDs and 12%are ECs.Each sample type has a proportion that interacts with upper-level flows.This result verifies that interaction can occur at any stage of development.We can also see that,although the proportion of EC samples is least,almost 60%of them experience interaction with upper-level flows.This ratio is significantly larger than for the other three types,showing that EC samples are more likely to interact with upper-level flow,while it is relatively rare to see interaction in the other three types.

    3.2.Spatial distribution

    The spatial distribution of interaction samples over the western North Pacific is shown in Fig.3a.The southwest to northeast orientationof the distribution from the South China Sea to the Bering Sea suggests that interactions generally occur after a TC’s recurvature to the northeast.The maximum frequencyof occurrenceof 0.5 per year per grid can be found aroundJapan.The interaction active regions are within(25?–45?N,125?–165?E),with two peaks at about 30?N and 40?N (Fig.3b),and 138?E and 151?E(Fig.3c),respectively.These results indicate that most interactions occur when TCs move into the midlatitudes,right under the influence of upper-level westerlies,but before dissipating.Notice that there are still many interaction samples(434)south of 20?N,over both the north South China Sea and east of the Philippines.More than 60%(276)are during the genesis or development stages (before reaching their peak intensities)when interaction occurs.

    4.Characteristics of upper-level flow patterns

    during interaction

    Previous studies(e.g.,Hanley et al.,2001)have already identified the most common upper-level flow,i.e.,a westerly trough pattern during interaction.Therefore,the interaction is also named TC–trough interaction.However,EFC is an abstract index and EFC larger than a certain threshold does not necessarily mean the flow would be characterized by a westerly trough.There could be other types of flow patterns with different synoptic signatures,although they may be related more or less to a“V”-shaped trough in the vicinity of a TC.Summarizing these features of the flow pattern will help to identify interactions using upper-level wind fields only.

    Fig.1.Distribution of 200-hPa EFC averaged over a 300–600 km radius.There were 21 685 samples in total from 1987 to 2011.

    Fig.2.Distributions of different TC types.There were 21 685 samples in total and 3890 interaction samples from 1987 to 2011.

    4.1.Extreme interaction

    Firstly,extreme interactions(extremely large magnitudes of EFC)from 1987 to 2011 were identified.Table 1 details the basic information of three negative and three positive extreme interaction samples.The three negative samples have EFC values of about-70 m s-1d-1.The three positive interactions show EFC values larger than 120 m s-1d-1,almosttwicethatofthenegativeinteractionsinmagnitude.The identified extreme positive EFC values are much larger than those reported previously in the literature[e.g.,~40 m s-1d-1,as shown by Wu and Cheng(1999)].Except for TC Andy(1989),which was in the TS stage,the other five samples were all in the EC stage.Low SST and large VWS are their common environmental features.It is noticeable that the negative extreme interaction samples could undergo filling(?p>0),while the positive ones could intensify further. Since these interaction samples were mainly in the EC stage, theirintensity(i.e.,minimumSLP)maynothavedependedas stronglyonSST and VWS as it wouldin the TS or TY stages, for example.Extremely large positive EFC from upper-level flows would possibly have been responsible for their intensification.

    Fig.3.(a)Spatial(interval:0.05 per square grid per year;9-point smoothing operator applied),(b)latitudinal and(c) longitudinal distributions of the interaction samples from 1987 to 2011.The plots are on a 1.5?×0.5?grid.

    Table 1.TC samples of extreme interactions.Environmental SST and VWS are averaged within 500 km and EFC is averaged within the 300–600 km radial band.

    Fig.4.200-hPa wind fields(m s-1;shaded for wind speed>50 m s-1)for(a–c)three TC samples experiencing extreme negative EFC values and(d–f)three TC samples experiencing extreme positive EFC values.Black dots show TC locations;thick solid lines indicate upper-layer troughs.

    Figure 4 shows the 200-hPa flows for three extreme negative interaction samples(Figs.4a–c)and three positive samples(Figs.4d–f).All six samples show that extreme EFC is associated with upper-level westerly troughs.For extreme negative interaction,TCs are right under the influence of northwesterly wind upstream of the trough axes and tend todissipate.In contrast,the three positive interactions cases (Figs.4d–f)were under the influence of southwesterly wind ahead of the troughs axis and retained a relatively high level of intensity compared to the negative cases(Table 1).Experiencingextremely large EFC,these TC samples were mostly close to westerly jet streams(shaded areas in Fig.4),or even inside the jet core regions.Although strong wind speed associated with the jets exerted large asymmetric forcing that resulted in large EFC,it also brought strong VWS to the TCs (Table 1).In such a strongly sheared environment,TCs cannot maintain their typical structure and usually undergo transition into ECs.

    4.2.All positive interaction

    According to the SEB theory,positive(negative)interaction would transport cyclonic(anticyclonic)angular momentum from the upper-level environment to TCs and enhance (reduce)their outflows.Due to its possible effects on TC intensification,previous studies have tended to concentrate on positive interaction.Following this line,we focused on upper-level flow patterns of all positive interaction samples. To extract the characteristics of the flow pattern from a large number of samples,the composite technique was used.Although the composite may smooth out some important characteristicsofindividualTCs,it emphasizesthesignaturesthat repeatedly appear.

    There were several considerations regarding the composite procedure.Firstly,as EFC is computed using stormrelative wind(i.e.,the translation speed of the TC is subtracted),it would be more meaningful to examine the stormrelative wind pattern.Secondly,to reduce the smoothing effect of the composite technique and to obtain more flow patterns,thesamplesweredividedintoeightgroupsaccordingto the storm-relativewind directions averaged within 500 km of TCs.It is worth noting that different synoptic patterns or different positions of troughs will be somehow smoothed in the composite.However,there is no best way to groupindividual samples into different flow patterns,and different composite criteriawouldalsosufferfromthesameproblem.Thepresent choice was shown to be effective in the sense of keeping flow features in the regions near the TC center and minimizing the blurring effect not too far from TCs.Besides,the most dominant flow type,the trough pattern,should be most frequently identified,in accordance with previous studies.Thirdly,different TCs may locate at different latitudes.Thus,traditional compositeanalysis using the latitude–longitudegriddedwind field at different latitudes may be compromised by the problem of the meridian converging at high latitudes,and thus blurring some flow features.Taking this into account,all the composites were performed in cylindrical coordinates on a sphere with respect to their origins(i.e.,the TCs’centers).

    Table 2 lists the statistics of the eight composites and Fig. 5 shows the correspondingcomposite flow patterns as well as localized EFC.Local EFC is defined similar to Eq.(1),except that no azimuthal average is taken so that the local EFC is also a function of azimuth.Not surprisingly,the southwest composite contains 1359 samples,which is almost half of all positive interaction samples(Table 2),indicating that classical positive interactions usually occur right under the southwesterly wind between a westerly trough and a ridge(Fig. 5h).Two westerly jets(black bold arrows)locate southwest and northeast of the TC’s center,respectively,and the TC’s axisymmetric outflow tends to enhance the northeast outflow jet but weaken the southwest inflow one.The southwest jet imports cyclonic eddy angular momentum to the TC,while the northeast jet exports anticyclonic eddy angular momentum from the TC.Both jets result in positive EFC channels (redareas),so that the300–600kmradialbandhasthe largest mean EFC value(26.9 m s-1d-1)among the eight composites.The west composite(Fig.5a),containing 983 samples,shows a strong westerly jet crossing the TC’s center.A smaller amplitudeof the synopticdisturbancealso indicates a smaller curvatureof the flow and larger wind speeds,as compared to the southwest composite.The south composite has 323 samples,which is the third highest amongthe eight composites.The correspondingflow pattern(Fig.5g)is similar to that of the southwest composite,except that the trough to the west of the TC center intensifies with its bottom extending to the south of the TC.The wind near the TC center is a uniform southerly,but weaker comparedto that of the southwest composite,leading to a relatively weak mean VWS(15.1 m s-1), butstill retainingalargevalueofmeanEFC(24.8ms-1d-1). The whole pattern resembles an“S”.This“S”-shaped flow has been reported to favor the intensification process of TCsduring the extratropical transition stage(Qian et al.,2011). However,from the statistical point of view,only 18.8%of the samples intensifiedand morethan43%of the samples decayed.The northwest composite(Fig.5b),which is the least among those composites(southwest,west,south and northwest)related to westerlies,contains only 204 samples.The flow featuresare quite similar,exceptthat the inflow jet shifts to the northwest and the outflow jet shifts to the southeast.

    Table 2.Characteristics of the eight composites for all positive interaction samples.Composites were selected according to the averaged 200-hPa wind direction within 500 km of TCs.Note that the mean SSTs were computed using samples that were over the ocean.Mean SST and VWS are averaged within 500 km and the mean EFC and inertial stabilityζaare averaged over the 300–600 km radial bands.

    Fig.5.Eight composites of 200-hPa storm-relative flow(arrows;m s-1)and local EFC(shaded,m s-1d-1) for all positive interaction samples.Black bold arrows indicate wind speeds>15 m s-1.Samples are grouped according to the averaged wind direction(shown in the lower-left corner of each panel)within 500 km of TCs. The composites were produced in cylindrical coordinates on a sphere,but the results were interpolated onto 1?×1?grids for a better visual.Black dots indicate TC centers and two black concentric circles in each panel show 300 km and 600 km radii from TC centers.The background map is shown for reference only.

    These four types,generally related to the upper-level westerlywind,containmoresamplesthantheremainingfour. A more interesting signature other than the trough in these composites is that there are two jets in the vicinity of the TC, like a dipole,resulting in an obvious azimuthal wavenumber-2 structure of eddy flux convergence.Their roles can be clearly identified by the local EFC in which no azimuthal average was calculated in Fig.5:the two jets bring in two local EFC maxima(red areas)with respect to the TC center, regardless of jet directions(i.e.,inflow or outflow).However, inside the trough regions(northwest quadrants between two jets),only negative local EFC(blue area)is found,and thus the troughs act to reduce the EFC when averaging within the 300–600kmradialband.Therefore,thisdouble-jetsignature, with the outflow anticyclonic curved and the inflow cyclonic curved,results in a largeEFC value rather than the traditional trough signature.These two jets cannot be ascribed to the interaction between the TC and trough because there would be such a flow pattern even without a TC.That is why their influence is also viewed as external(e.g.,Molinari and Vollaro, 1989).

    The remaining four composites of north,northeast,east and southeast(Figs.5c–f)have relatively fewer samples. Such few samples may not produce a statistically significant percentage of intensifying or weakening samples,especially for the east composite.It seems that these interactions occur at lower latitudes where the prevailing upper-level wind has an easterly component(Figs.5c–f)rather than synoptic waves embedded in westerlies.These TC samples are probably in their early stage as TDs because the underlying SST is relatively high(above 27?C,except that of the southeast composite).The previously mentioned double-jet signature also applies to these composites,exceptthat the jets are much weaker,but the traditional trough pattern can barely be identified in these four composites.

    Figure 6 shows the local inertial stability,i.e.,ζasuperposed by the composite wind fields.In terms of local inertial stability,the role of each jet,located in different quadrants, could be different because relatively weak inertial stability over a certain quadrant of the TC could minimize the energy expenditureandfacilitatetheformationofanoutflowchannel in that quadrant(Rappin et al.,2011).From Fig.6 we can see that all compositeζavalues are positive,although this may not be true for each individual sample.Thus,the composite wind flows are generally inertially stable.However, the azimuthal distribution ofζais quite uneven,showing obvious azimuthal wavenumber-1 asymmetry.In general,the jet cores separate theζainto two parts—a larger part that is more inertially stable and locates inside the troughs,and a smaller part that is less stable.As can be seen from Fig.6, therolesofthe twojets aredifferentintermsofinertialstability because the outflow jet is roughlywithin the less inertially stable region while the inflow jet is inside the more stable region.Given roughly the same environmental EFC forcing, the contribution of the outflow jet to the TC’s axisymmetric secondary circulation is more significant than that of the inflow jet due to weaker inertial stability.Besides,as mentioned previously,smaller axisymmetricζa,after azimuthal averaging,indicates less resistance to the environmentalEFC forcing and a stronger vortex response in the TC’s radial–verticalcirculation.Theseconditionsapplytothefourgroups of north,northeast,east,and southeast composites(Figs.6c–f),as theseinteractionsoccurat relativelylowlatitudes(Table 2),characterized by smaller planetary vorticity.

    It is worth mentioning that there are some similarities between the eight flow patterns identified here and those proposed in previous studies.For example,Holland and Merrill (1984)showed the manner in which a cyclone interacts with the subtropical westerlies and produces a poleward outflow channel(see their Fig.15,and notice that their illustration is in the Southern Hemisphere).They also emphasized the effect of the subtropical jet on enhancing the TC’s outflow, which is quite similar to Figs.5f and g shown here.Figures 5f and g are similar to Holland and Merrill’s(1984)illustration in the sense that the subtropical jet is located just north of the TC center and the TC keeps a certain distance away from the jet(heavy black arrows).However,the maximum radius is~1000 km in the composites.Therefore,the complete nature of the jet in Fig.5f cannot be shown.Notice that there are also jet maxima just to the north of the TC center in Fig.5h,but TCs are already inside the westerly jet core with large VWS,which is different from the picture of Holland and Merrill(1984).Another example was given by McTaggart-Cowan et al.(2006),who presented a conceptual model of tropical transition in a dipole-blockingenvironment (see the right-hand column of their Fig.2).Their flow pattern is quite similar to Figs.5f and g shown here,in which the anticyclone to the northeast of the TC center develops intensively so that it is gradually cut off from the westerlies and forms a cutoff blocking.Such a blocking flow pattern,as noted by McTaggart-Cowan et al.(2006),would facilitate a midlatitude vortex precursor transiting into a TC(known as tropical transition),whereas the similar flow pattern in Fig. 5g is likely to facilitate the transition of a TC into an EC (known as extratropical transition),depending on the accompanying VWS and underlying SST.A third example is the typical tropical upper-tropospheric trough(TUTT)flow pattern given by Fitzpatrick et al.(1995).The TUTT pattern (see their Fig.2),in which the axes of the shallow troughs lie in the northeast quadrants and also tilt northeast,resembles those of the north and northeast composites shown in Figs. 5c and d here.

    Fig.6.As in Fig.5 but for local absolute vorticityζa(shaded;10-5s-2),defined in Eq.(2),as a proxy for local inertial stability.

    The EFC criterion,as shown here,actually does not guarantee the occurrence of westerly troughs but reflects azimuthalshearsoftheflow crossingtheTC’s center,especially those induced by wavenumber-1 asymmetry.Curvatures ofthe flow(i.e.,cyclonic curved inflow and anticyclonic curved outflow)may eventually decide the sign of EFC.It was also found that extremely large EFC only occurs in those environments with large wind magnitude(e.g.,Fig.4),usually inside westerly jets.Therefore,we emphasize the feature of a double-jet signature of upper-level flows when interaction occurs because it is the curved jets,rather than troughs,that import large EFC from the environment to the TC(Fig.5), resulting in significant upper-level asymmetric forcing on the TC’s axisymmetric secondary circulation.

    5.Environmental factors and TC intensity change

    5.1.Relationship between EFC and other environmental factors

    Studies(e.g.,Wu and Cheng,1999;Chen et al.,2015) have shown that SST,VWS and EFC are three large-scale environmental factors affecting TC intensity.Their effects may overlap to result in a more favorable environment for TC strengthening,or be cancelled to lead to a less favorable one.Therefore,the relationship between EFC and other environmental factors was investigated,based on the 25-year dataset.

    According to the SEB theory(section 2.2),large EFC would favor TC intensification.However,from Table 2 we canseethatthereisnosimple“l(fā)argeEFC–largepercentageof intensifying samples(or small percentage of decaying samples)”relationship,primarily due to two reasons.One reason is that larger EFC is generally accompanied by greater inertial stability(see the last two columns in Table 2).Inertial stability is a measure of the resistance of vortex responses to eddy forcing.High inertial stability means that strong forcing will result in limited responses.As a TC intensifies,especially when doing so rapidly,it requires low inertial stability (e.g.,Rappin et al.,2011);the higher inertial stability shown here is partially responsible for reducing the positive effect of large EFC.Another reason is that TC intensity change is also controlled by two other environmental factors,i.e.,SST and VWS.Larger EFC corresponds to lower SST and larger VWS so that EFC’s positive effect on TC intensification is generally cancelled out by the negative effect of low SST and large VWS.All these facts can be explained by the variation in the latitude at which the TC is located,as a higher latitude usually means largerζa(larger planetary vorticity f),lower SST,larger wind speed,and thus larger VWS and EFC.

    Figure 7 shows scatter plots of EFC versus SST for TC samples grouped by different TC stages.For the TD,TS and TY types(Figs.7a–c),a large proportion of samples occur in a favorable environment of a warm ocean(SST>25?C). There are only a few samples whose underlying SSTs are lower than 20?C.The mean EFCs are slightly larger than 0,showing a skewness towards a positive EFC environment. For the EC type samples,as they have actually moved over cold water at higher latitudes,SST is much lower than those of TDs,with a mean of~16?C and large variance.Their mean EFC is close to 20 m s-1d-1,already exceeding the threshold of positive interaction,which means that ECs are essentially different from other TC types and large upperlevel EFC may be one of their characteristics.For all four types,it was found that the magnitude of EFC does not exceed 40 m s-1d-1if SST is higher than 25?C.However, extremely large EFCs occur frequently in lower SST environments.For the TS and TY samples,the highest EFC,at 80–120 m s-1d-1,was in the SST range of 15?C–20?C.For the EC samples,this was even more prominent,in that EFC couldreachabove120m s-1d-1whenSST was below15?C.

    The relationship between VWS and EFC(Fig.8)is the reverse of that between SST and EFC.Extremely large EFC usually occurs in a large VWS environment,but large VWS does not necessarily induce large EFC,as EFC also requires some curvature of the wind flow.According to the SEB theory,large EFC is usually expected to significantly enhance a TC’s axisymmetric outflow,strengthen the updraft near the TC center,and lead to TC intensification.Here,it is shown that extremely large EFC(>40 m s-1d-1)is usually accompanied by relatively low SST and strong VWS,so that the positive contribution of EFC to a TC’s intensification is cancelled out by the negative contribution of the low SST and strong VWS.Therefore,EFC larger than 40 m s-1d-1does not guarantee TC intensification.On the contrary,such an environment would destroy a TC’s warm-core structure and possibly lead to its transformation into an EC.

    5.2.Environmental factors and TC intensity change

    In order to investigate the relationship between environmental factors and TC intensity change,intensifying and weakening TC samples were selected.Generally,whether a TC is intensifying or weakening can be identified from the sign of pressure changes.Taking into account that the pressure recordsin the RSMC best track dataset are discrete,with aminimumintervalof2hPa,andthattherecordsmaycontain some uncertainties,records with|?p|>3 hPa were considered as the threshold for intensity change.

    Scatterplots oftheintensifying(?p<-3hPa)andweakening(?p>3 hPa)TC samples and their environmental factors are shown in Fig.9.For the intensifying samples(righthand column in Fig.9),there are obvious scatter clouds within regionsof SST>25?C(the circlein Fig.9b)andVWS<15 m s-1(the circle in Fig.9d).These clouds consist of samples from classic TC types,i.e.,TD,TS and TY,while the EC samples are generally located outside the clouds(see Figs.7d and 8d).For the classic TC types,environmentswith SST>25?C and VWS<15 m s-1are generally known as necessary conditions for their maintenance or development. Their EFC magnitudes seldom exceed 25 m s-1d-1.For the EC-type samples outside the dense clouds,it is surprising that all intensifying samples are in a positive EFC environment.Other factors,such as VWS or SST,place no obvious constraint on the distribution of these samples.This can be explainedbythe factthatECs arebaroclinicsystemsandgain energyprimarily from baroclinic instability,while classic TC types(TD,TS and TY)gain energy mainly from latent(dia-batic)heating.For weakening samples(left-hand column in Fig.9),the scatter clouds are not as concentrated as their intensifying counterparts,and there are no obvious features in selecting an environment.

    Fig.7.Scatter diagrams of SST(?C)averaged within 500 km versus 200-hPa EFC(m s-1d-1)averaged within 300–600 km for(a)TD,(b)TS,(c)TY and(d)EC samples.Samples of landing or near land(within 200 km) are excluded since no reliable SST can be obtained,yielding 15 999 samples.Black crosses indicate the mean EFC and SST for each panel.

    Since ECs are different from other TC types in response to their environments,we should exclude them for further study.It is also known that large EFC may not correlate with simultaneous TC intensity,but has a lead–lag impact on later intensity change(e.g.,Molinari and Vollaro,1990).Therefore,it is more appropriate to define an interaction event as three consecutive six-hourlysamples experiencingat least 10 m s-1d-1EFC.This definitiongroupsthose interactionsamples into events and excludes samples that only once or twice exceed the 10 m s-1d-1threshold,which has little persistenteffectonTC intensity.Sampleswith maximumsustained wind speed of<17.2 m s-1were not considered,so as to exclude most EC samples,since ECs are different from other types.This procedure also excluded TDs,as there are large uncertainties in intensity records of weak intensity samples (e.g.,Barcikowska et al.,2012).

    There were only 166 interaction events identified using the above constraints.The corresponding statistics of these events are listed in Table 3.The number of weakening events was twice as manyas intensifyingevents.Therefore,drawing the simple conclusionthat EFC>10m s-1d-1may facilitate TC intensification is still not possible.Other factors,such as SSTandVWS,seemmorefavorableforintensification.After excluding EC samples,the mean EFCs for both intensifying and weakening events remains at~25 m s-1d-1.This value can be viewed as a threshold,and if environmental EFC exceeds this value the positive effect of EFC would be offset by accompanying large VWS,as well as low SST,possibly leading to extratropical transition of TCs.

    Fig.8.As in Fig.7 but for VWS(m s-1).

    Table 3.Statistics of interaction events.Mean?p was estimated by the linear trend of SLP during the event using the least-squares fit method.

    It is known that TCs usually weaken when close to land or are already landing.Thus,we added another condition that interaction events near land,or that had already landed, shouldbeexcluded.Itis alsoknownthat,givenanunderlying SST,TCs have their maximum potential intensity(MPI).If a TC is closetoits MPI(usuallymeasuredassurfacemaximum wind speed),then any atmospheric environment is hostile to a TC’s intensity(Merrill,1988b).That is,there is no room for the TC to develop.We imposed the final constraint that TCs should be at least 30 m s-1away from their MPI,which was estimated using the empirical relation given by Zeng et al.(2007)for the western North Pacific as:

    where A=15.69 m s-1,B=98.03 m s-1,C=0.1806?C-1, Tsis SST,and T0=30.0?C.

    Fig.9.Scatter diagrams of(a,b)SST(?C)and(c,d)VWS(m s-1)averaged within 500 km versus 200-hPa EFC(m s-1d-1)averaged within 300–600 km for(a,c)weakening samples(?p>3 hPa)and(b,d)intensifying samples(?p<-3 hPa).TC samples landing or near land(within 200 km)are excluded since no reliable SST can be estimated.Black crosses indicate the mean VWS,SST and EFC for each panel.Black circles highlight the scatter clouds in panels(b,d).

    Table 3 shows that when the“no-landing”condition is added,the number of weakening events is greatly reduced from 111 to 67,and the rate of SLP change is also reduced from 5 hPa to 4.5 hPa per six hours.After the“havepotential”condition is applied,the number of weakening events drops to 26,roughlyhalf of that of intensifyingevents. Besides,the intensifying rate increases from 1.4 hPa to 1.7 hPa per six hours.However,even with these constraints,it is barely possible to draw the conclusion that higher EFC values will result in more intensifying events,because the mean EFC for an intensifying event is 16.2 m s-1d-1,which is relatively smaller than that of weakening events.Thus,regardless of these constraints,the intensifying events seem to be affected more by the SST and VWS rather than by EFC environments.As noted earlier,EFC exceeding the 25 m s-1d-1threshold is likely to induce very large VWS and low SST,and thus negative effects of the environment are more pronounced than positive ones,leading to weakening TCs or even deforming them into ECs.Here,EFC larger than this value tends to slow down the intensification rate or increase the possibility of filling a TC.So,we regard a moderate EFC environment,larger than 10 m s-1d-1but generally below 25 m s-1d-1,as a relatively good environment for TC intensification.

    6.Conclusions

    A preliminary statistical survey of the interactions between TCs and upper-level environmental flows over the western North Pacific is presented in this paper using a 25-year(1987–2011)best track dataset and ERA-Interim data. Interactions are defined as the magnitude of EFC exceeding 10 m s-1d-1.According to this definition,it was found that 17.7%of 21 685 six-hourly TC samples experienced interaction:14%positive and 3.7%negative,accordingto their EFC signs.These interactions were frequently identified around Japan,where TCs recurve into midlatitude westerlies.

    Six extreme interaction cases showed that negative(positive)interactions are characterized by upper-level northwest (southwest)cross-center wind associated with the upstream (downstream)flow of a trough.Extreme negative EFC values reach-70 m s-1d-1,whose magnitude is only half that of positive ones(120 m s-1d-1).These negative interactions,arising from slightly modifying the interaction definition,as compared to previous studies,also provides new features when EFC is negative.This kind of interaction also relates to upper-level troughs(Figs.4a–c).The main difference is that the TC locates upstream of the trough,rather than downstream as in traditionally defined(i.e.,positive)interaction.

    Composite analysis showed that positive interactions are characterized by an apparently curved upper-level crossing-center flow that consists of two jets.These two jets,identified in the present study through analyzing localized EFC, areamoreimportantsignaturethanthetraditionallyproposed trough pattern because large positive local EFC comes from these two jets rather than from the trough itself.Therefore, the interaction strength(e.g.,EFC value)largely depends on the jet strength as well as their curvature.It is also shown that the outflow jet,exporting anticyclonic eddy angular momentum,is more important than the inflow one,importing cyclonic momentum,as the axisymmetric outflow enhances the outflowjet strengthandlowers the adjacentinertial stability.That is why it is commonto observean enhancedoutflow channel in the northern part of a TC without seeing any clear inflow channel(e.g.,Rappin et al.,2011,Fig.17).

    Generally,according to the SEB theory,large positive EFC intheuppertropospheremeansthat largecyclonicangular momentumis transportedtowards the TC fromits vicinity and is thus favorable for TC axisymmetric outflow enhancement or even intensification.Previous studies(e.g.,DeMaria et al.,1993)have already pointed out that such a positive effect will be offset by accompanying strong VWS and low SST.Here,we suggested separating the interactions into extreme cases(EFC larger than 25 m s-1d-1)and moderate cases(10–25 m s-1d-1),according to Figs.7 and 8.It was found that extreme interactions generally occur at high latitudes,meaning the environmental SST is very low and the VWS is quite strong.Besides,the inertial stability is also largeathighlatitudes(Table2),whichmeansthatthepositive effect of extremely large EFC will be scaled down.During such cases,the TC tends to transform into an EC.Only those ECs after transformation,which are less dependenton underlying high SST and low VWS,can survive or intensify in extremely large EFC environments(Fig.9b).For moderate interactions that occur in the TS and TY stages,it is shown that their responses to EFC environments are quite different from that of the EC type.These traditional TCs can only survive in environments with small-to-moderate VWS,to maintain their warm core structures.The intensities of these TCs are also more sensitive to the underlying SST than EFC.Therefore,the“goodtrough/badtrough”issue still remains unclear for these traditional TCs.

    The EFC threshold of 25 m s-1d-1,empirically identified here,also distinguishes the interaction case of Hurricane Elena(1985)reported by Molinari and Vollaro(1989, 1990)as a canonical one,because this case only just reached the proposed upper limit of favorable interaction and then experienced a rapid intensification.It is also interesting to find cases over the western North Pacific that are similar to this canonical one,which may help us identify these“good”troughs.We selected TC cases with EFC larger than 20 m s-1d-1and VWS smaller than 10 m s-1for at least one sixhour snapshot.There were only 153 samples(5%of positive interactions)that met these requirements,including several consecutive records in a single TC.Six representative interacting TCs—two during the formation stage(cases 1 and 2), two during the mature stage(cases 3 and 4),and two during the dissipation stage(cases 5 and 6)—were selected(shown in Fig.10).As we can see,favorable interaction(cases 1–4)does not indicate a westerly trough in the vicinity of the TC(e.g.,cases 2 and 3).Even in the subtropics,flows with certain curvature can cause interaction.The enhanced EFC (green shaded areas)generally lead to subsequent intensifying processes.This is similar to the case of Hurricane Elena (1985),in which an EFC spike acted as a triggering effect for deepening the TC.Obviously,unfavorable interaction(cases 5 and 6)can be attributed to the increasing VWS and decreasing SST;although,at the interaction times(vertical gray lines),VWS does not exceed 10 m s-1.It also seems coincident that the 200-hPa flows in cases 5 and 6 resemble the non-intensifying pattern summarized by Merrill(1988b,Fig. 11),in which the anticyclones east of the TCs have closed streamlines.

    Based on the conclusions drawn from the present study, as well as these representative cases,we can gain some insights into the“good”trough issue.Favorable interactions for ECs should be characterized by very large environmental EFC(at least above 25 m s-1d-1),and the“l(fā)arger is better”conclusion may apply for ECs regardless of the associated VWS and underlying SST.However,for traditional TC types,such as TDs,TSs or TYs,favorable interaction occurs with only moderate EFC above 10 m s-1d-1,but generally below 25 m s-1d-1,in favor of low VWS.This requires that TCs should keep a certain distance from the disturbed westerly jet(usually with a trough embedded upstream)and not penetrate into the jet core regions(such as in extremeinteraction cases),as was schematically illustrated by Holland and Merrill(1984,Fig.15)and also exemplified by Rappin et al. (2011,Fig.17),so that their warm core structure will not be destroyed by persistently strong VWS.Besides,interaction should not last more than around two days,playing its role as a triggering effect for initiating wind-induced surface heat exchange(Emanuel,1986)or the eyewall replacement cycle (Willoughbyet al.,1982;Lerouxet al.,2013),as long-lasting interactionwill eventuallybringpersistentVWS,whichis destructive to a TC’s warm core.Finally,interaction should appear at low latitudes so that environmental inertial stability is relatively weak and thus the TC’s axisymmetric secondary circulation is more sensitive to the upper-level asymmetric forcing.

    The present study focused primarily on the dynamical effect of interaction in terms of EFC and inertial stability.In fact,a westerly trough approaching a TC will also import cold and dry air—also known as cold air intrusion.Such intrusion of cold and dry air will certainly alter the thermodynamic structure of a TC and exert some influence on its axisymmetric secondary circulation.Molinari and Vollaro (1990)found that the eddy heat flux induced by cold advection could contribute in the same direction as eddy momentum flux within 500 km,but opposite outside.Therefore,Molinari et al.(1995,1998)later analyzed interactions by adopting the Eliassen–Palm(EP)flux diagnostic,which has two components:eddy angular momentum flux(dynamical component)and eddy heat flux(thermodynamical component).As there is uncertainty as to whether the contribu-tions of these two components could be the same or opposite(e.g.,Molinari and Vollaro,1990;Qian et al.,2011),one way to address this issue is to use the EP flux diagnostic and its divergence.This diagnostic tool also requires the whole vertical structure of the TC,rather than a single upper-layer wind field,so that it will provide greater insights into TC–environment interaction.Using such a diagnostic in a statistical study of interaction,rather than EFC alone,may help identify“good”troughs and favorable interactions.

    Fig.10.Six selected interaction TCs(one row for each).The left-hand column shows the tracks for each case and 200-hPa flow fields(vectors;m s-1)valid at interaction times labelled in the titles.Black dots show the TC positions at those times.The right-hand column shows the temporal evolution of TC intensity i.e.,minimum SLP(black line;hPa),as well as the three environmental factors of SST(red line;?C),VWS (blue line;m s-1)and EFC(green line;m s-1d-1).Vertical gray lines in the right-hand column indicate the time shown in the titles of the left-hand column.The green shaded areas emphasize the interaction periods of interest with EFC larger than the 10 m s-1d-1threshold.

    Acknowledgements.This study was jointly sponsored by the National Natural Science Foundation of China(Grant Nos. 41205032,41405048 and 41375050),China Special Fund for MeteorologicalResearchin thePublicInterest(GrantNo. GYHY201406008),and the Guangdong Marine Disaster Emergency Response Technology Research Center(2012A032100004). The authors gratefully acknowledge the use of the High Performance Computing Cluster at the South China Sea Institute of Oceanology,Chinese Academy of Sciences.

    REFERENCES

    Barcikowska,M.,F.Feser,and H.von Storch,2012:Usability of best track data in climate statistics in the western North Pacific.Mon.Wea.Rev.,140,2818–2830.

    Bosart,L.F.,and J.A.Bartlo,1991:Tropical storm formation in a baroclinic environment.Mon.Wea.Rev.,119,1979–2013.

    Bosart,L.F.,C.S.Velden,W.E.Bracken,J.Molinari,and P.G. Black,2000:Environmental influences on the rapid intensification of Hurricane Opal(1995)over the Gulf of Mexico. Mon.Wea.Rev.,128,322–352.

    Bracken,W.E.,and L.F.Bosart,2000:The role of synopticscale flow during tropical cyclogenesis over the North Atlantic Ocean.Mon.Wea.Rev.,128,353–376.

    Challa,M.,and R.L.Pfeffer,1980:Effects of eddy fluxes of angular momentum on model hurricane development.J.Atmos. Sci.,37,1603–1618.

    Chan,J.C.L.,F.M.F.Ko,and Y.M.Lei,2002:Relationship between potential vorticity tendency and tropical cyclone motion.J.Atmos.Sci.,59,1317–1336.

    Chen,X.M.,Y.Q.Wang,and K.Zhao,2015:Synoptic flow patterns and large-scale characteristics associated with rapidly intensifying tropical cyclones in the South China Sea.Mon. Wea.Rev.,143,64–87.

    Choi,Y.,K.-S.Yun,K.-J.Ha,K.-Y.Kim,S.-J.Yoon,and J. C.L.Chan,2013:Effects of asymmetric SST distribution on straight-moving Typhoon Ewiniar(2006)and recurving Typhoon Maemi(2003).Mon.Wea.Rev.,141,3950–3967.

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

    DeMaria,M.,J.-J.Baik,and J.Kaplan,1993:Upper-level eddy angular momentum fluxes and tropical cyclone intensity change.J.Atmos.Sci.,50,1133–1147.

    DeMaria,M.,J.A.Knaff,and C.Sampson,2007:Evaluation of long-term trends in tropical cyclone intensity forecasts.Meteor.Atmos.Phys.,97,19–28.

    Eliassen,A.,1952:Slow thermally or frictionally controlled meridional circulation in a circular vortex.Astrophysica Norvegica,5,19–60.

    Elsberry,R.L.,G.J.Holland,H.Gerrish,M.DeMaria,C.P. Guard,and K.A.Emanuel,1992:Is there any hope for tropical cyclone intensity prediction?—A panel discussion.Bull. Amer.Meteor.Soc.,73,264–275.

    Emanuel,K.A.,1986:An air-sea interaction theory for tropical cyclones.Part I:Steady-state maintenance.J.Atmos.Sci.,43, 585–605.

    Emanuel,K.A.,2000:A statistical analysis of tropical cyclone intensity.Mon.Wea.Rev.,128,1139–1152.

    Erickson,C.O.,1967:Some aspects of the development of Hurricane Dorothy.Mon.Wea.Rev.,95,121–130.

    Fitzpatrick,P.J.,1997:Understanding and forecasting tropical cyclone intensity change with the typhoon intensity prediction scheme(TIPS).Wea.Forecasting,12,826–846.

    Fitzpatrick,P.J.,J.A.Knaff,C.W.Landsea,and S.V.Finley, 1995:Documentation of a systematic bias in the aviation model’s forecast of the Atlantic tropical upper-tropospheric trough:Implications for tropical cyclone forecasting.Wea. Forecasting,10,433–446.

    Hanley,D.E.,2002:The evolution of a hurricane-trough interaction from a satellite perspective.Wea.Forecasting,17,916–926.

    Hanley,D.,J.Molinari,and D.Deyser,2001:A composite study of the interactions between tropical cyclones and uppertropospheric troughs.Mon.Wea.Rev.,129,2570–2584.

    Holland,G.J.,and R.T.Merrill,1984:On the dynamics of tropical cyclone structural changes.Quart.J.Roy.Meteor.Soc., 110,723–745.

    Hoskins,B.J.,M.E.McIntyre,and A.W.Robertson,1985:On the use and significance of isentropic potential vorticitymaps. Quart.J.Roy.Meteor.Soc.,111,877-946.

    Kimball,S.K.,and J.L.Evans,2002:Idealized numerical simulations of hurricane-trough interaction.Mon.Wea.Rev.,130, 2210–2227.

    Leroux,M.-D.,M.Plu,D.Barbary,F.Roux,and P.Arbogast, 2013:Dynamical and physical processes leading to tropical cyclone intensification under upper-level trough forcing.J. Atmos.Sci.,70,2547–2565.

    Lewis,B.M.,and D.P.Jorgensen,1978:Study of the dissipation of Hurricane Gertrude(1974).Mon.Wea.Rev.,106,1288–1306.

    Li,Y.,L.S.Chen,and X.T.Lei,2006:Numerical study on impacts of upper-level westerly trough on the extratropical transitionprocess of Typhoon Winnie(1997).ActaMeteorologica Sinica,64,552–563.(in Chinese)

    Martin,J.D.,and W.M.Gray,1993:Tropical cyclone observation and forecasting with and without aircraft reconnaissance. Wea.Forecasting,8,519–532.

    McTaggart-Cowan,R.,L.F.Bosart,C.A.Davis,E.H.Atallah,J. R.Gyakum,and K.A.Emanuel,2006:Analysis of Hurricane Catarina(2004).Mon.Wea.Rev.,134,3029–3053.

    Merrill,R.T.,1988a:Characteristics of the upper-tropospheric environmental flow around hurricanes.J.Atmos.Sci.,45,1665–1677.

    Merrill,R.T.,1988b:Environmental influences on hurricane intensification.J.Atmos.Sci.,45,1678–1687.

    Molinari,J.,and D.Vollaro,1989:External influences on hurricane intensity.Part I:Outflow layer eddy angular momentum fluxes.J.Atmos.Sci.,46,1093–1104.

    Molinari,J.,and D.Vollaro,1990:External influences on hurricane intensity.Part II:Vertical structure and response of the hurricane vortex.J.Atmos.Sci.,47,1902–1918.

    Molinari,J.,and D.Vollaro,1993:Environmental controls on eye wall cycles and intensity change in Hurricane Allen(1980). Tropical Cyclone Disasters,J.Lighthill et al.,Eds.,Peking University Press,328–337.

    Molinari,J.,D.Vollaro,and F.Robasky,1992:Use of ECMWF operational analyses for studies of the tropical cyclone environment.Meteor.Atmos.Phys.,47,127–144.

    Molinari,J.,S.Skubis,and D.Vollaro,1995:External influences on hurricane intensity.Part III:Potential vorticity structure.J. Atmos.Sci.,52,3593–3606.

    Molinari,J.,S.Skubis,D.Vollaro,F.Alsheimer,and H.E. Willoughby,1998:Potential vorticity analysis of tropical cyclone intensification.J.Atmos.Sci.,55,2632–2644.

    Molinari,J.,P.Dodge,D.Vollaro,and K.L.Corbosiero,2006: Mesoscale aspects of the downshear reformation of a tropical cyclone.J.Atmos.Sci.,63,341–354.

    Montgomery,M.T.,and R.K.Smith,2014:Paradigms for tropical cyclone intensification.Australian Meteorological and Oceanographic Journal,64,37–66.

    Pfeffer,R.L.,and M.Challa,1981:A numerical study of the role of eddy fluxes of momentum in the development of Atlantic hurricanes.J.Atmos.Sci.,38,2393–2398.

    Qian,Y.K.,C.X.Liang,Q.Q.Liang,L.X.Lin,and Z.J.Yuan, 2011:On the forced tangentially-averaged radial-vertical circulation within vortices.Part II:The transformation of Tropical Storm Haima(2004).Adv.Atmos.Sci.,28,1143–1158, doi:10.1007/s00376-010-0060-x.

    Rappin,E.D.,M.C.Morgan,and G.J.Tripoli,2011:The impact ofoutflow environment ontropicalcyclone intensificationand structure.J.Atmos.Sci.,68,177–194.

    Rodgers,E.B.,S.W.Chang,J.Stout,J.Steranka,and J.-J.Shi, 1991:Satellite observations of variations in tropical cyclone convection caused by upper-tropospheric troughs.J.Appl. Meteor.,30,1163–1184.

    Rodgers,E.B.,W.S.Olson,V.M.Karyampudi,and H.F.Pierce, 1998:Satellite-derived latent heating distribution and environmental influences in Hurricane Opal(1995).Mon.Wea. Rev.,126,1229–1247.

    Sawyer,J.S.,1956:The vertical circulation at meteorological fronts and its relation to frontogenesis.Proc.Roy.Soc.London,234A,346–362.

    Shi,J.J.,S.Chang,and S.Raman,1997:Interaction between Hurricane Florence(1988)and an upper-tropospheric westerly trough.J.Atmos.Sci.,54,1231–1247.

    Smith,R.K.,and M.T.Montgomery,2015:Toward clarity on understanding tropical cyclone intensification.J.Atmos.Sci., 72,3020–3031.

    Sundqvist,H.,1970:Numerical simulation of the development of tropical cyclones with a ten-level model.Part I.Tellus,22, 359–390.

    Titley,D.W.,and R.L.Elsberry,2000:Large intensity changes in tropical cyclones:A case study of Supertyphoon Flo during TCM-90.Mon.Wea.Rev.,128,3556–3573.

    Willoughby,H.E.,J.A.Clos,and M.G.Shoreibah,1982:Concentric eye walls,secondary wind maxima,and the evolution of the hurricane vortex.J.Atmos.Sci.,39,395–411.

    Wu,C.-C.,and H.-J.Cheng,1999:An observational study of environmental influences on the intensity changes of Typhoons Flo(1990)and Gene(1990).Mon.Wea.Rev.,127,3003–3031.

    Wu,L.G.,and B.Wang,2000:A potential vorticitytendency diagnostic approach for tropical cyclone motion.Mon.Wea.Rev., 128,1899–1911.

    Yu,H.,and H.J.Kwon,2005:Effect of TC-trough interaction on the intensity change of two typhoons.Wea.Forecasting,20, 199–211.

    Zeng,Z.H.,Y.Q.Wang,and C.-C.Wu,2007:Environmental dynamical control of tropical cyclone intensity—An observational study.Mon.Wea.Rev.,135,38–59.

    Qian,Y.-K.,C.-X.Liang,Z.J.Yuan,S.Q.Peng,J.J.Wu,and S.H.Wang,2016:Upper-tropospheric environment–tropical cyclone interactions over the western North Pacific:A statistical study.Adv.Atmos.Sci.,33(5),614–631,

    10.1007/s00376-015-5148-x.

    13 June 2015;revised 16 September 2015;accepted 22 October 2015)

    ?Yu-Kun QIAN

    Email:qianyk@mail3.sysu.edu.cn

    亚洲精品成人av观看孕妇| 水蜜桃什么品种好| 一级毛片久久久久久久久女| 精品一区二区免费观看| 99热国产这里只有精品6| www.色视频.com| 欧美精品国产亚洲| 又粗又硬又长又爽又黄的视频| 欧美精品国产亚洲| 我要看日韩黄色一级片| 久久人人爽av亚洲精品天堂 | 麻豆国产97在线/欧美| 最近最新中文字幕免费大全7| 久久久久久久久大av| 777米奇影视久久| 国产视频首页在线观看| 听说在线观看完整版免费高清| 日本爱情动作片www.在线观看| 国精品久久久久久国模美| 18禁在线无遮挡免费观看视频| 老司机影院成人| 在线免费观看不下载黄p国产| 三级国产精品片| 免费看av在线观看网站| 国产欧美日韩一区二区三区在线 | 日韩欧美一区视频在线观看 | 一区二区三区乱码不卡18| 老司机影院成人| 国产成人aa在线观看| 特大巨黑吊av在线直播| 在线观看人妻少妇| 久久这里有精品视频免费| 欧美高清成人免费视频www| 大香蕉久久网| 欧美区成人在线视频| 永久免费av网站大全| 亚洲真实伦在线观看| 国产一区二区在线观看日韩| www.色视频.com| 亚洲自偷自拍三级| 97在线视频观看| 免费观看无遮挡的男女| 亚洲综合精品二区| 国产精品一区www在线观看| 欧美成人精品欧美一级黄| 成人亚洲精品av一区二区| 一级爰片在线观看| 亚洲在线观看片| 欧美变态另类bdsm刘玥| 欧美高清成人免费视频www| 免费大片黄手机在线观看| 欧美日韩综合久久久久久| 男的添女的下面高潮视频| 久久热精品热| 国产一区二区三区av在线| 亚洲成人中文字幕在线播放| 美女国产视频在线观看| 少妇的逼好多水| 深爱激情五月婷婷| 色播亚洲综合网| 久久韩国三级中文字幕| 亚洲最大成人手机在线| 在线 av 中文字幕| 免费在线观看成人毛片| 国产精品久久久久久精品电影小说 | 欧美变态另类bdsm刘玥| 黄色一级大片看看| 亚洲精品视频女| 日韩制服骚丝袜av| 丰满乱子伦码专区| 美女主播在线视频| 国产美女午夜福利| 两个人的视频大全免费| 亚洲国产精品国产精品| 国产成年人精品一区二区| 国产91av在线免费观看| 少妇高潮的动态图| 最近2019中文字幕mv第一页| 自拍欧美九色日韩亚洲蝌蚪91 | 久久这里有精品视频免费| 欧美激情国产日韩精品一区| 国产高清不卡午夜福利| 麻豆久久精品国产亚洲av| 欧美日韩视频精品一区| 亚洲精品色激情综合| 视频中文字幕在线观看| 性色avwww在线观看| 狂野欧美激情性bbbbbb| 欧美高清性xxxxhd video| av专区在线播放| 国产毛片a区久久久久| 免费观看av网站的网址| 久久久亚洲精品成人影院| 欧美日韩在线观看h| 99久久中文字幕三级久久日本| av在线天堂中文字幕| 性插视频无遮挡在线免费观看| 中文乱码字字幕精品一区二区三区| 91精品国产九色| 国产大屁股一区二区在线视频| 久久久欧美国产精品| 内射极品少妇av片p| av.在线天堂| 久久国产乱子免费精品| 国产欧美另类精品又又久久亚洲欧美| 久久久亚洲精品成人影院| 日本色播在线视频| 爱豆传媒免费全集在线观看| 深爱激情五月婷婷| 毛片女人毛片| 国产69精品久久久久777片| 国产亚洲av片在线观看秒播厂| 在线观看一区二区三区激情| 亚洲,欧美,日韩| 国产免费福利视频在线观看| 国产一区二区在线观看日韩| 美女国产视频在线观看| 深爱激情五月婷婷| 国产精品熟女久久久久浪| 极品教师在线视频| 亚洲三级黄色毛片| 看免费成人av毛片| 国产女主播在线喷水免费视频网站| 自拍偷自拍亚洲精品老妇| 丰满人妻一区二区三区视频av| 亚洲av一区综合| 亚洲欧美中文字幕日韩二区| 高清在线视频一区二区三区| 欧美老熟妇乱子伦牲交| 欧美日韩综合久久久久久| 亚洲自拍偷在线| 七月丁香在线播放| 一区二区三区四区激情视频| 最近2019中文字幕mv第一页| 国产探花在线观看一区二区| 伦理电影大哥的女人| 男女啪啪激烈高潮av片| 99久久中文字幕三级久久日本| 欧美zozozo另类| 国产视频内射| 最近最新中文字幕大全电影3| 看免费成人av毛片| 亚洲精品亚洲一区二区| 免费电影在线观看免费观看| 国产免费视频播放在线视频| 少妇 在线观看| 丝袜喷水一区| 内射极品少妇av片p| 精品人妻偷拍中文字幕| 亚洲av成人精品一二三区| 在线观看国产h片| 毛片女人毛片| 日韩一区二区三区影片| 精品久久久噜噜| 国产黄片美女视频| 免费看a级黄色片| 国产国拍精品亚洲av在线观看| 午夜激情久久久久久久| 啦啦啦在线观看免费高清www| 亚洲综合精品二区| 成年av动漫网址| 色5月婷婷丁香| 中文字幕av成人在线电影| 成人国产av品久久久| 欧美高清性xxxxhd video| 欧美日本视频| 人妻夜夜爽99麻豆av| 精品人妻视频免费看| 国产真实伦视频高清在线观看| 黑人高潮一二区| 亚洲婷婷狠狠爱综合网| 国模一区二区三区四区视频| 亚洲成人久久爱视频| 国产精品一区www在线观看| 亚洲久久久久久中文字幕| 成人高潮视频无遮挡免费网站| 免费看a级黄色片| 国产亚洲午夜精品一区二区久久 | 国产精品久久久久久久电影| 一区二区av电影网| 国产高清有码在线观看视频| 在线免费观看不下载黄p国产| 欧美潮喷喷水| 国产熟女欧美一区二区| 久久精品国产自在天天线| 日本爱情动作片www.在线观看| 久久精品久久久久久久性| 一级毛片aaaaaa免费看小| 成年女人看的毛片在线观看| 草草在线视频免费看| 丝袜美腿在线中文| 亚洲精品456在线播放app| 91精品一卡2卡3卡4卡| 欧美精品人与动牲交sv欧美| 一本色道久久久久久精品综合| 王馨瑶露胸无遮挡在线观看| 99久久人妻综合| 亚洲成色77777| 亚洲欧美成人综合另类久久久| 亚洲色图综合在线观看| 国产一级毛片在线| 日本免费在线观看一区| 国内少妇人妻偷人精品xxx网站| 亚洲av.av天堂| 精品久久久久久久末码| 成人亚洲精品av一区二区| 国产成人freesex在线| 少妇人妻一区二区三区视频| 激情五月婷婷亚洲| 午夜爱爱视频在线播放| 一级爰片在线观看| 成人亚洲精品一区在线观看 | 国产精品一区二区在线观看99| 97人妻精品一区二区三区麻豆| 18禁在线播放成人免费| 午夜激情福利司机影院| 国产在线男女| 久久精品人妻少妇| 久久久久久国产a免费观看| 熟女人妻精品中文字幕| 免费人成在线观看视频色| 国产高潮美女av| 久久人人爽人人爽人人片va| 中文乱码字字幕精品一区二区三区| 男女那种视频在线观看| 九九在线视频观看精品| 日韩制服骚丝袜av| 亚洲婷婷狠狠爱综合网| 免费黄网站久久成人精品| 黄色一级大片看看| 91aial.com中文字幕在线观看| 免费大片18禁| 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 成人毛片60女人毛片免费| 国产av不卡久久| 国产精品99久久99久久久不卡 | 午夜激情福利司机影院| 国产伦精品一区二区三区四那| 亚洲av中文av极速乱| a级毛片免费高清观看在线播放| 97精品久久久久久久久久精品| 亚洲精华国产精华液的使用体验| av在线蜜桃| 九色成人免费人妻av| 亚洲精品日本国产第一区| 日韩欧美一区视频在线观看 | 亚洲精品456在线播放app| 岛国毛片在线播放| 又大又黄又爽视频免费| 亚洲国产精品成人综合色| 成年av动漫网址| .国产精品久久| 亚洲一级一片aⅴ在线观看| 婷婷色综合www| 国产成人a区在线观看| 永久免费av网站大全| 久久久久久久午夜电影| 少妇人妻精品综合一区二区| 日韩在线高清观看一区二区三区| 男女边摸边吃奶| 99九九线精品视频在线观看视频| 美女高潮的动态| 特级一级黄色大片| 99热国产这里只有精品6| av国产久精品久网站免费入址| av黄色大香蕉| 免费av不卡在线播放| 亚州av有码| 麻豆成人午夜福利视频| 秋霞伦理黄片| 久久久午夜欧美精品| 一级二级三级毛片免费看| 免费看av在线观看网站| 欧美+日韩+精品| 99热国产这里只有精品6| 久久6这里有精品| 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线 | 最近中文字幕2019免费版| 亚洲av.av天堂| 日韩免费高清中文字幕av| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱码久久久久久按摩| 成人特级av手机在线观看| 成人二区视频| 国产成人精品福利久久| 日本wwww免费看| 毛片女人毛片| 国产在线男女| 久久精品国产自在天天线| 婷婷色综合www| 国产在线男女| 美女主播在线视频| 免费观看的影片在线观看| 久久久久久久亚洲中文字幕| 国内少妇人妻偷人精品xxx网站| 久久久久久国产a免费观看| 免费在线观看成人毛片| 在线精品无人区一区二区三 | 色吧在线观看| 午夜福利视频1000在线观看| 男人添女人高潮全过程视频| 男人舔奶头视频| 久久6这里有精品| 国产色爽女视频免费观看| 亚洲精品aⅴ在线观看| 青青草视频在线视频观看| 熟妇人妻不卡中文字幕| 欧美xxⅹ黑人| 在线观看国产h片| 成人亚洲精品av一区二区| 人体艺术视频欧美日本| 国产国拍精品亚洲av在线观看| 三级国产精品片| 久久鲁丝午夜福利片| 久久99热这里只有精品18| 精品久久久精品久久久| 久久这里有精品视频免费| av免费观看日本| 亚洲aⅴ乱码一区二区在线播放| 新久久久久国产一级毛片| 国产免费福利视频在线观看| 精品一区在线观看国产| 免费看日本二区| 一级毛片黄色毛片免费观看视频| 国产成人福利小说| 成人二区视频| 美女脱内裤让男人舔精品视频| 亚洲精品色激情综合| 久久久久国产精品人妻一区二区| 亚洲在久久综合| 日韩大片免费观看网站| 中文乱码字字幕精品一区二区三区| 国产成人一区二区在线| 精品国产三级普通话版| 高清欧美精品videossex| 亚洲精华国产精华液的使用体验| 亚洲欧洲国产日韩| 国产成人aa在线观看| 成人高潮视频无遮挡免费网站| 亚洲国产欧美在线一区| 韩国高清视频一区二区三区| 夫妻性生交免费视频一级片| 国产av不卡久久| 欧美 日韩 精品 国产| 我的老师免费观看完整版| 亚洲精品久久午夜乱码| 免费黄频网站在线观看国产| 2021天堂中文幕一二区在线观| 18禁在线无遮挡免费观看视频| 亚洲人成网站在线观看播放| 美女主播在线视频| 欧美日韩国产mv在线观看视频 | 国产伦理片在线播放av一区| 九九久久精品国产亚洲av麻豆| 水蜜桃什么品种好| 日韩,欧美,国产一区二区三区| 我的女老师完整版在线观看| 男人爽女人下面视频在线观看| 搞女人的毛片| 日韩三级伦理在线观看| 99热这里只有是精品在线观看| 一个人看的www免费观看视频| 丰满乱子伦码专区| 晚上一个人看的免费电影| 国产成人精品久久久久久| 深爱激情五月婷婷| tube8黄色片| 免费看a级黄色片| 日日摸夜夜添夜夜爱| 国产色婷婷99| 国产黄色免费在线视频| 国模一区二区三区四区视频| 午夜视频国产福利| 亚洲av成人精品一二三区| 国产精品99久久99久久久不卡 | 亚洲欧美日韩另类电影网站 | 国产色婷婷99| 婷婷色综合大香蕉| 久久久久精品性色| 成人高潮视频无遮挡免费网站| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 交换朋友夫妻互换小说| 在线观看免费高清a一片| 男女无遮挡免费网站观看| av免费观看日本| 黄色视频在线播放观看不卡| 久久精品国产亚洲av涩爱| 国产精品人妻久久久影院| 蜜桃久久精品国产亚洲av| 国产永久视频网站| 少妇猛男粗大的猛烈进出视频 | 国产国拍精品亚洲av在线观看| 亚洲国产成人一精品久久久| 国产精品久久久久久久电影| 听说在线观看完整版免费高清| 久久精品国产鲁丝片午夜精品| 大话2 男鬼变身卡| 亚洲,一卡二卡三卡| 午夜免费观看性视频| 亚洲国产av新网站| 亚洲欧美一区二区三区国产| 黄色日韩在线| 久久久久久久午夜电影| 极品少妇高潮喷水抽搐| av播播在线观看一区| 人妻少妇偷人精品九色| 一个人看的www免费观看视频| 久久久精品94久久精品| 久久久久国产精品人妻一区二区| 人体艺术视频欧美日本| 日韩成人av中文字幕在线观看| 精品视频人人做人人爽| 最新中文字幕久久久久| 秋霞在线观看毛片| 精品熟女少妇av免费看| 午夜福利视频精品| 久久亚洲国产成人精品v| 久久久久久久久大av| 嘟嘟电影网在线观看| 男女边摸边吃奶| 听说在线观看完整版免费高清| 免费看日本二区| 亚洲欧美一区二区三区黑人 | 亚洲婷婷狠狠爱综合网| 亚洲性久久影院| xxx大片免费视频| 欧美成人午夜免费资源| 精品人妻视频免费看| 最近2019中文字幕mv第一页| 免费观看av网站的网址| 亚洲,一卡二卡三卡| 日韩在线高清观看一区二区三区| 国产成年人精品一区二区| 免费黄网站久久成人精品| 只有这里有精品99| 男人爽女人下面视频在线观看| 人人妻人人澡人人爽人人夜夜| 欧美一区二区亚洲| 成年女人在线观看亚洲视频 | 亚洲av在线观看美女高潮| 国产一区二区三区av在线| 97精品久久久久久久久久精品| 欧美日韩一区二区视频在线观看视频在线 | 国产高清国产精品国产三级 | 成人一区二区视频在线观看| 亚洲欧美成人综合另类久久久| 亚洲成色77777| 久久99精品国语久久久| 好男人在线观看高清免费视频| 一个人看的www免费观看视频| 舔av片在线| 午夜免费观看性视频| xxx大片免费视频| 中文字幕人妻熟人妻熟丝袜美| av免费观看日本| 99热全是精品| 国内少妇人妻偷人精品xxx网站| 天美传媒精品一区二区| 久久这里有精品视频免费| 中文精品一卡2卡3卡4更新| 亚洲国产精品国产精品| 成人黄色视频免费在线看| 亚洲精品成人av观看孕妇| 欧美精品一区二区大全| 日韩欧美精品免费久久| 五月开心婷婷网| 国产综合精华液| 午夜免费男女啪啪视频观看| 欧美精品国产亚洲| av线在线观看网站| 精品人妻视频免费看| 三级国产精品片| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久精品电影小说 | 最近2019中文字幕mv第一页| 国产午夜福利久久久久久| 最近中文字幕高清免费大全6| 91精品一卡2卡3卡4卡| 男女边吃奶边做爰视频| 在线观看人妻少妇| www.色视频.com| 欧美成人午夜免费资源| 国产一区二区在线观看日韩| 热re99久久精品国产66热6| 神马国产精品三级电影在线观看| 男人爽女人下面视频在线观看| 亚洲电影在线观看av| 毛片一级片免费看久久久久| 久久久久久久久久成人| 国产成人福利小说| 欧美97在线视频| 国产一区亚洲一区在线观看| 在线亚洲精品国产二区图片欧美 | 国产av不卡久久| 国产高清三级在线| 久久久成人免费电影| 别揉我奶头 嗯啊视频| 亚洲精品成人久久久久久| 性色av一级| 自拍偷自拍亚洲精品老妇| 99热6这里只有精品| 久久久久久九九精品二区国产| 五月天丁香电影| 国产午夜福利久久久久久| 男插女下体视频免费在线播放| 人人妻人人澡人人爽人人夜夜| 久久99热这里只频精品6学生| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 久久久精品免费免费高清| 精品久久国产蜜桃| 久久这里有精品视频免费| 日韩av免费高清视频| 黄色欧美视频在线观看| 亚洲美女视频黄频| 狠狠精品人妻久久久久久综合| 亚洲精品日韩av片在线观看| 欧美丝袜亚洲另类| 中文字幕av成人在线电影| 国产精品99久久99久久久不卡 | 夜夜爽夜夜爽视频| 亚洲国产欧美在线一区| 国产高清三级在线| 黑人高潮一二区| 男女无遮挡免费网站观看| 男女国产视频网站| 国产精品国产av在线观看| 精品久久久久久久久av| 波多野结衣巨乳人妻| 丰满乱子伦码专区| 亚洲丝袜综合中文字幕| 大码成人一级视频| 男女无遮挡免费网站观看| 一边亲一边摸免费视频| 亚洲三级黄色毛片| 一级片'在线观看视频| 亚洲三级黄色毛片| 国产黄a三级三级三级人| 国产一区有黄有色的免费视频| 成人国产麻豆网| 亚洲av不卡在线观看| av在线蜜桃| 水蜜桃什么品种好| 边亲边吃奶的免费视频| 熟女人妻精品中文字幕| 成人毛片60女人毛片免费| 水蜜桃什么品种好| 亚洲精品亚洲一区二区| 水蜜桃什么品种好| av在线观看视频网站免费| 一级av片app| 看非洲黑人一级黄片| 亚洲精品乱码久久久v下载方式| 精品一区二区免费观看| 亚洲精品一区蜜桃| 一级毛片 在线播放| 黄色欧美视频在线观看| 亚洲国产欧美人成| 国产亚洲一区二区精品| freevideosex欧美| 亚洲人与动物交配视频| 欧美亚洲 丝袜 人妻 在线| 欧美丝袜亚洲另类| 国产真实伦视频高清在线观看| 一级av片app| 久久久久精品久久久久真实原创| 精品亚洲乱码少妇综合久久| 亚洲av.av天堂| 成人美女网站在线观看视频| 色哟哟·www| 尾随美女入室| 国产免费视频播放在线视频| 波野结衣二区三区在线| 国产亚洲av嫩草精品影院| 午夜福利网站1000一区二区三区| 在线免费十八禁| 国产亚洲5aaaaa淫片| 如何舔出高潮| 日日撸夜夜添| 午夜精品一区二区三区免费看| 日韩av免费高清视频| 午夜视频国产福利| 午夜亚洲福利在线播放| 久久久久久久久久久免费av| 18禁裸乳无遮挡免费网站照片| 在线看a的网站| 麻豆乱淫一区二区| 97人妻精品一区二区三区麻豆| 高清毛片免费看| 七月丁香在线播放| 成人欧美大片| 免费大片18禁| 亚洲av二区三区四区| 在现免费观看毛片| 婷婷色麻豆天堂久久| 欧美3d第一页| 少妇丰满av| 欧美一级a爱片免费观看看| a级毛片免费高清观看在线播放| 色5月婷婷丁香| 久久午夜福利片| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区国产| 深夜a级毛片| av黄色大香蕉| 色哟哟·www| 特级一级黄色大片| 丰满少妇做爰视频| 在线观看av片永久免费下载| a级毛色黄片| 草草在线视频免费看| 国产免费又黄又爽又色|