Footprint prediction of scalar fluxes from analytical solutions of the diffusion equation

The use of analytical solutions of the diffusion equation for footprint prediction is explored. Quantitative information about the footprint, i.e., the upwind area most likely to affect a downwind flux measurement at a given height z, is essential when flux measurements from different platforms, particularly airborne ones, are compared. Analytical predictions are evaluated against numerical Lagrangian trajectory simulations which are detailed in a companion paper (Leclerc and Thurtell, 1990). For neutral stability, the structurally simple solutions proposed by Gash (1986) are shown to be capable of satisfactory approximation to numerical simulations over a wide range of heights, zero displacements and roughness lengths. Until more sophisticated practical solutions become available, it is suggested that apparent limitations in the validity of some assumptions underlying the Gash solutions for the case of very large surface roughness (forests) and tentative application of the solutions to cases of small thermal instability be dealt with by semi-empirical adjustment of the ratio of horizontal wind to friction velocity. An upper limit of validity of these solutions for z has yet to be established.     

Diurnal variation of horizontal wind direction fluctuations in complex terrain at Geysers,Cal.

Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km² Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s^-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.     

Free convection scaling over the ocean derived from buoy measurements during gate

The structure of atmospheric turbulence in the surface layer over the open ocean is examined under conditions of local free convection. The raw data consist of profile and fluctuation measurements of wind and temperature as obtained from a meteorological buoy. For near neutral conditions and for waves running approximately along the wind direction, wave-induced wind fluctuations can be described by a simplified linear theory based on Miles (1957). In this case, the spectrum of wind velocity is given as the sum of two parts; for the turbulent part, the parameterization as obtained by Kaimal et al. (1972) applies, while the wave-induced part is parameterized using a simplification of Miles' linear theory. For cases of local free convection, the measurements of the vertical component of the wind velocity are well described by similarity theory; as expected, _w /(-uw)^1/2 is proportional to (- z/L)^1/3. In order to scale the longitudinal wind velocity component, it seems to be reasonable to extend the list of relevant parameters by the height of the mixed layer z _i. We obtain _u /(- uw)^1/2 (z/z _i)^1/3(- z/L)^1/3 with only a poor correlation coefficient of r = 0.6. Overall, the results of local free convection scaling obtained from direct measurements show good agreement with those obtained from profile measurements. A comparison between direct and indirect determination of turbulent fluxes of momentum shows an unexplained difference of about 20%. This discrepancy is mainly due to a gap in the uw-cospectrum at the swell frequency.     

Wind profile development above a locally adjusted sea surface

Results are presented from a numerical experiment of wind and shear stress profile development away from a shore line; the water surface is assumed to obey the Charnock-Ellison relation between surface roughness and friction velocity. In typical cases the upwind, land surface is rough relative to the sea and the velocity and shear stress results are qualitatively similar to those for flows from relatively rough to relatively smooth solid surfaces. In the present case, however, the downwind surface roughness and friction velocity vary with position and we find that wind profile development may play a significant role in the relationship between sea surface roughness and fetch.     

One-dimensional checkerboards and blending heights

We analyze the checkerboard problem of many alternating surfaces with different properties, on scales up to (say) 3,000 m. Power-law representations of the vertical profiles of mean wind speed and eddy diffusivity lead to solutions in terms of Kelvin and trigonometric functions.These solutions are used to determine blending heights (_*), where deviations from the mean of concentration, or of vertical flux density, fall to some small fraction, , of their value at the surface. Values of _*are important for regional and larger-scale meteorological models. In smaller scale micrometeorological studies, they may serve also as the top levels of surface boundary layers.An important result for both theoretical and experimental contexts is that deviations of flux persist with elevation much more strongly than those of concentration, so that, in general, _* should be based on flux rather than concentration. Representative values of _*, for = 0.05, are of order 5 and 30 m for surface pattern wavelengths of 10² and 10³ m, respectively. Values of _* are robust to changes in adopted power-law indices, and are independent of wind speed. Surface roughness has a mild but calculable effect.     

Estimates of diabatic wind speed profiles from near-surface weather observations

In this paper we analyse diabatic wind profiles observed at the 213 m meteorological tower at Cabauw, the Netherlands. It is shown that the wind speed profiles agree with the well-known similarity functions of the atmospheric surface layer, when we substitute an effective roughness length. For very unstable conditions, the agreement is good up to at least 200 m or z/L–7(z is height, L is Obukhov length scale). For stable conditions, the agreement is good up to z/L1. For stronger stability, a semi-empirical extension is given of the log-linear profile, which gives acceptable estimates up to ~ 100 m. A scheme is used for the derivation of the Obukhov length scale from single wind speed, total cloud cover and air temperature. With the latter scheme and the similarity functions, wind speed profiles can be estimated from near-surface weather data only. The results for wind speed depend on height and stability. Up to 80 m, the rms difference with observations is on average 1.1 m s^–1. At 200 m, 0.8 m s^–1 for very unstable conditions increasing to 2.1 m s^–1 for very stable conditions. The proposed methods simulate the diurnal variation of the 80 m wind speed very well. Also the simulated frequency distribution of the 80 m wind speed agrees well with the observed one. It is concluded that the proposed methods are applicable up to at least 100 m in generally level terrain.     

Lateral coherence of longitudinal wind components in strong winds

A combination of lateral coherence measurements of wind speed at five locations suggests that the decay constant is a monotonically increasing function of the ratio of separation to height, under neutral conditions.     

Wind characteristics in the lowest 340 m of the atmospheric boundary layer at Pune,India

Wind characteristics in the lowest 340 m (agl) of the atmospheric boundary layer at Pune (18°32N, 73°51E, 559 m ASL) have been investigated using the pilot balloon wind observations obtained during the summer monsoon seasons of 1976, 1979 and 1980. Variations in the zonal and meridional components of wind at the surface, 40, 150 and 340 m (agl) have been described. Frequency distribution of the zonal component showed an unimodal character at the surface, which gradually approached a multimodal character at 340 m. The distribution pattern of the wind components was, by and large, normal. Spectral analysis of the wind components showed that the spectral energy was predominantly shared by 5–7 and 2–3 day periods.     

Boundary-layer parameterizations for applied dispersion modeling over urban areas

Many applied dispersion models require the knowledge of boundary-layer parameters such as sensible heat flux,Q _H , friction velocity,u _*, and turbulent energy components, _w and _v . Formulas are suggested for calculating these parameters over a wide variety of types of ground surfaces, based on simple observations of wind speed near the ground and fractional cloud cover, and specification of constants such as roughness length, albedo, and soil moisture availability. Observations ofu _*,Q _H , _w , and _v during field experiments in St. Louis and Indianapolis are used to test the formulas for urban sites. Relative errors of about ±20% in the predictions are seen to occur whenu _*,Q _H , _w , and _v are large. However, when these quantities are small (e.g.,u _* < 0.2 m/s), the errors in the predictions are as large as the mean value of the quantity itself.In addition, it is concluded from studies of available field data and theories that the magnitude of _w is not well-known at elevations above about 100m during the late afternoon and night. Some simple parameterizations for _w . are suggested that are consistent with the observed steady decrease in ground-level concentration in the afternoon and the sudden increase in concentration that can occur a few hours after sunset due to wind shears associated with a low-level jet, for continuous plumes emitted from moderate to tall stacks.     

Flux Footprints Within and Over Forest Canopies

The characteristics of turbulence within a forest arespatially heterogeneous and distinct from thoseassociated with the surface boundary layer. Consequently, the size and probability distribution offlux footprints emanating from sources below aforest canopy have the potential to differ from thoseobserved above forests.A Lagrangian random walk model was used toinvestigate this problem since no analytical solutionof the diffusion equation exists. Model calculations suggest that spatialcharacteristics of flux footprint probabilitydistributions under forest canopies are muchcontracted, compared to those evaluated in the surfaceboundary layer. The key factors affecting thestatistical spread of the flux footprint, and theposition of the peak of its probability distribution,are horizontal wind velocity and the standarddeviations of vertical and horizontal velocityfluctuations. Consequently, canopies, which attenuatemean horizontal wind speed, or atmospheric conditions,which enhance vertical velocity fluctuations, willcontract flux footprint distributions mostly near thefloor of a forest. It was also found that theprobability distributions of the flux footprint arenarrower when horizontal wind velocity fluctuationsare considered, instead of the simpler case that considers only vertical velocity fluctuations and meanhorizontal wind velocity.     

On characteristics of wind direction fluctuations in the surface layer

It is shown that the ratio of standard deviation of lateral velocity to the friction velocity, /u _*, and therefore wind direction fluctuations, are sensitive to mesoscale terrain properties. Under neutral conditions, /u _* is almost 40% larger in rolling terrain than over a horizontal surface. In the lee of a low mountain, the fluctuations may be 2.5 times as strong as over horizontal terrain. In contrast, vertical velocity fluctuations are little influenced by mesoscale terrain features.Now with Air Weather Service, Offutt AFB, Omaha, Nebraska.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

The Value Of The β Coefficient In The Relaxed Eddy Accumulation Method In Terms Of Fourth-Order Moments

The commonly measured value of in the relaxed eddy accumulationmethod of about 0.56is shown to arise from the non-Gaussiannature of turbulence. Fourth-orderGram–Charlier functions forthe two-dimensional probability distributionsof variation in the horizontal component of wind velocityand concentrations of water vapour, carbondioxide and methane with respect to thevertical component of wind velocity are used to examinethe value of .An analytical solution for ispresented in terms of fourth-order moments.Under mean conditions, this solution givesa value for of0.557. Variation of is shown to be controlledprimarily by the ratio of the mean ofc'w3 (where c'is relevant to the entity of interest andw' is vertical component of windvelocity) to the correlationcoefficient between the entity concentrationand vertical component of wind velocity.     

On the memory of wind standard deviation for upstream roughness

Big eddies in the outer part of the atmospheric boundary layer contribute to the variance of the horizontal velocity fluctuations near the surface. Because of the slow adjustment of these eddies to new boundary conditions, they carry the roughness characteristics of a large upstream terrain. A scaling relation is proposed that accounts for the memory effects in the big eddies. It is concluded that the standard deviation of the horizontal wind ( _u ) measured at a given height is representative for the shear stress at greater height. This gives at least qualitative support to existing work where u is used for exposure correction of mean wind.     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Surface drag and turbulence over an inhomogeneous land surface

Data collected over an inhomogeneous semi-rural area are presented. The data are compared with previous surface-layer data to determine how representative the fixed-point flux measurements are of area averages. Departures from the standard surface-layer results are found to be relatively small (~10–20%), which supports the concept of a blending height above which the flow ceases to respond to variations in the underlying surface and becomes horizontally homogeneous.Effective roughness lengths are derived for different wind directions and the relationship between the effective roughness length and upwind surface is examined in the light of recent ideas on averaging surface roughness lengths. It is found that by averaging drag coefficients, realistic values of the effective roughness length can be calculated which are not very sensitive to the precise choice of the component roughness lengths.     

Averaging Intervals For Spectral Analysis Of Nonstationary Turbulence

We formulate a method for determining the smallest time interval Tover which a turbulence time series can be averaged to decompose it intoinstantaneous mean and random components. From the random part the method defines the optimal interval (or averaging window) AW over which this part should be averaged to obtain the instantaneous spectrum. Both T and AW vary randomly with time and depend on physical properties of the turbulence. T also depends on the accuracy of the measurements and is thus independent of AW. Interesting features of the method are its real-time capability and the non-equality between AW and T.     

Two-point velocity statistics over Lake Ontario

Cross-spectra between wind speeds on several masts in Lake Ontario have been analyzed. As previously predicted, coherence over water (small intensity of turbulence) between wind speeds measured on masts lined up with the wind appears to be larger than over land, and increases with decreasing Richardson numbers. As a result, in cold air over warm water, wind speed fluctuations are well predictable from upstream measurements. For large angles between the anemometer line and wind, the difference between coherence over land and water disappears. Furthermore, there is no significant difference in vertical coherence between water and land. When the wind is parallel to the anemometer line, small eddies travel, in agreement with Taylor's hypothesis, with the local mean wind speed. Large eddies travel significantly faster. Vertical phase delay increases with increasing Richardson number.     

Spatial Variability of Turbulent Fluxes in the Roughness Sublayer of an Even-Aged Pine Forest

The spatial variability of turbulent flow statistics in the roughness sublayer (RSL) of a uniform even-aged 14 m (= h) tall loblolly pine forest was investigated experimentally. Using seven existing walkup towers at this stand, high frequency velocity, temperature, water vapour and carbon dioxide concentrations were measured at 15.5 m above the ground surface from October 6 to 10 in 1997. These seven towers were separated by at least 100m from each other. The objective of this study was to examine whether single tower turbulence statistics measurements represent the flow properties of RSL turbulence above a uniform even-aged managed loblolly pine forest as a best-case scenario for natural forested ecosystems. From the intensive space-time series measurements, it was demonstrated that standard deviations of longitudinal and vertical velocities (_u, _w) and temperature (_T) are more planar homogeneous than their vertical flux of momentum (u_* ²) and sensible heat (H) counterparts. Also, the measured H is more horizontally homogeneous when compared to fluxes of other scalar entities such as CO₂ and water vapour. While the spatial variability in fluxes was significant (>15 %), this unique data set confirmed that single tower measurements represent the canonical structure of single-point RSL turbulence statistics, especially flux-variance relationships. Implications to extending the moving-equilibrium hypothesis for RSL flows are discussed. The spatial variability in all RSL flow variables was not constant in time and varied strongly with spatially averaged friction velocity u_*, especially when u_* was small. It is shown that flow properties derived from two-point temporal statistics such as correlation functions are more sensitive to local variability in leaf area density when compared to single point flow statistics. Specifically, that the local relationship between the reciprocal of the vertical velocity integral time scale (I_w) and the arrival frequency of organized structures (/h) predicted from a mixing-layer theory exhibited dependence on the local leaf area index. The broader implications of these findings to the measurement and modelling of RSL flows are also discussed.