A lagrangian random-walk model for simulating water vapor,CO2 and sensible heat flux densities and scalar profiles over and within a soybean canopy

An integrated canopy micrometeorological model is described for calculating CO₂, water vapor and sensible heat exchange rates and scalar concentration profiles over and within a crop canopy. The integrated model employs a Lagrangian random walk algorithm to calculate turbulent diffusion. The integrated model extends previous Lagrangian modelling efforts by employing biochemical, physiological and micrometeorological principles to evaluate vegetative sources and sinks. Model simulations of water vapor, CO₂ and sensible heat flux densities are tested against measurements made over a soybean canopy, while calculations of scalar profiles are tested against measurements made above and within the canopy. The model simulates energy and mass fluxes and scalar profiles above the canopy successfully. On the other hand, model calculations of scalar profiles inside the canopy do not match measurements.The tested Lagrangian model is also used to evaluate simpler modelling schemes, as needed for regional and global applications. Simple, half-order closure modelling schemes (which assume a constant scalar profile in the canopy) do not yield large errors in the computation of latent heat (LE) and CO₂ (F _c ) flux densities. Small errors occur because the source-sink formulation of LE andF _c are relatively insensitive to changes in scalar concentrations and the scalar gradients are small. On the other hand, complicated modelling frames may be needed to calculate sensible heat flux densities; the source-sink formulation of sensible heat is closely coupled to the within-canopy air temperature profile.     

Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

Modelling Vegetation-Atmosphere Co2 Exchange By A Coupled Eulerian-Langrangian Approach

A Eulerian-Lagrangian canopy microclimate model wasdeveloped with the aim of discerning physical frombiophysical controls of CO₂ and H₂O fluxes. The model couples radiation attenuation with mass,energy, and momentum exchange at different canopylevels. A unique feature of the model is its abilityto combine higher order Eulerian closure approachesthat compute velocity statistics with Lagrangianscalar dispersion approaches within the canopy volume. Explicit accounting for within-canopy CO₂,H₂O, and heat storage is resolved by consideringnon-steadiness in mean scalar concentration andtemperature. A seven-day experiment was conducted inAugust 1998 to investigate whether the proposedmodel can reproduce temporal evolution of scalar(CO₂, H₂O and heat) fluxes, sources andsinks, and concentration profiles within and above auniform 15-year old pine forest. The modelreproduced well the measured depth-averaged canopy surfacetemperature, CO₂ and H₂O concentrationprofiles within the canopy volume, CO₂ storageflux, net radiation above the canopy, and heat andmass fluxes above the canopy, as well as the velocitystatistics near the canopy-atmosphere interface. Implications for scaling measured leaf-levelbiophysical functions to ecosystem scale are alsodiscussed.     

The Turbulent Lagrangian Time Scale in Forest Canopies Constrained by Fluxes,Concentrations and Source Distributions

One-dimensional Lagrangian dispersion models, frequently used to relate in-canopy source/sink distributions of energy, water and trace gases to vertical concentration profiles, require estimates of the standard deviation of the vertical wind speed, which can be measured, and the Lagrangian time scale, T _L , which cannot. In this work we use non-linear parameter estimation to determine the vertical profile of the Lagrangian time scale that simultaneously optimises agreement between modelled and measured vertical profiles of temperature, water vapour and carbon dioxide concentrations within a 40-m tall temperate Eucalyptus forest in south-eastern Australia. Modelled temperature and concentration profiles are generated using Lagrangian dispersion theory combined with source/sink distributions of sensible heat, water vapour and CO₂. These distributions are derived from a multilayer Soil-Vegetation-Atmospheric-Transfer model subject to multiple constraints: (1) daytime eddy flux measurements of sensible heat, latent heat, and CO₂ above the canopy, (2) in-canopy lidar measurements of leaf area density distribution, and (3) chamber measurements of CO₂ ground fluxes. The resulting estimate of Lagrangian time scale within the canopy under near-neutral conditions is about 1.7 times higher than previous estimates and decreases towards zero at the ground. It represents an advance over previous estimates of T _L , which are largely unconstrained by measurements.     

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.     

The Effect of the Vertical Source Distribution on Scalar Statistics within and above a Forest Canopy

Little is known about in-canopy processes that may alter forest–atmosphere exchanges of trace gases and aerosols. To improve our understanding of in-canopy mixing, we use large-eddy simulation to study the effect of scalar source/sink distributions on scalar concentration moments, fluxes, and correlation coefficients within and above an ideal forest canopy. Scalars are emitted from: (1) the ground, (2) the canopy, and (3) both the ground and the canopy; a scalar is also deposited onto the canopy. All scalar concentration moments, fluxes, and correlation coefficients are affected by the source location/distribution, as is the scalar segregation intensity. We conclude that vertical source/sink distribution has a profound impact on scalar concentration profiles, fluxes, correlation coefficient, and scalar segregation.     

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     

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.     

Buoyancy and The Sensible Heat Flux Budget Within Dense Canopies

In contrast to atmospheric surface-layer (ASL) turbulence, a linear relationship between turbulent heat fluxes (F_T) and vertical gradients of mean air temperature within canopies is frustrated by numerous factors, including local variation in heat sources and sinks and large-scale eddy motion whose signature is often linked with the ejection-sweep cycle. Furthermore, how atmospheric stability modifies such a relationship remains poorly understood, especially in stable canopy flows. To date, no explicit model exists for relating F_T to the mean air temperature gradient, buoyancy, and the statistical properties of the ejection-sweep cycle within the canopy volume. Using third-order cumulant expansion methods (CEM) and the heat flux budget equation, a “diagnostic” analytical relationship that links ejections and sweeps and the sensible heat flux for a wide range of atmospheric stability classes is derived. Closure model assumptions that relate scalar dissipation rates with sensible heat flux, and the validity of CEM in linking ejections and sweeps with the triple scalar-velocity correlations, were tested for a mixed hardwood forest in Lavarone, Italy. We showed that when the heat sources (S_T) and F_T have the same sign (i.e. the canopy is heating and sensible heat flux is positive), sweeps dominate the sensible heat flux. Conversely, if S_T and F_T are opposite in sign, standard gradient-diffusion closure model predict that ejections must dominate the sensible heat flux.     

Footprints and Fetches for Fluxes over Forest Canopies with Varying Structure and Density

A stochastic trajectory model was used to estimate scalar fluxfootprints in neutral stabilityfor canopies of varying leaf area distributions andleaf area indices. An analytical second-order closure model wasused to predict mean wind speed, second moments and the dissipationrate of turbulent kinetic energy within a forest canopy.The influence of source vertical profile on the flux footprint wasexamined. The fetch is longer for surface sourcesthan for sources at higher levels in the canopy. In order tomeasure all the flux components, and thus the total flux, with adesired accuracy, sources were located at the forest floor in thefootprint function estimation. The footprint functions werecalculated for five observation levels above the canopy top. Itwas found that at low observation heights both canopy density andcanopy structure affect the fetch. The higher abovethe canopy top the flux is measured, the more pronounced is the effectof the canopy structure. The forest fetch for flux measurements isstrongly dependent on the required accuracy: The 90% flux fetchis greater by a factor of two or more compared to the 75% fetch. Theupwind distance contributing 75% of flux is as large as 45 timesthe difference between canopy height and the observation heightabove the canopy top, being even larger for low observationlevels.     

Coherent eddies and temperature structure functions for three contrasting surfaces. Part II: Renewal model for sensible heat flux

Sensible heat, latent heat, and other scalar fluxes cannot be measuredwithin short dense canopies, e.g., straw mulches, with standard approachessuch as eddy correlation, Bowen ratio-energy balance, aerodynamic, andvariance methods. However, recently developed surface renewal models, thatare based on the fact that most of the turbulent transfer within and abovecanopies is associated with large-scale coherent eddies, which are evidentas ramp patterns in scalar time series, offer a feasible solution. Wepresent a new air renewal model that calculates sensible heat flux atdifferent heights within and above a canopy from the average cubictemperature structure function, sampled at a moderate rate, and measuredaverage friction velocity. The model is calibrated and tested with datameasured above and within a Douglas-fir forest and above a straw mulch andbare soil. We show that the model describes half-hour variations ofsensible heat flux very well, both within the canopy and roughnesssublayers and in the inertial sublayer, for stable and unstable atmosphericconditions. The combined empirical coefficient that appears in the modelhas an apparently universal value of about 0.4 for all surfaces andheights, which makes application of the model particularly simple. Themodel is used to predict daytime and nighttime sensible heat flux profileswithin the straw mulch and within a small bare opening in the mulch.     

Dispersion of a Passive Scalar Within and Above an Urban Street Network

The transport of a passive scalar from a continuous point-source release in an urban street network is studied using direct numerical simulation (DNS). Dispersion through the network is characterized by evaluating horizontal fluxes of scalar within and above the urban canopy and vertical exchange fluxes through the canopy top. The relative magnitude and balance of these fluxes are used to distinguish three different regions relative to the source location: a near-field region, a transition region and a far-field region. The partitioning of each of these fluxes into mean and turbulent parts is computed. It is shown that within the canopy the horizontal turbulent flux in the street network is small, whereas above the canopy it comprises a significant fraction of the total flux. Vertical fluxes through the canopy top are predominantly turbulent. The mean and turbulent fluxes are respectively parametrized in terms of an advection velocity and a detrainment velocity and the parametrization incorporated into a simple box-network model. The model treats the coupled dispersion problem within and above the street network in a unified way and predictions of mean concentrations compare well with the DNS data. This demonstrates the usefulness of the box-network approach for process studies and interpretation of results from more detailed numerical simulations.     

Statistics of atmospheric turbulence within and above a corn canopy

Two three-dimensional split-film anemometers were used to measure turbulence statistics within and above a corn canopy. Normalised profiles of mean windspeed, root-mean-square velocity, momentum flux, and heat flux were constructed from half-hourly averages by dividing within-canopy measurements by the simultaneous canopy-top measurement. With the exception of the heat flux, these profiles showed consistent shape from day to day. Time series of the three velocity components were recorded on magnetic tape and subsequently analysed to obtain Eulerian time and length scales and the power spectrum of each component at several heights. The timescale was found to have a local minimum value at the top of the canopy. However the length scale L _wformed from the timescale and the root-mean-square vertical velocity varied with height as L _w 0.1 z. The power-spectra were non-dimensionalised to facilitate comparison of spectra at different heights and times. All spectra had -5/3 regions spanning at least two decades in frequency.     

Numerical analysis of flux footprints for different landscapes

Summary A model for the canopy – planetary boundary layer flow and scalar transport based on E- closure was applied to estimate footprint for CO₂ fluxes over different inhomogeneous landscapes. Hypothetical heterogeneous vegetation patterns – forest with clear-cuts as well as hypothetical heterogeneous relief – a bell-shaped valley and a ridge covered by forest were considered. The distortions of airflow caused by these heterogeneities are shown – the upwind deceleration of the flow at the ridge foot and above valley, acceleration at the crest and the flow separation with the reversed flow pattern at lee slopes of ridge and valley. The disturbances induce changes in scalar flux fields within the atmospheric surface layer comparing to fluxes for homogeneous conditions: at a fixed height the fluxes vary as a function of distance to disturbance. Correspondingly, the flux footprint estimated from model data depends on the location of the point of interest (flux measurement point) and may significantly deviate from that for a flat terrain. It is shown that proposed method could be used for the choice of optimal sensor position for flux measurements over complex terrain as well as for the interpretation of data for existing measurement sites. To illustrate the latter the method was applied for experimental site in Solling, Germany, taking into account the complex topography and vegetation heterogeneities. Results show that in certain situations (summer, neutral stratification, south or north wind) and for a certain sensor location the assumptions of idealized air flow structure could be used for measurement interpretation at this site, though in general, extreme caution should be applied when analytical footprint models are used in the interpretation of flux measurements over complex sites.     

Horizontal and Vertical Co2 Advection In A Sloping Forest

A system measuring the horizontal and vertical advection was devised and installed in a sloping forest at the Vielsalm site, Belgium. The measurements showed that under stable conditions a flow regime established below the canopy: air flowed horizontally along the slope and entrained the air above the canopy vertically. This movement occurs during stable nights characterised by strongly negative net radiation. It creates negative air concentration gradients in both the vertical and horizontal directions. The advection fluxes associated with these movements are opposite and of a similar order of magnitude. This implies that the horizontal advection cannot be ignored in the carbon budget equation at night. Unfortunately, the large variability of, and considerable uncertainty about, advection fluxes does not enable one to produce estimates of the source term from these equations. Advection measurement systems should be improved in order to enable such estimates to be made. Particular attention should be paid to the estimation of the vertical velocity above the canopy and to the vertical profiles of the horizontal velocity and horizontal CO₂ gradient below the canopy.     

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 ).     

Dawn‐to‐dusk evolution of air turbulence,temperature and sensible and latent heat fluxes above a forest canopy: Concepts,model and field comparisons

Abstract

Dawn‐to‐dusk evolution of air turbulence, sensible heat and latent heat above a forest during cloud‐free or near‐cloud‐free summer conditions is modelled by way of a system of differential equations. Temperatures in and above the canopy, near canopy‐top wind velocities, early morning leaf moisture (dew) and afternoon canopy ventilation (i.e. heat released from the canopy and from below the canopy) are included in the mathematical treatment. Computed results are compared with field data for atmospheric temperature and wind speed profiles up to 1200 m, within‐canopy temperature, and canopy‐level radiation, turbulent fluxes and wind speeds. Data were collected at a central New Brunswick mixed‐wood forest site dominated by spruce (Picea spp. ) and shade‐tolerant hardwoods for four representative summer days. It was found that the effective canopy temperature was not only affected by insolation, but also by the extent of canopy ventilation and the amount of dew on the foliage. The growth of the mixing layer was affected by canopy ventilation and by above‐canopy wind speeds. Model calculations closely simulated the meteorological observations.     

Estimating a Lagrangian Length Scale Using Measurements of CO2 in a Plant Canopy

Analytical Lagrangian equations capable of predicting concentration profiles from known source distributions offer the opportunity to calculate source/sink distributions through inverted forms of these equations. Inverse analytical Lagrangian equations provide a practical means of estimating source profiles using concentration and turbulence measurements. Uncertainty concerning estimates of the essentially immeasurable Lagrangian length scale ( ${\mathcal{L}}$ ), a key input, impedes the operational practicality of this method. The present study evaluates ${\mathcal{L}}$ within a corn canopy by using field measurements to constrain an analytical Lagrangian equation. Measurements of net CO₂ flux, soil-to-atmosphere CO₂ flux, and in-canopy profiles of CO₂ concentration provided the information required to solve for ${\mathcal{L}}$ in a global optimization algorithm for 30-min time intervals. For days when the canopy was a strong CO₂ sink, the optimization frequently located ${\mathcal{L}}$ profiles that follow a convex shape. A constrained optimization then fit the profile shape to a smooth sigmoidal equation. Inputting the optimized ${\mathcal{L}}$ profiles in the forward and inverse Lagrangian equations leads to strong correlations between measured and calculated concentrations and fluxes. Coefficients of the sigmoidal equation were specific to each 30-min period and did not scale with any measured variable. Plausible looking ${\mathcal{L}}$ profiles were associated with negative bulk Richardson number values. Once the canopy senesced, a simple eddy diffusivity profile sufficed to relate concentrations and sources in the analytical Lagrangian equations.     

Comparison of turbulence statistics within three boreal forest canopies

Three-dimensional sonic anemometers were used to measure velocities and temperatures within three natural boreal forest canopies. Vertical profiles of atmospheric turbulence statistics for a black spruce forest, a jack pine forest, and a trembling aspen forest, all located in southeastern Manitoba, were plotted and compared. The canopy structures were quite different, with total leaf-area indices of 2, 4 and 10, for the pine, aspen, and spruce forests, respectively.The profiles of the first and second moments differed among the canopies, where velocities decreased more rapidly in the top portions of the denser canopies. The velocity distributions were skewed and kurtotic within all canopies, and showed some differences among the canopies. Eulerian time scale profiles were generally similar among the canopies, and the vertical and streamwise time scale profiles were almost mirror images of each other. Eulerian length scale profiles showed some differences among canopies caused by differences in the velocity profiles. Ratios of vertical-to-horizontal time and length scales had a maximum in mid-canopy.Shear stress profiles were similar in the top parts of all canopies, and upward momentum fluxes were occasionally observed within the canopy trunk space. Countergradient heat fluxes were also observed sometimes. The countergradient fluxes and the skewed, kurtotic velocity distributions indicate the contribution of intermittent, large-scale eddies that are important for energy and mass transfer within canopies.