Flows through forest canopies in complex terrain

Recent progress on boundary layer flow within and above tall forest canopies in complex terrain is reviewed from the perspective of developing methods to interpret carbon dioxide fluxes from tower measurements in real terrain. Two examples of complex terrain are considered in detail: a forest edge, which exemplifies nonuniform forests, and hilly terrain, which can lead to drainage currents at night. Dynamical arguments show that, when boundary layer winds approach a forest edge, the mean wind adjusts on a length scale of approximately 3L(c), where L(c) is the canopy drag length scale, which depends inversely on the leaf area density of the forest. Over a further distance that also scales on L(c), turbulence in the flow adjusts, and the mixing and transport in the canopy approaches the homogeneous limit. Even low hills change the neutral flow within and above the forest canopy substantially. When the canopy is tall, pressure gradients drive flow up both the upwind and downwind slopes of the hill, leading to an ejection of air out of the top of the canopy just downwind of the crest. This flow at the crest can then advect scalar out of the top of the forest, leading to large variations in the flux of scalar across the hill. At night, when the air near the ground cools and becomes stably stratified, turbulence within the canopy can collapse, even when the flow above the canopy remains turbulent. This leads to a decoupling of the air motions within the canopy from those above. The air above the canopy can then continue to pass up and over the hill, as it does in the neutral case, but at the same time, air within the canopy drains down the hill slopes as drainage currents. These analyses will help us understand when flux towers are reliably measuring the net ecosystem exchange and suggest ways of correcting the flux tower data in more complex situations.     

Turbulence characteristics within sparse and dense canopies

Boundary layer interactions with canopies control various environmental processes. In the case of dense and homogeneous canopies, the so-called mixing layer analogy is most generally used. When the canopy becomes sparser, a transition occurs between the mixing layer and the boundary layer perturbed by interactions between element wakes. This transition has still to be fully understood and characterized. The experimental work presented here deals with the effect of the canopy density on the flow turbulence and involves an artificial canopy placed in a fully developed turbulent boundary layer. One and two-component velocity measurements are performed, both within and above the canopy. The influence of the spacing between canopy elements is studied. Longitudinal velocity statistical moments and Reynolds stresses are calculated and compared to literature data. For spacings greater than the canopy height, evidences of this transition are found in the evolution of the skewness factor, shear length scale and mixing length.     

Development of a rotorcraft dust-emission parameterization using a CFD model

Flights of rotorcraft over the desert floor can result in significant entrainment of particulate matter into the atmospheric boundary layer. Continuous or widespread operation can lead to local and regional impacts on visibility and air quality. To account for this pollutant source in air quality models, a parameterization scheme is needed that addresses the complex vertical distribution of dust ejected from the rotorcraft wake into the atmospheric surface layer. A method to parameterize the wind and turbulence fields and shear stress at the ground is proposed here utilizing computational fluid dynamics and a parameterized rotor model. Measurements taken from a full scale experiment of rotorcraft flight near the surface are compared to the simulation results in a qualitative manner. The simulation is shown to adequately predict the forward detachment length of the induced ground jet compared to the measured detachment lengths. However, the simulated ground vortex widths and vorticity deviate substantially from the measured values under a range of flight speeds. Results show that the method may be applicable for air quality modeling assuming slow airspeeds of the rotorcraft, with advance ratios of 0.005–0.02.     

A closure approach with multi-scale mixing length,and its application to wind and shear stress profiles within the atmospheric boundary layer

A method to determine flow specific first-order closure for the turbulent flux of momentum in the atmospheric boundary layer (ABL) is presented. This is based on the premise that eddy viscosity is a flow rather than a fluid property, and the physically more realistic assumption that the transfer of momentum and other scalar quantities in a turbulent flow takes place by a large, but finite number of length scales, than the often used single length scale, the ‘mixing length’. The resulting eddy viscosity is flow specific and when applied to the study of the ABL, yields the vertical profiles of shear stress and mean wind velocity in good agreement with observations. The method may be extended to other types of turbulent flows, however it should be recognized that each type of flow may yield a different eddy viscosity profile. Using the derived eddy viscosity the paper presents simple analytical solutions of the ABL equations to determine observationally consistent wind speed and shear stress profiles in the ABL for a variety of practical applications including air pollution modelling.     

Embedded large eddy simulation approach for pollutant dispersion around a model building in atmospheric boundary layer

In the present article, the potential of embedded large eddy simulation (ELES) approach to reliably predict pollutant dispersion around a model building in atmospheric boundary layer is assessed. The performance of ELES in comparison with large eddy simulation (LES) is evaluated in several ways. These include a number of qualitative and quantitative comparisons of time-averaged and instantaneous results with wind tunnel measurements supplemented by statistical data analyses using scatter plots and standard evaluation metrics. Results obtained by both LES and ELES approaches show very good agreement with the experiment. However, addition of turbulence to mean flow at Reynolds averaged Navier–Stokes (RANS)–LES interface in ELES approach not only increases the turbulence intensity, it also results in larger values of turbulent kinetic energy (TKE) as well as a shorter reattachment length in the wake region. Accordingly, higher levels of TKE predicted by ELES increase the local intensity of concentration leading to shorter plume shapes as compared with LES. In general, ELES shows better agreement with experiment on the surfaces of model building and also in the downstream wake region. In terms of computational costs, the CPU time required to obtain statistical values in ELES is about 49 % lower than that of LES and the number of iterations per time step is also reduced by 55 % as compared with LES.     

Flow structure and turbulence characteristics downstream of a spanwise suspended linear array

Laboratory experiments are conducted to quantify the mean flow structure and turbulence properties downstream of a spanwise suspended linear array in a uniform ambient water flow using Particle Tracking Velocimetry. Eighteen experimental scenarios, with four depth ratios (array depth to water column depth) of 0.35, 0.52, 0.78, and 0.95 and bulk Reynolds number (length scale is the array depth) from 11,600 to 68,170, are investigated. Three sub-layers form downstream of the array: (1) an internal wake zone, where the time-averaged velocity decreases with increasing distance downstream, (2) a shear layer which increases in vertical extent with increasing distance downstream of the array, and the rate of the increase is independent of the bulk Reynolds number or the depth ratio, and (3) an external wake layer with enhanced velocity under the array. The location of the shear layer is dependent on the depth ratio. The spatially averaged and normalized TKE of the wake has a short production region, followed by a decay region which is comparable to grid turbulence decay and is dependent on the depth ratio. The results suggest that the shear layer increases the transfer of horizontal momentum into the internal wake zone from the fluid outside of the array and that the turbulence in the internal wake zone can be modeled similarly to that of grid turbulence.     

A Simple Model to Determine In Situ Mixing Height Growth from Surface Measurements

A generic In Situ Mixing Height Growth (IMG) model, capable of predicting the real-time growth of the mixed layer and its diurnal evolution from routinely observed simple surface meteorological is developed. The algorithm for the determination of temporally growing daytime mixing height includes both convective and mechanical turbulence contributions. It accounts for neutral as well as height varying potential temperature gradients above the mixed layer. For thermally stable and mechanically dominated unstable night time Atmospheric Boundary Layer (ABL) the module uses similarity formulae based on the wind velocity [1]), the Monin—Obukhov length [2], and the Coriolis parameter. In the convective case simple slab model is integrated, based on initial lapse rate and the surface heat flux. The lapse rate is evaluated on the basis of vertical atmospheric stability, surface type and surface temperature. This differentiates the IMG model from other existing mixing height models that need initial measured lapse rate for calculation. IMG model is site specific as it calculates the radiative incoming heat flux depending on the solar declination estimates based on-site latitude and longitude. The IMG model is applied to calculate mixing height for India by using surface data (viz. wind speed, surface temperature, surface type) from 152 surface meteorological stations. Results have been evaluated with radiosonde mixing height data procured from 18 upper air stations. Sensitivity analysis of the model with respect to various parameters is performed. The model is formulated after reviewing presently available radiosonde mixing height data in India and can satisfactorily provide an alternative means of estimating mixing height for air pollution dispersion models.     

CO2 fluxes near a forest edge: a numerical study

In contrast with recent advances on the dynamics of the flow at a forest edge, few studies have considered its role on scalar transport and, in particular, on CO2 transfer. The present study addresses the influence of the abrupt roughness change on forest atmosphere CO2 exchange and contrasts the concentration and flux fields against those of a uniform forested surface. We use an atmospheric boundary layer two-equation closure model that accounts for the flow dynamics and vertical divergence of CO2 sources/sinks within a plant canopy. This paper characterizes the spatial variation of CO2 fluxes as a function of both sources/sinks distribution and the vertical structure of the canopy. Results suggest that the ground source plays a major role in the formation of wave-like vertical CO2 flux behavior downwind of a forest edge, despite the fact that the contribution of foliage sources/sinks changes monotonously. Such a variation is caused by scalar advection in the trunk space and reveals itself as a decrease or increase in vertical fluxes over the forest relative to carbon dioxide exchange of the underlying forest. The effect was more pronounced in model forests where the leaf area is concentrated in the upper part of the canopy. These results can be useful both for interpretation of existing measurements of net ecosystem exchange of CO2 (NEE) from flux towers in limited fetch conditions and in planning future CO2 transport experiments.     

Double-Arranged Mound-Mounted Shelterbelts: Influence of Porosity on Wind Reduction between the Shelters

Using a two-dimensional LDA-System, measurements in an atmospheric boundary layer wind tunnel have been made in the intermediate field of double-arranged, mound-mounted shelterbelts. The porosity of the shelterbelts ranged from 0% to 52%. For each porosity, seven different distances of the two windbreaks between 5 and 25 times the windbreak height were investigated. The efficiency of the windbreak systems is assessed by means of protection parameters considering the change of horizontal momentum on bodies immersed in the flow, of pedestrian comfort conditions and of momentum exchange in vertical direction. Furthermore recirculation zones and the change of kinetic energy are shown. The measurement results are given in the form of contour plots indicating flow quantities and protection volumes between the two windbreaks. The shelter efficiency of the investigated windbreak systems is compared to that of adequate single line windbreaks. It is shown that for certain windbreak configurations the wind-protected area in the intermediate field of consecutively arranged windbreaks becomes larger in comparison to a single windbreak.     

Physical modeling of complex airflows developing above rural terrains

The rural atmospheric boundary layer (ABL) flow was reproduced in a wind tunnel at three different simulation length scales to investigate possible effects of the simulation length scale on flow characteristics. Performance of truncated vortex generators developed for part-depth ABL wind-tunnel simulations was tested in rural terrain exposure against the full-size Counihan vortex generators. A procedure to design the ABL developing above rural type terrain has been described. The 1:395 and 1:236 simulations were created as full-depth simulations, i.e., wind characteristics throughout an entire ABL were reproduced in the wind tunnel. The 1:208 simulation was a part-depth simulation, i.e., only a lower 70% of the ABL was experimentally modelled. The projected scaled-up ABL thicknesses are 395, 354, and 416 m full-scale in the 1:395, 1:236, and 1:208 simulations, respectively. Experimental results show similar trends in all three configurations not depending on the simulation length scale factor. This clearly indicates a possibility to physically, in the wind tunnel, reproduce the same rural atmospheric airflows at different simulation length scales.     

Turbulence structure and coherent vortices in open-channel flows with wind-induced water waves

When wind-induced water waves appear over the free-surface flows such as natural rivers and artificial channels, large amounts of oxygen gas and heat are transported toward the river bed through the interface between water and wind layers. In contrast, a bed region is a kind of turbulent boundary layer, in which turbulence generation and its transport is promoted by the production of bed shear stress. In particular, coherent hairpin vortices, together with strong ejection events toward the outer part of the layer, promote mass and momentum exchanges between the inner and outer layers. It is inferred that such a near-bed turbulence may be influenced significantly by these air–water interfacial fluctuations accompanied with free-surface velocity shear and wind-induced water waves. However, these wind effects on the wall-turbulence structure are less understood. To address these exciting and challenging topics, we conducted particle imagery velocimetry (PIV) measurements in open-channel flows combined with air flows, and furthermore the present measured data allows us to investigate the effects of air–water interactions on turbulence structure through the whole depth region.     

Stratified flow interactions with a suspended canopy

Field observations of the interactions between a stratified flow and a canopy suspended from the free surface above a solid boundary are described and analysed. Data were recorded in and around the canopy formed by a large long-line mussel farm. The canopy causes a partial blockage of the water flow, reducing velocities in the upper water column. Deceleration of the approaching flow results in a deepening of isopycnals upstream of the canopy. Energy considerations show that the potential for an approaching stratified flow to be diverted beneath a porous canopy is indicated by a densimetric Froude number. Strong stratification or low-velocities inhibit vertical diversion beneath the canopy, instead favouring a horizontal diversion around the sides. The effect on vertical mixing is also considered with a shear layer generated beneath the canopy and turbulence generated from drag within the canopy. In the observations, stratification is shown to be of sufficient strength to limit the effectiveness of the first mixing process, while the turbulence within the canopy is likely to enhance vertical exchange. Velocity and temperature microstructure measurements are used to investigate the effect of the canopy on turbulent dissipation and show that dissipation is enhanced within the canopy.     

A Parameterization of the Roughness Length for the Air-Sea Interface in Free Convection

The response of the upper ocean to the parameterization of the roughness length z ₀ on the air side of the air-sea interface is studied using a one-dimensional mixed-layer model. In particular, it is shown that in the free convection limit when both the wind speed and the friction velocity approach zero, the familiar Charnock formula for the momentum roughness, which relies solely on wind generation, can be modified to account for contributions arising from the thermally generated turbulence. Therefore, a new parameterization is proposed for the momentum roughness length which extends the Charnock formula down to zero friction velocity. The value of a parameter which enters in the new formulation is determined by making use of exsisting free convection surface flux parameterizations. The effect of the new parameterization on the model performance is tested using data from the ocean weathership station Papa (OWS P), and data from the Long-Term Upper-Ocean Study (LOTUS) experiment. Simulations were carried out using a recently developed one-dimensional, second-order, turbulence closure scheme over diurnal as well as seasonal time scales. The findings suggest that the new momentum roughness parameterization improves the overall agreement between the observed and simulated sea-surface temperature (SST).     

The Structure of the Shear Layer in Flows over Rigid and Flexible Canopies

Flume experiments were conducted with rigid and flexible model vegetation to study the structure of coherent vortices (a manifestation of the Kelvin–Helmholtz instability) and vertical transport in shallow vegetated shear flows. The vortex street in a vegetated shear layer creates a pronounced oscillation in the velocity profile, with the velocity near the top of a model canopy varying by a factor of three during vortex passage. In turn, this velocity oscillation drives the coherent waving of flexible canopies. Relative to flows over rigid vegetation, the oscillation in canopy geometry has the effect of decreasing the amount of turbulent vertical momentum transport in the shear layer. Using a waving plant to determine phase in the vortex cycle, each vortex is shown to consist of a strong sweep at its front (during which the canopy is most deflected), followed by a weak ejection at its rear (when the canopy height is at a maximum). Whereas in unobstructed mixing layers the vortices span the entire layer, they encompass only 70% of the flexibly obstructed shear layer studied here.     

Some characteristics of the urban boundary layer above Rome, Italy, and applicability of Monin–Obukhov similarity

Air temperature and wind speed profiles measured during one year by means of a SODAR-RASS system located within a large park were examined for the urban boundary layer (UBL) over Rome, Italy. These data, combined with velocity and temperature measurements performed near the ground were used to analyze the vertical structure of the boundary layer and to estimate some turbulence parameters characterizing the surface layer. About 52,000 vertical profiles of wind speed and temperature were used for the analysis, allowing investigation for a large variety of stability conditions. First, friction velocity and Obukhov length were examined, showing clearly their dependence on the time of day and season. Second, the applicability of the Monin–Obukhov (MO) similarity theory—developed over rural terrain—was tested up to 200?m above ground level. For the wind speed profiles, the performance of the MO similarity degrades with both increasing height and stability, with maximum errors that are on the order of 300?% at 200?m for the most stable case. In contrast, for the air temperature the error always remains below 50?%.     

Numerical Simulation of Interaction of the Heavy Gas Cloud with the Atmospheric Surface Layer

The numerical time-dependent three-dimensional model [Kovalets, I.V. and Maderich, V.S.: 2001, Int. J. Fluid Mech. Res. 30, 410–429] of the heavy gas dispersion in the atmospheric boundary layer has been improved by parameterizing momentum and heat fluxes on the surface of Earth using Monin–Obukhov similarity theory. Three parameterizations of heat exchange with the surface of Earth were considered: (A) formula of Yaglom A.M. and Kader B.A. [1974, J. Fluid Mech. 62, 601–623] for forced convection, (B) interpolation formula for mixed convection and (C) similarity relationship for mixed convection [Kader, B.A. and Yaglom, A.M.: 1990, J. Fluid Mech. 212, 637–662]. Two case studies were considered. In the first study based on experiment of Zhu et al., J. Hazard Mater 62, 161–186], the interaction of an isothermal heavy gas plume with an atmospheric surface layer was simulated. It was found that stable stratification in the cloud essentially suppresses the turbulence in the plume, reducing the turbulent momentum flux by a factor of down to 1/5 in comparison with the undisturbed value. This reduction essentially influences velocities in the atmospheric boundary layer above the cloud, increasing the mean velocity by a factor of up to 1.3 in comparison with the undisturbed value. A simulation of cold heavy gas dispersion was carried out in the second case based on field experiment BURRO 8. It was shown that both forced and free convections under moderate wind speeds significantly influence the plume. The relative rms and bias errors of prediction the plume’s height were σ_H ≈ 30% and ɛ_H = − 10%, respectively, for parameterization B, while for A and C the errors were σ_H ≈ 80% and ɛ_H ≈ − 65%. It is therefore advised to use the simple parameterization B in dense gas dispersion models.     

Wind-induced pressure at a tunnel portal

In order to properly size the mechanical ventilation system of a tunnel, it is essential to estimate the wind-driven pressure difference that might rise between its two portals. In this respect, we explore here the pressure distribution over a tunnel portal under the influence of an incident atmospheric boundary layer and, in particular, its dependency on wind direction and on tunnel geometry. Reduced scale models of generic configurations of a tunnel portal are studied in an atmospheric wind tunnel. Pressure distributions over the front section of different open cavities are measured with surface taps, which allows us to infer the influence of the tunnel aspect ratio and wind direction on a pressure coefficient \(C_{P}\), defined as a spatially and time averaged non-dimensional pressure. Experiments reveal that the magnitude of the coefficient \(C_{P}\), as a function of the wind direction, is significantly influenced by the portal height-to-width ratio and almost insensitive to its length. The experimental data set is completed by hot-wire anemometry measurements providing vertical distribution of velocity statistics. The same configurations are simulated by numerically solving the Reynolds-averaged Navier–Stokes equations, adopting the standard \(k - \varepsilon\) turbulence model. Despite some discrepancies between numerical and experimental estimates of some flow parameters (namely the turbulent kinetic energy field), the numerical estimates of the pressure coefficients \(C_{P}\) show very good agreement with experimental data. The latter is also compared to the predictions of an analytical model, based on the estimate of a spatially averaged velocity within an infinitely long street canyon. The results of the model, which takes into account varying canyon aspect ratios, are in reasonable agreement with experimental data for all cases studied. Notably, its predictions are significantly better than those provided by the simple analytical relations usually adopted as a reference in tunnel ventilation studies.     

Subgrid-Scale Modeling of Reacting Scalar Fluxes in Large-Eddy Simulations of Atmospheric Boundary Layers

In large-eddy simulations of atmospheric boundary layer turbulence, the lumped coefficient in the eddy-diffusion subgrid-scale (SGS) model is known to depend on scale for the case of inert scalars. This scale dependence is predominant near the surface. In this paper, a scale-dependent dynamic SGS model for the turbulent transport of reacting scalars is implemented in large-eddy simulations of a neutral boundary layer. Since the model coefficient is computed dynamically from the dynamics of the resolved scales, the simulations are free from any parameter tuning. A set of chemical cases representative of various turbulent reacting flow regimes is examined. The reactants are involved in a first-order reaction and are injected in the atmospheric boundary layer with a constant and uniform surface flux. Emphasis is placed on studying the combined effects of resolution and chemical regime on the performance of the SGS model. Simulations with the scale-dependent dynamic model yield the expected trends of the coefficients as function of resolution, position in the flow and chemical regime, leading to resolution-independent turbulent reactant fluxes.     

Observations of the cross-wind velocity variance in the stable boundary layer

Nine tower datasets over grassland, brush rangeland, snow covered plain, the ocean, three different pine forests, an aspen forest and an urban site, are used to document the scale-dependence of the cross-wind velocity variance in the stable boundary layer. The turbulence velocity variance scales with the surface momentum flux, as reported in previous studies. Such scaling removes the stability dependence of the variance at a given site, and also removes most of the differences between sites. The scaling is more effective with use of a record-dependent averaging time for defining the turbulent fluctuations. The variable averaging time is the timescale associated with the gap region in the heat flux multiresolution cospectra. On scales larger than turbulence and less than a few hours (mesoscale), variations in the cross-wind velocity variance at a given site are not related to local variables such as the friction velocity. Possible exceptions include suppression of turbulence and mesoscale motions in well-defined drainage flows and enhancement of turbulence and mesoscale motions in stronger winds downstream of a ridge. Larger mesoscale variance is associated with complex terrain and forested sites compared to the more homogeneous sites in flat terrain with short or no vegetation. These differences between sites are related to the absence of a gap region in the velocity spectra at the complex terrain and forested sites. The observed probability distribution functions of the total variance and the mesoscale variance are documented for different averaging times, stability classes and site characteristics.