Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-Eddy Simulations. Part II: Forced Convection

Forced convection in a quasi-steady atmospheric boundary layer is investigated based on a large-eddy simulation (LES) model. The performed simulations show that in the upper portion of the mixed layer the dimensionless (in terms of mixed layer scales) vertical gradients of temperature, humidity, and wind velocity depend on the dimensionless height z/z _i and the Reech number Rn. The peak values of variances and covariances at the top of the mixed layer, scaled in terms of the interfacial scales, are functions of the interfacial Richardson number Ri. As a result expressions for the entrainment rates, in the case when the interfacial layer has a finite depth, and a condition for the presence of moistening or drying regimes in the mixed layer, are derived. Profiles of dimensionless scalar moments in the mixed layer are proposed to be expressed in terms of two empirical similarity functions F _m and F _i , dependent on dimensionless height z/z _i , and the interfacial Richardson number Ri. The obtained similarity expressions adequately approximate the LES profiles of scalar statistics, and properly represent the impact of stability, shear, and entrainment. They are also consistent with the parameterization proposed for free convection in the first part of this paper.     

Formulation of the thermal internal boundary layer in a mesoscale model

Mesoscale models using a non-local K-scheme for parameterization of boundary-layer processes require an estimate of the planetary boundary layer (PBL) height z _i at all times. In this paper, two-dimensional sea-breeze experiments are carried out to evaluate three different formulations for the advective contribution in the z _i prognostic equation of Deardorff (1974).Poor representation of the thermal internal boundary layer in the sea breeze is obtained when z _i is advected by the wind at level z _i . However, significantly better results are produced if the mean PBL wind is used for the advecting velocity, or if z _i is determined simply by checking for the first sufficiently stable layer above the ground.A Lagrangian particle model is used to demonstrate the effect of each formulation on plume dispersion by the sea breeze.     

Scaling Properties of Temperature Spectra and Heat-Flux Cospectra in the Surface Friction Layer Beneath an Unstable Outer Layer

Temperature variance and temperature power spectra in the unstable surface layer have always presented a problem to the standard Monin-Obukhov similarity model. Recently that problem has intensified with the demonstration by Smedman et al. (2007, Q J Roy Meteorol Soc 133: 37–51) that temperature spectra and heat-flux cospectra can have two distinct peaks in slightly unstable conditions, and by McNaughton et al. (2007, Nonlinear Process Geophys 14: 257–271) who showed that the wavenumber of the peak of temperature spectra in a convective boundary layer (CBL), closely above the surface friction layer (SFL), can be sensitive to the CBL depth, z _i. Neither the two-peak form at slight instability nor the dependence of peak position on z _i at large instability is compatible with the Monin-Obukhov model. Here we examine the properties of temperature spectra and heat-flux cospectra from between these extremes, i.e. from within the unstable SFL, in two experiments. The analysis is based on McNaughton’s model of the turbulence structure in the SFL. According to this model, heat is transported through most of the SFL by sheet plumes, created by the action of impinging outer eddies. The smallest and most effective of these outer eddies have sizes that scale on SFL depth, z _s. The z _s-scale eddies and plumes are organised within the overall convection pattern in the CBL, and in turn they organise the motion of smaller eddies within the SFL, whose sizes scale on height, z. The main experimental results are: (1) the peak amplitudes of the temperature spectra in the SFL are collapsed with a scaling factor (z_sz)^1/3e_o^2/3{(z_{\rm s}z)^{1/3}\varepsilon_{\rm o}^{2/3}} divided by the square of the surface temperature flux, where e_o{\varepsilon_{\rm o}} is the dissipation rate of turbulent energy in the outer CBL (above the SFL); (2) the peak wavenumbers of the temperature spectra are collapsed with the mixed length scale (z _i z _s)^1/2; (3) the peak wavenumbers of the heat-flux cospectra are collapsed with the doubly-mixed length scale (z _i z _s)^1/4 z ^1/2; (4) for z/z _s < 0.03, the peak in the cospectrum is replaced by another peak at a wavenumber about a magnitude larger. This peak’s position scales on z; (5) all these findings are consistent with the observations of Smedman et al.     

Velocity and Surface Shear Stress Distributions Behind a Rough-to-Smooth Surface Transition: A Simple New Model

A simple new model is proposed to predict the distribution of wind velocity and surface shear stress downwind of a rough-to-smooth surface transition. The wind velocity is estimated as a weighted average between two limiting logarithmic profiles: the first log law, which is recovered above the internal boundary-layer height, corresponds to the upwind velocity profile; the second log law is adjusted to the downwind aerodynamic roughness and local surface shear stress, and it is recovered near the surface, in the equilibrium sublayer. The proposed non-linear form of the weighting factor is equal to ln(z/z ₀₁)/ln(δ _i /z ₀₁), where z, δ _i and z ₀₁ are the elevation of the prediction location, the internal boundary-layer height at that downwind distance, and the upwind surface roughness, respectively. Unlike other simple analytical models, the new model does not rely on the assumption of a constant or linear distribution for the turbulent shear stress within the internal boundary layer. The performance of the new model is tested with wind-tunnel measurements and also with the field data of Bradley. Compared with other existing analytical models, the proposed model shows improved predictions of both surface shear stress and velocity distributions at different positions downwind of the transition.     

A numerical study of daily transitions in the convective boundary layer

We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z _m (defined by the maximum of the potential temperature gradient) and its standard deviation σ _m , the mixed-layer height z _i (defined by the minimum of the potential temperature flux), entrainment velocity W _e, and the entrainment flux H _i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w_*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ_*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.     

The characteristics of turbulent velocity components in the surface layer under convective conditions

It is proposed that the ratios of the standard deviations of the horizontal velocity components to the friction velocity in the surface layer under convective conditions depend only onz _i /L wherez _i is the height of the lowest inversion andL is the Monin-Obukhov length. This hypothesis is tested by using observations from several data sets over uniform surfaces and appears to fit the data well. Empirical curves are fitted to the observations which have the property that at largez _i /-L, the standard deviations become proportional tow _*, the convective scaling velocity.Fluctuations of vertical velocity obtained from the same experiments scale withz/L, wherez is the height above the surface, in good agreement with Monin-Obukhov theory.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Direct numerical simulation of a vigorously heated low Reynolds number convective boundary layer

Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The results are obtained by directly solving the three-dimensional time-dependen Navier-Stokes equations with the Boussinesq buoyancy approximation, resolving all relevant scales of motion (no turbulence closure is needed). The flow is capped by a stable temperature inversion and heated from below at a rate that produces an inversion-height to Obukhov-length ratio −z_i/L_* = 32. Temperature and velocity variance profiles are found to agree well with those from an earlier vigorously heated under-resolved computation at higher Reynolds number, and with experimental data of Deardorff and Willis (Boundary-Layer Meteorol., 32: 205–236, 1985). Significant helicity is found in the layer, and helical convection patterns of the scale of the inversion height are observed.     

Measurements of spray at an air-water interface

The spray content in the surface boundary layer above an air—water interface was determined by a series of measurements at various feteches and wind speeds in a laboratory facility. The droplet flux density N(z) can be described in terms of the scaling flux density N_* and von Karman constant K throguh the equation, N(z)/N_* = −(1/K) ln(z/z_0d) where z is height above the mean water level and z_0d is the droplet boundary layer thickness. N_* is given by a unique relationship in terms of the roughness Reynolds number u_*σ/ν where σ is the root-mean-square surface displacement. Spray inception occurred for u_* 0.3. The dominant mode of spray generation in the present and most other laboratory tests, as well as in available field data, appears to be bubble bursting.     

On the Time Evolution of the Turbulent Kinetic Energy Spectrum for Decaying Turbulence in the Convective Boundary Layer

Our focus is the time evolution of the turbulent kinetic energy for decaying turbulence in the convective boundary layer. The theoretical model with buoyancy and inertial transfer terms has been extended by a source term due to mechanical energy and validated against large-eddy simulation data. The mechanical effects in a boundary layer of height z _i at a convective surface-layer height z = 0.05z _i are significant in the time evolution of the vertical component of the spectrum, i.e. they enhance the decay time scale by more than an order of magnitude. Our findings suggest that shear effects seem to feedback to eddies with smaller wavenumbers, preserving the original shape of the spectrum, and preventing the spectrum from shifting towards shorter wavelengths. This occurs in the case where thermal effects only are considered.     

Mean Flow and Turbulence Characteristics in an Urban Roughness Sublayer

In this study, a detailed model of an urban landscape has been re-constructed inthe wind tunnel and the flow structure inside and above the urban canopy has beeninvestigated. Vertical profiles of all three velocity components have been measuredwith a Laser-Doppler velocimeter, and an extensive analysis of the measured meanflow and turbulence profiles carried out. With respect to the flow structure inside thecanopy, two types of velocity profiles can be distinguished. Within street canyons,the mean wind velocities are almost zero or negative below roof level, while closeto intersections or open squares, significantly higher mean velocities are observed.In the latter case, the turbulent velocities inside the canopy also tend to be higherthan at street-canyon locations. For both types, turbulence kinetic energy and shearstress profiles show pronounced maxima in the flow region immediately above rooflevel.Based on the experimental data, a shear-stress parameterization is proposed, inwhich the velocity scale, u_s, and length scale, z_s, are based on the level and magnitude of the shear stress peak value. In order to account for a flow region inside the canopy with negligible momentum transport, a shear stress displacement height, d_s, is introduced. The proposed scaling and parameterization perform well for the measured profiles and shear-stress data published in the literature.The length scales derived from the shear-stress parameterization also allowdetermination of appropriate scales for the mean wind profile. The roughnesslength, z₀, and displacement height, d₀, can both be described as fractions of the distance, z_s - d_s, between the level of the shear-stress peak and the shear-stress displacement height. This result can be interpreted in such a way that the flow only feels the zone of depth z_s - d_s as the roughness layer. With respect to the lower part of the canopy (z < d_s) the flow behaves as a skimming flow. Correlations between the length scales z_s and d_s and morphometric parameters are discussed.The mean wind profiles above the urban structure follow a logarithmic windlaw. A combination of morphometric estimation methods for d₀ and z₀ with wind velocity measurements at a reference height, which allow calculation of the shear-stress velocity, u_*, appears to be the most reliable and easiest procedure to determine mean wind profile parameters. Inside the roughnesssublayer, a local scaling approach results in good agreement between measuredand predicted mean wind profiles.     

On Transpiration and Soil Moisture Content Sensitivity to Soil Hydrophysical Data

Sensitivity of evapotranspiration E and root zone soil moisture content θ to the parameterization of soil water retention Ψ(θ) and soil water conductivity K(Ψ), as well as to the definition of field capacity soil moisture content, is investigated by comparing Psi1-PMSURF and Theta-PMSURF models. The core of PMSURF (Penman–Monteith Surface Fluxes) consists of a 3-layer soil moisture prediction module based on Richard’s equation in combination with the Penman–Monteith concept for estimating turbulent heat fluxes. Psi1- PMSURF and Theta-PMSURF differ only in the parameterization of the moisture availability function F_ma. In Psi1,F_ma is parameterized by using Ψ(θ) and K(Ψ) hydrophysical functions; in Theta, F_ma is parameterized by using hydrophysical parameters: the field capacity θ_f and wilting point θ_w soil moisture contents. Both Psi1 and Theta are based on using soil hydrophysical data, that is, there is no conceptual difference between them in the parameterization of E even if in Psi1F_ma depends on 12 parameters, while in Theta only on two soil/vegetation parameters. Sensitivity tests are performed using the Cabauw dataset. Three soil datasets are used: the vG (van Genuchten), CH/vG (Clapp and Hornberger/van Genuchten) and CH/PILPS (Clapp and Hornberger/Project for Intercomparison of Land-surface Parameterization Schemes) datasets. The vG dataset is used in van Genuchten’s parameterization, while in Clapp and Hornberger’s the CH/vG and CH/PILPS datasets are used. It is found that the consistency of soil hydrophysical data in the simulation of transpiration is quite important. The annual sum of E obtained by Psi1E^Psi1, differs from the annual sum of E obtained by Theta, E^Theta, because of the inconsistency between the fitting parameters of Ψ(θ) and K(Ψ) and the θ_f, and not because of the differencies in the parameterization of F_ma. Further, θ_f can be estimated not only on the basis of using soil hydrophysical functions (the θ_f so obtained is θ^Soil_f) but also on the basis of analysing the transpiration process (the θ_fso obtained is θ^tr_f). θ^tr_f values estimated from the condition E^Theta ≈ E^Psi1 are in acceptable accordance with the θ^Soil_f values proposed by Wösten and co-workers. The results are useful in optimizing the parameterization of transpiration in land-surface schemes.     

A comparison of turbulence statistics in the surface layer over plant canopies with those over several other surfaces

Turbulence statistics, including higher order moments, in the surface layer over plant canopies were compared with those observed over several different surfaces, using a nondimensional height (z – d)/z ₀: The values of (z – d)/z ₀extend over a very wide range from 10 over plant canopies to 10⁷ over the ocean. Several properties such as intensities of turbulence and skewness factors show a remarkable height-dependency in the air layer below (z – d)/z ₀ = 10², which is supposed to be much influenced by the underlying surface. In that layer, some peculiar phenomena, such as a downward energy transport and positive flux of shear stress, are frequently observed.     

Can Water Vapour Raman Lidar Resolve Profiles of Turbulent Variables in the Convective Boundary Layer?

High-resolution water vapour measurements made by the Atmospheric Radiation Measurement (ARM) Raman lidar operated at the Southern Great Plains Climate Research Facility site near Lamont, Oklahoma, U.S.A. are presented. Using a 2-h measurement period for the convective boundary layer (CBL) on 13 September 2005, with temporal and spatial resolutions of 10 s and 75 m, respectively, spectral and autocovariance analyses of water vapour mixing ratio time series are performed. It is demonstrated that the major part of the inertial subrange was detected and that the integral scale was significantly larger than the time resolution. Consequently, the major part of the turbulent fluctuations was resolved. Different methods to retrieve noise error profiles yield consistent results and compare well with noise profiles estimated using Poisson statistics of the Raman lidar signals. Integral scale, mixing-ratio variance, skewness, and kurtosis profiles were determined including error bars with respect to statistical and sampling errors. The integral scale ranges between 70 and 130 s at the top of the CBL. Within the CBL, up to the third order, noise errors are significantly smaller than sampling errors and the absolute values of turbulent variables, respectively. The mixing-ratio variance profile rises monotonically from ≈0.07 to ≈3.7 g² kg⁻² in the entrainment zone. The skewness is nearly zero up to 0.6 z/z _i , becomes −1 around 0.7–0.8 z/z _i , crosses zero at about 0.95 z/z _i , and reaches about 1.7 at 1.1 z/z _i (here, z is the height and z _i is the CBL depth). The noise errors are too large to derive fourth-order moments with sufficient accuracy. Consequently, to the best of our knowledge, the ARM Raman lidar is the first water vapour Raman lidar with demonstrated capability to retrieve profiles of turbulent variables up to the third order during daytime throughout the atmospheric CBL.     

A suggested refinement for O'Brien's convective boundary layer Eddy exchange coefficient formulation

With observational data collected and interpreted by Crane et al. (1977), the adequacy of the O'Brien polynomial to represent the exchange profile of heat and pollution in a convective boundary layer is examined and a refinement suggested. Also, it is shown that the height of the surface layer, h=0.04 z _i , developed by Blackadar and Tennekes (1968) for a neutrally stratified boundary layer (with z _z =0.25u _*/f) appears to be equally valid for the convective boundary layer where z _i , defined as the top of the mixed layer, is used.     

On the parameterization of precipitation in warm clouds

The Sundqvist parameterization for warm rain production by autoconversion processes as the function of cloud liquid water mixing ratio m is tested by defining a realistic ‘driving’ profile m(z) for a maritime low, warm stratocumulus cloud, and comparing with various recent observations. The results show that the parameterization is acceptable, especially after tuning its rain collection constant C₁. It is somewhat sensitive to the vertical resolution of the host model, though. Extending the calculations by considering typical cloud and raindrop size spectra, extra variables such as drizzle amounts and droplet effective radii (forced by the bulk Sundqvist rain rates) could be estimated by numerical integration. Also, these seem to agree fairly well with the available observations.     

地形重力波拖曳参数化对热带气旋强度和路径预报影响的研究

本文在GRAPES_TMM（Global/Regional Assimilation and Prediction System for Tropical Mesoscale Model）——中国南海台风模式版（面向南海和东南亚）中发展和引进了KA95（Kim and Arakawa,1995）地形重力波拖曳参数化方案（GWDO）,并对2012年主要的9个登陆台风进行了试验对比研究,考察了不同标准Richardson数（Ri_c）的GWDO试验对台风路径和强度预报的影响。结果表明,在引入地形重力波拖曳参数化过程后,模式对台风登陆时路径和强度的预报能力均要有提高,对台风预报时长越长,GWDO的影响也更为显著。对双台风“SAOLA”和“DAMREY”试验结果表明,GWDO对台风外围距台风中心150 km的对流层中下层风速减弱较为明显,减弱了GRAPES区域模式对台风强度预报偏强的现象,对台风强度长时间预报改善更为明显。不同标准Ri_c对重力波拖曳力的计算较为敏感,当Ri_c取1.0时,动能迅速的在低层被频散,能量无法有效地上传;Ri_c取0.25时,大部分的能量在中高层被频散。总的来说,Ri_c取0.75时对台风路径和强度预报改进更为显著,其结果可为业务预报提供指导意义。     

Momentum Absorption in Rough-Wall Boundary Layers with Sparse Roughness Elements in Random and Clustered Distributions

A wind-tunnel experiment has been used to investigate momentum absorption by rough surfaces with sparse random and clustered distributions of roughness elements. An unusual (though longstanding) method was used to measure the boundary-layer depth δ and friction velocity u _* and thence to infer the functional relationship z ₀/h = f(λ) between the normalised roughness length z ₀/ h and the roughness density λ (where z ₀ is the roughness length and h the mean height of the roughness elements). The method for finding u _* is based on fitting the velocity defect in the outer layer to a functional form for the dimensionless velocity-defect profile in a canonical zero-pressure-gradient boundary layer. For the conditions investigated here, involving boundary layers over sparse roughness with strong local heterogeneity, this velocity-defect-law method is found to be more robust than several alternative methods for finding u _* (uw covariance, momentum integral and slope of the logarithmic velocity profile).The experimental results show that, (1) there is general agreement in the relationship z ₀/h = f(λ) between the present experiment with random arrays and other wind-tunnel experiments with regular arrays; (2) the main effect of clustering is to increase the scatter in the z ₀/h = f(λ) relationship, through increased local horizontal heterogeneity; (3) this scatter obscures any trend in the z ₀/h = f(λ) relationship in response to clustering; and (4) the agreement between the body of wind-tunnel data (taken as a whole) and field data is good, though with scatter for which it is likely that a major contribution stems from local horizontal heterogeneity in the field.     

The impact of the PBL scheme and the vertical distribution of model layers on simulations of Alpine foehn

Summary This paper investigates the influence of the planetary boundary-layer (PBL) parameterization and the vertical distribution of model layers on simulations of an Alpine foehn case that was observed during the Mesoscale Alpine Programme (MAP) in autumn 1999. The study is based on the PSU/NCAR MM5 modelling system and combines five different PBL schemes with three model layer settings, which mainly differ in the height above ground of the lowest model level (z ₁). Specifically, z ₁ takes values of about 7 m, 22 m and 36 m, and the experiments with z ₁ = 7 m are set up such that the second model level is located at z = 36 m. To assess if the different model setups have a systematic impact on the model performance, the simulation results are compared against wind lidar, radiosonde and surface measurements gathered along the Austrian Wipp Valley. Moreover, the dependence of the simulated wind and temperature fields at a given height (36 m above ground) on z ₁ is examined for several different regions. Our validation results show that at least over the Wipp Valley, the dependence of the model skill on z ₁ tends to be larger and more systematic than the impact of the PBL scheme. The agreement of the simulated wind field with observations tends to benefit from moving the lowest model layer closer to the ground, which appears to be related to the dependence of lee-side flow separation on z ₁. However, the simulated 2 m-temperatures are closest to observations for the intermediate z ₁ of 22 m. This is mainly related to the fact that the simulated low-level temperatures decrease systematically with decreasing z ₁ for all PBL schemes, turning a positive bias at z ₁ = 36 m into a negative bias at z ₁ = 7 m. The systematic z ₁-dependence is also observed for the temperatures at a fixed height of 36 m, indicating a deficiency in the self-consistency of the model results that is not related to a specific PBL formulation. Possible reasons for this deficiency are discussed in the paper. On the other hand, a systematic z ₁-dependence of the 36-m wind speed is encountered only for one out of the five PBL schemes. This turns out to be related to an unrealistic profile of the vertical mixing coefficient. Correspondence: Günther Z?ngl, Meteorologisches Institut der Universitat München, 80333 München, Germany     

A Modelling Study of Flux Imbalance and the Influence of Entrainment in the Convective Boundary Layer

It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u _*/w _*, where u _* and w _* are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z _i , where z _i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.     

Sublayer-Stanton numbers of heat and matter for aerodynamically smooth surfaces: basic considerations and evaluation

Summary Sublayer-Stanton numbers, B_i, of heat and matter for the interfacial sublayer over aerodynamically smooth surfaces determined for forced convective conditions by elementary and numerical integration are reviewed and evaluated. The results are based on Roth’s modified Heisenberg model for the spectral energy transfer in the equilibrium range under locally isotropic conditions and the approaches of Reichardt, Elser, Deissler, van Driest, Rannie, Sheppard, and Spalding for the normalized eddy diffusivity K_m/v. The results substantiate that with the exception of Sheppard’s K_m/v-approach all formulations are appropriate to provide sublayer-Stanton numbers with a sufficient degree of accuracy. From a theoretical point of view the K_m/v-relationships of Roth, Reichardt, van Driest, and Spalding are to be preferred, when a turbulent Prandtl number Pr_t = 1 is presumed. Since within the framework of mesoscale meteorological modelling numerical integration techniques would consume too much CPU-time because of the large number of near-wall grid points, a parameterization formula for the sublayer-Stanton number is proposed and evaluated. Compared to the B_i ⁻¹-results obtained by numerical integration, this kind of parameterization leads to a relative error of less than 5 percent for roughness Reynolds numbers, η_r, ranging from 30 to 600. Received January 2001 Revised November 5, 2001