Turbulence structure in a deciduous forest

Three-dimensional wind velocity components were measured at two levels above and at six levels within a fully-leafed deciduous forest. Greatest shear occurs in the upper 20% of the canopy, where over 70% of the foliage is concentrated. The turbulence structure inside the canopy is characterized as non-Gaussian, intermittant and highly turbulent. This feature is supported by large turbulence intensities, skewness and kurtosis values and by the large infrequent sweeps and ejections that dominate tangential momentum transfer. Considerable day/night differences were observed in the vertical profiles of the mean streamwise wind velocity and turbulence intensities since the stability of the nocturnal boundary layer dampens turbulence above and within the canopy.     

Turbulence in an almond orchard: Spatial variations in spectra and coherence

The spatial variability of turbulence in a fully-leafed almond orchard was studied. Two three-dimensional sonic anemometers were used to measure turbulence spectra and coherence at different vertical and lateral separations inside the canopy. Peak frequencies of the horizontal velocity components, normalized by local horizontal wind speed, are greater in the canopy crown than in the trunkspace. Peak-normalized frequencies for the vertical velocity power spectra are similar in the canopy crown and in the subcanopy trunkspace. Spectral slopes in the inertial subrange are more negative than those predicted with Kolmogorov's -2/3 theory. It is thought that the foliage elements act to short-circuit the eddy cascade. Lateral separation of the instruments in the subcanopy trunkspace has little effect on the shape of the velocity spectra. On the other hand, lateral and vertical velocity coherences between spatially separated sensors are low inside the canopy. These low coherences are due to the Eulerian length scales being of the same order of magnitude as the separation distances of the anemometers. Phase angles between velocity components are about zero for small separation distances. When the two instruments are separated by 9 m and one instrument is positioned in a row while the other is between two rows, vertical velocities are about 180 deg out of phase and the streamwise velocities are about 40 to 60 deg out of phase. These data support the contention that preferred differences occur between within- and between-row wind flow regimes.     

Wave-like wind fluctuations observed in the stable surface layer over a plant canopy(2)

Observations of wind velocity and air temperature fluctuations were made in the nocturnal surface inversion layer over a sorghum field. Wave-like fluctuations of temperature and wind velocity with a period of 15–20 min were observed for about 2 hours, 3 to 5 hours before sunrise. Wave-like fluctuations of temperature were observed in the air layer above and within a plant canopy and were most noticeable at the top of the plant canopy. Spectral analysis of temperature and wind velocity fluctuations reveals a separation of energy into wave-like and turbulent fluctuations. Cospectral analysis shows that for both momentum and heat, vertical transports are partitioned almost equally in the frequency ranges characteristic of wave-like and turbulent fluctuations. This suggests wave- turbulence interactions at low frequencies in the air layer near a plant canopy.     

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.

     

Mean Flow Near Edges and Within Cavities Situated Inside Dense Canopies

A streamfunction-vorticity formulation is used to explore the extent to which turbulent and turbulently inviscid solutions to the mean momentum balance explain the mean flow across forest edges and within cavities situated inside dense forested canopies. The turbulent solution is based on the mean momentum balance where first-order closure principles are used to model turbulent stresses. The turbulently inviscid solution retains all the key terms in the mean momentum balance but for the turbulent stress gradients. Both exit and entry versions of the forest edge problem are explored. The turbulent solution is found to describe sufficiently the bulk spatial patterns of the mean flow near the edge including signatures of different length scales reported in canopy transition studies. Next, the ‘clearing inside canopy’ or the so-called ‘cavity’ problem is solved for the inviscid and turbulent solutions and then compared against flume experiments. The inviscid solution describes the bulk flow dynamics in much of the zones within the cavity. In particular, the solution can capture the correct position of the bulk recirculation zone within the cavity, although with a weaker magnitude. The inviscid solution cannot capture the large vertical heterogeneity in the mean velocity above the canopy, as expected. These features are better captured via the first-order closure representation of the turbulent solution. Given the ability of this vorticity formulation to capture the mean pressure variations and the mean advective acceleration terms, it is ideal for exploring the distributions of scalars and roughness-induced flow adjustments on complex topography.     

A spatial length scale analysis of turbulent temperature and velocity fluctuations within and above an orchard canopy

Turbulent flow within and above an almond orchard was measured with three-dimensional wind sensors and fine-wire thermocouple sensors arranged in a horizontal array. The data showed organized turbulent structures as indicated by coherent asymmetric ramp patterns in the time series traces across the sensor array. Space-time correlation analysis indicated that velocity and temperature fluctuations were significantly correlated over a transverse distance more than 4m. Integral length scales of velocity and temperature fluctuations were substantially greater in unstable conditions than those in stable conditions. The coherence spectral analysis indicated that Davenport's geometric similarity hypothesis was satisfied in the lower frequency region. From the geometric similarity hypothesis, the spatial extents of large ramp structures were also estimated with the coherence functions.     

A spectral and lag-correlation analysis of turbulence in a deciduous forest canopy

The processes influencing turbulence in a deciduous forest and the relevant length and time scales are investigated with spectral and cross-correlation analysis. Wind velocity power spectra were computed from three-dimensional wind velocity measurements made at six levels inside the plant canopy and at one level above the canopy. Velocity spectra measured within the plant canopy differ from those measured in the surface boundary layer. Noted features associated with the within-canopy turbulence spectra are: (a) power spectra measured in the canopy crown peak at higher wavenumbers than do those measured in the subcanopy trunkspace and above the canopy; (b) peak spectral values collapse to a relatively universal value when scaled according to a non-dimensional frequency comprised of the product of the natural frequency and the Eulerian time scale for vertical velocity; (c) at wavenumbers exceeding the spectral peak, the slopes of the power spectra are more negative than those observed in the surface boundary layer; (d) Eulerian length scales decrease with depth into the canopy crown, then increase with further depth into the canopy; (e) turbulent events below crown closure are more correlated with turbulent events above the canopy than are those occurring in the canopy crown; and (f) Taylor's frozen eddy hypothesis is not valid in a plant canopy. Interactions between plant elements and the mean wind and turbulence alter the processes that produce, transport and remove turbulent kinetic energy and account for the noted observations.     

Mean Flow and Turbulence Characteristics in an Urban Roughness Sublayer

In this study, a detailed model of an urban landscape has been re-constructed inthe wind tunnel and the flow structure inside and above the urban canopy has beeninvestigated. Vertical profiles of all three velocity components have been measuredwith a Laser-Doppler velocimeter, and an extensive analysis of the measured meanflow and turbulence profiles carried out. With respect to the flow structure inside thecanopy, two types of velocity profiles can be distinguished. Within street canyons,the mean wind velocities are almost zero or negative below roof level, while closeto intersections or open squares, significantly higher mean velocities are observed.In the latter case, the turbulent velocities inside the canopy also tend to be higherthan at street-canyon locations. For both types, turbulence kinetic energy and shearstress profiles show pronounced maxima in the flow region immediately above rooflevel.Based on the experimental data, a shear-stress parameterization is proposed, inwhich the velocity scale, u_s, and length scale, z_s, are based on the level and magnitude of the shear stress peak value. In order to account for a flow region inside the canopy with negligible momentum transport, a shear stress displacement height, d_s, is introduced. The proposed scaling and parameterization perform well for the measured profiles and shear-stress data published in the literature.The length scales derived from the shear-stress parameterization also allowdetermination of appropriate scales for the mean wind profile. The roughnesslength, z₀, and displacement height, d₀, can both be described as fractions of the distance, z_s - d_s, between the level of the shear-stress peak and the shear-stress displacement height. This result can be interpreted in such a way that the flow only feels the zone of depth z_s - d_s as the roughness layer. With respect to the lower part of the canopy (z < d_s) the flow behaves as a skimming flow. Correlations between the length scales z_s and d_s and morphometric parameters are discussed.The mean wind profiles above the urban structure follow a logarithmic windlaw. A combination of morphometric estimation methods for d₀ and z₀ with wind velocity measurements at a reference height, which allow calculation of the shear-stress velocity, u_*, appears to be the most reliable and easiest procedure to determine mean wind profile parameters. Inside the roughnesssublayer, a local scaling approach results in good agreement between measuredand predicted mean wind profiles.     

Relationship between turbulent flux and variance in the urban canopy

The relationship between fluxes and variances has been widely studied in the frame of the similarity theory. However, the dynamics of urban morphology should be highlighted in the urban canopy. Therefore, the relationship between fluxes and variances of turbulent quantities (wind velocity, temperature, carbon dioxide, and water vapor concentration) without stability parameters has been investigated with data collected from the Beijing 325-m meteorological tower in the urban canopy. Though the similarity theory is available for the statistics of wind velocity and temperature, urban land use plays an important role in the linear relationship between momentum fluxes and velocity variances, while the linear relationship between sensible heat fluxes and temperature variances is influenced at large by stratification. As for carbon dioxide and water vapor, the linear relationship between their fluxes and variances is a better way compared to the similarity theory to describe their motion, though the plots are a little scattered owing to their fugitive source.     

Influence of non-flat terrain and wind direction shear on canopy turbulence

We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk _w and r _uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.     

一方程模型在植被冠层中的应用

利用建立的一方程模式对植被气象场和湍流场进行了模拟计算。结果表明：在植被叶面积密度最大值处，由植被叶面积产生的阻力使风速急剧减小，而后平缓接近地面风速值。Reynolds应力从植被顶部向下剧烈减小，在Z／Hc值为0．4～0．6之间衰减迅速；冠层下部风速小且切变弱，湍流未能充分发展，因而动量输送甚微。由于在植被内部环境中存在着湍流通量的辐散或辐合的现象，湍流强度从植被底部开始由下至上逐渐增强，在冠层上方则基本保持不变。     

Daytime Turbulence Statistics above a Steep Forested Slope

Six levels of simultaneously sampled ultrasonic data are used to analyse the turbulence structure within a mixed forest of 13 m height on a steep slope (35°) in an alpine valley. The data set is compared to other studies carried out over forests in more ideal, flat terrain. The analysis is carried out for 30-min mean data, joint probability distributions, length scales and spectral characteristics.Thermally induced upslope winds and cold air drainage lead to a wind speed maximum within the trunk space. Slope winds are superimposed on valley winds and the valley-wind component becomes stronger with increasing height. Slope and valley winds are thus interacting on different spatial and time scales leading to a quite complex pattern in momentum transport that differs significantly from surface-layer characteristics. Directional shear causes lateral momentum transports that are in the same order or even larger than the longitudinal ones. In the canopy, however, a sharp attenuation of turbulence is observed. Skewed distributions of velocity components indicate that intermittent turbulent transport plays an important role in the energy distribution.Even though large-scale pressure fields lead to characteristic features in the turbulent structure that are superimposed on the canopy flow, it is found that many statistical properties typical of both mixing layers and canopy flow are observed in the data set.     

Turbulent Structures in a Pine Forest with a Deep and Sparse Trunk Space: Stand and Edge Regions

Forested landscapes often exhibit large spatial variability in vertical and horizontal foliage distributions. This variability may affect canopy-atmosphere exchanges through its action on the development of turbulent structures. Here we investigate in neutral stratification the turbulent structures encountered in a maritime pine forest characterized by a high, dense foliated layer associated with a deep and sparse trunk space. Both stand and edge regions are considered. In situ measurements and the results of large-eddy simulations are used and analyzed together. In stand conditions, far from the edge, canopy-top structures appear strongly damped by the dense crown layer. Turbulent wind fluctuations within the trunk space, where the momentum flux vanishes, are closely related to these canopy-top structures through pressure diffusion. Consequently, autocorrelation and spectral analyses are not quite appropriate to characterize the vertical scale of coherent structures in this type of canopy, as pressure diffusion enhances the actual scale of structures. At frequencies higher than those associated with canopy-top structures, wind fluctuations related to wake structures developing behind tree stems are observed within the trunk space. They manifest themselves in wind velocity spectra as secondary peaks in the inertial subrange region, confirming the hypothesis of spectral short-cuts in vegetation canopies. In the edge region specific turbulent structures develop just below the crown layer, in addition to canopy-top structures. They are generated by the wind shear induced by the sub-canopy wind jet that forms at the edge. These structures provide a momentum exchange mechanism similar to that observed at the canopy top but in the opposite direction and with a lower magnitude. They may develop as in plane mixing-layer flows, with some perturbations induced by canopy-top structures. Wake structures are also observed within the trunk space in the edge region.     

A note on vertical coherence of streamwise turbulence inside and above a model plant canopy

Measurements of longitudinal turbulent velocity were made at pairs of levels inside and above a model plant canopy in the wind tunnel. It was found that above approximately the zero-plane displacement level, the coherence and phase results were similar in many respects to atmospheric data, but that deviations from this behaviour appeared deeper in the model canopy.     

An Analytical Model for Mean Wind Profiles in Sparse Canopies

Existing analytical models for mean wind profiles within canopies are applicable only in dense canopy scenarios, where all momentum is absorbed by canopy elements and, hence, the effect of the ground on turbulent mixing is not important. Here, we propose a new analytical model that can simulate mean wind profiles within sparse canopies under neutral conditions. The model adopts a linearized canopy-drag parametrization and a first-order turbulence closure scheme taking into account the effects of both the ground and canopy elements on turbulent mixing. The resulting wind profile within a sparser canopy appears to be more like a logarithmic form, with the no-slip condition at the ground being satisfied. The analytical solution converges exactly to the standard surface-layer logarithmic wind profile in the case of zero canopy density (i.e., no-canopy scenario) and tends to be an exponential wind profile for a dense canopy; this feature is unique compared with existing analytical models for canopy wind profiles. Results from the new model are in good agreement with those from laboratory experiments and numerical simulations.     

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.     

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.     

Particle Image Velocimetry Measurements of Turbulent Flow Within Outdoor and Indoor Urban Scale Models and Flushing Motions in Urban Canopy Layers

We investigate the spatial characteristics of urban-like canopy flow by applying particle image velocimetry (PIV) to atmospheric turbulence. The study site was a Comprehensive Outdoor Scale MOdel (COSMO) experiment for urban climate in Japan. The PIV system captured the two-dimensional flow field within the canopy layer continuously for an hour with a sampling frequency of 30 Hz, thereby providing reliable outdoor turbulence statistics. PIV measurements in a wind-tunnel facility using similar roughness geometry, but with a lower sampling frequency of 4 Hz, were also done for comparison. The turbulent momentum flux from COSMO, and the wind tunnel showed similar values and distributions when scaled using friction velocity. Some different characteristics between outdoor and indoor flow fields were mainly caused by the larger fluctuations in wind direction for the atmospheric turbulence. The focus of the analysis is on a variety of instantaneous turbulent flow structures. One remarkable flow structure is termed ‘flushing’, that is, a large-scale upward motion prevailing across the whole vertical cross-section of a building gap. This is observed intermittently, whereby tracer particles are flushed vertically out from the canopy layer. Flushing phenomena are also observed in the wind tunnel where there is neither thermal stratification nor outer-layer turbulence. It is suggested that flushing phenomena are correlated with the passing of large-scale low-momentum regions above the canopy.     

植被内部及其上方湍流场的数值模拟

植被内部及其上方的湍流流场对于了解植被与大气之间的动量、热量和质量交换过程极其重要。本文把高阶湍流封闭模型的Reynolds应力方程模型(RSM)应用于植被湍流的计算,得到了风速、湍流动能、Reynolds应力及能量耗散率的垂直分布,与现场观测数据比较,甚为满意。