Length Scales of Scalar Diffusion in the Convective Boundary Layer: Laboratory Observations

A laboratory study of scalar diffusion in the convective boundary layer has found results that are consistent with a 1999 large-eddy simulation (LES) study by Jonker, Duynkerke and Cuijpers. For bottom-up and top-down scalars (introduced as ‘infinite’ area sources of passive tracer at the surface and inversion, respectively) the dominant length scale was found to be much larger than the length scale for density fluctuations, the latter being equal to the boundary-layer depth h. The variance of the normalized passive scalar grew continuously with time and its magnitude was about 3–5 times larger for the top-down case than for the bottom-up case. The vertical profiles of the normalized passive scalar variance were found to be approximately constant through the convective boundary layer (CBL) with a value of about 3–8c_*² for bottom-up and 10–50c_*² for top-down diffusion. Finally, there was some evidence of a minimum in the variance and dominant length scale for scalar flux ratios (top-down to bottom-up flux) close to −0.5. All these convection tank results confirm the LES results and support the hypothesis that there is a distinct difference in behaviour between the dynamic and passive variables in the CBL.     

An evaluation of neutral and convective planetary boundary-layer parameterizations relative to large eddy simulations

This paper compares a number of one-dimensional closure models for the planetary boundary layer (PBL) that are currently in use in large-scale atmospheric models. Using the results of a large-eddy simulation (LES) model as the standard of comparison, the PBL models are evaluated over a range of stratifications from free convective to neutral and a range of surface shear stresses. Capping inversion strengths for the convective cases range from weakly to strongly capped. Six prototypical PBL models are evaluated in this study, which focuses on the accuracy of the boundary-layer fluxes of momentum, heat, and two passive scalars. One scalar mimics humidity and the other is a top-down scalar entrained into the boundary layer from above. A set of measures based on the layer-averaged differences of these fluxes from the LES solutions is developed. In addition to the methodological framework and suite of LES solutions, the main result of the evaluation is the recognition that all of the examined PBL parameterizations have difficulty reproducing the entrainment at the top of the PBL, as given by the LES, in most parameter regimes. Some of the PBL models are relatively accurate in their entrainment flux in a subset of parameter regimes. The sensitivity of the PBL models to vertical resolution is explored, and substantive differences are observed in the performance of the PBL models, relative to LES, at low resolution typical of large scale atmospheric models.     

Large-eddy simulations of thermally forced circulations in the convective boundary layer. Part I: A small-scale circulation with zero wind

We have conducted large-eddy simulations (LES) of the atmospheric boundary layer with surface heat flux variations on a spatial scale comparable to the boundary layer depth.We first ran a simulation with a horizontally homogeneous heat flux. In general the results are similar to those of previous large-eddy simulations. The model simulates a field of convective eddies having approximately the correct velocity and spatial scales, and with the crucial property that kinetic energy is transported vigorously upwards through the middle levels. However, the resolved temperature variance is only about half what is observed in the laboratory or the atmosphere. This deficiency — which is shared by many other large-eddy simulations — has dynamic implications, particularly in the pressure/temperature interaction terms of the heat flux budget. Recent simulations by other workers at much higher resolution than ours appear to be more realistic in this respect.The surface heat flux perturbations were one-dimensional and sinusoidal with a wavelength equal to 1.3 times the boundary-layer depth. The mean wind was zero. Results were averaged over several simulations and over time. There is a mean circulation, with ascent over the heat flux maxima (vertical velocity ~0.1w _*) and descent over the heat flux minima. Turbulence is consistently stronger over the heat flux maxima. The horizontal velocity variance components (calculated with respect to the horizontal average) become unequal, implying that convective eddies are elongated parallel to the surface heat flux perturbations.A consideration of the budgets for temperature and velocity suggests several simplifying concepts.The research reported in this paper was conducted while the first author was on study leave at Colorado State University.     

Sensitivity of simulated salinities in a three-dimensional ocean general circulation model to vertical mixing of destabilizing surface fluxes

 We present simulations performed with a three dimensional global ocean general circulation model which show that simulated salinities and amounts of convective mixing are very sensitive to vertical mixing of surface buoyancy fluxes. If, as usual, surface buoyancy fluxes are placed entirely in the topmost model level, our model produces excessive convective mixing in the Southern Ocean. This results in poor stimulated salinity in the Southern Ocean. In this simulation, we assume, as usual, that both surface buoyancy forcing and vertical mixing are homogeneous within each grid cell. If, on the other hand, destabilizing surface fluxes are instantaneously mixed into the subsurface ocean, the model produces much less convective mixing and much more realistic salinities. The vertical mixing of surface buoyancy fluxes performed in this simulation is equivalent to assuming that those fluxes affect only a small fraction of each grid cell, and cause vertical mixing only in that limited area. Our interpretation of these results is that the usual assumption that both surface buoyancy forcing and vertical mixing are uniform within each grid cell has a detrimental effect on model results; these results could be significantly improved by good parametrizations which treat the horizontal inhomogeneity of surface buoyancy forcing and of vertical mixing. Received: 25 February 1998 / Accepted: 9 September 1998     

A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows

A long-standing problem in large-eddy simulations (LES) of the planetary boundary layer (PBL) is that the mean wind and temperature profiles differ from the Monin-Obukhov similarity forms in the surface layer. This shortcoming of LES has been attributed to poor grid resolution and inadequate sub-grid-scale (SGS) modeling. We study this deficiency in PBL LES solutions calculated over a range of shear and buoyancy forcing conditions. The discrepancy from similarity forms becomes larger with increasing shear and smaller buoyancy forcing, and persists even with substantial horizontal grid refinement. With strong buoyancy forcing, however, the error is negligible.In order to achieve better agreement between LES and similarity forms in the surface layer, a two-part SGS eddy-viscosity model is proposed. The model preserves the usual SGS turbulent kinetic energy formulation for the SGS eddy viscosity, but it explicitly includes a contribution from the mean flow and a reduction of the contributions from the turbulent fluctuations near the surface. Solutions with the new model yield increased fluctuation amplitudes near the surface and better correspondence with similarity forms out to a distance of 0.1–0.2 times the PBL depth, i.e., a typical surface-layer depth. These results are also found to be independent of grid anisotropy. The new model is simple to implement and computationally inexpensive.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Footprint estimation for scalar flux measurements in the atmospheric surface layer

The flux footprint is the contribution, per unit emission, of each element of a surface area source to the vertical scalar flux measured at height z _m ; it is equal to the vertical flux from a unit surface point source. The dependence of the flux footprint on crosswind location is shown to be identical to the crosswind concentration distribution for a unit surface point source; an analytic dispersion model is used to estimate the crosswind-integrated flux footprint. Based on the analytic dispersion model, a normalized crosswind-integrated footprint is proposed that principally depends on the single variable z/z _m , where z is a measure of vertical dispersion from a surface source. The explicit dependence of the crosswind-integrated flux footprint on downwind distance, thermal stability and surface roughness is contained in the dependence of z on these variables. By also calculating the flux footprint with a Lagrangian stochastic dispersion model, it is shown that the normalized flux footprint is insensitive to the analytic model assumption of a self-similar vertical concentration profile.The National Center for Atmospheric Research is funded by the National Science Foundation.     

The Influence of Thermally-Induced Mesoscale Circulations on Turbulence Statistics Over an Idealized Urban Area Under a Zero Background Wind

The influence of mesoscale circulations induced by urban-rural differential surface sensible heat flux and roughness on convective boundary-layer (CBL) flow statistics over an isolated urban area has been examined using large-eddy simulation (LES). Results are analyzed when the circulations influence the entire urban area under a zero background wind. For comparison, the CBL flow over an infinite urban area with identical urban surface characteristics under the same background meteorological conditions is generated as a control case (without circulations). The turbulent flow over the isolated urban area exhibits a mix of streaky structure and cellular pattern, while the cellular pattern dominates in the control case. The mixed-layer height varies significantly over the isolated urban area, and can be lower near the edge of the urban area than over the rural area. The vertical profiles of turbulence statistics over the isolated urban area vary horizontally and are dramatically different from the control case. The turbulent kinetic energy (TKE) sources include wind shear, convergence, and buoyancy productions, compared to only buoyancy production in the control case. The normalized vertical velocity variance is reduced compared to the control case except in the central urban area where it is little affected. The low-level flow convergence is mainly responsible for the enhanced horizontal velocity variance in the central urban area, while wind shear is responsible for the additional local maximum of the horizontal velocity variance near the middle of the CBL outside the central area. Parameterizations in the prognostic equation for TKE used in mesoscale models are evaluated against the LES results over the isolated urban area. We also discuss conditions under which the urban-induced circulations occur and when they may affect the entire urban area. Given that urban-induced circulations can influence the entire urban area within hours for an urban area of a realistic size, it is inappropriate to directly apply empirical relations of turbulence statistics derived under horizontally-homogenous flow conditions to an urban area.     

Evaluating Source Area Contributions from Aircraft Flux Measurements Over Heterogeneous Land Using Large-Eddy Simulation

The estimation of spatial patterns in surface fluxes from aircraft observations poses several challenges in the presence of heterogeneous land cover. In particular, the effects of turbulence on scalar transport and the different behaviour of passive (e.g. water vapour) versus active (e.g. temperature) scalars may lead to large uncertainties in the source area/flux-footprint estimation for sensible (H) and latent (LE) heat-flux fields. This study uses large-eddy simulation (LES) of the land–atmosphere interactions to investigate the atmospheric boundary-layer (ABL) processes that are likely to create differences in airborne-estimated H and LE footprints. We focus on 32~m altitude aircraft flux observations collected over a study site in central Oklahoma during the Southern Great Plains experiment in 1997 (SGP97). Comparison between the aircraft data and traditional model estimates provide evidence of a difference in source area for turbulent sensible and latent heat fluxes. The LES produces reasonable representations of the observed fluxes, and hence provides credible evidence and explanation of the observed differences in the H and LE footprints. Those differences can be quantified by analyzing the change in the sign of the spatial correlation of the H and LE fields provided by the LES model as a function of height. Dry patterns in relatively moist surroundings are able to generate strong, but localized, sensible heating. However, whereas H at the aircraft altitude is still in phase with the surface, LE presents a more complicated connection to the surface as the dry updrafts force a convergence of the surrounding moist air. Both the observational and LES model evidence support the concept that under strongly advective conditions, H and LE measured at the top of the surface layer (≈50 m) can be associated with very different upwind source areas, effectively contradicting surface-layer self-similarity theory for scalars. The results indicate that, under certain environmental conditions, footprint models will need to predict differing source area/footprint contributions between active (H) and passive (LE) scalar fluxes by considering land-surface heterogeneity and ABL dynamics.     

Large-Eddy Simulation of Katabatic Flows

A large-eddy simulation model with rotated coordinates and an open boundary is used to simulate the characteristics of katabatic flows over simple terrain. Experiments examine the effects of cross winds on the development of the slope-flow boundary layer for a steep (20°) slope and the role of drainage winds in preventing turbulence collapse on a gentle slope (1°). For the steep flow cases, comparisons between model average boundary-layer velocity, temperature deficit, and turbulence kinetic energy budget terms and tower observations show reasonable agreement. Results for different cross slope winds show that as the cross slope winds increase, the slope flow deepens faster and behaves more like a weakly stratified, sheared boundary layer. Analysis of the momentum budget shows that near the surface the flow is maintained by a balance between downslope buoyancy forcing and vertical turbulence flux from surface drag. Above the downslope jet, the turbulence vertical momentum flux reverses sign and acceleration of the flow by buoyancy is controlled by horizontal advection of slower moving ambient air. The turbulence budget is dominated by a balance between shear production and eddy dissipation, however, buoyancy and pressure transport both are significant in reducing the strength of turbulence above the jet. Results from the gentle slope case show that even a slight terrain variation can lead to significant drainage winds. Comparison of the gentle slope case with a flat terrain simulation indicates that drainage winds can effectively prevent the formation of very stable boundary layers, at least near the top of sloping terrain.     

The Effect of Chessboard Variability of the Surface Fluxes on the Aggregated Turbulence Fields in a Convective Atmospheric Surface Layer

To what degree the variability of surface features can be identified in the turbulent signals observed in the atmospheric boundary layer is still an unresolved problem. This was investigated by conducting an analytical experiment for a one-dimensional 'chessboard'-type surface-flux distribution on the basis of local free convection scaling. The results showed that, due to their nonlinear dependency on the surface fluxes, the dimensionless gradients of the mean quantities and the dimensionless standard deviations are altered by the surface-flux variability. Furthermore, passive scalars, such as humidity, are considerably more sensitive to surface variability than the main active scalar, temperature. However, the response of the gradients of the mean quantities is fairly negligible in the range of variability studied herein as compared to that of the standard deviations, which were found to be more sensitive to the surface-flux variability. In addition, the phase difference between the active and the passive scalar flux distribution strongly affects the passive scalar turbulence. This dissimilarity between passive and active scalars, or between passive scalars when their source distributions are different, brings into question the use of variance methods for the measurement of a scalar flux, such as evaporation, over variable surfaces. The classical Bowen ratio method, which depends on the validity of the Reynolds analogy for the vertical gradients of the mean quantities, was shown to be relatively more robust. However, under conditions of strong surface variability, it can also be expected to fail.     

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.     

Dispersion of a Passive Scalar Within and Above an Urban Street Network

The transport of a passive scalar from a continuous point-source release in an urban street network is studied using direct numerical simulation (DNS). Dispersion through the network is characterized by evaluating horizontal fluxes of scalar within and above the urban canopy and vertical exchange fluxes through the canopy top. The relative magnitude and balance of these fluxes are used to distinguish three different regions relative to the source location: a near-field region, a transition region and a far-field region. The partitioning of each of these fluxes into mean and turbulent parts is computed. It is shown that within the canopy the horizontal turbulent flux in the street network is small, whereas above the canopy it comprises a significant fraction of the total flux. Vertical fluxes through the canopy top are predominantly turbulent. The mean and turbulent fluxes are respectively parametrized in terms of an advection velocity and a detrainment velocity and the parametrization incorporated into a simple box-network model. The model treats the coupled dispersion problem within and above the street network in a unified way and predictions of mean concentrations compare well with the DNS data. This demonstrates the usefulness of the box-network approach for process studies and interpretation of results from more detailed numerical simulations.     

Large Wind Farms and the Scalar Flux over an Heterogeneously Rough Land Surface

The influence of surface heterogeneities extends vertically within the atmospheric surface layer to the so-called blending height, causing changes in the fluxes of momentum and scalars. Inside this region the turbulence structure cannot be treated as horizontally homogeneous; it is highly dependent on the local surface roughness, the buoyancy and the horizontal scale of heterogeneity. The present study analyzes the change in scalar flux induced by the presence of a large wind farm installed across a heterogeneously rough surface. The change in the internal atmospheric boundary-layer structure due to the large wind farm is decomposed and the change in the overall surface scalar flux is assessed. The equilibrium length scale characteristic of surface roughness transitions is found to be determined by the relative position of the smooth-to-rough transition and the wind turbines. It is shown that the change induced by large wind farms on the scalar flux is of the same order of magnitude as the adjustment they naturally undergo due to surface patchiness.     

Scaling Variances of Scalars in a Convective Boundary Layer Under Different Entrainment Regimes

For the presentation and analysis of atmospheric boundary-layer (ABL) data, scales are used to non-dimensionalise the observed quantities and independent variables. Usually, the ABL height, surface sensible heat flux and surface scalar flux are used. This works well, so long as the absolute values of the entrainment ratio for both the scalar and temperature are similar. The entrainment ratio for temperature naturally ranges from −0.4 to −0.1. However, the entrainment ratio for passive scalars can vary widely in magnitude and sign. Then the entrainment flux becomes relevant as well. The only customary scalar scale that takes into account both the surface flux and the entrainment flux is the bulk scalar scale, but this scale is not well-behaved for large negative entrainment ratios and for an entrainment ratio equal to −1. We derive a new scalar scale, using previously published large-eddy simulation results for the convective ABL. The scale is derived under the constraint that scaled scalar variance profiles are similar at those heights where the variance producing mechanisms are identical (i.e., either near the entrainment layer or near the surface). The new scale takes into account that scalar variance in the ABL is not only related to the surface flux of that scalar, but to the scalar entrainment flux as well. Furthermore, it takes into account that the production of variance by the entrainment flux is an order of magnitude larger than the production of variance by the surface flux (per unit flux). Other desirable features of the new scale are that it is always positive (which is relevant when scaling standard deviations) and that the scaled variances are always of order 1–10.     

Dissimilarity of Scalar Transport in the Convective Boundary Layer in Inhomogeneous Landscapes

A land-surface model (LSM) is coupled with a large-eddy simulation (LES) model to investigate the vegetation-atmosphere exchange of heat, water vapour, and carbon dioxide (CO₂) in heterogeneous landscapes. The dissimilarity of scalar transport in the lower convective boundary layer is quantified in several ways: eddy diffusivity, spatial structure of the scalar fields, and spatial and temporal variations in the surface fluxes of these scalars. The results show that eddy diffusivities differ among the three scalars, by up to 10–12%, in the surface layer; the difference is partly attributed to the influence of top-down diffusion. The turbulence-organized structures of CO₂ bear more resemblance to those of water vapour than those of the potential temperature. The surface fluxes when coupled with the flow aloft show large spatial variations even with perfectly homogeneous surface conditions and constant solar radiation forcing across the horizontal simulation domain. In general, the surface sensible heat flux shows the greatest spatial and temporal variations, and the CO₂ flux the least. Furthermore, our results show that the one-dimensional land-surface model scheme underestimates the surface heat flux by 3–8% and overestimates the water vapour and CO₂ fluxes by 2–8% and 1–9%, respectively, as compared to the flux simulated with the coupled LES-LSM.     

Simulating Dispersion in the Evening-Transition Boundary Layer

We investigate dispersion in the evening-transition boundary layer using large-eddy simulation (LES). In the LES, a particle model traces pollutant paths using a combination of the resolved flow velocities and a random displacement model to represent subgrid-scale motions. The LES is forced with both a sudden switch-off of the surface heat flux and also a more gradual observed evolution. The LES shows ‘lofting’ of plumes from near-surface releases in the pre-transition convective boundary layer; it also shows the subsequent ‘trapping’ of releases in the post-transition near-surface stable boundary layer and residual layer above. Given the paucity of observations for pollution dispersion in evening transitions, the LES proves a useful reference. We then use the LES to test and improve a one-dimensional Lagrangian Stochastic Model (LSM) such as is often used in practical dispersion studies. The LSM used here includes both time-varying and skewed turbulence statistics. It is forced with the vertical velocity variance, skewness and dissipation from the LES for particle releases at various heights and times in the evening transition. The LSM plume spreads are significantly larger than those from the LES in the post-transition stable boundary-layer trapping regime. The forcing from the LES was thus insufficient to constrain the plume evolution, and inclusion of the significant stratification effects was required. In the so-called modified LSM, a correction to the vertical velocity variance was included to represent the effect of stable stratification and the consequent presence of wave-like motions. The modified LSM shows improved trapping of particles in the post-transition stable boundary layer.     

Comparison of Kansas Data with High-resolution Large-eddy Simulation Fields

The surface layer of an atmospheric boundary layer (ABL) is most accessible to field measurements and hence its ensemble-mean structure has been well established. The Kansas field measurements were the first detailed study of this layer, providing numerous benchmark statistical profiles for a wide range of stability states. Large-eddy simulation (LES), in contrast, is most suitable for studying the mixed layer of the ABL where the energy-containing range of the vertical velocity field is well resolved. In the surface layer, typical large-eddy simulations barely resolve the energy-containing vertical-velocity fields and hence do not provide sufficient data for a detailed analysis.We carried out a nested-mesh simulation of a moderately convective ABL (-z_i/L = 8) in which the lower 6% of the boundary layer had an effective grid resolution of 512³. We analyze the LES fields above the 6th vertical grid level (z = 23 m) where the vertical velocity field has a well formed inertial subrange, for a detailed comparison with the Kansas results. Various terms in the budgets of turbulent kinetic energy, temperature variance, Reynolds stress, temperature flux, and some higher order moments are compared. The agreement is generally quite good; however, we do observe certain discrepancies, particularly in the terms involving pressure fluctuations.     

Scalar Fluxes Near a Tall Building in an Aligned Array of Rectangular Buildings

Scalar dispersion from ground-level sources in arrays of buildings is investigated using wind-tunnel measurements and large-eddy simulation (LES). An array of uniform-height buildings of equal dimensions and an array with an additional single tall building (wind tunnel) or a periodically repeated tall building (LES) are considered. The buildings in the array are aligned and form long streets. The sensitivity of the dispersion pattern to small changes in wind direction is demonstrated. Vertical scalar fluxes are decomposed into the advective and turbulent parts and the influences of wind direction and of the presence of the tall building on the scalar flux components are evaluated. In the uniform-height array turbulent scalar fluxes are dominant, whereas the tall building produces an increase of the magnitude of advective scalar fluxes that yields the largest component. The presence of the tall building causes either an increase or a decrease to the total vertical scalar flux depending on the position of the source with respect to the tall building. The results of the simulations can be used to develop parametrizations for street-canyon dispersion models and enhance their capabilities in areas with tall buildings.     

Turbulent Statistics of Neutrally Stratified Flow Within and Above a Sparse Forest from Large-Eddy Simulation and Field Observations

Turbulent statistics of neutrally stratified shear-driven flow within and above a sparse forest canopy are presented from a large-eddy simulation (LES) and compared with those from observations within and above a deciduous forest with similar height and foliage density. First- and second-order moments from the LES agree with observations quite well. Third-order moments from the LES have the same sign and similar vertical patterns as those from the observations, but the LES yields smaller magnitudes of such higher-order moments. Turbulent spectra and cospectra from the LES agree well with observations above the forest. However, at the highest frequencies, the LES spectra have steeper slopes than observations. Quadrant and conditional analyses of the LES resolved-scale flow fields also agree with observations. For example, both LES and observation find that sweeps are more important than ejections for the transport of momentum within the forest, while inward and outward interaction contributions are both small, except near the forest floor. The intermittency of the transport of momentum and scalar increases with depth into the forest. Finally, ramp structures in the time series of a passive scalar at multiple levels within and above the forest show similar features to those measured from field towers. Two-dimensional (height-time cross-section) contours of the passive scalar and wind vectors show sweeps and ejections, and the characteristics of the static pressure perturbation near the ground resemble those deduced from field tower-based measurements. In spite of the limited grid resolution (2 m × 2 m × 2 m) and domain size (192 m × 192 m × 60 m) used in this LES, we demonstrate that the LES is capable of resolving the most important characteristics of the turbulent flow within and above a forest canopy.