Turbulent exchange above a pine forest II. Organized structures

Data from two 100-min runs on the turbulent atmospheric wind, temperature and humidity fields above a pine forest have been analysed using conditional sampling techniques. With the aid of the temperature time series, ramp events were identified and all fields were averaged in order to remove smaller scale turbulence and random low-frequency turbulence, and to map the organized structures revealed in this way. It is shown that the turbulent fluxes of momentum, heat and humidity, determined from these ramp-events, in fact constitute a large part of the total fluxes, about 90% during the actual events.On average, the duration of and the time between ramps were found to be about 35 s and 100 s, respectively. These periods were also found in the uw-, wT- and wq-cospectra, not only in the data analysed in this paper, but in the major part of all day-time data from this forest site (about 80 runs of 60 or 100 min). This indicates that the ramp events described arc a rather common phenomenon. The organized structures also have remarkable resemblance to near-wall structures observed in the laboratory.     

Turbulence measurements above a pine forest

Eddy fluxes of momentum, sensible and latent heat, and turbulence spectra measured over the Thetford Forest during 10 days in the Spring of 1973 are described. The measured total heat flux (H + E) for 122 20-min periods agreed closely on average with independent estimates from an energy balance method. There was evidence that the energy balance data gave small systematic overestimates of available energy during the hours before noon, compensated by slight underestimates for the remainder of the day. A comparison of measured wind speeds and friction velocities in neutral stability confirmed the validity of the aerodynamic method for estimating momentum fluxes at heights of a few roughness lengths above the canopy. In stable conditions the log-linear wind profileU = (u _*/k)(ln ((z -d)/z _o) + (z -d -z _o)/L) with = 3.4 ± 0.4 provided a good fit to the data. Spectra in unstable conditions were generally more sharply peaked than those measured by other workers over smoother terrain: differences were less marked in the case of vertical velocity in stable conditions. Temperature spectra in these stable conditions showed high energy at relatively low wavenumbers, andwT cospectra showed a cospectral gap; both of these results were associated with an intermittent sawtooth structure in the temperature fluctuations.Now at the Meteorological Office, Bracknell     

Fluxes of gases and particles above a deciduous forest in wintertime

Eddy-correlation measurements of the vertical fluxes of ozone, carbon dioxide, fine particles with diameter near 0.1 m, and particulate sulfur, as well as of momentum, heat and water vapor, have been taken above a tall leafless deciduous forest in wintertime. During the experimental period of one week, ozone deposition velocities varied from about 0.1 cm s^–1 at night to more than 0.4 cm s^-1 during the daytime, with the largest variations associated primarily with changes in solar irradiation. Most of the ozone removal took place in the upper canopy. Carbon dioxide fluxes were directed upward due to respiration and exhibited a strong dependence on air temperature and solar heating. The fluxes were approximately zero at air temperatures less than 5 °C and approached 0.8 mg m^–2 s^–1 when temperatures exceeded 15 °C during the daytime. Fine-particle deposition rates were large at times, with deposition velocities near 0.8 cm s^–1 when turbulence levels were high, but fluxes directed upward were found above the canopy when the surface beneath was covered with snow. Diffusional processes seemed to dominate fine-particle transfer across quasilaminar layers and subsequent deposition to the upper canopy. Deposition velocities for particulate sulfur were highly variable and averaged to a value small in magnitude as compared to similar measurements taken previously over a pine forest in summer.     

Monoterpene concentrations in and above a forest of scots pine

Diurnal and vertical ambient air measurements of the monoterpenes have been made in and above a Scots pine (Pinus sylvestris) forest of central Sweden, within the boreal northern coniferous biome. Sampling was done with Tenax TA, and analysis by GC and ion trap detection. Daytime mixing ratios were on the order of tenths of a ppbv from the forest floor to the top of the forest, and a factor of 2 or 3 lower above the forest. Mixing ratios at night were at the ppbv level, highest near the forest floor and the crown, and decreased with height above the forest. The highest total concentration observed was 8 ppbv inside the forest at 3 am (GMT). The average terpene composition was 3-carene 32%, -pinene 29%, limonene 18%, -pinene 10%, -phellandrene 7%, camphene 5%, and sabinene at less than 2%. The 3-carene/-pinene ratio varied with wind direction and speed, relative humidity, and wet/dry vegetation, but not with ozone or NO₂ concentration, solar radiation, or temperature. Variations in the observed terpene composition at the sampling site are mainly caused by the influence of other vegetation in the vicinity of the site. It would seem that wet Scots pine emits more 3-carene relative to -pinene than does dry pine.     

Turbulent kinetic energy budgets from a large-eddy simulation of airflow above and within a forest canopy

The output of a large-eddy simulation was used to study the terms ofthe turbulent kinetic energy (TKE) budget for the air layers above andwithin a forest. The computation created a three-dimensional,time-dependent simulation of the airflow, in which the lowest third ofthe domain was occupied by drag elements and heat sources to representthe forest. Shear production was a principal source of TKE in theupper canopy, diminishing gradually above tree-top height and moresharply with depth in the canopy. The transfer of energy to subgridscales (dissipation) was the main sink in the upper part of the domainbut diminished rapidly with depth in the canopy. Removal ofresolved-scale TKE due to canopy drag was extremely important,occurring primarily in the upper half of the forest where the foliagedensity was large. Turbulent transport showed a loss at the canopytop and a gain within the canopy. These general features have beenfound elsewhere but uncertainty remains concerning the effects ofpressure transport. In the present work, pressure was calculateddirectly, allowing us to compute the pressure diffusion term. Wellabove the canopy, pressure transport was smaller than, and opposite insign to, the turbulent transport term. Near the canopy top andbelow, pressure transport acted in concert with turbulent transport toexport TKE from the region immediately above and within the uppercrown, and to provide turbulent energy for the lower parts of theforest. In combination, the transport terms accounted for over half ofthe TKE loss near the canopy top, and in the lowest two-thirds of thecanopy the transport terms were the dominant source terms in thebudget. Moreover, the pressure transport was the largest source ofturbulent kinetic energy in the lowest levels of the canopy, beingparticularly strong under convective conditions. These resultsindicate that pressure transport is important in the plant canopyturbulent kinetic energy budget, especially in the lowest portion ofthe stand, where it acts as the major driving force for turbulentmotions.     

Analysis of long time series of long-wave radiation fluxes above a pine forest

Summary ?The long-term variations of upward terrestrial (E) and downward atmospheric (A) long-wave radiation fluxes above a pine stand in the southern part of the upper Rhine valley plain are analysed based on monthly mean values from 27 years of monitoring. Equivalent blackbody temperatures of the canopy and the atmosphere are calculated and compared to air temperatures at nearby sites. Based on 324 monthly values each, correlations between A and E as well as A and global radiation G are analysed. Only the former are highly correlated. Monthly mean values of long-wave radiation A can be calculated from air temperature, water vapour pressure and cloud cover. The long-term yearly average of equivalent blackbody temperature of the canopy is 0.6 K lower than kinetic air temperature at nearby sites. Only small, insignificant increases of both blackbody temperatures and air temperature are found. Despite the strong forest growth, it is surprising that the ratio of canopy emission temperature to air temperature did not change significantly. Relationships between the changes of general atmospheric circulation and equivalent blackbody temperatures of the canopy point to a strong dependence on shifts of general atmospheric flow. Received February 24, 2000; revised April 18, 2002; accepted July 20, 2002     

Turbulent transport within and above a maize canopy

Hot-wire anemometers were used to measure air temperature and the three velocity components of the wind within and above a maize canopy. From digitized anemometer outputs, correlation coefficients for vertical heat flux and turbulent momentum transfer were calculated. A comparison of these coefficients with profiles of mean wind speed and mean temperature indicates that the main features of the turbulence may be explained in terms of the usual mixing-length theory. Instantaneous records of heat and momentum flux, however, indicate the existence of other competing turbulent mechanisms due to the unsteady, non-equilibrium nature of the turbulent flow. Regimes of flow dominated by mechanical and/or thermal mixing are indicated. Spectral results show that high shear and turbulent intensity levels as well as the presence of the maize leaves and stalks as vortex-shedding surfaces complicate the energy transfer mechanism. An energy balance between radiation and convection reveals that the energy budget is primarily a balance between solar radiation and the flux of latent heat.Contribution of the Sibley School of Mechanical and Aerospace Engineering, Cornell University, in cooperation with the Agricultural Research Service, U.S. Department of Agriculture, Ithaca, N.Y., U.S.A. and the Cornell University Agricultural Experiment Station. Department of Agronomy Series No. 1116.Sibley School of Mechanical and Aerospace Engineering, Cornell University; U.S. Department of Agriculture, Gainesville, Florida Section for Estuary and Fjord Studies, River and Harbour Laboratory, Technical University of Norway, Trondheim, Norway; State Univ. of New York at Buffalo; and U.S. Department of Agriculture and Cornell University; respectively.     

Flux agreement above a Scots pine plantation

Summary The surface energy exchange of 12m high Scots pine plantation at Hartheim, Germany, was measured with a variety of methods during a 11-day period of fine weather in mid-May 1992. Net radiation and rate of thermal storage were measured with conventional net radiometers, soil heat flux discs and temperature-based storage models. The turbulent fluxes discussed in this report were obtained with an interchanging Bowen ratio energy budget system (BREB, at 14 m), two one-propeller eddy correlation systems (OPEC systems 1 and 2 at 17m), a 1-dimensional sonic eddy correlation system (SEC system 3) at 15 m, all on one low tower, and a 3-dimensional sonic eddy correlation system (SEC system 22) at 22 m on the high tower that was about 46 m distant. All systems measured sensible and latent heat (H and LE) directly, except for OPEC systems 1 and 2 which estimated LE as a residual term in the surface energy balance. Closure of turbulent fluxes from the two SEC systems was around 80% for daytime and 30% for night, with closure of 1-dimensional SEC system 3 exceeding that of 3-dimensional SEC system 22. The night measurements of turbulent fluxes contained considerable uncertainty, especially with the BREB system where measured gradients often yielded erroneous fluxes due to problems inherent in the method (i.e., computational instability as Bowen's ratio approaches –1). Also, both eddy correlation system designs (OPEC and SEC) appeared to underestimate |H| during stable conditions at night. In addition, both sonic systems (1- and 3-dimensional) underestimated |LE| during stable conditions. The underestimate of |H| at night generated residual estimates of OPEC LE containing a phantom dew error that erroneously decreased daily LE totals by about 10 percent. These special night problems are circumvented here by comparing results for daytime periods only, rather than for full days. To summarize, turbulent fluxes on the low tower from OPEC system 2 and the adjacent SEC system 3 were in reasonable agreement, while the BREB system appeared to overestimate H and underestimate LE; H and LE measured by SEC system 22 on the high tower were lower than from OPEC and SEC3 on the low tower. The turbulent flux measurements tended to converge, but the data exhibit unexplained differences between days, between systems, and between locations.With 7 Figures     

Turbulent characteristics over a rough natural surface part I: Turbulent structures

A comprehensive experiment for turbulence and profiles has been carried out over a very rough natural mallee bushland. In Part I, the large inherent structures of the turbulence are presented. Short-period space-time correlations of wind and temperature indicate that periodic fluctuations pervade the lower surface layer above the canopy. These are associated with the appearance of large turbulent structures. Prominent peaks in the spectra are identified with the arrival of these structures, which show that they have a relatively large length scale over rough surfaces. These findings suggest a different mechanism for the enhancement of turbulent diffusivity over rough surfaces from that by the simple generation of turbulence by plant wakes.Analysis via conditional sampling of the turbulent components has revealed that bursts make almost the same contributions as gusts (at the measurement height of 8.4 m above ground) for either momentum or temperature under neutral conditions; for unstable conditions and high hole size H (see Equation (3)), gusts dominate for momentum while bursts dominate for heat transfer. This implies that the mechanisms for momentum and heat transfer are different.     

Turbulent transport of scalar quantities within and above a paddy field

An experiment was conducted to study turbulent transport processes of scalar quantities within and above a rice plant canopy. A sonic anemometer-thermometer and a Lyman- humidiometer were used to measure the turbulent fluxes of sensible and latent heat and related turbulence statistics within a paddy field. The sensible and latent heat fluxes measured at two heights within and above the plant canopy showed that the upper layer of this plant canopy was an active source region and that the source strength of sensible and latent heat depended on the solar radiation and physiology of rice plants. Analysis of joint probability distributions of w and T and of w and q within this plant canopy showed that downdrafts were remarkably efficient for upward transport of sensible and latent heat in the daytime. The vertical fluxes of temperature and humidity variance were also divergent from the upper layer of plant canopies. The power spectra of temperature and humidity within the plant canopy decreased rapidly in the high frequency range, compared with the - 2/3 law relationship of nS(n) vs log n observed above plant canopies.     

Subfilter-scale Fluxes over a Surface Roughness Transition. Part I: Measured Fluxes and Energy Transfer Rates

A wind-tunnel experiment was designed and carried out to study the effect of a surface roughness transition on subfilter-scale (SFS) physics in a turbulent boundary layer. Specifically, subfilter-scale stresses are evaluated that require parameterizations and are key to improving the accuracy of large-eddy simulations of the atmospheric boundary layer. The surface transition considered in this study consists of a sharp change from a rough, wire-mesh covered surface to a smooth surface. The resulting magnitude jump in aerodynamic roughnesses, M = ln(z ₀₁/z ₀₂), where z ₀₁ and z ₀₂ are the upwind and downwind aerodynamic surface roughnesses respectively, is similar to that of past experimental studies in the atmospheric boundary layer. The two-dimensional velocity fields used in this study are measured using particle image velocimetry and are acquired at several positions downwind of the roughness transition as well as over a homogeneous smooth surface. Results show that the SFS stress, resolved strain rate and SFS transfer rate of resolved kinetic energy are dependent on the position within the boundary layer relative to the surface roughness transition. A mismatch is found in the downwind trend of the SFS stress and resolved strain rate with distance from the transition. This difference of behaviour may not be captured by some eddy-viscosity type models that parameterize the SFS stress tensor as proportional to the resolved strain rate tensor. These results can be used as a benchmark to test the ability of existing and new SFS models to capture the spatial variability SFS physics associated with surface roughness heterogeneities.     

Turbulent exchange over a surface with chessboard-type inhomogeneities

The results of an experiment carried out in the summer of 1988 in the forest-steppe region of the USSR are given. The aim of the experiment was to study turbulent exchange between the atmosphere and a non-homogeneous underlying surface consisting of a patchwork quilt of agricultural lands sown with different crops (hereafter called chessboard-type inhomogeneities). Simultaneous measurements of turbulent fluxes of heat, momentum and humidity by the eddy-correlation methods and also the measurements of the mean wind velocity and temperature are described. The values of turbulent fluxes calculated from the profile measurements are in good agreement with those obtained from the eddy-correlation measurements. Thus it may be concluded that the universal functions of Monin-Obukhov similarity theory (determined from measurements over a flat and homogeneous underlying surface) pertain also to a non-homogeneous (chessboard-type) underlying surface.     

Spectral Analysis of Turbulent Aerosol Fluxes by Fourier Transform,Wavelet Analysis,and Multiresolution Decomposition

Aerosol number fluxes are spectrally analyzed using fast Fourier transform analysis, wavelet analysis and multiresolution decomposition. All three methods yield similar spectral features in general, although a detailed evaluation of the cospectra shows some differences, e.g. due to different resolutions in the time and frequency domains. Wavelet analysis yields aerosol flux estimates with a high time resolution that can be used to assess the flux variability. Multiresolution decomposition has been applied successfully to evaluate cospectra of the aerosol number flux, the buoyancy flux and the momentum flux of three 1-day datasets from diverse environments. For all scalars and all environments, the dimensionless frequency (f) of the cospectral peak was found between $f = 0.1$ and 0.2. In addition, the cospectral gap time scale of the aerosol number flux was found between 100 and 1,000 s. Thus, in this study several spectral features such as the dominant time scale and the cospectral gap time scale of aerosol number fluxes are similar to buoyancy fluxes. However, the shape of aerosol number flux cospectra often deviates from buoyancy and momentum flux cospectra, especially at very small and at very large time scales.     

Resistance of a pine forest to ozone uptake

The surface resistance of a pine forest to ozone uptake is calculated from aircraft data obtained in southern New Jersey. The resistances are 130 and 300 s m^-1 for the afternoon and early evening, respectively, in good agreement with estimates used by Galbally and Roy (1980) in their calculation of global surface destruction of ozone.     

Large-eddy simulation of turbulent flow above and within a forest

A large-eddy simulation has been performed of an atmospheric surface layer in which the lower third of the domain is occupied by a drag layer and heat sources to represent a forest. Subgridscale processes are treated using second-order closure techniques. Lateral boundaries are periodic, while the upper boundary is a frictionless fixed lid. Mean vertical profiles of wind velocity derived from the output are realistic in their shape and response to forest density. Similarly, vertical profiles of Reynolds stress, turbulent kinetic energy and velocity skewness match observations, at least in a qualitative sense. The limited vertical extent of the domain and the artificial upper boundary, however, cause some departures from measured turbulence profiles in real forests. Instantaneous turbulent velocity and scalar fields are presented which show some of the features obtained by tower instrumentation in the field and in wind tunnels, such as the vertical coherence of vertical velocity and the slope of structures revealed by temperature patterns.     

Vertical divergence of fogwater fluxes above a spruce forest

Two almost identical eddy covariance measurement setups were used to measure the fogwater fluxes to a forest ecosystem in the “Fichtelgebirge” mountains (Waldstein research site, 786 m a.s.l.) in Germany. During the first experiment, an intercomparison was carried out with both setups running simultaneously at the same measuring height on a meteorological tower, 12.5 m above the forest canopy. The results confirmed a close agreement of the turbulent fluxes between the two setups, and allowed to intercalibrate liquid water content (LWC) and gravitational fluxes. During the second experiment, the setups were mounted at a height of 12.5 and 3 m above the canopy, respectively. For the 22 fog events, a persistent negative flux divergence was observed with a greater downward flux at the upper level. To extrapolate the turbulent liquid water fluxes measured at height z to the canopy of height h_c, a conversion factor 1/[1+0.116(z−h_c)] was determined. For the fluxes of nonvolatile ions, no such correction is necessary since the net evaporation of the fog droplets appears to be the primary cause of the vertical flux divergence. Although the net evaporation reduces the liquid water flux reaching the canopy, it is not expected to change the absolute amount of ions dissolved in fogwater.     

Turbulent dispersion in a stable layer with a quadratic exchange coefficient

Dispersion from a line source into a stable boundary layer of thickness l is analysed through solution of the diffusion equation assuming an exchange coefficient K(z) (1 – z/l) ² and wind profiles u(z) z ⁿ, with n = 0,1. Estimates of ground-level concentrations are made by developing analytic formulae where this is possible. A general method of solution using Laplace transformation and Green's function techniques is developed as an alternative to the eigenfunction expansion method discussed previously (Robson, 1983).