对流边界层中的一个双高斯型PDF扩散模式

对流边界层（ＣＢＬ）中的污染扩散是非高斯型的。本文在下列三个假设下建立了双高斯型几率分布函数（ＰＤＦ）模式：１．对流边界层任一确定高度铅直速度Ｗ的几率分布函数ｐｗ由两个高斯分布迭加而成；２．从污染源释放的粒子具有源高的铅直速度几率分布，且其轨迹是线性的；３．粒子在地面的反射为全反射，在混合层高度Ｚｉ为全反射或有部分吸收。然后分析了三个高度上铅直速度Ｗ的一些统计特征量，比较了由ＰＤＦ模式计算的横风向积分浓度和Ｌａｍｂ的数值模拟［１－３］，Ｄｅａｒｄｏｒｆｆ的水槽模拟［４－６］结果，并用美国ＣＯＮＤＯＲＳ计划的外场试验资料［７］对ＰＤＦ模式进行了验证，结果均相当一致。     

云南安宁市辐射雾大气扩散参数研究

利用三向风标探测资料，分析了安宁县辐射雾中的扩散参数，结果表明：平均风速对水平风向脉动标准差起着重要作用，是影响水平湍流扩散的关键因子；雾中的横向和铅直方向扩散参数分别相当于ＰａｓｑｕｉｌｌＢ类和Ｄ类稳定度下的扩散参数，表明在辐射雾出现期间，虽然铅直方向上的扩散能力不强，但有较强的水平扩散能力。     

基于Weierstrass-Mandelbrot函数的分形风速脉动仿真

实际风速脉动普遍具有自相似分形特征,而传统的谐波合成法和线性滤波法仿真的风速脉动均不具有自相似分形特征.因此,基于随机型Weierstrass-Mandelbrot函数,设计了一种能够仿真自相似风速脉动的方案.其中,表征风速脉动自相似特征的重要参数分形维度可与湍流惯性区能谱的幂指数建立联系.将该方案仿真的风速脉动与实际风速脉动一些重要的统计特征,如功率谱和概率密度函数等,进行了比较,结果表明提出的新方案能有效仿真风速脉动的中高频变化及其概率分布特征.     

中尺度数值模式中湿过程参数化的敏感性试验

本文利用一个三维中尺度模式，分别采用两种不同的湿过程处理方案ＦＰＡ和ＮＣＡ模拟同一个洋面冷锋个例，以考察模拟结果对湿过程参数化方案的敏感性。两种方案的区别在于前一种使用了次网格湿对流参数化而后者没有。通过比较分别使用ＦＰＡ和ＮＣＡ方案所进行的模拟ＦＰＳ和ＮＣＳ，发现冷锋结构在两者间出现较大差异。从８５０ hPa等压面|Δθe|等值线图所显示的结构看，ＮＣＳ中冷锋呈现为一长一短两条锋带，而ＦＰＳ中冷锋仅呈现为一条锋带。在相应的冷锋横向垂直剖面中，ＮＣＳ中存在两个相邻的锋面垂直环流，而ＦＰＳ中只有一个。另外，ＮＣＳ中对流不稳定的冷锋区出现垂直运动过分发展的现象，而ＦＰＳ中不存在这个问题。通过与卫星云图比较，我们注意到，ＮＣＳ中尽管出现了垂直运动过分发展的现象，但相对于ＦＰＳ而言，其模拟的冷锋结构与实际较为接近，成功地模拟出云图上显示的双冷锋结构。ＮＣＳ中垂直运动过分发展的一个重要原因是模式中采用的静力平衡近似。ＦＰＳ中次网格对流参数化方案的使用，通过减少甚至消除对流不稳定度，一方面使垂直运动的过分发展倾向受到约束，另一方面也可能使在对流不稳定层结下的锋区环境中本应出现的中尺度结构失去了形成的机会。     

下曳气流在暴洪中的作用

介绍了Ｓｐｅｎｃｅｒ Ｐ Ｌ和Ｓｔｅｎｓｒｕｄ Ｄ Ｊ关于下曳气流在暴洪中作用的模拟结果。Ｓｐｅｎｃｅｒ和Ｓｔｅｎｓｒｕｄ利用Ｋａｉｎ－Ｆｒｉｓｃｈ对流参数化方案及其３种修改方案即最大降水效率方案（ＭＰＥ）、无下曳气流方案（ＮＤＤ）和滞后下曳气流方案（ＤＤＤ）对暴洪个例进行模拟试验。结果表明：利用种修改了的对流参数化方案模拟出的２４ｈ最大降水量的结果与未修改的最大降水量均比实测的小，但是ＤＤＤ所模     

二阶闭合模式在污染气象学中的某些应用

本文采用二阶闭合的湍流边界层模式,进行一系列数值试验以模拟边界层中连续线源的扩散状况。试验表明:无论在稳定的或不稳定的边界层中,高源的扩散能力都低于低源;在稳定层中,粗糙地表上的大气扩散能力高于光滑表面;在相同风速和地表净辐射情况下,粗糙表面上的大气扩散能力反而低于光滑表面;对流边界层中存在反梯度输送,因而K理论的应用受到限制。试验还表明,修正的Kazanski-Monin参数可能比Monin-Obukhov长度更能反映大气的扩散能力。     

边界层大气环境与空气污染问题的风洞流体模拟实验

在大气环境风洞中,采用人工形成边界层的方法,建立一定厚度的模拟边界层,测试确定边界层特征量,如粗糙度长z。风廓线指数n等。实验揭示了形成大气边界层的物理机制。用烟云照相法对模拟形成的高架点源排放烟流的抬升高度相大气扩散参数作了实验测定,并采用波动烟流法进行数据处理。讨论了排烟速度Vs和环境风速u对烟气动力抬升的影响;此外还对阶跃下垫面条件和障碍物后旋涡区中的扩散进行了分析。     

中尺度对流雨带飞机观测资料分析和数值模拟研究

本文使用MM5中尺度数值模式，采用不同分辨率(90、45、30和15Km)相同的Reisner混合相微物理显示方案，及Betts—Miller和Grell的对流参数化方案，对2003年7月8日吉林省境内中尺度对流雨带天气过程进行了数值模拟试验。结果表明：不同尺度的天气系统，需要不同分辨率、不同参数化方案的中尺度模式进行模拟，提高模式的分辨率可以增强对含有对流的中小尺度天气系统的模拟和预报能力；模式的水平分辨率和对流参数化方案对模拟强降水中心有重要影响；当MM5模式分辨率提高时，模拟的细化给降水分布、降水强度带来一些改进，但主要雨区内也出现了一些虚假预报中心水平扩散；低层正涡度区，高层负涡度区的配合，对低层辐合上升运动有利、能促进对流发展；从散度场的分布特征看。低层850hPa的水平辐合区域强弱与降水的大小有密切关系，强水平辐合区与强降水区有很好的对应关系。因此，在MM5模式业务化时应该根据天气特点来选择模式分辨率和对流参数化方案，以使模拟结果更接近实况。     

对流能量计算及强对流天气落区预报技术研究

文章分析了两种典型的大气湿绝热过程及其处理方法，对大气对流能量参数的计算技术进行了研究。结合实际个例，利用可逆饱和绝热过程，对包含液态水重力拖曳作用的修正对流有效位能(MCAPE)和修正下沉对流有效位能(MDCAPE)进行了定量计算。文章结合数值模式输出探空分析，预报不稳定和对流能量的区域分布，在此基础上建立了综合多指标叠套强对流天气落区预报方法，用MM5 及国家气象中心T106模式输出及诊断产品预报强对流天气落区，并检验强对流落区预报技术。     

SAMIL模式中Tiedtke 积云对流方案对热带降水模拟的影响

目前, 大多数全球耦合模式及大气环流模式在降水模拟中普遍存在不同程度的“热带偏差”问题, 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室所发展的全球大气环流谱模式SAMIL-R42L26也存在这一现象, 主要表现在SPCZ (南太平洋辐合区) 降水过强且过分东伸、 赤道附近降水偏少等方面。本文通过修改SAMIL中的积云对流方案有效地削弱了这一偏差, 并进一步探讨其原因, 发现对流方案修改后, 改变了对流层低层至地面的温度分布状态, 进而影响了风速及散度场的模拟, 最终通过垂直速度的调整反作用于对流过程。比较修改前后对流过程云底质量通量, 发现修改后的方案主要通过削弱浅对流来提高热带降水的模拟性能。     

ON THE PDF MODELS FOR ATMOSPHERIC DIFFUSION IN BOUNDARY LAYER

In this paper, taking its turbulent exchange coefficient as a function of the Lagrangian timescale and standard variance of the turbulence in atmosphere, the atmospheric dispersion PDFmodels are obtained on the basis of atmospheric diffusion K-theory. In the model the statistics ofwind speed are directly used as its parameters instead of classic dispersion parameters. The bi-Gaussian PDF is derived in convective boundary layer (CBL), from the statistics of verticalvelocity in both of the downdraft and updraft regions that are investigated theoretically in the otherpart of this paper. Giving the driven parameters of the CBL (including the convective velocity scalew* and the mixing depth h_i) and the time-averaged wind speed at release level, the PDF model isable to simulate the distribution of concentration released at any levels in the CBL. The PDF'ssimulations are fairly consistent with the measurements in CONDORS experiment or the resultsbrought out by some numerical simulations.     

Flux Footprints in the Convective Boundary Layer: Large-Eddy Simulation and Lagrangian Stochastic Modelling

We investigated the flux footprints of receptors at different heights in the convective boundary layer (CBL). The footprints were derived using a forward Lagrangian stochastic (LS) method coupled with the turbulent fields from a large-eddy simulation model. Crosswind-integrated flux footprints shown as a function of upstream distances and sensor heights in the CBL were derived and compared using two LS particle simulation methods: an instantaneous area release and a crosswind linear continuous release. We found that for almost all sensor heights in the CBL, a major positive flux footprint zone was located close to the sensor upstream, while a weak negative footprint zone was located further upstream, with the transition band in non-dimensional upwind distances −X between approximately 1.5 and 2.0. Two-dimensional (2D) flux footprints for a point sensor were also simulated. For a sensor height of 0.158 z _i, where z _i is the CBL depth, we found that a major positive flux footprint zone followed a weak negative zone in the upstream direction. Two even weaker positive zones were also present on either side of the footprint axis, where the latter was rotated slightly from the geostrophic wind direction. Using CBL scaling, the 2D footprint result was normalized to show the source areas and was applied to real parameters obtained using aircraft-based measurements. With a mean wind speed in the CBL of U = 5.1 m s⁻¹, convective velocity of w _* = 1.37 m s⁻¹, CBL depth of z _i = 1,000 m, and flight track height of 159 m above the surface, the total flux footprint contribution zone was estimated to range from about 0.1 to 4.5 km upstream, in the case where the wind was perpendicular to the flight track. When the wind was parallel to the flight track, the total footprint contribution zone covered approximately 0.5 km on one side and 0.8 km on the other side of the flight track.     

Vertical velocity field in the convective boundary layer as observed with an acoustic Doppler sodar

The morning development of the daytime convective boundary layer (CBL) during fine weather has been observed with an acoustic Doppler sodar of the C.R.P.E. In particular, the vertical profile of the vertical velocity third-order statistic W* ³ has been obtained. This quantity is a maximum near 0.3z _I where z _I, is the height of the CBL. The histogram of vertical velocity in the CBL shows a relationship between W ³ and the convective velocity W _* and is useful for convective plume determination.     

On the Determination of the Neutral Drag Coefficient in the Convective Boundary Layer

Based on the idea that free convection can be considered as a particular case of forced convection, where the gusts driven by the large-scale eddies are scaled with the Deardorff convective velocity scale, a new formulation for the neutral drag coefficient, C_Dn, in the convective boundary layer (CBL) is derived. It is shown that (i) a concept of C_Dn can still be used under strongly unstable conditions including a pure free-convection regime even when no logarithmic portion in the velocity profile exists; (ii) gustiness corrections must be applied for rational calculations of C_Dn; and (iii) the stratification function used in the derivation of C_Dn should satisfy the theoretical free-convection limit. The new formulation is compared with the traditional relationship for C_Dn, and data collected over the sea (during the Tropical Ocean-Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) and the San Clemente Ocean Probing Experiment (SCOPE)) and over land (during the BOREX-95 experiment) are used to illustrate the difference between the new and traditional formulations. Compared to the new approach, the traditional formulation strongly overestimates C_Dn and z_o in the CBL for mean wind speed less than about 2 m s^-1. The new approach also clarifies several contradictory results from earlier works. Some aspects related to an alternate definition of the neutral drag coefficient and the wind speed and the stress averaging procedure are considered.     

Estimation of the Exchange Coefficient of Heat During Low Wind Convective Conditions

For the first time, the exchange coefficient of heat C_H has been estimated from eddy correlation of velocity and virtual temperature fluctuations using sonic anemometer measurements made at low wind speeds over the monsoon land atJodhpur (26°18' N, 73°04' E), a semi arid station. It shows strong dependence on wind speed, increasing rapidly with decreasing wind speed, and scales according to a power law C_H = 0.025U₁₀ ^-0.7 (where U₁₀ is the mean wind speed at 10-m height). A similar but more rapid increase in the drag coefficient C_Dhas already been reported in an earlier study. Low winds (<4 m s^-1) are associated with both near neutral and strong unstable situations. It is noted that C_H increases with increasing instability. The present observations best describe a low wind convective regime as revealed in the scaling behaviour of drag, sensible heat flux and the non-dimensional temperature gradient. Neutral drag and heat cofficients,corrected using Monin–Obukhov (M–O) theory, show a more uniform behaviour at low wind speeds in convective conditions, when compared with the observed coefficients discussed in a coming paper.At low wind convective conditions, M-O theory is unable to capture the observed linear dependence of drag on wind speed, unlike during forced convections. The non-dimensional shear inferred from the present data shows noticeable deviations from Businger's formulation, a forced convection similarity. Heat flux is insensitive to drag associated with weak winds superposed on true free convection. With heat flux as the primary variable, definition of new velocity scales leads to a new drag parameterization scheme at low wind speeds during convective conditionsdiscussed in a coming paper.     

The computation of the friction velocity u * and the temperature scale T * from temperature and wind velocity profiles by least-square methods

A method is given to calculate the surface layer parameters: u _* (friction velocity) and T _* (temperature scale) from wind speed and temperature profiles.The problem is formulated as a minimization of a least-square function, which is constructed from the difference between the measured profiles and the well-known Kansas profile relations.The wind speed and temperature profiles are treated simultaneously in this procedure. All the available wind speed and temperature measurements are used in order to reduce the effect of measurement errors.Estimates of the goodness of fit and confidence limits on the estimated parameters are discussed.The method has been applied to data obtained during experiments in a wide variety of conditions: Project Prairie Grass, experiments over Lake Flevo and experiments at the meteorological tower at Cabauw, the last two in the Netherlands.     

A general radar equation for the bistatic acoustic sounder

A radar equation which can be applied to any bistatic acoustic sounder is derived. The equation reduces to the expression normally used for the special case of the monostatic sounder. Numerical results using this equation are given for specific acoustic sounders, including the relative contributions to the scattering from the temperature and velocity parameters (C _T ² and C _V ²), the effect of wind speed on the scattered intensity and the measurement of horizontal wind velocity.Now with the Department of Electrical and Electronic Engineering, Portsmouth Polytechnic, Anglesea Road, Portsmouth PO1 3DJ.     

Length Scales of Scalar Diffusion in the Convective Boundary Layer: Laboratory Observations

A laboratory study of scalar diffusion in the convective boundary layer has found results that are consistent with a 1999 large-eddy simulation (LES) study by Jonker, Duynkerke and Cuijpers. For bottom-up and top-down scalars (introduced as ‘infinite’ area sources of passive tracer at the surface and inversion, respectively) the dominant length scale was found to be much larger than the length scale for density fluctuations, the latter being equal to the boundary-layer depth h. The variance of the normalized passive scalar grew continuously with time and its magnitude was about 3–5 times larger for the top-down case than for the bottom-up case. The vertical profiles of the normalized passive scalar variance were found to be approximately constant through the convective boundary layer (CBL) with a value of about 3–8c_*² for bottom-up and 10–50c_*² for top-down diffusion. Finally, there was some evidence of a minimum in the variance and dominant length scale for scalar flux ratios (top-down to bottom-up flux) close to −0.5. All these convection tank results confirm the LES results and support the hypothesis that there is a distinct difference in behaviour between the dynamic and passive variables in the CBL.     

Dispersion of passive tracer in the atmospheric convective boundary layer with wind shears: a review of laboratory and numerical model studies

Summary Paper reviews recent laboratory and numerical model studies of passive gaseous tracer dispersion in the atmospheric convective boundary layer (CBL) with surface and elevated wind shears. Atmospheric measurement data used for validation of these two model techniques are briefly discussed as well. A historical overview is given of laboratory studies of dispersion in the atmospheric CBL. Model studies of tracer dispersion in two CBL types, the (i) non-steady, horizontally homogeneous CBL and (ii) quasi-stationary, horizontally heterogeneous CBL, are reviewed. The discussion is focused on the dispersion of non-buoyant plume emitted from a point source located at different elevations within the CBL. Approaches towards CBL modeling employed in different laboratory facilities (water tanks and wind tunnels) are described. The reviewed numerical techniques include Large Eddy Simulation (LES) and Lagrangian modeling. Numerical data on dispersion in the sheared CBL is analyzed in conjunction with experimental results from wind-tunnel CBLs.     

The Role of Shear in the Morning Transition Boundary Layer

We use large-eddy simulation (LES) to better define the early stages of the morning transition boundary layer. Previous LES studies relating to the morning transition boundary layer focus on the role of the entraining convective boundary layer (CBL). By using a combination of different domain sizes and grid lengths, the full evolution from the stable boundary layer (SBL) to the CBL is modelled here. In the early stages of the morning transition the boundary layer is shown to be a combination of a shallow mixed layer capped by a significant shear driven stable boundary layer (the so-called mixed CBL–SBL state). The mixed CBL–SBL state is the key to understanding the sensitivity to shear. Turbulent kinetic energy budgets also indicate that it is shear driven. The negative flux from the mixed CBL–SBL state extends much further above the minimum than is typically found for the CBL later in the day, and the depth of penetration scales as w _m /N _i , where w _m is the combined friction and convective velocity scale and N _i the static stability at the inversion top.