• <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片在线观看| 91午夜精品亚洲一区二区三区 | 一区二区三区高清视频在线| 精品不卡国产一区二区三区| 欧美不卡视频在线免费观看| 亚洲av五月六月丁香网| 一个人看的www免费观看视频| a级一级毛片免费在线观看| 脱女人内裤的视频| 九色国产91popny在线| 久久久久久久亚洲中文字幕 | 国产精品一及| 亚洲,欧美精品.| 丰满人妻一区二区三区视频av| 亚洲,欧美精品.| 精品久久久久久久久av| а√天堂www在线а√下载| 午夜免费成人在线视频| 国产精品人妻久久久久久| 欧美最黄视频在线播放免费| 国产欧美日韩精品一区二区| 中文亚洲av片在线观看爽| 男人舔女人下体高潮全视频| 啪啪无遮挡十八禁网站| 最新在线观看一区二区三区| 9191精品国产免费久久| 亚洲精品在线观看二区| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩福利视频一区二区| 99精品在免费线老司机午夜| 1000部很黄的大片| 亚洲熟妇熟女久久| 国产亚洲av嫩草精品影院| 亚洲男人的天堂狠狠| 内射极品少妇av片p| 成人av一区二区三区在线看| 变态另类丝袜制服| 成人性生交大片免费视频hd| 可以在线观看毛片的网站| 久久精品91蜜桃| 久久99热这里只有精品18| 亚洲av熟女| 精品日产1卡2卡| 熟妇人妻久久中文字幕3abv| 精品人妻视频免费看| 乱人视频在线观看| 欧美黄色淫秽网站| 老司机深夜福利视频在线观看| av天堂在线播放| 久久久久久国产a免费观看| 简卡轻食公司| 赤兔流量卡办理| 三级国产精品欧美在线观看| 国产国拍精品亚洲av在线观看| 精品久久久久久久末码| 亚洲avbb在线观看| 国产精品1区2区在线观看.| 欧美bdsm另类| 亚洲精品影视一区二区三区av| 亚洲av不卡在线观看| 久久精品国产亚洲av涩爱 | 亚洲第一电影网av| 亚洲经典国产精华液单 | 99热这里只有是精品50| 亚洲一区二区三区色噜噜| 十八禁网站免费在线| 国产淫片久久久久久久久 | 国产精品久久久久久精品电影| 99久久无色码亚洲精品果冻| 亚洲av电影不卡..在线观看| 亚洲av五月六月丁香网| 国产精品免费一区二区三区在线| 国产精品久久视频播放| 国产大屁股一区二区在线视频| 久久久久免费精品人妻一区二区| 亚洲国产精品999在线| 淫秽高清视频在线观看| 色视频www国产| 91在线观看av| 午夜影院日韩av| 美女xxoo啪啪120秒动态图 | 在线观看一区二区三区| 长腿黑丝高跟| 偷拍熟女少妇极品色| 亚洲五月天丁香| 久久亚洲精品不卡| 好男人电影高清在线观看| 久久国产乱子伦精品免费另类| 免费大片18禁| 国产久久久一区二区三区| 国产乱人伦免费视频| 久久亚洲精品不卡| 亚洲av免费高清在线观看| 亚洲欧美日韩高清在线视频| 日韩人妻高清精品专区| 亚洲狠狠婷婷综合久久图片| 免费看日本二区| 免费看光身美女| 精品久久久久久久久亚洲 | 欧美中文日本在线观看视频| 99久久九九国产精品国产免费| 男女做爰动态图高潮gif福利片| 中亚洲国语对白在线视频| 国产av一区在线观看免费| 夜夜看夜夜爽夜夜摸| 亚洲精品色激情综合| 亚洲成av人片免费观看| 国产精品电影一区二区三区| 日韩精品青青久久久久久| 国产三级黄色录像| 夜夜夜夜夜久久久久| 欧美激情国产日韩精品一区| 亚洲美女视频黄频| 国产精华一区二区三区| 精品国内亚洲2022精品成人| av中文乱码字幕在线| 黄色视频,在线免费观看| 欧美黄色淫秽网站| 又紧又爽又黄一区二区| 国产一区二区激情短视频| 少妇人妻精品综合一区二区 | 天堂√8在线中文| 欧美成狂野欧美在线观看| 国产精品电影一区二区三区| 2021天堂中文幕一二区在线观| 免费看美女性在线毛片视频| 午夜福利成人在线免费观看| 在线观看免费视频日本深夜| 亚洲精品色激情综合| 亚洲精品影视一区二区三区av| 国产毛片a区久久久久| 真人做人爱边吃奶动态| 搡老熟女国产l中国老女人| 国产探花极品一区二区| 国产一级毛片七仙女欲春2| 一级黄色大片毛片| 国产爱豆传媒在线观看| 1000部很黄的大片| 欧美在线黄色| 免费人成视频x8x8入口观看| 久久精品国产亚洲av天美| 又紧又爽又黄一区二区| 欧美日韩中文字幕国产精品一区二区三区| 色视频www国产| 国产白丝娇喘喷水9色精品| 国产欧美日韩一区二区三| 禁无遮挡网站| 午夜日韩欧美国产| 午夜a级毛片| 久久久久国内视频| netflix在线观看网站| 日韩欧美精品v在线| 制服丝袜大香蕉在线| 亚洲久久久久久中文字幕| 成年女人看的毛片在线观看| 露出奶头的视频| 亚洲第一欧美日韩一区二区三区| 波野结衣二区三区在线| 免费高清视频大片| 免费看日本二区| 在线天堂最新版资源| 精品久久国产蜜桃| 精品日产1卡2卡| 动漫黄色视频在线观看| 国产一区二区激情短视频| 久久草成人影院| h日本视频在线播放| 88av欧美| 极品教师在线视频| 女同久久另类99精品国产91| or卡值多少钱| 在线播放国产精品三级| 91午夜精品亚洲一区二区三区 | 美女cb高潮喷水在线观看| 中文字幕高清在线视频| 国内揄拍国产精品人妻在线| 亚洲成人久久爱视频| 欧美精品啪啪一区二区三区| 五月玫瑰六月丁香| 首页视频小说图片口味搜索| 在线a可以看的网站| 少妇的逼好多水| 国产视频内射| 99久久精品一区二区三区| 1000部很黄的大片| 久久久久久久久中文| 久久精品国产亚洲av涩爱 | 校园春色视频在线观看| 婷婷亚洲欧美| 悠悠久久av| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 一本综合久久免费| 啦啦啦韩国在线观看视频| 波野结衣二区三区在线| 九九久久精品国产亚洲av麻豆| 国模一区二区三区四区视频| 午夜福利在线观看吧| 亚洲一区高清亚洲精品| 免费在线观看日本一区| 日本熟妇午夜| 亚洲综合色惰| 青草久久国产| 成人美女网站在线观看视频| 动漫黄色视频在线观看| 久久久久久久午夜电影| 夜夜夜夜夜久久久久| 欧美+日韩+精品| 性插视频无遮挡在线免费观看| 九九在线视频观看精品| 日本 av在线| 岛国在线免费视频观看| 一级作爱视频免费观看| 亚洲人成电影免费在线| 级片在线观看| 欧美不卡视频在线免费观看| 久久99热这里只有精品18| 精品久久久久久久末码| 在线看三级毛片| 色综合婷婷激情| av中文乱码字幕在线| 日韩欧美免费精品| 内地一区二区视频在线| av福利片在线观看| ponron亚洲| 久久久精品欧美日韩精品| 国产中年淑女户外野战色| 国产在线男女| 国内少妇人妻偷人精品xxx网站| 18禁黄网站禁片免费观看直播| 亚洲av美国av| 久久亚洲真实| 久久久久性生活片| 亚洲在线自拍视频| 日日摸夜夜添夜夜添av毛片 | 看黄色毛片网站| 亚洲精品在线美女| 亚洲av.av天堂| 97热精品久久久久久| 精品人妻熟女av久视频| 又黄又爽又刺激的免费视频.| 国产亚洲精品av在线| 大型黄色视频在线免费观看| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 三级毛片av免费| 久久精品国产99精品国产亚洲性色| 午夜视频国产福利| 亚洲成人免费电影在线观看| 欧美黄色片欧美黄色片| 国产精华一区二区三区| av专区在线播放| 免费搜索国产男女视频| 三级毛片av免费| 国产av一区在线观看免费| 亚洲三级黄色毛片| 精品一区二区三区视频在线观看免费| 久久久成人免费电影| 日韩av在线大香蕉| 特大巨黑吊av在线直播| 久久国产精品影院| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 好男人在线观看高清免费视频| 亚洲人与动物交配视频| 啦啦啦观看免费观看视频高清| 99视频精品全部免费 在线| 99久久精品热视频| 夜夜爽天天搞| 国产成人av教育| 级片在线观看| 免费av不卡在线播放| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲| 色尼玛亚洲综合影院| 人人妻,人人澡人人爽秒播| 桃红色精品国产亚洲av| 久久国产精品影院| 天天躁日日操中文字幕| 国产精品人妻久久久久久| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 国产高清视频在线播放一区| 99riav亚洲国产免费| 国产精品99久久久久久久久| 综合色av麻豆| av专区在线播放| 人人妻人人澡欧美一区二区| 国产精品,欧美在线| 真人一进一出gif抽搐免费| 免费一级毛片在线播放高清视频| eeuss影院久久| 三级毛片av免费| 国产极品精品免费视频能看的| 欧美zozozo另类| 免费av观看视频| 哪里可以看免费的av片| 国产高清三级在线| 丰满人妻一区二区三区视频av| 午夜福利成人在线免费观看| 18禁黄网站禁片免费观看直播| 精品无人区乱码1区二区| 亚洲一区二区三区色噜噜| 久久亚洲真实| 网址你懂的国产日韩在线| 亚洲av日韩精品久久久久久密| 成人美女网站在线观看视频| 深夜精品福利| 直男gayav资源| av国产免费在线观看| 国产人妻一区二区三区在| 亚洲 欧美 日韩 在线 免费| 18美女黄网站色大片免费观看| 老女人水多毛片| 欧美精品啪啪一区二区三区| 欧美高清成人免费视频www| 在线看三级毛片| 18+在线观看网站| 老熟妇仑乱视频hdxx| 我要搜黄色片| 欧美成人免费av一区二区三区| 欧美绝顶高潮抽搐喷水| 国产三级中文精品| 午夜福利免费观看在线| 久9热在线精品视频| av在线老鸭窝| 在线a可以看的网站| 国产亚洲欧美在线一区二区| 国产精品电影一区二区三区| 国产精品三级大全| 国产精品电影一区二区三区| 久久精品国产清高在天天线| bbb黄色大片| 亚洲国产精品合色在线| 日日摸夜夜添夜夜添小说| 亚洲色图av天堂| 亚洲av.av天堂| 一进一出抽搐动态| 亚洲,欧美,日韩| АⅤ资源中文在线天堂| 亚洲人成网站在线播| 91午夜精品亚洲一区二区三区 | 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 日韩欧美国产在线观看| 最好的美女福利视频网| 欧美性猛交黑人性爽| 欧美色欧美亚洲另类二区| 亚洲成av人片在线播放无| 少妇被粗大猛烈的视频| 亚洲五月婷婷丁香| 久久久久久久久中文| 美女大奶头视频| 午夜免费激情av| 午夜日韩欧美国产| 国内揄拍国产精品人妻在线| 久99久视频精品免费| 欧美不卡视频在线免费观看| 国产精品一及| 桃色一区二区三区在线观看| 婷婷亚洲欧美| 麻豆国产97在线/欧美| 精品久久久久久成人av| 国产乱人视频| 亚洲黑人精品在线| 国产在线男女| 免费电影在线观看免费观看| 欧美3d第一页| 成人精品一区二区免费| 亚洲av电影在线进入| 村上凉子中文字幕在线| 最新中文字幕久久久久| www.www免费av| 真人做人爱边吃奶动态| 亚洲国产精品sss在线观看| 乱人视频在线观看| 老司机午夜十八禁免费视频| 99精品久久久久人妻精品| 亚洲中文字幕日韩| 国产精品三级大全| 成人亚洲精品av一区二区| 国产精品影院久久| 在线观看午夜福利视频| 真人做人爱边吃奶动态| 宅男免费午夜| 亚洲人与动物交配视频| 亚洲国产日韩欧美精品在线观看| 亚洲精品日韩av片在线观看| 男女之事视频高清在线观看| 午夜久久久久精精品| 国产精品久久久久久亚洲av鲁大| 亚洲欧美清纯卡通| 亚洲,欧美,日韩| 国产伦精品一区二区三区视频9| 麻豆av噜噜一区二区三区| 成人永久免费在线观看视频| 欧美bdsm另类| 国产中年淑女户外野战色| 小说图片视频综合网站| 亚洲av熟女| 亚洲中文日韩欧美视频| 波多野结衣高清无吗| 欧美成人a在线观看| 看片在线看免费视频| 九九久久精品国产亚洲av麻豆| 亚洲国产色片| 黄色一级大片看看| 999久久久精品免费观看国产| 赤兔流量卡办理| 在线免费观看的www视频| 日韩欧美在线二视频| 中文字幕人成人乱码亚洲影| 久久久国产成人免费| 变态另类成人亚洲欧美熟女| 成人午夜高清在线视频| 免费在线观看日本一区| 久久久久久久久大av| 老司机午夜福利在线观看视频| 自拍偷自拍亚洲精品老妇| 一进一出抽搐动态| 国产一级毛片七仙女欲春2| 亚洲精品影视一区二区三区av| 亚洲,欧美,日韩| 国产一区二区在线av高清观看| 淫妇啪啪啪对白视频| 国产精品久久久久久精品电影| 国产精品久久久久久久久免 | av在线观看视频网站免费| 精品一区二区三区视频在线| av天堂中文字幕网| 淫秽高清视频在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 在线观看美女被高潮喷水网站 | 无遮挡黄片免费观看| 国产麻豆成人av免费视频| 国产精品综合久久久久久久免费| 97热精品久久久久久| 国产黄片美女视频| 国产伦人伦偷精品视频| 老女人水多毛片| 国产一区二区三区视频了| 99久久无色码亚洲精品果冻| 91在线观看av| 国产黄a三级三级三级人| 白带黄色成豆腐渣| 欧美激情久久久久久爽电影| 69av精品久久久久久| h日本视频在线播放| 俺也久久电影网| 亚洲中文字幕一区二区三区有码在线看| 国产精品嫩草影院av在线观看 | 久久热精品热| 极品教师在线视频| 97超视频在线观看视频| 麻豆一二三区av精品| 日本五十路高清| 国产精品久久久久久久久免 | 91麻豆av在线| 丰满人妻一区二区三区视频av| 美女被艹到高潮喷水动态| 嫩草影院精品99| 日本熟妇午夜| 久久性视频一级片| 国产男靠女视频免费网站| 亚洲avbb在线观看| 美女高潮喷水抽搐中文字幕| netflix在线观看网站| 99国产精品一区二区蜜桃av| 成人特级黄色片久久久久久久| 亚洲人与动物交配视频| 久久精品国产自在天天线| 久久人妻av系列| 综合色av麻豆| 亚洲经典国产精华液单 | 亚洲国产精品sss在线观看| 国产av一区在线观看免费| 九色成人免费人妻av| 特级一级黄色大片| 亚洲欧美清纯卡通| av中文乱码字幕在线| 亚洲不卡免费看| 天堂√8在线中文| 久久久久久久久久成人| 一本久久中文字幕| av中文乱码字幕在线| 久久精品国产亚洲av天美| 99在线人妻在线中文字幕| 国产成+人综合+亚洲专区| 在线观看午夜福利视频| 波多野结衣高清无吗| 免费黄网站久久成人精品 | 欧美日韩国产亚洲二区| 性插视频无遮挡在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 性欧美人与动物交配| 色吧在线观看| 久久久久久久久中文| 免费搜索国产男女视频| ponron亚洲| 老司机福利观看| 国产精品一区二区三区四区久久| 搞女人的毛片| 久久久久九九精品影院| 久久欧美精品欧美久久欧美| 国产精品永久免费网站| 90打野战视频偷拍视频| 欧美最新免费一区二区三区 | 亚洲中文字幕一区二区三区有码在线看| aaaaa片日本免费| 国产主播在线观看一区二区| 免费在线观看日本一区| 岛国在线免费视频观看| 国产高潮美女av| 欧美成人免费av一区二区三区| 五月玫瑰六月丁香| 日韩国内少妇激情av| 亚洲av免费高清在线观看| 国产真实伦视频高清在线观看 | 精品一区二区三区视频在线观看免费| 国产69精品久久久久777片| 国产免费av片在线观看野外av| 欧美潮喷喷水| а√天堂www在线а√下载| 亚洲一区二区三区色噜噜| 男人狂女人下面高潮的视频| 日本精品一区二区三区蜜桃| 国产成+人综合+亚洲专区| 日本一二三区视频观看| 精品人妻视频免费看| 亚洲在线自拍视频| 日日干狠狠操夜夜爽| 悠悠久久av| 欧美另类亚洲清纯唯美| 国产精品日韩av在线免费观看| 亚洲色图av天堂| 国产麻豆成人av免费视频| 婷婷六月久久综合丁香| 国产在线男女| 成人无遮挡网站| 日韩欧美国产一区二区入口| 级片在线观看| 国产精品女同一区二区软件 | 亚洲第一欧美日韩一区二区三区| 波多野结衣高清作品| 人妻制服诱惑在线中文字幕| 1024手机看黄色片| 久久精品久久久久久噜噜老黄 | 日韩欧美国产一区二区入口| 级片在线观看| 制服丝袜大香蕉在线| 精品一区二区三区人妻视频| 午夜福利成人在线免费观看| 国产精品日韩av在线免费观看| 国产精品1区2区在线观看.| 18禁黄网站禁片午夜丰满| 久久婷婷人人爽人人干人人爱| 日韩欧美在线乱码| 成人无遮挡网站| 色在线成人网| 精品久久久久久久久久免费视频| 制服丝袜大香蕉在线| 身体一侧抽搐| 免费黄网站久久成人精品 | av福利片在线观看| 夜夜夜夜夜久久久久| 天堂动漫精品| 免费无遮挡裸体视频| 很黄的视频免费| 97超视频在线观看视频| 亚洲欧美日韩东京热| 麻豆国产av国片精品| 五月伊人婷婷丁香| 小蜜桃在线观看免费完整版高清| www.www免费av| 美女大奶头视频| 国产欧美日韩精品一区二区| 亚洲黑人精品在线| av天堂在线播放| 欧美一级a爱片免费观看看| 国产高清有码在线观看视频| 91午夜精品亚洲一区二区三区 | 亚洲一区二区三区色噜噜| bbb黄色大片| 91av网一区二区| 国产 一区 欧美 日韩| 最近最新中文字幕大全电影3| 久久99热这里只有精品18| 国产在线精品亚洲第一网站| 一个人免费在线观看的高清视频| 国产毛片a区久久久久| 欧美黑人巨大hd| 亚洲 欧美 日韩 在线 免费| 日韩亚洲欧美综合| 99国产极品粉嫩在线观看| 国产v大片淫在线免费观看| 九九久久精品国产亚洲av麻豆| av在线老鸭窝| 露出奶头的视频| 99久久九九国产精品国产免费| 给我免费播放毛片高清在线观看| 日韩精品青青久久久久久| 欧美中文日本在线观看视频| 日本 av在线| 啦啦啦韩国在线观看视频| 日韩欧美免费精品| 久久国产精品影院| 国产精品美女特级片免费视频播放器| 中文字幕人成人乱码亚洲影| 国产精品人妻久久久久久| 真人一进一出gif抽搐免费| 久久99热这里只有精品18|