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

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

A Homogeneous Langevin Equation Model, Part ii: Simulation of Dispersion in the Convective Boundary Layer

We present a Lagrangian stochastic model of vertical dispersion in the convective boundary layer (CBL). This model is based on a generalized Langevin equation that uses the simplifying assumption that the skewed vertical velocity probability distribution is spatially homogeneous. This approach has been shown to account for two key properties of CBL turbulence associated with large-scale coherent turbulent structures: skewed vertical velocity distributions and long velocity correlation time. A 'linear-skewed' form of the generalized Langevin equation is used, which has a linear (in velocity) deterministic acceleration and a skewed random acceleration. 'Reflection' boundary conditions for selecting a new velocity for a particle that encounters a boundary were investigated, including alternatives to the standard assumption that the magnitudes of the particle incident and reflected velocities are positively correlated. Model simulations were tested using cases for which exact, analytic statistical properties of particle velocity and position are known, i.e., well-mixed spatial and velocity distributions. Simulations of laboratory experiments of CBL dispersion show that (1) the homogeneous linear-skewed Langevin equation model (as well as an alternative 'nonlinear-Gaussian' Langevin equation model) can simulate the important aspects of dispersion in the CBL, and (2) a negatively-correlated-speed reflection boundary condition simulates the observed dispersion of material near the surface in the CBL significantly better than alternative reflection boundary conditions. The homogeneous linear-skewed Langevin equation model has the advantage that it is computationally more efficient than the homogeneous nonlinear-Gaussian Langevin equation model, and considerably more efficient than inhomogeneous Langevin equation models.     

Dispersion Coefficients for Gaussian Puff Models

The Gaussian distribution is a good approximation for transient (instantaneously released) puff concentration distributions within a short period of time after release. Artificial neural network (ANN) models for puff dispersion coefficients were developed, based on observations from field experiments covering a wide range of meteorological conditions (in March, May, August and November). Their average predictions were in very good agreement with measurements, having high correlation coefficients (r > 0.99). A non-linear multi-variable regression model for dispersion coefficients was also developed, under the assumption that puff dispersion coefficients increase with time, and follow power laws. Both ANN-based and multi-regression non-linear models were able to use easily measured atmospheric parameters directly, without the necessity of predefining the Pasquill stability category. Predictions of ANN-based and multi-regression-based Gaussian puff models were compared with those of Gaussian puff models using Slade’s dispersion coefficients and COMBIC, a sophisticated model based on Gaussian distributions. Predictions from our two new models showed better agreement with concentration measurements than the other Gaussian puff models, by having a much higher fraction within a factor of two of measured values, and lower normalized mean square errors.     

On the influence of non-Gaussianity on turbulent transport

Non-Gaussianity effects, first of all the influence of the third and fourth moments of the velocity probability density function, have to be assessed for higher-order closure models of turbulence and Lagrangian modelling of turbulent dispersion in complex flows. Whereas the role and the effects of the third moments are relatively well understood as essential for the explanation of specific observed features of the fully developed convective boundary layer, there are indications that the fourth moments may also be important, but little is known about these moments. Therefore, the effects of non-Gaussianity are considered for the turbulent motion of particles in non-neutral flows without fully developed convection, where the influence of the fourth moments may be expected to be particularly essential. The transport properties of these flows can be characterized by a diffusion coefficient which reflects these effects. It is shown, for different vertical velocity distributions, that the intensity of turbulent transport may be enhanced remarkably by non-Gaussianity. The diffusion coefficient is given as a modification of the Gaussian diffusivity, and this modifying factor is found to be determined to a very good approximation by the normalized fourth moment of the vertical velocity distribution function. This provides better insight into the effect of fourth moments and explains the varying importance of third and fourth moments in different flows.     

Proper Orthogonal Decomposition of Mesoscale Vertical Velocity in the Convective Boundary Layer

The proper orthogonal decomposition technique is used to analyze wind-profiler observations from the surface to the top of the convective boundary layer. The mesoscale thermal structures are identified by decomposing the vertical velocity measurements into a sequence of characteristic modes with random coefficients. Results show that the first two dominant modes contribute over 85 % to the average kinetic energy. These most energetic modes also show mixed-layer similarity, which indicates that the non-local static instability plays a major role in determining the structure of the most energetic modes. Reconstruction of the vertical wind profiles by the first two dominant modes shows that they represent the most significant thermal structures. The probability distributions of the random coefficients related to these first two dominant modes are also analyzed and found to be Gaussian.     

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.     

Non-Gaussian Lagrangian Stochastic Model for Wind Field Simulation in the Surface Layer

Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality, gene flow(through pollen drift), and plant disease transmission(spore dispersion). Although Lagrangian stochastic(LS) models describe stochastic wind behaviors, such models assume that wind velocities follow Gaussian distributions. However,measured surface-layer wind velocities show a strong skewness and kurtosis. This paper presents an improved model, a non-Gaussian LS model, which inco...     

Turbulent velocity distributions and implied trajectory models

Well-mixed, first-order Lagrangian stochastic (LS) particle trajectory models are derived from several idealized (“toy”) turbulent velocity distributions, and their performance is compared against the observations of Project Prairie Grass, i.e., the case of a continuous point source of tracer near the ground, in the horizontally homogeneous and neutrally stratified surface layer. Although in a context of limited information a Gaussian distribution is the preferred choice, and although the Gaussian corresponds to the simplest of this set of LS models (namely, the Langevin equation), models stemming from other velocity distributions give similar, albeit distinguishable, predictions.     

On the Formulation of Lagrangian Stochastic Models of Scalar Dispersion Within Plant Canopies

Lagrangian stochastic models, quadratic in velocity and satisfying the well-mixed condition for two-dimensional Gaussian turbulence, are used to make predictions of scalar dispersion within a model plant canopy. The non-uniqueness associated with satisfaction of the well-mixed condition is shown to be non-trivial (i.e. different models produce different predictions for scalar dispersion). The best agreement between measured and predicted mean concentrations of scalars is shown to be obtained with a small sub-class of optimal models. This sub-class of optimal models includes Thomson's model (J. Fluid Mech. 180, 529–556, 1987), the simplest model that satisfies the well-mixed condition for Gaussian turbulence, but does not include two other models identified recently as being in optimal agreement with the measured spread of tracers in a neutral boundary layer. It is therefore demonstrated that such models are not universal, i.e. applicable to a wide range of flows without readjustment of model parameters. Predictions for scalar dispersion in the model plant canopy are also obtained using the model of Flesch and Wilson (Boundary-Layer Meteorol. 61, 349–374, 1992). It is shown that, when used with a Gaussian velocity distribution or a maximum-missing-information velocity distribution, which accounts for the measured skewness and kurtosis of velocity statistics, the agreement between predictions obtained using the model of Flesch and Wilson and measurements is as good as that obtained using Thomson's model.     

Measurements of Turbulence and Dispersion in Three Idealized Urban Canopies with Different Aspect Ratios and Comparisons with a Gaussian Plume Model

Water-tunnel measurements of velocity, turbulence and scalar concentration for three model urban canopies with aspect ratios A _r of building height-to-width of 0.25, 1 and 3 are presented. The measurements for the canopies with A _r = 1 and 3 are new, while the measurements for A _r = 0.25 were previously published. A passive scalar was continuously released from a near-ground point source, and the concentration was measured at several distances from the source and at different heights above the ground. Plume spreads, concentration and distance from the source were non-dimensionalized using length, time and velocity scales reflecting the geometry of the buildings. The scaling collapses the data for all aspect ratios and is valid when the vertical extent of the plume is smaller than the canopy height. The observed plume spreads are compared with analytical relations, which predict linear growth in both transverse and vertical directions. The observed mean concentration is compared with a Gaussian dispersion model that predicts a ?2 power-law decay with distance from the source.     

Airborne measurements of the vertical flux of ozone in the boundary layer

A fast-response chemiluminescent ozone sensor was mounted in an aircraft instrumented for air motion and temperature measurements. Measurements of the vertical flux of ozone by the eddy correlation technique were obtained after correcting for time delay and pressure sensitivity in the ozone sensor output. The observations were taken over eastern Colorado for two days in April, one a morning and the other an afternoon flight. Since the correlation coefficient of ozone and vertical velocity is small compared to, for example, temperature and vertical velocity in the lower part of the convective boundary layer, an averaging length of the order of 100 km was required to obtain a reasonably accurate estimate of the ozone flux. The measured variance of ozone appeared to be too large, probably mainly due to random noise in the sensor output, although the possibility of the production of ozone fluctuations by chemical reactions cannot be dismissed entirely. Terms in the budget equation for ozone were estimated from the aircraft measurements and the divergence of the ozone flux was found to be large compared to the flux at the surface divided by the boundary-layer height.Part of this work was carried out while a visitor at Risø National Laboratory, Denmark.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Wind-tunnel Modelling of Dispersion from a Scalar Area Source in Urban-Like Roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ _p = λ _f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w _T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness.     

Simulation of Time Series of Concentration Fluctuations in Atmospheric Dispersion Using a Correlation-distortion Technique

Short-duration fluctuations in the concentration of airborne substances can be important in a variety of atmospheric dispersion problems, especially when assessing the risks posed by harmful materials. This paper discusses a simulation technique for generating individual realisations of fluctuating concentration time series in dispersing plumes based on target probability distributions and spectral functions. The scheme uses a correlation-distortion approach to simulate these time series. Gaussian processes with modified spectral characteristics are generated and then transformed to yield non-Gaussian processes with the desired spectral characteristics. The simulation scheme is initially developed for a single receptor, and is then extended to model pairs of correlated time series at two receptors. In fact, the modelling technique can be generalised to an arbitrary number of receptors and this provides, in principal, an approach that is applicable to a wide class of similar problems (such as the modelling of instantaneous puff releases or the response of line-of-sight detection systems). The simulation technique is illustrated using observations made during recent field experiments, conducted both in the United Kingdom and in the U.S.A., investigating the short-range dispersion of a passive tracer.     

Wave-Modified Flux and Plume Dispersion in the Stable Boundary Layer

The effects of a pressure jump and a following internal gravity wave on turbulence and plume diffusion in the stable planetary boundary layer are examined. The pressure jump was accompanied by a sudden increase in turbulence and plume dispersion. The effects of wave perturbations on turbulence statistics are analysed by calculating fluxes and variances with and without the wave signal for averaging times ranging from 1 to 30 min. The wave signals are obtained using a band-pass filter. It is shown that second-order turbulence quantities calculated without first subtracting the wave perturbations from the time are greater than those calculated when the wave signal is separated from the turbulence. Estimates of the vertical dispersion of an elevated tracer plume in the stable boundary layer are made using an elastic backscatter lidar. Plume dispersion observed 25 m downwind of the source increases rapidly with the arrival of the flow disturbances. Measured plume dispersion and plume centreline height correlate with the standard deviation of the vertical velocity but not with the wave signal.     

Multi-beam techniques for deriving wind fields from airborne doppler radars

Summary Two techniques for deriving horizontal and vertical air motions using vertically scanning airborne Doppler radar data are presented and discussed. These techniques make use of the scanning ability of the NOAA P-3 tail-mounted radar antenna to view a region of space from at least two vantage points during a straight-line flight track. The scanning methodology is termed the Fore/Aft Scanning Technique or FAST because the antenna is alternately scanning forward and then aft of the flight track. The major advantages of FAST over flying two quasi-orthogonal flight tracks with the antenna scanning normal to the flight track are that the data are collected in roughly half the time and the aircraft does not have to execute a right-angle turn. However, accuracy of the resulting wind field is compromised slightly because the beam intersection angle is reduced from 90° to about 50°. The reduction of area covered because of large drift angles is also discussed.A three-dimensional wind field can be constructed using the dual-Doppler equations from FAST data using the two radial velocity estimates and vertical integration of the continuity equation with a boundary condition of no vertical motion at cloud top and the Earth's surface. To keep errors in the calculated winds acceptably small, the elevation angles are typically restricted to ±45° from the horizontal to minimize contamination of the horizontal wind by terminal fallspeeds.A different, and perhaps more believable vertical velocity, can be derived using a second technique that utilizes two (or more) airborne Doppler radar equipped aircraft each using FAST to observe the echo-top vertical velocity at common point (e.g., two aircraft flying parallel flight, paths, or by using an L-shaped flight track with a single aircraft). This technique results in 4 (or more) radial velocity estimates at each point (hence is called the quad-Doppler technique). Horizontal winds can be derived using either an overdetermined three-equation solution or an overdetermined dual-Doppler solution, whichever is more accurate. For the calculation of vertical velocity a new approach is proposed that utilizes the overdetermined triple-Doppler solution for vertical particle motion near cloud top, minus an estimate of terminal fallspeeds, as a top boundary condition for the downward vertical divergence integration to derive vertical air velocity elsewhere in the domain. In addition, this approach allows measurements at steep elevation angles allowing for more depth of coverage for a given range.To show the utility of the method, analyses of data collected using FAST are compared to conventional dual-Doppler-derived wind fields constructed from data collected simultaneously by S-band ground-based Doppler radars. An example of the quad-Doppler technique is also presented from the recently completed Tropical Oceans/Global Atmospheres Coupled Ocean/Atmosphere Response Experiment (TOGA/COARE). Comparisons of quad-Doppler vertical velocity are made with in-situ derived vertical air motions collected by the NASA DC-8 to judge the quality of the approach.With 9 Figures     

全球海洋CFC-11吸收对传输速度的敏感性

本文在中国科学院大气物理研究所发展的全球海洋模式(LICOM)中使用五个不同的海气交换的气体 传输速度公式对CFC-11(一氟三氯甲烷)在海洋中的分布和吸收做了模拟。讨论了不同气体传输速度的差异, 发现差异最大的两个公式得到的全球年平均传输速度相差81%。对CFC-11的海表浓度分布、海气通量、水柱总量、海水累积吸收量以及在大洋断面的垂直浓度分布进行了分析。分析结果显示, 使用Liss and Merlivat (1986) 的传输速度公式的试验在海气通量和海洋吸收总量的模拟上均小于其他试验, Nightingale et al. (2000)、Ho et al. (2006) 和Sweeney et al. (2007) 等的公式虽然全球年平均值相近, 但在高风速地区Nightingale et al. (2000) 公式的传输速度要小于后两者, 导致了使用该公式的试验模拟结果在主吸收区和存储区的强度比后二者偏小。Wanninkhof (1992) 的公式在形式上与Ho et al. (2006) 以及Sweeney et al. (2007) 的公式一致, 但在系数上存在差别, 这使得模拟的水柱总量在南大洋的分布明显好于其他试验, 尽管其最大值仍比观测资料略小。在海洋累积吸收量的计算上, 使用Wanninkhof (1992) 传输速度公式得到的模拟结果比观测资料小8%左右。计算了Liss and Merlivat (1986) 和Wanninkhof (1992) 的传输速度公式模拟的单年吸收量相对差, 其总体上一直保持持续下降的趋势, 到2007年仅为2%。从该相对差变化趋势看, 在最初的前10年, 海气CFC-11交换通量对海气交换传输速度的敏感性更强, 而在更长时间的模拟上, 海洋对CFC-11的吸收则更依赖于物理模式的通风速率。通过对CFC-11垂直断面分布的分析可知, 不同的传输速度在主要吸收区的不同导致了一定的垂直分布差异。基于本文的结果可以认为Wanninkhof (1992) 的海气气体传输速度公式更适合本模式对CFC-11的模拟。     

Validation of very high cloud droplet number concentrations in air masses transported thousands of kilometres over the ocean

The microstructure of orographic clouds related to the aerosol present was studied during the second Aerosol Characterisation Experiment (ACE‐2). Very high cloud droplet number concentrations (almost 3000 cm⁻³) were observed. These high concentrations occurred when clouds formed on a hill slope at Tenerife in polluted air masses originating in Europe that had transported the order of 1000 km over the Atlantic Ocean. The validity of the measured droplet number concentrations was investigated by comparing with measurements of the aerosol upstream of the cloud and cloud interstitial aerosol. Guided by distributions of the ratios between the measurements, three criteria of typically 30% in maximum deviation were applied to the measurements to test their validity. Agreement was found for 88% of the cases. The validated data set spans droplet number concentrations of 150–3000 cm⁻³. The updraught velocity during the cloud formation was estimated to 2.2 m s⁻¹ by model calculations, which is typical of cumuliform clouds. The results of the present study are discussed in relation to cloud droplet number concentrations previously reported in the literature. The importance of promoting the mechanistic understanding of the aerosol/cloud interaction and the use of validation procedures of cloud microphysical parameters is stressed in relation to the assessment of the indirect climatic effect of aerosols.     

On the use of forward scatter techniques in the study of turbulent stratified layers in the troposhere

This review paper considers the potentials of forward scatter radio techniques as a diagnostic tool in the study of stratified turbulent layers. Seven classes of experimental techniques have been considered (beamswinging, multifrequency, vertical and horizontal field strength correlation, coupling loss, pulse delay and bandwidth measurements). All these techniques with the exception of bandwidth and delay experiments measure directly or indirectly the irregularity spectrum of refractive index. Bandwidt hand delay distributions are determined by the vertical distribution of the refractivity irregularities. Thus for the purpose of obtaining information about a layered structure where our interest is focused on the vertical refractivity distribution rather than on the averaged spectral properties of the structure, the last two categories of experiments appear to represent a powerful tool. If one is interested in the spectral properties, multifrequency and beamswinging experiments are of great potential value whereas other experiments such as that measuring coupling loss are largely influenced by pure refraction effects.     

超声波传感器雪深测量与人工观测对比试验分析

通过对2010年2、3月佳木斯国家基准气候观测站人工观测与超声波雪深传感器测量获取数据对比、分析不同天气条件下（高风速、低风速、低温度）两种方法获得的数据差值的变化情况。结果表明：观测结果差异的大小主要受温度和风速两个因素的影响。低温环境下温度会影响超声波的行程时间使得超声波传感器的测量精度受到影响,通过温度补偿的方法对雪深进行订正可以提高超声波传感器的测量精度,两者的差值较小,且观测资料有较好的可比性;高风速时,风速会影响超声波脉冲,使其偏离传感器的下方,同时超声波的速度也会受到影响,从而影响测量结果,使得两者的差值较大。