Experiments on scalar dispersion within a model plant canopy part I: The turbulence structure

This is the first of a series of three papers describing experiments on the dispersion of trace heat from elevated line and plane sources within a model plant canopy in a wind tunnel. Here we consider the wind field and turbulence structure. The model canopy consisted of bluff elements 60 mm high and 10 mm wide in a diamond array with frontal area index 0.23; streamwise and vertical velocity components were measured with a special three-hot-wire anemometer designed for optimum performance in flows of high turbulence intensity. We found that:

1. The momentum flux due to spatial correlations between time-averaged streamwise and vertical velocity components (the dispersive flux) was negligible, at heights near and above the top of the canopy.
2. In the turbulent energy budget, turbulent transport was a major loss (of about one-third of local production) near the top of the canopy, and was the principal gain mechanism lower down. Wake production was greater than shear production throughout the canopy. Pressure transport just above the canopy, inferred by difference, appeared to be a gain in approximate balance with the turbulent transport loss.
3. In the shear stress budget, wake production was negligible. The role of turbulent transport was equivalent to that in the turbulent energy budget, though smaller.
4. Velocity spectra above and within the canopy showed the dominance of large eddies occupying much of the boundary layer and moving downstream with a height-independent convection velocity. Within the canopy, much of the vertical but relatively little of the streamwise variance occurred at frequencies characteristic of wake turbulence.
5. Quadrant analysis of the shear stress showed only a slight excess of sweeps over ejections near the top of the canopy, in contrast with previous studies. This is a result of improved measurement techniques; it suggests some reappraisal of inferences previously drawn from quadrant analysis.

     

Multiple windbreaks: An aeolean ensemble

Near-neutral measurements of the turbulent wind field within and above a sequence of 15 parallel windbreaks on a flat pastoral site are presented. The windbreak fences each had a porosity of 60% and were equally-spaced at 6 times their height (h = 2 m). The following conclusions seem justified for wind directions within 10 ° of the normal to the array:

1. Above the windbreaks (2h), mean windspeeds first decreased and then increased asymptotically to a value in equilibrium with the new surface roughness. At 0.5h, windspeeds exhibited a slow increase down the entire array.
2. Reflecting differences in approach flows, the drag on the initial fence was almost twice that on barriers farther downstream. This reduction in momentum extraction per windbreak was associated with an elevation in the zero-plane displacement to a level equal to 0.8h.
3. At positions well-removed from the initial fences, mean windspeeds were reduced throughout the entire region below shelter height. In this region, the flow became increasingly dominated by downward moving air with velocities much greater than the local average. The zone of reduced turbulence was small, extending only 2h downstream of a barrier at a height of 0.25h. This corresponded with the region excluded from smoke trails released at the top of windbreaks.
4. An approximate TKE budget mid-way between windbreaks 7 and 8 suggests that shear and wake production peak near z = h and that production is balanced by dissipation and vertical transport components. Advective and inertial interaction terms are negligible at this midway position but are likely to be major sources of TKE closer to the windbreak. Local equilibrium is attained above z = 1.5h implying the existence of a constant-stress layer.

The measurements show the practical difficulty of simultaneously reducing both mean windspeeds and turbulence levels with repeated windbreaks at conventional spacings for horticultural applications.     

Atmospheric turbulence within and above a Douglas-fir stand. Part I: Statistical properties of the velocity field

This is the first of two papers reporting the results of a study of the turbulence regimes and exchange processes within and above an extensive Douglas-fir stand. The experiment was conducted on Vancouver Island during a two-week rainless period in July and August 1990. The experimental site was located on a 5^o slope. The stand, which was planted in 1962, and thinned and pruned uniformly in 1988, had a (projected) leaf area index of 5.4 and a heighth=16.7 m. Two eddy correlation units were operated in the daytime to measure the fluctuations in the three velocity components, air temperature and water vapour density, with one mounted permanently at a height of 23.0m (z/h=1.38) and the other at various heights in the stand with two to three 8-hour periods of measurement at each level. Humidity and radiation regimes both above and beneath the overstory and profiles of wind speed and air temperature were also measured. The most important findings are:

1. A marked secondary maximum in the wind speed profile occurred in the middle of the trunk space (aroundz/h=0.12). The turbulence intensities for the longitudinal and lateral velocity components increased with decreasing height, but the intensity for the vertical velocity component had a maximum atz/h=0.60 (middle of the canopy layer). Magnitudes of the higher order moments (skewness and kurtosis) for the three velocity components were higher in the canopy layer than in the trunk space and above the stand.
2. There was a 20% reduction in Reynolds stress fromz/h=1.00 to 1.38. Negative Reynolds stress or upward momentum flux perisistently occurred atz/h=0.12 and 0.42 (base of the canopy), and was correlated with negative wind speed gradients at the two heights. The longitudinal pressure gradient due to the land-sea/upslope-downslope circulations was believed to be the main factor responsible for the negative Reynolds stress.
3. Momentum transfer was highly intermittent. Sweep and ejection events dominated the transfer atz/h=0.60, 1.00 and 1.38, with sweeps playing the more important role of the two atz/h=0.60 and 1.00 and the less important role atz/h=1.38. But interaction events were of greater magnitude than sweep and ejection events atz/h=0.12 and 0.42.

     

Wind profiles at the Boulder Tower

Analysis of wind profiles at the Boulder Tower (BAO) leads to these conclusions:

1. The variation of roughness with wind direction found earlier is confirmed. Roughness lengths measured on the tower are larger than those measured close to the surface.
2. The profiles and measurements of Reynolds stress are consistent with a von-Karman constant of 0.35.
3. The form φ_m=(1?15z/L)^-1/3 fits best in the range -0.6 < z/L < 0. In the range 0 < z/L < 0.5, θ _m ～ 1 + 4.7z/L provides a good fit to the observations. For z/L < 0.1, φ _m also depends on h, the thickness of the PBL. For z/L < -0.6, Φ _m approaches the constant 0.5, in contrast to all previous suggestions. For larger stabilities, the upper level is usually not in the surface layer, and wind ratios become independent of z/L.
4. With snow cover, the effective roughness diminishes to about 1 cm, even for directions for which the roughness length without snow is large.
5. Estimation of winds at 100 or 150 m from information near the surface is best for similarity theory provided that the ratio of height to Monin-Obukhov L is less than 0.1. For larger z/L, simple power laws seem more appropriate.

     

The internal boundary layer — A review

A review is given of relevant work on the internal boundary layer (IBL) associated with:

1. Small-scale flow in neutral conditions across an abrupt change in surface roughness,
2. Small-scale flow in non-neutral conditions across an abrupt change in surface roughness, temperature or heat/moisture flux,
3. Mesoscale flow, with emphasis on flow across the coastline for both convective and stably stratified conditions.

The major theme in all cases is on the downstream, modified profile form (wind and temperature), and on the growth relations for IBL depth.     

A comparison of ABL heights inferred routinely from lidar and radiosondes at noontime

The height of the atmospheric boundary layer (ABL) obtained with lidar and radiosondes is compared for a data set of 43 noon (12.00 GMT) cases in 1984. The data were selected to represent the synoptic circulation types appropriately. Lidar vertical profiles at 1064 nm were used to obtain three estimates for the ABL height (h _lid), based on the first gradient in the back-scatter profile, namely, at the beginning, middle and top of the gradient. The boundary-layer height obtained with the radiosondes (h _s) was determined with the dry-parcel-intersection method in unstable conditions. As a first guess for near-neutral and stable conditions, the height of the first significant level in the potential temperature profile was taken. Overall, the boundary-layer thickness estimates agree surprisingly well (regression lineh _lidb=h_s:cc.=0.93 and the standard error=121 m). However, in 10% of the cases, the lidar estimate was significantly lower (difference>400 m) than the routinely inferredh _s. These outliers are discussed separately. For stable conditions, an estimate of ABL height (h _N) is also made based on the friction velocity and the Brunt-Väisälä frequency. The agreement betweenh _Nandh _lidbis good. Discrepancies between the two methods are caused by:

1. rapid growth of the boundary layer arround the measurement time;
2. the presence of a deep entrainment layer leading to a large zone in which quantities are not well mixed;
3. a large systematic error of 100–200 m in the estimate of boundary-layer height obtained from the radiosonde due to the way that profiles are recorded, as a series of significant points.

     

Atmospheric boundary layer research at Cabauw

At Cabauw, The Netherlands, a 213 m high mast specifically built for meteorological research has been operational since 1973. Its site, construction, instrumentation and observation programs are reviewed. Regarding analysis of the boundary layer at Cabauw, the following subjects are discussed:

- terrain roughness;

- Monin-Obukhov theory in practice;

- the structure of stable boundary layers;

- observed evolution of fog layers;

- inversion rise and early morning entrainment;

- use of the geostrophic wind as a predictor for wind profiles;

- height variation of wind climate statistics;

- air pollution applications: long range transport and short range dispersion;

- dependence of sound wave propagation on boundary-layer structure;

- testing of weather and climate models.

     

‘Application of a New Spectral Theory of Stably Stratified Turbulence to the Atmospheric Boundary Layer over Sea Ice’

A new spectral closure model of stably stratified turbulence is used to develop a K–ε model suitable for applications to the atmospheric boundary layer. This K–ε model utilizes vertical viscosity and diffusivity obtained from the spectral theory. In the ε equation, the Coriolis parameter-dependent formulation of the coefficient C ₁ suggested by Detering and Etling is generalized to include the dependence on the Brunt-Väisälä frequency, N. The new K–ε model is tested in simulations of the ABL over sea ice and compared with observations from BASE as simulated in large-eddy simulations by Kosovic and Curry, and observations from SHEBA.     

Bulk characteristics of heat transfer in the unstable,baroclinic atmospheric boundary layer

Field data for the unstable, baroclinic, atmospheric boundary layer over land and over the sea are considered in the context of a general similarity theory of vertical heat transfer. The dependence of δθ/θ_* upon logarithmic functions of h _c z _T and stability (through the similarity function C) is clearly demonstrated in the data. The combined data support the conventional formulation for the heat transfer coefficient δθ/θ_* when,

1. the surface scaling length is z _T (« z ₀), the height at which the surface temperature over land is obtained by extrapolation of the temperature profile
2. the height scale is taken as the depth of convective mixing h _c
3. the temperature profile equivalent of the von Karman constant is taken as 0.41
4. areal average, rather than single point, values of δθ are employed in strongly baroclinic conditions. No significant effect of baroclinity or the height scale ratio as proposed in the general theory is found. Variations in C about a linear regression relation against stability are most probably due to uncertainties in the areal surface temperature and to experimental errors in general temperature measurements.

     

Experiments on scalar dispersion within a model plant canopy part II: An elevated plane source

This paper describes a wind-tunnel experiment on the dispersion of trace heat from an effectively planar source within a model plant canopy, the source height being h _s = 0.80 h _c , where h _c is the canopy height. A sensor assembly consisting of three coplanar hot wires and one cold wire was used to make simultaneous measurements of the temperature and the streamwise and vertical velocity components. It was found that:

1. The thermal layer consisted of two parts with different length scales, an inner sublayer (scaling with h _s and h _c ) which quickly reached streamwise equilibrium downstream of the leading edge of the source, and an outer sublayer which was self-preserving with a length scale proportional to the depth of the thermal layer.
2. Below 2h _c , the vertical eddy diffusivity for heat from the plane source (K _HP ) was substantially less than the far-field limit of the corresponding diffusivity for heat from a lateral line source at the same height as the plane source. This shows that dispersion from plane or other distributed sources in canopies is influenced, near the canopy, by turbulence ‘memory’ and must be considered as a superposition of both near-field and far-field processes. Hence, one-dimensional models for scalar transport from distributed sources in canopies are wrong in principle, irrespective of the order of closure.
3. In the budgets for temperature variance, and for the vertical and streamwise components of the turbulent heat flux, turbulent transport was a major loss between h _s and h _c and a principal gain mechanism below h _s , as also observed in the budgets for turbulent energy and shear stress.
4. Quadrant analysis of the vertical heat flux showed that sweeps and ejections contributed about equal amounts to the heat flux between h _s and h _c , though among the more intense events, sweeps were dominant. Below h _s , almost all the heat was transported by sweeps.

     

On the desertification of the Sahel zone

Ground based measurements which were carried out in the Northern Sahel in southern Tunisia showed the following results:

1. The albedo difference between ground and protected land is about 10%, half of the amount Charney (1975) used in his model.
2. Bare soil is always warmer during times of bright sunshine than vegetated soil, which is in agreement with Jackson and Idso (1975). Temperature differences in excess of the 10 °C were observed between plants and the surrounding soil.
3. For bare soil, the surface temperature increases with declining albedo. However the opposite holds true for plants. Here, when lowering the albedo, a decrease in temperature was found.
4. In a sand dune field, the surface temperature depends strongly on the exposure. Surface temperature differences of 8 °C were observed for slopes of different exposures for measurements carried out around noon.

     

Turbulence decay in stratified and homogeneous marine layers

A spectral approach is applied to shear-induced turbulence in stratified layers. A system of spectral equations for stationary balance of turbulent energy and temperature variances was deduced in the vicinity of the local shear scale L_U = (ε/U_Z³)^1/2. At wavenumbers between the inertial-convective (k^−5/3) and wak turbulence (k⁻³) subranges, additional narrow spectral intervals—‘production’ subranges—may appear (E k⁻¹, E_T k⁻²). The upper boundary of these subranges is determined as L_U, and the lower boundaries as L_R (ε/U_ZN²)^1/2(χ/T_Z²). It is shown that the scale L_U is a unique spectral scale that is uniform up to a constant value for every hydrophysical field. It appears that the spectral scale L_U is equivalent to the Thorpe scale L_Th for the active turbulence model. Therefore, if turbulent patches are generated in a background of permanent mean shear, a linear relation between temperature and mass diffusivities exists. In spectral terms, the fossil turbulence model corresponds to the regime of the Boldgiano-Obukhov buoyancy subrange (E k^−11/5, E_T k^−7/5). During decay the buoyancy subrange is expanded to lower and higher wavenumbers. At lower wavenumbers the buoyancy subrange is bounded by L_** = 3(χ^1/2/N^1/2T_Z), which is equivalent to the Thorpe scale L_Th. In such a transition regime only, when the viscous dissipation rate is removed from the set of main turbulence parameters, the Thorpe scale does not correlate with the buoyancy scale L_N ε^1/2/N^3/2 and fossil turbulence is realized. Oceanic turbulence measurements in the equatorial Pacific near Baker Island confirm the main ideas of the active and fossil turbulence models.     

Intercomparison of turbulence data measured by SODAR and sonic anemometers

The capability of SODAR to measure the mean wind field in the lower boundary layer is well known and documented. Therefore, mean wind data are easily obtainable by means of the SODAR-technique, and are used to simulate the transport of pollutants after their release into the atmosphere. But when calculating the diffusion of pollutants, information about atmospheric turbulence is needed, too. In principle, a SODAR can measure turbulence data like the standard deviation of the vertical wind speed or horizontal wind direction. But when measuring turbulence data with a SODAR, one is beset by a host of limitations like volume sampling, spatial and temporal separation of sampling volume, attenuation of the acoustic waves and the slow speed of sound. Therefore, successful turbulence measurements with SODAR are not numerous and little is known about the quality of these data. In this context an intercomparison between a REMTECH-SODAR and a sonic anemometer mounted at the 100 m level of our meteorological tower was performed in summer 1990 at the Kernforschungszentrum Karlsruhe. The intercomparison is in two parts:

1. Half hour mean values of the standard deviation of the vertical wind speed are intercompared by scatter plots and by a linear regression and correlation analysis.
2. During 7 periods, 2 hours each, and covering atmospheric stabilities from unstable to slightly stable, the instantaneous vertical wind speeds were measured by both instruments and spectra were calculated.

The intercomparison demonstrates that DOPPLER-SODAR sounding is a reliable technique to determine besides the mean field, also athmospheric turbulence data like Sigma(w).     

Vertical cross-spectra of horizontal velocity components at the Boulder observatory

Cross-spectra between horizontal wind components at different levels of the Boulder Atmospheric Observatory (BAO) tower lead to the following conclusions:

1. Davenport's hypothesis is satisfied that coherence decays exponentially with the ratio of vertical separation to horizontal wave length, at least to very small values of coherence.
2. The decay coefficients increase with z/L for z/L < 0.5. For larger stabilities, irregular fluctuations with periods of order 10–20 min have considerable vertical coherence. Results at BAO are quite consistent with those elsewhere.
3. Eddy slopes in vertical planes increase with wind shear up to a point where the slope (horizontal delay over vertical separation) is just above 2. Beyond that point, the systematic increase of slopes with shear ceases. Since wind shear decreases upward, slopes tend to decrease upward. Slopes for lateral components are significantly larger than those for u-components.

     

Bare-ground surface heat and water exchanges under dry conditions: Observations and parameterization

A simplified land-surface parameterization is tested against bare-soil data collected during the EFEDA experiment conducted in Spain in June 1991. A complete data set, made up of soil properties as well as hydrological and atmospheric measurements, is described and discussed. The 11-day data set is characterized by very dry conditions and high surface temperatures during the day. Large values of sensible and soil heat fluxes and small values of surface evaporation (≈1 mm/day) were observed. This data set was modelled, leading to the following conclusions:

1. In the model, the parameterization provides values of the soil thermal properties and subsequently of the predicted soil heat fluxes which are overestimated when compared with the observations.
2. Following the literature, a value of the ratio between the roughness lengths for momentumZ _oand heatZ _ohof close to 10 for fairly homogeneous areas of bare soil and vegetation is used. This value leads to a fair prediction of the surface temperature. If the roughness lengths were taken to be equal, as is often assumed in atmospheric modelling, a poorer prediction results.
3. Finally, the vapor phase transfer mode is found dominant close to the surface and a modified parameterization including this effect is proposed. It allows a fair prediction of both surface evaporation and near-surface water content.

     

The scavenging of nitrate by clouds and precipitation

A model with spectral microphysics was developed to describe the scavenging of nitrate aerosol particles and HNO₃ gas. This model was incorporated into the dynamic framework of an entraining air parcel model with which we computed the uptake of nitrate by cloud drops whose size distribution changes with time because of condensation, collision-coalescence and break-up. Significant differences were found between the scavenging behavior of nitrate and our former results on the scavenging behavior of sulfate. These reflect the following chemical and microphysical differences between the two systems:

1. nitrate particles occur in a larger size range than sulfate particles.
2. HNO₃ has a much greater solubility than SO₂ and is taken up irreversibly inside the drops in contrast to SO₂.
3. nitric acid in the cloud water is formed directly on uptake of HNO₃ gas whereas on uptake of SO₂ sulfuric acid is formed only after the reaction with oxidizing agents such as e.g., H₂O₂ or O₃.
4. nitrate resulting from uptake of HNO₃ is confined mainly to small drops, whereas sulfate resulting from uptake of SO₂ is most concentrated in the largest, oldest drops, which have had the greatest time for reaction.

Sensitivity studies showed that the nitrate concentration of small drops is significantly affected by the mass accommodation coefficient.     

Numerical study on the bulk heat transfer coefficient for a variety of vegetation types and densities

The relationship between the geometrical structure of a canopy layer and the bulk transfer coefficient was investigated using a numerical canopy model. The following results were obtained:

1. The bulk transfer coefficients for momentum and heat, C _M and C _H , change with non-dimensional canopy density C _* each has a maximum.
2. The value of C _M is always larger than the value of C _H for a canopy with c _m > c _h , c _m and c _h being the drag coefficient and the heat transfer coefficient of an individual canopy element, respectively.
3. The value of C _* at which C _H has its maximum value is larger than the value of C _* at which C _M has its maximum. Therefore, the reciprocal of the sublayer Stanton number b _h ^?1 ranges between 50 and 65 for C _* around 0.1 while it ranges between 0 and 30 for C _* < 10^?2 and C _* > 2 (when c _m = 0.5).
4. The value of B _H ^?1 in the present study is consistent with most available observations, except for canopies of medium density (when C _* is around 0.1) for which no observational value has been obtained.

     

Experiments on scalar dispersion within a model plant canopy,part III: An elevated line source

An experiment is reported in which heat was released as a passive tracer from an elevated lateral line source within a model plant canopy, with h _s = 0.85 h _c (h _s and h _c being the source and canopy heights, respectively). A sensor assembly consisting of three coplanar hot wires and one cold wire was used to measure profiles of mean temperature % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikamaana% aabaGaeqiUdehaaiaacMcaaaa!390C!\[(\overline \theta )\], temperature variance (Σ_θ ²), vertical and streamwise turbulent heat fluxes, and third moments of wind and temperature fluctuations. Conclusions were:

1. Despite the very heterogeneous flow within the canopy, the observed dispersive heat flux (due to spatial correlation between time-averaged temperature and vertical velocity) was small. However, there is evidence from the plume centroid (which was lower than h _s at the source) of systematic recirculating motions within the canopy.
2. The ratio % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS% baaSqaaiabeI7aXjaab2gacaqGHbGaaeiEaaqabaGccaGGVaWaa0aa% aeaacqaH4oqCaaWaaSbaaSqaaiaab2gacaqGHbGaaeiEaaqabaaaaa!41DF!\[\sigma _{\theta {\text{max}}} /\overline \theta _{{\text{max}}} \] (of maximum values on vertical profiles) decreased from 1 near the source to an asymptotic value of 0.4 far downstream, in good agreement with previous experimental and theoretical work for concentration fluctuations in the surface layer well above the canopy.
3. The eddy diffusivity for heat from the line source (K _HL ) increased, downstream of the source, to a nearly constant ‘far-field’ vertical profile. Within the canopy, the far-field K _HL was an order of magnitude larger than K _HP , the equivalent diffusivity for a plane source; well above the canopy, the two were equal. The time scale defined by (far-field K _HL )/(vertical velocity variance) was independent of height within the canopy.
4. Budgets for temperature variance, vertical heat flux and streamwise heat flux are remarkably similar to the equivalent budgets for an elevated line source in the surface layer well above the canopy, except in the lower part of the canopy in the far field, where vertical transport is much more important than in the surface layer.
5. A random flight simulation of the mean height and depth of the temperature plume was generally in good agreement with experiment. However, details of the temperature and streamwise turbulent heat flux profiles were not correct, suggesting that the model formulation needs to be improved.

     

Outcomes and policy recommendations from the IPCC/AFOS working group on climate change response strategies and emission reductions

The most promising response strategies are

-sustainable practices in agriculture to improve productivity on existing arable land especially in developing countries to meet the food requirements of a still rising population

-sustainable practices in forestry both in tropical forests as well as in forests of temperate and boreal zones, in the latter case to achieve sufficient fast adaption to climate change.

     

Problems of atmospheric shear flows and their laboratory simulation

Shear flows generated by movement of the atmosphere near the earth's surface are accompanied by complexities not ordinarily encountered in the treatment of turbulent boundary layers. Problems arising from the following physical features are considered:

1. thermal stratification;
2. surface roughness in the form of forests and cities;
3. non-uniformity of surface roughness and/or temperature (leading to 3-dimensional turbulent boundary layers);
4. surface irregularities in the form of hilly and mountainous topography.

The complex nature of atmospheric shear flows has stimulated efforts to study their characteristics in the laboratory under controlled conditions. Accordingly, questions of similarity between the laboratory and the atmospheric flows for both mean and turbulent quantities arise. Similarity criteria, or appropriate scaling relationships, are discussed. Wind tunnels designed for investigations related to atmospheric shear flows are described. These facilities are shown to have a capability for simulating such flows for a wide range of the physical features listed above.