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

    Diagnostic Study of Global Energy Cycle of the GRAPES Global Model in the Mixed Space-Time Domain

    2014-12-14 06:58:00ZHAOBin趙濱andZHANGBo張博
    Journal of Meteorological Research 2014年4期
    關(guān)鍵詞:張博

    ZHAO Bin(趙濱)and ZHANG Bo(張博)

    National Meteorological Center,Beijing 100081

    Diagnostic Study of Global Energy Cycle of the GRAPES Global Model in the Mixed Space-Time Domain

    ZHAO Bin(趙濱)and ZHANG Bo?(張博)

    National Meteorological Center,Beijing 100081

    Some important diagnostic characteristics for a model’s physical background are reflected in the model’s energy transport,conversion,and cycle.Diagnosing the atmospheric energy cycle is a suitable way towards understanding and improving numerical models.In this study,formulations of the“Mixed Space-Time Domain”energy cycle are calculated and the roles of stationary and transient waves within the atmospheric energy cycle of the Global-Regional Assimilation and Prediction System(GRAPES)model are diagnosed and compared with the NCEP analysis data for July 2011.Contributions of the zonal-mean components of the energy cycle are investigated to explain the performance of numerical models.

    Mixed Space-Time Domain energy cycle,energy reservoir,energy conversion,stationary wave, transient wave,GRAPES model

    1.Introduction

    Atmospheric systems such as cyclones and anticyclones are measured in terms of their kinetic energy (KE)(Storch et al.,2012),and those systems that are intensifying or weakening are often defined as gaining or losing KE(Luo,1994).Therefore,knowledge regarding the sources and sinks of KE in this context is very important(Li and Zhu,1995;Gao et al.,2006).

    The total energy of the whole atmosphere will remain constant under adiabatic motion and the only source and sink of KE should be the available potential energy(APE).The APE can be considered as the difference between the total potential energy and minimum potential energy,and we can state that APE should vanish when the system’s distribution becomes universally horizontal and stable.However,the general motion of the atmosphere is not adiabatic;the presence of friction should alter KE directly,and conversion between APE and KE should be considered as the generation and dissipation energy within the whole life cycle of the atmospheric motions.

    Supported by the National Nature Science Foundation of China(41305091)and China Meteorological Administration Special Fund for Numerical Prediction(GRAPES).

    ?Corresponding author:zhaob@cma.gov.cn.

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

    Lorenz(1955)developed a life cycle for energy conversions and defined a series of formulations to represent such conversions from APE to KE.He subdivided APE and KE into mean and eddy forms.The eddy form is the deviation from the mean form and the mean form can convert into the eddy form by eddytransport of sensible heat from low to high latitudes. In his energy cycle process,the general circulation can be characterized as a conversion of zonal APE,which is generated by low-latitude heating and high-latitude cooling,to eddy APE,then to eddy kinetic energy and to zonal mean KE.The conversion between the two forms of energy involves horizontal and vertical transport of momentum and sensible heat.The dissipation and generation terms within the life cycle can barely be estimated directly,and they are only a few balance terms in the budget equations.

    This life cycle can be used to distinguish different contributions from different synoptic-scale and general circulation processes.It should give an intuitive judgment for evaluating the sources of different model performances.With regard to the Hadley cell,this largescale system involves a conversion from mean APE to mean KE(Dickinson,1969;Wu,1987),and for the barotropic atmosphere,the rate of transport between the eddy and mean KE is given by the product of the eddy transport if the momentum and gradient of mean angular rotation are both taken in the north-south direction.Concerning the baroclinic atmosphere(Stone, 1978),by the eddy transport of sensible heat from low to high latitudes,the mean APE converts into eddy APE and the eddy conversion from APE to KE is given by vertical motion and temperature within the latitudinal circle.Thus,it is a result of warm air rising and cold air sinking(Stein,1986;Huang and Vincent, 1998).In terms of the real atmosphere,the entire energy life cycle should display a complete conversion from APE to KE.

    Oort(1964,1983)used the classical Lorenz energy cycle concept to estimate the reservoirs of APE and KE,together with related sources and sinks for the Northern Hemisphere.Steinheimer et al.(2008) estimated the conversion of Lorenz theory both in gridscale processes and subgridscale processes obtained from parameterization schemes.Their results showed that the subgridscale processes contributed significantly to the Lorenz energy cycle and the total dissipation terms of the subgridscale were more intense than what all earlier gridscale estimates had indicated.

    Lorenz formulas can be introduced in the space domain,and we can further acquire the energy value in the mixed space and time domain,the so-called“Mixed Space-Time Domain”energy cycle.The eddy energy is subdivided into stationary and transient waves.The stationary form is the time mean and a result of diabatic and orographic forcing.The transient form is the departure from the time mean and a result of baroclinic instability of zonal mean flow(Simmons and Hoskins,1978,1980).

    Ulbrich and Speth(1991)expanded the Lorenz classical process to give a series of detailed formulations for the“Mixed Space-Time Domain”energy cycle(Arpe et al.,1986),and examined the role of stationary and transient waves within the atmospheric energy cycle with the ECMWF data for winter and summer.Their results showed that all terms of the energy cycle related to stationary waves reveal a predominance of the planetary scale while the transient waves are governed by synoptic-scale waves.

    In the present study,an approach based on the“Mixed Space-Time Domain”energy cycle is adopted.The formulations of the energy cycle are calculated,and the role of stationary and transient waves within the atmospheric energy cycle from the Global-Regional Assimilation and Prediction System (GRAPES)is firstly diagnosed and compared with NCEP Final Operational Global Analysis(FNL)data for July 2011.Three main kinds of considerable processes(planetary scale,barotropic conversion,and baroclinic conversion)are examined separately in order to investigate in detail the characteristics and features of energy cycle and its component contribution. The zonal-mean components of the energy cycle are investigated to diagnose the performance of numerical integrations.The forecast results at different lead times are used to explain the reason for deteriorating forecast performance.

    2.Method

    By following Ulbrich and Speth(1991)’s energy cycle framework(illustrated in Fig.1),all the energy reservoirs and conversions are calculated in this study with Ulbrich and Speth’s detailed formulations. Zonal available potential energy(AZ)and zonal kinetic energy(KZ)are zonal APE and KE.Stationary eddy available potential energy(ASE)and transient eddy available potential energy(ATE)are stationary and transient eddy APE,and the stationary and transient eddy KE is named KSE and KTE,respectively. Besides energy reservoirs,some terms of the energy cycle,such as CZ(conversion from AZ to KZ),CAS (conversion from AZ to ASE),CAT(conversion from AZ to ATE),CES(conversion from ASE to KSE), CET(conversion from ATE to KTE),CKS(conversion from KZ to KSE),CKT(conversion from KZ to KTE),and the nonlinear conversions between stationary and transient terms(CATE and CKTE)all play important roles in the global energy cycle process and can be regarded as an important benchmark to estimate the basic model performance.

    3.Data

    Fig.1.Diagram of the global atmospheric energy cycle in the Mixed Space-Time Domain.Arrows indicate orientation of conversions corresponding to the definitions of parameters(from Ulbrich and Speth,1991).Note:G denotes generation of potential energy;D denotes dissipation of kinetic energy.

    The GRAPES model daily forecast data for July 2011 with four different lead times of 24,72,120,and 168 h are used.The model resolution is selected as 0.5°and it is initialized with global 1200 UTC analyses of the NCEP FNL data.Archived GRAPES model data consist of geopotential height,temperature,specific humidity,zonal and meridional wind,and vertical wind for 29 pressure levels(1000,962.5,925,887.5, 850,800,750,700,650,600,550,500,450,400,350, 300,275,250,225,200,175,150,125,100,70,50, 30,20,and 10 hPa).NCEP FNL data are selected as the analysis data for comparison and are interpolated from 1.0°to 0.5°grid,and vertically from 26 to 29 pressure levels.NCEP global forecast(not analysis)data are also selected for the same time period and the same horizontal and vertical resolutions for comparison.

    4.Results

    4.1 The global energy cycle for July 2011

    As a preliminary work for the energy cycle analysis,it is necessary to evaluate the performance of the GRAPES model and provide an intuitive impression of the quality of the FNL data,which serve as the comparative observation.In the classical treatment of statistical verification for different model forecasts, root-mean-square-error(RMSE)is a predominant index used as a measure of model forecast performance. Temperature and wind velocity are the main variables involved in the calculation of energy cycle and conversion terms.Figure 2 depicts the global averaged RMSE for temperature and wind velocity against the FNL analysis data and NCEP BUFR(Binary Universal Form of the Representation of meteorological data) sounding data.The latter dataset is chosen as the observation,which is only distributed between 850 and 50 hPa.Figure 2a shows that the maxima of temperature RMSE are mostly concentrated within 200 hPa,and with increasing lead time,RMSE grows from 1.59 K(24 h)to 3.80 K(168 h).Meanwhile,in Fig. 2c,the maxima of wind velocity RMSE are located within 250 hPa,and with increasing lead time,RMSE rises from 6.14 m s-1(24 h)to 16.16 m s-1(168 h). The RMSEs of GRAPES model outputs relative to the FNL data in Fig.2b and 2d present a similar pattern to that in Figs.2a and 2c.Temperature RMSE at different lead times in Fig.2c is consistently smaller

    than that in Fig.2a(e.g.,at 200 hPa,1.09 K for 24 h and 3.53 K for 168 h).Moreover,wind velocity RMSE is smaller in Fig.2d than in Fig.2b at 120-h lead time and increases to 17.00 m s-1(168 h)at 250 hPa.This is possibly due to the GRAPES model being initiated with the FNL data.With respect to upper levels(above 50 hPa),where there are no comparative observations,large temperature and wind velocity RMSE maxima are observed to intensify with lead time.This will lead to strengthening potential energy and KE.

    Fig.2.Global averaged RMSE for GRAPES produced (a,b)temperature(K)against(a,c)NCEP BUFR observation and(b,d)FNL analysis in July with different lead times.

    Following Ulbrich and Speth’s work(1991),we have calculated the globally averaged energy cycle of the GRAPES model with different lead times and the results are shown in Fig.3.It is confirmed that the GRAPES model has the capability to reproduce the main features of the global energy cycle as compared with NCEP analysis data.AZ is converted into ASE and ATE;the stationary waves cannot be neglected compared with the transient ones,and ASE and ATE have about the same value.The nonlinear conversion between the two eddy available potential energy terms, i.e.,CATE,plays an important role in the global energy cycle and is directed from the stationary to the

    transient reservoir of APE.It can be deduced that the damping of stationary temperature by horizontal transient fluxes of sensible heat is an important process in the global general circulation.

    Fig.3.The global atmospheric energy cycle of GRAPES in July 2011.Various energy components(in boxes)are in J m-2,while conversions between the components are in W m-2.Numbers at the top indicate values based on NCEP FNL data,and 24,72,120,and 168 h are the different lead times from 1 to 7 days.The dotted frame refers to the planetary-scale(Hadley cell)process,the dashed frame refers to the barotropic process,and the dot-dashed frame refers to the baroclinic process.

    With increased forecast lead time,AZ becomes larger,which reflects the meridional temperature gradient between high and low latitudes becoming steeper. This enhances zonal baroclinic processes, which increases the meridional heat flux and enhances the conversion from AZ to eddy potential energy,especially the conversion of the transient term,CAT.The stationary conversion of APE,i.e.,CAS,is just 1/6 of the magnitude of the transient conversion of CAT.

    The zonal KE(i.e.,KZ)has a similar value to the sum of stationary and transient eddy KE(KSE and KTE),and the transient term is three times larger than the stationary term,while there is almost no global net conversion between the two eddy KE terms. Barotropic stationary and transient conversions(CKS and CKT)are directed from eddy KE to zonal KE. The zonal planetary-scale conversion CZ(from AZ to KZ)is about 1.5 times larger in GRAPES than in the NCEP analysis.

    4.2 Planetary-scale processes

    The globally averaged energy values are presented in Fig.3,which gives us an intuitive impression of how the global energy cycle works.It is also a useful diagnosis tool to analyze the source of the bias between different model forecasts.However,globally averaged values can barely depict the contributions and interrelations of special energy terms.More needs to be done to reveal detailed features.We further separate the framework of the energy cycle into three major processes:planetary-scale,barotropic,and baroclinic atmospheric processes.

    With regard to planetary-scale processes,largescale systems such as the Hadley cell provide a conversion from zonal APE(AZ)to zonal KE(KZ).The AZ and KZ increase with the increasing lead time.Figure 4 shows the zonal mean temperature distribution of NCEP FNL analysis and the bias of the GRAPES forecast relative to FNL data at different lead times.As expected,the maximum overestimation is located at low levels near the equator and at high levels near the poles,and the overestimation becomes stronger with increasing lead time.This indicates that the strong warm bias in the low latitudes increases the meridional temperature gradient,which leads to stronger AZ.Concerning the CZ term,which represents the conversion between AZ and KZ,this reflects the zonal features of the vertical wind and temperature.When warm air rises and cold air sinks,CZ>0 and shows Hadley cell features.In contrast,when CZ<0,it shows the features of the Ferrel cell.In Fig.3,it presents positive values of CZ,so it reflects a Hadley cell feature in planetary-scale processes.

    4.3 Barotropic conversion

    In this section,the barotropic conversion and KE are diagnosed.As discussed in Section 4.1,zonal KE (KZ)is the sum of the stationary and transient eddy KE,while KSE is about 1/3 of KTE and there is almost no net global conversion between the two eddy terms.The zonal distribution of KTE(Fig.5)shows that the maxima in the Southern Hemisphere are displaced northward at the location of the zonal mean jet (50°S)and the maxima in the Northern Hemisphere

    are displaced southward in the opposite latitudes. With increasing lead time,the strong western flow decreases in the upper troposphere(near 300 hPa)and stratosphere(near the top of the model),which leads to a decreasing of the global averaged energy value.

    Fig.4.(a)Zonal averaged temperature distribution of FNL and(b)-(e)biases of GRAPES relative to FNL with different lead times(1-7 days,i.e.,D+1 to D+7).

    Fig.5.Transient eddy kinetic energy(KTE)in July(J m-2Pa-1).(a)FNL data,and(b)24-,(c)72-,(d)120-, and(e)168-h GRAPES forecast.

    The stationary eddy KE(KSE)displays a very different pattern to the transient one(Fig.6).The comparatively weak maxima are presented at the locations of 30°S,15°N,and 45°N,separately,and are associated with the zonal mean jets in both hemispheres. Another maximum appears in high latitudes of the Southern Hemisphere,which is related to the polar jet.With increasing leading time,the four weaker maxima decrease significantly.A band of high energy values presents in the Southern Hemisphere over 100 hPa due to the presence of a large zonal wind jet at the upper 30-hPa level,which causes a strengthening of KSE.

    As mentioned before,the nonlinear conversion be tween KSE and KTE is directed from KSE to KTE, and it is difficult to determine the symbol of the global averaged value because it is approximately zero.However,the zonal mean distribution of transient eddy momentum against or along the direction of the stationary eddy momentum gradient can be explained. The zonal mean distribution of the nonlinear conversion CKTE shows a pair of strong maxima in Fig.7. Comparing the location with respect to the zonal mean

    contribution of KSE and KTE,we find that the maximum of the nonlinear conversion is displaced near the location of the maximum KTE and minimum KSE. This means that the strong nonlinear conversion can be considered as a physical mechanism forcing the atmosphere to a more uniform state.Figure 7 shows that as the lead time increases,the maximum nonlinear conversion value becomes weaker,which indicates a weakening of the local contribution.

    Fig.6.As in Fig.5,but for the stationary eddy kinetic energy(KSE)in July(J m-2Pa-1).

    Fig.7.Distributions of the conversion term CKTE(from KTE to KSE;10-6W m-2Pa-1)in July.(a)FNL data, and(b)24-,(c)72-,(d)120-,and(e)168-h GRAPES forecast.

    CKTE is a very important nonlinear conversion term and its positive or negative value indicates the direction of transfer between transient waves and stationary waves.Moreover,the strength of KSE can be used as a basis to determine the jet intensity changes.When CKTE is positive,heat transport is transferred from transient waves to stationary waves, which strengths KSE.Conversely,if CKTE is negative,heat transport is directed into transient waves, resulting in a weakening of KSE.The result obtained by Ulbrich and Speth(1991)showed that a strengthening of KSE corresponded to local jet maxima while a weakening prevailed over the jet entrance regions.

    Figure 7 shows the conversion term CKTE from FNL analysis and GRAPES forecast,while the CKTE in the Southern Hemisphere is depicted in Fig.8. It is seen that in the FNL analysis,positive values are mainly found downstream of the jet maxima over Southwest Pacific,while negative values are distributed in the jet entrance regions over Southeast Atlantic where wind velocity is relatively weak.The distribution of positive and negative local values and jet regions verifies Ulbrich’s conclusion.Compared with the FNL analysis,the strength of the local jet in the GRAPES forecast becomes gradually weaker with increased lead time,while the corresponding CKTE decreases significantly,together leading to the dissipation of KSE.Until 7 days,the positive values over the local jet maxima show a clear decreasing trend,while the corresponding jet is significantly weaker than in the FNL analysis.It can be concluded that CKTE can be used to determine stationary wave changes,and act as an indicator of the jet intensity features.

    4.4 Baroclinic conversion

    For the baroclinic conversion process,the global averaged values reported in Fig.3 indicate how this process is working.With increased lead time,the conversion from zonal APE to stationary eddy APE remains constant.As a result of the meridional temperature gradient increasing,meridional transport of sensible heat enhances,which causes an increase in the conversion CAT,transferring energy from zonal to transient eddy APE.

    In order to compare details of the process associated with baroclinic conversion,the zonal contributions of ASE and ATE are plotted in Fig.9.Larger ASE values appear in both the Southern and Northern Hemispheres.In the Northern Hemisphere,the largest ASE is found below 850 hPa between 20°and 60°N,and the second largest ASE is seen between 400 and 200 hPa around 30°N.However,in the Southern Hemisphere,the largest ASE occurs between 850 and 700 hPa and further south of 60°S,with the maximum value greater than 30 J m-2Pa-1.In contrast to ASE,the zonal mean distribution of ATE has a similar structure to KTE(Fig.5).The maximum value of transient eddy potential energy ATE is only about half that of ASE.With increased lead time,the maximum value of ASE increases and ATE decreases slightly.

    In terms of the zonal mean distribution of ASE, Fig.10 shows that the maxima of ASE are associated with the linear conversion term CAS.The zonal contribution of stationary conversion is negative over 200-400 hPa at 30°N,corresponding to similar conversion rates of ASE maximizing at a similar location. The same coincidence can be seen in high latitudes of the Southern Hemisphere.The zonal mean contribution to the conversion from zonal to transient eddy APE,i.e.,CAT,is determined by the meridional temperature gradient.Comparing CAT(Fig.10)with ATE(Fig.9)reveals that a coincidence in the maximum energy value with conversion terms is apparent.The conversion value increases over low levels in the Southern Hemisphere with increasing lead time. This is determined by the decreasing of vertical stability and,simultaneously,the decreasing of transient sensible heat transport.The coincidence in the distribution of CAS is found to be very similar to that observed for the transient term with respect to midtropospheric values.This indicates that the baroclinic

    process associated with the stationary wave is comparable to that associated with the transient wave.

    Fig.8.Wind velocity(m s-1)at 250 hPa in July 2011 in the Southern Hemisphere.Shaded areas depict the local contributions to CKTE greater than 50×10-6W m-2Pa-1(dark)and lower than-50×10-6W m-2Pa-1(light).(a) FNL data,and(b)24-,(c)72-,(d)120-,and(e)168-h GRAPES forecast.

    Fig.9.(a1-e1)Stationary eddy available potential energy(ASE)and(a2-e2)transient eddy available energy(ATE)(J m-2Pa-1)in July.(a1-a2)FNL data,and(b1-b2)24-,(c1-c2)72-,(d1-d2)120-,and(e1-e2)168-h forecast.

    Fig.10.As in Fig.9,but for the conversion terms CAS(a1-e1)and CAT(a2-e2)(10-6W m-2Pa-1)in July.

    The nonlinear conversion CATE is directed from ASE to ATE.Local zonal mean values of this nonlinear conversion are depicted in Fig.11.The maximum contribution to the global averaged energy is located at lower levels over Antarctica,where the maximum value of ASE presents.It is considered that this kind of energy conversion does not originate from the zonal reservoir of AZ,but from the stationary eddy APE’s nonlinear interactions.The maximum ASE above the troposphere at 30°N is not associated with the strong nonlinear conversions.The negative values of CATE is seen to intensify over Antarctica with increased

    lead time. The nonlinear conversion is thought as the horizontal heat transport by transient waves(Lau, 1982)and a negative value means that the heat transport is directly from warm to cold regions.

    Fig.11.As in Fig.7,but for CATE(ATE to ASE;10-6W m-2Pa-1)in July.

    As the maximum zonal mean CATE appears at 700 hPa in the Southern Hemisphere,we choose to analyze the nonlinear conversion CATE at 700-hPa in the Southern Hemisphere(Fig.12).With regard to the FNL analysis,some negative and positive values of CATE are located around Antarctica,coinciding with regions of relatively warm and cold air.We superimpose temperature on CATE,and see that the distribution coincides with the horizontal heat transport by the transient waves.When heat is transported against the local temperature gradient by transient waves,meaning that heat transport is directed into the warm region,the CATE value is positive.In contrast,a negative value indicates that the heat transport is directed into the cold region.Compared with FNL data,the GRAPES forecast data show a larger temperature gradient and increasing negative values.

    4.5 Comparison with NCEP-GFS

    According to the aforementioned analysis of the atmospheric energy cycle,the main characteristics of the energy cycle is now used to aid in traditional verification of model performance.Real forecast data with the same resolution as that of FNL data from the NCEP Global Forecast Systems(NCEP-GFS)are used to compare the features of the energy cycle with the GRAPES model.The same 31 days of analysis are carried out and monthly averages are considered.

    Firstly,the anomaly correlation coefficient(ACC) is calculated for geopotential height at 500 hPa. We use wavelet analysis and our verification procedure include a separation of the waves with different zonal wave numbers.The long waves of the planetary scale(zonal wavenumbers 0-3),synoptic scale waves(zonal wavenumbers 4-9),and small-scale waves (zonal wavenumbers 10-20)are distinguished separately.We assume that fast moving waves are synoptic or small scale,and slowly moving ones are predominantly planetary scale(Kung,1988).

    In Fig. 13,the separation into the different wavenumber groups is shown in order to compare the performance of NCEP forecast and GRAPES forecast. With respect to NCEP forecast results,the ACC of GRAPES falls behind with different lead times.The bias at 168 h is approximately 0.03 and all verification indices are significant at the 95%confidence level. The planetary-scale contribution is roughly similar to the aggregate performance. However,the synoptic scale and shorter scale represent some differences:for the synoptic scale(wavenumbers 4-9),the confidence interval becomes larger at 168 h,meaning that the standard deviation of ACC bias also becomes larger. Therefore,the bias can barely satisfy the 95%confidence level.For the shorter scale(wavenumbers 10-20),the bias of model performance becomes more statistically insignificant after 48 h.We conclude that significant bias is reflected on the large and meso scales, and is roughly equal on the smaller scale.

    Fig.12.Temperature(K)at 700 hPa in July 2011 in the Southern Hemisphere.Shaded areas depict local contributions to CATE greater than 50×10-6W m-2Pa-1(dark)and lower than-50×10-6W m-2Pa-1(light).(a)FNL data,and (b)24-,(c)72-,(d)120-,and(e)168-h forecast.

    Fig.13.Monthly mean of 500-hPa geopotential height ACC for GRAPES and NCEP-GFS with different lead times of 1-7 days.(a)Total statistical score,(b)planetary scale(wavenumbers 0-3),(c)synoptic scale(wavenumbers 4-9), and(d)short scale(wavenumbers 10-20).Note:“w.r.t.”means“with regard to.”

    Figure 14 shows the energy values for the global integrals of NCEP-GFS with different lead times.In contrast to the GRAPES model,it is generally accepted that most energy values calculated by baroclinic and barotropic processes are roughly similar in July.The primary difference lies in the conversion of zonal APE to zonal KE,which is a result of the planetary-scale processes.AZ of NCEP-GFS is larger than in GRAPES,which means that NCEP-GFS has a steeper meridional temperature gradient between the high and low latitudes.The conversion term CZ of NCEP-GFS is smaller than in GRAPES,which is mainly due to the simulated weak vertical wind.Due to the hydrostatic assumption,NCEP-GFS simulates weaker vertical air motion than GRAPES,leading to the nonlinear conversion CZ becoming smaller.All of these features will seriously affect the Hadley cell simulation.KZ of GRAPES is larger than that of NCEP-GFS,mainly because of the zonal mean wind amplitude being too large.It should be noted at this point that the difference in model performance between GRAPES and NCEP-GFS is mainly reflected in the large-scale process,and such a conclusion is consistent with the results obtained via traditional verification in Fig.13.Furthermore,detailed characteristics should be investigated and compared with a full analysis of unique energy cycle processes.

    Fig.14.The global atmospheric energy cycle in the Mixed Space-Time Domain for NCEP-GFS in July 2011.Various energy components(in boxes)are in J m-2,while conversions between the components are in W m-2.Numbers at the top indicate values based on NCEP FNL data,and 24,72,120,and 168 h are the different lead times from 1 to 7 days. The dashed box indicates the large difference of NCEP-GPS with the GRAPES model.

    5.Summary

    An investigation on the atmospheric energy cycle was carried out by using GRAPES forecast data,with a focus on the role of model performance at different lead times.Three main atmospheric-scale processes were diagnosed separately and the characteristics of stationary and transient eddy energy terms and their conversions were calculated and presented.The results show that:

    (1)The GRAPES model has the capability to reproduce the main characteristics of the global energy cycle as compared with NCEP FNL analysis data.

    (2)ASE and ATE have approximately the same values and the nonlinear conversion is directed from the stationary to the transient.

    (3)Barotropic conversions(CKS and CKT)are directed from eddy kinetic energy to zonal kinetic energy,and the kinetic energy of stationary eddy is about 1/3 of that of the transient eddy.

    (4)With increasing forecast lead time,AZ becomes larger and the zonal conversion CZ(from AZ to KZ)in GRAPES is around 1.5 times larger than in the NCEP analysis.

    (5)In contrast with traditional verification,the energy cycle diagnosis can help to derive statistical scores and identify the source of model differences.

    This study on the diagnosis of the energy cycle is at a preliminary stage in terms of operational application.There are some areas of work that need to be improved and some problems that should be noted. For example,only short-term integrations of forecast data were collected for analysis in order to acquire stable performances of energy values.Low-resolution and long-term integrations should be run for“seasonal forecasts.”Moreover,we have recognized that diagnosis of the energy cycle can be achieved with a similar conclusion to traditional verification methods,but how to combine traditional verification methods with this energy cycle diagnosis for model evaluation remains to be examined.These problems should be discussed in future investigations.

    Arpe,K.,C.Brancovic,E.Oriol,et al.,1986:Variability in time and space of energetics from a long series of atmospheric data produced by ECMWF.Beitr. Phys.Atmos.,59,321-355.

    Dickinson,R.E.,1969:Theory of planetary wave-zonal flow interaction.J.Atmos.Sci.,26,73-81.

    Gao Hui,Chen Longxun,He Jinhai,et al.,2006:Characteristics of zonal propagation of atmospheric kinetic energy at equatorial region in Asia.Acta Meteor. Sinica,20,86-94.

    Huang,H.-J.,and D.G.Vincent,1988:Daily spectral energy conversions of the global circulation during 10-27 January 1979.Tellus,40A,37M9.

    Kung,E.C.,1988:Spectral energetics of the general circulation and time spectra of transient waves during the FGGE year.J.Climate,1,5-19.

    Lau,N.C.,and A.H.Oort,1982:A comparative study of observed Northern Hemisphere circulation statistics based on GFDL and NMC analyses.Part II: Transient eddy statistics and the energy cycle.Mon. Wea.Rev.,110,889-906.

    Li Qingquan and Zhu Qiangen,1995:Analysis on the source and sink of kinetic energy of atmospheric 30-60-day period oscillation and the probable causes. Acta Meteor.Sinica,9,420-431.

    Lorenz,E.N.,1955:Available potential energy and the maintenance of the general circulation.Tellus,7, 157-167.

    Luo Zhexian,1994:Effect of energy dispersion on the structure and motion of tropical cyclone.Acta Meteor.Sinica,8,51-59.

    Oort,A.H.,1964:On estimates of the atmospheric energy cycle.Mon.Wea.Rev.,92,483-493.

    —-,1983:Global Atmospheric Circulation Statistics 1958-1973,Volume 14,NOAA Professional Paper. U.S.Department of Commerce/NOAA,180.

    Simmons,A.J.,and B.J.Hoskins,1978:The life cycles of some nonlinear baroclinic waves.J.Atmos.Sci., 35,414-432.

    —-,and—-,1980:Barotropic influences on the growth and decay of nonlinear baroclinic waves.J.Atmos. Sci.,37,1679-1684.

    Stein, C., 1986: Das zeitliche Zusammenwirken baroklinerEnergieumwandlungen durch grogrfiumige Wellen in derAtmosph/ire. Diplomarbeit, Institut fiir Geophysik und Meteorologie der Universit/it zu K61n.

    Steinheimer,M.,M.Hantel,and P.Bechtold,2008:Convection in Lorenz’s global energy cycle with the ECMWF model.Tellus,60A,1001-1022.

    Stone,P.H.,1978:Baroclinic adjustment.J.Atmos. Sci.,35,561-571.

    Storch,J.S.,C.Eden,and I.Fast,2012:An estimate of the Lorenz energy cycle for the world ocean based on the STORM/NCEP simulation.J.Phys.Oceanogr., 42,2185-2205.

    Ulbrich,U.,and P.Speth,1991:The global energy cycle of stationary and transient atmospheric waves: Results from ECMWF analyses.Meteor.Atmos. Phys.,45,125-138.

    Wu Rongsheng,1987: Energy,energy flux and Lagrangian of Rossby wave. Acta Meteor. Sinica, 1,143-150.

    :Zhao Bin and Zhang Bo,2014:Diagnostic study of global energy cycle of the GRAPES global model in the mixed space-time domain.J.Meteor.Res.,28(4),592-606,

    10.1007/s13351-014-3072-0.

    (Received October 16,2013;in final form June 17,2014)

    The results show that the GRAPES model has the capability to reproduce the main features of the global energy cycle as compared with the NCEP analysis.Zonal available potential energy(AZ)is converted into stationary eddy available potential energy(ASE)and transient eddy available potential energy(ATE),and ASE and ATE have similar values.The nonlinear conversion between the two eddy energy terms is directed from the stationary to the transient.AZ becomes larger with increased forecast lead time,reflecting an enhancement of the meridional temperature gradient,which strengthens the zonal baroclinic processes and makes the conversion from AZ to eddy potential energy larger,especially for CAT(conversion from AZ to ATE).The zonal kinetic energy(KZ)has a similar value to the sum of the stationary and transient eddy kinetic energy.Barotropic conversions are directed from eddy to zonal kinetic energy.The zonal conversion from AZ to KZ in GRAPES is around 1.5 times larger than in the NCEP analysis.The contributions of zonal energy cycle components show that transient eddy kinetic energy(KTE)is associated with the Southern Hemisphere subtropical jet and the conversion from KZ to KTE reduces in the upper tropopause near 30°S.The nonlinear barotropic conversion between stationary and transient kinetic energy terms(CKTE) is reduced predominantly by the weaker KTE.

    猜你喜歡
    張博
    過(guò) 橋
    對(duì)善舉的尊重
    項(xiàng)目學(xué)習(xí)背景下高中信息技術(shù)課堂教學(xué)有效實(shí)施的探索
    師道·教研(2021年8期)2021-09-22 12:12:03
    集成創(chuàng)新發(fā)揮優(yōu)勢(shì) 自主研發(fā)服務(wù)造紙
    造紙信息(2019年7期)2019-09-10 11:33:18
    心里住著彼得·潘
    不孝之子
    一道高考試題的四次拓展
    張博:一個(gè)任性Girl的折騰史
    不孝之子
    基因組編輯技術(shù)研究前沿與挑戰(zhàn)
    遺傳(2015年10期)2015-01-03 06:31:58
    日本免费一区二区三区高清不卡| 一个人免费在线观看电影| 午夜老司机福利剧场| 91精品国产九色| av专区在线播放| 成年女人看的毛片在线观看| 亚洲自偷自拍三级| 国产老妇伦熟女老妇高清| 国产单亲对白刺激| a级一级毛片免费在线观看| 国产高清不卡午夜福利| 一区福利在线观看| 好男人视频免费观看在线| 51国产日韩欧美| 中文字幕久久专区| 人人妻人人澡人人爽人人夜夜 | 欧洲精品卡2卡3卡4卡5卡区| 岛国毛片在线播放| 男女下面进入的视频免费午夜| 欧美性感艳星| 非洲黑人性xxxx精品又粗又长| 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 亚洲18禁久久av| 高清毛片免费观看视频网站| 男女那种视频在线观看| 我的老师免费观看完整版| 色吧在线观看| 久久九九热精品免费| 国产乱人偷精品视频| 大香蕉久久网| 99久久精品热视频| 高清毛片免费观看视频网站| 国产色婷婷99| 亚洲av中文av极速乱| 国产一区二区激情短视频| 久久久国产成人免费| 久久久久久久亚洲中文字幕| 国产熟女欧美一区二区| 亚洲av成人av| 1024手机看黄色片| 天堂av国产一区二区熟女人妻| 色哟哟哟哟哟哟| 亚洲一级一片aⅴ在线观看| 国模一区二区三区四区视频| 亚洲人成网站高清观看| 久久久久久九九精品二区国产| 亚洲国产欧洲综合997久久,| 亚洲最大成人手机在线| 免费人成在线观看视频色| 身体一侧抽搐| 亚洲无线在线观看| 人妻系列 视频| 午夜免费男女啪啪视频观看| 只有这里有精品99| 亚洲av熟女| 女的被弄到高潮叫床怎么办| 老熟妇乱子伦视频在线观看| 爱豆传媒免费全集在线观看| 乱系列少妇在线播放| 99热精品在线国产| 99国产极品粉嫩在线观看| 国产视频首页在线观看| 青春草国产在线视频 | 亚洲精品色激情综合| 性色avwww在线观看| 99久久成人亚洲精品观看| av福利片在线观看| 九色成人免费人妻av| 1024手机看黄色片| 赤兔流量卡办理| 日韩一本色道免费dvd| 精品久久久久久久久av| 久久久国产成人精品二区| 久久久久久久久久成人| 色综合亚洲欧美另类图片| 成人无遮挡网站| 中文字幕人妻熟人妻熟丝袜美| 在线观看免费视频日本深夜| 99久久中文字幕三级久久日本| 久久久久久久午夜电影| h日本视频在线播放| 精品一区二区三区视频在线| 精品欧美国产一区二区三| 搡老妇女老女人老熟妇| 深爱激情五月婷婷| 日韩高清综合在线| 春色校园在线视频观看| 日本免费a在线| 男女视频在线观看网站免费| 深爱激情五月婷婷| 美女内射精品一级片tv| 日韩 亚洲 欧美在线| 欧美3d第一页| 亚洲,欧美,日韩| 97人妻精品一区二区三区麻豆| 嫩草影院新地址| 亚洲人成网站在线播| 国产黄片美女视频| 美女黄网站色视频| 国产视频内射| 国产精品一区二区三区四区免费观看| 国产精品,欧美在线| 亚洲av免费高清在线观看| 欧美一级a爱片免费观看看| 成人毛片60女人毛片免费| 亚洲欧洲日产国产| 久久99热6这里只有精品| 国产伦在线观看视频一区| 五月伊人婷婷丁香| 国产成人影院久久av| 观看免费一级毛片| 精品人妻偷拍中文字幕| 嫩草影院精品99| 国产一级毛片在线| av.在线天堂| 国产亚洲91精品色在线| 久久精品久久久久久久性| 免费看光身美女| 26uuu在线亚洲综合色| 六月丁香七月| 欧美潮喷喷水| 国产一区二区在线观看日韩| 午夜福利视频1000在线观看| 亚洲无线观看免费| 美女国产视频在线观看| 国产精品一区www在线观看| 午夜免费男女啪啪视频观看| 久久久久久久午夜电影| 国产精品伦人一区二区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧洲国产日韩| 乱系列少妇在线播放| a级一级毛片免费在线观看| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 亚洲av不卡在线观看| 波野结衣二区三区在线| 成人特级av手机在线观看| 亚洲成a人片在线一区二区| 亚洲在久久综合| 成人国产麻豆网| 亚洲欧美清纯卡通| 小蜜桃在线观看免费完整版高清| 91aial.com中文字幕在线观看| 成人午夜高清在线视频| 免费人成在线观看视频色| 欧美日本亚洲视频在线播放| 男人的好看免费观看在线视频| 国产探花在线观看一区二区| 国产三级在线视频| 男插女下体视频免费在线播放| 国产精品电影一区二区三区| 国产高清激情床上av| 亚洲四区av| 亚洲人成网站在线观看播放| 国内揄拍国产精品人妻在线| 成人性生交大片免费视频hd| 大香蕉久久网| 亚洲精品国产成人久久av| 99久国产av精品| 人体艺术视频欧美日本| 一级二级三级毛片免费看| 亚洲av电影不卡..在线观看| 男的添女的下面高潮视频| 国产精品精品国产色婷婷| 欧洲精品卡2卡3卡4卡5卡区| 99国产极品粉嫩在线观看| 免费搜索国产男女视频| 黄色欧美视频在线观看| 26uuu在线亚洲综合色| 给我免费播放毛片高清在线观看| 久久精品国产亚洲av涩爱 | 免费人成在线观看视频色| 亚洲欧美精品专区久久| 好男人在线观看高清免费视频| 国产午夜精品一二区理论片| 亚洲综合色惰| 麻豆精品久久久久久蜜桃| 少妇高潮的动态图| a级毛片a级免费在线| 亚洲av电影不卡..在线观看| 亚洲精品自拍成人| 日韩欧美精品免费久久| 国产精品99久久久久久久久| 又粗又硬又长又爽又黄的视频 | 身体一侧抽搐| 亚洲aⅴ乱码一区二区在线播放| 男女那种视频在线观看| 日韩中字成人| 久久精品夜色国产| 美女脱内裤让男人舔精品视频 | 国产成人精品一,二区 | 成人美女网站在线观看视频| 神马国产精品三级电影在线观看| 性插视频无遮挡在线免费观看| 中文字幕久久专区| 国产不卡一卡二| 99久国产av精品国产电影| 亚洲性久久影院| 日韩欧美一区二区三区在线观看| 久久99热6这里只有精品| 欧美日韩综合久久久久久| 国产一区二区在线av高清观看| 啦啦啦啦在线视频资源| 不卡一级毛片| 亚洲av第一区精品v没综合| 天堂影院成人在线观看| 欧美日韩在线观看h| 男人的好看免费观看在线视频| 国产亚洲精品av在线| 99热网站在线观看| 国产亚洲欧美98| 亚洲人成网站在线播放欧美日韩| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 国产伦精品一区二区三区视频9| 精品一区二区三区视频在线| 亚洲国产日韩欧美精品在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲五月天丁香| 亚洲av中文字字幕乱码综合| 深夜精品福利| 亚洲第一电影网av| 成人三级黄色视频| 性色avwww在线观看| 乱码一卡2卡4卡精品| 波多野结衣巨乳人妻| 国产日本99.免费观看| 国产成人精品一,二区 | 狂野欧美激情性xxxx在线观看| 国内精品宾馆在线| 99久久成人亚洲精品观看| 国产精品一区二区三区四区免费观看| 国产精品久久久久久亚洲av鲁大| 中文字幕熟女人妻在线| 欧美色欧美亚洲另类二区| 亚洲成av人片在线播放无| 国产精品综合久久久久久久免费| 国产亚洲91精品色在线| 亚洲av免费高清在线观看| 午夜精品一区二区三区免费看| 精品熟女少妇av免费看| 99久久中文字幕三级久久日本| 久久久久九九精品影院| 国产精品99久久久久久久久| 久久久久久久午夜电影| 午夜精品一区二区三区免费看| 黄片wwwwww| 秋霞在线观看毛片| 不卡一级毛片| 男插女下体视频免费在线播放| 日本熟妇午夜| 青青草视频在线视频观看| 亚洲三级黄色毛片| 色吧在线观看| 在线播放国产精品三级| 国产一区二区三区av在线 | 中文欧美无线码| 久久久国产成人免费| h日本视频在线播放| 黄色视频,在线免费观看| 成人综合一区亚洲| 丰满人妻一区二区三区视频av| 免费看日本二区| 国产乱人视频| 国产精品一区www在线观看| 小说图片视频综合网站| 久久6这里有精品| 日韩成人伦理影院| a级毛片免费高清观看在线播放| 99久久成人亚洲精品观看| 国产午夜精品久久久久久一区二区三区| 麻豆国产av国片精品| 久久久久久久久久久免费av| 欧美xxxx黑人xx丫x性爽| 中文字幕久久专区| 亚洲一级一片aⅴ在线观看| 老司机福利观看| 又粗又爽又猛毛片免费看| 99riav亚洲国产免费| 不卡视频在线观看欧美| 在线观看美女被高潮喷水网站| 99久久精品一区二区三区| av免费在线看不卡| 成人鲁丝片一二三区免费| 22中文网久久字幕| 久久久久网色| 国产一级毛片七仙女欲春2| 丝袜美腿在线中文| 免费看光身美女| 舔av片在线| 国产探花在线观看一区二区| 亚洲18禁久久av| 高清在线视频一区二区三区 | 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类| 亚洲最大成人中文| 此物有八面人人有两片| 亚洲最大成人中文| 神马国产精品三级电影在线观看| av视频在线观看入口| 性色avwww在线观看| 国产成人精品一,二区 | or卡值多少钱| 欧美另类亚洲清纯唯美| 成人亚洲欧美一区二区av| 九九久久精品国产亚洲av麻豆| 日本一二三区视频观看| 免费一级毛片在线播放高清视频| 日韩欧美精品v在线| 欧美性猛交╳xxx乱大交人| 国产精品一二三区在线看| 18+在线观看网站| 欧美bdsm另类| 性色avwww在线观看| 精品国内亚洲2022精品成人| 欧美一区二区精品小视频在线| 最近视频中文字幕2019在线8| 亚洲一区高清亚洲精品| 免费看a级黄色片| 三级毛片av免费| 成人毛片60女人毛片免费| 免费观看的影片在线观看| 日本成人三级电影网站| 男女视频在线观看网站免费| 亚洲av二区三区四区| 日本爱情动作片www.在线观看| 日日摸夜夜添夜夜添av毛片| 中国国产av一级| 亚洲人成网站在线观看播放| 国产视频首页在线观看| 日韩成人伦理影院| 在线观看美女被高潮喷水网站| 亚洲成av人片在线播放无| 一级毛片电影观看 | 国产免费男女视频| 一个人免费在线观看电影| 中国美白少妇内射xxxbb| 亚洲av成人av| 中文字幕免费在线视频6| 亚洲中文字幕一区二区三区有码在线看| 国产视频首页在线观看| 蜜桃久久精品国产亚洲av| 成人毛片a级毛片在线播放| 亚洲熟妇中文字幕五十中出| 国产精品嫩草影院av在线观看| 少妇的逼水好多| 黄片无遮挡物在线观看| 国产成年人精品一区二区| 国产一级毛片在线| 一进一出抽搐gif免费好疼| 一级二级三级毛片免费看| 国产在视频线在精品| 99久久久亚洲精品蜜臀av| 一级毛片久久久久久久久女| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐gif免费好疼| 少妇猛男粗大的猛烈进出视频 | 午夜福利在线观看免费完整高清在 | 人妻少妇偷人精品九色| 变态另类成人亚洲欧美熟女| 综合色av麻豆| 尤物成人国产欧美一区二区三区| 男女啪啪激烈高潮av片| 一区二区三区四区激情视频 | 蜜臀久久99精品久久宅男| 嫩草影院新地址| 大又大粗又爽又黄少妇毛片口| 在线观看一区二区三区| 久久韩国三级中文字幕| 国产精品不卡视频一区二区| 国产蜜桃级精品一区二区三区| 三级毛片av免费| 国产免费男女视频| 九色成人免费人妻av| 麻豆精品久久久久久蜜桃| 免费观看在线日韩| 亚洲18禁久久av| 一级黄片播放器| 寂寞人妻少妇视频99o| 亚洲美女视频黄频| 人妻系列 视频| 国产伦一二天堂av在线观看| 日韩欧美 国产精品| 少妇猛男粗大的猛烈进出视频 | 少妇人妻精品综合一区二区 | 高清日韩中文字幕在线| 国产成人a∨麻豆精品| 国产又黄又爽又无遮挡在线| 亚洲成人精品中文字幕电影| 色综合站精品国产| 99热这里只有精品一区| 三级经典国产精品| 午夜a级毛片| 久久久久九九精品影院| 久久久久久久久久久丰满| 日韩一区二区视频免费看| 美女大奶头视频| 好男人在线观看高清免费视频| 色综合色国产| 亚洲自偷自拍三级| 免费看av在线观看网站| 欧美三级亚洲精品| 中国国产av一级| 午夜久久久久精精品| 亚洲成人中文字幕在线播放| 乱人视频在线观看| 精品久久久久久久久久免费视频| 在线观看免费视频日本深夜| 久久6这里有精品| 久久久a久久爽久久v久久| 亚洲美女视频黄频| 亚洲人与动物交配视频| 国产伦在线观看视频一区| 99热只有精品国产| 中文字幕免费在线视频6| 国产女主播在线喷水免费视频网站 | 一本精品99久久精品77| 亚洲精品乱码久久久久久按摩| 99riav亚洲国产免费| 亚洲人成网站在线观看播放| 国产一级毛片七仙女欲春2| 深夜a级毛片| 午夜福利在线观看吧| 午夜视频国产福利| 在线免费观看的www视频| 午夜福利成人在线免费观看| 久久精品人妻少妇| 亚洲欧美清纯卡通| 少妇高潮的动态图| 一级毛片久久久久久久久女| 美女被艹到高潮喷水动态| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 给我免费播放毛片高清在线观看| 久久久久性生活片| 亚洲欧美精品自产自拍| 性插视频无遮挡在线免费观看| 国产伦在线观看视频一区| 国产国拍精品亚洲av在线观看| 九九爱精品视频在线观看| 国产极品天堂在线| 日本欧美国产在线视频| 国产精品.久久久| 深爱激情五月婷婷| 美女被艹到高潮喷水动态| 97超视频在线观看视频| 国产亚洲av片在线观看秒播厂 | 欧美一区二区精品小视频在线| 99热全是精品| 久久精品夜夜夜夜夜久久蜜豆| 观看免费一级毛片| 日韩强制内射视频| 国产精品.久久久| 禁无遮挡网站| 蜜臀久久99精品久久宅男| 精品日产1卡2卡| www.色视频.com| 精品国内亚洲2022精品成人| 久久综合国产亚洲精品| 国产高清不卡午夜福利| 一本久久中文字幕| 久久欧美精品欧美久久欧美| 亚洲成a人片在线一区二区| 日日啪夜夜撸| 色播亚洲综合网| 午夜激情福利司机影院| 国产精品一区二区性色av| 久99久视频精品免费| 日本免费一区二区三区高清不卡| 麻豆成人午夜福利视频| 日本免费一区二区三区高清不卡| 一进一出抽搐动态| 日韩一本色道免费dvd| 美女 人体艺术 gogo| 日本免费一区二区三区高清不卡| 久久国内精品自在自线图片| 国语自产精品视频在线第100页| 免费大片18禁| 人人妻人人看人人澡| 国产成人福利小说| 国产色爽女视频免费观看| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 中文欧美无线码| 美女 人体艺术 gogo| 日韩欧美 国产精品| 久久九九热精品免费| 在线天堂最新版资源| 国产精品久久视频播放| 特级一级黄色大片| 深爱激情五月婷婷| 97超碰精品成人国产| 国产一区二区亚洲精品在线观看| 尤物成人国产欧美一区二区三区| 搡女人真爽免费视频火全软件| 国产精品蜜桃在线观看 | 淫秽高清视频在线观看| 黄色一级大片看看| 一本一本综合久久| 久久国产乱子免费精品| 亚洲在线自拍视频| 日韩成人av中文字幕在线观看| 精品不卡国产一区二区三区| 秋霞在线观看毛片| 国产成人91sexporn| 亚洲人成网站在线观看播放| 性插视频无遮挡在线免费观看| 欧美一区二区精品小视频在线| 极品教师在线视频| 国产成人精品婷婷| 国内精品一区二区在线观看| 久久久久网色| 黄片wwwwww| 99热精品在线国产| 丝袜美腿在线中文| 久久热精品热| 亚洲图色成人| 欧美性猛交黑人性爽| 黄色一级大片看看| 熟女人妻精品中文字幕| 国内揄拍国产精品人妻在线| 国产激情偷乱视频一区二区| 女人十人毛片免费观看3o分钟| 身体一侧抽搐| 人妻久久中文字幕网| 一边亲一边摸免费视频| 少妇高潮的动态图| 国产黄色视频一区二区在线观看 | 毛片一级片免费看久久久久| 国产爱豆传媒在线观看| 高清毛片免费看| 丰满乱子伦码专区| 成人亚洲欧美一区二区av| 少妇熟女欧美另类| www.色视频.com| 日韩欧美在线乱码| 亚洲va在线va天堂va国产| 免费看a级黄色片| 国产 一区 欧美 日韩| 日韩强制内射视频| 熟女电影av网| 性欧美人与动物交配| 少妇被粗大猛烈的视频| 久久精品国产亚洲av天美| 亚州av有码| 一个人看视频在线观看www免费| 在线观看午夜福利视频| 91久久精品电影网| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久大精品| 人妻制服诱惑在线中文字幕| 天天躁日日操中文字幕| 成人特级黄色片久久久久久久| 免费看av在线观看网站| 五月玫瑰六月丁香| 中文字幕av在线有码专区| 少妇裸体淫交视频免费看高清| 人人妻人人澡人人爽人人夜夜 | 深爱激情五月婷婷| 永久网站在线| 成熟少妇高潮喷水视频| 精品午夜福利在线看| 日本熟妇午夜| 熟女人妻精品中文字幕| www.色视频.com| av在线亚洲专区| 久久午夜亚洲精品久久| 五月伊人婷婷丁香| 国产激情偷乱视频一区二区| 中国美女看黄片| 最好的美女福利视频网| 午夜精品一区二区三区免费看| 欧美激情在线99| 又粗又爽又猛毛片免费看| 午夜久久久久精精品| 观看美女的网站| 国产女主播在线喷水免费视频网站 | 少妇裸体淫交视频免费看高清| 国内精品美女久久久久久| 校园春色视频在线观看| 综合色av麻豆| 精品熟女少妇av免费看| 午夜福利成人在线免费观看| 久久久久久大精品| av.在线天堂| 成熟少妇高潮喷水视频| 日韩强制内射视频| 国产午夜福利久久久久久| 特大巨黑吊av在线直播| 亚洲无线在线观看| 免费观看在线日韩| 一卡2卡三卡四卡精品乱码亚洲| 91久久精品国产一区二区三区| 亚洲中文字幕日韩| 成熟少妇高潮喷水视频| 亚洲最大成人av| 国产一级毛片七仙女欲春2| 精品久久久久久成人av| 青春草亚洲视频在线观看| 色哟哟哟哟哟哟| 亚洲图色成人| 欧美日韩乱码在线| 99久久久亚洲精品蜜臀av| 高清在线视频一区二区三区 | 国产亚洲av嫩草精品影院| 看非洲黑人一级黄片| 亚洲成av人片在线播放无| 老司机影院成人| 色综合亚洲欧美另类图片| 国产高清不卡午夜福利| 成人鲁丝片一二三区免费| 男女边吃奶边做爰视频| 国产精品久久久久久久电影| 99久国产av精品|