On the turbulent Prandtl number in the stable atmospheric boundary layer

This study focuses on the behaviour of the turbulent Prandtl number, Pr _t , in the stable atmospheric boundary layer (SBL) based on measurements made during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA). It is found that Pr _t increases with increasing stability if Pr _t is plotted vs. gradient Richardson number, Ri; but at the same time, Pr _t decreases with increasing stability if Pr _t is plotted vs. flux Richardson number, Rf, or vs. ζ = z/L. This paradoxical behaviour of the turbulent Prandtl number in the SBL derives from the fact that plots of Pr _t vs. Ri (as well as vs. Rf and ζ) for individual 1-h observations and conventional bin-averaged values of the individual quantities have built-in correlation (or self-correlation) because of the shared variables. For independent estimates of how Pr _t behaves in very stable stratification, Pr _t is plotted against the bulk Richardson number; such plots have no built-in correlation. These plots based on the SHEBA data show that, on the average, Pr _t decreases with increasing stability and Pr _t < 1 in the very stable case. For specific heights and stabilities, though, the turbulent Prandtl number has more complicated behaviour in the SBL.     

Wind-Tunnel Study Of Atmospheric Stable Boundary Layers Over A Rough Surface

Wind-tunnel simulations of theatmospheric stable boundary layer (SBL) developedover a rough surface were conducted by using athermally stratified wind tunnel at the Research Institutefor Applied Mechanics (RIAM), Kyushu University. Thepresent experiment is a continuation of the workcarried out in a wind tunnel at Colorado StateUniversity (CSU), where the SBL flows were developed over asmooth surface. Stably stratified flows were createdby heating the wind-tunnel airflow to a temperature ofabout 40–50°and by cooling the test-section floor toa temperature of about 10°. To simulate therough surface, a chain roughness was placed over thetest-section floor. We have investigated the buoyancyeffect on the turbulent boundary layer developed overthis rough surface for a wide range of stability,particularly focusing on the turbulence structure andtransport process in the very stable boundary layer.The present experimental results broadly confirm theresults obtained in the CSU experiment with the smoothsurface, and emphasizes the following features: thevertical profiles of turbulence statistics exhibitdifferent behaviour in two distinct stability regimes with weak and strong stability,corresponding to the difference in the verticalprofiles of the local Richardson number. The tworegimes are separated by the critical Richardsonnumber. The magnitudes in turbulence intensities andturbulent fluxes for the weak stability regime aremuch greater than those of the CSU experiments becauseof the greater surface roughness. For the very stableboundary layer, the turbulent fluxes of momentum andheat tend to vanish and wave-like motions due to theKelvin–Helmholtz instability and the rolling up andbreaking of those waves can be observed. Furthermore,the appearance of internal gravity waves is suggestedfrom cross-spectrum analyses.     

The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer

Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin–Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a ‘critical value’ of about 0.20–0.25, the inertial subrange associated with the Richardson–Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson–Kolmogorov energy cascade weakens; therefore, the applicability of local Monin–Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri _cr and Rf < Rf _cr. However, it is found that Rf _cr  =  0.20–0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin–Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson–Kolmogorov cascade) have been filtered out.     

Characteristics of intermittent turbulence in the upper stable boundary layer over Greenland

The experiment IGLOS (Investigation of the Greenland Boundary Layer Over Summit) was conducted in June and July 2002 in the central plateau of the Greenland inland ice. The German research aircraft Polar2, equipped with the turbulence measurement system Meteopod, was used to investigate turbulence and radiation flux profiles near research station “Summit Camp”. Aircraft measurements are combined with measurements of radiation fluxes and turbulent quantities made from a 50 m tower at Summit Camp operated by Eidgenössische Technische Hochschule (ETH) Zürich. During all six flight missions, well-developed stable boundary layers were found. Even in high-wind conditions, the surface inversion thickness did not exceed roughly 100 m. The turbulent height of the stable boundary layer (SBL) was found to be much smaller than the surface inversion thickness. Above the surface layer, significant turbulent fluxes occurred only intermittently in intervals on the order of a few kilometres. Turbulent event fraction in the upper SBL shows the same dependence on gradient Richardson number as reported for near-surface measurements. Clear-air longwave radiation divergence was always found to contribute significantly to the SBL heat budget. In low-wind cases, radiative cooling even turned out to be dominant.     

Analysis of Atmospheric Boundary Layer Height Characteristics over the Arctic Ocean Using the Aircraft and GPS Soundings

Utilizing aircraft sounding data collected from the Surface Heat Budget of the Arctic Ocean (SHEBA, 1998) campaign, the authors evaluated commonly-used profile methods for Arctic ABL height estimation by validating against the’true’ABL height from aircraft sounding profiles, where ABL height is defined as the top of the layer with significant turbulence. Furthermore, the best performing method was used to estimate ABL height from the one-year GPS soundings obtained during SHEBA (October 1997-October 1998). It was found that the temperature gradient method produces a reliable estimate of ABL height. Additionally, the authors determined optimal threshold values of temperature gradient for stable boundary layer (SBL) and convective boundary layer (CBL) to be 6.5 K/100 m and 1.0 K/100 m, respectively. The maximum ABL height during the year was 1150 m occurred in May. Median values of Arctic ABL height in May, June, July, and August were 400 m, 430 m, 180 m, and 320 m, respectively. Arctic ABL heights are clearly higher in the spring than in the summer.     

Nocturnal Boundary-Layer Regimes

This study analyzes turbulence data collected over a grassland site in the nocturnal boundary layer. Examination of the dependence of the nocturnal boundary layer on stability suggests three regimes: a) the weakly stable case, b) a transition stability regime where many of the variables change rapidly with increasing stability and c) the very stable case. The value of z/L where the downward heat flux is a maximum defines the stability boundary between the weakly stable and transition regimes, where L is the Obukhov length. In the present analysis, the downward heat flux reaches a maximum at z/L approximately equal to 0.05 for 10 m, although comparison with other data indicates that this is not a universal value. For weaker stability, the heat flux decreases with decreasing z/L due to weaker temperature fluctuations. In the transition stability regime, the heat flux decreases rapidly with increasing stability due to restriction of vertical velocity fluctuations by the increasing stratification.For weakly stable conditions, the variances scale according to Monin-Obukhov similarity theory. For very stable conditions, the variances are contaminated by non-turbulent horizontal motions and do not follow the scaling laws. An alternative length scale based on variances is developed which explains more of the variance of the transfer coefficients compared to the Obukhov length.     

A Study of the Stable Boundary Layer Based on a Single-Column K-Theory Model

We document numerical experiments with a single-column, high-resolution model of the stable boundary layer. The model resolves the logarithmic layer, and does not require inverting the Monin–Obukhov similarity functions in order to calculate the surface fluxes. The turbulence closure is based on the K-theory approach, with a new form of stability functions of the Richardson number, evaluated by using the Surface Heat Budget of the Arctic Ocean (SHEBA) and the Cooperative Atmosphere-Surface Exchange Study (CASES-99) data. A comparison with two, high-resolution large-eddy simulation models shows very good agreement. The reported numerical experiments test the effects of shear, surface cooling, the Coriolis parameter, subsidence, and baroclinicity. The time evolution of the drag coefficient, the heat-transfer coefficient, and the cross-isobar angle is also evaluated.     

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.     

Stratified Atmospheric Boundary Layers

Various features of different stability regimes of the stable boundary layer are discussed. Traditional layering is examined in terms of the roughness sublayer, surface layer, local similarity, z-less stratification and the region near the boundary-layer top. In the very stable case, the strongest turbulence may be detached from the surface and generated by shear associated with a low level jet, gravity waves or meandering motions. In this case, similarity theory and the traditional concept of a boundary-layer break down. The elevated turbulence may intermittently recouple to the surface. Inability to adequately measure turbulent fluxes in very stable conditions limits our knowledge of this regime.     

SHEBA flux–profile relationships in the stable atmospheric boundary layer

Measurements of atmospheric turbulence made during the Surface Heat Budget of the Arctic Ocean Experiment (SHEBA) are used to examine the profile stability functions of momentum, φ _m , and sensible heat, φ _h , in the stably stratified boundary layer over the Arctic pack ice. Turbulent fluxes and mean meteorological data that cover different surface conditions and a wide range of stability conditions were continuously measured and reported hourly at five levels on a 20-m main tower for 11 months. The comprehensive dataset collected during SHEBA allows studying φ _m and φ _h in detail and includes ample data for the very stable case. New parameterizations for φ _m (ζ) and φ _h (ζ) in stable conditions are proposed to describe the SHEBA data; these cover the entire range of the stability parameter ζ = z/L from neutral to very stable conditions, where L is the Obukhov length and z is the measurement height. In the limit of very strong stability, φ _m follows a ζ ^1/3 dependence, whereas φ _h initially increases with increasing ζ, reaches a maximum at ζ ≈ 10, and then tends to level off with increasing ζ. The effects of self-correlation, which occur in plots of φ _m and φ _h versus ζ, are reduced by using an independent bin-averaging method instead of conventional averaging.     

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.     

Early Warning Signals for Regime Transition in the Stable Boundary Layer: A Model Study

The evening transition is investigated in an idealized model for the nocturnal boundary layer. From earlier studies it is known that the nocturnal boundary layer may manifest itself in two distinct regimes, depending on the ambient synoptic conditions: strong-wind or overcast conditions typically lead to weakly stable, turbulent nights; clear-sky and weak-wind conditions, on the other hand, lead to very stable, weakly turbulent conditions. Previously, the dynamical behaviour near the transition between these regimes was investigated in an idealized setting, relying on Monin–Obukhov (MO) similarity to describe turbulent transport. Here, we investigate a similar set-up, using direct numerical simulation; in contrast to MO-based models, this type of simulation does not need to rely on turbulence closure assumptions. We show that previous predictions are verified, but now independent of turbulence parametrizations. Also, it appears that a regime shift to the very stable state is signaled in advance by specific changes in the dynamics of the turbulent boundary layer. Here, we show how these changes may be used to infer a quantitative estimate of the transition point from the weakly stable boundary layer to the very stable boundary layer. In addition, it is shown that the idealized, nocturnal boundary-layer system shares important similarities with generic non-linear dynamical systems that exhibit critical transitions. Therefore, the presence of other, generic early warning signals is tested as well. Indeed, indications are found that such signals are present in stably stratified turbulent flows.     

Modelling of the Evolving Stable Boundary Layer

A single-column model of the evolving stable boundary layer (SBL) is tested for self-similar properties of the flow and effects of ambient forcing. The turbulence closure of the model is diagnostic, based on the K-theory approach, with a semi-empirical form of the mixing length, and empirical stability functions of the Richardson number. The model results, expressed in terms of local similarity scales, are universal functions, satisfied in the entire SBL. Based on similarity expression, a realizability condition is derived for the minimum allowable turbulent heat flux in the SBL. Numerical experiments show that the development of “horse-shoe” shaped, fixed-elevation hodographs in the interior of the SBL around sunrise is controlled by effects imposed by surface thermal forcing.     

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.     

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.     

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.     

黄土高原复杂地形受中尺度运动影响的稳定边界层湍流特征

利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University,简称SACOL)2008年12月观测资料,研究了稳定边界层湍流特征.使用涡动相关资料研究湍流通量时,定义湍流的平均时间τ内的中尺度运动是造成湍流统计量变化范围大的主要原因,稳定情形? τ取几十秒至几分钟.对梯度理查森数大于0.3的强稳定情形的湍流尺度分解(MRD)谱分析表明,感热通量在112.4～449.9 s存在谱隙,尺度大于谱隙的中尺度运动造成了通量观测资料离散性大,甚至有支配性影响.动量通量的谱隙在112.4～224.9 s之间.弱风时,中尺度运动的影响更大,垂直风速标准差以0.1的比率随中尺度风速变化;垂直风速标准差同广义风速表现出很好的相关性,并随着广义风速消失而消失.三维风速标准差与摩擦速度呈很好的线性关系,垂直、水平、横风风速的无量纲标准差分别为1.35、2.54、2.21.对湍流动能的研究发现,在梯度理查森数大于0.3的条件下,仍然存在连续的湍流.以湍动能为依据,分析了湍流的平稳时间长度,其长度随稳定度变化而变化,2008年12月7～11日从133.5 s变化到856.2 s,湍流平稳时间长度反映了中尺度运动的发生频率.     

Aircraft-Based Measurements Of Turbulence Structures In The Katabatic Flow Over Greenland

Turbulence structures in the katabatic flow in the stable boundary layer (SBL) over the ice sheet are studied for two case studies with high wind speeds during the aircraft-based experiment KABEG (Katabatic wind and boundary layer front experiment around Greenland) in the area of southern Greenland. The aircraft data allow the direct determination of turbulence structures in the katabatic flow. For the first time, this allows the study of the turbulence structure in the katabatic wind system over the whole boundary layer and over a horizontal scale of 80 km.The katabatic flow is associated with a low-level jet (LLJ), with maximum wind speeds up to 25 m s^-1. Turbulent kinetic energy (TKE) and the magnitude of the turbulent fluxes show a strong decrease below the LLJ. Sensible heat fluxes at the lowest level have values down to -25 W m^-2. Latent heat fluxes are small in general, but evaporation values of up to +13 W m^-2 are also measured. Turbulence spectra show a well-defined inertial subrange and a clear spectral gap around 250-m wavelength. While turbulence intensity decreases monotonously with height above the LLJ for the upper part of the slope, high spectral intensities are also present at upper levels close to the ice edge. Normalized fluxes and variances generally follow power-law profiles in the SBL.Terms of the TKE budget are computed from the aircraft data. The TKE destruction by the negative buoyancy is found to be very small, and the dissipation rate exceeds the dynamical production.     

Turbulence Structure of Stable Boundary Layers with a Near-Linear Temperature Profile

By using a thermally stratified wind tunnel, we have successfullysimulated stably stratified boundary layers (SBL), in which the meantemperature increases upward almost linearly. We have investigated the flow structure and the effects of near-linearstable stratification on the transfer of momentum and heat. Thevertical profiles of turbulence quantities exhibit different behaviour in two distinct stability regimes of the SBLflows with weak and strong stability. For weak stability cases, theturbulent transfer of momentum and heat is basically similar to that for neutral turbulent boundary layers, although it is weakenedwith increasing stability. For strong stability cases, on the other hand,the time-mean transfer is almost zero over the whole boundary-layer depth.However, the instantaneous turbulent transfer frequently occurs in bothgradient and counter-gradient directions in the lower part of the boundary layer. This is due to the Kelvin–Helmholtz (K–H) shear instability and therolling up and breaking of K–H waves. Moreover, the internal gravity wavesare observed in the middle and upper parts of all stable boundary layers.