北京城市化发展对大气边界层特性的影响

利用中国科学院大气物理研究所大气边界层物理与大气化学国家重点实验室的北京325 m气象塔1993年～2003年夏季 (7月～9月) 的观测资料, 统计分析了各年的风速与温度廓线分布特征。统计分析结果表明, 随着城市化的发展, 相对风速有逐年减小的趋势, 并且越靠近地面, 相对风速的减小越明显, 这反映了城市建筑对近地面层空气流动的摩擦作用。对风速廓线进行线性拟合得到风速随高度的垂直递增率, 发现无论是100 m以下的近地面层还是较高层, 风速的垂直递增率都随城市化发展存在逐年增大的趋势, 表明粗糙下垫面的影响已经向高层扩展。根据温度廓线计算了各年的温度垂直递减率, 发现其有增大的趋势, 这表明城市化发展对边界层热力结构同样有显著影响。本文还依据统计整理得到的近中性层结下的风速廓线资料, 利用莫宁-奥布霍夫相似理论计算了下垫面的空气动力学参数, 结果表明, 地表粗糙度、 零平面位移随着城市化发展皆有明显增加的趋势。同时, 分析了各空气动力学参数与平均风速及无量纲风速的关系。其中, 摩擦速度和平均风速二者基本成正相关, 且摩擦速度随平均风速的增大而增大的趋势越发明显。本文研究结果对研究城市化发展对区域大气边界层结构、 气候和环境影响有参考意义, 可为城市大气边界层模式和区域气候模式提供参数化依据。     

井周钻孔补砾工艺在修井中的应用

提出了通过井周钻孔恢复旧井填砾层，提高管井涌水量的方法，通过在某电厂供水管井中的实践，证明效果良好。     

深圳两次大雾天气过程对比分析及预报启示

2005年2月23-25日和2006年3月6日深圳分别出现了一次大雾天气。从天气学角度对两次大雾过程的形成原因、特点进行对比分析，分析表明：大雾天气需在一定的形势场中出现并维持，近地面层气象要素场的变化会促进大雾的形成、维持和消失，而近地面层风场或温度场的改变除了与其环境、当地气候特征有关外，与大的形势场是分不开的。通过对比，对深圳大雾天气的生消和维持机制有一定了解，对预报本地大雾天气有指示作用，也为其它地区特别是沿海地区雾的预报提供借鉴。     

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.     

Atmospheric boundary-layer structure observed during a haze event due to forest-fire smoke

During a haze event in Baltimore, U.S.A. from July 6 to 8, 2002, smoke from forest fires in the Québec region (Canada), degraded air quality and impacted upon local climate, decreasing solar radiation and air temperature. The smoke particles in and above the atmospheric boundary layer (ABL) served as a tracer and provided a unique opportunity to investigate the ABL structure, especially entrainment. Elastic backscatter lidar measurements taken during the haze event distinctly reveal the downward sweeps (or wisps) of smoke-laden air from the free atmosphere into the ABL. Visualisations of mechanisms such as dry convection, the entrainment process, detrainment, coherent entrainment structures, and mixing inside the ABL, are presented. Thermals overshooting at the ABL top are shown to create disturbances in the form of gravity waves in the free atmosphere aloft, as evidenced by a corresponding ripple structure at the bottom of the smoke layer. Lidar data, aerosol ground-based measurements and supporting meteorological data are used to link free atmosphere, mixed-layer and ground-level aerosols. During the peak period of the haze event (July 7, 2002), the correlation between time series of elastic backscatter lidar data within the mixed layer and the scattering coefficient from a nephelometer at ground level was found to be high (R=0.96 for z =324 m, and R=0.89 for z=504 m). Ground-level aerosol concentration was at a maximum about 2 h after the smoke layer intersected with the growing ABL, confirming that the wisps do not initially reach the ground.     

An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer

Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.     

A Large-Eddy Simulation and Lagrangian Stochastic Study of Heavy Particle Dispersion in the Convective Boundary Layer

Large-eddy simulation and Lagrangian stochastic dispersion models were used to study heavy particle dispersion in the convective boundary layer (CBL). The effects of various geostrophic winds, particle diameters, and subgrid-scale (SGS) turbulence were investigated. Results showed an obvious depression in the vertical dispersion of heavy particles in the CBL and major vertical stratification in the distribution of particle concentrations, relative to the passive dispersion. Stronger geostrophic winds tended to increase the dispersion of heavy particles in the lower CBL. The SGS turbulence, particularly near the surface, markedly influenced the dispersion of heavy particles in the CBL. For reference, simulations using passive particles were also conducted; these simulation results agreed well with results from previous convective tank experiments and numerical simulations.     

Maximal Overlap Wavelet Statistical Analysis With Application to Atmospheric Turbulence

Statistical tools based on the maximal overlap discrete wavelet transform (MODWT) are reviewed, and then applied to a dataset of aircraft observations of the atmospheric boundary layer from the tropical eastern Pacific, which includes quasi-stationary and non-stationary segments. The wavelet methods provide decompositions of variances and covariances, e.g. fluxes, between time scales that effectively describe a broadband process like atmospheric turbulence. Easily understood statistical confidence bounds are discussed and applied to these scale decompositions, and results are compared to Fourier methods for quasi-stationary turbulence. The least asymmetric LA(8) wavelet filter yields coefficients that exhibit better uncorrelatedness across scales than the Haar filter and is better suited for decomposition of broadband turbulent signals. An application to a non-stationary segment of our dataset, namely vertical profiles of the turbulent dissipation rate, highlights the flexibility of wavelet methods.     

Turbulent Kinetic Energy Dissipation in the Surface Layer

We estimated the turbulent kinetic energy (TKE) dissipation rate for thirty-two 1-h intervals of unstable stratification covering the stability range 0.12 ≤ −z/L ≤ 43 (z/L is the ratio of instrument height to the Obukhov length), by fitting Kolmogorov’s inertial subrange spectrum to streamwise spectra observed over a desert flat. Estimated values are compatible with the existence of local equilibrium, in that the TKE dissipation rate approximately equalled the sum of shear and buoyant production rates. Only in the neutral limit was the turbulent transport term in the TKE budget measured to be small.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).