Similarity of scalars under stable conditions

The question whether two different scalars have the same behaviour in the surface layer under stable conditions is investigated. Similarity of two scalars is defined in terms of the equality of their corresponding dimensionless Monin-Obukhov similarity functions. Previous theoretical and experimental results concerning the issue are briefly reviewed: they are found to be contradictory. An analytical derivation of the square of the correlation coefficient between two scalars is obtained based on the correlation structure of the turbulent dissipation functions for stable conditions, when it can be assumed that the divergence of the vertical transport of scalar variance/covariance is negligible. The resulting expression elucidates some earlier conflicting results, and helps to establish the equality of the similarity functions for all scalars in the stable surface layer. A statistical analysis in the time domain is also performed using temperature and humidity turbulence data measured in nocturnal stable conditions during FIFE-89. Our results, both from the analytical derivation and the statistical analysis of turbulence data, confirm that under validity of the Monin-Obukhov similarity theory assumptions, the corresponding similarity functions for temperature and humidity are equal to within the statistical uncertainty of the measurements. An important consequence is that the eddy diffusivities of temperature and humidity are also equal.Now at SIMEPAR — Sistema Meteorológico do Paraná and UFPR — Universidade Federal do Paraná, Caixa Postal 318 CEP 80001-970, Curitiba PR Brazil.     

A Study of Correlations of Scalar Quantities in the Atmospheric Surface Layer

Standard deviations for vertical velocity and scalar quantities, such as temperature, T, and specific humidity, q, were analyzed on the basis of Monin-Obukhov (M-O) similarity theory in the atmospheric surface layer. The correlation coefficient between scalar quantities T and q, R_Tq, was derived from the similarity functions and can be expressed as the ratio of B_T/B_q (B_T≤ B_q), where parameter B is the value of the normalized standard deviation of any scalar quantity at neutral conditions.     

Use of variance techniques to measure dry air-surface exchange rates

The variances of fluctuations of scalar quantities can be measured and interpreted to yield indirect estimates of their vertical fluxes in the atmospheric surface layer. Strong correlations among scalar fluctuations indicate a similarity of transfer mechanisms, which is utilized in some of the variance techniques. The ratios of the standard deviations of two scalar quantities, for example, can be used to estimate the flux of one if the flux of the other is measured, without knowledge of atmospheric stability. This is akin to a modified Bowen ratio approach. Other methods such as the normalized standard-deviation technique and the correlation-coefficient technique can be utilized effectively if atmospheric stability is evaluated and certain semi-empirical functions are known. In these cases, iterative calculations involving measured variances of fluctuations of temperature and vertical wind velocity can be used in place of direct flux measurements. For a chemical sensor whose output is contaminated by non-atmospheric noise, covariances with fluctuations of scalar quantities measured with a very good signal-to-noise ratio can be used to extract the needed standard deviation. Field measurements have shown that many of these approaches are successful for gases such as ozone and sulfur dioxide, as well as for temperature and water vapor, and could be extended to other trace substances. In humid areas, it appears that water vapor fluctuations often have a higher degree of correlation to fluctuations of other trace gases than do temperature fluctuations; this makes water vapor a more reliable companion or reference scalar. These techniques provide some reliable research approaches but, for routine or operational measurement, they are limited by the need for fast-response sensors. Also, all variance approaches require some independent means to estimate the direction of the flux.This research has been funded as part of the National Acid Precipitation Assessment Program by the U.S. Environmental Protection Agency through IAGDW89930069-01 to the U.S. Department of Energy.     

Convective scaling of the average dissipation rate of temperature variance in the atmospheric surface layer

The flux of sensible heat from the land surface is related to the average rate of dissipation of temperature fluctuations in the atmospheric surface layer through the temperature variance budget equation. In many cases it is desirable to estimate the heat flux from measurement or inference of the dissipation rate. Here we study how the dissipation rate scales with atmospheric stability, using three inertial range methods to calculate the dissipation rate: power spectra, second order structure functions, and third order structure functions. Experimental data are analyzed from a pair of field experiments, during which turbulent fluctuations of velocity and temperature were measured over a broad range of neutral and unstable atmospheric flows. It is shown that the temperature dissipation rate scales with a single convective power law continuously from near-neutral to strongly unstable stratification. The dissipation scaling is found to nearly match production in the near-neutral region, but to be consistently lower than production in the more convective regimes. The convective scaling is shown to offer a simplified means of computing sensible heat flux from the dissipation rate of temperature variance.Also at Johns Hopkins University, Baltimore, MarylandAlso at Los Alamos National Laboratory, Los Alamos, New Mexico.     

Radiative Effects on Temperature in the Stable Surface Layer

The interaction between longwave radiation and temperature fluctuations plays a role in the dissipation of temperature variance. This interaction is most easily described by spectral models of atmospheric turbulence and a spectral radiative dissipation function which gives the intensity of the damping at each radiative wavelength and wavenumber k. We have used a Corrsin–Pao closure for the spectral budgets of turbulent kinetic energy and temperature to study the coupling of radiation to turbulence. The spectral radiative dissipation function and a related integral have been fitted by analytical approximations with the correct asymptotic behavior. This resulted in a simple analytical formula for the dimensionless temperature spectrum as a function of Monin-Obukhov stability, and a new dimensionless parameter describing the relative importance of radiation in the temperature spectral budget. The radiative effects both on the temperature spectrum and on the dimensionless temperature variance can then be calculated. Based on typical values of the radiative dimensionless parameters for the surface layer, we conclude that radiative dissipation is probably negligible there.     

On the correlation between turbulent velocity and temperature derivatives in the atmospheric surface layer

Velocity and temperature derivatives were obtained at a height of 4 m in the atmospheric surface layer above land. With the assumption of local isotropy, these measurements are used to obtain some statistics of the turbulent energy and scalar dissipation fields. These statistics include the variances of the logarithms of the scalar and velocity dissipation fields and the correlation between these logarithms. When used in conjunction with the hypotheses for fluctuations in turbulent dissipation rates of Obukhov and Kolmogorov, the statistics suggest that the dependence of the flatness factor of temperature derivative on the turbulent Rynolds number R is not as large as that which had been previously reported in the literature. The experimental data indicate a R ^0.5 dependence for the kurtosis of the temperature derivative and a R ^0.15 variation for the strain rate-scalar dissipation correlation.     

On Monin–Obukhov Similarity In The Stable Atmospheric Boundary Layer

Atmospheric measurements from several field experiments have been combined to develop a better understanding of the turbulence structure of the stable atmospheric boundary layer. Fast response wind velocity and temperature data have been recorded using 3-dimensional sonic anemometers, placed at severalheights (1 m to 4.3 m) above the ground. The measurements wereused to calculate the standard deviations of the three components of the windvelocity, temperature, turbulent kinetic energy (TKE) dissipation andtemperature variance dissipation. These data were normalized and plottedaccording to Monin–Obukhov similarity theory. The non-dimensional turbulencestatistics have been computed, in part, to investigate the generalapplicability of the concept of z-less stratification for stable conditions. From the analysis of a data set covering almost five orders ofmagnitude in the stability parameter = z/L (from near-neutral tovery stable atmospheric stability), it was found that this concept does nothold in general. It was only for the non-dimensional standard deviation oftemperature and the average dissipation rate of turbulent kinetic energythat z-less behaviour has been found. The other variables studied here(non-dimensional standard deviations of u, v, and w velocity components and dissipation of temperature variance) did not follow the concept of z-less stratification for the very stable atmospheric boundary layer. An imbalance between production and dissipation of TKE was found for the near-neutral limit approached from the stable regime, which matches with previous results for near-neutral stability approached from the unstable regime.     

Temperature structure in the atmospheric surface layer

Production, transport and dissipation terms in the temperature variance equation have been measured in the atmospheric surface layer. The transport term is, within the experimental uncertainty, negligible. The dissipation term, determined by assuming local isotropy, is approximately equal to production under near-neutral conditions. For moderately unstable conditions, the ratio of production to dissipation is 1.4. The resulting imbalance in the budget is attributed to the inequality between the three components of the dissipation term. The Kolmogorov constant for temperature is found to be about 0.8.     

Some Basic Properties of the Surrogate Subgrid-Scale Heat Flux in the Atmospheric Boundary Layer

The development of improved subgrid-scale (SGS) models for large-eddy simulation of scalar transport in the atmospheric boundary layer requires an improved understanding of basic properties of the SGS fluxes. High frequency atmospheric wind speed and temperature data sampled at a height of 1.7 m are used to measure SGS heat fluxes and dissipation of temperature variance, by means of one-dimensional filtering and invoking Taylor's hypothesis. Conditional averaging is used to isolate interesting features of the SGS signals, and to relate them to the large-scale characteristics of the flow, such as the presence of coherent structures. Both mean and conditionally averaged SGS quantities are compared with those obtained using a standard eddy-diffusivity model. Within the limitations imposed by the one-dimensional data analysis, we observe that the model appears unable to reproduce important features of the real signals, such as the negative dissipation of temperature variance associated with strong negative resolved temperature gradients due to the ejection of warm air under unstable atmospheric stability conditions.     

Vertical Structure in the Marine Atmospheric Boundary Layer and its Implication for the Inertial Dissipation Method

The structure of the marine atmospheric boundarylayer and the validity ofMonin–Obukhov similarity theory over the seahave been investigated using longterm measurements. Three levels of turbulencemeasurements (at 10 m, 18 mand 26 m) at Östergarnsholm in themiddle of the Baltic Sea have beenanalysed. The results show that turbulentparameters have a strong dependenceon the actual height due to wave influence.The wind profile and thus thenormalised wind gradient are very sensitiveto wave state. The lower part of theboundary layer can be divided into three heightlayers, a wave influenced layerclose to the surface, a transition layer andan undisturbed ordinary surfacelayer; the depth of the layers is determinedby the wave state. This heightstructure can, however, not be found for thenormalised dissipation, which is onlya function of the stability, except duringpronounced swell where the actualheight also has to be accounted for. Theresults have implications for the heightvariation of the turbulent kinetic energy(TKE) budget. Thus, the imbalancebetween production and dissipation willalso vary with height according to thevariation of wave state. This, in turn,will of course have strong implicationsfor the inertial dissipation method, inwhich a parameterisation of the TKEbudget is used.     

Fractal properties of temperature fluctuations in the convective surface layer

Temperature fluctuations in a convective surface layer were investigated. Box counting analysis was performed to investigate fractal properties of surfaces of constant temperature and was performed on sets of points obtained by setting thresholds on detrended records. Results indicate that surfaces of constant temperature have fractal properties for thresholds far from the mean. Estimated fractal dimensions of one-dimensional cuts through these surfaces varied between 0.23 and 0.66, increasing with threshold value approaching the mean temperature. For thresholds close to the mean, no fractal behavior was found. Asymmetry in results for thresholds above and below the mean temperature was attributed to the asymmetry between updrafts and downdrafts in the convective surface layer.The temperature dissipation rate (TD) was also investigated. It was found to be strongly intermittent with large fluctuations of the intermittency exponent. Moments were analyzed in order to investigate multifractal properties of TD. Results indicate scaling in the range of 50–1000 (where is the Kolmogorov scale) and multifractal properties resembling those observed for passive scalar dissipation in laboratory flows.     

Accuracy of moments of velocity and scalar fluctuations in the atmospheric surface layer

A detailed accuracy analysis is presented for moments, up to order four, of both velocity (horizontal u and vertical w) and scalar (temperature and humidity q) fluctuations, as well as of the products uw, w and wq, in the atmospheric surface layer. The high-order moments and integral time scales required for this analysis are evaluated from data obtained at a height of about 5 m above the ocean surface under stability conditions corresponding to Z/L \- –0.05. Measured moments and probability density functions of some of the individual fluctuations show departures from Gaussianity, but these are sufficiently small to enable good estimates to be obtained using Gaussian instead of measured moments. For the products, the assumption of joint Gaussianity for individual fluctuations provides a reasonable, though somewhat conservative, estimate for the integration times required. The concept of Reynolds number similarity implies that differences in integration time requirements for flows at different Reynolds numbers arise exclusively from differences in integral time scales. A first approximation to the integral time scales relevant to atmospheric flows is presented.     

The Temperature–Humidity Covariance in the Marine Surface Layer: A One-dimensional Analytical Model

An analytical model that predicts how much of the temperature–humidity covariance within the marine atmospheric surface layer (ASL) originates just above the ASL and just near the surface is proposed and tested using observations from the Risø Air Sea Experiment (RASEX). The model is based on a simplified budget for the two-scalar covariance that retains three basic terms: production, dissipation, and vertical transport. Standard second-order closure formulations are employed for the triple moments and the dissipation terms, and the canonical mixing length for the closure model is assumed linear with height (z) from the surface. Despite the poor performance of the gradient–diffusion closure in reproducing the measured triple moment, the overall covariance model was shown to be sufficiently robust to these assumptions. One of the main findings from the analytical treatment is the origin of the asymmetry in how the top and bottom boundary conditions affect the two-scalar covariance in the ASL. The analytical model reveals that ‘bottom-up’ boundary-condition variations scale with \(z^{-\sqrt{a}}\) , while ‘top-down’ variations scale with \(z^{\sqrt{a}}\) , where a is a constant that can be derived from similarity and closure constants. The genesis of this asymmetry stems from the flux-transport term but is modulated by the dissipation, and persists even in the absence of any inhomogeneity in the local production function. It is shown that the local production function acts to adjust the relative proportions of these two boundary conditions with weights that vary with the Obukhov length. The findings here do not provide ‘finality’ to the discussions on the covariance between humidity and temperature or the role of entrainment in modulating the turbulence within the ASL. Rather, they are intended to guide new hypotheses about interpretations of existing field data and identify needs for future field and numerical experiments.     

Quantifying the Influence of Random Errors in Turbulence Measurements on Scalar Similarity in the Atmospheric Surface Layer

The influence of random errors in turbulence measurements on scalar similarity for temperature, water vapour, \(\hbox {CO}_{2}\), and \(\hbox {NH}_{3}\) is investigated using two eddy-covariance datasets collected over a lake and a cattle feedlot. Three measures of scalar similarity, namely, the similarity constant in the flux–variance relationship, the correlation coefficient between two scalars and the relative transport efficiency, are examined. The uncertainty in the similarity constant \(C_{s}\) in the flux–variance relationship resulting from random errors in turbulence measurements is quantified based on error propagation analyses and a Monte-Carlo sampling method, which yields a distribution instead of a single value for \(C_{s}\). For different scalars, the distributions of \(C_{s}\) are found to significantly overlap, implying that scalars are transported similarly under strongly unstable conditions. The random errors in the correlation coefficients between scalars and the relative transport efficiencies are also quantified through error propagation analyses, and they increase as the atmosphere departs from neutral conditions. Furthermore, the correlation coefficients between three scalars (water vapour, \(\hbox {CO}_{2}\), and \(\hbox {NH}_{3}\)) are statistically different from unity while the relative transport efficiencies are not, which highlights the difference between these two measures of scalar similarity. The results suggest that uncertainties in these measures of scalar similarity need to be quantified when using them to diagnose the existence of dissimilarity among different scalars.     

Numerical simulations of the tropical air-sea planetary boundary layer

A one-dimensional numerical model of the planetary boundary layer was used to investigate thermal and kinetic energy budgets. The simulation experiments were based on two sets of data. The first set was based on a ‘typical’ June with climatological data extracted for the oceanic region slightly northeast of Barbados. The second set used data from the third phase of project BOMEX, for approximately the same area and time of year as the first set. Comparison with observations of three simulated elements (viz., sea surface temperature and wind and humidity at 6 m) which are important in determining the near-interface energy transports shows that:

1. the model is capable of realistic simulations of both ‘typical’ conditions, and conditions for a specific four-day period;
2. the model is capable of realistically simulating the differences between prevailing values of these parameters in the two cases (‘typical’ and specific four-day period).

The simulated interface fluxes are those of incoming and outgoing short- and long-wave radiation; transmitted radiation at -0.5 m in the ocean, sensible heat transfer into the ocean and air, and latent heat flux of evaporation. Comparison with observational analyses shows that the diurnal variations in net radiation and heat storage in the mixed layer are realistically simulated. The simulated values of evaporation are consistent with other estimates for both ‘typical’ conditions and specific conditions during this four-day period. The rate of heat storage varies between +51 and -37 percent of the diurnal maximum incoming radiation, and the evaporation varies between +16% and -13% of this term. The non-dimensional transfer coefficients (C _D, C_T, C_q) computed from the model show general agreement with the coefficients calculated from observations in the simulated region (Pondet al., 1971). The simulated vertical profiles of temperature are in general agreement with observed profiles, except in the uppermost portions of the atmospheric boundary layer where deviations of approximately 1.5C occur. Simulated vertical profiles of wind speed are generally consistent with observed profiles, with the largest deviations appearing to be of the order of 0.5 m s^-1. Simulated vertical profiles of the eddy fluxes of sensible heat, water vapor, and momentum are generally consistent with Bunker's (1970) aircraft-based measurements of these quantities. The time averages of these simulated profiles show regular decreases with height, while simulated profiles for specific hours of the day show intermediate maxima and minima, which are also seen in the measured profiles. The vertically integrated kinetic energy budgets of the modelled atmospheric layer are presented through the four terms of the kinetic energy budget, viz., the upper and the lower boundary drags, dissipation, and potential-to-kinetic conversion. The dominant terms in the atmospheric energy budgets are the production and dissipation terms, with kinetic energy being exported both to the overlying atmospheric layer and to the underlying oceanic layer at rates of about 2 to 6% of the production, respectively. Comparisons between the climatological and BOMEX simulations are presented. The vertically integrated humidity budgets are presented for the two simulation experiments. Under ‘typical’ conditions, the humidity budget reveals an upper boundary flux of about +29% of the lower boundary flux with the vertically integrated advective flux being -59% of the lower flux. For the specific four-day simulation, the upper boundary flux and advection are about +28 and -70%, respectively, of the lower boundary flux.     

A Velocity–Dissipation Lagrangian Stochastic Model for Turbulent Dispersion in Atmospheric Boundary-Layer and Canopy Flows

An extended Lagrangian stochastic dispersion model that includes time variations of the turbulent kinetic energy dissipation rate is proposed. The instantaneous dissipation rate is described by a log-normal distribution to account for rare and intense bursts of dissipation occurring over short durations. This behaviour of the instantaneous dissipation rate is consistent with field measurements inside a pine forest and with published dissipation rate measurements in the atmospheric surface layer. The extended model is also shown to satisfy the well-mixed condition even for the highly inhomogeneous case of canopy flow. Application of this model to atmospheric boundary-layer and canopy flows reveals two types of motion that cannot be predicted by conventional dispersion models: a strong sweeping motion of particles towards the ground, and strong intermittent ejections of particles from the surface or canopy layer, which allows these particles to escape low-velocity regions to a high-velocity zone in the free air above. This ejective phenomenon increases the probability of marked fluid particles to reach far regions, creating a heavy tail in the mean concentration far from the scalar source.     

THE EJECTION-SWEEP CHARACTER OF SCALAR FLUXES IN THE UNSTABLE SURFACE LAYER

In the atmospheric surface layer, it is widely accepted that ejection andsweep eddy motions, typically associated with coherent structures, areresponsible for much of the land-surface evaporation, sensible heat, andmomentum fluxes. The present study analyzes the ejection-sweep propertiesusing velocity and scalar fluctuation measurements over tall natural grassand bare soil surfaces. It is shown that momentum ejections and sweeps occurat equal frequencies (D eject D sweep 0.29) irrespective of surfaceroughness length or atmospheric stability conditions. Also, their magnitudesare comparable to values reported from open channel velocity measurements (Dsweep 0.33; D eject : 0.30). The scalar D eject is constant andsimilar in magnitude to the momentum D eject( 0.29) over both surfacesand for a wide range of atmospheric stability conditions, in contrast to thescalar D sweep. The scalar sweep frequency is shown to depend on the scalarskewness for the dynamic convective and free convective sublayers, but isidentical to D eject for the dynamic sublayer. The threshold scalar skewnessat which the D sweep dependence occurs is 0.25, in agreement with theaccepted temperature skewness value at near-neutral conditions. In contrastto a previous surface-layer experiment, this investigation demonstrates thatthe third-order cumulant expansion method (CEM) reproduces the measuredrelative flux contribution of ejections and sweeps (S0) for momentumand scalars at both sites. Furthermore, a linkage between S0 and thescalar variance budget is derived via the third-order CEM in analogy tomomentum. It is shown that S0 can be related to the flux divergenceterm and that such a relationship can be estimated from surface-layersimilarity theory, and the three sublayer model of Kader and Yaglom andproposed similarity functions.     

Detection of long-term coherent exchange over spruce forest using wavelet analysis

Summary This study presents a discussion of a method for automated and quasi-online analysis of coherent structures using wavelet transform. The method is optimised for rapid processing of vector and scalar variables obtained over tall vegetation. It has been designed to assess long-term statistics of coherent structures, as it is applicable over a wide range of atmospheric conditions. Data of artificial and real turbulent signals are used to perform the analysis and to evaluate the presented method.Different wavelet functions are used for filtering the original signals, determining characteristic time scales, and detecting individual coherent structures. On this basis, statistics of temporal separation of coherent structures and phase shift between different variables can be calculated.Background turbulence and spikes are found to be efficiently removed without changing the shape, particularly the sharp localised gradients, of coherent structures. The determined peak in the calculated wavelet variance spectrum is observed to correspond very well to characteristic event durations and to satisfy the definition of coherent structures present in vector and scalar variables. The detection algorithm was successful in analysing data covering a wide range of atmospheric conditions. Detected individual coherent structures provide a parallel temporal pattern for scalar variables, but a phase shift between scalar and vector components.     

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.     

Boundary Conditions for Scalar (Co)Variances over Heterogeneous Surfaces

The problem of boundary conditions for the variances and covariances of scalar quantities (e.g., temperature and humidity) at the underlying surface is considered. If the surface is treated as horizontally homogeneous, Monin–Obukhov similarity suggests the Neumann boundary conditions that set the surface fluxes of scalar variances and covariances to zero. Over heterogeneous surfaces, these boundary conditions are not a viable choice since the spatial variability of various surface and soil characteristics, such as the ground fluxes of heat and moisture and the surface radiation balance, is not accounted for. Boundary conditions are developed that are consistent with the tile approach used to compute scalar (and momentum) fluxes over heterogeneous surfaces. To this end, the third-order transport terms (fluxes of variances) are examined analytically using a triple decomposition of fluctuating velocity and scalars into the grid-box mean, the fluctuation of tile-mean quantity about the grid-box mean, and the sub-tile fluctuation. The effect of the proposed boundary conditions on mixing in an archetypical stably-stratified boundary layer is illustrated with a single-column numerical experiment. The proposed boundary conditions should be applied in atmospheric models that utilize turbulence parametrization schemes with transport equations for scalar variances and covariances including the third-order turbulent transport (diffusion) terms.