Radiative Processes in the Stable Boundary Layer: Part II. The Development of the Nocturnal Boundary Layer

The interaction between radiation and turbulence in the stable boundary layer over land is explored using an idealized model, with a focus on the surface layer after the evening transition. It is shown that finer vertical resolution is required in transitional boundary layers than in developed ones. In very light winds radiative cooling determines the temperature profile, even if similarity functions without a critical Richardson number are used; standard surface similarity theory applied over thick layers then yields poor forecasts of near-surface air temperatures. These points are illustrated with field data. Simulations of the developing nocturnal boundary layer are used to explore the wider role of radiation. Comparatively, radiation is less significant within the developed stable boundary layer than during the transition; although, as previous studies have found, it remains important towards the top of the stable layer and in the residual layer. Near the ground, reducing the surface emissivity below one is found to yield modest relative radiative warming rather than intense cooling, which reduces the potential importance of radiation in the developed surface layer. The profile of the radiative heating rate may be strongly dependent on other processes, leading to quite varied behaviour.     

Stable Atmospheric Boundary Layer Over a Uniform Slope: Some Theoretical Concepts

Theoretical arguments are developed to describe the effects of a uniform slope on the development of the stably stratified atmospheric boundary layer (SBL). A maximum sustainable surface buoyancy flux exists for the SBL overlying a uniform, non-sloping surface. In this study it is shown that the SBL overlying a uniform shallow slope (with gradient of the order of 1:1000) also supports a maximum sustainable buoyancy flux, B _max, but that the value of B _max is influenced by the gradient of the slope, γ. It is demonstrated that in the limit γ → 0, results for the SBL over a horizontal surface are recovered.     

Characteristics of Submeso Winds in the Stable Boundary Layer

The characteristics of submeso motions in the stable boundary layer are examined using observations from networks of sonic anemometers with network sizes ranging from a few hundred metres to 100 km. This study examines variations on time scales between 1 min and 1 h. The analysis focuses on the behaviour of the spectra of the horizontal kinetic energy, the ratios of the three velocity variances, their kurtosis, the dependence of horizontal variability on time scale, and the inter-relationship between vertical vorticity, horizontal divergence and deformation. Motions on larger time and space scales in the stable boundary layer are found to be nearly two-dimensional horizontal modes although the ratio of the vorticity to the divergence is generally on the order of one and independent of scale. One exception is a small network where stronger horizontal divergence is forced by a decrease in surface roughness. The horizontal variability, averaged over 1 h, appears to be strongly influenced by surface heterogeneity and increases with wind speed. In contrast, the time dependence of the horizontal structure on time scales less than one hour tends to be independent of wind speed for the present datasets. The spectra of the horizontal kinetic energy and the ratio of the crosswind velocity variance to the along-wind variance vary substantially between networks. This study was unable to isolate the cause of such differences. As a result, the basic behaviour of the submeso motions in the stable boundary layer cannot be generalized into a universal theory, at least not from existing data.     

Lagrangian Stochastic Modeling of Dispersion in the Stable Boundary Layer

Lagrangian stochastic models are well-suited for modeling dispersion in the stable boundary layer, especially in complex terrain. This note briefly describes the formulations and application of a Lagrangian stochastic model to predict dispersion of tracers released within nocturnal drainage flows.     

Observations of Solitary Waves in the Stable Atmospheric Boundary Layer

Large amplitude, isolated, wave-like phenomena have been observed in the lowest 40 m of the strongly stably stratified atmospheric boundary layer overlying a coastal Antarctic ice shelf. The waves only occur when prevailing wind speeds are low. They always propagate from over the land, with phase speeds exceeding the local mean wind speed. They have wavelengths of the order of 200 m. Several examples are described and a summary of the statistical properties of these waves events is presented.     

Turbulence Characteristics of the Stable Boundary Layer Over a Mid-Latitude Glacier. Part I: A Combination of Katabatic and Large-Scale Forcing

Profile and eddy-correlation (heights of 4 and 10 m) measurements performed on the Pasterze glacier (Austria) are used to study the characteristics of the stable boundary layer under conditions of katabatic and large-scale forcing. We consider cases where large-scale forcing results in a downslope (or following) ambient wind. The analysis of averaged spectra and cospectra reveals low frequency perturbations that have a large influence on the variances of temperature and horizontal wind components and also alter the cospectra of momentum and sensible heat flux. Only the spectrum of the vertical wind speed is comparable to universal spectra. The low frequency perturbations occur as brief intermittent events and result in downward entrainment of ambient air thereby producing enhanced downward sensible heat fluxes and downward as well as upward momentum fluxes with various magnitudes and timescales. After the variances were high pass filtered, the normalised standard deviations of wind speed and temperature compare favourably to findings in the literature within the range 0>z/L>0.5. For larger z/L they deviate as a result of an increased influence from low frequency perturbations and thus non-stationarity. In line with this, the turbulent kinetic energy budget (at 4 m height) indicates that production (shear) is in balance with destruction (buoyancy and dissipation) within the range 0>z/L>0.3. Non-dimensional gradients of wind speed within the range 0>z/L>0.3 have a slope of about 3.5. The scatter for the dimensionless temperature gradient is quite large, and the slope is comparable to that for wind speed gradients. For z/L>0.3 the imbalance in the turbulent kinetic energy budget grows and non-dimensional gradients for wind speed and temperature deviate considerably from accepted values as a result of increased non-stationarity. Average roughness lengths for momentum and sensible heat flux derived from wind speed and temperature profiles are respectively 1 × 10^-3 m and 6 × 10^-5 m, consistent with the literature. The ratio (z_0h/z_0m) compares to those predicted by surface renewal models. A variation of this ratio with the roughness Reynolds number is not indicated by our data.     

Local Scales of Turbulence in the Stable Boundary Layer

Local, gradient-based scales, which contain the vertical velocity and temperature variances, as well as the potential temperature gradient, but do not include fluxes, are tested using data collected during the CASES-99 experiment. The observations show that the scaling based on the temperature variance produces relatively smaller scatter of empirical points. The resulting dimensionless statistical moments approach constant values for sufficiently large values of the Richardson number Ri. This allows one to derive predictions for the Monin–Obukhov similarity functions φ _m and φ _h , the Prandtl number Pr and the flux Richardson number Rf in weak turbulence regime.     

Local Similarity in the Stable Boundary Layer and Mixing-Length Approaches: Consistency of Concepts

In stably stratified flows vertical movement of eddies is limited by the fact that kinetic energy is converted into potential energy, leading to a buoyancy displacement scale z _B . Our new mixing-length concept for turbulent transport in the stable boundary layer follows a rigid-wall analogy, in the sense that we assume that the buoyancy length scale is similar to neutral length scaling. This implies that the buoyancy length scale is: ℓ _B  = κ _B z _B , with κ _B ≈ κ, the von Karman constant. With this concept it is shown that the physical relevance of the local scaling parameter z/Λ naturally appears, and that the α coefficient of the log-linear similarity functions is equal to c/κ ², where c is a constant close to unity. The predicted value α ≈ 1/κ ² = 6.25 lies within the range found in observational studies. Finally, it is shown that the traditionally used inverse linear interpolation between the mixing length in the neutral and buoyancy limits is inconsistent with the classical log-linear stability functions. As an alternative, a log-linear consistent interpolation method is proposed.     

Turbulence Characteristics Of The Stable Boundary Layer Over A Mid-Latitude Glacier. Part Ii: Pure Katabatic Forcing Conditions

Observations obtained over a glacier surface in a predominantlykatabatic flow and with a distinctwind maximum below 13-m height are presented. The data werecollected using a 13-m high profilemast and two sonic anemometers (at about 2.5-m and 10-m heights).The spectra at frequencies belowthat of the turbulence range appear to deviate considerably fromthe curves obtained by Kaimal andco-workers during the 1968 Kansas experiment. The characteristicsof these deviations are compared tothe observations of others in surface-layers disturbed by anykind of large-scale outer-layer (orinactive) turbulence. In our case the disturbances arelikely to be induced by the highmountain ridges that surround the glacier. Moreover, the deviationsobserved in the cospectra seemto result from an, as yet, unspecified interaction between theinactive outer-layer turbulenceand the local surface-layer turbulence. Near the distinctwind maximum turbulence production ceasedwhile turbulence itself did not, probably the result ofturbulence transport from other levels. Consequently, we studied thelocal similarity relations using _w instead of u_* as an alternative velocity scale. Wellbelow the wind maximum, and for relatively low stability(0< Ri_g <0.2), the flow behaves accordingto well established local-scaling similarity relationshipsin the stable boundary layer. For higherstability (Ri_g > 0.2), and near or above the wind maximum, the boundary-layer structure conforms tothat of z-less stratification suggesting that the eddy sizeis restricted by the local stability ofthe flow. In line with this we observed that the sensibleheat fluxes relate remarkably well to thelocal flow parameters.     

A Note on Wave-Turbulence Interaction and the Possibility of Scaling the Very Stable Boundary Layer

Success in scaling the very stable boundary layer has been mixed. One possible reason for this is that wave-turbulence interaction can result in significant energy coupling between turbulence within the planetary boundary layer and the mean flow outside the PBL. Three regimes are described where wave-turbulence interaction forms a significant part of the dynamics. In the first, energy and momentum flows are confined within the PBL and do not preclude local scaling. In the other two, which involve topographically forced waves and propagating waves respectively, energy/momentum transfers across the PBL make success in local scaling unlikely.     

Intermittent Turbulence Associated with a Density Current Passage in the Stable Boundary Layer

Using the unprecedented observational capabilities deployed duringthe Cooperative Atmosphere-Surface Exchange Study-99 (CASES-99),we found three distinct turbulence events on the night of 18October 1999, each of which was associated with differentphenomena: a density current, solitary waves, and downwardpropagating waves from a low-level jet. In this study, we focus onthe first event, the density current and its associatedintermittent turbulence. As the cold density current propagatedthrough the CASES-99 site, eddy motions in the upper part of thedensity current led to periodic overturning of the stratifiedflow, local thermal instability and a downward diffusion ofturbulent mixing. Propagation of the density current induced asecondary circulation. The descending motion following the head ofthe density current resulted in strong stratification, a sharpreduction in the turbulence, and a sudden increase in the windspeed. As the wind surge propagated toward the surface, shearinstability generated upward diffusion of turbulent mixing. Wedemonstrate in detail that the height and sequence of the localthermal and shear instabilities associated with the dynamics ofthe density current are responsible for the apparent intermittentturbulence.     

宜昌冬季两次降水过程大气边界层的观测分析

利用宜昌2007年12月10-25日的加密观测资料,分析了两次低值系统经过宜昌时大气边界层的温湿风廓线结构及其日变化特征。结果表明：位温廓线具有明显的日变化特征,对流边界层在白天出现和发展,其高度可达600m,而稳定边界层在夜间出现和发展,其高度可达300m,降水会抑制对流边界层和稳定边界层的发展;湿度廓线结构及其日变化与对流边界层的发展有关,总体上湿度随高度减小,贴近地面的薄层湿度随高度减小较快,而混合层内湿度随高度变化较小,出现降水时,近地层的湿度有明显增加,大气边界层内湿度随高度快速平稳减小;风速廓线结构比较复杂,总体上风速随高度增大,在大气边界层低层有时会出现一个风速极大值,风速廓线没有明显的日变化特征,大气边界层内风向变化较大,但以偏东风为主。     

Large-Eddy Simulation of the Stable Atmospheric Boundary Layer using Dynamic Models with Different Averaging Schemes

Large-eddy simulation (LES) of a stable atmospheric boundary layer is performed using recently developed dynamic subgrid-scale (SGS) models. These models not only calculate the Smagorinsky coefficient and SGS Prandtl number dynamically based on the smallest resolved motions in the flow, they also allow for scale dependence of those coefficients. This dynamic calculation requires statistical averaging for numerical stability. Here, we evaluate three commonly used averaging schemes in stable atmospheric boundary-layer simulations: averaging over horizontal planes, over adjacent grid points, and following fluid particle trajectories. Particular attention is focused on assessing the effect of the different averaging methods on resolved flow statistics and SGS model coefficients. Our results indicate that averaging schemes that allow the coefficients to fluctuate locally give results that are in better agreement with boundary-layer similarity theory and previous LES studies. Even among models that are local, the averaging method is found to affect model coefficient probability density function distributions and turbulent spectra of the resolved velocity and temperature fields. Overall, averaging along fluid pathlines is found to produce the best combination of self consistent model coefficients, first- and second-order flow statistics and insensitivity to grid resolution.     

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.     

Airborne Measurements in the Stable Boundary Layer over the Larsen Ice Shelf,Antarctica

We present aircraft measurements of boundary-layer structure and surface turbulent fluxes from a flight over the Larsen Ice Shelf, Antarctica. Warm advection, associated with föhn flow, led to the formation of a stable boundary layer over the ice shelf, with a well-defined low-level jet at the top of the surface inversion. The strong shear associated with the jet kept the gradient Richardson number small and maintained a turbulent boundary layer over a depth of at least 600 m. The net surface energy balance amounted to 52 Wm^?2, equivalent to a melt rate of 13 mm water per day, with net radiation (48 Wm^?2) making the largest contribution to melt. The contribution from the sensible heat flux (13 Wm^?2) was largely balanced by an upwards latent heat flux (?9 Wm^?2). These measurements provide insight into the processes that control surface melt rates in an area that has experienced recent rapid warming and deglaciation.     

Surface-Layer Characteristics in the Stable Boundary Layer with Strong and Weak Winds

An extensive meteorological dataset obtained from the plumevalidation experiment conducted by the Electric Power Research Institute (EPRI) atKincaid during 1980–1981 is analysed for studying the characteristic differences in thesurface-layer parameters in strong and weak wind stable conditions. The surface-layerparameters are computed using the similarity functions _m and _h proposed byBeljaars and Holtslag. The weak winds are characterized using the geostrophic wind speedas well as the wind speed at the 10-m level. The surface fluxes are found to be finitein weak wind conditions.Empirical formulations for the eddy diffusivities of momentum(K_M) and heat (K_H), and drag (C_D) and heat exchange (C_H) coefficients, as powerlaw functions of the bulk Richardson number (Ri_B), are proposed under both strong andweak wind conditions. Results are close to those based on observations taken from the IndianInstitute of Technology low wind diffusion experiment, the Land surface processes experiment,the Hanford diffusion experiment, the Cabauw field experiment and the Cooperative Atmospheric SurfaceExchange Study 1999 (CASES-99) experiment. In addition, the fluxes obtained fromthe proposed empirical relations are in good agreement with those based on similarity theory as wellas the turbulence measurements taken from the CASES-99 experiment.     

The Stable Atmospheric Boundary Layer over an Antarctic ice Sheet

Turbulence measurements up to 11-m height and longterm profile measurements up to 45-m height performed at the German Neumayer Station in Antarctica are used to investigate different components of turbulence closure schemes of the stable boundary layer. The results confirm the linear relationships for the universal functions of momentum and heat exchange in the stability range z/L < 0.8 ... 1, whereas the local scaling approach should be used above the surface layer. Furthermore, boundary-layer heights below 50 m are frequently observed at this coastal Antarctic site, mainly due to the influence of stability above the boundary layer. It is shown that the inclusion of this stability into parametrization relations is necessary to provide realistic equilibrium heights of the stable boundary layer. Two relations, based on different physical approaches, were successfully applied for the parametrization of the equilibrium height.     

自由对流与稳定层结边界层风廓线的解析表达和边界层顶抽吸速度

基于近年来对自由对流和稳定边界层湍流交换特征的研究,求解边界层运动方程,得到这两种层结下边界层风的解析表达式.所得廓线与边界层特性参数符合观测特征.还求出了这二种层结下边界层顶抽吸速度的解析表达及其与某些参数的关系.结果表明,抽吸速度与层结有关,其特征可从物理上加以解释.     

Resolution Sensitivity and Scaling of Large-Eddy Simulations of the Stable Boundary Layer

Large-eddy simulations (LES) of the continuously turbulent quasi-equilibrium stable boundary layer (SBL) are conducted with grid lengths in the range of 12.5 m to 2 m, in order to explore resolution sensitivity, and determine at what point grid convergence occurs. The structure of the mean potential temperature, winds, and turbulent fluxes varies significantly over this resolution range. The highest resolution simulations show a significant degree of convergence. The dimensionless momentum diffusivity asymptotes to a value of 0.06, corresponding to a limiting flux Richardson number of 0.15.Using the converged simulations, some scaling hypotheses underpinning first-order and second-order closure models are revisited. The effective Richardson number stability functions of the LES are compared with the forms often used in numerical weather prediction (NWP). The mixing implied by the LES is less than that used in NWP. The commonly used similarity profiles for heat and momentum fluxes, and the scalings for dissipation and pressure covariances are compared with the LES. This information could provide guidance for the next generation of SBL parametrization schemes.