A wind tunnel study of air flow in waving wheat: Single-point velocity statistics

We analyse single-point velocity statistics obtained in a wind tunnel within and above a model of a waving wheat crop, consisting of nylon stalks 47 mm high and 0.25 mm wide in a square array with frontal area index 0.47. The variability of turbulence measurements in the wind tunnel is illustrated by using a set of 71 vertical traverses made in different locations, all in the horizontally-homogeneous (above-canopy) part of the boundary layer. Ensemble-averaged profiles of the statistical moments up to the fourth order and profiles of Eulerian length scales are presented and discussed. They are consistent with other similar experiments and reveal the existence of large-scale turbulent coherent structures in the flow. The drag coefficient in this canopy as well as in other reported experiments is shown to exhibit a characteristic height-dependency, for which we propose an interpretation. The velocity spectra are analysed in detail; within and just above the canopy, a scaling based on fixed length and velocity scales (canopy height and mean horizontal wind speed at canopy top) is proposed. Examination of the turbulent kinetic energy and shear stress budgets confirms the role of turbulent transport in the region around the canopy top, and indicates that pressure transport may be significant in both cases. The results obtained here show that near the top of the canopy, the turbulence properties are more reminiscent of a plane mixing layer than a wall boundary layer.     

Effects of Unstable Thermal Stratification on Vertical Fluxes of Heat and Momentum in Urban Areas

A Reynolds-averaged Navier–Stokes microscale model is used for the simulation of the effect of unstable thermal stratification on the flow within an aligned configuration of building-like cubes as used in Santiago et al. (Urban Clim 9:115–133, 2014). The spatially-averaged results show increased dispersive fluxes, turbulent length scales and sectional drag coefficient. An extension of K-theory is presented to parametrize the sum of the turbulent and dispersive fluxes, and the length scale and drag coefficient increases are parametrized as functions of the ratio of buoyant and inertial forces. This approach improves the results of urban canopy parametrization simulations inside and above the urban canyon and represents the first attempt to account for the dispersive fluxes and the effect of solar radiation on the flow.     

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.     

A Simple Model for Spatially-averaged Wind Profiles Within and Above an Urban Canopy

This paper deals with the modelling of the flow in the urban canopy layer. It critically reviews a well-known formula for the spatially-averaged wind profile, originally proposed by Cionco in 1965, and provides a new interpretation for it. This opens up a number of new applications for modelling mean wind flow over the neighbourhood scale. The model is based on a balance equation between the obstacle drag force and the local shear stress as proposed by Cionco for a vegetative canopy. The buildings within the canopy are represented as a canopy element drag formulated in terms of morphological parameters such as λ _f and λ _p (the ratios of plan area and frontal area of buildings to the lot area). These parameters can be obtained from the analysis of urban digital elevation models. The shear stress is parameterised using a mixing length approach. Spatially-averaged velocity profiles for different values of building packing density corresponding to different flow regimes are obtained and analysed. The computed solutions are compared with published data from wind-tunnel and water-tunnel experiments over arrays of cubes. The model is used to estimate the spatially-averaged velocity profile within and above neighbourhood areas of real cities by using vertical profiles of λ _f .     

A Dynamic Urban Canopy Parameterization for Mesoscale Models Based on Computational Fluid Dynamics Reynolds-Averaged Navier–Stokes Microscale Simulations

Urban canopy parameterizations (UCPs) are necessary in mesoscale modelling to take into account the effects of buildings on wind and turbulent structures. This study is focused on the dynamical part of UCPs. The main objective is twofold: first, computing important UCP input parameters (turbulent length scales and the sectional drag coefficient) by means of Reynolds-averaged Navier–Stokes (RANS) simulations of turbulent flow; and second, comparing UCP variables with spatially-averaged properties obtained from RANS simulations for the same configurations. The results show the importance of using a suitable parameterization of the drag force for different packing densities. An urban canopy parameterization that is a compromise between simplicity and accuracy is proposed. This scheme accounts for the variation of drag coefficients with packing densities, and has a parameterization of turbulent length scales. The technique adopted ensures that, at least for the simple configurations studied, the urban canopy parameterization gives values of spatially-averaged variables similar to those computed from a more complex simulation, such as RANS that resolves explicitly the flow around buildings.     

Three-dimensional simulations and wind-tunnel experiments on airflow over isolated forest stands

Prediction of windthrow risk to individual or groups of retained trees in harvested stands requires an improved understanding of canopy airflow dynamics. Large-eddy simulations were used to simulate wind-tunnel experiments in two and three dimensions to compare with observations for model validation and to address parameter space considerations for the design of subsequent retention pattern experiments. The three-dimensional simulations were similar to the observed wind-tunnel data for the statistical profiles for but there were greater differences in skewness and kurtosis. These results were obtained using a common leaf-area drag formulation without either skin friction or speed dependent drag that enables scaling with U ₀ (ambient wind speed) and h (height of the canopy). This scaling results in a single non-dimensional parameter h/h _c where h _c (x, y, z) is the momentum range resulting from the canopy drag. The validity of the model scaling was tested using two-dimensional simulations. The irrotational component of the flow (potential flow) was found to be important when defining vertical domain limitations and has significant implications for time dependent flow (i.e. turbulent conditions) when considering retention pattern design. The sudden onset of drag associated with the isolated stand presents some unexpected challenges. The horizontal scales of the shearing instabilities were simulated in two dimensions and found to range between 2h for early times to 7h for later times. The early-time horizontal scales are in the range of logical retention pattern scales and as such need to be taken into account as part of the parameter space, i.e. a range of retention pattern lengths need consideration.     

Statistics of atmospheric turbulence within a natural black spruce forest canopy

Turbulence statistics were measured in a natural black-spruce forest canopy in southeastern Manitoba, Canada. Sonic anemometers were used to measure time series of vertical wind velocity (w), and cup anemometers to measure horizontal wind speed (s), above the canopy and at seven different heights within the canopy. Vertical profiles were measured during 25 runs on eight different days when conditions above the canopy were near-neutral.Profiles of s and of the standard deviation ( _w ) of w show relatively little scatter and suggest that, for this canopy and these stability conditions, profiles can be predicted from simple measurements made above the canopy. Within the canopy, a negative skewness and a high kurtosis of the w-frequency distributions indicate asymmetry and the persistence of large, high-velocity eddies. The Eulerian time scale is only a weak function of height within the canopy.Although w-power spectra above the canopy are similar to those in the free atmosphere, we did not observe an extensive inertial subrange in the spectra within the canopy. Also, a second peak is present that is especially prominent near the ground. The lack of the inertial subrange is likely caused by the presence of sources and sinks for turbulent kinetic energy within our canopy. The secondary spectral peak is probably generated by wake turbulence caused by form drag on the wide, horizontal spruce branches.     

Exploring the Effects of Microscale Structural Heterogeneity of Forest Canopies Using Large-Eddy Simulations

The Regional Atmospheric Modeling System (RAMS)-based Forest Large-Eddy Simulation (RAFLES), developed and evaluated here, is used to explore the effects of three-dimensional canopy heterogeneity, at the individual tree scale, on the statistical properties of turbulence most pertinent to mass and momentum transfer. In RAFLES, the canopy interacts with air by exerting a drag force, by restricting the open volume and apertures available for flow (i.e. finite porosity), and by acting as a heterogeneous source of heat and moisture. The first and second statistical moments of the velocity and flux profiles computed by RAFLES are compared with turbulent velocity and scalar flux measurements collected during spring and winter days. The observations were made at a meteorological tower situated within a southern hardwood canopy at the Duke Forest site, near Durham, North Carolina, U.S.A. Each of the days analyzed is characterized by distinct regimes of atmospheric stability and canopy foliage distribution conditions. RAFLES results agreed with the 30-min averaged flow statistics profiles measured at this single tower. Following this intercomparison, two case studies are numerically considered representing end-members of foliage and midday atmospheric stability conditions: one representing the winter season with strong winds above a sparse canopy and a slightly unstable boundary layer; the other representing the spring season with a dense canopy, calm conditions, and a strongly convective boundary layer. In each case, results from the control canopy, simulating the observed heterogeneous canopy structure at the Duke Forest hardwood stand, are compared with a test case that also includes heterogeneity commensurate in scale to tree-fall gaps. The effects of such tree-scale canopy heterogeneity on the flow are explored at three levels pertinent to biosphere-atmosphere exchange. The first level (zero-dimensional) considers the effects of such heterogeneity on the common representation of the canopy via length scales such as the zero-plane displacement, the aerodynamic roughness length, the surface-layer depth, and the eddy-penetration depth. The second level (one-dimensional) considers the normalized horizontally-averaged profiles of the first and second moments of the flow to assess how tree-scale heterogeneities disturb the entire planar-averaged profiles from their canonical (and well-studied planar-homogeneous) values inside the canopy and in the surface layer. The third level (three-dimensional) considers the effects of such tree-scale heterogeneities on the spatial variability of the ejection-sweep cycle and its propagation to momentum and mass fluxes. From these comparisons, it is shown that such microscale heterogeneity leads to increased spatial correlations between attributes of the ejection-sweep cycle and measures of canopy heterogeneity, resulting in correlated spatial heterogeneity in fluxes. This heterogeneity persisted up to four times the mean height of the canopy (h _c ) for some variables. Interestingly, this estimate is in agreement with the working definition of the thickness of the canopy roughness sublayer (2h _c –5h _c ).     

Streamwise Versus Spanwise Spacing of Obstacle Arrays: Parametrization of the Effects on Drag and Turbulence

A Reynolds-averaged Navier–Stokes model is used to investigate the evolution of the sectional drag coefficient and turbulent length scales with the layouts of aligned arrays of cubes. Results show that the sectional drag coefficient is determined by the non-dimensional streamwise distance (sheltering parameter), and the non-dimensional spanwise distance (channelling parameter) between obstacles. This is different than previous approaches that consider only plan area density $(\lambda _\mathrm{p})$ . On the other hand, turbulent length scales behave similarly to the staggered case (e. g. they are function of $\lambda _\mathrm{p}$ only). Analytical formulae are proposed for the length scales and for the sectional drag coefficient as a function of sheltering and channelling parameters, and implemented in a column model. This approach demonstrates good skill in the prediction of vertical profiles of the spatially-averaged horizontal wind speed.     

Spatio-temporal Variability of Surface-layer Turbulent Fluxes Over the Bay of Bengal and Arabian sea During the ICARB Field Experiment

We report the spatio-temporal variability of surface-layer turbulent fluxes of heat, moisture and momentum over the Bay of Bengal (BoB) and the Arabian Sea (AS) during the Integrated Campaign for Aerosols, gases Radiation Budget (ICARB) field experiment. The meteorological component of ICARB conducted during March – May 2006 onboard the oceanic research vessel Sagar Kanya forms the database for the present study. The bulk transfer coefficients and the surface-layer fluxes are estimated using a modified bulk aerodynamic method, and then the spatio-temporal variability of these air-sea interface fluxes is discussed in detail. It is observed that the sensible and latent heat fluxes over the AS are marginally higher than those over the BoB, which we attribute to differences in the prevailing meteorological conditions over the two oceanic regions. The values of the wind stress, sensible and latent heat fluxes are compared with those obtained for the Indian Ocean Experiment (INDOEX) period. The variation of drag coefficient (C _D ), exchange coefficients of sensible heat and moisture (C _H = C _E ) and neutral drag coefficient (C _DN ) with wind speed is also discussed.

---

     

Aerodynamic roughness of vegetated surfaces

Available experimental results indicate that as the density of roughness elements over a horizontally homogeneous surface is varied, the roughness length, z ₀, varies in a manner that exhibits a maximum at intermediate density values. In an attempt to explain this behaviour, the available analytical solutions for the wind profile inside dense homogeneous canopies were reviewed. The review indicated that the variation of z ₀ with density depends on the interrelationship between the leaf density, a, and the mixing length, l. In view of this finding, a numerical model was devised based on a simple rule for constructing mixing-length profiles in the canopy. The rule states that the actual value of l is the maximum possible under the two constraints: l l _i and ¦dl/dz¦ k, where k is the von Karman constant and the intrinsic mixing length, l _i, is a function of the local internal structure of the canopy. The model which ensures a smooth transition from dense to thin canopy, was used to reproduce the observed maximum of z ₀. The model is also capable of handling vertically non-homogeneous canopies.     

Nocturnal Intermittent Coupling Between the Interior of a Pine Forest and the Air Above It

A 1-year set of measurements of CO₂ and energy turbulent fluxes above and within a 25-m pine forest in southern Brazil is analyzed. The study focuses on the coupling state between two levels and its impact on flux determination by the eddy-covariance method. The turbulent series are split in their typical temporal scales using the multiresolution decomposition, a method that allows proper identification of the time scales of the turbulent events. Initially, four case studies are presented: a continually turbulent, a continually calm, a calm then turbulent, and an intermittent night. During transitions from calm to turbulent, large scalar fluxes of opposing signs occur at both levels, suggesting the transference of air accumulated in the canopy during the stagnant period both upwards and downwards. Average fluxes are shown for the entire period as a function of turbulence intensity and a canopy Richardson number, used as an indicator of the canopy coupling state. Above the canopy, CO₂ and sensible heat fluxes decrease in magnitude both at the neutral and at the very stable limit, while below the canopy they increase monotonically with the canopy Richardson number. Latent heat fluxes decrease at both levels as the canopy air becomes more stable. The average temporal scales of the turbulent fluxes at both levels approach each other in neutral conditions, indicating that the levels are coupled in that case. Average CO₂ fluxes during turbulent periods that succeed very calm ones are appreciably larger than the overall average above the canopy and smaller than the average or negative within the canopy, indicating that the transfer of air accumulated during calm portions at later turbulent intervals affects the flux average. The implications of this process for mean flux determination are discussed.     

étude expérimentale de la couche limite de montagne

The characteristics of the boundary layer over complex terrain (Lannemezan - lat.: 43.7° N and, long.: 0.7 ° E) are analyzed for various scales, using measurements obtained during the COCAGNE Experiment. In this first part, the dynamic characteristics of the flow are studied with respect to atmospheric stability and the relief at small (~20 km) and medium scales (~100 km). These relief scales depend on the topographical profile of the Lannemezan Plateau along the dominant axis of the wind (E-W) and the Pyrénées Mountains located at the south of the experimental site. The terrain heterogeneities have a standard deviation of ~48 m and a wavelength of ~2 km.The averaged vertical profiles of wind speed and direction over the heterogeneous terrain are analyzed. The decrease of wind speed within the boundary layer is greater than over flat terrain (WANGARA Experiment). However, a comparison between ETTEX (complex terrain) and COCAGNE vertical wind speed profiles shows good agreement during unstable conditions. In contrast, during neutral conditions a more rapid increase with normalized height is found with COCAGNE than with ETTEX and WANGARA data. The vertical profiles of wind direction reveal an influence of the Pyrénées Mountains on the wind flow. The wind rotation in the BL is determined by the geostrophic wind direction-Pyrénées axis angle (negative deviation) as the geostrophic wind is connected with the Mountain axis.When the geostrophic wind does not interact with the Pyrénées axis, the mean and turbulent wind flow characteristics (drag coefficient C _D, friction velocity u _*) depend on the topography of the plateau. When the wind speed is strong (>6 m s ^-1), an internal boundary layer is generated from the leading edge of the Plateau.     

Numerical Simulations of the Microscale Heterogeneities of Turbulence Observed on a Complex Site

We simulate the microscale heterogeneities of turbulent variables observed at a complex site for different wind directions. The atmospheric computational fluid dynamics (CFD) results are compared with an ensemble of 36 months of data collected at the experimental site SIRTA “Site Instrumental de Recherche par Télédétection Atmosphérique”, located near Paris (France) in a semi-urban environment. The experimental data show that the normalized turbulent kinetic energy (TKE) k/U ² (where k is TKE and U is the wind speed) at 10-m height, for two different locations, is highly dependent on wind direction and strongly influenced by trees. These measurements show a strong increase of the normalized TKE downstream of the forest canopies with a large variability within the 36-month period in part due to the variation of the tree foliage. The numerical simulations are carried out using the CFD code Code_Saturne with the standard k?ε closure, in neutral stratification. The buildings are taken into account explicitly in the mesh and the forested areas are modelled with two approaches: the classical roughness wall law and a drag porosity. A comparison has been performed between the calculated values and the median of measured values of the normalized TKE and the normalized friction velocity, for each wind sector of 10°. A very good agreement is obtained with the drag porosity model, whereas the classical roughness law leads to a strong underestimation downstream of the forested areas. However, this large improvement of the results using the drag porosity model can only be obtained with a refinement of the grid, especially in forested areas, and an accurate land-use map.     

Turbulent exchange above a pine forest,I: Fluxes and gradients

Measurements of gradients of wind, temperature and humidity and of the corresponding turbulent fluxes have been carried out over a sparse pine forest at Jädra»s in Sweden. In order to ascertain that correct gradient estimates were obtained, two independent measuring systems were employed: one system with sensors at 10 fixed levels on a 51 m tower and another with reversing sensors for temperature and humidity, covering the height interval 23 to 32 m. Turbulent fluxes were measured at three levels simultaneously. Data from three field campaigns: in June 1985, June 1987 and September 1987 have been analyzed. The momentum flux is found on the average to be virtually constant from tree top level, at 20 to 50 m. The average fluxes of sensible and latent heat are not so well behaved. The ratio of the non-dimensional gradients of wind and temperature to their corresponding values under ideal conditions (low vegetation) are both found to be small immediately above the canopy (about 0.3 for temperature and 0.4 for wind). With increasing height, the ratios increase, but the values vary substantially with wind direction. The ratios are not found to vary systematically with stability (unstable stratification only studied). The ratio of the non-dimensional humidity gradient to the corresponding non-dimensional potential temperature gradient (equivalent to k _h /k _w ) is found to be unity for (z – d)/L _v less than about –0.1 and about 1.4 for near neutral stratification, but the scatter of the data is very large.     

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.     

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.     

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.     

Development and Validation of a Lagrangian Probability Density Function Model of Horizontally-Homogeneous Turbulence Within and Above Plant Canopies

Second- and third-order turbulence closure models have met with mixed success when applied to the prediction of turbulent flows within and above plant canopies and model predictions are typically no better than those obtained using simpler two-equation (k-e){(k-\varepsilon)} models. This is because the local gradient diffusion approximation (a crucial model requirement) cannot represent accurately turbulent transport that is dominated by the presence of ejections and sweeps whose length scales are comparable with the canopy height. To make progress, turbulent transport must be treated without approximation, as in Lagrangian probability density function (PDF) models. This study is the first to develop and to validate a PDF model of horizontally-homogeneous canopy flow. The model relies upon a prescribed length-scale that has been used elsewhere in the modelling of turbulent flows. Model predictions compare favourably with measurements of neutrally stratified turbulent flows within and above canopies of mature corn and forested eucalypt trees.