A Density Correction for Lagrangian Particle Dispersion Models

Current Lagrangian particle dispersion models, used to simulate the dispersion of passive tracers in the turbulent planetary boundary layer (PBL), assume that the density is constant within the PBL. In deep PBLs, where the density at the boundary-layer top may be lower by more than 20% than at the surface, this assumption leads to errors in the tracer concentrations on the order of 10%. In the presence of a vertical wind shear, this also leads to inaccurate calculations of the horizontal tracer transport. To remove this deficiency, a Langevin equation is presented that contains a density correction term. The effect of the density correction is studied using data from a large-scale tracer experiment. It is found that for this experiment, the main effect of the density correction is an increase in the surface tracer concentrations, whereas the horizontal tracer transport patterns remain largely unaffected.     

An eddy-resolving simulation of the diurnal variation of fair-weather convection and tracer transport

Tracer transport in the atmosphere is controlled not only by synoptic-scale to mesoscale weather disturbances but also by microscale boundary-layer processes especially under fair-weather conditions. The present study investigates numerically the diurnal variation of boundary-layer convection and cumulus clouds and their role in transporting tracers by conducting high-resolution simulations that explicitly resolve turbulent eddies. The transport of dust aerosols in a desert area under two distinct stability conditions is specifically examined. Convection plays a significant role in transporting dust upward; in other words, the vertical depth of the dust transport is critically determined by the depth of convection. Deep convection is effective in transporting dust into the free atmosphere. The early morning stratification strongly regulates the temporal evolution and the vertical growth of convection and therefore the amount of tracer emission and transport. A sensitivity to model resolution of O (1 km) in a cloud-resolving simulation range is also examined. A proper parameterization for activating microscale convection is required for representing the diurnal variation of convection and tracer transport.     

Idealised Simulations of Daytime Pollution Transport in a Steep Valley and its Sensitivity to Thermal Stratification and Surface Albedo

Numerical simulations of tracer transport in an idealised, east-west aligned valley are performed with the Regional Atmospheric Modeling System (RAMS), both two-dimensional and three-dimensional. The results are qualitatively consistent with wintertime observations in the Austrian Inn Valley. The simulations show an asymmetry in wind circulation and tracer distribution between the valley sidewalls according to the orientation of the slope with respect to the sun. Two-dimensional sensitivity experiments are run to investigate the influence of vertical inhomogeneities in thermal stratification and vegetation coverage on the slope-wind circulation and therewith the tracer transport. It is shown that a transition to a layer of higher stability or to a region with higher surface albedo causes a reduction of the mass flux in the upslope-wind layer and due to mass continuity a quasi-horizontal transport out of the slope-wind layer.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part I: A Parametrization-Oriented Conditional Sampling in Large-Eddy Simulations

A conditional sampling based on the combination of a passive tracer emitted at the surface and thermodynamic variables is proposed to characterise organized structures in large-eddy simulations of cloud-free and cloudy boundary layers. The sampling is evaluated against more traditional sampling of dry thermals or clouds. It enables the characterization of convective updrafts from the surface to the top of the boundary layer (or the top of cumulus clouds), describing in particular the transition from the sub-cloud to the cloud layer, and retrieves plume characteristics, entrainment and detrainment rates, variances and fluxes. This sampling is used to analyze the contribution of boundary-layer thermals to vertical fluxes and variances.     

热带深对流云对CO、NO、NOx和O3的垂直输送作用

利用2005年11月至2006年2月ACTIVE (Aerosol and Chemical Transport in tropIcal conVEction) 外场试验期间在澳大利亚北部达尔文地区取得的CO、O3、NO和NOx飞机探测资料, 并结合HYSPLIT后向轨迹模式结果, 分析这几种气体成分在对流卷云砧内外的分布情况, 并探讨热带深对流云对于污染气体的垂直输送作用。分析结果显示, 在孤立对流云卷云砧中, 云砧内部O3、NO、NOx浓度均大于云外; 而CO则不同, 只有在近地面浓度高时才如此, 在近地面浓度较小时, 卷云砧内部的浓度反而小于云外。进一步分析造成这两类气体分布差异的原因, 发现CO主要借助深对流云将对流层下层以及对流云周围环境中的CO夹卷并动力垂直输送到对流云顶部卷云砧中, 而对于O3、NO和NOx来说, 除了上述作用以外, 还可能与对流云内部其他物理机制(如闪电), 造成新的O3、NO和NOx有关, 这些新生气体随着风暴内部强烈的上升气流被最终输送进云砧中。     

半透明云风矢量高度算法中代表运动像元的使用

半透明云风矢量高度指定是卫星风矢量算法的重要部分,需要使用来自半透明云体的辐射和云下背景辐射两个变量。为了更精确地获得来自半透明云体的辐射,采用像元对追踪相关系数的贡献和红外亮温两个变量进行统计分析,将追踪图像块中的像元分为3个部分:小贡献像元、暖区段大贡献像元和冷区段大贡献像元;对暖区段、冷区段内大贡献像元的特征进行考察表明,冷区段的大贡献像元更能代表追踪图像块运动,称为代表运动像元。选择代表运动像元参与计算来自半透明云体的辐射,对半透明云风矢量高度指定算法进行改进,在该算法中分别使用和不使用代表运动像元,计算FY-2气象卫星风矢量,将结果与欧洲中期天气预报中心（ECMWF）分析场进行对比表明,在半透明云高度指定算法中使用代表运动像元,FY-2气象卫星风矢量误差明显降低。     

Large-Eddy Simulation of air Pollution Dispersion in the Nocturnal Cloud-Topped Atmospheric Boundary Layer

Effects of stratocumulus clouds on the dispersion of contaminants are studied in the nocturnal atmospheric boundary layer. The study is based on a large-eddy simulation (LES) model with a bulk parametrization of clouds. Computations include Lagrangian calculations of atmospheric dispersion of a passive tracer released from point sources at various heights above the ground. The results obtained show that the vertical diffusion is non-Gaussian and depends on the location of a source in the boundary layer.     

大涡模式垂直分辨率在夹卷及示踪物垂直传输模拟结果分析中的作用

基于敦煌野外观测资料和大涡模式,研究了垂直方向不同尺度湍涡对夹卷及示踪物垂直传输的影响,明确了模式垂直分辨率在模拟结果分析中的作用。结果表明:垂直方向上小尺度湍涡对夹卷作用贡献更大,小尺度湍涡较多时夹卷层相对更暖,而夹卷层厚度、夹卷强度和风速变化受垂直方向湍涡尺度影响较小。当垂直分辨率为50 m时,越往夹卷层上部,上升气流和下沉气流分布较多且强度较大;分辨率为10、20和30 m时,夹卷层各高度垂直速度、位温和示踪物浓度分布较接近。另外,垂直方向湍涡尺度对示踪物垂直传输高度影响不大,而对示踪物的空间分布有一定作用。当大尺度湍涡较多且强度较强时,越有利于将高浓度的示踪物向上传输。综合考虑到模式采用较高分辨率模拟时产生的噪音及计算时间等问题,认为模式采用30 m的垂直分辨率,既能较好地模拟出夹卷层平均结构特征,又能模拟出夹卷层湍流的精细分布,是较为理想的选择。     

A numerical study of the marine stratocumulus cloud layer

A one-dimensional grid-level model including longwave radiative transfer and a level-4 second-order turbulent transfer closure which contains prognostic equations for turbulent quantities, is used to study the physics and dynamics of inversion-capped marine stratocumulus clouds.A set of numerical experiments had been performed to examined the role of sea surface temperature, large-scale vertical velocity, wind speed, and vertical wind shear in the formation and the structure of low-level clouds. For a given sea surface and geostrophic wind speed, stratocumulus clouds can grow higher with smaller large-scale subsidence as less dry air entrains into the cloud. Clouds grow higher with higher sea surface temperature for a given geostrophic wind speed and large-scale subsidence as a result of enhanced moist convection. In high wind speeds, the entire cloud deck is lifted up because of larger surface energy flux. In the budget studies of the turbulent kinetic energy (TKE), the buoyancy term is a major source term when the wind speed and the vertical shear are small across the inversion top. When the wind speed and the vertical wind shear across the inversion top become large, the mixed layer is decoupled into a cloud and a subcloud layer. In the TKE budget studies, the shear generation term becomes an important term in the budgets of the TKE and the variance of vertical velocity.     

对流降水云中大气垂直运动反演及个例试验

为深入认识对流降水云结构及动力特征,基于降水频段调频连续波5520 MHz垂直指向雷达(VPR-CFMCW),使用地面至15 km高度的反射率因子及径向速度,建立对流降水云中大气垂直运动的反演方法,分析对流垂直结构及大气垂直运动随高度分布的演变特征.对在广东龙门测站探测的2019年4月20-22日前汛期4次对流降水进行...     

Vertical Structures of Convective and Stratiform Clouds in Boreal Summer over the Tibetan Plateau and Its Neighboring Regions

Cloud is essential in the atmosphere, condensing water vapor and generating strong convective or large-scale persistent precipitation. In this work, the relationships between cloud vertical macro- or microphysical properties, radiative heating rate, and precipitation for convective and stratiform clouds in boreal summer over the Tibetan Plateau (TP) are analyzed and compared with its neighboring land and tropical oceans based on CloudSat/CALIPSO satellite measurements and TRMM precipitation data. The precipitation intensity caused by convective clouds is twofold stronger than that by stratiform clouds. The vertical macrophysics of both cloud types show similar features over the TP, with the region weakening the precipitation intensity and compressing the cloud vertical expansion and variation in cloud top height, but having an uplift effect on the average cloud top height. The vertical microphysics of both cloud types under conditions of no rain over the TP are characterized by lower-level ice water, ice particles with a relatively larger range of sizes, and a relatively lower occurrence of denser ice particles. The features are similar to other regions when precipitation enhances, but convective clouds gather denser and larger ice particles than stratiform clouds over the TP. The atmospheric shortwave (longwave) heating (cooling) rate strengthens with increased precipitation for both cloud types. The longwave cooling layer is thicker when the rainfall rate is less than 100 mm d^?1, but the net heating layer is typically compressed for the profiles of both cloud types over the TP. This study provides insights into the associations between clouds and precipitation, and an observational basis for improving the simulation of convective and stratiform clouds over the TP in climate models.     

半透明云风矢量高度算法中云下背景辐射的估计

半透明云风矢量高度指定是卫星风矢量算法的重要部分,需要来自半透明云体的辐射和云下背景辐射两个变量。云下背景辐射发生在云层下面,未被卫星直接观测到,为了在半透明云风矢量高度指定算法中更精确地获得云下背景辐射,使用风矢量附近无云区的红外/水汽散点图,估计云下背景辐射。分析表明:在追踪区域里存在无云区的情况下,追踪区的最高红外亮温可代表红外通道的背景辐射;而水汽通道的背景辐射,却在红外亮温高值区段内水汽亮温相对较低区段。追踪区内找不到无云区时应扩大搜索范围,找到无云区后可估计云下背景辐射。在半透明云风矢量高度指定算法中使用云下背景辐射估计的改进算法前后,计算FY-2气象卫星进行风矢量,并将结果与欧洲中期天气预报中心（ECMWF）分析场进行对比表明,在半透明风矢量高度指定算法中使用云下背景辐射估计,FY-2气象卫星风矢量误差明显降低。     

风切变对边界层对流影响的大涡模拟研究

利用"西北干旱区陆-气相互作用野外观测实验"加密观测期间在敦煌站的观测资料以及大涡模式,模拟了对流边界层的发展,以及示踪物从混合层向残留层传输的时空变化。模拟的对流边界层的结构及演变特征与实测结果基本一致。进一步通过有风切变和无风切变的敏感性数值试验,研究了风切变对垂直速度、位温和示踪物浓度的水平分布以及示踪物传输高度的影响。研究结果表明,在有风切变的试验中(甚至风切变仅存在于近地层中),对流边界层的增长加强,而且示踪物被传输的高度也较高。与浮力驱动的对流边界层相比,由浮力和风切变共同驱动的边界层中上升气流较弱而下沉气流较强,但前者的上升气流与下沉气流的分布在垂直方向上更为倾斜。由于夹卷作用的增强,浮力和风切变共同驱动的对流边界层较浮力驱动的对流边界层暖。在夹卷层,浮力和风切变共同驱动的边界层对流的上升气流和下沉气流都比浮力驱动的边界层对流中的强,而且垂直速度的概率密度函数分布也较对称,其位温和示踪物浓度的概率密度函数分布也比浮力驱动的边界层中的平直。对湍流动能收支的分析也表明风切变对湍流动能有重要影响,尤其对夹卷层中的湍流动能切变产生项影响较大。示踪物浓度的概率密度函数垂直分布显示,浮力驱动的边界层中示踪物浓度随高度变化较小,而浮力和风切变共同驱动的边界层中示踪物浓度随高度递减,但是示踪物传输的高度比较高。     

Organization of tropical convection in a GCM with varying vertical resolution; implications for the simulation of the Madden-Julian Oscillation

 Experiments using a GCM with two different vertical resolutions show differences in the amount of variability in the tropical upper tropospheric zonal wind component associated with the Madden-Julian Oscillation (MJO). The GCM with lower vertical resolution shows very little variability in this quantity whereas when the vertical resolution is doubled in the free troposphere, the GCM produces variability which is of the same strength as observations. However, the eastward propagation of an enhanced convective region from the Indian Ocean into the west Pacific is not well represented in either simulation of this atmospheric GCM. A water-covered or “aqua-planet” version of the same GCM is used to investigate the behaviour of tropical convection when the vertical resolution is doubled. When the vertical resolution is increased, the spectrum of tropical cloud types changes from a bimodal distribution with peaks representing shallow cumulus and deep cumulonimbus clouds to a trimodal distribution with a third peak in mid-troposphere near the melting level. Associated with periods when these mid-level congestus clouds are dominant, the detrainment from these clouds significantly moistens the mid-troposphere. The appearance of these congestus clouds is shown to be partly due to improved resolution of the freezing level and the convective processes occurring at this level. However, due to the way in which convective detrainment is parametrized in this model, the vertical profile becomes rather noisy and this too contributes to the change in the nature of the convective clouds. The resulting cloud distribution more closely resembles observations, particularly during the suppressed phase of the MJO when cumulus congestus is the dominant cloud type. Received: 17 April 2000 / Accepted: 30 November 2000     

北京城市化对一次降雪过程影响的数值模拟研究

利用中尺度数值预报（Weather Research and Forecast,WRF）模式,针对2018年3月17日05—17时（北京时）北京地区的一次降雪过程模拟分析了城市化对降雪的主要影响机制。结果表明,城市化使得北京五环以内降雪量减少,降雨量增加,这主要是由于城市化低层增温效应加强了雪的融化过程,产生混合型降水,距离市中心越近越容易发生混合型降水。城市化对降雪的总降水量和降水的时、空分布也存在一定的影响。降水初期,城市化造成的“城市干热岛”效应不利于水汽的水平和垂直输送,不利于云的形成,地面总降水量减小。随着降水过程的发展,部分冰相粒子融化,使近地面水汽增多,“城市热岛效应”的热力抬升作用有利于水汽的垂直输送和云的发展,部分云滴或水汽抬升进入云中,增强冷云过程,使雪和霰粒子含量增大,地面总降水量增加。城市化产生的“城市效应”对低层大气温度和云微物理过程产生影响,而云微物理过程的非绝热过程反过来又影响低层大气温度和大气层结,影响能量和水汽输送,进而对云和地面降水产生影响。      

基于探空云识别方法的云垂直结构分布特征

云的垂直结构特征作为云重要的宏观特征之一,直接决定了云的类型,进而通过发射和吸收辐射的方式影响着地气系统的能量收支平衡,因此对云垂直结构特征的研究一直都是云物理研究的一个重要方向。作为观测云垂直结构特征的一种方式,探空气球通过获取沿路径方向高分辨率的廓线信息,采用一定反演方法从而能够较为准确的识别云的垂直结构。本文即利用我国业务布网探空站的观测资料,采用相对湿度阈值法识别云垂直结构,并同激光云高仪、“风云四号”静止卫星和毫米波云雷达对识别的云结构特征量进行了一致性检验。在此基础上,统计分析了2015～2017年单层、两层和三层云的垂直结构分布特征、日变化和季节变化特征以及全国区域分布特征,结果表明:（1）整体分布上,单层云在垂直方向上出现的高度范围介于多层云的高度范围内,并且随着云层数的增加,云在垂直方向上更为伸展,即高层云越高,低层云越低;（2）在日变化中,中午单层和多层云中最低层云的云底高度均高于早晨,而夜间单层和多层云中最高层云的云顶高度则高于早晨和中午,同时中间层云厚的变化要小于最上层和最下层云厚的变化;（3）在季节变化中,夏季云量较其他季节更多,云体发展也更为深厚,表明温暖的大气条件更有利于云的形成和发展;（4）我国云垂直结构分布特征具有明显的纬向变化趋势,从以青藏高原为中心的西南地区的云底较高云体较薄的云,逐步过渡到以东南沿海地区为中心的云底较低云体较为深厚的云,表明不同地形和气候带的差异与不同云类型的分布直接相关。     

The vertical transport of air pollutants by convective clouds. Part III: Transport features of different cloud systems

The vertical transport features of gaseous pollutants, with a negative exponent profile of concentration, by dif-ferent types of convective cloud systems are numerically investigated by using a two-dimensional, reactionless convective cloud transport model. The results show that an isolated, weak storm is able to pump pollutant gas out PBL and transport it to the mid-troposphere, whereas a deep, intense thunderstorm can very efficiently transport air pollutants up to the mid and upper troposphere and laterally spread with anvil, forming an extensive concentration surge layer at altitude of ten-odd kilometers altitude. Each type of convective transport results in concentration re-duction in PBL. In a wind shear environment the transport efficiency of deep thunderstorm significantly increases and the pollutants enter into clouds on the downshear side at low-level and spread downwind in anvil layer. On the other hand, for a cumulus cloud with plenty of liquid water, the gas dissolution effect is increased, and the irreversible aqueous reactions, in extreme, may significantly weaken the vertical transports of pollutant gases even with solubility coefficients no more than 103 M atm-1.     

ATTILA: atmospheric tracer transport in a Lagrangian model

The model ATTILA has been developed to treat the global-scale transport of passive trace species in the atmosphere within the framework of a general circulation model (GCM). ATTILA runs online within the GCM ECHAM4 and advects the centroids of 80.000 to 190.000 constant mass air parcels. Each trace constituent is thereby represented by a mass mixing ratio in each parcel. ATTILA contains state-of-the-art parameterizations of convection, turbulent boundary layer mixing and inter-parcel transport, and provides an algorithm to map the tracer concentrations from the trajectories to the ECHAM model grid. The transport characteristics of ATTILA are evaluated against observations and the standard semi-Lagrangian transport scheme of ECHAM by two experiments. (1) We simulate the distribution of the short-lived tracer radon (²²²Rn) in order to examine fast vertical transport over continents, and long-range transport from the continents to remote areas. (2) We simulate the distribution of radiocarbon (¹⁴C) from nuclear weapon tests in order to examine upper tropospheric and stratospheric transport characteristics. Contrary to the semi-Lagrangian scheme, ATTILA shows a greatly reduced meridional transport in the upper troposphere and lower stratosphere, and a reduced downward flux from the stratosphere to the troposphere, especially in mid-latitudes. Since ATTILA is a numerically non-diffusive scheme, it is able to maintain steep gradients, which compare better to the observations than the rather smooth gradients produced by the semi-Lagrangian scheme.     

边界层对流对示踪物抬升和传输影响的大涡模拟研究

利用"西北干旱区陆气相互作用野外观测实验"加密观测期间敦煌站的实测资料以及大涡模式, 通过一系列改变地表热通量和风切变的敏感性数值试验, 分析了地表热通量和风切变对边界层对流的强度、形式, 以及对对流边界层结构和发展的影响。模拟结果显示风切变一定, 增大地表热通量时, 由于近地层湍流运动增强, 向上输送的热量也较多, 使对流边界层变暖增厚, 而且边界层对流的强度明显增强, 对流泡发展的高度也较高。当地表热通量一定, 增大风切变时, 由于风切变使夹卷作用增强, 将逆温层中的暖空气向下卷入混合层中, 使对流边界层增暖增厚, 但是对流泡容易破碎, 对流的强度也较弱。另外通过在模式近地层释放绝对浓度为100的被动示踪物方法, 用最小二乘法定量地分析了地表热通量和风切变分别与示踪物抬升效率和传输高度的关系。分析结果表明, 风切变小于10.5×10-3 s-1时, 增大地表热通量加强了上层动量的下传, 使示踪物的抬升效率也线性增大;地表热通量小于462.5 W m-2时, 增大风切变减弱了边界层对流的强度, 从而使示踪物的抬升效率减弱。当风切变一定时, 示踪物的平均传输高度随地表热通量增加而增大, 而地表热通量一定, 只有风切变大于临界值时, 示踪物平均传输高度才随风切变的增加而增大, 而临界风速的大小由地表热通量决定。     

Precipitation Structures of a Thermal Convective System Happened in the Central Western Subtropical Pacific Anticyclone

In this paper, characteristics of precipitating clouds in a thermal convective system (TCS) occurred in the southeastern mainland of China at 15:00 BT (Beijing time) on August 2, 2003 in the central western subtropical Pacific anticyclone (WSPA) is studied by using TRMM tropical rainfallmeasuring mission PR (precipitution radar) and IR Infrared radiation measurements. The precipitating cloud structures in both horizontal and vertical, relationship among storm top, cloud top, and surface rain rate are particularly analyzed. Results show that a strong ascending air at 500 hPa and a strong convergence of moisture flux at 850 hPa in the central WSPA supply necessary conditions both in dynamics and moisture for the happening of the TCS precipitation. The TRMM PR observation shows that the horizontal scale of the most TCS precipitating clouds is about 30-40 km, their averaged vertical scale is above 10 km, and the maximum reaches 17.5 km. The maximum rain rate near surface of those TCS clouds is beyond 50 mm h-1. The mean rain profile of the TCS clouds shows that its maximum rain rate at 5 km altitude is 1 km lower than the estimated freezing level of the environment. Compared with the mesoscale convective system (MCS) of "98.7.20", both systems have the same altitude of the maximum rain rate displayed from both mean rain profiles, but the TCS is much deeper than the MCS. From the altitude of the maximum rain rate to near surface, profiles show that rain rate reducing in the TCS is faster than that in the MCS, which implies a strong droplet evaporation process occurring in the TCS. Relationship among cloud top, storm top, and surface rain rate analysis indicates a large variation of cloud top when storm top is lower. On the contrary, the higher the storm top, the more consistent both cloud top and storm top. And, the larger the surface rain rate, the higher and more consistent for both cloud top and storm top. At the end, results expose that area fractions of non-precipitating clouds and clear sky are 86% and 2%, respectively. The area fraction of precipitating clouds is only about 1/8 that of non-precipitating clouds.