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

    Modification of Cumulus Convection and Planetary Boundary Layer Schemes in the GRAPES Global Model

    2015-11-21 11:23:29LIUKun劉琨CHENQiying陳起英andSUNJian孫健
    Journal of Meteorological Research 2015年5期
    關(guān)鍵詞:劉琨孫健

    LIU Kun(劉琨),CHEN Qiying(陳起英),and SUN Jian(孫?。?/p>

    1 National Meteorological Center,China Meteorological Administration,Beijing 100081

    2 Numerical Weather Prediction Center of China Meteorological Administration,Beijing 100081

    Modification of Cumulus Convection and Planetary Boundary Layer Schemes in the GRAPES Global Model

    LIU Kun1,2?(劉琨),CHEN Qiying1,2(陳起英),and SUN Jian1,2(孫?。?/p>

    1 National Meteorological Center,China Meteorological Administration,Beijing 100081

    2 Numerical Weather Prediction Center of China Meteorological Administration,Beijing 100081

    Cumulus convection is a key linkage between hydrological cycle and large-scale atmospheric circulation. Cumulus parameterization scheme is an important component in numerical weather and climate modeling studies.In the Global/Regional Assimilation and Prediction Enhanced System(GRAPES),turbulent mixing and diffusion approach is applied in its shallow convection scheme.This method overestimates the vertical transport of heat and moisture fluxes but underestimates cloud water mixing ratio over the region of stratocumulus clouds.As a result,the simulated low stratocumulus clouds are less than observations. To overcome this problem,a mass flux method is employed in the shallow convection scheme to replace the original one.Meanwhile,the deep convection scheme is adjusted correspondingly.This modification is similar to that in the US NCEP Global Forecast System(GFS),which uses the simplified Arakawa Schubert Scheme(SAS).The planetary boundary layer scheme(PBL)is also revised by considering the coupling between the PBL and stratocumulus clouds.With the modification of both the cumulus and PBL schemes,the GRAPES simulation of shallow convective heating rate becomes more reasonable;total amounts of stratocumulus clouds simulated over the eastern Pacific and their vertical structure are more consistent with observations;the underestimation of stratocumulus clouds simulated by original schemes is less severe with the revised schemes.Precipitation distribution in the tropics becomes more reasonable and spurious precipitation is effectively suppressed.The westward extension and northward movement of the western Pacific subtropical high simulated with the revised schemes are more consistent with Final Operational Global Analysis(FNL)than that simulated with the original schemes.The statistical scores for the global GRAPES forecast are generally improved with the revised schemes,especially for the simulation of geopotential height in the Northern Hemisphere and winds in the tropics.Root mean square errors(RMSEs)decrease in the lower and upper troposphere with the revised schemes.The above results indicate that with the revised cumulus and PBL schemes,model biases in the tropics decrease and the global GRAPES performance is greatly improved.

    GRAPES model,cumulus parameterization,boundary layer diffusion,mass flux

    1.Introduction

    Cumulus convection is an important process that links hydrologic cycle and large-scale atmospheric circulation.Cumulus parameterization scheme is one of the most important physical schemes in numerical weather and climate models.Particularly in the tropics where cumulus convections are active,large amounts of precipitation and accompanied latent heat release play a key role in tropical atmospheric circulation and global energy balance.Large latent heat release in the tropics also exerts great impact on midhigh latitude wind fields by providing kinetic energy to these regions.Cumulus parameterization scheme is a key component and also a research focus in numerical weather and climate modeling studies.Appropriate cumulus scheme can well describe impacts of several kilometer sub-grid scale processes on large-scale cir-culation,and thus is a key issue in numerical model development.Cumulus scheme includes two types:deep convection scheme and shallow convection scheme,according to the depth of convective clouds.Shallow convection helps for the vertical diffusion of moisture and heat fluxes,while the latent heat released during shallow convection can affect the development of deep convection via re-condensation of atmospheric moisture.This process subsequently affects large-scale circulation over the entire tropics and extratropics,making great contributions to the maintenance of the thermodynamic structure of the whole atmosphere(von Salzen et al.,2005;Nie and Kuang,2012).In addition,the cloud base of shallow convection is often located at the top of the planetary boundary layer(PBL).Therefore,successful simulation of stratocumulus clouds,which is related to the PBL structure,is highly dependent on appropriate shallow convection scheme(Li and Wang,2010;Liu et al.,2010;Yang et al.,2010).

    Supported by the National Natural Science Foundation of China(41305090),National Science and Technology Support Program of China(2012BAC22B02),and China Meteorological Administration Special Public Welfare Research Fund(GYHY201406005).

    ?Corresponding author:kliu@cams.cma.gov.cn.

    ?The Chinese Meteorological Society and Springer-Verlag Berlin Heidelberg 2015

    The simplified Arakawa-Schubert(SAS)cumulus parameterization scheme(Pan and Wu,1995)is currently applied in the global GRAPES.The deep convection in this scheme is a bulk convection scheme,which is based on the concept of mass flux proposed by Arakawa and Schubert(1974)but with some modification.This scheme can well describe essential physical interaction between cumulus convection and largescale circulation,especially the circulation in the tropics(Deng et al.,2008;Xu et al.,2014).

    However,shallow convection scheme in the SAS still follows the method proposed by Tiedtke et al.(1988)to describe vertical transport of moisture and heat fluxes over stratifically unstable regions.This method describes the vertical transport of moisture and heat fluxes based on turbulent diffusion,which can effectively eliminate convective available potential energy(CAPE)over the areas where shallow convection occurs.While it can describe the physical processes involved in shallow convection,this method also overestimates shallow convection and large biases can be found at the top of the shallow convection.In this study,the GRAPES global model is initialized by using the Final Operational Global Analysis(FNL)analysis from NCEP and continuously runs for eight days. The results of the first two days are discarded due to the influence of initial fields.Simulated heating rate for the 3rd day of the simulation is analyzed here.The results are shown in Fig.1.Zonally average convective heating rate decreases abruptly over the tropics at around 800 hPa(Fig.1a).Comparing Fig.1a with Fig.1b,it is found that the decrease in convective heating is attributed to significant cooling at the top of the shallow convection.Over the areas where shallow convection is active,vertical profiles of convective heating rate at two single points are shown in Figs. 1c and 1d,respectively.One point is located to the east of the inter-tropical convergence zone(ITCZ)at 2°S,100°W,and the other is located inside the South Pacific convergence zone(SPCZ)at 5°S,120°E.Cold bias is found at the top of the simulated shallow convection at both points.The unrealistic simulation of shallow convection is mainly attributed to the turbulent mixing and diffusion method used in the shallow convection scheme of GRAPES,which leads to abnormal warming at low levels and cooling at upper levels.

    Mass flux method has been widely employed in many numerical models to describe cumulus convection(Arakawa and Schubert,1974;Tiedtke,1989;Kain and Fritsh,1990;Zhang and McFarlane,1995;Liu and Ding,2002).The physical processes involved in deep and shallow convection are essentially the same,although the convective cloud depth is different between the two types of convection.Therefore,it is reasonable to apply the mass flux method to shallow convection simulation.

    In order to overcome the problem in the simulation of convective heating rate and optimize the cumulus parameterization scheme,the turbulent mixing and diffusion method currently used in the GRAPES for shallow convection is replaced by the mass flux method.The mass flux method has been applied in the shallow convection scheme of the NCEP GFS operational model(Han and Pan,2011),which provides a good reference for the present study.In addition,some physical processes that are more realistically described in shallow convection scheme,for example,the overshooting and the convection triggering mechanismthat highly depends on the vertical velocity of the environment,are added to the deep convection scheme to keep consistency between deep and shallow convection schemes.These changes are expected to provide more reasonable simulations.

    Fig.1.Zonal mean convective heating rate(K day-1)simulated by the global GRAPES using original convective scheme.(a)Total convection;(b)vertical profile of heating rate at 2°S,100°W;(c)shallow convection only;and(d)vertical profile of heating rate at 5°S,120°E.

    In addition to the modification of shallow convection scheme,in this study we also modify the PBL scheme,which is closely related to shallow convection. At present,most studies about PBL scheme focus on how to better describe the nonlocal effects of large eddies(Troen and Mahrt,1986;Hong and Pan,1996)and how to better represent the strong entrainment at the top of the PBL(Siebesma et al.,2007).Lock et al.(2000)raised the question of coupling between the PBL and stratocumulus.The radiative cooling of the stratocumulus clouds can result in entrainment of dry and cold air mass at the top of the stratocumulus clouds,leading to the so-called cloud-top entrainment instability.Such kinds of instability can be used to describe the coupling-decoupling process between the PBL and stratocumulus.

    The PBL scheme is also modified in this study to match the changes in cumulus convection schemes. The global GRAPES PBL scheme is the same as that used in the medium range forecast model(MRF),which describes large eddies in the PBL by considering the counter-gradient diffusion(Hong and Pan,1996).However,this scheme does not consider the radiative cooling of stratocumulus clouds and strong entrainment at the top of the PBL.These processesmake the PBL more instable by intensifing turbulent mixing and vertical diffusion.It is necessary to consider these processes in the PBL scheme in order to better simulate the coupling between the stratocumulus clouds and PBL processes.In this study,the MRF PBL scheme is modified to include these processes.

    The simulation of stratocumulus cloud is an important index to evaluate whether the cumulus convection,shallow convection,and PBL processes are realistically simulated.Figure 2a shows the monthly average total cloud cover fraction simulated by the GRAPES.The model is initialized on a daily basis using the FNL and continuously runs for eight days in October 2012.October is chosen for this experiment mainly because boreal autumn is the season of maximum stratocumulus.The 5th day forecast of cloud fraction of each run is extracted to produce monthly average cloud cover fraction.Model results are compared with the Moderate-resolution Imaging Spectroradiometer(MODIS)satellite observations(Acker and Leptoukh,2007).It is found that the total stratocumulus clouds over the eastern Pacific are underestimated and the region of clouds shifts westward.

    Further analysis reveals that the underestimation of stratocumulus clouds is attributed to inappropriate simulation of shallow convection.Meanwhile,interaction between PBL processes and stratocumulus clouds is not considered.In the present study,the shallow convection scheme in the GRAPES is modified based on the mass flux method;the deep convection scheme is also modified to be consistent with the shallow convection scheme;physical processes in the PBL are coupled with stratocumulus clouds to make the PBL simulation more reasonable.

    This paper is organized as follows. Section 2 briefly introduces the cumulus scheme and PBL optimization.Simulations of cumulus convection and PBL processes using the modified cumulus and PBL schemes are discussed in Section 3.Summary and discussion are presented in Section 4.

    2.Modification of the cumulus and PBL schemes in the GRAPES global model

    2.1 Convection

    The above analysis clearly indicates that one major problem in the cumulus scheme of GRAPES is that turbulent mixing is applied in the shallow convection scheme,which results in unrealistic cooling at the top of the shallow convection.In this study,the mass fluxmethod proposed by Han and Pan(2011)is applied in the shallow convection scheme of the GRAPES global model.Compared with the turbulent mixing method used previously,the mass flux method can better describe the shallow convection process.

    Fig.2.(a)Monthly mean total cloud fraction for October 2012 simulated by the global GRAPES using the old convective scheme and(b)the corresponding observation from MODIS.

    2.1.1 Shallow convection scheme

    The mass-flux based SAS is employed in the GRAPES model to describe deep convective process. In order to keep consistency between the deep and shallow convection schemes in the GRAPES,the same bulk mass flux method is utilized in the shallow convection scheme and the same entrainment ratio is applied in both the deep and shallow convection schemes. However,compared with deep convection,shallow convection is more closely associated with the PBL processes.Therefore,constraints related to the surface condition must be considered in the closure assumption for the shallow convection.A brief introduction is given below.

    In a bulk mass flux scheme,a single updraftdowndraft branch of the air flow is used to represent an ensemble of cumulus elements within a grid box. The cloud model used in the new shallow convection scheme in this study is expressed as:

    where M represents mass flux;ε is the entrainment rate;δ is the detrainment rate;s is the dry static energy;q and qlare mixing ratio of the environment and updraft,respectively;qsis the saturated humidity;and overbar denotes the horizontal mean on resolvable scale;c is the condensation;r is the precipitation;and L is the latent heat of vaporization.Equations(2)and(3)can be merged to form Eq.(5),since moist static energy can be written as hd=s+Lq.

    Based on their large eddy studies(LES),Siebesma and Cuijpers(1995)suggested that for shallow convection,the turbulent entrainment rate at cloud base is around 2×10-3m-1.In the present study,based on previous results from LES experiments(Grant and Brown,1999;Siebesma et al.,2003)and observation(Raga et al.,1990),turbulent entrainment rate satisfies the following empirical function,

    where ceis empirical constant.Since the entrainment rate generally is larger in shallow convection than in deep convection,ceis set to 0.3 based on the results from several of our sensitivity experiments.

    Shallow convection is closely associated with surface processes and the closure assumption for shallow convection is related to surface properties.Following Grant(2001),the cloud base mass flux of shallow convection is written as:

    where mbis the cloud base mass flux,α is coefficient and set to 0.03,ρ is the air density,and w?is convective-scale vertical velocity,which is defined as

    where g is the gravitational acceleration,θ0is potential temperature at the lowest model level,h is the PBL height,andis the surface virtual kinematic heat flux.

    Identification of the shallow convection starting level is similar to that of deep convection,which is defined as the level of maximum moist static energy between surface and the level of p/ps=0.7.Note that shallow convection starting level is confined below the boundary layer height.The level of free convection(LFC)is defined as the level of cloud base,and the level of cloud top is limited below the level of p/ps= 0.7.Since the major difference between deep and shallow convection is the cloud depth,a deep convection is identified if the cloud depth is larger than 150 hPa.Otherwise,it is identified as a shallow convection.In addition,the effect of convective overshooting is considered for both shallow and deep convection.Convective triggering mechanism is also slightly changed. Detailed description is given next.

    2.1.2 Deep convection scheme

    The modification of deep convection scheme is part of the improvement of convection schemes in this study.Consistency between shallow and deep convection schemes is considered to make it more reasonable. Other modifications include improvement in triggering mechanism,convective entrainment rate,limitation of maximum mass flux at cloud base,consideration of the effect of overshooting,etc.

    In the original scheme,deep convection is triggered once the undiluted parcel is lifted upward from the convective starting level to the LFC level and the depth between the two levels must be less than 150 hPa.The triggering function is a rough estimate of convective inhibition(CIN),which can effectively determine the occurring of convection under the current resolution of global model.However,the empirical constant of 150 hPa is arbitrary.In this study,the empirical constant in the trigger function of the new convection scheme is determined based on grid mean vertical velocity,and the value ranges between 120 and 180 hPa.With this new trigger function,impacts of vertical velocity on the development of cumulus convection are considered,which is more physically reasonable.In addition,the old scheme only considers convective triggering under the condition of no entrainment.In the revised scheme,convective triggering is re-determined after consideration of entrainment.The height of cloud base between the diluted plumes with entrainment and detrainment and the undiluted ones should be less than 25 hPa,an addition for the convective triggering.This method is thus more realistic and avoids too much dry air mass being entrained into the updraft area.

    A realistic entrainment ratio has great impacts on the determination of cloud tops,and efficiently improves simulation of the Madden-Julian oscillation(MJO).Based on Bechtold et al.(2008),the entrainment rate in the present study is expressed as:

    where ε0is the turbulent entrainment rate,which is equivalent to ε in Eq.(6);d1is truncation coefficient,with its magnitude around 10-4;qsand qsbare the saturated mixing ratio at each cloud level and at the cloud base,respectively;u is relatively humidity;and F0and F1are dimensionless variables.

    To avoid unrealistically stronger convection,an arbitrary limitation of maximum mass flux(Mbmax= 0.1 kg m-2s-1)is imposed at the cloud base of the original scheme.Note that this limitation is not appropriate for global models under current resolution(50 km).For this reason,it is adjusted to depend on the vertical resolution and time step,and the Courant Friedrichs Lewy criterion(CFL)is used to determine it as:

    where mbmaxis the maximum mass flux at the cloud base,Δp is the pressure depth of model level of cloud base,and Δt is time step.With this definition,the arbitrary value in the old scheme is removed and strong convection is allowed. Furthermore,it is also correlated with vertical resolution and time step of the model,making it a more reasonable option.

    The level of cloud top in the original scheme is defined as the level of neutral buoyancy.However,neutral buoyancy condition only indicates that acceleration is zero at this level whereas convective overshooting still exists.In this study,the height of overshooting is identified above the neutral buoyancy level where the upward parcel lifted from the neutral buoyancy level costs 10%of CAPE.

    2.2 Boundary layer diffusion

    The boundary layer scheme is also revised in this study since it is highly correlated with shallow conve-ction scheme.One major purpose of this study is to improve the stratocumulus cloud simulation in the eastern Pacific.Stratocumulus clouds develop under combined effects of shallow convection and boundary layer diffusion.In order to better simulate stratocumulus clouds,it is necessary to improve not only shallow convection scheme but also vertical PBL diffusion scheme.

    The current boundary layer scheme used in the global GRAPES is MRF boundary layer diffusion scheme(Troen and Mahrt,1986;Hong and Pan,1996),which is a K-closure nonlocal PBL scheme.In the present study,the cloud-top-driven vertical diffusion is added to the new boundary layer diffusion scheme based on the scheme proposed by Lock et al.(2000)for the coupling between the PBL and the stratocumulus clouds.The turbulent fluxes of boundary layer satisfy the following equation:

    where κ is Von Kalman constant(κ=0.4);hbis the cloud top of stratocumulus;zbis the stratocumulus cloud base,which is the first level where the cloudtop-driven eddy diffusion should be calculated;and VScis the vertical velocity caused by entrainment at the cloud top,which is defined by

    where Vradand Vbrare vertical velocities due to radiative cooling and buoyancy,respectively,and Vbris negligible in the study.Variable Vradis defined by

    where ΔR is the radiative net flux from cloud base to cloud top;cpis specific heat at constant pressure level. The cloud-top entrainment flux is given by

    where c is a constant set to 0.2 based on Moeng et al.(1999).In addition,the cloud-top entrainment instability(CTEI)must satisfy the following condition(Deardorff,1980;Randall,1980)

    where Δθeand Δqtare changes in the equivalent potential temperature and total water content from the cloud starting level across to the cloud top,and c1is constant,which equals 0.7(MacVean and Mason,1990).

    where Km/h(Ri)is momentum/heat turbulent diffusivity;Riis Richardson number;l is mixing length; s the absolute value of the vertical wind shear; and fm/h(Ri)is the stability function of momentum/heat,which is highly dependent on Richardson number.Variable fh(Ri)is modified under the stable condition and expressed as

    while for unstable condition,they are defined as

    Background turbulent diffusivity is also modified,but their impact is negligible,so we will not introduce details here.

    3.Evaluation of the new convection and PBL schemes

    The governing equations of the GRAPES global forecast system are the fully compressible,nonhydrostatic equations describing the complete dynamics of the atmosphere.A hydrostatic reference state of atmosphere is applied.It adopts a semi-implicit and semi-Lagrangian time integration scheme.The spatial discretization uses horizontally staggered Arakawa-C grid and the vertical coordinate is terrain-following with Charney-Philips vertical variation configuration(Chen et al.,2008;Xue and Chen,2008).The horizontal resolution of GRAPES is 50 km and there are 60 levels in the vertical,with the model top at around 36 km.Important physical schemes used in the GRAPES include the RRTMG longwave and shortwave radiation scheme(Mlawer et al.,1997),SAS cumulus parameterization,MRF PBL scheme,WSM6 microphysics scheme(Hong and Lim,2006),gravity wave drag scheme(Kim and Arakawa,1995;Lott and Miller 1997),and CoLM land surface scheme(Dai et al.,2003).

    3.1 Evaluation of the new convection scheme

    Considering the close relation between cumulus convection and PBL physics,both the convection scheme and the boundary layer diffusion scheme have been modified.The impact of the modifications in each scheme is considered separately in order to identify their unique impact.

    Simulation of the total convective heating rate becomes more realistic with the modified cumulus scheme except for an abrupt low bias at around 800 hPa(Fig.3).A weak cooling tendency is found nearthe surface,which may be attributed to rainfall evaporation in shallow convection.Note that there are no constraints for precipitation imposed in the modified shallow convection scheme.This is based on the fact that shallow convection does induce drizzle sometimes in the real world.However,shallow convectioninduced precipitation is very small in the GRAPES(less than 2 mm;figure omitted).The total convective heating rate is larger in the modified scheme than in the old scheme,and the location of maximum heating rate shifts from 450-700 to 350-650 hPa.The level of maximum apparent heat source(Q1)from the modified scheme is more consistent with that proposed by Yanai et al.(1973).

    Fig.3.Zonal mean of the(a,b)total and(c,d)shallow convective heating rate(K day-1)simulated by the GRAPES global model using(a,c)the revised schemes and(b,d)the original schemes.

    The simulated vertical heating profile at some single points shows similar improvement.In order to ensure that the simulations with the modified and the original cumulus schemes are carried out under similar environmental conditions,the vertical profiles of convective heating rate at all model points are extracted from the 24-h forecast outputs.The results at the point(2°S,100°W),which is to the east of the ITCZ,indicate that with the original scheme,shallow convective heating largely occurs at lower levels(950 hPa)while a cooling close to 3 K day-1occurs at 900 hPa. However,with the modified scheme,shallow convective heating becomes dominant and the largest heating occurs at around 925 hPa,a little higher than that for the original scheme.At the SPCZ point(2°S,120°E),shallow convective heating mainly occurs at around 1000 hPa with the original scheme,and a cooling center of about 0.5 K day-1occurs at 700 hPa.This is apparently not realistic.With the modified scheme,shallow convective warming occurs at almost all the levels except those close to the ocean surface,and the convective warming center lifts from the previous 1000 hPa to around 925 hPa(Fig.4b).These results indicate that the revised cumulus scheme can improve the model simulation and produce more reasonable results.

    3.2 Necessity to improve the boundary layer diffusion scheme

    One focus of the present study is to improve the shallow convection scheme.Does this imply any necessity for the improvement of the PBL scheme?Figure 5c is the same as Fig.2a,except that it shows the monthly mean total cloud fraction simulated by GRAPES with the revised cumulus scheme.Comparing Fig.5c with Fig.2a,F(xiàn)ig.5c shows clearly that total cumulus clouds simulated over the eastern Pacific have increased and stratocumulus clouds now cover a larger area.Meanwhile,total clouds over the ITCZ,which are often associated with deep convection,also increase in Fig.5c,probably due to intensified detrainment of high clouds generated by strengthened deep convection.The results in Fig.5c are more consistent with MODIS observations,compared with Fig. 2a,suggesting that the revised cumulus scheme has improved the simulation of both shallow and deep convection.However,a maximum cloud cover area parallel to the ITCZ appears over the tropical southern Pacific in Fig.5c,indicating errors in simulation of a“double ITCZ”.When both the revised cumulus and PBL schemes are applied,the“double ITCZ”error is effectively removed(Fig.5a).Large amounts of clouds with cloud cover fraction higher than 90%appear over the coastal region of Peru in the eastern Pacific and South America,which agrees better to MODIS observations.

    Further analysis indicates that without the revised PBL scheme,turbulent mixing and diffusion at the top of the stratocumulus clouds is not considered,which directly leads to overestimation of stratocumulus clouds.The“double ITCZ”error in the simulation can further affect the simulation of the entire tropical circulation system.Hence,it is not sufficient to solely revise the cumulus scheme in the GRAPES.With modification of both the cumulus and PBL schemes,the“double ITCZ”error is overcome while the high cloud simulation over the ITCZ improves and is more consistent with MODIS observation.Therefore,it is necessary to consider the coupling between the PBL processes and stratocumulus clouds.

    To diagnose the vertical properties of stratocumulus clouds,a vertical cross-section of cloud-water mixing ratio along 5°S is shown in Fig.6.During the boreal autumn,stratocumulus clouds appear relativelylow due to the low inversion layer over the eastern Pacific(Yang et al.,2010).In the present study,we identify the height of cumulus clouds according to cloud top pressure.MODIS observations indicate that the cloud tops are below 800 hPa around 5°S to the east of 90°W.Between 140°and 120°W,the top of stratocumulus clouds is even below 900 hPa.With the original cumulus and PBL schemes,the simulated top of stratocumulus clouds is almost always above 800 hPa.In the regions to the west of 120°W,cloud water mixing ratio is underestimated whereas the simulated cloud top is around 700-800 hPa.The above errors in the simulation are attributed to unrealistically strong vertical transport of cloud water,which is caused by overestimation of shallow convection in the model simulation.Such an overestimation in shallow convection and vertical transport of cloud water completely violates the coupling relations between the PBL processes and cumulus convective activities.With the revised PBL and cumulus schemes,the simulated height of stratocumulus clouds is basically below 800 hPa,and the monthly average height of stratocumulus cloudsover 120°-140°W is around 850 hPa.The base of the simulated stratocumulus clouds is largely below 925 hPa with the old schemes,whereas it is located around 950 hPa with the revised schemes.Apparently,the revised schemes produce a more realistic simulation compared with the original schemes. In addition,in the simulation with the revised schemes,the maximum cloud water mixing ratio occurs over 80°-85°W,corresponding to the cold offshore current area over the ocean.Figure 5d shows clearly a maximum center over this region,while in the simulation with the old schemes,this centercompletely disappears.In general,the vertical properties of stratocumulus clouds are more reasonably simulated with the revised schemes and more consistent with observations.

    Fig.4.The vertical convective heating rate profile(K day-1)at(a)2°S,100°W and(b)5°S,120°E simulated with the revised cumulus scheme(green line)and original cumulus scheme(black line),respectively.

    Fig.5.Monthly mean total cloud cover fraction from simulations and MODIS observations.(a)Simulation with the modified PBL and cumulus schemes,(b)simulation with the original schemes,(c)simulation with the modified convection scheme only,and(d)MODIS observations.

    Fig.6.(a)Monthly mean cloud top pressure(hPa)for October 2012 from MODIS data,(b)vertical cross-section of simulated cloud water mixing ratio with the revised cumulus and PBL schemes along 5°S,and(c)the same as(b)but with the original schemes.

    3.3 Evaluation of the GRAPES global model performance with the modifiedPBLand cumulus schemes

    In Section 3.2,we have evaluated simulations of convective heating rate and stratocumulus cloud properties over the eastern Pacific,which is a region of large amounts of stratocumulus clouds.Impacts of the revised schemes have been illustrated.In this section,we will further investigate how the revised cumulus and PBL schemes affect simulations of global precipitation and large-scale circulation.The capability of the global GRAPES will also be assessed.

    Modification of cumulus and PBL schemes is a focus of the present study,and their impacts on the GRAPES model performance will be evaluated based on the GRAPES simulations over the tropical region of the Northern Hemisphere during summertime.All experiments are conducted for July 2013.The GRAPES is initialized on a daily basis by using the FNL analysis and runs for 8 days.The monthly average results are used for further analysis.

    We are concerned about the precipitation simulation over the tropical region since it directly reflects the impact of the revised schemes.Figure 7 shows the precipitation simulation with the original and the revised scheme as well as the observation from the Tropical Rainfall Measuring Mission(TRMM).With the revised schemes,the spurious precipitation over the Indian Peninsula and the Bay of Bengal is effectively suppressed.Precipitation simulation is also improved over the warm pool region,where the precipitation underestimation is severe with the original schemes.Compared with the original schemes,the revised schemes can produce precipitation simulation that is more consistent with the TRMM observation(Kummerow et al.,1998).

    Figure 8 presents time evolution of the summertime western Pacific subtropical high,including its westward extension and northward movement.The results are compared with the FNL analysis.In general,the simulated subtropical high becomes stronger with the revised schemes than with the original schemes,particularly in its westward extension.During 17-25 July 2013,the simulated 5880-gpm contour extends to nearby 120°E with the revised schemes,whereas it is not clear at all with the original schemes.However,the subtropical high intensity is still underestimated in the GRAPES simulation even with the revised schemes,indicating that further analysis should be done to explore the reasons behind such an underestimation.It is probably due to inconsistencies between physical schemes and/or problems in the coupling between the model dynamics and physics.

    The GRAPES evaluation tools(GET),which is a system used to verify the operational GRAPES,is employed to evaluate results of batch tests of GRAPES in the present study and to obtain objective assessment of the revised schemes.Figure 9a shows the correlation coefficient of 500-hPa geopotential height anomalies between GRAPES simulation and FNL analysis for July 2013.This month is taken as a typical summer month in the Northern Hemisphere.The red boxes in Figs.9a and 9b indicate the range of 95%confidence level.The red lines represent results with the revised cumulus and PBL schemes,and the black lines for control run with the original schemes.Figure 9a shows that the correlation coefficient of 500-hPa geopotential height anomalies becomes significantly large after the 4th day of forecast and remains large until the 8th day of forecast.This indicates that with the revised cumulus and PBL schemes,the GRAPES forecast is valid for 7 days in advance.In addition,the root mean square errors(RMSEs)of wind fields in the tropics with the revised schemes and the original schemes also indicate that the revised schemes can improve simulation of winds,and the improvements are statistically significant.The RMSEs of winds can be reduced by approximately 1 m s-1.

    Statistical analysis(figure omitted)indicates that with the revised schemes,the simulation of geopotential height and wind fields are improved over the Northern Hemisphere and the tropical region,and the improvements are statistically significant.

    4.Conclusions and discussion

    The cumulus and PBL schemes are two important physical schemes in numerical modeling study.They play a key role in maintaining the heat,moisture,and momentum fluxes in the earth system.The heating of the atmosphere caused by latent heat release during the cumulus convection and the radiative cooling of stratocumulus clouds over the eastern Pacific are two important processes in numerical weather forecast.In the original GRAPES,a turbulent mixing and diffusion method is used to describe shallow convection. This method overestimates shallow convection in the PBL and produces unreasonable convective heating rate in the vertical.In this study,both the shallow convection and PBL schemes are modified to realistically simulate related physical processes in the atmosphere.Batch experiments have been conducted by using both the revised and the original schemes. Diagnostic analyses of the results from these experiments indicate that the modified cumulus and PBL schemes have significantly improved the performance of the global GRAPES.Major results are as follows.

    (1)The original shallow convection scheme in theglobal GRAPES generates unrealistic cooling at the top of the shallow convection.As a result,the total convective heating rate becomes unrealistically small at the top of shallow convection.Meanwhile,stratocumulus clouds in the coastal region of Peru in eastern Pacific are underestimated while the cloud top is abnormally higher than observation.

    Fig.8.Evolution of the western Pacific subtropical high in July 2013.(a,c,e)Time-longitude cross-sections of the WPSH simulated with the revised and the original schemes respectively and from FNL analysis,and(b,d,f)time-latitude cross-sections of the WPSH simulated with the revised and the original schemes and from FNL.

    (2)In the present study,a mass flux method is applied in the shallow convection scheme to replace the original method that is based on turbulent diffusion.The closure assumption for shallow convection considers the surface properties,and deep convection is also adjusted to become more reasonable.This change ensures that deep and shallow convection schemes are consistent with each other. The coupling between the PBL processes and stratocumulus clouds,i.e.,the CTEI effect,is also well considered in the present study.Under the condition of stable stratification that is favorable for the formation of stratocumulus clouds,the stability function is adjusted and becomes more realistic.

    Fig.9.(a)500-hPa geopotential height anomaly correlation coefficient in the Northern Hemisphere and(b)root mean square error(RMSE)of winds at 850 hPa over the tropics.Red lines represent results from the revised schemes and black lines from the original schemes.

    (3)With only modification of convection scheme,the simulation of convective heating rate becomes more reasonable but the stratocumulus clouds over the eastern Pacific are overestimated and cover an unrealistically large area.Hence,it is necessary to modify the PBL scheme simultaneously.With modification of both the cumulus and PBL schemes,stratocumulus clouds are well simulated and the vertical structure of the stratocumulus clouds is consistent with observations.

    (4)With the modified cumulus and PBL schemes,spurious precipitation over the tropics is suppressed and the underestimation of precipitation over the warm pool region is improved.These modifications also increase the 500-hPa geopotential height anomaly correlation coefficient and effectively extend the valid forecast period.Circulations at the tropics are also improved with the revised schemes.The statistical scorecard indicates that the performance of the global GRAPES with the revised schemes has been significantly improved. The simulations of geopotential height and wind fields are improved over the Northern Hemisphere and the tropical region.These results demonstrate that the modification of cumulus and PBL schemes in the present study could contribute to the development and improvement of the global GRAPES.

    Acknowledgments.The authors would like to express their gratitude to Drs.Pan Hua-lu and Shen Xueshun for helpful discussion and comments.We are also grateful to Dr.Zhao Bin for supporting us in running the GET evaluation system.

    Acker,J.G.,and G.Leptoukh,2007:Online analysis enhances use of NASA earth science data.Eos.Trans. Amer.Geophys.Union,88,14-17.

    Arakawa,A.,and W.H.Schubert,1974:Interaction of a cumulus cloud ensemble with the large-scale environment,Part I.J.Atmos.Sci.,31,674-701.

    Bechtold,P.,M.K¨ohler,T.Jung,et al.,2008:Advances in simulating atmospheric variability with the ECMWF model:From synoptic to decadal timescales.Quart.J.Roy.Meteor.Soc.,134,1337-1351.

    ChenDehui,XueJishan,YangXuesheng,etal.,2008:New generation of multi-scale NWP system(GRAPES):General scientific design.Chin.Sci. Bull.,53,3433-3445.

    Dai,Y.J.,X.B.Zeng,R.E.Dickinson,et al.,2003:The common land model.Bull.Amer.Meteor.Soc.,84,1013-1023.

    Deardorff,J.W.,1980:Cloud top entrainment instability.J.Atmos.Sci.,37,131-147.

    Deng Hua,Xue Jishan,Xu Haiming,et al.,2008:Study of different cumulus parameterization schemes of GRAPES-Meso model in simulation of convection provocation.J.Trop.Meteor.,24,327-334.(in Chinese)

    Grant,A.L.,and A.R.Brown,1999:A similarity hypothesis for shallow-cumulus transports.Quart.J. Roy.Meteor.Soc.,125,1913-1936.

    Grant,A.L.M.,2001:Cloud-base fluxes in the cumuluscapped boundary layer.Quart.J.Roy.Meteor. Soc.,127,407-421.

    Han,J.,and H.L.Pan,2011:Revision of convection and vertical diffusion schemes in the NCEP global forecast system.Wea.Forecasting,26,520-533.

    Hong,S.Y.,and H.L.Pan,1996:Nonlocal boundary layer vertical diffusion in a medium-range forecast model.Mon.Wea.Rev.,124,2322-2339.

    Hong,S.-Y.,and J.-O.J.Lim,2006:The WRF singlemoment 6-class microphysics scheme(WSM6).J. Korean Meteor.Soc.,42,129-151.

    Kain,J.S.,and J.M.Fritsh,1990:A one-dimensional entraining/detraining plume model and its application in convective parameterization.J.Atmos.Sci.,47,2784-2802.

    Kim,Y.-J.,and A.Arakawa,1995:Improvement of orographic gravity wave parameterization using a mesoscale gravity-wave model.J.Atmos.Sci.,52,1875-1902.

    Kummerow,C.,W.Barnes,T.Kozu,et al.,1998:The tropical rainfall measuring mission(TRMM)sensor package.J.Atmos.Ocean.Technol.,15,809-817.

    Li Lijuan and Wang Bin,2010:Influence of two different convection schemes on the radiative energy budget in GAMIL1.0.Acta Meteor.Sinica,24,318-327.

    Liu Kun,Liu Yimin,and Wu Guoxiong,2010:The impacts of the modified Tiedtke cumulus convective parameterization scheme on the tropical rainfall simulation in SAMIL model. Chinese J.Atmos. Sci.,34,163-174.(in Chinese)

    Liu Yiming and Ding Yihui,2002:Modified mass flux cumulus convective parameterization scheme and its simulation experiment.Part I:Mass flux scheme and its simulation of the 1991 flood event. Acta Meteor.Sinica,16,37-49.

    Lock,A.P.,A.R.Brown,M.R.Bush,et al.,2000:A new boundary layer mixing scheme.Part I:Scheme description and single-column model tests. Mon. Wea.Rev.,128,3187-3199.

    Lott,F(xiàn).,and M.J.Miller,1997:A new subgrid-scale orographic drag parameterization:Its formulation and testing.Quart.J.Roy.Meteor.Soc.,123,101-127. MacVean,M.K.,and P.J.Mason,1990:Cloud-top entrainment instability through small-scale mixing and its parameterization in numerical models.J. Atmos.Sci.,47,1012-1030.

    Mlawer,E.J.,S.J.Taubman,P.D.Brown,et al.,1997:Radiative transfer for inhomogeneous atmospheres:RRTM,a validated correlated-k model for the longwave.J.Geophys.Res.,102,16663-16682.

    Moeng,C.H.,P.P.Sunllivan,and B.Stevens,1999:Including radiative effects in an entrainment rate formula for buoyancy-driven PBLs.J.Atmos.Sci.,56,1031-1049.

    Nie,J.,and Z.M.Kuang,2012:Responses of shallow cumulus convection to large-scale temperature and moisture perturbations:A comparison of large-eddy simulations and a convective parameterization based on stochastically entraining parcels.J.Atmos.Sci.,69,1936-1956.

    Pan,H.-L.,and W.-S.Wu,1995:Implementing a Mass Flux Convective Parameterization Package for the NMC Medium-Range Forecast Model.NMC Office Note 409,Washington,DC,1-40.

    Raga,G.B.,J.B.Jensen,and M.B.Baker,1990:Characteristics of cumulus band clouds off the coast of Hawaii.J.Atmos.Sci.,47,338-356.

    Randall,D.A.,1980:Conditional instability of the first kind upside-down.J.Atmos.Sci.,37,125-130.

    Siebesma,A.P.,and J.W.M.Cuijpers,1995:Evaluation of parametric assumptions for shallow cumulus convection.J.Atmos.Sci.,52,650-666.

    Siebesma,A.P.,C.S.Bretherton,A.Brown,et al.,2003:A large eddy simulation intercomparison study of shallow cumulus convection.J.Atmos.Sci.,60,1201-1219.

    Siebesma,A.P.,P.M.M.Soares,and J.Teixeira,2007:A combined eddy-diffusivity mass flux approach for the convective boundary layer.J.Atmos.Sci.,64,1230-1248.

    Tiedtke,M.,1989:A comprehensive mass flux scheme for cumulus parameterization in large-scale models. Mon.Wea.Rev.,117,1779-1800.

    Tiedtke,M.,W.A.Heckley,and J.Slingo,1988:Tropical forecasting at ECMWF:The influence of physical parametrization on the mean structure of forecasts and analyses.Quart.J.Roy.Meteor.Soc.,114,639-664.

    Troen,I.,and L.Mahrt,1986:A simple model of the atmospheric boundary layer:Sensitivity to surface evaporation.Bound.-Layer Meteor.,37,129-148.

    von Salzen,K.,N.A.McFarlane,and M.Lazare,2005: The role of shallow convection in the water and energy cycles of the atmosphere.Climate Dyn.,25,671-688.

    Xu Daosheng,Chen Zitong,Dai Guangfeng,et al.,2014:The influence of an improved cumulus parameterization scheme on typhoon forecast from GRAPES model.J.Trop.Meteor.,30,210-218.(in Chinese)

    Xue Jishan and Chen Dehui,2008:Scientific Design and Application of Numerical Prediction System GRAPES.Science Press,Beijing,334-335.(in Chinese)

    Yanai,M.,S.Esbensen,and J.-H.Chu,1973:Determination of bulk properties of tropical cloud clusters from large-scale heat and moisture budgets.J.Atmos.Sci.,30,611-627.

    Yang,F(xiàn).L.,H.L.Pan,S.K.Krueger,et al.,2010:Evaluation of the NCEP Global Forecast System at the ARM SGP site.Mon.Wea.Rev.,134,3668-3690.

    Zhang,G.J.,and N.A.McFarlane,1995:Sensitivity of climate simulations to the parameterization of cumulus convection in the Canadian Climate Centre general circulation model.Atmos.-Ocean,33,407-446.

    Liu Kun,Chen Qiying,and Sun Jian,2015:Modification of cumulus convection and planetary boundary layer schemes in the GRAPES global model.J.Meteor.Res.,29(5),806-822,

    10.1007/s13351-015-5043-5.

    (Received May 10,2015;in final form August 16,2015)

    猜你喜歡
    劉琨孫健
    A pressure-calibration method of wavelength modulation spectroscopy in sealed microbial growth environment*
    There Is No Perfect Solution
    西部論叢(2017年7期)2017-11-20 06:59:42
    量臀定碼
    三月三(2017年6期)2017-07-01 08:13:42
    量臀定碼
    三月三(2017年6期)2017-07-01 07:31:42
    枕戈
    也有會(huì)吹的
    三月三(2016年12期)2016-12-27 18:09:10
    Modelling hydrodynamic processes in tidal stream energy extraction*
    漫畫四幅
    安身立命的資本
    安身立命的資本
    婷婷亚洲欧美| 成年版毛片免费区| 日韩人妻高清精品专区| 国产主播在线观看一区二区| 一夜夜www| 午夜激情福利司机影院| 国产aⅴ精品一区二区三区波| 在线十欧美十亚洲十日本专区| 身体一侧抽搐| 黄片小视频在线播放| 国产黄a三级三级三级人| 精品人妻视频免费看| 国产乱人伦免费视频| 国产精品电影一区二区三区| 国产在线男女| 波多野结衣巨乳人妻| 老鸭窝网址在线观看| 国产精品伦人一区二区| 少妇熟女aⅴ在线视频| 日韩精品中文字幕看吧| 欧美一区二区亚洲| 国产白丝娇喘喷水9色精品| 亚洲成人久久性| 免费av毛片视频| 亚洲七黄色美女视频| 成人性生交大片免费视频hd| 亚洲成av人片在线播放无| 一进一出抽搐gif免费好疼| 在线免费观看的www视频| 性插视频无遮挡在线免费观看| 特大巨黑吊av在线直播| 亚洲国产日韩欧美精品在线观看| 亚洲精品一区av在线观看| 欧美精品国产亚洲| 日韩 亚洲 欧美在线| 怎么达到女性高潮| 中文字幕av成人在线电影| 国产精品,欧美在线| 欧美日本亚洲视频在线播放| 一区福利在线观看| 网址你懂的国产日韩在线| 欧美激情国产日韩精品一区| 色5月婷婷丁香| 最近最新中文字幕大全电影3| 99久久精品国产亚洲精品| 校园春色视频在线观看| 午夜福利成人在线免费观看| 美女高潮的动态| 国产精品人妻久久久久久| 永久网站在线| www.999成人在线观看| 亚洲五月婷婷丁香| 看免费av毛片| 欧美又色又爽又黄视频| 天美传媒精品一区二区| 免费黄网站久久成人精品 | 宅男免费午夜| 成人午夜高清在线视频| 亚洲精品456在线播放app | 国产欧美日韩精品亚洲av| 免费看光身美女| 国产亚洲精品av在线| 婷婷亚洲欧美| 两个人的视频大全免费| 亚洲精品成人久久久久久| 少妇的逼好多水| 精品欧美国产一区二区三| 亚洲人成网站高清观看| 日韩人妻高清精品专区| 69人妻影院| 女人被狂操c到高潮| 亚洲无线观看免费| 男人舔女人下体高潮全视频| 国模一区二区三区四区视频| 高潮久久久久久久久久久不卡| 两个人的视频大全免费| 少妇裸体淫交视频免费看高清| 美女高潮的动态| 久久精品国产亚洲av涩爱 | 免费av观看视频| 久久久久精品国产欧美久久久| 国产aⅴ精品一区二区三区波| 99久久精品一区二区三区| 搞女人的毛片| 久久国产乱子伦精品免费另类| 日韩欧美免费精品| 精品午夜福利视频在线观看一区| 国内精品一区二区在线观看| 国产高清激情床上av| 日本 欧美在线| 乱码一卡2卡4卡精品| 亚洲人成网站在线播放欧美日韩| 女生性感内裤真人,穿戴方法视频| 久9热在线精品视频| 亚洲av成人精品一区久久| 国产午夜福利久久久久久| 国产国拍精品亚洲av在线观看| 能在线免费观看的黄片| 欧美绝顶高潮抽搐喷水| 国产精品,欧美在线| 给我免费播放毛片高清在线观看| 桃色一区二区三区在线观看| 欧美xxxx性猛交bbbb| 中国美女看黄片| 99久久无色码亚洲精品果冻| 好男人在线观看高清免费视频| 亚洲欧美精品综合久久99| 亚洲人成网站高清观看| 成人精品一区二区免费| av在线观看视频网站免费| 国产精品女同一区二区软件 | 国产亚洲精品综合一区在线观看| 人人妻人人澡欧美一区二区| 嫩草影院新地址| 成人亚洲精品av一区二区| 婷婷精品国产亚洲av| 我的女老师完整版在线观看| 中文字幕精品亚洲无线码一区| 少妇的逼水好多| 欧美日本亚洲视频在线播放| 亚洲五月天丁香| 国产 一区 欧美 日韩| 三级男女做爰猛烈吃奶摸视频| 村上凉子中文字幕在线| 日本黄色视频三级网站网址| 国产精品一区二区三区四区免费观看 | 精品福利观看| 九九热线精品视视频播放| 亚洲综合色惰| 亚洲内射少妇av| 亚洲欧美日韩卡通动漫| 久9热在线精品视频| 日本三级黄在线观看| 国产伦在线观看视频一区| 亚洲精品日韩av片在线观看| 性插视频无遮挡在线免费观看| 又黄又爽又刺激的免费视频.| 99久久久亚洲精品蜜臀av| 欧美性猛交黑人性爽| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 国内毛片毛片毛片毛片毛片| 日韩欧美国产在线观看| 亚洲一区二区三区色噜噜| 国产一区二区三区在线臀色熟女| 中文在线观看免费www的网站| 国产精品综合久久久久久久免费| 丰满乱子伦码专区| 日韩亚洲欧美综合| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩东京热| 中文字幕av在线有码专区| 婷婷精品国产亚洲av在线| 深夜精品福利| 亚洲七黄色美女视频| 国产精品久久电影中文字幕| 国产伦在线观看视频一区| 我的老师免费观看完整版| 精品乱码久久久久久99久播| a级毛片免费高清观看在线播放| 美女cb高潮喷水在线观看| 天堂网av新在线| 亚洲第一区二区三区不卡| 他把我摸到了高潮在线观看| 日韩欧美精品免费久久 | 国产精品亚洲一级av第二区| 亚洲成人中文字幕在线播放| 极品教师在线免费播放| 性欧美人与动物交配| av黄色大香蕉| 欧美三级亚洲精品| 少妇人妻一区二区三区视频| 一个人免费在线观看电影| 亚洲精品456在线播放app | 中文资源天堂在线| 亚洲国产高清在线一区二区三| 少妇丰满av| 制服丝袜大香蕉在线| 99久久精品国产亚洲精品| 午夜日韩欧美国产| 国产乱人视频| 国产国拍精品亚洲av在线观看| 美女高潮的动态| 精品不卡国产一区二区三区| 国产三级黄色录像| 高清在线国产一区| 欧美又色又爽又黄视频| 欧美成狂野欧美在线观看| netflix在线观看网站| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清专用| 亚洲国产精品合色在线| 国产精品日韩av在线免费观看| 国产v大片淫在线免费观看| 国产午夜福利久久久久久| 超碰av人人做人人爽久久| 国产aⅴ精品一区二区三区波| 午夜福利在线在线| 久久天躁狠狠躁夜夜2o2o| 成人精品一区二区免费| 国产精品久久久久久精品电影| 午夜精品久久久久久毛片777| 少妇熟女aⅴ在线视频| 亚洲熟妇中文字幕五十中出| 丰满人妻熟妇乱又伦精品不卡| 成人国产一区最新在线观看| 欧美黑人巨大hd| 日日夜夜操网爽| 18禁黄网站禁片免费观看直播| 青草久久国产| 少妇高潮的动态图| 亚洲狠狠婷婷综合久久图片| 精品一区二区三区av网在线观看| 最新中文字幕久久久久| 三级毛片av免费| 黄色日韩在线| 免费看美女性在线毛片视频| 看片在线看免费视频| 在线观看一区二区三区| 可以在线观看的亚洲视频| 精品福利观看| 国产精品自产拍在线观看55亚洲| 欧美潮喷喷水| 亚洲av美国av| 成人av在线播放网站| 99久国产av精品| 露出奶头的视频| 久久久色成人| www.熟女人妻精品国产| 国产成人福利小说| 日韩高清综合在线| 日韩欧美国产在线观看| 国产一区二区三区视频了| 十八禁网站免费在线| 又爽又黄无遮挡网站| av在线蜜桃| 91久久精品国产一区二区成人| 99久久精品一区二区三区| 亚洲美女视频黄频| 国产69精品久久久久777片| 999久久久精品免费观看国产| 日本免费a在线| 脱女人内裤的视频| 国产精品一及| 欧美另类亚洲清纯唯美| 美女大奶头视频| 久99久视频精品免费| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| 欧美乱色亚洲激情| 国产精品一区二区三区四区久久| 午夜久久久久精精品| 亚洲国产精品999在线| 国产人妻一区二区三区在| 亚洲av五月六月丁香网| 少妇裸体淫交视频免费看高清| 精品国内亚洲2022精品成人| av欧美777| 国产亚洲欧美98| 999久久久精品免费观看国产| 日韩av在线大香蕉| 久久久久久久久久成人| 九九热线精品视视频播放| 九色国产91popny在线| 久久国产乱子伦精品免费另类| 国产成人影院久久av| 国产精品乱码一区二三区的特点| 国产精品爽爽va在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 国产白丝娇喘喷水9色精品| 91在线观看av| 久久人妻av系列| 看十八女毛片水多多多| 国产蜜桃级精品一区二区三区| 久久性视频一级片| 亚洲在线观看片| 91午夜精品亚洲一区二区三区 | 男人舔女人下体高潮全视频| 欧美日韩综合久久久久久 | 夜夜爽天天搞| 亚洲精品乱码久久久v下载方式| 看十八女毛片水多多多| 美女黄网站色视频| 男人和女人高潮做爰伦理| 久久久久久久久中文| 人人妻人人澡欧美一区二区| 51国产日韩欧美| 变态另类成人亚洲欧美熟女| 老鸭窝网址在线观看| 亚洲,欧美,日韩| 国产亚洲欧美在线一区二区| 亚洲三级黄色毛片| 美女高潮喷水抽搐中文字幕| av黄色大香蕉| 亚洲综合色惰| 亚洲男人的天堂狠狠| av天堂在线播放| 亚洲欧美日韩高清专用| 欧美日韩黄片免| 制服丝袜大香蕉在线| 亚洲天堂国产精品一区在线| 国模一区二区三区四区视频| 国产单亲对白刺激| 国语自产精品视频在线第100页| 欧美午夜高清在线| 久久久色成人| 韩国av一区二区三区四区| 色综合亚洲欧美另类图片| 非洲黑人性xxxx精品又粗又长| 变态另类成人亚洲欧美熟女| or卡值多少钱| 亚洲精品乱码久久久v下载方式| 日本免费一区二区三区高清不卡| 成人特级av手机在线观看| 内地一区二区视频在线| 国产av在哪里看| 久久久久九九精品影院| 中文资源天堂在线| 精品乱码久久久久久99久播| 日本撒尿小便嘘嘘汇集6| 欧美乱色亚洲激情| 亚洲av成人不卡在线观看播放网| 亚洲熟妇中文字幕五十中出| 欧美zozozo另类| 女人被狂操c到高潮| 国产亚洲精品久久久久久毛片| 国产熟女xx| 全区人妻精品视频| 久久99热这里只有精品18| 禁无遮挡网站| 小蜜桃在线观看免费完整版高清| 深夜a级毛片| 国产亚洲av嫩草精品影院| 9191精品国产免费久久| 国产伦精品一区二区三区视频9| 男人的好看免费观看在线视频| 国产色婷婷99| 三级国产精品欧美在线观看| 欧美黑人欧美精品刺激| 熟女人妻精品中文字幕| 免费av观看视频| 天堂√8在线中文| 亚洲在线自拍视频| 深夜精品福利| 99精品在免费线老司机午夜| 免费高清视频大片| 国产aⅴ精品一区二区三区波| 一进一出好大好爽视频| 久久久久精品国产欧美久久久| 国内揄拍国产精品人妻在线| 亚洲人成网站在线播| 少妇丰满av| 午夜福利成人在线免费观看| 亚洲美女视频黄频| 中文字幕人妻熟人妻熟丝袜美| av欧美777| 精品久久久久久成人av| 精华霜和精华液先用哪个| 青草久久国产| 在线观看舔阴道视频| 国产精品一区二区性色av| 嫩草影院精品99| 美女高潮的动态| 亚洲 欧美 日韩 在线 免费| 偷拍熟女少妇极品色| 男女那种视频在线观看| 色在线成人网| 男插女下体视频免费在线播放| 小说图片视频综合网站| 久久性视频一级片| 成人av一区二区三区在线看| 成人性生交大片免费视频hd| 日韩欧美国产一区二区入口| 在线看三级毛片| 久久精品人妻少妇| 中国美女看黄片| 精品久久久久久久人妻蜜臀av| 亚洲真实伦在线观看| 性色av乱码一区二区三区2| 国产 一区 欧美 日韩| 欧美区成人在线视频| 久久精品综合一区二区三区| 校园春色视频在线观看| 国内少妇人妻偷人精品xxx网站| 啪啪无遮挡十八禁网站| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女| 欧美绝顶高潮抽搐喷水| 欧美激情在线99| 人妻制服诱惑在线中文字幕| 啪啪无遮挡十八禁网站| 在线天堂最新版资源| 国产亚洲av嫩草精品影院| 亚洲中文字幕一区二区三区有码在线看| 久久午夜亚洲精品久久| 又粗又爽又猛毛片免费看| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 免费电影在线观看免费观看| 欧美最黄视频在线播放免费| 嫩草影院精品99| 国产精品99久久久久久久久| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 青草久久国产| 91午夜精品亚洲一区二区三区 | 亚洲不卡免费看| 久久精品国产亚洲av香蕉五月| 尤物成人国产欧美一区二区三区| 波多野结衣高清作品| 中文字幕久久专区| 亚洲狠狠婷婷综合久久图片| 757午夜福利合集在线观看| 欧美成人a在线观看| 女同久久另类99精品国产91| 别揉我奶头~嗯~啊~动态视频| 国产高清视频在线播放一区| 精品无人区乱码1区二区| 亚洲自拍偷在线| 一卡2卡三卡四卡精品乱码亚洲| 久久婷婷人人爽人人干人人爱| www.熟女人妻精品国产| 久久6这里有精品| 91狼人影院| 久久6这里有精品| 91狼人影院| 欧美日韩中文字幕国产精品一区二区三区| 我要看日韩黄色一级片| 美女cb高潮喷水在线观看| 亚洲,欧美,日韩| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av天美| 无遮挡黄片免费观看| 日韩 亚洲 欧美在线| 身体一侧抽搐| 久久天躁狠狠躁夜夜2o2o| 天天一区二区日本电影三级| 免费高清视频大片| 成人欧美大片| 99国产综合亚洲精品| 日韩精品中文字幕看吧| 精品人妻熟女av久视频| 精品一区二区三区人妻视频| 无人区码免费观看不卡| 精品久久久久久久久亚洲 | 免费av观看视频| 亚洲七黄色美女视频| 变态另类成人亚洲欧美熟女| av女优亚洲男人天堂| 99国产极品粉嫩在线观看| 精品午夜福利视频在线观看一区| 精品乱码久久久久久99久播| 国产激情偷乱视频一区二区| 国产高清有码在线观看视频| 欧美日韩综合久久久久久 | 亚洲欧美精品综合久久99| 国产精品久久久久久精品电影| 国产色婷婷99| 免费观看人在逋| av天堂中文字幕网| 熟女人妻精品中文字幕| 美女免费视频网站| 又黄又爽又刺激的免费视频.| 久久久久久久午夜电影| 国产探花极品一区二区| 97热精品久久久久久| 国产高潮美女av| 欧美丝袜亚洲另类 | 色播亚洲综合网| 亚洲第一区二区三区不卡| 麻豆成人午夜福利视频| 欧美一区二区精品小视频在线| 欧美丝袜亚洲另类 | a在线观看视频网站| ponron亚洲| 亚洲精品456在线播放app | 亚洲最大成人手机在线| 国产爱豆传媒在线观看| 亚洲一区高清亚洲精品| 国产精品一及| 99久久99久久久精品蜜桃| 91在线精品国自产拍蜜月| 在线天堂最新版资源| 国产成人av教育| 男人舔女人下体高潮全视频| 中文字幕久久专区| 嫩草影院新地址| 在线天堂最新版资源| 精品一区二区三区视频在线| 亚洲精品日韩av片在线观看| 亚洲人成网站高清观看| 中文字幕熟女人妻在线| 人妻制服诱惑在线中文字幕| 欧美成狂野欧美在线观看| 最好的美女福利视频网| 日韩欧美国产在线观看| 亚洲狠狠婷婷综合久久图片| av福利片在线观看| 午夜激情欧美在线| 久久精品国产亚洲av香蕉五月| 成人欧美大片| 91久久精品国产一区二区成人| 精品一区二区三区人妻视频| 精品无人区乱码1区二区| 国产精品av视频在线免费观看| 成人欧美大片| 亚洲黑人精品在线| 精品人妻偷拍中文字幕| 亚洲人成网站在线播| 天堂√8在线中文| 99国产精品一区二区三区| 久久性视频一级片| 久久国产乱子伦精品免费另类| 乱码一卡2卡4卡精品| 久久久精品欧美日韩精品| 黄色一级大片看看| 欧美色视频一区免费| 日本与韩国留学比较| 日韩欧美国产在线观看| 哪里可以看免费的av片| 成人国产综合亚洲| 精品人妻熟女av久视频| 亚洲美女视频黄频| 少妇的逼水好多| 日本 欧美在线| 嫩草影视91久久| 一级av片app| 美女黄网站色视频| 欧美+日韩+精品| 青草久久国产| 美女免费视频网站| 成人一区二区视频在线观看| 免费人成视频x8x8入口观看| 国产乱人视频| 国产av一区在线观看免费| or卡值多少钱| 亚洲成人免费电影在线观看| 首页视频小说图片口味搜索| 天堂动漫精品| 99国产精品一区二区三区| 真人一进一出gif抽搐免费| 丁香欧美五月| 夜夜躁狠狠躁天天躁| 日本a在线网址| 国内精品美女久久久久久| 女人被狂操c到高潮| 日本精品一区二区三区蜜桃| 久久久久久久久久成人| 最近在线观看免费完整版| 一级a爱片免费观看的视频| 免费观看精品视频网站| 啪啪无遮挡十八禁网站| 国产极品精品免费视频能看的| 国产人妻一区二区三区在| av天堂在线播放| 丝袜美腿在线中文| 成年女人毛片免费观看观看9| 亚洲国产精品sss在线观看| 国产v大片淫在线免费观看| 国产蜜桃级精品一区二区三区| 国产成人a区在线观看| a级毛片a级免费在线| 久久久久国产精品人妻aⅴ院| 男插女下体视频免费在线播放| 嫩草影视91久久| 99热这里只有是精品在线观看 | 国模一区二区三区四区视频| 九九在线视频观看精品| 麻豆成人午夜福利视频| 亚洲精品一区av在线观看| 精品日产1卡2卡| 草草在线视频免费看| 亚洲人成网站高清观看| 色综合站精品国产| 亚洲无线观看免费| 日本 欧美在线| 91久久精品电影网| 一级av片app| 宅男免费午夜| 色综合婷婷激情| 99国产综合亚洲精品| 欧美激情久久久久久爽电影| 最近最新免费中文字幕在线| 国产老妇女一区| 久久久国产成人免费| 波多野结衣高清无吗| 国产成人欧美在线观看| 久久久久久久久大av| a级一级毛片免费在线观看| 国产高清有码在线观看视频| 999久久久精品免费观看国产| 在线免费观看不下载黄p国产 | 免费电影在线观看免费观看| 亚洲中文日韩欧美视频| 久久久久久久精品吃奶| 亚洲一区二区三区色噜噜| av女优亚洲男人天堂| av天堂在线播放| 中出人妻视频一区二区| 国产不卡一卡二| 午夜久久久久精精品| 精品久久国产蜜桃| 国产精品亚洲美女久久久| 免费人成视频x8x8入口观看| 国产精品1区2区在线观看.| 91午夜精品亚洲一区二区三区 | 久久午夜亚洲精品久久| 三级国产精品欧美在线观看| 很黄的视频免费| 婷婷色综合大香蕉| 国产色爽女视频免费观看| 十八禁网站免费在线| 欧美高清成人免费视频www| 亚洲av一区综合|