论边界层中的大气扩散PDF模式

基于大气扩散Ｋ理论，用作为风速脉动均方差和拉氏时间尺度函数的湍流交换系数，得到了直接利用风速脉动几率密度而不用扩散参数的大气扩散ＰＤＦ模式。分别研究了对流边界层上升气流区与下降区垂直速度的统计特征，求得双正态ＰＤＦ模式。在给定ＣＢＬ自身参数如对流特征速度ｗ＊，顶高ｈｉ和源高度上的平均风速时，该模式计算出的无量纲浓度分布与室内外测试结果一致。     

牙克石地区大气扩散试验分析

利用牙克石电厂野外平衡球观测资料分别计算了不同稳定度下的大气扩散参数,并与其他结果进行了比较,揭示了牙克石地区大气湍流扩散特征。试验结果仅适用于远距离的高架源排烟计算。     

沿海地区大气边界层扩散特征探讨

利用辽东湾近地层三轴风速仪资料分别计算不同稳定度的大气扩散参数及湍流强度。并将不同情况下的大气扩散参数及湍流强度进行了比较。结果表明 :在地形相对平坦的辽东湾 ,在不同的稳定度下 ,σy 和σz值均略大于同级布里格斯公式的计算值。在不稳定条件 (C、B类 )下 ,湍强值最大。     

沿海地区大气边界层扩散特征探讨

利用辽东湾近地层三轴风速仪资料分别计算不同稳定度的大气扩散参数及湍流强度。并将不同情况下的大气扩散参数及湍流强度进行了比较。结果表明：在地形相对平坦的辽东湾,在不同的稳定度下,σy和σx值均略大于同级布里格斯公式的计算值。在不稳定条件(C、B类)下,湍强值重大。     

辽东湾地区大气边界层扩散特征探讨

利用辽东湾地区近地层三轴风速仪资料,分别计算了不同稳定度下的大气扩散参数及湍流强度。并对不同情况下的大气扩散参数及湍流强度进行了比较。结果表明,在地形相对平坦的辽东湾地区,在不同的稳定度下,σy和σx值均略大于同级Briggs公式的计算值,但小于地形相对复杂的本溪地区。在不稳定条件(C、B类)下,湍强值最大。     

大气边界层热量输送的非局地多尺度湍流理论及试验研究

该文在大气边界层多尺度湍流理论和湍流穿越理论的基础上，建立了大气边界层感热通量的非局地多尺度湍流计算方法，并选用1998年5～8月份第二次青藏高原大气科学试验（TIPEX）的大气边界层综合观测基地昌都站观测资料，分别用感热通量的非局地多尺度湍流计算方法和经典相似理论计算感热通量。对计算结果的对比分析表明，用感热通量的非局地多尺度湍流计算方法能计算出感热通量的逆梯度输送，也能计算出青藏高原近中性层结条件下明显的感热通量。此外，多尺度湍流理论还有待进一步发展完善，提高感热通量的计算精度。     

玛纳斯电厂冬季大气扩散的试验研究

在进行大气环境质量的评价和预测时,对大气扩散规律的分析是十分重要的.目前普遍采用的方法之一是用观测平衡气球轨迹的方法计算扩散参数.本文对1989年冬季,在玛纳斯电厂用平衡气球轨迹法进行的大气湍流扩散试验资料作了分析计算,并与P-G(帕斯圭尔一吉福德)结果进行了比较,提出了厂区不同稳定度条件下的大气湍流扩散参数σ_y和σ_z,为研究和摸清该区大气污染物的扩散能力提供了依据.     

利用测温声雷达回波图估算大气扩散的一种参数化方法

本文提出了一种基于湍流扩散统计理论的利用声雷达回波图估算大气扩散参数的新方法。这种新方法,仅需要从声雷达回波图上读出边界层结构的信息并取得常规地面测风资料,即可利用参数化公式进行计算,不需进行Pasquill稳定度分类。 所应用的参数化公式用外场观测资料进行了检验,表明公式计算结果与观测值比较符合。这种方法计算得到的扩散参数也与P-G扩散曲线进行了比较。     

不同下垫面的大气边界层湍流扩散特征

利用海洋、戈壁、城郊和城市4种不同下垫面上的湍流观测资料,运用多尺度湍流概念和方法计算分析了风速谱和扩散参数等湍流特征量。结果表明:下垫面的粗糙度对大气边界层湍流特征影响很大,随着粗糙度的增大,速度谱的峰值频率向高频方向移动,扩散参数增大。     

城市扩散参数的综合测定与研究

利用常州市城市边界层风、温场及湍流场的综合观测资料,对该地区大气扩散稀释能力进行了全面系统的分析,总结出一组适合平原地区中等城市的扩散参数计算公式。研究表明,城市贴地层及高层的湍流扩散特征有明显差异,应该用不同的扩散参数计算近地面面源和高架点源的浓度分布。在此基础上,首次提出了按高度层次划分的城市扩散参数体系。实际监测资料的验证结果,表明该组公式的模式计算效果甚好。     

Along-wind dispersion of puffs released in a built-up urban area

The SF₆ gas tracer observations for puffs released near the ground during the Joint Urban 2003 (JU2003) urban dispersion experiment in Oklahoma City have been analysed. The JU2003 observations, at distances of about 100–1,100 m from the source, show that, at small times, when the puff is still within the built-up downtown domain, the standard deviation of the concentration time series, σ_t, is influenced by the initial puff spread due to buildings near the source and by hold-up in the wakes of large buildings at the sampler locations. This effect is parameterised by assuming an initial σ_to of about 42 s, leading to a comprehensive similarity formula: σ_t = 42 + 0.1t. The second term, 0.1t, is consistent with an earlier similarity relation, σ_t = 0.1t, derived from puff observations in many experiments over rural terrain. The along-wind dispersion coefficient, σ_x, is assumed to equal σ_t u, in which u is the puff speed calculated as the distance from the source to the sampler, x, divided by the time after the release that the maximum concentration is observed at the sampler. σ_x can be expressed as σ_x = σ_xo + 0.14x, with the initial σ_xo of 45 m. This initial σ_xo agrees with the suggestion of an initial plume spread of about 40 m, made by McElroy and Pooler from analysis of the 1960s’ St. Louis urban dispersion experiment. The puff speeds, u, are initially only about 20% of the observed wind speed, averaged over about 80 street-level and rooftop anemometers in the city, but approach the mean observed wind speed as the puffs grow vertically. The scatter in the σ_t data is about ± a factor of two or three at any given travel time. The maximum σ_t is about 250 s, and the maximum duration of the puff over the sampler, Dt, sometimes called the retention time, is about 1,100 s or 18 min for these puffs and distances.     

Study on horizontal relative diffusion in the troposphere and lower stratosphere

The behaviour of relative diffusion theory and Gifford’s random-force theory for long-range atmospheric diffusion is examined. When a puff scale is smaller than the Lagrangian length scale, 2KTL, an accelerative relative diffusion region exists, i.e., σy∝t3/2. While the puff diffusion enters a two-dimensional turbulence region, in which the diffusion scale is larger than 500 km, or time scale is larger than 1 day, divergence and convergence are main cause of horizontal diffusion. Between the two above-mentioned regimes, diffusion deviation is given byσy=2KTL. The large-scale horizontal relative diffusion parameters were obtained by analyzing the data of radioactive cloud width collected in air nuclear tests.     

Large-Eddy Simulation of Flow and Pollutant Transports in and Above Two-Dimensional Idealized Street Canyons

A large-eddy simulation (LES) model, using the one-equation subgrid-scale (SGS) parametrization, was developed to study the flow and pollutant transport in and above urban street canyons. Three identical two-dimensional (2D) street canyons of unity aspect ratio, each consisting of a ground-level area source of constant pollutant concentration, are evenly aligned in a cross-flow in the streamwise direction x. The flow falls into the skimming flow regime. A larger computational domain is adopted to accurately resolve the turbulence above roof level and its influence on the flow characteristics in the street canyons. The LES calculated statistics of wind and pollutant transports agree well with other field, laboratory and modelling results available in the literature. The maximum wind velocity standard deviations σ _i in the streamwise (σ _u ), spanwise (σ _v ) and vertical (σ _w ) directions are located near the roof-level windward corners. Moreover, a second σ _w peak is found at z ≈ 1.5h (h is the building height) over the street canyons. Normalizing σ _i by the local friction velocity u _*, it is found that σ _u /u _* ≈ 1.8, σ _v /u _* ≈ 1.3 and σ _w /u _* ≈ 1.25 exhibiting rather uniform values in the urban roughness sublayer. Quadrant analysis of the vertical momentum flux u′′w′′ shows that, while the inward and outward interactions are small, the sweeps and ejections dominate the momentum transport over the street canyons. In the x direction, the two-point correlations of velocity R _v,x and R _w,x drop to zero at a separation larger than h but R _u,x (= 0.2) persists even at a separation of half the domain size. Partitioning the convective transfer coefficient Ω _T of pollutant into its removal and re-entry components, an increasing pollutant re-entrainment from 26.3 to 43.3% in the x direction is revealed, suggesting the impact of background pollutant on the air quality in street canyons.     

Spectra of temperature derivatives in the atmospheric surface layer

Measured spectra ofθ _x,θ _y,θ _z, the derivatives of temperature in streamwise, lateral and vertical directions, respectively, indicate that the spectral densities ofθ _z andθ _y are nearly similar but significantly different from the spectral density ofθ _x. The high-frequency parts of the three spectra satisfy, in a qualitative sense, local isotropy requirements. In the high-frequency end of the inertial subrange, the relative behaviour of spectra ofθ _x,θ _y andθ _z is also consistent with local isotropy.     

2003年淮河流域大水期间体积降水量的研究

针对 2 0 0 3年梅雨期淮河流域大水提出了体积降水量的概念及其计算方法 ,在计算出淮河流域和其各子流域逐日和总体积降水量的基础上 ,将体积降水量与水文站的水位和流量进行了对比分析。结果表明 :水位对累计体积降水量有较好的即时响应 ,流量对体积降水量有较好的延迟响应。另外 ,讨论了流域体积降水量的预报问题 ,用国内外数值天气预报产品和中央气象台指导预报产品 ,进行了流域体积降水量预报试验 ,并对预报结果进行了检验分析。结果表明 :在目前天气预报水平条件下 ,利用数值天气预报和中央气象台业务预报产品制作体积降水量是可行的 ,能够延长洪水预报的预见期 ;而且 ,数值天气预报产品在预报体积降水量方面有明显的优势 ,因此可以直接利用数值天气预报产品进行体积降水量的客观预报 ,为防汛工作提供重要的依据。     

Variance Method to Determine Turbulent Fluxes of Momentum And Sensible Heat in The Stable Atmospheric Surface Layer

Evidence is presented that in the stable atmospheric surface layer turbulent fluxes of heat and momentum can be determined from the standard deviations of longitudinal wind velocity and temperature, σ _u and σ _T respectively, measured at a single level. An attractive aspect of this method is that it yields fluxes from measurements that can be obtained with two-dimensional sonic anemometers. These instruments are increasingly being used at official weather stations, where they replace the standard cup anemometer–wind vane system. With methods such as the one described in this note, a widespread, good quality, flux network can be established, which would greatly benefit the modelling community. It is shown that a ‘variance’ dimensionless height (ζ _σ) defined from σ _u and σ _T is highly related to the ‘conventional’ dimensionless stability parameter ζ=z/L, where z is height and L is the Obukhov length. Empirical functions for ζ _σ are proposed that allow direct calculation of heat and momentum fluxes from σ _u and σ _T. The method performs fairly well also during a night of intermittent turbulence.     

Markov Chain Simulations of Vertical Dispersion in the Neutral Surface Layer for Surface and Elevated Releases

Vertical dispersion in the neutral surface layer is investigated using a Markov Chain simulation procedure. The conceptual basis of the procedure is discussed and computation procedures outlined. Wind and turbulence parameterizations appropriate to the neutral surface layer are considered with emphasis on the Lagrangian time scale. Computations for a surface release are compared with field data. Good agreement is found for the variation of surface concentration and cloud height to distances 500 m downwind of the source. The functional form of the vertical concentration profile is examined and an exponential with exponent ∼1.6 is found to give the best fit with simulations. For elevated releases, it is demonstrated that an initial dip of the mass mean height from the simulation can be normalized for various release heights using a non-dimensionalized downwind coordinate incorporating advective wind speed and wind shear. The vertical distribution standard deviation (σ_z), as employed in Gaussian models, shows a fair degree of independence with source height but close examination reveals an optimum source height for maximum σ_z at a given downwind distance,x. This source height increases with downwind distance. Also the simulations indicate that vertical wind shear is more important than vertical variation of Lagrangian time scale close to the source, with a reverse effect farther downwind.     

A Scheme for Scalar Exchange in the Urban Boundary Layer

We present a scheme for parameterising scalar transfer in the urban boundary layer, which is divided into an inertial layer and a roughness layer. The latter is further divided into a shear layer and a canyon layer. In the inertial layer, scalar transfer is determined by turbulence related to canyon macroscopic features, while in the roughness layer, it is determined by shear-generated turbulence, canyon vortex and vortex-generated turbulence. We first describe a conceptual model for the canyon flow and the aerodynamic resistance network, and then estimate the resistances from the point of view of drag partition and vortex advection. The results are compared with the measurements from wind-tunnel experiments. It is found that for small canyon aspect ratio, σ_c, the transfer velocity increases with σ_c, reaching a maximum at around σ_c=0.5 and then decreases with σ_c. We also show that the scheme is not sensitive to adjustable parameters     

Relating Urban Surface-layer Structure to Upwind Terrain for the Salford Experiment (Salfex)

Profiles of wind and turbulence over an urban area evolve with fetch in response to surface characteristics. Sodar measurements, taken on 22 April 2002 during the Salford Experiment in the UK (Salfex), are here related to upstream terrain. A logarithmic layer up to z = 65m was observed in all half-hour averaged profiles. Above this height the profile showed a different vertical gradient, suggesting a change in surface cover upstream. The drag coefficient varied by a factor of two over only a 20° direction change. Turbulence intensity (σ _x ) for each wind component (x) decreased with height, but the ratio suggested an underestimate of σ _u compared to previous results. Mean urban and suburban cover fraction within the source area for each height decreased sharply between z = 20 and 50m, increasing slightly above. The near-convergence of cover fractions thus occured for source areas of minimum length ≈ 2,200 m. In comparison, the mean length scale of heterogeneity L _P was calculated from surface cover data to be 1,284 m, and the corresponding mean blending height h _b was 175 m. Finally, the mean streamline angle, α, was negative and the magnitude decreased with height. An exponential fit to α for z ≤ 65m gave an e-folding height scale of 159 m. A simple relationship between this height scale and L _P was assumed, giving L _P ≈ 1,080 m, which is in reasonable agreement with the estimate from surface cover type. The results suggest that more emphasis is required on modelling and measuring surface-layer flow over heterogeneous urban canopies.     

Sea-breeze Front Variations in Space and Time

Summary  This paper is a contribution to experimental meteorology: A sea-breeze front was investigated by aircraft observations and thorough numerical analysis using an unprecedented number of runs crossing the same front within a timespan of . The 33 runs were flown in a situation of offshore geostrophic wind of 5 m/s in 1000 hPa and with the strategy of obtaining information on the four-dimensional field (t=time, x=cross-coastal coordinate, y=coast-parallel coordinate, z=height): 9 runs in x-direction (and reverse) at different heights to yield x,z-cross-sections of the observed meteorological quantities (specific humidity q, potential temperature Θ and the components u, v and w of the wind velocity), assuming a frozen structure in time; the next 7 runs again in x-direction but all at the same level and on the same track to yield x,t-diagrams of the same quantities in order to study the temporal changes compared to those with x and z; the next 10 runs as a zig-zagging flight track crossing the front but drifting in y-direction, all at the same height, in order to obtain the y-dependency; andfinally 7 runs for another x,z-cross-sectional analysis, which can be compared to that evaluated from the runs at the beginning of the mission. The paper describes the 4-dimensional dependencies in detail. Pure x-variations at constant z are expressed by VCM low-pass filtered space series (VCM=variance conserving multiresolution, according to Howell and Mahrt, 1994). The x,z-analyses are similar to those in Kraus et al. (1990) and Finkele et al. (1995) verifying these results. The comparison of the x,z-studies gained from the data at the beginning and at the end of the mission show how the sea-breeze frontal area changes its structure. The fluctuations (in time) revealed by the low-pass filtered x,t-runs (same track and same height) are smaller than the contour intervals chosen in the x,z-cross-sections. This shows, that the single runs, from which the x,z-cross-sections are constructed, reliably and significantly contribute to the interpolated structure. The paper also demonstrates the overall development of the front within the 3^1/2 h of continuous observation. The x,y-fields demonstrate that the y-dependency of the various quantities is generally one order of magnitude smaller than the x-dependency and that the assumption of negligible y-dependency holds in the first order of approximation for a fairly homogeneous coast. Convective disturbances of a horizontal scale of 1 to 4 km at the landward side of the front, embedded in the offshore flow and bouncing against the landward propagating sea-breeze front, considerably contribute to variations of the frontal propagation speed and of the frontal shape and also to changes of the parameters with the along-frontal coordinate y. Received April 24, 1998 Revised November 3, 1998