Effect of Baroclinicity on Vortex Axisymmetrization.Part Ⅰ：Barotropic Basic Vortex

The barotropic and baroclinic disturbances axisymmetrized by the barotropic basic vortex are examined in an idealized modeling framework consisting of two layers.Using a Wentzel-Kramers-Brillouin approach,the radial propagation of a baroclinic disturbance is shown to be slower than a barotropic disturbance,resulting in a slower linear axisymmetrization for baroclinic disturbances.The slower-propagating baroclinic waves also cause more baroclinic asymmetric kinetic energy to be transferred directly to the barotropic symmetric vortex than from barotropic disturbances,resulting in a faster axisymmetrization process in the nonlinear baroclinic wave case than in the nonlinear barotropic wave case.     

Effect of Baroclinicity on Vortex Axisymmetrization.Part Ⅱ：Baroclinic Basic Vortex

The effect of baroclinicity on vortex axisymmetrization is examined within a two-layer dynamical model.Three basic state vortices are constructed with varying degrees of baroclinicity：（i） barotropic,（ii） weak baroclinic,and （iii） strong baroclinic.The linear and nonlinear evolution of wavenumber-2 baroclinic disturbances are examined in each of the three basic state vortices.The results show that the radial propagating speed of the vortex Rossby wave at the lower level is larger with the stronger baroclinicity,resulting in a faster linear axisymmetrization process in the stronger baroclinic vortex.It is found that the nonlinear axisymmetrization process takes the longest time in the strongest baroclinic vortex among the three different basic vortices due to the weaker kinetic energy transfer from asymmetric to symmetric circulations at the lower level.A major finding in this study is that the same initial asymmetric perturbation can have different effects on symmetric vortices depending on the initial vortex baroclinicity.In numerical weather prediction models,this implies that there exists a sensitivity of the subsequent structural and intensity change solely due to the specification of the initial vertical shear of the tropical cyclone vortex.     

Effect of Initial Vortex Intensity Correction on Tropical Cyclone Intensity Prediction:A Study Based on GRAPES_TYM

Predicting the intensity of tropical cyclones(TCs)is challenging in operational weather prediction systems,partly due to the difficulty in defining the initial vortex.In an attempt to solve this problem,this study investigated the effect of initial vortex intensity correction on the prediction of the intensity of TCs by the operational numerical prediction system GRAPES_TYM(Global and Regional Assimilation and Prediction System_Typhoon Model)of the National Meteorological Center of the China Meteorological Administration.The statistical results based on experiments using data for major TCs in 2018 show that initial vortex intensity correction can reduce the errors in mean intensity for up to 120-h integration,with a noticeable decrease in the negative bias of intensity and a slight increase in the mean track error.The correction leads to an increase in the correlation coefficient of Vmax(maximum wind speed at 10-m height)for the severe typhoon and super typhoon stages.Analyses of the errors in intensity at different stages of intensity(including tropical storms,severe tropical storms,typhoons,severe typhoons,and super typhoons)show that vortex intensity correction has a remarkable positive influence on the prediction of super typhoons from 0 to 120h.Analyses of the errors in intensity for TCs with different initial intensities indicate that initial vortex correction can significantly improve the prediction of intensity from 24 to 96 h for weak TCs(including tropical storms and severe tropical storms at the initial time)and up to 24 h for strong TCs(including severe typhoons and super typhoons at the initial time).The effect of the initial vortex intensity correction is more important for developing TCs than for weakening TCs.     

Study on Instability in Baroclinic Vortex Symmetric Disturbance under Effect of Nonuniform Environmental Parameters

With the aid of the baroclinicity parameter M2,inertial instability parameter F2 and the stratification instability paramter N2 as the slowly varying function both spatially and temporally,an energetic equation is derived of symmetric perturbation waves in baroclinic vortices in the framework of progressively changing wavetrain theory,or WKB,alongside the examination of effects of these parameters upon the vortex disturbance.     

The Solitary Wave of Barotropic Atmosphere on a Sphere

In this paper, we have investigated large scale disturbance in barotropic atmosphere on a sphere. It demonstrates that: considering nonlinear effects of interaction between waves in barotropic vorticity equation, the wave packet of the disturbance is governed by the famous equation-nonlinear Schrodinger equation. For the solitary wave, two factors are very important: one is spherical effect of the disturbance and the other is meridional shear of blocking high. In comparison with the results of local Cartesian coordinates, the former factor is the individuality of spherical soliton.     

Impact of Cloud Microphysical Processes on the Simulation of Typhoon Rananim near Shore.Part Ⅰ: Cloud Structure and Precipitation Features

By using the Advanced Regional Eta-coordinate Model (AREM),the basic structure and cloud features of Typhoon Rananim are simulated and verified against observations.Five sets of experiments are designed to investigate the effects of the cloud microphysical processes on the model cloud structure and precipitation features.The importance of the ice-phase microphysics,the cooling effect related to microphysical characteristics change,and the influence of terminal velocity of graupel are examined.The results indicate that the cloud microphysical processes impact more on the cloud development and precipitation features of the typhoon than on its intensity and track.Big differences in the distribution pattern and content of hydrometeors,and types and amount of rainfall occur in the five experiments,resulting in different heating and cooling effects.The largest difference of 24-h rain rate reaches 52.5 mm h-1.The results are summarized as follows:1) when the cooling effect due to the evaporation of rain water is excluded,updrafts in the typhoon's inner core are the strongest with the maximum vertical velocity of-19 Pa s-1 and rain water and graupel grow most dominantly with their mixing ratios increased by 1.8 and 2.5 g kg-1,respectively,compared with the control experiment; 2) the melting of snow and graupel affects the growth of rain water mainly in the spiral rainbands,but much less significantly in the eyewall area; 3) the warm cloud microphysical process produces the smallest rainfall area and the largest percentage of convective precipitation (63.19％),while the largest rainfall area and the smallest percentage of convective precipitation (48.85％) are generated when the terminal velocity of graupel is weakened by half.     

An Updated Coupled Model for Land-Atmosphere Interaction. Part Ⅰ: Simulations of Physical Processes

A new two-way land-atmosphere interaction model (R42_AVIM) is fulfilled by coupling the spectral at- mospheric model (SAMIL_R42L9) developed at the State Key Laboratory of Numerical Modeling for Atmo- spheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sci- ences (LASG/IAP/CAS) with the land surface model, Atmosphere-Vegetation-Interaction-Model (AVIM). In this coupled model, physical and biological components of AVIM are both included. Climate base state and land surface physical fluxes simulated by R42_AVIM are analyzed and compared with the results of R42_SSIB [which is coupled by SAMIL_R42L9 and Simplified Simple Biosphere (SSIB) models]. The results show the performance of the new model is closer to the observations. It can basically guarantee that the land surface energy budget is balanced, and can simulate June-July-August (JJA) and December-January- February (DJF) land surface air temperature, sensible heat flux, latent heat flux, precipitation, sea level pressure and other variables reasonably well. Compared with R42_SSIB, there are obvious improvements in the JJA simulations of surface air temperature and surface fluxes. Thus, this land-atmosphere coupled model will offer a good experiment platform for land-atmosphere interaction research.     

Dynamical and Thermal Problems in Vortex Development and Movement. Part I: A PV–Q View

Based on the Lagrangian change equation of vertical vorticity deduced from the equation of threedimensional Ertel potential vorticity(PV e),the development and movement of vortex are investigated from the view of potential vorticity and diabatic heating(PV-Q).It is demonstrated that the asymmetric distribution in the vortex of the non-uniform diabatic heating in both vertical and horizontal can lead to the vortex’s development and movement.The theoretical results are used to analyze the development and movement of a Tibetan Plateau(TP) vortex(TPV),which appeared over the TP,then slid down and moved eastward in late July 2008,resulting in heavy rainfall in Sichuan Province and along the middle and lower reaches of the Yangtze River.The relative contributions to the vertical vorticity development of the TPV are decomposed into three parts:the diabatic heating,the change in horizontal component of PV e(defined as PV 2),and the change in static stability θ z.The results show that in most cases,diabatic heating plays a leading role,followed by the change in PV 2,while the change of θ z usually has a negative impact in a stable atmosphere when the atmosphere becomes more stable,and has a positive contribution when the atmosphere approaches neutral stratification.The intensification of the TPV from 0600 to 1200 UTC 22 July 2008 is mainly due to the diabatic heating associated with the precipitation on the eastern side of the TPV when it uplifted on the up-slope of the northeastern edge of the Sichuan basin.The vertical gradient of diabatic heating makes positive(negative) PV e generation below(above) the maximum of diabatic heating;the positive PV e generation not only intensifies the low-level vortex but also enhances the vertical extent of the vortex as it uplifts.The change in PV e due to the horizontal gradient of diabatic heating depends on the vertical shear of horizontal wind that passes through the center of diabatic heating.The horizontal gradient of diabatic heating makes positive(negative) PV e generation on the right(left) side of the vertical shear of horizontal wind.The positive PV e generation on the right side of the vertical shear of horizontal wind not only intensifies the local vertical vorticity but also affects direction of movement of the TPV.These diagnostic results are in good agreement with the theoretic results developed from the PV-Q view.     

Development and Validation of an Evaporation Duct Model.Part Ⅰ:Model Establishment and Sensitivity Experiments

Based on the Coupled Ocean-Atmospheric Response Experiment（COARE）bulk algorithm and the Naval Postgraduate School（NPS）model,a universal evaporation duct（UED）model that can flexibly accommodate the latest improvements in component（such as stability function,velocity roughness,and scalar roughness）schemes for different stratification and wind conditions,is proposed in this paper.With the UED model,the sensitivity of the model-derived evaporation duct height（EDH）to stability function（Ψ）,ocean wave effect under moderate to high wind speeds,and scalar roughness length parameterization,is investigated,and relative contributions of these factors are compared.The results show that the stability function is a key factor influencing the simulated EDH values.Under unstable conditions,the EDH values from stability functions of Fairall et al.（1996）and Hu and Zhang（1992）are generally higher than those from others;while under stable conditions,unreasonably high EDHs can be avoided by use of the stability functions of Hu and Zhang（1992）and Grachev et al.（2007）.Under moderate to high wind speeds,the increase in velocity roughness length z₀ due to consideration of the true ocean wave effect acts to reduce modeled EDH values;this trend is more pronounced under stable conditions.Although the scalar roughness length parameterization has a minor effect on the model-derived EDH,a positive correlation is found between the scalar roughness length z_0qand the model-derived EDH.     

Monthly Extended Predicting Experiments with Nonlinear Regional Prediction.PartⅠ:Prediction of Zonal Mean Flow

Systematic errors have recently been founded to be distinct in the zonal mean component forecasts, which account for a large portion of the total monthly-mean forecast errors. To overcome the difficulty of numerical model, the monthly pentad-riean nonlinear dynamic regional prediction models of the zonal mean geopotential height at 200, 300, 500, and 700 hPa based on a large number of historical data (NCEP/NCAR reanalysis data) were constituted by employing the local approximation of the phase space reconstruction theory and nonlinear spatio-temporal series prediction method. The 12-month forecast experiments of 1996 indicated that the results of the nonlinear model are better than those of the persistent, climatic prediction, and T42L9 model either over the high- and mid-latitude areas of the Northern and Southern Hemispheres or the tropical area. The root-mean-square of the monthly-mean height of T42L9 model was considerably decreased with a change of 30.4%, 26.6%, 82.6%, and 39.4%, respectively, over the high- and mid-latitudes of the Northern Hemisphere, over the high- and mid-latitudes of the Southern Hemisphere, over the tropics and over the globe, and also the corresponding anomaly correlation coefficients over the four areas were respectively increased by 0.306-0.312, 0.304-0.429, 0.739-0.746, and 0.360-0.400 (averagely a relative change of 11.0% over the globe) by nonlinear correction after integration, implying that the forecasts given by nonlinear model include more useful information than those of T42L9 model.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅰ: Retrieval of Microphysical Elements and Vertical Velocity

Assuming that cloud reaches static state in the warm microphysical processes, water vapor mixing ratio(qv), cloud water mixing ratio (qc), and vertical velocity (w) can be calculated from rain water mixing ratio (qr)- Through relation of Z-qr, qr can be retrieved by radar reflectivity factor (Z). Retrieval results indicate that the distributions of mixing ratios of vapor, cloud, rain, and vertical velocity are consistent with radar images, and the three-dimensional spatial structure of the convective cloud is presented. Treating q,v saturated at the echo area, the retrieved qr is about 0.1 g kg-1, qc is always less than 0.3 g kg-1, w is usually below 0.5 m s-1, and rain droplet terminal velocity (vr) is around 5.0 m s-1 in the place where radar reflectivity factor is about 25 dBz; in the place where echo is 45 dBz, the retrieved qr and qc are always about 3.0 g kg-1, w is greater than 5.0 m s-1, and vr is around 7.0 m s-1. In the vertical, the maximum updraft velocity is greater than 3.0 m s-1 at the height of around 5.0 kin, the maximum cloud water content is about 3.0 g kg-1 above 5 km and the maximum rain water content is about 3.0 g kg-1 below 6 kin. Due to the assumption that the cloud is in static state, there will be some errors in the retrieved variables within the clouds which axe rapidly growing or dying-out, and in such cases, more sophisticated radar data control technique will help to improve the retrieval results.     

Meridional Wind Stress Anomalies over Tropical Pacific and the Onset of El Nino. Part Ⅰ: Data Analysis

对热带太平洋海表经向风应力异常与EI Nino事件之间的关系进行了诊断分析。结果表明，超前的经向风应力距平场与NINO3区(150°-90°W，5°S-5°N)的海面温度异常(SSTA)有显著的超前相关，这种相关性在超前6个月甚至更早一些就有显示。利用奇异值分解方法分析超前的经向风应力距平场与太平洋海表温度异常场之间的耦合模，结果表明对应于赤道中东太平洋的海面温度异常升高，大气风应力场在超前6个月甚至更早的时候，在赤道中东太平洋表现为辐合的经向异常风应力场，即赤道以北为北风异常应力，赤道以南为南风异常应力。这种耦合模的时间系数与NINO3 SSTA指数所表示的El Nino事件有很好的对应关系，表明这种耦合模反映的正是超前的经向风应力异常与El Nino事件所对应的海表温度异常之间的相关模态。通过与热带西太平洋纬向风应力异常的比较，赤道中东太平洋辐合的经向风应力异常与El Nino事件发生的同样具有重要的联系。     

Impact of Cloud Microphysical Processes on the Simulation of Typhoon Rananim near Shore.Part Ⅱ:Typhoon Intensity and Track

The impact of cloud microphysical processes on the simulated intensity and track of Typhoon Rananim is discussed and analyzed in the second part of this study.The results indicate that when the cooling effect due to evaporation of rain water is excluded,the simulated 36-h maximum surface wind speed of Typhoon Rananim is about 7 m s-1 greater than that from all other experiments; however,the typhoon landfall location has the biggest bias of about 150 km against the control experiment.The simulated strong outer rainbands and the vertical shear of the environmental flow are unfavorable for the deepening and maintenance of the typhoon and result in its intensity loss near the landfall.It is the cloud microphysical processes that strengthen and create the outer spiral rainbands,which then increase the local convergence away from the typhoon center and prevent more moisture and energy transport to the inner core of the typhoon.The developed outer rainbands are supposed to bring dry and cold air mass from the middle troposphere to the planetary boundary layer (PBL).The other branch of the cold airflow comes from the evaporation of rain water itself in the PBL while the droplets are falling.Thus,the cut-off of the warm and moist air to the inner core and the invasion of cold and dry air to the eyewall region are expected to bring about the intensity reduction of the modeled typhoon.Therefore,the deepening and maintenance of Typhoon Rananim during its landing are better simulated through the reduction of these two kinds of model errors.     

Numerical Simulation Study on the Impacts of Tropospheric and CO_2 Concentration Changes on Winter Wheat. Part I: Model Description

Ozone is well documented as the air pollutant most damaging to agricultural crops and other plants. It is reported that tropospheric O3 concentration increases rapidly in recent 20 years. Evaluating and predicting impacts of ozone concentration changes on crops are drawing great attention in the scientific community. In China, main study method about this filed is controlled experiments, for example, Open Top Chambers. But numerical simulation study about impacts of ozone on crops with crop model was developed slowly, what is more, the study about combined impacts of ozone and carbon dioxide has not been reported. The improved agroecosystem model is presented to evaluate simultaneously impacts of tropospheric O3 and CO2 concentration changes on crops in the paper by integrating algorithms about impacts of ozone on photosynthesis with an existing agroecosystem biogeochemical model (named as DNDC). The main physiological processes of crop growth (phenology, leaf area index, photosynthesis, respiration, assimilated allocation and so on) in the former DNDC are kept. The algorithms about impacts of ozone on photosynthesis and winter wheat leaf are added in the modified DNDC model in order to reveal impacts of ozone and carbon dioxide on growth, development, and yield formation of winter wheat by coupling the simulation about impacts of carbon dioxide on photosynthesis of winter wheat which exists in the former DNDC. In the paper, firstly assimilate allocation algorithms and some genetic parameters (such as daily thermal time of every development stage) were modified in order that DNDC can be applicable in North China. Secondly impacts of ozone on crops were simulated with two different methods- one was impacts of ozone on light use efficiency, and the other was direct effects of ozone on leaves photosynthesis. The latter simulated results are closer to experiment measurements through comparing their simulating results. At last the method of direct impacts of ozone on leaf growth is adopted and the coefficients about impacts of ozone on leaf growth and death are ascertained. Effects of climate changes, increasing ozone, and carbon dioxide concentration on agroecosystem are tried to be simulated numerically in the study which is considered to be advanced and credible.