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.     

Simplification and Validation of a Spectral-Tensor Model for Turbulence Including Atmospheric Stability

A spectral-tensor model of non-neutral, atmospheric-boundary-layer turbulence is evaluated using Eulerian statistics from single-point measurements of the wind speed and temperature at heights up to 100 m, assuming constant vertical gradients of mean wind speed and temperature. The model has been previously described in terms of the dissipation rate \(\epsilon \), the length scale of energy-containing eddies \(\mathcal {L}\), a turbulence anisotropy parameter \(\varGamma \), the Richardson number Ri, and the normalized rate of destruction of temperature variance \(\eta _\theta \equiv \epsilon _\theta /\epsilon \). Here, the latter two parameters are collapsed into a single atmospheric stability parameter z / L using Monin–Obukhov similarity theory, where z is the height above the Earth’s surface, and L is the Obukhov length corresponding to \(\{Ri,\eta _\theta \}\). Model outputs of the one-dimensional velocity spectra, as well as cospectra of the streamwise and/or vertical velocity components, and/or temperature, and cross-spectra for the spatial separation of all three velocity components and temperature, are compared with measurements. As a function of the four model parameters, spectra and cospectra are reproduced quite well, but horizontal temperature fluxes are slightly underestimated in stable conditions. In moderately unstable stratification, our model reproduces spectra only up to a scale \(\sim \) 1 km. The model also overestimates coherences for vertical separations, but is less severe in unstable than in stable cases.     

A Hierarchy of Energy- and Flux-Budget (EFB) Turbulence Closure Models for Stably-Stratified Geophysical Flows

Here we advance the physical background of the energy- and flux-budget turbulence closures based on the budget equations for the turbulent kinetic and potential energies and turbulent fluxes of momentum and buoyancy, and a new relaxation equation for the turbulent dissipation time scale. The closure is designed for stratified geophysical flows from neutral to very stable and accounts for the Earth’s rotation. In accordance with modern experimental evidence, the closure implies the maintaining of turbulence by the velocity shear at any gradient Richardson number Ri, and distinguishes between the two principally different regimes: “strong turbulence” at ${Ri \ll 1}$ typical of boundary-layer flows and characterized by the practically constant turbulent Prandtl number Pr _T; and “weak turbulence” at Ri > 1 typical of the free atmosphere or deep ocean, where Pr _T asymptotically linearly increases with increasing Ri (which implies very strong suppression of the heat transfer compared to the momentum transfer). For use in different applications, the closure is formulated at different levels of complexity, from the local algebraic model relevant to the steady-state regime of turbulence to a hierarchy of non-local closures including simpler down-gradient models, presented in terms of the eddy viscosity and eddy conductivity, and a general non-gradient model based on prognostic equations for all the basic parameters of turbulence including turbulent fluxes.     

How do Stability Corrections Perform in the Stable Boundary Layer Over Snow?

We assess sensible heat-flux parametrizations in stable conditions over snow surfaces by testing and developing stability correction functions for two alpine and two polar test sites. Five turbulence datasets are analyzed with respect to, (a) the validity of the Monin–Obukhov similarity theory, (b) the model performance of well-established stability corrections, and (c) the development of new univariate and multivariate stability corrections. Using a wide range of stability corrections reveals an overestimation of the turbulent sensible heat flux for high wind speeds and a generally poor performance of all investigated functions for large temperature differences between snow and the atmosphere above (>10 K). Applying the Monin–Obukhov bulk formulation introduces a mean absolute error in the sensible heat flux of \(6\,\hbox {W m}^{-2}\) (compared with heat fluxes calculated directly from eddy covariance). The stability corrections produce an additional error between 1 and \(5\,\hbox {W m}^{-2}\), with the smallest error for published stability corrections found for the Holtslag scheme. We confirm from previous studies that stability corrections need improvements for large temperature differences and wind speeds, where sensible heat fluxes are distinctly overestimated. Under these atmospheric conditions our newly developed stability corrections slightly improve the model performance. However, the differences between stability corrections are typically small when compared to the residual error, which stems from the Monin–Obukhov bulk formulation.     

On the Velocity Gradient in Stably Stratified Sheared Flows. Part 1: Asymptotic Analysis and Applications

We give a new derivation of the familiar linear relation for the dimensionless velocity gradient in the stably stratified surface layer and provide physical and empirical grounds for its universal applicability in stationary homogeneous turbulence over the whole range of static stabilities from Ri =  0 to very large Ri. Combining this relation with the budget equation for the turbulent kinetic energy we obtain the “equilibrium formulation” of the turbulent dissipation length scale, and recommend it for use in turbulence closure models.     

Outlier Problem in Evaluating Similarity Functions in the Stable Atmospheric Boundary Layer

The gradient-based similarity approach removes turbulent fluxes as governing parameters and replaces them with vertical gradients of mean wind speed and potential temperature. As a result, the gradient Richardson number, Ri, appears as a stability parameter instead of the Monin–Obukhov stability parameter z/L (L is the Obukhov length). The gradient-based scaling is more appropriate for moderate and very stable conditions when the gradients are large and their errors are relatively small whereas z/L becomes ambiguous in these conditions because turbulent fluxes are small. However, the gradient-based formulation is faced with a problem related to the influence of Ri outliers: outliers with high values of Ri can exist in conditions that are really near-neutral. These outliers are mapped into the very stable range in plots in which Ri is the independent variable and may lead to spurious dependencies for bin-averaged data (spurious bin-averaging). This effect is quite large for functions that are steep for the gradient-based scaling. The present study uses the Surface Heat Budget of the Arctic Ocean (SHEBA) data to examine the problem and proposes two methods, conditional analysis and independent binning, to limit the influence of outliers on bin-averaging. A disadvantage of the conditional analysis is associated with eliminating outliers based on criteria that could be considered as subjective. The independent bin-averaging method does not have this disadvantage, but the scatter of the bin-averaged points is higher than for the conditional analysis, rendering data analysis and interpretation difficult.     

Estimation of atmospheric surface layer parameters and numerical simulation using MM5 at Mangalore, West Coast of India

Atmospheric surface layer meteorological observations obtained from 20-m-high meteorological tower at Mangalore, situated along the west coast of India are used to estimate the surface layer scaling parameters of roughness length (z _o) and drag coefficient (C _D), surface layer fluxes of sensible heat and momentum. These parameters are computed using the simple flux–profile relationships under the framework of Monin–Obukhov (M–O) similarity theory. The estimated values of z _o are higher (1.35–1.54 m) than the values reported in the literature (>0.4–0.9 m) probably due to the undulating topography surrounding the location. The magnitude of C _D is high for low wind speed (<1.5 m s^?1) and found to be in the range 0.005–0.03. The variations of sensible heat fluxes (SHF) and momentum fluxes are also discussed. Relatively high fluxes of heat and momentum are observed during typical days on 26–27 February 2004 and 10–11 April 2004 due to the daytime unstable atmospheric conditions. Stable or near neutral conditions prevail after 1700 h IST with negative SHF. A mesoscale model PSU/NCAR MM5 is run using a high-resolution (1 km) grid over the study region to examine the influence of complex topography on the surface layer parameters and the simulated fluxes are compared with estimated values. Spatial variations of the frictional velocity (u _*), C _D, surface fluxes, planetary boundary layer (PBL) height and surface winds are noticed according to the topographic variations in the simulation.     

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin–Obukhov similarity theory is estimated as \(\beta \approx 5.7\) in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations.     

Intermittent Bursting of Turbulence in a Stable Boundary Layer with Low-level Jet

The atmospheric stable boundary layer (SBL) with a low-level jet is simulated experimentally using a thermally stratified wind tunnel. The turbulence structure and flow characteristics are investigated by simultaneous measurements of velocity and temperature fluctuations and by flow visualization. Attention is focused on the effect of strong wind shear due to a low-level jet on stratified boundary layers with strong stability. Occasional bursting of turbulence in the lower portion of the boundary layer can be found in the SBL with strong stability. This bursting originates aloft away from the surface and transports fluid with relatively low velocity and temperature upward and fluid with relatively high velocity and temperature downward. Furthermore, the relationship between the occurrence of turbulence bursting and the local gradient Richardson number (Ri) is investigated. The Ri becomes larger than the critical Ri, Ri_cr = 0.25, in quiescent periods. On the other hand, the Ri number becomes smaller than Ri_cr during bursting events.     

An Efficient Non-iterative Bulk Parametrization of Surface Fluxes for Stable Atmospheric Conditions Over Polar Sea-Ice

In climate and weather prediction models the near-surface turbulent fluxes of heat and momentum and related transfer coefficients are usually parametrized on the basis of Monin–Obukhov similarity theory (MOST). To avoid iteration, required for the numerical solution of the MOST equations, many models apply parametrizations of the transfer coefficients based on an approach relating these coefficients to the bulk Richardson number \(Ri_{b}\). However, the parametrizations that are presently used in most climate models are valid only for weaker stability and larger surface roughnesses than those documented during the Surface Heat Budget of the Arctic Ocean campaign (SHEBA). The latter delivered a well-accepted set of turbulence data in the stable surface layer over polar sea-ice. Using stability functions based on the SHEBA data, we solve the MOST equations applying a new semi-analytic approach that results in transfer coefficients as a function of \(Ri_{b}\) and roughness lengths for momentum and heat. It is shown that the new coefficients reproduce the coefficients obtained by the numerical iterative method with a good accuracy in the most relevant range of stability and roughness lengths. For small \(Ri_{b}\), the new bulk transfer coefficients are similar to the traditional coefficients, but for large \(Ri_{b}\) they are much smaller than currently used coefficients. Finally, a possible adjustment of the latter and the implementation of the new proposed parametrizations in models are discussed.     

On the Scale-dependence of the Gradient Richardson Number in the Residual Layer

We present results of a technique for examining the scale-dependence of the gradient Richardson number, Ri, in the nighttime residual layer. The technique makes use of a series of high-resolution, in situ, vertical profiles of wind speed and potential temperature obtained during CASES-99 in south-eastern Kansas, U.S.A. in October 1999. These profiles extended from the surface, through the nighttime stable boundary layer, and well into the residual layer. Analyses of the vertical gradients of both wind speed, potential temperature and turbulence profiles over a wide range of vertical scale sizes are used to estimate profiles of the local Ri and turbulence structure as a function of scale size. The utility of the technique lies both with the extensive height range of the residual layer as well as with the fact that the sub-metre resolution of the raw profiles enables a metre-by-metre ‘sliding’ average of the scale-dependent Richardson number values over hundreds of metres vertically. The results presented here show that small-scale turbulence is a ubiquitous and omnipresent feature of the residual layer, and that the region is dynamic and highly variable, exhibiting persistent turbulent structure on vertical scales of a few tens of metres or less. Furthermore, these scales are comparable to the scales over which the Ri is less than or equal to the critical value of Ri _c of 0.25, although turbulence is also shown to exist in regions with significantly larger Ri values, an observation at least consistent with the concept of hysteresis in turbulence generation and maintenance. Insofar as the important scale sizes are comparable to or smaller than the resolution of current models, it follows that, in order to resolve the observed details of small Ri values and the concomitant turbulence generation, future models need to be capable of significantly higher resolutions.     

Energy- and Flux-Budget Turbulence Closure Model for Stably Stratified Flows. Part II: The Role of Internal Gravity Waves

We advance our prior energy- and flux-budget (EFB) turbulence closure model for stably stratified atmospheric flow and extend it to account for an additional vertical flux of momentum and additional productions of turbulent kinetic energy (TKE), turbulent potential energy (TPE) and turbulent flux of potential temperature due to large-scale internal gravity waves (IGW). For the stationary, homogeneous regime, the first version of the EFB model disregarding large-scale IGW yielded universal dependencies of the flux Richardson number, turbulent Prandtl number, energy ratios, and normalised vertical fluxes of momentum and heat on the gradient Richardson number, Ri. Due to the large-scale IGW, these dependencies lose their universality. The maximal value of the flux Richardson number (universal constant ≈0.2–0.25 in the no-IGW regime) becomes strongly variable. In the vertically homogeneous stratification, it increases with increasing wave energy and can even exceed 1. For heterogeneous stratification, when internal gravity waves propagate towards stronger stratification, the maximal flux Richardson number decreases with increasing wave energy, reaches zero and then becomes negative. In other words, the vertical flux of potential temperature becomes counter-gradient. Internal gravity waves also reduce the anisotropy of turbulence: in contrast to the mean wind shear, which generates only horizontal TKE, internal gravity waves generate both horizontal and vertical TKE. Internal gravity waves also increase the share of TPE in the turbulent total energy (TTE = TKE + TPE). A well-known effect of internal gravity waves is their direct contribution to the vertical transport of momentum. Depending on the direction (downward or upward), internal gravity waves either strengthen or weaken the total vertical flux of momentum. Predictions from the proposed model are consistent with available data from atmospheric and laboratory experiments, direct numerical simulations and large-eddy simulations.     

Integral Turbulence Statistics Over a Central European City Centre

Atmospheric measurements over 5 years (2005–2010) at two sites in ?ód?, central Poland have been analyzed to develop a better understanding of turbulence in urban areas. Fast response wind velocity, temperature, humidity and CO₂ concentration were measured using sonic anemometers and gas analyzers, placed on narrow masts at 37 and 42 m above the ground. The measurements were used to calculate standard deviations of each parameter, and were then normalized according to local Monin–Obukhov similarity theory and plotted as a function of stability parameter ζ = z′/L. Results for the wind components show typical scaling with a power law with exponent ±1/3 in the free convection limit, and that approaches a constant value close to neutral stratification. For stable conditions, the constant value in the neutral limit remains the same for stability parameters lower than 0.1–0.2, then increases. The normalized standard deviation of temperature fits the ?1/3 law in the free convection limit, approaching a constant value within a stable limit. However, it exhibits hyperbolic characteristics for close to neutral stratification. The normalized standard deviations for humidity and CO₂ concentration exhibit scaling similar to the wind components in the unstable regime and remain constant in the stable domain. The results for the wind components and for temperature are in the range of various functions found in other studies. The absolute values for humidity and CO₂ concentration seem to be slightly higher, but only single examples of such investigations can be found in the literature.     

A Validation of CloudSat and CALIPSO's Temperature,Humidity, Cloud Detection,and Cloud Base Height over the Arctic Marine Cryosphere

Abstract

This study presents a semi‐analytic non‐iterative solution for the Monin‐Obukhov similarity equations under unstable surface conditions. The solution is represented in terms of the non‐dimensional Monin‐Obukhov stability parameter z/L. This parameter is given as a function of the bulk Richardson number and other surface parameters including the heat and momentum roughness lengths which are generally assumed to be different in this formulation. The proposed formulations give results that are both quantitatively and qualitatively consistent with the fully iterated numerical solution for a wide range of surface parameters.     

Air–Sea Interaction Features in the Baltic Sea and at a Pacific Trade-Wind Site: An Inter-comparison Study

A systematic comparison of wind profiles and momentum exchange at a trade wind site outside Oahu, Hawaii and corresponding data from the Baltic Sea is presented. The trade wind data are to a very high degree swell dominated, whereas the Baltic Sea data include a more varied assortment of wave conditions, ranging from a pure growing sea to swell. In the trade wind region swell waves travel predominantly in the wind direction, while in the Baltic, significant cross-wind swells are also present. Showing the drag coefficient as a function of the 10-m wind speed demonstrates striking differences for unstable conditions with swell for the wind-speed range 2 m s^?1 < U ₁₀ < 7 m s^?1, where the trade-wind site drag values are significantly larger than the corresponding Baltic Sea values. In striking contrast to this disagreement, other features studied are surprisingly similar between the two sites. Thus, exactly as found previously in Baltic Sea studies during unstable conditions and swell, the wind profile in light winds (3 m s^?1) shows a wind maximum at around 7–8 m above the water, with close to constant wind speed above. Also, for slightly higher wind speeds (4 m s^?1 < U ₁₀ < 7 m s^?1), the similarity between wind profiles is striking, with a strong wind-speed increase below a height of about 7–8 m followed by a layer of virtually constant wind speed above. A consequence of these wind-profile features is that Monin–Obukhov similarity is no longer valid. At the trade-wind site this was observed to be the case even for wind speeds as high as 10 m s^?1. The turbulence kinetic energy budget was evaluated for four cases of 8–16 30- min periods at the trade-wind site, giving results that agree very well with corresponding figures from the Baltic Sea.     

Relating Mean Radiosounding Profiles to Surface Fluxes for the Very Stable Boundary Layer

A dataset collected during a measurement campaign in the middle of the Po Valley, Italy, is used to investigate the boundary-layer structure in stable conditions. Empirical formulations for temperature and wind profiles derived from Monin–Obukhov similarity theory are used as regression curves to fit radiosounding profiles in the lower half of the boundary-layer. The best fitting parameters of the regression are then compared to the surface turbulent fluxes as measured by a co-located sonic anemometer. This comparison shows significant discrepancies and supports earlier results showing that surface fluxes, in the limit of high stability, are not adequate scalings for mean profiles. The most evident differences are found for cases for which the bulk Richardson number turns out to be quite large. One of the practical consequences is that boundary-layer height diagnostic formulations that mainly rely on surface fluxes are in disagreement with those obtained by inspecting the thermodynamic profiles recorded during the radiosounding ascent. Moreover the incorrect scaling of similarity profiles in stable conditions leads to the erroneous diagnosis of 2-m air temperatures used in numerical weather prediction validation.     

Helicopter-Borne Flux Measurements in the Nocturnal Boundary Layer Over Land – a Case Study

This case study introduces measurements of turbulent fluxes in a nocturnal boundary layer in North Germany with the new helicopter-borne turbulence measurement system HELIPOD, a detailed data analysis and examination in regard of systematic errors of the instrument, and some comparison with local similarity theory and experiments of the past, in order to confirm the occurrence of small vertical turbulent fluxes. The examined nocturnal boundary layer offered excellent conditions to analyse the quality of the measurement system. In this connection, a detailed look at a strong ground-based inversion disclosed small turbulent fluxes with a spectral maximum at ten metres wavelength or less, embedded in intermittent turbulence. For verification of these fluxes, the measurements were compared with well established results from past experiments. Local similarity theory was applied to calculate dimensionless variances of the turbulent quantities, which were found in good agreement with other observations. Since shear and stratification varied significantly on the horizontal flight legs due to global intermittency, a method was developed to determine vertical gradients on a horizontal flight pattern, by use of small fluctuations of the measurement height. With these locally determined gradients, gradient transport theory became applicable and the turbulent diffusivities for heat and momentum, the Richardson number, and the flux Richardson number were estimated within isolated strong turbulent outbursts. Within these outbursts the flux Richardson number was found between 0.1 and 0.2. The functional relationship between the gradient Richardson number and the turbulent Prandtl number agreed well with observations in past experiments and large eddy simulation. The impact of the stratification on the vertical turbulent exchange, as already described for the surface layer using Monin–Obukhov similarity, was analogously observed in the very stably stratified bulk flow when local scaling was applied.     

Stable Boundary-Layer Scaling Regimes: The Sheba Data

Turbulent and mean meteorological data collected at five levels on a 20-m tower over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are analyzed to examine different regimes of the stable boundary layer (SBL). Eleven months of measurements during SHEBA cover a wide range of stability conditions, from the weakly unstable regime to very stable stratification. Scaling arguments and our analysis show that the SBL can be classified into four major regimes: (i) surface-layer scaling regime (weakly stable case), (ii) transition regime, (iii) turbulent Ekman layer, and (iv) intermittently turbulent Ekman layer (supercritical stable regime). These four regimes may be considered as the basic states of the traditional SBL. Sometimes these regimes, especially the last two, can be markedly perturbed by gravity waves, detached elevated turbulence (‘upside down SBL’), and inertial oscillations. Traditional Monin–Obukhov similarity theory works well in the weakly stable regime. In the transition regime, Businger–Dyer formulations work if scaling variables are re-defined in terms of local fluxes, although stability function estimates expressed in these terms include more scatter compared to the surface-layer scaling. As stability increases, the near-surface turbulence is affected by the turning effects of the Coriolis force (the turbulent Ekman layer). In this regime, the surface layer, where the turbulence is continuous, may be very shallow (< 5 m). Turbulent transfer near the critical Richardson number is characterized by small but still significant heat flux and negligible stress. The supercritical stable regime, where the Richardson number exceeds a critical value, is associated with collapsed turbulence and the strong influence of the earth’s rotation even near the surface. In the limit of very strong stability, the stress is no longer a primary scaling parameter.     

Classification of Vertical Wind Speed Profiles Observed Above a Sloping Forest at Nighttime Using the Bulk Richardson Number

Wind speed profiles above a forest canopy relate to scalar exchange between the forest canopy and the atmosphere. Many studies have reported that vertical wind speed profiles above a relatively flat forest can be classified by a stability index developed assuming wind flow above a flat plane. However, can such a stability index be used to classify vertical wind speed profiles observed above a sloping forest at nighttime, where drainage flow often occurs? This paper examines the use of the bulk Richardson number to classify wind speed profiles observed above a sloping forest at nighttime. Wind speed profiles above a sloping forest were classified by the bulk Richardson number Ri _B . Use of Ri _B distinguished between drainage flow, shear flow, and transitional flow from drainage flow to shear flow. These results suggest that Ri _B is useful to interpret nighttime CO₂ and energy fluxes above a sloping forest. Through clear observational evidence, we also show that the reference height should be high enough to avoid drainage-flow effects when calculating Ri _B .