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 qv 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 km, 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 km. Due to the assumption that the cloud is in static state, there will be some errors in the retrieved variables within the clouds which are rapidly growing or dying-out, and in such cases, more sophisticated radar data control technique will help to improve the retrieval results.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅱ: Numerical Simulation Experiments

Microphysics elements and vertical velocity retrieved were incorporated using the nudging method into the initial data assimilation of GRAPES (Global/Regional Assimilation and Prediction System) model.Simulation experiments indicated that nudging technique was effective in forcing the model forecast gradually consistent to the observations, yielding the thermodynamically and dynamically balanced analysis field. As viewed from the simulation results, water vapor is vital to precipitation, and it is a governing factor for the amount and duration of precipitation. The initial cloud water, rain water, and vertical velocity determine the strength distribution of convection and precipitation at the beginning time of forecast; the horizontal wind field steers the motion of the mesoscale weather system embedded in and impacts the position of precipitation zone to a large extent. The simulation experiments show that the influence of the initial retrieval data on prediction weakens with the increase of forecast time, and within the first hour of forecast, the retrieval data have an important impact on the evolution of the weather system, but its influence becomes trivial after the first three hours. Changing the nudging coefficient and the integral time-spacing of numerical model will bring some influences to the results. Herein only one radar reflectivity was used, the radar observations did not cover the whole model domain, and some empirical parameters were used in the retrieval method, therefore some differences still lie between simulation and observation to a certain extent, and further studies on several aspects are expected.     

Application of Radar Reflectivity Factor in Initializing Cloud-Resolving Mesoscale Model. Part Ⅱ： Numerical Simulation Experiments

Microphysics elements and vertical velocity retrieved were incorporated using the nudging method into the initial data assimilation of GRAPES （Global/Regional Assimilation and Prediction System） model. Simulation experiments indicated that nudging technique was effective in forcing the model forecast gradually consistent to the observations, yielding the thermodynamically and dynamically balanced analysis field. As viewed from the simulation results, water vapor is vital to precipitation, and it is a governing factor for the amount and duration of precipitation. The initial cloud water, rain water, and vertical velocity determine the strength distribution of convection and precipitation at the beginning time of forecast; the horizontal wind field steers the motion of the mesoscale weather system embedded in and impacts the position of precipitation zone to a large extent. The simulation experiments show that the influence of the initial retrieval data on prediction weakens with the increase of forecast time, and within the first hour of forecast, the retrieval data have an important impact on the evolution of the weather system, but its influence becomes trivial after the first three hours. Changing the nudging coefficient and the integral time-spacing of numerical model will bring some influences to the results. Herein only one radar reflectivity was used, the radar observations did not cover the whole model domain, and some empirical parameters were used in the retrieval method, therefore some differences still lie between simulation and observation to a certain extent, and further studies on several aspects are expected.     

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.     

Validation of the Atmospheric Infrared Sounder Retrieval Products over China and Their Application in Numerical Model

The atmospheric infrared sounder (AIRS) instrument onboard Aqua Satellite is a high spectral resolution infrared sounder. In recent years, AIRS has gradually become the primary method of atmospheric vertical observations. To examine the validation of AIRS retrieval products (V3.0) over China, the AIRS surface air temperature retrievals were compared with the ground observations obtained from 540 meteorological stations in July 2004 and January 2005, respectively. The sources of errors were considerably discussed. Based on the error analysis, the AIRS retrieved surface air temperature products were systemi-cally corrected. Moreover, the AIRS temperature and humidity profile retrievals were compared with T213 numerical forecasting products. Because T213 forecasting products are not the actual atmospheric states, to further verify the validation, the AIRS temperature and humidity profile products were assimilated into the MM5 model through the analysis nudging. In this paper, the case on February 14, 2005 in North China was simulated in detail. Then, we investigated the effects of AIRS retrievals on snowfall, humidity field, vertical velocity field, divergence field, and cloud microphysical processes. The major results are: (1) the errors of AIRS retrieved surface air temperature products are largely systematic deviations, for which the influences of terrain altitude and surface types are the major reasons; (2) the differences between the AIRS atmospheric profile retrievals and T213 numerical prediction products in temperature are generally less than 2 K, the differences in relative humidity are generally less than 25%; and (3) the AIRS temperature and humidity retrieval products can adjust the model initial field, and thus can improve the capacity of snowfall simulation to some extent.     

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.     

The Vertical Transport of Air Pollutants by Convective Clouds. Part I: A Non-Reactive Cloud Transport Model

A convective cloud transport model, without chemical processes, is developed by joining a set of concentration conservative equations into a two-dimensional, slab-symmetric and fully elastic numerical cloud model, and a numerical experiment is completed to simulate the vertical transport of ground-borne, inert gaseous pollutant by deepthunderstorm. The simulation shows that deep convective storm can very effectively transport high concentrated pollutant gas from PBL upward to the upper troposphere in 30 to 40 minutes, where the pollutant spreads laterally outward with strong anvil outflow, forming an extensive high concentration area. Meanwhile, relatively low concentration areas are formed in PBL both below and beside the cloud, mainly caused by dynamic pumping effect and sub-cloud downdraft flow. About 80% of the pollutant gas transported to the upper troposphere is from the layer below 1.5 km AGL (above ground level).     

Coamps<Superscript>®</Superscript>-Les: Model Evaluation and Analysis of Second-and Third-Moment Vertical Velocity Budgets

The Naval Research Laboratory Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS^®) has been extended to perform as a large-eddy simulation (LES) model. It has been validated with a series of boundary-layer experiments spanning a range of cloud nighttime, and includes a nighttime stratocumulus case, a trade wind cumulus layer, shallow cumulus convection over land, and a mixed regime consisting of cumulus clouds under broken stratocumulus. COAMPS-LES results are in good agreement with other models for all the cases simulated. Exact numerical budgets for the vertical velocity second\((\overline{w^{'2}})\) and third moment\((\overline{w^{'3}})\) have been derived for the stratocumulus and trade wind cumulus cases. For the\(\overline{w^{'3}}\) budget in the stratocumulus, the buoyancy contribution from the updraughts and downdraughts largely cancel each other due to their similar magnitudes but opposite signs. In contrast, for the cumulus layer, the negative buoyancy contribution from the environmental downdraughts is negligible and the positive contribution from the updraughts completely dominates due to the conditional instability in the environment. As a result,\(\overline{w^{'3}}\) is significantly larger in the cumulus than in the stratocumulus layer.     

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.     

Applications of AMSR-E Measurements for Tropical Cyclone Predictions Part Ⅰ: Retrieval of Sea Surface Temperature and Wind Speed

Existing satellite microwave algorithms for retrieving Sea Surface Temperature(Sst)and wind(SSW)are applicable primarily for non-raining cloudy conditions.With the launch of the Earth Observing System (EOS)Aqua satellite in 2002,the Advanced Microwave Scanning Radiometer(AMSR-E)onboard provides some unique measurements at lower frequencies which are sensitive to ocean surface parameters under ad-verse weather conditions.In this study,a new algorithm is developed to derive SST and SSW for hurricane predictions such as hurricane vortex analysis from the AMSR-E measurements at 6.925 and 10.65 GHz.In the algorithm,the effects of precipitation emission and scattering on the measurements are properly taken into account.The algorithm performances are evaluated with buoy measurements and aircraft dropsonde data.It is found that the root mean square (RMS) errors for SST and SSW are about 1.8K and 1.9m s(-1),respectively,when the results are compared with the buoy data over open oceans under precipitating clouds (e.g.,its liquid water path is larger than 0.5 mm),while they are 1.1 K for SST and 2.0 ms(-1)for SSW,respectively,when the retrievals are validated against the dropsonde measurements over warm oceans.These results indicate that our newly developed algorithm catl provide some critical surface information for trop-ical cycle predictions.Currently,this newly developed algorithm has been implemented into the hybrid variational scheme for the hurricane vortex analysis to provide predictions of SST and SSW fields.     

Retrieval of Tropospheric CO Profiles Using Correlation Radiometer.Ⅱ: Effects of Other Gases and the Retrieval in Cloudy Atmosphere

l.Intr0ducti0nAschemeforretrievingCOprofilesinthetroPOspherehasbeendiscussedinthepreviousPartI(WuandGille,partI)foragascorrelationradiometerMeasurementsofPollutionintheTroPOsphere(M0PITT)workingatthe4.6pmwaveband.Thebasicequationshavebeendiscussed.TheuseofthewidebandsignalforprovidingthesurfacetemPeratureandanaddi-tionalchannelforthecolumnarC0amounttoimprovetheretrievalinthenearsurfacelayerhasbeentested.Itisfoundthaterrorsinthetemperatureprofilemayincreaseerrorsinthere-trievedprofiles…     

Application of Linear Thermodynamics to the Atmospheric System. Part Ⅰ: Linear Phenomenological Relations and Thermodynamic Property of the Atmospheric System

从一般的热力学原理或其它自然原理对唯象关系所强加的限制，能够演绎出大气系统的一系列热力学性质。利用非平衡态线性热力学导出了湍流K闭合理论中湍流交换系数同唯象系数的关系，从理论上证明大气系统热量湍流输送同水泡之间存在交叉耦合，还导出了湍流强度同速度和位温梯度的关系，从而证明速度和位温空间分布的非均匀性是湍流之源。并证明湍流强度定理，不可压缩气体和各向同性湍流大气中，湍流强度正比于速度与位温梯度的标积。进而证明大气涡旋定理，位温的切变将导致涡旋运动或各种环流运动，速度涡度等于速度同位温相对梯度的矢积。展现了线性热力学在大气系统的应用前景。     

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.     

Mesoscale and Microphysical Characteristics of a Double Rain Belt Event in SouthChina on May 10–13, 2022

A second rain belt sometimes occurs ahead of a frontal rain belt in the warm sector over coastal South China,leading to heavy precipitation. We examined the differences in the mesoscale characteristics and microphysics of thefrontal and warm sector rain belts that occurred in South China on May 10–13, 2022. The southern rain belt occurred in anenvironment with favorable mesoscale conditions but weak large-scale forcing. In contrast, the northern rain belt wasrelated to low-level horizontal shear and the surface-level front. The interaction between the enhanced southeasterly windsand the rainfall-induced cold pool promoted the persistent growth of convection along the southern rain belt. The con vective cell propagated east over the coastal area, where there was a large temperature gradient. The bow-shaped echo inthis region may be closely related to the rear-inflow jet. By contrast, the initial convection of the northern rain belt wastriggered along the front and the region of low-level horizontal shear, with mesoscale interactions between the enhancedwarm-moist southeasterly airflow and the cold dome associated with the earlier rain. The terrain blocked the movement ofthe cold pool, resulting in the stagnation of the frontal convective cell at an early stage. Subsequently, a meso-γ-scalevortex formed during the rapid movement of the convective cell, corresponding to an enhancement of precipitation. Therepresentative raindrop spectra for the southern rain belt were characterized by a greater number and higher density ofraindrops than the northern rain belt, even though both resulted in comparable hourly rainfalls. These results help us betterunderstand the characteristics of double rain belts over South China.