The Role of Shear in the Morning Transition Boundary Layer

We use large-eddy simulation (LES) to better define the early stages of the morning transition boundary layer. Previous LES studies relating to the morning transition boundary layer focus on the role of the entraining convective boundary layer (CBL). By using a combination of different domain sizes and grid lengths, the full evolution from the stable boundary layer (SBL) to the CBL is modelled here. In the early stages of the morning transition the boundary layer is shown to be a combination of a shallow mixed layer capped by a significant shear driven stable boundary layer (the so-called mixed CBL–SBL state). The mixed CBL–SBL state is the key to understanding the sensitivity to shear. Turbulent kinetic energy budgets also indicate that it is shear driven. The negative flux from the mixed CBL–SBL state extends much further above the minimum than is typically found for the CBL later in the day, and the depth of penetration scales as w _m /N _i , where w _m is the combined friction and convective velocity scale and N _i the static stability at the inversion top.     

Observations Of The Morning Transition Of The Convective Boundary Layer

The morning transition between the stable nocturnal situation and the daytime convective boundary layer (CBL) is of interest both for basic understanding and for initializing prognostic models. While the morning growth phase of the CBL has been studied in detail, relatively little has been published on the transition itself. In this paper, conventional observations of surface temperature, humidity, and turbulent fluxes,and data from a meteorological tower, are combined with measurements of the onset of convection by boundary-layer wind profilers to explore the timing and behaviour of the transition period. The transition is defined here as the period between sunrise and the time at which the depth ofconvection reaches about 200 m AGL. Diagnostic relationships based on surface heat flux, the temperature difference between 2 m and 200 m, and bulk Richardson number are explored. The transition is foundto be enabled by surface heating relaxing the surface stability, while the warming of the layerbetween 2 m and 200 m is in large part due to shear-driven entrainment.     

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.     

Intermittency and the Exchange of Scalars in the Nocturnal Surface Layer

The determination of nocturnal surface fluxes in low wind conditions is a major problem for micrometeorological studies. The eddy correlation technique, extensively used in field measurements, becomes inappropriate if not enough turbulent activity exists. At the same time, the phenomenon of turbulence intermittency is responsible for the existence of localized events of short duration within which a large fraction of the total nighttime scalar exchange occurs. The scalar flux within a certain intermittent event varies considerably depending on the window used for the flux calculation. In many cases, events with very different time durations occur in the same night, and therefore, the proper determination of the surface flux would require averaging within data windows of different sizes for each event. In this work, the surface exchanges of temperature, moisture and carbon dioxide are analysed at a micrometeorological tower at southern Brazil. Intermittent turbulence is a common occurrence at the location. The analysis shows that the fluxes vary with turbulence intensity and the estimation technique. A variable-window size method for flux estimation is suggested and shown to cause an increase in the magnitude of the nocturnal surface fluxes     

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.     

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.     

Simulation of a Stably Stratified Atmospheric Boundary Layer Using One-Dimensional Turbulence

One-dimensional turbulence (ODT) is a single-column simulation in which vertical motions are represented by an unsteady advective process, rather than their customary representation by a diffusive process. No space or time averaging of mesh-resolved motions is invoked. Molecular-transport scales can be resolved in ODT simulations of laboratory-scale flows, but this resolution of these scales is prohibitively expensive in ODT simulations of the atmospheric boundary layer (ABL), except possibly in small subregions of a non-uniform mesh.Here, two methods for ODT simulation of the ABL on uniform meshes are described and applied to the GABLS (GEWEX Atmospheric Boundary Layer Study; GEWEX is the Global Energy and Water Cycle Experiment) stable boundary-layer intercomparison case. One method involves resolution of the roughness scale using a fixed eddy viscosity to represent subgrid motions. The other method, which is implemented at lower spatial resolution, involves a variable eddy viscosity determined by the local mesh-resolved flow, as in multi-dimensional large-eddy simulation (LES). When run at typical LES resolution, it reproduces some of the key high-resolution results, but its fidelity is lower in some important respects. It is concluded that a more elaborate empirically based representation of the subgrid physics, closely analogous to closures currently employed in LES of the ABL, might improve its performance substantially, yielding a cost-effective ABL simulation tool. Prospects for further application of ODT to the ABL, including possible use of ODT as a near-surface subgrid closure framework for general circulation modeling, are assessed.     

Tethered Balloon Observations Of The Nocturnal Stable Boundary Layer In A Valley

Detailed measurements of profiles of mean and turbulent variablesthrough the nocturnal stable boundary layer over a valley arepresented. Two nights of data are analysed in detail, one with only aweakly stable boundary layer and one with a strongly stable boundarylayer. The weakly stratified night shows high levels of turbulence inwhich the flow remains attached to the valley and the boundary layeracts as a single coherent layer. On the strongly stratified night, twoflow regimes are identified: attached flow, as on the weaklystratified night, and decoupled flow in which the air in the valleybecomes so dynamically stable that there is no turbulent interactionwith the mean flow aloft. Because the valley is sloping, the decoupledlower stagnant air then forms a drainage current. We find that theFroude number evaluated at the hill height, F_H = U(H)/N(H) H,diagnoses the flow regime: when F_H = 2, the flow remainsattached and when F_H 2 the flow in the valley becomesdecoupled from the flow aloft. The dynamics of the flow regimes areshown to be understandable in terms of the gradient Richardson number,which indicates the turbulent mixing. We show that the gradientRichardson number is a key parameter in diagnosing each flow regime.     

Large-Eddy Simulation of air Pollution Dispersion in the Nocturnal Cloud-Topped Atmospheric Boundary Layer

Effects of stratocumulus clouds on the dispersion of contaminants are studied in the nocturnal atmospheric boundary layer. The study is based on a large-eddy simulation (LES) model with a bulk parametrization of clouds. Computations include Lagrangian calculations of atmospheric dispersion of a passive tracer released from point sources at various heights above the ground. The results obtained show that the vertical diffusion is non-Gaussian and depends on the location of a source in the boundary layer.     

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.     

Large-Eddy Simulation Of The Stably Stratified Planetary Boundary Layer

In this work, we study the characteristics of a stably stratifiedatmospheric boundary layer using large-eddy simulation (LES).In order to simulate the stable planetary boundary layer, wedeveloped a modified version of the two-part subgrid-scalemodel of Sullivan et al. This improved version of themodel is used to simulate a highly cooled yet fairly windy stableboundary layer with a surface heat flux of(W)_o = -0.05 m K s^-1and a geostrophic wind speed of U_g = 15 m s^-1.Flow visualization and evaluation of the turbulencestatistics from this case reveal the development ofa continuously turbulent boundary layer with small-scalestructures. The stability of the boundary layercoupled with the presence of a strong capping inversionresults in the development of a dominant gravity wave atthe top of the stable boundary layer that appears to be relatedto the most unstable wave predicted by the Taylor–Goldsteinequation. As a result of the decay of turbulence aloft,a strong-low level jet forms above the boundary layer.The time dependent behaviour of the jet is compared with Blackadar'sinertial oscillation analysis.     

Study of Coherent Structures and Estimation of the Pressure Transport Terms for the Nocturnal Stable Boundary Layer

During the CASES-99 field experiment, three quartz-based microbarographs were installed on the 58-m main tower at the Central Site. These devices measuredabsolute pressure with temperature compensated output at a resolution better than 0.2 Pa and a sampling frequency of 2 s^-1 during the whole campaign. This sampling rate is not adequate to compute turbulent pressure fluxes with the classic averaging method, but the wavelet transform allows flux estimations at a wide range of scales. The resolution of the devices is suitable to study pressure perturbations such as internal gravity waves. The night period of the Intensive Operational Period number 6 (IOP6), where wave-like structures were present, is chosen to illustrate the method. A complete wavelet analysis of pressure recordsand data from sonic anemometers located at the same heights in the tower is performed. Wavelet methods make it possible to identify the relevant scales in the flowand to study the vertical structure of pressure perturbations, including coherent structures and small-scale motions.A study of a simplified turbulence kinetic energy budget equation is made and the contribution of the pressure terms is discussed.     

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.     

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.     

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.     

Different Turbulent Regimes and Vertical Turbulence Structures of the Urban Nocturnal Stable Boundary Layer

Turbulence in the nocturnal boundary layer(NBL) is still not well characterized, especially over complex underlying surfaces. Herein, gradient tower data and eddy covariance data collected by the Beijing 325-m tower were used to better understand the differentiating characteristics of turbulence regimes and vertical turbulence structure of urban the NBL. As for heights above the urban canopy layer(UCL), the relationship between turbulence velocity scale(V_TKE) and wind speed(V) was con...     

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.     

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.