Conditional analysis of temperature and humidity microfronts and ejection/sweep motions within and above a deciduous forest

A conditional sampling technique using a multilevel scheme was applied to the detection of temperature and humidity microfronts and organized ejection/sweep motions under different atmospheric stabilities. Data were obtained with seven triaxial sonic anemometer/thermometers and three Lyman-alpha hygrometers within and above a deciduous forest. Both temperature and humidity microfronts were identified in unstable cases, but only humidity microfronts could be detected under neutral conditions. Inverted temperature ramps occurred under slightly stable conditions. Occasionally, wave-like patterns appeared within the canopy, seemingly coupled with inverse ramps occurring above the forest. The frequency of occurrence of scalar microfronts appears to have no clear dependence on atmospheric stability, and averages 74–84 s per cycle with a mode of about 50 s per cycle. However, the strength of ejections and sweeps, shown by the vertical velocity averaged within structures, was reduced by increasing atmospheric stability. Structures identified under different stabilities show many similarities in their patterns of scalar ramps, and associated velocity and surface pressure. Profiles of short-term averaged longitudinal velocity at different times during the microfront passage show that the air within the canopy was retarded and an intensified shear above the canopy occurred prior to the passage of the microfront. Results from the present conditional analysis strongly suggest an important role of shear instability in the formation of canopy coherent structure.     

Estimating Heat Sources And Fluxes In Thermally Stratified Canopy Flows Using Higher-Order Closure Models

Over the past two decades, several inverse methods have been proposed to estimatescalar source and sink strengths from measured mean concentration profiles withinthe canopy volume (hereafter termed the `inverse' problem). These inverse methodscommonly assumed neutral atmospheric stability conditions for the entire canopyvolume. For non-neutral conditions, atmospheric stability corrections in inverseschemes were limited to adjusting the integral time scale or other flow statistics tomatch well-established surface-layer similarity relations above the canopy. Suchstability corrections do not explicitly consider the local stability effects within thecanopy volume. Currently, there is no satisfactory inverse scheme that explicitlyaccounts for local atmospheric stability for canopy turbulence. A Eulerian inversemethod that explicitly accounts for local atmospheric stability within the canopy isdeveloped using second-order closure principles. Field testing the method is conductedusing temperature measurements from two field experiments collected in an even-ageduniform loblolly pine forest. It is demonstrated that by accounting for local atmospheric stability in the inversion scheme, the agreement between modelled sensible heat flux calculations and measurements improve by 60% for stable conditions, 10% for near-neutral conditions and 20% for unstable conditions     

Stability Dependence of Canopy Flows over a Flat Larch Forest

The dependence on atmospheric stability of flow characteristics adjacent to a very rough surface was investigated in a larch forest in Japan. Micrometeorological measurements of three-dimensional wind velocity and air temperature were taken at two heights above the forest, namely 1.7 and 1.2 times the mean canopy height h. Under near-neutral and stable conditions, the observed turbulence statistics suggest that the flow was likely to be that of the atmospheric surface layer (ASL) at 1.7h, and of the roughness sublayer (RSL) at 1.2h. However, in turbulence spectra, canopy-induced large coherent motions appeared clearly at both heights. Even under strongly stable conditions, the large-scale motions were retained at 1.2h, whereas they were overwhelmed by small-scale motions at 1.7h. This phenomenon was probably due to the enhanced contribution of the ASL turbulence associated with nocturnal decay of the RSL depth, because the small-scale motions appeared at frequencies close to the peak frequencies of well-known ASL spectra. This result supports the relatively recent concept that canopy flow is a superimposition of coherent motions and the ASL turbulence. The large-scale motions were retained in temperature spectra over a wider region of stability compared to streamwise wind spectra, suggesting that a canopy effect extended higher up for temperature than wind. The streamwise spacing of dominant eddies according to the plane mixing-layer analogy was only valid in a narrow range at near neutral, and it was stabilised at nearly half its value under stable conditions.     

An Analytical one-Dimensional Second-Order Closure Model of Turbulence Statistics and the Lagrangian Time Scale Within and Above Plant Canopies of Arbitrary Structure

An analytical one-dimensional second-order closure model is developed to describe the within canopy velocity variances, turbulent intensities, dissipation rates, Lagrangian time scale and Lagrangian far field diffusivities for vegetation canopies of arbitrary structure and density. The model incorporates and extends the model of momentum transfer developed by Massman (1997) and the model of within canopy velocity variances developed by Weil (unpublished) from the second-order closure model of Wilson and Shaw (1977). Model predictions of within and above canopy velocity variances, turbulent intensities, dissipation rates and the Lagrangian time scale are in reasonable agreement with previously measured or estimated values for these parameters. The present model suggests that the Lagrangian time scale and the far field diffusivity could be strongly dependent upon foliage structure and density through the foliage effects on the velocity variances. A simple formulation for the Lagrangian time scale at canopy height is derived from model results. Taken as a whole, the present model may provide a relatively simple way to incorporate turbulence parameters into models of soil/canopy/atmosphere mass transfer.     

Turbulence Statistics Measurements in a Northern Hardwood Forest

Tower-based turbulence measurements were collected in and over a mixed hardwood forest at the University of Michigan BiologicalStation (UMBS) UMBSflux site in the northern summerof 2000. Velocity and temperature fluctuations were measured at five levels within the canopy (up to the canopy height, H = 21.4 m), using one- and three-dimensional sonic anemometers and fine-wire thermocouples. Six additional thermocouples were distributed over the canopy-layer depth. Three-dimensional velocities and sonic temperatures were also measured above the canopy at 1.6H and at 2.15H on the AmeriFlux tower located at the UMBSflux site. Vertical profiles of buoyancy flux, mean horizontal velocity, Reynolds stress, and standard deviation and skewness of velocity components were calculated. The analysis of these measurements aims at a multi-layer parameterization framework of turbulence statistics forimplementation in Lagrangian stochastic models. Turbulence profiles and power spectra above the canopy were analyzed in the context of Monin-Obukhov similarity theory (MOST) and Kolmogorov theory, as determined by stability at the top level (2.15H), to assess the extent to which surface scaling is valid as the canopy top is approached. Velocity spectra were computed to explore the potential of estimating the viscous dissipation rate, and results show that the high frequency range of the spectra above the canopy exhibits the roll-off predicted by Kolmogorov theory. Similarly, velocity standard deviations above the canopy converge to MOST predicted values toward the top level, and spectral peaks shift with stability, as expected. Within the canopy, both turbulence statistics profiles and spectral distributions follow the general known characteristics inside forests.     

Flow over high roughness elements

The mean velocity and longitudinal turbulence-intensity distributions inside the zone of and above high roughness elements were investigated experimentally. This was accomplished by using a model forest canopy. The results indicate that the flow may be divided into transition and fully-developed flow regions, followed by a short adjustment distance near the downstream terminus of the rough boundary. The transition region has a strong effect on the flow characteristics within and above the layer of roughness elements. Generally, a similar qualitative variation for both velocity and turbulence was found inside and above the roughness zone, whose influence extends to more than three times the roughness height.Investigation of the modified universal logarithmic law for describing the velocity variation above the roughness zone revealed that both of the so-called similarity parameters, i.e., friction velocity and roughness length, are not local constants. On the contrary, for a given flow and local conditions they vary drastically with height. It is suspected that this is due to the fact that the classical assumption of constant shear stress throughout the boundary layer or significant portions of it is not satisfied in the case of roughness elements many times greater in height than the thickness of the viscous wall zone.     

Turbulence Spectra And Dissipation Rates Above And Within A Forest Canopy

Three velocity componentsand temperature were measured usingthree-dimensional sonic anemometers/thermometers attwo levels, above and within a forest canopy, in theChangbai Mountains of northeast China. Turbulencespectral structure, local isotropy anddissipation rates above and within the forest canopywere calculated using the eddy correlation method.Results show that the normalized turbulent spectralcurves have -2/3 slopes in the inertial subrange.While the shapes of the spectra are in good agreementwith the Kansas flat terrain results, the atmosphericturbulence is anisotropic above the forest canopy. Dueto breaking down of large eddies by the foliage,branches and trunks, the spectral peak frequencies forvelocity and temperature are higher withinthan above the forest canopy. Compared withmeasurements from previous studies over flat terrain,the velocity and temperature spectra above andinside the forest canopy appear to shift toward higherfrequencies. The turbulence is approximately isotropicin the inertial subrange within the forest canopy, and isanisotropic above the forest canopy. The turbulentkinetic energy and heat energy dissipation rates aboveand inside the forest canopy are much larger thanthose obtained by Kaimal and Hogstrom over grasslandand grazing land. The distinct features in the resultsof the present experiment may be attributed to thedynamic forcing caused by the rough surface of the forestcanopy.     

Turbulent Flux Measurements Above and Below the Overstory of a Boreal Aspen Forest

Turbulent flux measurements both above and beneath the canopy of a boreal aspen forest are described. Velocity skewness showed that, beneath the aspen canopy, turbulence was dominated by intermittent, downward penetrating gusts. Eulerian horizontal length scales calculated from integration of the autocorrelation function or spectral peaks were 9.0 and 1.4 times the mean aspen height of 21.5 m respectively. Above-canopy power spectral slopes for all velocity components followed the -2/3 power law, whereas beneath-canopy slopes were closer to -1 and showed a spectral short cut in the horizontal and vertical components. Cospectral patterns were similar both above and beneath the canopy. The Monin–Obukhov similarity function for the vertical wind velocity variance was a well-defined function of atmospheric stability, both above and beneath the canopy. Nocturnal flux underestimation and departures of this similarity function from that expected from Monin–Obukhov theory were a function of friction velocity. Energy balance closure greater than 80% was achieved at friction velocities greater than 0.30 and 0.10 m s-1, above and below the aspen canopy, respectively. Recalculating the latent heat flux using various averaging periods revealed a minimum of 15 min were required to capture 90% of the 30-min flux. Linear detrending reduced the flux at shorter averaging periods compared to block averaging. Lack of energy balance closure and erratic flux behaviour led to the recalculation of the latent and sensible heat fluxes using the ratio of net radiation to the sum of the energy balance terms.     

Bulk Transfer Relations for the Roughness Sublayer

In the roughness sublayer (RSL), Monin–Obukhov surface layer similarity theory fails. This is problematic for atmospheric modelling applications over domains that include rough terrain such as forests or cities, since in these situations numerical models often have the lowest model level located within the RSL. Based on empirical RSL profile functions for momentum and scalar quantities, and scaling the height with the RSL height z _*, we derive a simple bulk transfer relation that accounts for RSL effects. To verify the validity of our approach, these relations are employed together with wind speed and temperature profiles measured over boreal forest during the BOREAS experimental campaign to estimate momentum and heat fluxes. It is demonstrated that, when compared with observed flux values, the inclusion of RSL effects in the transfer relations yields a considerable improvement in the estimated fluxes.     

Influence of foliar density and thermal stability on profiles of Reynolds stress and turbulence intensity in a deciduous forest

Observations were made of turbulence in an extensive deciduous forest on level terrain using a vertical array of seven three-dimensional sonic anemometer/thermometers within and above the canopy. Data were collected through the period of leaf fall and over a range of thermal stabilities. A bulk canopy drag coefficient was nearly independent of the density of the forest but decreased greatly with the onset of nocturnal stability. The depth of penetration of momentum into the forest increased with leaf fall but, again, was greatly curtailed by stable conditions. Turbulent velocities decreased with increasing depth in the forest but relative turbulence intensities increased to mid-canopy levels. Leaf density influenced turbulence levels but not as strongly as did thermal stability. Thermal effects were adequately described by the single parameter h/L, where h is the canopy height and L is the Monin-Obukhov length. The longitudinal and vertical velocity correlation coefficient was larger in magnitude than expected in the upper layers of the forest but decreased to a small value in the lowest layers where the Reynolds stress was small. The ratio _w /u ^*, where u ^* is the local friction velocity, reflected changes in the uw correlation, becoming smaller than usual in the upper canopy layers. It is believed that these effects result from the intermittent, spatially coherent structures that are responsible for a large fraction of the momentum flux to the forest.     

Are Urban-Canopy Velocity Profiles Exponential?

Using analyses of data from extant direct numerical simulations and large-eddy simulations of boundary-layer and channel flows over and within urban-type canopies, sectional drag forces, Reynolds and dispersive shear stresses are examined for a range of roughness densities. Using the spatially-averaged mean velocity profiles these quantities allow deduction of the canopy mixing length and sectional drag coefficient. It is shown that the common assumptions about the behaviour of these quantities, needed to produce an analytical model for the canopy velocity profile, are usually invalid, in contrast to what is found in typical vegetative (e.g. forest) canopies. The consequence is that an exponential shape of the spatially-averaged mean velocity profile within the canopy cannot normally be expected, as indeed the data demonstrate. Nonetheless, recent canopy models that allow prediction of the roughness length appropriate for the inertial layer’s logarithmic profile above the canopy do not seem to depend crucially on their (invalid) assumption of an exponential profile within the canopy.     

Spatial Variability of Both Turbulent Fluxes and Temperature Profiles in an Urban Roughness Layer

The spatial variability of both turbulent flow statistics in the roughness sublayer (RSL) and temperature profiles within and above the canopy layer (CL) were investigated experimentally in a densely built-up residential area in Tokyo, Japan. Using five towers with measuring devices, each tower isolated from the others by at least 200 m, we collected high-frequency measurements of velocity and temperature at a height z=1.8 z _H, where z _H, the mean building height in the area, is 7.3 m. Also, temperature profiles were measured from z=0.4 to 1.8 z _H. The ‘areal mean’ geometric parameters that were obtained for the areas within 200 m of each tower were fairly homogeneous among the tower sites. The main results are as follows: (1) The spatial variability of all RSL turbulent statistics, except the sensible heat flux, was comparable to that reported in a pine forest. Also, the variability decreased with increasing friction velocity. (2) The spatial variability of the RSL sensible heat flux was larger than that reported in a pine forest. Also, the variability depended on the time of the day and became larger in the morning. The difference among the sites was well related to the areal fraction of vegetation. (3) The spatial variability of the CL temperature profile depended on the time of the day and became larger in the morning. Nevertheless, the spatial standard deviation of CL temperature was always below 0.7 K. (4) It is suggested that the “warming-up” process in the morning when heat storage is dominant increases the spatial variation of RSL sensible heat flux and CL temperature according to the local properties around each tower and the variation decreases once there is further convective mixing in the midday     

The sensitivity of the water balance of a wet multilayer model pine canopy to variations in meteorological input

A multilayer canopy model of a pine forest is used to investigate the sensitivity of the water balance of the wet canopy to variations in meteorological input. The multilayer model does not take into account large-scale eddies, which are now considered to be of importance to canopy transport. It does, however, provide realistic simulations of wet canopy water balance and often predicts interception loss rates higher than those predicted by a unilayer model for the same meteorological input. Stable layers both within and above the canopy are often simulated during rainfall events, and these may help to spontaneously generate large-scale eddies or waves within forest canopies. The sensitivity study for a wet canopy suggests that low vapour pressure deficits and low wind speeds are associated with unstable surface conditions, and increasing values of both variables are associated with decreasing canopy drainage values and increasing evaporative losses. Low short- or long-wave radiation inputs are associated with stable surface conditions, and increasing values of both variables are associated with decreasing canopy drainage values and increasing evaporative losses. Increasing temperature is associated with increasing surface stability and increasing canopy drainage and decreasing evaporative losses. In real situations the tendency for increasing temperature to cause surface stability and decreased evaporative loss is probably compensated by the opposite effects of increasing short- or long-wave radiation. The model simulations suggest that wet forest canopies may be better ventilated at low temperatures, if other meteorological conditions are constant.     

The influence of buoyancy on third-order turbulent velocity statistics within a deciduous forest

The influence of atmospheric stability on the behaviour of the third moment of flow velocities observed inside a deciduous forest canopy is examined. Results suggest that buoyancy plays a dominant role in dictating the magnitude of gusts observed inside tall vegetation. Furthermore, an examination of the turbulence recorded throughout leaf fall inside the same forest indicates that larger velocity skewnesses are observed inside a canopy in full leaf than inside a sparse canopy. The behaviour of the measured terms in the non-dimensionalized rate equation of the third moment of canopy flow velocities is also examined. Turbulent diffusion and turbulence gradient interaction terms are largest in stable conditions in the upper canopy layer while these are most important in unstable conditions in the lower canopy layer. In all stability regimes, the turbulent diffusion term is the main source of skewness. The turbulence gradient interaction term, the residual and buoyant production terms all contribute to destroy skewness in stable conditions.     

Measurements and Langevin simulations of mean tracer concentration fields downwind from a circular line source inside an alfalfa canopy

Two Langevin simulations of trajectories of marked fluid elements in inhomogenous turbulence, where the Lagrangian length and vertical velocity scales are height dependent, were compared with field data. A CO₂ tracer was released from a circular line source and the concentration profiles were measured for diffusion distances of 50 and 100 cm inside and above an alfalfa canopy.One of the simulations, suggested by Wilson et al. (1983), biases the vertical velocities by adding a mean upward drift. The second simulation proposed here by-passes this difficulty by reflecting marked particles according to a probability calculated from the gradient in vertical velocity variance between the beginning and the end of each step. This simulation also makes use of a constant time-scale within the canopy, following preliminary results from a turbulence experiment within a forest (Leclerc, 1987).Comparing the results of these simulations with the field data shows that the simulation proposed by Wilson et al. (1983) does not correctly reproduce the difusion for the larger fetch in systems exhibiting strong gradients in vertical velocity variance. Instead, the modelled plumes exhibit a bulge at the source height whereas the field data show smooth profiles. In addition, the modelled plumes overestimate the vertical spread of the plumes, which is possibly due to the inadequacy of the approach in severely inhomogeneous systems. In contrast, the results from the tracer experiments indicate that the diffusion can be better reproduced with the use of a reflection probability calculated at each step. The discrepancies between the experimental results and the simulation using a reflection probability are attributed to stability effects.     

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.     

Influence of Source/Sink Distributions on Flux–Gradient Relationships in the Roughness Sublayer Over an Open Forest Canopy Under Unstable Conditions

The flux–gradient relationships in the unstable roughness sublayer (RSL) over an open canopy of black spruce forest were examined using long-term observations from an instrumented tower. The observed gradients normalised with the surface fluxes and height above the zero-plane displacement showed differences from a universal function established in the surface layer. The magnitude of differences was not constant throughout the year even at the same observation height. Also the magnitude of the differences was different for each scalar, and scalar similarity in the context of the flux–gradient relationship did not always hold. The variation of the differences was explained by the relative contribution of overstorey vegetation to the total flux from the entire ecosystem. This suggests that a mismatch of the vertical source/sink distributions between scalars leads to a different strength of the near-field dispersion effect for each scalar, and this resulted in inequality of eddy diffusivity among scalars in the RSL. An empirical method that predicts the magnitude of differences is proposed. With this method, it is possible to estimate the eddy diffusivity of scalars provided that the relative contribution of overstorey vegetation to the total flux from the ecosystem is known. Also this method can be used to estimate the eddy diffusivity for scalars whose primary sources are at ground level, such as methane and nitrous oxide.     

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.     

Estimation of Scalar Source/Sink Distributions in Plant Canopies Using Lagrangian Dispersion Analysis: Corrections for Atmospheric Stability and Comparison with a Multilayer Canopy Model

Source/sink distributions of heat, water vapour andCO₂ within a rice canopy were inferred using aninverse Lagrangian dispersion analysis and measuredmean profiles of temperature, specific humidity andCO₂ mixing ratio. Monin–Obukhov similarity theorywas used to account for the effects of atmosphericstability on _w(z), the standard deviation ofvertical velocity and _L(z), the Lagrangian timescale of the turbulence. Classical surface layer scaling was applied in the inertial sublayer (z > z_ruf)using the similarity parameter = (z - d)/L, where z is height above ground, d is the zero plane displacementheight for momentum, L is the Obukhov length,and z_ruf 2.3h_c, where h_c iscanopy height. A single length scale h_c, was usedfor the stability parameter 3 = h_c/L in the height range 0.25 < z/h_c < 2.5. This choice is justified by mixing layer theory, which shows that within the roughness sublayer there is one dominant turbulence length scaledetermined by the degree of inflection in the windprofile at the canopy top. In the absence of theoretical or experimental evidence for guidance,standard Monin–Obukhov similarity functions, with = h_c/L, were used to calculate the stabilitydependence of _w(z) and _L(z) in the roughness sublayer. For z/h_c < 0.25 the turbulence length and time scales are influenced by the presence of the lowersurface, and stability effects are minimal. With theseassumptions there was excellent agreement between eddycovariance flux measurements and deductions from theinverse Lagrangian analysis. Stability correctionswere particularly necessary for night time fluxes whenthe atmosphere was stably stratified.The inverse Lagrangian analysis provides a useful toolfor testing and refining multilayer canopy models usedto predict radiation absorption, energy partitioningand CO₂ exchanges within the canopy and at thesoil surface. Comparison of model predictions withsource strengths deduced from the inverse analysisgave good results. Observed discrepancies may be dueto incorrect specification of the turbulent timescales and vertical velocity fluctuations close to theground. Further investigation of turbulencecharacteristics within plant canopies is required toresolve these issues.     

Analysis of Low-Frequency Turbulence Above Tall Vegetation Using a Doppler Sodar

This study applies acoustic sounding to observe coherent structures in the roughness sublayer (RSL) above tall vegetated surfaces. Data were collected on 22 days during two separate field experiments in summer 2003. A quality control scheme was developed to ensure high data quality of the collected time series. The data analysis was done using both discrete and continuous wavelet transform. The flow in the RSL was found to be a superposition of dynamic Kelvin–Helmholtz instabilities and convective mixing. The characteristic time scales for coherent structures resulting from the dynamic instabilities were observed to be approximately 20–30 s while thermal eddies have much larger time scales of 190–210 s. The degree of vertical coherency in the RSL increases with the flow evolving from neutral to near-convective conditions. This increase in the degree of organisation is attributed to the evolution of attached thermal eddies. The coherent structures resulting from instabilities were found to be present throughout the RSL but do not contribute to the increased vertical coherency. An alternative conceptual approach for the definition of the RSL is proposed, which yields its maximum vertical extent to five times the canopy height.