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.     

Joint probability analysis of momentum and heat fluxes at a deciduous forest

Turbulent transport processes for momentum and scalar quantities are examined by a joint probability distribution analysis using data observed within and above a deciduous forest. Characteristics of transport processes in the frequency domain were also analyzed using Tukey's procedure. The results confirm that sweep phenomena prevail within and at the top of a tall plant canopy and that downdrafts are more effective for vertical transport of momentum and scalar quantities. On the other hand, updrafts become more efficient for vertical transport in the daytime at levels about twice treetop height. The results show that within the forest, the sweep phenomenon prevails over a wide frequency range, while above the forest, prevalence of the ejection phenomenon is limited to low frequencies. It is again noted that the plant canopy plays an important role in the sweep-ejection cycle as well as in turbulent transport processes.     

On the downward transport of turbulent kinetic energy in the surface layer over plant canopies

Downward fluxes of turbulent kinetic energy have been frequently observed in the air layer just above plant canopies. In order to investigate the mechanism for such downward transport, analysis of observational data is attempted. Height-dependency of turbulent kinetic energy flux and turbulence statistics including higher order moments is represented as a function of a non-dimensional height z/H, where z is an observational height and H an average height of plant canopies. Downward fluxes and non-Gaussianity of wind velocity fluctuations are predominant just above plant canopies and decrease with increasing height. The downward flux is closely related to the high intensity of turbulence and the non-Gaussianity of wind velocity fluctuations, especially with a positive skewness in the longitudinal wind and a negative skewness in the vertical wind. The analysis method of conditional sampling and averaging is applied to the present observations. The results show that the predominance of the intermittent inrush phase over the intermittent ejection phase leads to the above-mentioned non-Gaussianity. Finally, a simple explanation is given in order to interpret the turbulent flow structure in the air layer near the plant canopies, which is associated with the downward energy transport process.     

Turbulence measurements during mistral winds with a 1-dimensional sonic anemometer

Fluctuations in the vertical wind velocity and air temperature were measured with a 1-dimensional sonic anemometer and fine thermocouple over a flat agricultural site in the Rhone Valley, France. Strong Mistral winds with speeds up to 20 m s^–1 kept atmospheric conditions very close to neutral and ensured stationarity. Friction velocities estimated both by eddy correlation (sonic plus Gill Bivane) and inertialdissipation (sonic only) methods agreed within 1 and 5 % respectively of traditional profile measurements over the measured range of 0.2 to 1.2 m s^–1. The coefficient of eddy transport for heat exceeded that of momentum by a factor of 1.38 (± 0.05), a result almost identical to that obtained in the Kansas experiment (Businger et al., 1971). For - 0.15 >= z/L >= 0.05, the ratio _w /u _* was 1.69 and 1.34 for unstable and stable conditions, respectively. For ¦z/L¦ >= 0.05, the ratio /T _* was 1.40 independent of whether neutrality was approached from either stable or unstable conditions.     

Observed generation of an atmospheric gravity wave by shear instability in the mean flow of the planetary boundary layer

Observations of a single boundary-layer event — the generation of an atmospheric gravity wave by an unstable shear flow at Haswell, Colorado on November 12, 1971 — are briefly described and discussed. The observations were made using: (a) an acoustic echo sounder, (b) anemometers mounted at two fixed levels on a 150-m tower, (c) an anemometer and a thermometer mounted on a movable carriage on the tower, and (d) a microbarograph array, including one microbarograph mounted atop the tower. The wave phase velocity (–3.5–4.0 m s^–1) was found to equal the wind velocity in the middle of the shear flow, as assumed by other authors. The wave-associated vertical fluxes of momentum and energy measured just above the wave critical layer were estimated to be –5 dyn cm^–2 and –800 erg cm^–2 s^–1, respectively. These are large values. The annual average vertical flux of momentum at temperate and high latitudes is –0.25 dyn cm^–2, while the average kinetic energy dissipation rate in a unit column of atmosphere is –5 × 10³ erg cm^–2 s^–1. If the region of wave generation was itself propagating horizontally, its propagation velocity was large compared with the horizontal phase speed of the small-scale waves generated. Wave generation appeared to occur over an area large compared with the size of the microbarograph array (i.e., 2 km).     

A comparison of turbulence statistics in the surface layer over plant canopies with those over several other surfaces

Turbulence statistics, including higher order moments, in the surface layer over plant canopies were compared with those observed over several different surfaces, using a nondimensional height (z – d)/z ₀: The values of (z – d)/z ₀extend over a very wide range from 10 over plant canopies to 10⁷ over the ocean. Several properties such as intensities of turbulence and skewness factors show a remarkable height-dependency in the air layer below (z – d)/z ₀ = 10², which is supposed to be much influenced by the underlying surface. In that layer, some peculiar phenomena, such as a downward energy transport and positive flux of shear stress, are frequently observed.     

Vertical transport of turbulent kinetic energy in the surface layer over a paddy field

Turbulent kinetic energy and its vertical flux were measured at two heights over a paddy field. The vertical transport of turbulent kinetic energy was always downward right above the paddy field and was frequently downward at higher levels within a few metres above the crop. Contributions to the downward transport arise mainly from the turbulent kinetic energy of horizontal wind velocity components. It is shown from the analysis of probability distributions that appreciable transport takes place intermittently in a few large downward bursts and that these downdrafts are efficient for downward energy transport.In the budget of turbulent kinetic energy, the flux divergence term and the energy dissipation term are the main loss terms under unstable conditions. These terms increase in magnitude with increasing instability. Buoyant production is insufficient to balance these losses. The imbalance term involving the pressure-work term is probably one of the main energy sources in unstable conditions.     

Detection and analysis of coherent structures in urban turbulence

Summary The continuous wavelet transform provides a suitable tool to visualize the vertical structure of turbulence and to detect coherent structures in turbulent time series. This is demonstrated with a simple example of an artificially ramp structured time series. In this study turbulence data, i.e. the fluctuations of the horizontal wind components u′ and v′, the vertical component w′ and temperature T′, sampled with 20.83 Hz and measured simultaneously at three levels (z/h=1.5, 2.1 and 3.2, with z as the sensor height and h the height of the roughness elements) over an urban canopy in the inner city of Basel, Switzerland, are analyzed. The detection of the coherent structures was performed using the Mexican hat wavelet and the zero-crossing method. The analysis for unstable conditions shows that organized structures (ejection-sweep cycles) cover about 45% of the total run time. A conditional average from a total of 116 detected ejection-sweep sequences during 7 hours was calculated over a time window of 100 seconds. This dominating time scale was derived from peak frequencies of the wavelet spectra as well as from the Fourier spectra. It is shown that the normalized amplitudes of fluctuations of temperature and longitudinal wind speed during the events are largest at the lowest measurement level just above the canopy and decrease with increasing distance from the roughness elements. A comparison of related studies over different non-urban surfaces (mainly forests) shows that the shape of conditionally averaged ejection-sweep sequences is very similar for all canopies, however, the dominating time scale in general increases the rougher the surface is and the higher the roughness elements are.     

Tilt errors and precipitation effects in trivane measurements of turbulent fluxes over open water

For 390 ten-minute samples of turbulent flux, made with a trivane above a lake, the vertical alignment is determined within 0.1 ° through azimuth-dependent averaging. One degree of instrumental misalignment is found to produce an average tilt error of 9 ± 4% for momentum flux, and 4 ± 2% for heat flux. The tilt error in the vertical momentum flux depends mainly ons _u/u_*, and cannot be much diminished with impunity by high-pass pre-filtering of the turbulence signals. The effects of rain on trivane measurements of vertical velocity are shown to be negligible at high wind speeds, and adaptable to correction in any case.The normalized vertical velocity variance,s _w/u_*, appears to be proportional to the square root ofz/L for unstable stratification. For a wind speed range of 2 to 15 m s^–1, the eddy correlation stresses measured at 4- and 8-m heights can be reasonably well estimated by using a constant drag coefficientC _d=1.3 X 10^-3, while cup anemometer profile measurements give an overestimate of eddy stress at high wind speeds. A good stress estimate is also obtained from the elevation variance; it is suggested that trivane measurement of this variance might be made from a mobile platform, e.g., a moderately stabilized spar buoy.     

Seasonal and diurnal variations of coherent structures over a deciduous forest

Coherent structures in turbulent flow above a midlatitude deciduous forest are identified using a wavelet analysis technique. Coupling between motions above the canopy (z/h=1.5, whereh is canopy height) and within the canopy (z/h=0.6) are studied using composite velocity and temperature fields constructed from 85 hours of data. Data are classified into winter and summer cases, for both convective and stable conditions. Vertical velocity fluctuations are in phase at both observation levels. Horizontal motions associated with the structures within the canopy lead those above the canopy, and linear analysis indicates that the horizontal motions deep in the canopy should lead the vertical motions by 90°. On average, coherent structures are responsible for only about 40% of overall turbulent heat and momentum fluxes, much less than previously reported. However, our large data set reveals that this flux fraction comes from a wide distribution that includes much higher fractions in its upper extremes. The separation distanceL _s between adjacent coherent structures, 6–10h, is comparable to that obtained in previous observations over short canopies and in the laboratory. Changes in separation between the summer and winter (leafless) conditions are consistent withL _s being determined by a local horizontal wind shear scale.     

Structure Function Analysis of Two-Scale Scalar Ramps. Part II: Ramp Characteristics and Surface Renewal Flux Estimation

Ramp features in the turbulent scalar field are associated with turbulent coherent structures, which dominate energy and mass fluxes in the atmospheric surface layer. Although finer scale ramp-like shapes embedded within larger scale ramp-like shapes can readily be perceived in turbulent scalar traces, their presence has largely been overlooked in the literature. We demonstrate the signature of more than one ramp scale in structure functions of the turbulent scalar field measured from above bare ground and two types of short plant canopies, using structure-function time lags ranging in scale from isotropic to larger than the characteristic coherent structures. Spectral analysis of structure functions was used to characterize different scales of turbulent structures. By expanding structure function analysis to include two ramp scales, we characterized the intermittency, duration, and surface renewal flux contribution of the smallest (i.e., Scale One) and the dominant (i.e., Scale Two) coherent structure scales. The frequencies of the coherent structure scales increase with mean wind shear, implying that both Scale One and Scale Two are shear-driven. The embedded Scale One turbulent structure scale is ineffectual in the surface-layer energy and mass transport process. The new method reported here for obtaining surface renewal-based scalar exchange works well over bare ground and short canopies under unstable conditions, effectively eliminating the α calibration for these conditions and forming the foundation for analysis over taller and more complex surfaces.     

The Influence of Local Stability on Heat and Momentum Transfer within Open Canopies

Eddy-covariance data have been analyzed to investigate the influence of local stability on heat transfer within open canopies. The flux–gradient relationship for heat is derived from the temperature variance equation, and the stability dependence of the flux–gradient relationship is examined and discussed. The results indicate that the strong stability dependence of the nondimensional standard deviation of temperature, and the small contributions of turbulent transport to the temperature variance, lead to a strong stability dependence of the nondimensional temperature gradient within open canopies. Quadrant analysis and hole size analysis were performed for momentum and heat fluxes in the subcanopy, and the results indicate that the contribution of each quadrant to the total flux depends on both the local stability and canopy depth. The intermittency of the turbulent flux does not show a clear dependence on local stability. As the contribution of ejections to the heat flux increases, the vertical flux of the temperature variance changes sign from negative to positive, leading to small temperature variance transport in unstable conditions. Multi-resolution analysis indicates that heat and momentum are transported with different dominant time scales in very unstable conditions, suggesting a different role of local buoyancy in heat and momentum transfer.     

NUMERICAL SIMULATIONS OF TURBULENT FLOW IN AND ABOVE PLANT CANOPIES

The turbulent flow in and above plant canopies is of fundamental importance to the understanding oftransport processes of momentum,heat and mass between plant canopies and atmosphere,and to microme-teorology.The Reynolds stress equation model(RSM)has been applied to calculate the turbulence in cano-pies in this paper.The calculated mean wind velocity profiles,Reynolds stress,turbulent kinetic energy andviscous dissipation rate in a corn canopy and a spruce forest are compared with field observed data and withWilson's and Shaw's model.The velocity profiles and Rynolds stress calculated by both models are in goodagreement,and the length scale of turbulence appears to be similar.     

Daytime variations of ozone eddy fluxes to maize

The vertical fluxes of ozone, momentum and heat in the atmospheric surface layer have been measured by eddy correlation above both mature and senescent maize canopies. Aerodynamic formulae are applied to find that the bulk canopy surface resistancer _c to ozone uptake and destruction varies between 4.0 and 0.5 s cm⁻¹ during the daytime. Apparently, surface properties tend to control the removal of ozone at the surface of the earth. For a lush canopy, the stomatal diffusion resistance is the most important property, while changes in surface temperature have little effect. Destruction at the soil and exterior plant surfaces appears to account for 20–50% of the total loss if leaf mesophyll resistances are assumed to be very small. Free water at leaf surfaces may at times inhibit ozone removal by both senescent and healthy plants.     

Drag coefficients and turbulence spectra within three boreal forest canopies

Atmospheric turbulence was measured within a black spruce forest, a jack pine forest, and a trembling aspen forest, located in southeastern Manitoba, Canada. Drag coefficients (C _d ) varied little with height within the pine and aspen canopies, but showed some height dependence within the dense spruce canopy. A constant C _d of 0.15, with the measured momentum flux and velocity profiles, gave good estimates of leaf-area-index (LAI) profiles for the pine and aspen canopies, but underestimated LAI for the spruce canopy.Velocity spectra were scaled using the Eulerian integral time scales and showed a substantial inertial subrange above the canopies. In the bottom part of the canopies, the streamwise and cross-stream spectra showed rapid energy loss whereas the vertical spectra showed an apparent energy gain, in the region where the inertial subrange is expected. The temperature spectra showed an inertial subrange with the expected -2/3 slope at all heights. Cospectra of momentum and heat flux had slopes of about -1 in much of the inertial subrange. Possible mechanisms to explain some of the spectral features are discussed.     

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

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

Measurements Of Fine-Scale Structure At The Top Of Marine Stratocumulus

During the Dynamics and Chemistry of the MarineStratocumulus (DYCOMS) experiment in July–August 1985, the NCAR Electra aircraft flew a series of flight legs just at the top of the marinestratocumulus cloud decks that cap the mixed layer off the coast of southernCalifornia. Because of the corrugated structure of the cloud-top, the aircraft, which was flown at a nearly constant level and adjusted only to maintain its altitude at the average cloud-top height, was alternately within and above the clouds – roughly half the time in each domain. These legs were used to examine the structure of the cloud-top by compositing the segments on either side of the cloud/clear-air interface, which was identified by the transitions of liquid water measured by the Forward Scattering Spectrometer Probe (either increasing or decreasing) through a threshold of 0.04 × 10^-3 kg m^-3.An equivalent vertical distance (EVD) from the cloud-top was obtained from the horizontal flight legs by estimating the average slope of the cloud-top from the cloud-top radiation temperature. The results show that a near discontinuity occurs in variables across cloud top over an EVD of 0.3 m, but that above this, the air has already been modified by boundary-layer air. Thus, cloud-top is not the limit of mixing of boundary-layer air. This mixing may extend to tens of metres or more. The bulk Richardson number in the vicinity of cloud-top increases from near zero within the cloud to about 1.2 at an EVD of 3–6 m above cloud. Fluctuations of the three velocity components within cloud are nearly equal; above cloud the vertical component structure function is about half the horizontal components. The scalar structure functions are about an order of magnitude higher above cloud than in cloud. The structure parameters of temperature and humidity measured just below cloud-top agree reasonably well with predicted values based on a previously-developed model for the clear convective boundary layer. Above cloud, the scalar structure parameters are much larger, but their interpretation is questionable, since this region does notcontain isotropic turbulence.     

Scalar Transport over Forested Hills

Numerical simulations of scalar transport in neutral flow over forested ridges are performed using both a 1.5-order mixing-length closure scheme and a large-eddy simulation. Such scalar transport (particularly of CO₂) has been a significant motivation for dynamical studies of forest canopy–atmosphere interactions. Results from the 1.5-order mixing-length simulations show that hills for which there is significant mean flow into and out of the canopy are more efficient at transporting scalars from the canopy to the boundary layer above. For the case with a source in the canopy this leads to lower mean concentrations of tracer within the canopy, although they can be very large horizontal variations over the hill. These variations are closed linked to flow separation and recirculation in the canopy and can lead to maximum concentrations near the separation point that exceed those over flat ground. Simple scaling arguments building on the analytical model of Finnigan and Belcher (Q J Roy Meteorol Soc 130:1–29, 2004) successfully predict the variations in scalar concentration near the canopy top over a range of hills. Interestingly this analysis suggests that variations in the components of the turbulent transport term, rather than advection, give rise to the leading order variations in scalar concentration. The scaling arguments provide a quantitative measure of the role of advection, and suggest that for smaller/steeper hills and deeper/sparser canopies advection will be more important. This agrees well with results from the numerical simulations. A large-eddy simulation is used to support the results from the mixing-length closure model and to allow more detailed investigation of the turbulent transport of scalars within and above the canopy. Scalar concentration profiles are very similar in both models, despite the fact that there are significant differences in the turbulent transport, highlighted by the strong variations in the turbulent Schmidt number both in the vertical and across the hill in the large-eddy simulation that are not represented in the mixing-length model.     

Accuracy of the relaxed eddy-accumulation technique,evaluated using CO2 flux measurements

A system capable of measuring the fluxes of trace gases was developed. It is based on a simpler version of the eddy-accumulation technique (EA), known as the relaxed eddy-accumulation technique (REA). It accumulates air samples associated with updrafts and downdrafts at a constant flow rate in two containers for later analysis of the trace gas mean concentration. The flux integration is based on the durations of updraft and downdraft events, rather than on the vertical wind velocity (W) as is the case for EA and eddy-correlation (EC) techniques. The flux, calculated by the REA technique, is equal to the difference in the mean concentration of the trace gas of interest between the upward and downward moving eddies, multiplied by the standard deviation of the vertical wind velocity and an empirical coefficient. CO₂ fluxes measured for 162 half-hour periods over a soybean field by both EC and REA techniques showed excellent agreement (coefficient of determination,R ²=0.92). The slope (0.985) and the intercept (–0.042 mg m^–2 s^–1) were not significantly different from 1 and 0, respectively, at the 5% level; and the standard error of estimate was 0.074 mg m^–2 s^–1. It is also shown that the empirical coefficient can be calculated from either latent or sensible heat fluxes. A model describing the effect on this empirical coefficient of not sampling aroundW equal to zero is proposed.Centre for Land and Biological Resources Research Contribution No. 92-212.     

Atmospheric Stability Effects on Wind Fields and Scalar Mixing Within and Just Above a Subalpine Forest in Sloping Terrain

Air temperature T _a , specific humidity q, CO₂ mole fraction χ _c , and three-dimensional winds were measured in mountainous terrain from five tall towers within a 1 km region encompassing a wide range of canopy densities. The measurements were sorted by a bulk Richardson number Ri _b . For stable conditions, we found vertical scalar differences developed over a “transition” region between 0.05 < Ri _b < 0.5. For strongly stable conditions (Ri _b > 1), the vertical scalar differences reached a maximum and remained fairly constant with increasing stability. The relationships q and χ _c have with Ri _b are explained by considering their sources and sinks. For winds, the strong momentum absorption in the upper canopy allows the canopy sublayer to be influenced by pressure gradient forces and terrain effects that lead to complex subcanopy flow patterns. At the dense-canopy sites, soil respiration coupled with wind-sheltering resulted in CO₂ near the ground being 5–7 μmol mol⁻¹ larger than aloft, even with strong above-canopy winds (near-neutral conditions). We found Ri _b -binning to be a useful tool for evaluating vertical scalar mixing; however, additional information (e.g., pressure gradients, detailed vegetation/topography, etc.) is needed to fully explain the subcanopy wind patterns. Implications of our results for CO₂ advection over heterogenous, complex terrain are discussed.