Turbulent kinetic energy budgets from a large-eddy simulation of airflow above and within a forest canopy

The output of a large-eddy simulation was used to study the terms ofthe turbulent kinetic energy (TKE) budget for the air layers above andwithin a forest. The computation created a three-dimensional,time-dependent simulation of the airflow, in which the lowest third ofthe domain was occupied by drag elements and heat sources to representthe forest. Shear production was a principal source of TKE in theupper canopy, diminishing gradually above tree-top height and moresharply with depth in the canopy. The transfer of energy to subgridscales (dissipation) was the main sink in the upper part of the domainbut diminished rapidly with depth in the canopy. Removal ofresolved-scale TKE due to canopy drag was extremely important,occurring primarily in the upper half of the forest where the foliagedensity was large. Turbulent transport showed a loss at the canopytop and a gain within the canopy. These general features have beenfound elsewhere but uncertainty remains concerning the effects ofpressure transport. In the present work, pressure was calculateddirectly, allowing us to compute the pressure diffusion term. Wellabove the canopy, pressure transport was smaller than, and opposite insign to, the turbulent transport term. Near the canopy top andbelow, pressure transport acted in concert with turbulent transport toexport TKE from the region immediately above and within the uppercrown, and to provide turbulent energy for the lower parts of theforest. In combination, the transport terms accounted for over half ofthe TKE loss near the canopy top, and in the lowest two-thirds of thecanopy the transport terms were the dominant source terms in thebudget. Moreover, the pressure transport was the largest source ofturbulent kinetic energy in the lowest levels of the canopy, beingparticularly strong under convective conditions. These resultsindicate that pressure transport is important in the plant canopyturbulent kinetic energy budget, especially in the lowest portion ofthe stand, where it acts as the major driving force for turbulentmotions.

[1]  J. Deardorff.  Stratocumulus-capped mixed layers derived from a three-dimensional model , 1980 .

[2]  Steven A. Orszag,et al.  Pseudospectral approximation to two-dimensional turbulence , 1973 .

[3]  R. Shaw,et al.  The Influence of Atmospheric Stability on the Budgets of the Reynolds Stress and Turbulent Kinetic Energy within and above a Deciduous Forest , 1990 .

[4]  T. Meyers,et al.  Testing of a higher-order closure model for modeling airflow within and above plant canopies , 1986 .

[5]  Peter A. Coppin,et al.  Experiments on scalar dispersion within a model plant canopy part II: An elevated plane source , 1986 .

[6]  N. Wilson,et al.  A Higher Order Closure Model for Canopy Flow , 1977 .

[7]  M. Raupach,et al.  Experiments on scalar dispersion within a model plant canopy part I: The turbulence structure , 1986 .

[8]  John C. Wyngaard,et al.  Evaluation of turbulent transport and dissipation closures in second-order modeling , 1989 .

[9]  Ulrich Schumann,et al.  Large-eddy simulation of turbulent flow above and within a forest , 1992 .

[10]  T. Meyers,et al.  The budgets of turbulent kinetic energy and Reynolds stress within and above a deciduous forest , 1991 .

[11]  C. Moeng.  A Large-Eddy-Simulation Model for the Study of Planetary Boundary-Layer Turbulence , 1984 .

[12]  U. Schumann,et al.  Large‐eddy simulation of a neutrally stratified boundary layer: A comparison of four computer codes , 1994 .

[13]  Manabu Kanda,et al.  Organized structures in developing turbulent flow within and above a plant canopy, using a Large Eddy Simulation , 1994 .

[14]  G. McBean,et al.  The Vertical Transports of Kinetic Energy by Turbulence and Pressure in the Boundary Layer , 1975 .

[15]  John C. Wyngaard,et al.  Spectral analysis of large-eddy simulations of the convective boundary layer , 1988 .

[16]  R. Shaw,et al.  Influence of foliar density and thermal stability on profiles of Reynolds stress and turbulence intensity in a deciduous forest , 1988 .

[17]  Retrieval of turbulent pressure fluctuations at the ground surface beneath a forest , 1990 .

[18]  B. Amiro,et al.  Pressure Fluctuations during Coherent Motions and Their Effects on the Budgets of Turbulent Kinetic Energy and Momentum Flux within a Forest Canopy , 1994 .

[19]  P. Mason.  Atmospheric boundary layer flows: Their structure and measurement , 1995 .

[20]  Ulrich Schumann,et al.  Large-Eddy Simulation of the Convective Boundary Layer: A Comparison of Four Computer Codes , 1993 .

[21]  On Pressure-Driven Wind in Deep Forests , 1989 .

[22]  E. Shaughnessy,et al.  Microscale pressure fluctuations in a mature deciduous forest , 1983 .

[23]  X. J. Zhang,et al.  Evidence of pressure-forced turbulent flow in a forest , 1992 .

[24]  Turbulent Shear Flows 8 , 1993 .

[25]  Tilden P. Meyers,et al.  Modelling the plant canopy micrometeorology with higher-order closure principles , 1987 .

[26]  J. Finnigan,et al.  A wind tunnel study of air flow in waving wheat: Two-point velocity statistics , 1994 .

[27]  P. Mason.  Large‐eddy simulation: A critical review of the technique , 1994 .

[28]  J. Finnigan,et al.  A wind tunnel study of air flow in waving wheat: Single-point velocity statistics , 1994 .

[29]  T. Maitani,et al.  Estimates of velocity-pressure and velocity-pressure gradient interactions in the surface layer over plant canopies , 1985 .

[30]  P. Mason,et al.  Large-Eddy Simulation of the Convective Atmospheric Boundary Layer , 1989 .

[31]  John D. Wilson,et al.  A second-order closure model for flow through vegetation , 1988 .

[32]  M. Raupach,et al.  Averaging procedures for flow within vegetation canopies , 1982 .