Richardson number profiles through shear instability wave regions observed in the lower planetary boundary layer

Remote sensing of the lower planetary boundary layer in the vicinity of a meteorological tower on many occasions reveals the existence of shear instability (Kelvin-Helmholtz) waves. In general, such waves are found within shallow strata which are marked by strong thermal stability and large vertical wind shear. The independent and concurrent measurements of the vector wind and temperature, made on a 152-m high tower, allow the construction of wind and temperature profiles. From such measurements, the Richardson number profile is constructed as well as the instability regime according to Drazin's criterion. The results show that regions of shear-instability waves as depicted by the remote sensor (an acoustic sounder) agree well with Drazin's instability regime, and that within such regions the Richardson number is indeed 0.25.     

A bulk model for the atmospheric planetary boundary layer

The integrated momentum and thermodynamic equations through the planetary boundary layer (PBL) are solved numerically to predict the mean changes of wind and potential temperature from which surface fluxes are computed using bulk transfer coefficients of momentum and heat. The second part of the study involves a formulation and testing of a PBL height model based on the turbulent energy budget equation where turbulent fluxes of wind and heat are considered as the source of energy. The model exhibits capability of predicting the PBL height development for both stable and unstable regimes of observed conditions. Results of the model agree favourably with those of Deardorff's (1974a) and Tennekes' (1973) models in convective conditions.Contribution number 396.     

A study of gravity waves in the planetary boundary layer by acoustic sounding

Gravity waves associated with stably stratified layers were observed in the planetary boundary layer at two locations in France. Using an array of three monostatic Doppler acoustic sounders, the wave speeds and directions were found. A quantitative study of the waves has been performed to determine their physical characteristics. Based on their dynamics, a classification into three types was possible. Most were found to be trapped beneath a critical level.     

On the instability of a planetary boundary layer with Rossby-number similarity

The formation of longitudinal vortex rolls in the planetary boundary layer (PBL) is investigated by means of perturbation analysis. The method is the same as that used by previous authors who have investigated the instability of a laminar Ekman layer. To study the instability of the turbulent boundary layer of the atmosphere, vertical profiles are needed of the eddy viscosity and of the two components of the basic flow. These profiles have been obtained by a numerical PBL-model; they are universal for zz ₀. (However, the stability equations are not completely universal, i.e., independent of the external parameters). For each thermal stratification, expressed by the internal stratification parameter , one has a set of three consistent profiles.The numerical solution of the stability equations gives the critical values and the perturbation growth rates as functions of thermal stratification and of the surface Rossby number Ro₀. This is in contrast to the case of a laminar Ekman layer, where the instability depends on a Reynolds number only, which makes atmospheric applications difficult. The most unstable perturbations are found for Ro₀ = 1 × 10⁶ approximately, which is in the range of surface Rossby numbers observed in the atmosphere. However, considering vortex rolls oriented in the direction of the surface stress, the instability is found to be universal, i.e., independent of the external parameters combined in the surface Rossby number. With respect to thermal stratification, the results show that the instability of the perturbations increases with increasing static stability.     

On the instability of a planetary boundary layer with Rossby number similarity

In Part I (1975), a linear stability analysis with respect to the formation of longitudinal vortex rolls was given for aturbulent boundary layer of the atmosphere. However, that analysis investigated the effect of inflection point instability only; therefore it is applicable only to the case of neutral stratification. In Part II presented here, the analysis is extended to include the combined effect of inflection point instability and instability due to heating from below. In contrast with the result in Part I, the main result is that in considering both these effects, longitudinal vortex rolls can develop only if the boundary layer has an unstable stratification. Another important result is that the structure of developing vortex rolls and their growth-rates are universal in a boundary layer with Rossby-number similarity, i.e., they are independent of any external parameter. The same is true for the orientation of the vortex rolls: the angle between the axis of the rolls and the surface stress is independent of external parameters. The only quantity which is not universal is the phase speed, which indicates the speed with which the rolls move in a direction perpendicular to the vortex axis; this phase speed depends on the geostrophic wind and on the roughness-length. Paper presented at the XIIIth Biennual Fluid Dynamics Symposium 5–10 September, 1977, Olsztyn, Poland.     

Structure of mean winds and turbulence in the planetary boundary layer over rural terrain

A detailed analysis has been carried out of the temporal and spatial structure of mean winds and turbulence in the neutrally-stable planetary boundary layer over typically rural terrain. The data were obtained from a horizontal array of tower-mounted propeller anemometers (z = 11 m) during a five-hour period for which the mean wind direction was virtually perpendicular to the main span of the array. Various turbulence characteristics have been obtained for all three components of velocity and have been compared with idealized models for such a flow and with some of the other available atmospheric results.Considerable tower-to-tower and block-to-block variability has been observed in many of the measured results, particularly in those for the horizontal-component integral scales. Surface shear stress, roughness length and turbulence intensities were in good agreement with expected values for such a site. Power spectra for all components displayed significantly more energy at middle and lower frequencies than that observed by Kaimal et al. (1972) over flat, relatively featureless terrain. This is felt to be a result of the generally rougher gross features of the terrain in the present case and has led to the development of a modified version of the Kaimal-spectral model which fits the observed data better than either the original Kaimal model or the von Kármán model. It is suggested that it may in future be possible to represent power spectra over a wide range of terrain types by using such a modified spectral model.Integral scales of turbulence were calculated by three different techniques and in most cases displayed a strong dependence on the technique used. Averaged values of scale showed reasonable agreement with most of the available atmospheric data and with the values suggested by ESDU (1975). The anticipated elongation of turbulent eddies in the longitudinal direction was confirmed for all three velocity components, although it was found to be not as large as some other observations.     

Effects of turbulent dispersion of atmospheric balance motions of planetary boundary layer

New Reynolds' mean momentum equations including both turbulent viscosity and dispersion are used to analyze atmospheric balance motions of the planetary boundary layer. It is pointed out that turbulent dispersion with r 0 will increase depth of Ekman layer, reduce wind velocity in Ekman layer and produce a more satisfactory Ekman spiral lines fit the observed wind hodograph. The wind profile in the surface layer including tur-bulent dispersion is still logarithmic but the von Karman constant k is replaced by k1 = 1 -2/k, the wind increasesa little more rapidly with height.     

The role of radiation in an inversion-capped planetary boundary layer

In an inversion-capped planetary boundary layer (PBL), the structure of the turbulent fluxes as well as the height of the inversion are determined by the interaction of turbulent mixing in the PBL, large-scale subsidence above the PBL and radiational cooling. Here the sensitivity of the inversion height and of the turbulent fluxes due to radiational processes is investigated with the aid of a three-layered model for a well mixed PBL. For an example of the Trade-Wind region, the inversion height (i.e., the difference between surface pressure and pressure at the inversion level) varies between 46 and 257 mb and the surface flux of moist static energy between 417 and 99 W m^-2, if the (mean) radiative net flux divergence for both the inversion and the well-mixed layer is changed over a reasonable range of values. None of the parameterization schemes existing in the literature is able to describe these radiational effects in an appropriate way. This is due to the fact that these parameterizations are either not or not flexibly enough linked to the thermodynamical model parameter. Therefore the demand for an adequate parameterization of the radiational influence in a well-mixed PBL under a subsidence inversion is obvious.     

Scaling turbulence in the planetary boundary layer

Two different approaches to scaling turbulence in the planetary boundary layer over Lake Ontario are investigated. The height up to the inversion was found to be the appropriate scaling height while u. for near‐neutral and w* for unstable conditions were the appropriate scaling velocities. The results were in general agreement with the numerical models of Deardorff (1972) and Wyngaard, Cote, and Rao (1974).     

Turbulence structure in the planetary boundary layer

This review of the last three years of progress in the understanding of wind profiles and the structure of turbulence in the planetary boundary layer is divided into three parts. The first part, by N. E. Busch, deals with the atmospheric surface layer below 30 m. It is shown that the Monin-Oboukhov similarity hypotheses fail at low frequencies and large wave-lengths, probably due to mesoscale influences. Also, it is suggested that the neutral surface layer is a poor reference state in some respects, because the structure of turbulence in unstable conditions is quite different from that in stable stratification. The second part, by H. Tennekes, is concerned with the intermittency of the dissipative structure of turbulence and its effects on the velocity and temperature structure functions. It is shown that the modified Kolmogorov-Oboukhov theory, which attempts to explain the consequences of the dissipative intermittency, is unable to predict the shape of the temperature structure functions. The third part of this review, by H. A. Panofsky, deals with wind profiles and turbulence structure above 30 m. It is shown that between 30 and 150 m, surface-layer formulas can be used, if such mesoscale effects as changes of terrain roughness are taken into account where needed. Experimental data on turbulence above 150 m are quite sparse; some of the current scaling laws that can be used in this region are described.     

大气边界层中Kelvin-Helmholtz波对降水的作用

在夜间稳定边界层大气中,北京受燕山山脉地形的影响,重力流一般在300米以下的高度内传播,波的传播速度基本等于切变层的平均风速.在一定的高度上存在大风轴,并形成低层风的强切变,大风中心一般在100米高度附近,强度约5米/秒左右. 本文讨论了Helmholtz波不稳定性对降水的影响,在Ri≤0.25的临界值时,与重力流相伴随有时可有小雨,当大尺度盛行气流为东风时可出现大雨.     

Three-dimensional numerical study of the height and mean structure of a heated planetary boundary layer

The heated boundary layer for DAY 33 of the Wangara data of southeast Australia (Clarke et al., 1971) is studied numerically with a three-dimensional model using 64000 grid points within a volume 5 km on a side and 2 km deep. Subgrid-scale transport equations were utilized in place of eddy-coefficient formulations. The rate of growth of the mixed layer is examined and parameterized, and the vertical profiles of heat flux, moisture flux and momentum fluxes are examined. The momentum boundary layer is found to coincide essentially with the mixed layer, and to grow with the latter during the hours of solar heating of the surface.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Equilibrium theory of the planetary boundary layer with an inversion lid

The maintenance of an elevated inversion in steady flow above a cold, rotating surface is shown to be possible for a certain range of the buoyancy number bfV _g, where b is the buoyant acceleration appropriate to the density deficiency of the fluid above the inversion, f is Coriolis parameter and V _gis geostrophic velocity (so that fV _gis also horizontal pressure gradient in kinematic units). The height of the inversion lid is determined by a balance of surface stress and buoyancy, in a way which may be deduced from laboratory experiments. With the aid of such empirical evidence a theory is constructed for the layer below the inversion lid. The cross-isobar angle of ground-level stress is found to increase with the buoyancy number, to a limiting value of 90, by which time the inversion descends to the ground. Under typical conditions, a temperature difference of order 10C is necessary to eliminate the possibility of an equilibrium, elevated inversion lid and reduce ground level wind stress to a vanishingly small value.Woods Hole Oceanographic Institution Contribution #3011On leave from the University of Waterloo, Ontario     

Properties of pressure fluctuations in an atmospheric boundary layer

From experimental data it is shown that there is a relationship between measured wind speed and associated micropressure levels. It is also shown that turbulent eddies associated with the wind vanish within distances of five to six times their size, and that these eddies exhibit normal dispersion. These results confirm, at the atmospheric scales, results previously observed only in model experiments.     

The vertical structure of the atmospheric planetary boundary layer in undisturbed trade winds over the Atlantic Ocean

During the Atlantic Expedition 1965 and the Atlantic Tradewind Experiment (ATEX) 1969, shipborne aerological measurements were obtained in order to investigate the thermodynamical and kinematic structure of the planetary boundary layer in low latitudes. Under suppressed convection, the subdivision of this layer into five sub-regions was found to be a rather permanent feature. Enhanced cumulus convection effects a smoothing of the thermodynamical discontinuities and leads sometimes to a destruction of the trade inversion. Due to the surface pressure distribution and the thermal wind distribution in the lower atmosphere, the actual wind speed and direction are nearly constant with height below cloud base. In the cloud layer up to the inversion, the wind speed generally decreases and the air flow tends to become more zonal.     

尺度自适应大气边界层参数化改进及其对一次海雾的数值模拟研究

大气边界层湍流运动是地球大气运动最重要的能量输送过程之一.当数值模式分辨率接近活跃含能湍涡长度尺度时,湍流运动被部分解析,被称为"灰色区域",传统的边界层方案不适合此时模式湍流问题的描述.为了提高模式边界层方案在包括"灰色区域"的不同网格尺度上的描述能力,适应不同分辨率模式的需要,在雷诺平均湍流理论基础上,修正Mell...     

Numerical model of a non-steady atmospheric planetary boundary layer,based on similarity theory

A numerical model of a non-stationary atmospheric planetary boundary layer (PEL) over a horizontally homogeneous flat surface is derived on the basis of similarity theory. The two most typical turbulence regimes are reproduced: one corresponding to a convectively growing PBL and another corresponding to a stably stratified quasi-equilibrium PBL. The PBL is treated as a unit. Hence, its evolution is characterized by temporal changes of intrinsic parameters, such as the PBL depth, the vertically averaged values of potential temperature, specific humidity and components of wind velocity, the near-surface values of heat, water vapor and momentum fluxes. The internal structure of the PBL is considered self-similar. This allows one to represent the interaction between the air flow and the underlying surface by means of universal heat/mass transfer and resistance laws. Numerical experiments on the diurnal variations of meteorological fields in the lower 2 km layer confirm the ability of the model to reproduce the main features of the phenomena, known from observations.On leave from: Insitute of Limnology, Academy of Sciences of the Russia, SU-196199, St. Petersburg, Russia.     

Simulation of three-dimensional wind flow over complex terrain in the atmospheric boundary layer

A three-dimensional model for wind prediction over rough terrain has been developed for practical use. It is a compromise between hydrodynamic and objective wind models. The proposed model includes: (1) a statistical model to predict the wind velocity and potential temperature at anemometer height at observing stations, (2) the drainage wind model expressed by Prandtl's analytic solution for the slope wind, (3) the Businger-Dyer surface-layer formulation which considers the surface energy budget and (4) the model for three-dimensional boundary-layer solutions to the stationary flow. In this model, mass consistency is guaranteed by using flow fields that satisfy the continuity equation. Model predictions show good agreement with the observations.     

Parametrization of orographical effects in the planetary boundary layer

Studies of the influence of orography on the dynamics of atmospheric processes usually assume the following relation as a boundary condition at the surface of the Earth, or at the top of the planetary layer: $$w = u\frac{{\delta z_0 }}{{\delta x}} + v\frac{{\delta z_0 }}{{\delta y}}$$ where u, v and w are the components of wind velocity along the x, y and z axes, respectively, and z ₀ = z₀(x, y) is the equation of the Earth's orography. We see that w, and consequently the influence of orography on the dynamics of atmospheric processes, depend on the wind (u, v) and on the slope of the obstacle (δz ₀/δx, δz₀/δy). In the present work, it is shown that the above relation for w is insufficient to describe the influence of orography on the dynamics of the atmosphere. It is also shown that the relation is a particular case of the expression: $$\begin{gathered} w_h = \left| {v_g } \right|\left[ {a_1 (Ro,s)\frac{{\delta z_0 }}{{\delta x}} + a_2 (Ro,s)\frac{{\delta z_0 }}{{\delta y}}} \right] + \hfill \\ + \frac{{\left| {v_g } \right|^2 }}{f}\left[ {b_1 (Ro,s)\frac{{\delta ^2 z_0 }}{{\delta x^2 }} + b_2 (Ro,s)\frac{{\delta ^2 z_0 }}{{\delta y^2 }} + b_3 (Ro,s)\frac{{\delta ^2 z_0 }}{{\delta x\delta y}}} \right] \hfill \\ \end{gathered} $$ where ¦vv _g¦ is the strength of the geostrophic wind, a ₁, a₂, b₁, b₂, b₃ are functions of Rossby number Ro and of the external stability parameter s. The above relation is obtained with the help of similarity theory, with a parametrization of the planetary boundary layer. Finally, the authors show that a close connection exists between the effects described by the above expression and cyclo- and anticyclogenesis.     

Velocity profiles,the resistance law and the dissipation rate of mean flow kinetic energy in a neutrally and stably stratified planetary boundary layer

The logarithmic + polynomial approximation is suggested for vertical profiles of velocity components in a planetary boundary layer (PBL) at neutral and stable stratification. The resistance law functions A and B are determined on the basis of this approximation, using integral relations derived from the momentum equations, the Monin-Obukhov asymptotic formula for the wind profile in a stably stratified near-surface layer and the known expressions for the PBL depth. This result gives a realistic and convenient method for calculating the surface friction velocity and direction and the total dissipation rate of mean flow kinetic energy in terms of geostrophic velocity, buoyancy flux at the surface, the roughness parameter and the Coriolis parameter. In the course of these derivations a review is given of current views on the main problems of the neutral and stable PBL.