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

    Dominant Cloud Microphysical Processes of a Torrential Rainfall Event in Sichuan,China

    2015-02-24 03:39:46HUANGYongjieandCUIXiaopeng
    Advances in Atmospheric Sciences 2015年3期

    HUANG Yongjieand CUI Xiaopeng

    1Key Laboratory of Cloud-Precipitation Physics and Severe Storms,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing100029

    2University of Chinese Academy of Sciences,Beijing100049

    3Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science&Technology,Nanjing210044

    Dominant Cloud Microphysical Processes of a Torrential Rainfall Event in Sichuan,China

    HUANG Yongjie1,2and CUI Xiaopeng?1,3

    1Key Laboratory of Cloud-Precipitation Physics and Severe Storms,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing100029

    2University of Chinese Academy of Sciences,Beijing100049

    3Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science&Technology,Nanjing210044

    High-resolution numerical simulation data of a rainstorm triggering debris fow in Sichuan Province of China simulated by theWeather Researchand Forecasting(WRF)Model wereused tostudythedominant cloud microphysical processes of the torrential rainfall.The results showed that:(1)In the strong precipitation period,particle sizes of all hydrometeors increased, and mean-mass diametersof graupel increased themost signifcantly,ascompared withthose intheweakprecipitation period; (2)The terminal velocity of raindrops was the strongest among all hydrometeors,followed by graupel’s,which was much smaller than that of raindrops.Differences between various hydrometeors’terminal velocities in the strong precipitation period were larger than those in the weak precipitation period,which favored relative motion,collection interaction and transformation between the particles.Absolute terminal velocity values of raindrops and graupel were signifcantly greater than those of air upward velocity,and the stronger the precipitation was,the greater the differences between them were; (3)The orders of magnitudes of the various hydrometeors’sources and sinks in the strong precipitation period were larger than those in the weak precipitation period,causing a difference in the intensity of precipitation.Water vapor,cloud water, raindrops,graupel and their exchange processes played a major role in the production of the torrential rainfall,and there were two main processes via which raindrops were generated:abundant water vapor condensed into cloud water and,on the one hand,accretion of cloud water by rain water formed rain water,while on the other hand,accretion of cloud water by graupel formed graupel,and then the melting of graupel formed rain water.

    torrential rainfall,Sichuan,cloud microphysical processes,numerical simulation

    1. Introduction

    Orographictorrentialrainfallcancausefoods,landslides, debris fows and other natural disasters,constituting a great threat to the lives and properties of people,and representing a great challenge to mountain weather forecasters.The Sichuan Basin,with its complex terrain,is located east of the Tibetan Plateau,north of the Yunnan–Guizhou Plateau,and south of the Qinling Heights.It is a storm-proneplace,where debris fows,landslides and other secondary geological disasters occur frequently.Following two major earthquakes—Wenchuanearthquake(2008)andYa’anearthquake(2013)—the situation has worsened because the soil of Sichuan’s complex terrain surrounding these earthquake zones has become much looser,meaning frequent torrential rainfall ismore likely to lead to debris fows and other disasters.Therefore,understandingthe physicalprocesses responsibleforthe developmentof orographictorrentialrainfall is crucial forthe improvement of operational forecasts in this area.

    Sawyer(1956)revealed three important factors for orographic rainfall:larger-scale atmospheric circulation;the interaction of the ambient fow by terrain;and cloud microphysical processes.Smith and Barstad(2004)developed a linear steady-state theory for orographic rainfall and derived a linear precipitation model that can represent many complex processes in orographic rainfall using a relatively simple equation(Crochet et al.,2007;Smith and Evans,2007). However,this linear model cannot include fow blocking, gravity wave breaking,and other non-linear processes,and usesrelativelysimplemicrophysics.Comparedwiththissimple linear model,high-resolution mesoscale models,such as the Weather Research and Forecasting(WRF)model,can better refect nonlinear dynamics,thermodynamics,and de-tailed microphysical processes in complex orographic rainfall,which has been intensively applied to studies of severe weather and torrential rainfall(Ge et al.,2008;Maussion et al.,2010;Smith et al.,2010;Wang et al.,2011;He and Li, 2013;Wang et al.,2013;Li et al.,2014).Luet al.(2009)simulated and analyzed several local heavy rainfall events over the western Sichuan Basin during August 2003 using the AdvancedRegionalEta-coordinateModel(AREM),andshowed that heavyrainfall events overthis area are highly sensitive to the initial local water vaporconditions.Chen et al.(2010)exploredamethodtoimprovetheaccuracyofSichuantorrential rainfall forecasts based on Bayesian decision theory.Chen and Li(2013)used satellite,radar and routine meteorological data to analyze a mesoscale convectivesystem(MCS)and its effects on short-term heavy rainfall in the Sichuan Basin during July 2012.Wang et al.(2013)simulated and analyzed a rainstorm process infuenced by the terrain in Sichuan using the WRF model and carried out a diagnostic analysis using dynamical parameters.Li et al.(2014)analyzed observations and a WRF simulation of a rainstorm process that triggered debris fows,and their results indicated that this precipitation process had a short duration,strong intensity and local characteristics,and the heavy rainfall occurred mainly near steep terrain in the Wenchuan–Ya’an earthquake fault zone.Their study(Li et al.,2014)also noted that,under favorable largescale circulations,the developmentof a southwest vortex and its interaction with strong wind at the western edge of the subtropical high enhanced the lifting effects of steep terrain and caused heavy rainfall.Because of the combined effects of complex dynamic,thermodynamic and cloud microphysical processes over complex terrain,more accurate quantitative precipitation forecasts(QPFs)remain a huge challenge. Studies on orographictorrential rainfall in Sichuan,in particular,and studies on cloudmicrophysicalprocesses associated with torrential rainfallhave seldom been conducted,although some numerical sensitivity experiments of cloud microphysical parameterization schemes have been carried out(Li and He,2011;Liao et al.,2012;Ma et al.,2012).Overall,however,a number of open questions remain.For example,what is the key cloud microphysical process for the formation of orographic torrential rainfall?And what are the main raindrop sources for the production of torrential rainfall?

    Li et al.(2014)carried out a numerical simulation of a rainstorm in Sichuan using the WRF model.They discussed the role of large-scale circulation evolution and development of the mesoscale system in the emergence of the rainstorm over complex terrain from macroscopic and dynamic perspectives.However,since the formation of a rainstorm is related to both macroscopic dynamic and cloud microphysical processes,as discussed above,it is necessary to also examine the dominant cloud microphysical processes involved in the torrential rainfall event.To address this,we used similar schemes and the same model as Li et al.(2014)to carry out a high-resolution numerical simulation of the orographic rainfall event.However,as an extension to the work of Li et al.(2014),in our study cloud microphysical data were output and analyzed to investigate the above aspect.The numerical simulation and cloud microphysical scheme are briefy described in section 2,the results are presented in section 3, and a summary is given in section 4.

    2. Numerical simulationandcloudmicrophysical scheme

    2.1.Numerical simulation

    A torrential rainfall event in Sichuan Province of China from 1800 UTC 17 August to 0000 UTC 20 August 2010 was simulated using the same WRF model setup as in Li et al.(2014).The simulation was integrated for 54 hours,using three nests with resolutions of 27,9,and 3 km.In the 3 km domain,only the cloud microphysical parameterization scheme was used[for simulation verifcation and other details,see Li et al.(2014)],and the cloud microphysical data (including cloud microphysical conversion terms)were output and analyzed in this study.

    2.2.Cloud microphysical scheme

    The cloud microphysical scheme used was the Milbrandt 2-mom scheme(Milbrandt and Yau,2005a,b),which includes water vapor and six cloud species:water vapor(Qv); cloud water(Qc);rain water(Qr);cloud ice(Qi);snow(Qs); graupel(Qg);and hail(Qh).Figure 1 shows the microphysical fowchart for the Milbrandt 2-mom scheme[for detailed microphysical conversion processes,see Milbrandt and Yau (2005b)].The tendencies for the mixing ratios are

    whereSQv,SQc,SQr,SQi,SQs,SQgandSQhare microphysical source/sink terms of water vapor and cloud species,respectively.The notation for the terms involving two interacting categories is denoted by QAByx(also see Milbrandt and Yau, 2005b),whereQis theprognosticvariablemass mixingratio, AB represents microphysical processes(QCL for collection, QCN for conversion,QFZ for freezing,QIM for ice multiplication,QML formelting,QNU fornucleation,QSH forshedding,QVD for diffusional growth),and the subscript“yx”indicates that mass is being transferred from categoryytox[x,y∈(v,c,r,i,s,g,h)denote various water species].δxryis a determination coeffcient meaning that categoryx(s,i,g) and raindrops collide to form categoryy(s,g,h).For example,if raindrops and snow collide to form graupel,thenδsrg=1,δsrs=0 andδsrh=0(for details,see Milbrandt and Yau,2005b).

    3. Results

    3.1.Precipitation and cloud hydrometeors distribution

    Figure 2 shows the distributions of 6-hourly accumulated rainfall and 6-hourlyaveragedand verticallyintegratedaIntegral calculation method:where,Qx,vintrepresents vertically integrated mixing ratio of category x(Qx),PtandPsarepressure at the top and the surface of the atmosphere.mixing ratios of hydrometeors in the weak(0000–0600 UTC 18 August 2010)and strong(1800UTC 18 August to 0000UTC 19 August 2010)precipitation periods.The mixing ratio of hail was negligibly small in this rainfall event and is not analyzed in the following discussions.In the strong precipitation period,the area and intensity of precipitation and cloud hydrometeors increased compared with those in the weak precipitation period(Fig.2).In the weak precipitationperiod,6-hourlyaccumulatedrainfall didnot exceed40 mm,compared to the maximum 6-hourly accumulated rainfall of over 100 mm in the strong precipitationperiod(Figs.2a and b).As for the 6-hourly averaged and vertically integrated mixing ratios ofhydrometeors,themaximumofthemixingratiosincreased from less than 0.09×10-3kg kg-1to more than 0.27×10-3kg kg-1for cloud water(Figs.2c and d),from less than 0.06×10-3kg kg-1to morethan0.3×10-3kgkg-1forrain water(Figs.2e and f),from less than 0.06×10-3kg kg-1to morethan 0.14×10-3kg kg-1forcloudice(Figs.2gand h), from less than 0.2×10-3kg kg-1to more than 0.7×10-3kg kg-1for snow(Figs.2i and 2j),and from less than 0.04× 10-3kg kg-1to more than 0.6×10-3kg kg-1for graupel (Figs.2k and l).In the strong precipitation period,the centers of the rain water mixing ratio coincided well with rainfall centers as well as cloud water,snow and graupel,indicating that the evolutions of both liquid and solid hydrometeors may have important contributions to precipitation development.Compared with the weak precipitation period,in the strong precipitation period,not only liquid hydrometeors(cloud water and rain water)increased signifcantly(themaximumgrowthof cloud water contentwas more than three times greater,and the maximum growth of rain water content was more than fve times greater),but also solid hydrometeors(cloudice,snowandgraupel)grewsignifcantly(themaximum growth of cloud ice content was more than two times greater,the maximum growth of snow was more than three times greater,and the maximum growth of graupel was more than ffteen times greater),and their centers were consistent with strong rainfall centers,showing that cold cloud processes developed remarkably and ice-phase processes contributed much more to the rainfall in the strong precipitation period.

    3.2.Cloud microphysical characteristics

    To further study the differences in cloud microphysical characteristics between the weak and strong precipitation periods,we focused our analysis on the region of(28°–33°N, 100°–106°E)in the 3 km domain.Grids whose 6-hourly accumulated rainfall was greater than 35 mm(256 grids in total)in the weak precipitation period(0000–0600 UTC 18August 2010)and grids whose 6-hourly accumulated rainfall was greater than 80 mm(228 grids in total)in the strong precipitation period(1800 UTC 18 August to 0000 UTC 19 August 2010)were chosen for analysis.

    Figure3 shows compositesof verticalproflesofhydrometeor mixing ratios(Figs.3a and b),number concentrations (Figs.3c and d),mean-mass diameters(Figs.3e and f),and terminal velocities and air vertical velocities(Figs.3g and h)in the weak and strong precipitation periods.In the Milbrandt2-momscheme,the mean-mass diameterDmxforeach hydrometeorcategoryxwas computed by

    whereρis the air density,cxanddxare parameters for the mass–diameter relations[see Milbrandt and Yau(2005a)for specifc values],Qxis the mixing ratio,andNxis the total number concentration.Tables 1 and 2 show the vertical accumulations of the above profles of mixing ratios and number concentrations,respectively.In the weak precipitation period,graupel was the most abundant(3.9036),the mixing ratios of cloud water(2.5644)and rain water(2.3247) were similar,and cloud ice(0.6282)and snow(0.7426)were the least prevalent(Fig.3a,Table 1).In the strong precipitation period,the abundance of the main hydrometeors increasedsignifcantly(Fig.3b).Theratiosofthemixingratios of cloud water,rain water,cloud ice,snow and graupel between the weak andstrong precipitationperiods were 1.6170,2.5864,5.1195,9.0687 and 1.7009,respectively(Table 1). A larger rain water mixing ratio(6.0126)corresponded with the strong precipitation.Higher snow(6.7344)and graupel (6.6396)mixing ratios may have made an important contribution to the development of the cloud system and strong precipitation;compared with the weak precipitation period, the higher cloud water mixing ratio(4.1467)in the strong precipitation period may also have played an important role in the strong surface precipitation(Fig.3b,Table 1).Furthermore,as shown in Figs.3a and b,and the vertically stratifed accumulationresults in Tables 1 and 2,the cloud water above and below the zero layer possessed similar magnitudes,and abundantsuper-cooledcloud water may havemade an important contribution to the heavy rainfall.Rain water(cloud ice, snow and graupel)was mainly distributes below(above)the zero layer,and the coalescence of cloud water and rain water below the zero layer should also plays an important role in strong rainfall.There was also some graupel distributed below the zero degree layer,whose melting may have made an important contribution to the growth of rain water and the enhancement of surface rainfall.

    As shown in Figs.3e and f,in the strong precipitation period,particle sizes of all hydrometeors increased compared with those in the weak precipitation period,and mean-mass diameters of graupel increased the most signifcantly(maximum increased from about 0.9 mm to more than 1.9 mm), while the increase in the mean-mass diameters of cloud ice was not signifcant.In Eq.(8),the mean-mass diameterDmxis not only correlated with the number concentrationNx,but also with the mixing ratioQx.Although the number concentrations of all hydrometeors increased(Figs.3c and d,Table 2),their mixingratios increased as well(Figs.3a and b,Table 1).Therefore,the ratio of the mixing ratioQxto the number concentrationNxforeachhydrometerincreased(butwithdistinctly different magnitudes),indicating that the mean-mass diameter increased,thereby increasing the mean particle size (Figs.3e and f).The increasing of graupel particle sizes may be partially related to the accretion of cloud water by graupel to form graupel,which causes the graupel mixing ratio to increase while the graupel number concentration remains roughly the same.Similarly,the increasing of rain water particle sizes may be partially related to the accretion of cloud water by rain water to form rain water.Besides,the increasing of particle sizes of all hydrometeors is favorable for the occurrence of collision between particles,especially cloud water with rain water and graupel.

    Meanwhile,Figs.3c and d show that the height of the number concentration maximum of cloud water and rain water in the strong precipitation period was higher than in the weak precipitation period.This may have been related to stronger vertical upward motion in the strong precipitation period(Figs.3g and h).By analyzing the vertical profles of vertical velocity of air and various hydrometeors’terminal velocities in the two periods(Figs.3g and h),we can see thatthe rain water’s terminal velocity was the strongest among all hydrometeors in both the strong(strongest value was close to-6.5 m s-1)and weak(strongest value was close to-4.5 m s-1)precipitationperiods.Graupel’sterminalvelocitywas second strongest,but its strongest value was only about-2 m s-1(in the strong precipitation period),far below that of rain water.In the strong precipitation period,the differences among the various hydrometeors’terminal velocities became larger compared to in the weak precipitation period,which would have been more benefcial to relative motion,coalescence and transformation between particles.It is worth noting that the absolute values of the rain water’s and graupel’s terminal velocities were signifcantly greater than that of air vertical velocity,and the stronger the precipitation was,the larger the differences between them were(Figs.3g and h). This was the result of the development of the cloud microphysicalprocessesduringtherainstorm,whichmayhavepresaged the gradual demise of the storm cloud system and the gradual end of the rainstorm process.

    Table 1.Vertical cumulative mixing ratios of cloud hydrometeors(Qv,Qc,Qr,Qi,Qs,Qg,Qh,units:10-3kg kg-1)throughout the whole layer,above and below the zero layer,in the weak precipitation period(0000–0600 UTC 18 Aug 2010)and the strong precipitation period (1800 UTC 18 Aug–0000 UTC 19 Aug 2010).The ratio is(1818–1900)/(1800–1806).

    Table 2.Vertical cumulative number concentrations of cloud hydrometeors(Nc,Nr,Ni,Ns,Ng,Nh,units:kg-1)throughout the whole layer, above and below the zero layer,in the weak precipitation period(0000–0600 UTC 18 Aug 2010)and the strong precipitation period(1800 UTC 18 Aug–0000 UTC 19 Aug 2010).The ratio is(1818–1900)/(1800–1806).

    3.3.Cloud microphysical conversion processes

    The microphysical sources/sinks and conversion processes of hydrometeors in the weak and strong precipitation periods are shown in Fig.4.The hydrometeors’main microphysical sources/sinks were similar in both periods,but the sources/sinks in the strong precipitation period were signifcantly stronger than those in the weak precipitation period,ultimately resulting in differences in the intensity of precipitation.Furthermore,the orders of magnitude of the sources/sinks(cloud microphysical conversion rates)of water vapor,cloud water,rain water and graupel were much larger than those of cloud ice and snow,indicating that water vapor,cloud water,rain water,graupel and their conversion processes played a major role during the development of the rainstorm.Condensation of water vapor to form cloud water(QVDvc)was the main source term of cloud water,while accretion of cloud water by rain water(QCLcr)and the accretion of cloud water by graupel(QCLcg)were the main sink terms of cloud water(Figs.4c and d).In the budget of graupel,the collection of cloud water by graupel(QCLcg)was a main source term of graupel,while the melting of graupel to form rain water(QMLgr)was a major consumption term of graupel(Figs.4k and l).The graupel production from the process of graupel–rain collision(QCLgrg)required the expenditure of graupel(QCLgr)and rain water(QCLrg).Therefore,the net growth of graupel by graupel–rain collision and coalescence was(QCLgrg-QCLgr).Meanwhile,since absolute values of QCLgrgand QCLgrwere similar(Figs.4k and l),values of(QCLgrg-QCLgr)were small,meaning that the collection of rain water by graupel to form graupel was not a main source/sink term of graupel(Figs.4k and l).As for rain water,the accretionof cloudwater byrain water(QCLcr) and the melting of graupel to form rain water(QMLgr)were the main source terms(Figs.4e and f),which caused the rapid growth of rain water(Fig.4f)and eventually enhanced surface precipitation in the strong precipitation period(Fig. 2b).

    The above analysis of hydrometeors’microphysical sources/sinks is summarizedin Fig.5.The solid thick arrows (Fig.5)represent the most important microphysical conversion processes in this rainstorm.The main cloud microphysical conversion processes related to precipitation in the two periods were the same,shown as solid thick arrows in Fig.5, and proportions of these processes to total water vapor loss (WVLbWVL,Water Vapor Loss is defned asWVL=max(QVDvc,0)+max(QVDvr,0)+max(QNUvi,0)+max(QVDvi,0)+max(QVDvs,0)+max(QVDvg,0)+max(QVDvh,0),Colle and Zeng,2004)in each period were roughly equal(see the values on the right-hand sides of the conversion process notations in Fig.5).However,due to different total water vapor loss rates in the strong(WVL=15.73)and weak(WVL=7.38)precipitation periods,the microphysical conversion rates between these two periods were signifcantlydifferent(seevaluesin parenthesesin Fig.5),resulting in signifcant differences in precipitation intensity between these two periods.In the strong precipitation period(Fig. 5b),there were two main pathways to generate rain water: QVDvc→QCLcrand QVDvc→QCLcg→QMLgr,which can be described as follows:Abundant water vapor condensedintocloudwater(QVDvc)and,ontheonehand,accretion of cloud water by rain water(QCLcr)formed rain water, while on the other hand,accretion of cloud water by graupel(QCLcg)formed graupel,and then the melting of graupel formed rain water(QMLgr).As discussed above,the net growthof graupelby graupel–rainwater collision and coalescence(QCLgrg-QCLgr)was not a major source of graupel and was also not a major sink of rain water(Fig.5).Therefore,the net result of(QCLgrg-QCLgr)contributed little to raindrop growth and enhancement of surface precipitation. In Fig.5,we can see another distinct difference between the two periods:In the strong precipitation period,depositional growth of snow(QVDvs),melting of snow to form rain water(QMLsr),and autoconversion(aggregation)of snow to form graupel(QCNsg)strengthened,which played a role in the growth of rain water(Fig.5b).Due to the accumulation of abundant moisture in the strong precipitation period, water vapor deposition to form ice-phase particles,such as snow(QVDvs),was more likely to occur in the upper levels. Meanwhile,in the strong precipitation period,mixing ratios of cloud water and snow were signifcantly larger than those in the weak precipitation period(Figs.3a and b),so the accretion of cloud water by snow(QCLcs)strengthened(Fig. 5).The increasing of snow also had an active effect on the autoconversion of snow to form graupel.Besides,from Fig. 3b,we can see that the vertical profle of the snow mixing ratio extended below the zero layer,thereby enhancing the process of snow melting into rain water.However,the orders of magnitudes of these conversion processes associated with snow were signifcantly smaller than those of the mainmicrophysical conversion processes mentioned above.

    4. Summary

    In this study,by using high-resolution numerical simulation data of an orographic torrential rainfall process in Sichuan simulated by the WRF model,differences of cloud microphysical characteristics and conversion processes between weak and strong precipitation periods were analyzed to investigatethedominantcloudmicrophysicalprocesses responsible for the development of the torrential rainfall.The major results can be summarized as follows:

    (1)In the strong precipitation period,the area and intensity of precipitation and cloud hydrometeors increased compared with those in the weak precipitation period,and the centers of rain water,cloud water and graupel were collo-cated with strong rainfall centers,showing that both liquid and solid hydrometeors may play an important role in surface precipitation.

    (2)Composite analysis of cloud hydrometeors revealed that,in the strong precipitation period,the cloud hydrometeors increased signifcantly.Abundantsuper-cooledcloud water may have made an important contribution to heavy rainfall.The coalescence of cloud water and rain water below the zero layer and the melting of graupel may also have played important roles in the growth of rain water and the enhancement of surface rainfall.In the strong precipitation period, particle sizes of all hydrometeors increased,especially graupel’s.Theincreasingofparticlesizes ofall hydrometeorswas favorable to the occurrenceof collision between particles,especially cloud water with rain water and graupel.

    (3)The terminal velocity of raindrops was the strongest among all hydrometeors in both the weak and strong precipitation periods,followed by that of graupel.Differences of terminal velocities among the various hydrometeors in the strong precipitation periodwere larger than those in the weak precipitationperiod,which was favorable for the relative motion,collection interaction and exchange between particles. Absolute values of terminal velocities of raindrops and graupel were signifcantly greater than those of air updrafts.Thestronger the precipitation was,the greater the differences between them were,which may imply a decay of the rainfall.

    (4)In both periods,the budgets of cloud microphysics weresimilar,buttheordersofmagnitudesofthecloudmicrophysical conversion terms in the strong precipitation period were signifcantly larger than those in the weak precipitation period,resulting in a difference in the intensity of surface rainfall.Besides,the orders of magnitude of the conversion rates of water vapor,cloudwater,rain water and graupelwere much larger than those of cloud ice and snow,showing that water vapor,cloud water,rain water and graupel and their conversionprocesses may have played a major role in the development of the torrential rainfall event.

    (5)In the rainstorm process,there were two main pathways for the generation of raindrops:QVDvc→QCLcrand QVDvc→QCLcg→QMLgr.In detail,abundant water vapor condensed into cloud water(QVDvc)and,on the one hand, the accretion of cloud water by rain water(QCLcr)formed rainwater,whileontheotherhand,theaccretionofcloudwater by graupel(QCLcg)formed graupel,and then the melting of graupel formed rain water(QMLgr).Furthermore,the net growthofgraupelby graupel–rainwater collision andcoalescence(QCLgrg-QCLgr)contributedlittle to raindropgrowth and enhancement of surface precipitation.

    Our own previous analyses and studies have shown that the type of orographic torrential rainfall discussed in this study is very common in summer in the Sichuan region. Thus,although only one case has been analyzed in this paper,the conclusions may have a certain degree of universality,which is of great signifcance as we strive for an in-depth and comprehensive understanding of the mechanisms underpinning orographic torrential rainfall in the Sichuan region. In future work,we will analyze more cases to make the conclusions more robust and statistically signifcant.

    Acknowledgements.This work was supported by the Key Research Program of the Chinese Academy of Sciences(Grant No. KZZD-EW-05-01)and the National Basic Research Program of China(973 Program)(Grant No.2014CB441402).

    REFERENCES

    Chen,C.P.,H.Z.Feng,and J.Chen,2010:Application of Sichuan heavy rainfall ensemble prediction probability products based on Bayesian method.Meteorological Monthly,36(5),32–39. (in Chinese)

    Chen,Y.R.,and Y.Q.Li,2013:Characteristics of mesoscale convective system and its effects on short-time severe rainfall in Sichuan Basin during 21–22 July 2012.Meteorological Monthly,39(7),848–860.(in Chinese)

    Colle,B.A.,and Y.G.Zeng,2004:Bulk microphysical sensitivities within the MM5 for orographic precipitation.Part II: Impact of barrier width and freezing level.Mon.Wea.Rev., 132,2802–2815.

    Crochet,P.,T.J′ohannesson,T.J′onsson,O.Siguresson,H. Bj¨ornsson,F.P′alsson,and I.Barstad,2007:Estimating the spatial distribution of precipitation in Iceland using a linear model of orographic precipitation.J.Hydrometeor.,8,1285–1306.

    Ge,J.J.,W.Zhong,N.Du,and H.C.Lu,2008:Numerical simulationandanalysis of Sichuanrainstormunder terraininfuence.Scientia Meteorologica Sinica,28(2),176–183.(in Chinese)

    He,Y.,and G.P.Li,2013:Numerical experiments on infuence of Tibetan Plateau on persistent heavy rain in South China.Chinese J.Atmos.Sci.,37(4),933–944.(in Chinese)

    Li,A.T.,and H.R.He,2011:Impact of different cloud microphys-ical parameterization schemes on the numeric simulation results of“8.8”rainstorm process in Zhouqu.Meteorology and Disaster Reduction Research,34(3),9–16.(in Chinese)

    Li,Q.,X.P.Cui,and J.Cao,2014:The observational analysis and numerical simulation of a heavy rain event in Sichuan Province.Chinese J.Atmos.Sci.,doi:10.3878/j.issn.1006-9895.1401.13255.(in Chinese)

    Liao,J.B.,X.M.Wang,B.C.Xia,T.J.Wang,and Z.M.Wang, 2012:The effects of different physics and cumulus parameterization schemes in WRF on heavy rainfall simulation in PRD.J.Trop.Meteor.,28(4),461–470.(in Chinese)

    Lu,P.,R.C.Yu,and T.J.Zhou,2009:Numerical simulation on the sensitivity of heavy rainfall over the western Sichuan Basin to initial water vapor condition.Chinese J.Atmos.Sci.,33(2), 241–250.(in Chinese)

    Ma,Y.Z.,C.G.Lu,and S.T.Gao,2012:The effects of different microphysical schemes in WRF on a heavy rainfall in North China during 18–19 August 2010.Chinese J.Atmos. Sci.,36(4),835–850.(in Chinese)

    Maussion,F.,D.Scherer,R.Finkelnburg,J.Richters,W.Yang, and T.D.Yao,2010:WRF simulation of a precipitation event over the Tibetan Plateau,China—An assessment using remote sensing and ground observations.Hydrology and Earth System Sciences Discussions,15,1795–1817.

    Milbrandt,J.A.,and M.K.Yau,2005a:A multimoment bulk microphysics parameterization.Part I:Analysis of the role of the spectral shape parameter.J.Atmos.Sci.,62,3051–3064.

    Milbrandt,J.A.,and M.K.Yau,2005b:A multimoment bulk microphysics parameterization.Part II:A proposed threemoment closure and scheme description.J.Atmos.Sci.,62, 3065–3081.

    Sawyer,J.S.,1956:The physical and dynamical problems of orographic rain.Weather,11(12),375–381.

    Smith,B.L.,S.E.Yuter,P.J.Neiman,and D.E.Kingsmill,2010: Water vapor fuxes and orographic precipitation over Northern California associated with a landfalling atmospheric river.Mon.Wea.Rev.,138,74–100.

    Smith,R.B.,and I.Barstad,2004:A linear theory of orographic precipitation.J.Atmos.Sci.,61,1377–1391.

    Smith,R.B.,and J.P.Evans,2007:Orographic precipitation and water vapor fractionation over the Southern Andes.J.Hydrometeor.,8,3–19.

    Wang,C.X.,S.T.Gao,L.Liang,and Y.Z.Ma,2013:Diagnostic analysis of dynamical parameters for Sichuan rainstorm infuenced by terrain.Chinese J.Atmos.Sci.,37(5),1099–1110. (in Chinese)

    Wang,W.,J.Liu,and X.J.Cai,2011:Impact of mesoscale gravity waves on a heavy rainfall event in the east side of the Tibetan Plateau.Transactions of Atmospheric Sciences,34(6),737–747.(in Chinese)

    :Huang,Y.J.,and X.P.Cui,2015:Dominant cloud microphysical processes of a torrential rainfall event in Sichuan, China.Adv.Atmos.Sci.,32(3),389–400,

    10.1007/s00376-014-4066-7.

    (Received 28 March 2014;revised 5 June 2014;accepted 10 July 2014)

    ?Corresponding author:CUI Xiaopeng

    Email:xpcui@mail.iap.ac.cn

    国模一区二区三区四区视频| 国产一区二区三区av在线| 成人综合一区亚洲| 日韩亚洲欧美综合| 精品午夜福利在线看| 身体一侧抽搐| 欧美97在线视频| 免费黄色在线免费观看| av一本久久久久| 深夜a级毛片| 亚洲经典国产精华液单| av在线老鸭窝| 51国产日韩欧美| 久久ye,这里只有精品| 免费观看a级毛片全部| 久久鲁丝午夜福利片| 亚洲精品一区蜜桃| 男人和女人高潮做爰伦理| 精品久久国产蜜桃| 高清日韩中文字幕在线| 欧美激情久久久久久爽电影| 国产日韩欧美在线精品| 精品亚洲乱码少妇综合久久| 少妇熟女欧美另类| 亚洲四区av| 好男人在线观看高清免费视频| 成人黄色视频免费在线看| 国产成人a区在线观看| 久久久色成人| 午夜老司机福利剧场| av一本久久久久| 在线看a的网站| 另类亚洲欧美激情| 国产爽快片一区二区三区| 少妇熟女欧美另类| 国语对白做爰xxxⅹ性视频网站| 日韩不卡一区二区三区视频在线| 大陆偷拍与自拍| av在线亚洲专区| 在线观看一区二区三区激情| 97精品久久久久久久久久精品| 国精品久久久久久国模美| 欧美日韩视频精品一区| 亚洲经典国产精华液单| 18+在线观看网站| 亚洲精品第二区| 69av精品久久久久久| 欧美激情在线99| 欧美性感艳星| 精品久久久久久久久亚洲| www.av在线官网国产| 亚洲精品一区蜜桃| 亚洲国产精品国产精品| 水蜜桃什么品种好| 性插视频无遮挡在线免费观看| 一边亲一边摸免费视频| 亚洲天堂av无毛| 两个人的视频大全免费| 亚洲最大成人中文| 精品一区二区三卡| 中文字幕制服av| 欧美成人一区二区免费高清观看| 亚洲伊人久久精品综合| 成人毛片a级毛片在线播放| 久久久精品免费免费高清| 亚洲av不卡在线观看| 乱系列少妇在线播放| 一二三四中文在线观看免费高清| 免费电影在线观看免费观看| 国产极品天堂在线| 免费观看av网站的网址| 国产伦在线观看视频一区| 干丝袜人妻中文字幕| 欧美成人a在线观看| a级毛片免费高清观看在线播放| 夫妻午夜视频| 久久人人爽av亚洲精品天堂 | 波多野结衣巨乳人妻| 欧美成人午夜免费资源| 亚洲国产色片| 久久97久久精品| 国产精品三级大全| 亚洲美女视频黄频| 国产精品国产三级专区第一集| 中文字幕av成人在线电影| 丝瓜视频免费看黄片| 在线免费观看不下载黄p国产| 男人爽女人下面视频在线观看| 国产精品99久久99久久久不卡 | 一二三四中文在线观看免费高清| 久久人人爽人人爽人人片va| av又黄又爽大尺度在线免费看| 亚洲怡红院男人天堂| 欧美xxⅹ黑人| 在线a可以看的网站| 精品99又大又爽又粗少妇毛片| 久久久久国产精品人妻一区二区| 亚洲av免费高清在线观看| 狂野欧美激情性xxxx在线观看| 亚洲精品乱久久久久久| 婷婷色综合大香蕉| 亚洲精品乱码久久久久久按摩| 中文资源天堂在线| 男女无遮挡免费网站观看| 久久久欧美国产精品| 久久久久精品性色| 久久午夜福利片| 美女xxoo啪啪120秒动态图| 夫妻午夜视频| 夫妻午夜视频| 狂野欧美激情性xxxx在线观看| 国产精品久久久久久精品电影| 人人妻人人澡人人爽人人夜夜| 在线观看一区二区三区激情| av线在线观看网站| 在线观看人妻少妇| 熟女人妻精品中文字幕| 亚洲精品久久午夜乱码| 在线观看国产h片| 欧美另类一区| 日韩三级伦理在线观看| 国产精品一及| 日日撸夜夜添| 国产精品久久久久久精品电影| 午夜日本视频在线| 成人毛片a级毛片在线播放| 99久久九九国产精品国产免费| 免费av观看视频| 青春草视频在线免费观看| 亚洲av在线观看美女高潮| 免费看日本二区| 乱系列少妇在线播放| 亚洲av免费高清在线观看| 女人十人毛片免费观看3o分钟| 亚洲最大成人av| 少妇的逼好多水| xxx大片免费视频| 精品久久国产蜜桃| 国产一区二区亚洲精品在线观看| av国产免费在线观看| 久久精品国产自在天天线| 日韩制服骚丝袜av| 国产精品秋霞免费鲁丝片| 大片电影免费在线观看免费| 久久精品人妻少妇| 大陆偷拍与自拍| h日本视频在线播放| 欧美日韩综合久久久久久| 国内精品美女久久久久久| 亚洲欧美精品自产自拍| 国产成人午夜福利电影在线观看| 80岁老熟妇乱子伦牲交| 免费黄频网站在线观看国产| 成人美女网站在线观看视频| av专区在线播放| 欧美变态另类bdsm刘玥| 久久久a久久爽久久v久久| 国产日韩欧美亚洲二区| 精品久久久久久久人妻蜜臀av| 成人欧美大片| 大话2 男鬼变身卡| 80岁老熟妇乱子伦牲交| 不卡视频在线观看欧美| 国产伦理片在线播放av一区| 插逼视频在线观看| 日本免费在线观看一区| 波野结衣二区三区在线| 日本免费在线观看一区| 中文字幕av成人在线电影| 国产一区二区三区av在线| 99热国产这里只有精品6| 成人国产av品久久久| 99热这里只有是精品在线观看| 春色校园在线视频观看| 国产亚洲av嫩草精品影院| 男人和女人高潮做爰伦理| 国产淫片久久久久久久久| 狂野欧美激情性bbbbbb| 简卡轻食公司| 熟女人妻精品中文字幕| 男人和女人高潮做爰伦理| 简卡轻食公司| 在线免费观看不下载黄p国产| 日韩一区二区三区影片| 国产爽快片一区二区三区| 免费av观看视频| 国产精品国产三级专区第一集| 91久久精品国产一区二区三区| 精品人妻熟女av久视频| 国产色爽女视频免费观看| 国产免费一级a男人的天堂| 男女边吃奶边做爰视频| 亚洲精品成人av观看孕妇| 波多野结衣巨乳人妻| 久久久午夜欧美精品| 一级毛片电影观看| 在线免费观看不下载黄p国产| 内射极品少妇av片p| av国产免费在线观看| 嫩草影院入口| 中文字幕人妻熟人妻熟丝袜美| 在线看a的网站| 夫妻午夜视频| 人妻系列 视频| 免费看不卡的av| 欧美三级亚洲精品| 成人黄色视频免费在线看| 日韩av在线免费看完整版不卡| 边亲边吃奶的免费视频| 国产爱豆传媒在线观看| 天天躁夜夜躁狠狠久久av| 男女边摸边吃奶| 国产高潮美女av| 国产乱人视频| 边亲边吃奶的免费视频| 1000部很黄的大片| 22中文网久久字幕| 亚洲精品视频女| 免费看a级黄色片| 自拍偷自拍亚洲精品老妇| 午夜福利在线在线| 男女边摸边吃奶| 只有这里有精品99| 国产真实伦视频高清在线观看| 成人黄色视频免费在线看| 亚洲av男天堂| 超碰av人人做人人爽久久| 嫩草影院新地址| 成年人午夜在线观看视频| 欧美3d第一页| 别揉我奶头 嗯啊视频| 97超碰精品成人国产| av福利片在线观看| 神马国产精品三级电影在线观看| av免费在线看不卡| 精品一区在线观看国产| 亚洲av一区综合| 精品少妇久久久久久888优播| 日韩免费高清中文字幕av| 色视频www国产| 极品少妇高潮喷水抽搐| 中文资源天堂在线| www.av在线官网国产| 国产免费福利视频在线观看| 日日啪夜夜撸| 草草在线视频免费看| 中国国产av一级| 男人和女人高潮做爰伦理| 久久99热6这里只有精品| 亚洲精品成人av观看孕妇| 99热这里只有精品一区| 一级毛片aaaaaa免费看小| 不卡视频在线观看欧美| 日本熟妇午夜| 岛国毛片在线播放| 国产精品福利在线免费观看| 有码 亚洲区| 精品久久久久久久末码| 日韩 亚洲 欧美在线| 免费黄网站久久成人精品| 国产精品国产av在线观看| 十八禁网站网址无遮挡 | 亚洲精品一区蜜桃| 伦精品一区二区三区| 国产乱人视频| 少妇丰满av| 欧美成人一区二区免费高清观看| 久久久久精品性色| 女的被弄到高潮叫床怎么办| 丝瓜视频免费看黄片| 伊人久久精品亚洲午夜| 亚洲精品视频女| 国产国拍精品亚洲av在线观看| 99九九线精品视频在线观看视频| 丰满少妇做爰视频| 一级黄片播放器| av在线app专区| 尾随美女入室| 国产有黄有色有爽视频| 黄片无遮挡物在线观看| 五月玫瑰六月丁香| 69av精品久久久久久| 亚洲精品第二区| 老师上课跳d突然被开到最大视频| 亚洲性久久影院| 亚洲av日韩在线播放| 2022亚洲国产成人精品| 国产亚洲91精品色在线| 亚洲国产最新在线播放| 韩国高清视频一区二区三区| 又黄又爽又刺激的免费视频.| 国产精品99久久久久久久久| 欧美三级亚洲精品| 日韩中字成人| 亚洲欧洲日产国产| 国产精品久久久久久av不卡| 少妇人妻精品综合一区二区| 国产91av在线免费观看| 亚洲人成网站在线观看播放| 美女脱内裤让男人舔精品视频| 亚洲成人精品中文字幕电影| 国产国拍精品亚洲av在线观看| 卡戴珊不雅视频在线播放| 亚洲精品成人av观看孕妇| 国产在线男女| 国产成人精品久久久久久| 亚洲国产精品专区欧美| 欧美成人a在线观看| 亚洲精品中文字幕在线视频 | 看十八女毛片水多多多| 亚洲欧美日韩无卡精品| 国产黄a三级三级三级人| 国产在视频线精品| 狠狠精品人妻久久久久久综合| 97超碰精品成人国产| 91午夜精品亚洲一区二区三区| av又黄又爽大尺度在线免费看| 综合色丁香网| a级毛色黄片| 天堂网av新在线| 九九在线视频观看精品| 少妇人妻 视频| 嫩草影院入口| 狠狠精品人妻久久久久久综合| 国产精品一区二区三区四区免费观看| 亚洲精品乱码久久久久久按摩| 日本wwww免费看| 色播亚洲综合网| 在线亚洲精品国产二区图片欧美 | 国产 一区精品| 国产精品秋霞免费鲁丝片| 免费观看性生交大片5| 狂野欧美白嫩少妇大欣赏| 人人妻人人看人人澡| 成人毛片a级毛片在线播放| 七月丁香在线播放| 国产免费一区二区三区四区乱码| 欧美成人午夜免费资源| 人人妻人人澡人人爽人人夜夜| 菩萨蛮人人尽说江南好唐韦庄| 亚洲婷婷狠狠爱综合网| 国产乱人偷精品视频| 亚洲性久久影院| 一区二区三区精品91| 亚洲最大成人av| 色视频www国产| 女人久久www免费人成看片| 制服丝袜香蕉在线| 成人国产av品久久久| 麻豆国产97在线/欧美| 在现免费观看毛片| 久久这里有精品视频免费| 国产免费一级a男人的天堂| 国产精品一及| 王馨瑶露胸无遮挡在线观看| 中文字幕免费在线视频6| 国产探花极品一区二区| 2021天堂中文幕一二区在线观| 日本猛色少妇xxxxx猛交久久| 18禁裸乳无遮挡免费网站照片| 内地一区二区视频在线| 免费av不卡在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产黄片视频在线免费观看| 国产精品久久久久久久久免| 少妇高潮的动态图| 中国三级夫妇交换| 简卡轻食公司| kizo精华| 午夜日本视频在线| 色吧在线观看| 亚洲,一卡二卡三卡| 一级av片app| 国产毛片在线视频| eeuss影院久久| 在线免费十八禁| 青春草视频在线免费观看| 日本色播在线视频| 嫩草影院入口| 日日啪夜夜爽| 在线观看av片永久免费下载| 18禁在线播放成人免费| 久久99热6这里只有精品| 18+在线观看网站| 免费看a级黄色片| 国产综合精华液| 国产久久久一区二区三区| 人体艺术视频欧美日本| 高清av免费在线| 人妻系列 视频| 十八禁网站网址无遮挡 | 日日摸夜夜添夜夜爱| 尤物成人国产欧美一区二区三区| 欧美三级亚洲精品| 欧美少妇被猛烈插入视频| a级毛色黄片| 免费大片黄手机在线观看| 亚洲伊人久久精品综合| 好男人视频免费观看在线| 又黄又爽又刺激的免费视频.| 国产女主播在线喷水免费视频网站| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 亚洲色图av天堂| 国产成人免费观看mmmm| 国产片特级美女逼逼视频| 国产精品人妻久久久久久| 久久精品国产自在天天线| 欧美成人a在线观看| 天堂网av新在线| 日韩av在线免费看完整版不卡| 青春草亚洲视频在线观看| 亚洲不卡免费看| 大话2 男鬼变身卡| 偷拍熟女少妇极品色| 久久久久久久久久久丰满| 成人亚洲精品一区在线观看 | 人妻 亚洲 视频| 国产精品成人在线| 69人妻影院| av在线app专区| 中文天堂在线官网| 视频中文字幕在线观看| 九九爱精品视频在线观看| 国产女主播在线喷水免费视频网站| 久久影院123| 欧美成人午夜免费资源| 两个人的视频大全免费| 偷拍熟女少妇极品色| 亚洲国产成人一精品久久久| 三级男女做爰猛烈吃奶摸视频| 亚洲av免费在线观看| 国产黄色视频一区二区在线观看| 午夜视频国产福利| 欧美日韩视频精品一区| 美女内射精品一级片tv| 亚洲av成人精品一区久久| 中文天堂在线官网| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影 | 国产真实伦视频高清在线观看| 午夜福利视频精品| av在线老鸭窝| 亚洲精品乱码久久久久久按摩| 久久影院123| 亚洲久久久久久中文字幕| 日本黄大片高清| 亚洲色图综合在线观看| 男女那种视频在线观看| 亚洲av欧美aⅴ国产| 亚洲aⅴ乱码一区二区在线播放| 国产一区有黄有色的免费视频| 亚洲真实伦在线观看| 亚洲高清免费不卡视频| 少妇的逼水好多| 日韩制服骚丝袜av| 在线免费观看不下载黄p国产| 白带黄色成豆腐渣| 黑人高潮一二区| 狂野欧美白嫩少妇大欣赏| 日韩国内少妇激情av| 精品一区二区免费观看| 久久99热6这里只有精品| 全区人妻精品视频| 校园人妻丝袜中文字幕| 国产一区有黄有色的免费视频| 亚洲欧美日韩无卡精品| 美女cb高潮喷水在线观看| 小蜜桃在线观看免费完整版高清| 国产午夜精品久久久久久一区二区三区| 男插女下体视频免费在线播放| 亚洲性久久影院| 亚洲自偷自拍三级| 中文字幕制服av| 亚洲精品视频女| 国产免费一区二区三区四区乱码| 人人妻人人看人人澡| 美女国产视频在线观看| 久久久久性生活片| 国产男女超爽视频在线观看| av一本久久久久| 伊人久久国产一区二区| 男人舔奶头视频| 国产亚洲av嫩草精品影院| 乱系列少妇在线播放| 高清毛片免费看| 亚洲激情五月婷婷啪啪| 久久午夜福利片| 亚洲精华国产精华液的使用体验| 日韩制服骚丝袜av| 亚洲精品第二区| 亚洲成色77777| 久久影院123| videos熟女内射| 看免费成人av毛片| 777米奇影视久久| 精品99又大又爽又粗少妇毛片| 天天躁夜夜躁狠狠久久av| 亚洲婷婷狠狠爱综合网| 亚洲无线观看免费| 麻豆国产97在线/欧美| 极品少妇高潮喷水抽搐| 久久久久精品久久久久真实原创| 国产精品国产三级国产专区5o| 午夜日本视频在线| 国产人妻一区二区三区在| av福利片在线观看| 国产黄片视频在线免费观看| 另类亚洲欧美激情| 黑人高潮一二区| 亚洲国产成人一精品久久久| 91久久精品电影网| 精品久久久久久久末码| 中文资源天堂在线| 日本wwww免费看| 人人妻人人爽人人添夜夜欢视频 | 一级毛片久久久久久久久女| 欧美成人精品欧美一级黄| 成人毛片60女人毛片免费| 久久6这里有精品| 国产成人一区二区在线| 麻豆国产97在线/欧美| 国产精品麻豆人妻色哟哟久久| 精品午夜福利在线看| 国产淫语在线视频| 丰满乱子伦码专区| 久久亚洲国产成人精品v| av在线天堂中文字幕| kizo精华| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲精品久久久com| 精品国产一区二区三区久久久樱花 | 免费黄色在线免费观看| 亚洲av日韩在线播放| 性色avwww在线观看| 少妇高潮的动态图| 偷拍熟女少妇极品色| 久热这里只有精品99| 亚洲国产欧美在线一区| 国产色爽女视频免费观看| 国产精品av视频在线免费观看| 狠狠精品人妻久久久久久综合| 搡女人真爽免费视频火全软件| 内地一区二区视频在线| 草草在线视频免费看| 99热这里只有精品一区| 日本与韩国留学比较| 免费在线观看成人毛片| 国产午夜精品久久久久久一区二区三区| 春色校园在线视频观看| 久久久久久久久久人人人人人人| 国产成人免费无遮挡视频| 亚洲在久久综合| 视频中文字幕在线观看| 麻豆成人av视频| 久久精品综合一区二区三区| 久久久久久久亚洲中文字幕| 日韩制服骚丝袜av| 亚洲图色成人| 国产精品秋霞免费鲁丝片| 亚洲久久久久久中文字幕| 婷婷色综合大香蕉| 亚洲精品中文字幕在线视频 | 少妇人妻久久综合中文| 成人亚洲精品av一区二区| 亚洲天堂av无毛| 女人被狂操c到高潮| 一区二区三区精品91| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 午夜爱爱视频在线播放| 我的老师免费观看完整版| 99久久精品国产国产毛片| 高清午夜精品一区二区三区| 能在线免费看毛片的网站| 亚洲成人av在线免费| 日韩欧美精品免费久久| 久久久久久久久久人人人人人人| 亚洲国产最新在线播放| 国产精品一区二区在线观看99| 天天躁日日操中文字幕| 亚洲成色77777| 男男h啪啪无遮挡| 成年版毛片免费区| 欧美zozozo另类| 国产成人精品福利久久| 一区二区三区四区激情视频| 国产成人一区二区在线| 深爱激情五月婷婷| 熟女av电影| av在线天堂中文字幕| 国产伦理片在线播放av一区| 青春草亚洲视频在线观看| 日韩成人伦理影院| 成人漫画全彩无遮挡| 亚洲无线观看免费| 亚洲av一区综合| 亚洲在久久综合| 嫩草影院新地址| 成人高潮视频无遮挡免费网站| 日本爱情动作片www.在线观看| 在线精品无人区一区二区三 | 国产一区亚洲一区在线观看| 尤物成人国产欧美一区二区三区| 婷婷色麻豆天堂久久| 毛片女人毛片| 18禁裸乳无遮挡免费网站照片| 老司机影院毛片| 精品久久久精品久久久| 日本免费在线观看一区| 日本欧美国产在线视频| 亚洲电影在线观看av| 一本一本综合久久| 边亲边吃奶的免费视频|