The HartX-synthesis: An experimental approach to water and carbon exchange of a Scots pine plantation

Summary In May 1992 during the interdisciplinary measurement campaign HartX (Hartheim eXperiment), several independent estimates of stand water vapor flux were compared at a 12-m high Scots pine (Pinus silvestris) plantation on a flat fluvial terrace of the Rhine close to Freiburg, Germany. Weather during the HartX period was characterized by ten consecutive clear days with exceptionally high input of available energy for this time of year and with a slowly shifting diurnal pattern in atmospheric variables like vapor pressure deficit. Methods utilized to quantify components of stand water flux included porometry measurements on understory graminoid leaves and on pine needles and three different techniques for determining individual tree xylem sap flow. Micrometeorological methods included eddy covariance and eddy covariance energy balance techniques with six independent systems on two towers separated by 40 m. Additionally, Bowen ratio energy balance estimates of water flux were conducted and measurements of the gradients in water vapor, CO₂, and trace gases within and above the stand were carried out with an additional, portable 30 m high telescoping mast.Biologically-based estimates of overstory transpiration were obtained by up-scaling tree sap flow rates to stand level via cumulative sapwood area. Tree transpiration contributed between 2.2 and 2.6 mm/day to ET for a tree leaf area index (LAI) of 2.8. The pine stand had an understory dominated by sedge and grass species with overall average LAI of 1.5. Mechanistic canopy gas exchange models that quantify both water vapor and CO₂ exchange were applied to both understory and tree needle ecosystem compartments. Thus, the transpiration by graminoid species was estimated at approximately 20% of total stand ET. The modelled estimates for understory contribution to stand water flux compared well with micrometeorologically-based determinations. Maximum carbon gain was estimated from the canopy models at approximately 425 mmol/(m²day) for the tree needles and at 100 mmol/(m²day) for the understory. Carbon gain was suggested by the modelling analysis to remain relatively constant during the HartX period, while water use efficiency in carbon fixation increased with decreasing vapor pressure deficit. Biologically- and micrometeorologically-based estimates of stand water flux showed good general agreement with variation of up to 20% that reflects both errors due to the inherent assumptions associated with different methods as well as natural spatial variability in fluxes. The various methods support a reliable estimate of average ET from this homogeneous canopy during HartX of about 2.6 mm/day (a maximum of about 3.1 mm/day) with an insignificant decreasing trend in correlation with decreasing vapor pressure deficit and possibly soil moisture.Findings during HartX were embedded in local scale heterogeneity with greater roughness over the forest and much higher ET over the surrounding agricultural fields which results in weak but clearly existant circulation patterns. A variety of measurements were continued after the HartX campaign. They allow us to extend our findings for six months with changing environmental conditions, including shortage of soil moisture. Hydrological estimates of soil water extractions and micrometeorological estimates of ET by the one-propeller eddy covariance (OPEC) system were in very good agreement, supporting the use of this robust eddy covariance energy balance technique for long-term monitoring.With 5 Figures     

Comparisons of xylem sap flow and water vapour flux at the stand level and derivation of canopy conductance for Scots pine

Summary Simultaneous measurements of xylem sap flow and water vapour flux over a Scots pine (Pinus sylvestris) forest (Hartheim, Germany), were carried out during the Hartheim Experiment (HartX), an intensive observation campaign of the international programme REKLIP. Sap flow was measured every 30 min using both radial constant heating (Granier, 1985) and two types of Cermak sap flowmeters installed on 24 trees selected to cover a wide range of the diameter classes of the stand (min 8 cm; max 17.5 cm). Available energy was high during the observation period (5.5 to 6.9 mm.day^–1), and daily cumulated sap flow on a ground area basis varied between 2.0 and 2.7 mm day^–1 depending on climate conditions. Maximum hourly values of sap flow reached 0.33 mm h^–1, i.e., 230 W m^–2.Comparisons of sap flow with water vapour flux as measured with two OPEC (One Propeller Eddy Correlation, University of Arizona) systems showed a time lag between the two methods, sap flow lagging about 90 min behind vapour flux. After taking into account this time lag in the sap flow data set, a good agreement was found between both methods: sap flow = 0.745^* vapour flux,r ² = 0.86. The difference between the two estimates was due to understory transpiration.Canopy conductance (g _c ) was calculated from sap flow measurements using the reverse form of Penman-Monteith equation and climatic data measured 4 m above the canopy. Variations ofg _c were well correlated (r ² = 0.85) with global radiation (R) and vapour pressure deficit (vpd). The quantitative expression forg _c =f (R, vpd) was very similar to that previously found with maritime pine (Pinus pinaster) in the forest of Les Landes, South Western France.With 6 Figures     

Soil moisture variation and plant water stress at the Hartheim Scots pine plantation

Summary During two measurement campaigns in 1992 (the Hartheim Experiment HartX- and an additional experiment in autumn), measurements of soil moisture were carried out in aPinus sylvestris stand at Hartheim on the Oberrhein. Several methods were used to determine soil water status. They were compared in terms of suitability for estimating stand evapotranspiration (ET) via soil water depletion. Measurements of tree water potential suggested that conductance of the trees was affected by soil water depletion during the period of the HartX campaign in spring 1992. We interpret the observations to indicate a lesser influence of soil water availability on tree transpiration during the autumn experiment.Eddy correlation and xylem sapflow measurements provided reference ET values with which to compare the stand ET calculations based on soil moisture measurements. Profile measurements of soil moisture showed that particularly in springtime when the lower soil layers are saturated with water, the water transport from depths below the major rooting zone is a very important factor affecting evaluation of stand ET. Decreases in soil water storage may be determined best with permanently installed soil moisture sensors such as used in tensiometric or TDR measurements that provide high resolution of changes over time.With 8 Figures     

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

Swaying of trees in response to wind and thinning in a stand of Scots pine

Observations have been made of the windspeed, wind direction, and tree movement at the edge and 20 m within a stand of Scots pine (Pinus sylvestris L.) close to 11 m in height. The spectra of windspeed near canopy top, together with the output of accelerometers and video observations of tree movement at mid-crown, were compared in the same stand prior and two years after first thinning. Furthermore, the transfer of wind energy into tree movement was investigated by calculating the mechanical transfer function (H _m ² ) between the wind spectrum (S _uu) and the tree's response (S _yy), i.e. H _m ² = S_yy/S_uu. Trees were found to behave like damped harmonic oscillators. They reacted to sudden increases in windspeed, reached their greatest displacement during the first cycle, and then returned to their rest position under the influence of damping. The spectral peak frequencies in S _yy and in H _m ²coincided with the estimated natural sway frequency of trees. Response in the second mode was, however, also evident, especially within the unthinned stand. The periodogram plots showed a consistent trend of a marked decrease in the response of the tree to increase in frequency. Almost no difference in the wind energy transfer, i.e. peak frequencies and peak width, and damping of the system was found between Scot pine at 2700 and 1500 stems per hectare. However, along the stand edge tree movement was greater than within the stand indicating greater wind energy transfer and damping of the system along the stand edge than within the stand.     

Model-based estimates of water loss from “patches” of the understory mosaic of the Hartheim Scots pine plantation

Summary During the Hartheim Experiment (HartX) 1992 conducted in the Upper Rhine Valley, Germany, we estimated water vapor flux from the understory and the forest floor by several methods. At the vegetation patch level, direct estimates were made with small weighing lysimeters, and water loss was scaled-up to the stand level based on vegetation patchtype distribution. At the leaf level, transpiration flux was determined with a CO₂/H₂O porometer for the dominant understory plant species,Brachypodium pinnatum, Carex alba, andCarex flacca. Measured leaf transpiration was scaled-up to patch level with a canopy light interception and leaf gas exchange model, and then to stand level as in the case of lysimeter data, but with further consideration of patchtype leaf area index (LAI). On two days, total understory latent heat flux was estimated by eddy correlation methods below the tree canopy.The understory vegetation was subdivided into five major patch-types which covered 62% of the ground area and resulted in a cumulative LAI of approx. 1.54 when averaged over total stand ground area and compared to the average tree canopy LAI of 2.8. The remaining 38% of ground area was unvegetated bare soil and/or covered by moss (mainly byScleropodium purum) or litter. The evapotranspiration from the understory and unvegetated areas equaled approx. 20% of total forest stand transpiration during the HartX period. The understory vegetation transpired about 0.4 mm d^–1 (13%) estimated over the period of May 13 to 21, whereas evaporation from moss and soil patches amounted 0.23 mm d^–1 (7.0%). On dry, sunny days, total water vapor flux below the tree canopy exceeded 0.66 mm d^–1. Using the transpiration rates derived from the GAS-FLUX model together with estimates of evaporation from moss and soil areas and a modified application of the Penman-Monteith equation, the average daily maximum conductance of the understory and the forest floor was 1.7 mm s^–1 as compared to 5.5 mm s^–1 for the tree canopy.With 6 Figures     

Aerodynamic Scaling for Estimating the Mean Height of Dense Canopies

We used an aerodynamic method to objectively determine a representative canopy height, using standard meteorological measurements. The canopy height may change if the tree height is used to represent the actual canopy, but little work to date has focused on creating a standard for determining the representative canopy height. Here we propose the ‘aerodynamic canopy height’ h _a as the most effective means of resolving the representative canopy height for all forests. We determined h _a by simple linear regression between zero-plane displacement d and roughness length z ₀, without the need for stand inventory data. The applicability of h _a was confirmed in five different forests, including a forest with a complex canopy structure. Comparison with stand inventory data showed that h _a was almost equivalent to the representative height of trees composing the crown surface if the forest had a simple structure, or to the representative height of taller trees composing the upper canopy in forests with a complex canopy structure. The linear relationship between d and z ₀ was explained by assuming that the logarithmic wind profile above the canopy and the exponential wind profile within the canopy were continuous and smooth at canopy height. This was supported by observations, which showed that h _a was essentially the same as the height defined by the inflection point of the vertical profile of wind speed. The applicability of h _a was also verified using data from several previous studies.     

Energy balance and transpiration from tussock grassland in New Zealand

The energy balance was measured for the dry canopy of narrow-leaved snow tussock (Chionochloa rigida), and measurements of transpiration were obtained from a large weighing lysimeter.Typical maximum summer transpiration rates of 0.21–0.43 mmhr^-1 (140–290 W m^-2) were recorded. The latent heat flux accounted for less than 40% of net radiation. The estimated value of the bulk stomatal resistance (r _ST) for 29 days was 158 s m^-1, and the decoupling parameter () was 0.17. Transpiration rates were not driven directly by net radiation, but were closely linked to the size of the regional saturation deficit imposed at the level of the canopy by efficient overhead mixing, and were constrained by a large bulk stomatal resistance. A linear relationship between r _ST and the saturation deficit is proposed as a realistic method for estimating transpiration for water yield studies of tussock catchments.     

Water vapor transfer from a vegetated surface: A numberical study of bulk transfer coefficients and canopy resistances

The semi-analytical model outlined in previous studies (Massman, 1987a, b) to describe momentum and heat exchange between the atmosphere and vegetated surfaces is extended to include water vapor exchange. The methods employed are based on one-dimensional turbulent diffusivities and use numerical solutions to the steady-state diffusion equation. The model formulates stomatal response as a function of vapor pressure deficit and the within-canopy profile of mean photosynthetically-active radiation (PAR). It is then used to assess the influence that foliage structure, density, and sheltering can have upon the bulk transfer coefficient, kB _v ^-1, and the canopy resistance. A general analytical formulation of the canopy resistance based on the mean within-canopy profile of PAR is proposed and found to agree with the model's solutions for canopy resistance to within a few percent.     

Influence of leaf water potential on diurnal changes in CO2 and water vapour fluxes

Mass and energy fluxes between the atmosphere and vegetation are driven by meteorological variables, and controlled by plant water status, which may change more markedly diurnally than soil water. We tested the hypothesis that integration of dynamic changes in leaf water potential may improve the simulation of CO₂ and water fluxes over a wheat canopy. Simulation of leaf water potential was integrated into a comprehensive model (the ChinaAgrosys) of heat, water and CO₂ fluxes and crop growth. Photosynthesis from individual leaves was integrated to the canopy by taking into consideration the attenuation of radiation when penetrating the canopy. Transpiration was calculated with the Shuttleworth-Wallace model in which canopy resistance was taken as a link between energy balance and physiological regulation. A revised version of the Ball-Woodrow-Berry stomatal model was applied to produce a new canopy resistance model, which was validated against measured CO₂ and water vapour fluxes over winter wheat fields in Yucheng (36°57′ N, 116°36′ E, 28 m above sea level) in the North China Plain during 1997, 2001 and 2004. Leaf water potential played an important role in causing stomatal conductance to fall at midday, which caused diurnal changes in photosynthesis and transpiration. Changes in soil water potential were less important. Inclusion of the dynamics of leaf water potential can improve the precision of the simulation of CO₂ and water vapour fluxes, especially in the afternoon under water stress conditions.     

A Measurement Based Sensitivity Analysis of the Penman-Monteith Actual Evapotranspiration Model for Crops of Different Height and in Contrasting Water Status

Summary This paper presents a study of the sensibility of the Penman-Monteith evapotranspiration model to climatic (available energy and vapour pressure deficit) and parametric (aerodynamic and canopy resistances, r _a and r _c respectively) factors in a semi-arid climate, for crops in contrasting water status (well irrigated and under water stress) and of different heights. Three experiments were carried out in southern Italy on reference grass (≈ 0.1 m), grain sorghum (≈ 1 m) and sweet sorghum (≈ 3 m). For this analysis the sensitivity coefficients, taken as hourly means, were evaluated during the growth season when the crops completely covered the soil. The relative errors on evapotranspiration were also evaluated for r _a and r _c . The results showed that, for reference grass, available energy and aerodynamic resistance play a major role. For crops under water stress the most important term to evaluate is canopy resistance. For a tall crop, as sweet sorghum, the role of the vapour pressure deficit is fundamental, both when the crop is in good water status and under water stress. Received July 14, 1997 Revised February 5, 1998     

Calibration of a land-surface model using data from primary forest sites in Amazonia

Summary A land-surface model (MOSES) was tested against observed fluxes of heat, water vapour and carbon dioxide for two primary forest sites near Manaus, Brazil. Flux data from one site (called C14) were used to calibrate the model, and data from the other site (called K34) were used to validate the calibrated model. Long-term fluxes of water vapour at C14 and K34 simulated by the uncalibrated model were good, whereas modelled net ecosystem exchange (NEE) was poor. The uncalibrated model persistently underpredicted canopy conductance (g _c ) from mid-morning to mid-afternoon due to saturation of the response to solar radiation at low light levels. This in turn caused a poor simulation of the diurnal cycles of water vapour and carbon fluxes. Calibration of the stomatal conductance/photosynthesis sub-model of MOSES improved the simulated diurnal cycle of g _c and increased the diurnal maximum NEE, but at the expense of degrading long-term water vapour fluxes. Seasonality in observed canopy conductance due to soil moisture change was not captured by the model. Introducing realistic depth-dependent soil parameters decreased the amount of moisture available for transpiration at each depth and led to the model experiencing soil moisture limitation on canopy conductance during the dry season. However, this limitation had only a limited effect on the seasonality in modelled NEE.     

Climate change impacts on dryland cropping systems in the Central Great Plains,USA

Agricultural systems models are essential tools to assess potential climate change (CC) impacts on crop production and help guide policy decisions. In this study, impacts of projected CC on dryland crop rotations of wheat-fallow (WF), wheat-corn-fallow (WCF), and wheat-corn-millet (WCM) in the U.S. Central Great Plains (Akron, Colorado) were simulated using the CERES V4.0 crop modules in RZWQM2. The CC scenarios for CO₂, temperature and precipitation were based on a synthesis of Intergovernmental Panel on Climate Change (IPCC 2007) projections for Colorado. The CC for years 2025, 2050, 2075, and 2100 (CC projection years) were super-imposed on measured baseline climate data for 15–17 years collected during the long-term WF and WCF (1992–2008), and WCM (1994–2008) experiments at the location to provide inter-annual variability. For all the CC projection years, a decline in simulated wheat yield and an increase in actual transpiration were observed, but compared to the baseline these changes were not significant (p > 0.05) in all cases but one. However, corn and proso millet yields in all rotations and projection years declined significantly (p < 0.05), which resulted in decreased transpiration. Overall, the projected negative effects of rising temperatures on crop production dominated over any positive impacts of atmospheric CO₂ increases in these dryland cropping systems. Simulated adaptation via changes in planting dates did not mitigate the yield losses of the crops significantly. However, the no-tillage maintained higher wheat yields than the conventional tillage in the WF rotation to year 2075. Possible effects of historical CO₂ increases during the past century (from 300 to 380 ppm) on crop yields were also simulated using 96 years of measured climate data (1912–2008) at the location. On average the CO₂ increase enhanced wheat yields by about 30%, and millet yields by about 17%, with no significant changes in corn yields.     

Aerodynamic and canopy resistance of short-rotation forest in relation to leaf area index and climate

The aerodynamic and canopy resistances of a willow short-rotation stand were estimated during the course of a growing season on the basis of micrometeorological measurements. The normalized roughness length (z ₀/h) decreased from about 0.10 at a leaf area index of one to 0.05 at a leaf area index of seven. This implies that the aerodynamic resistance at peak leaf area index is more than twice the value at zero leaf area index, all other variables unchanged. The canopy resistance depended strongly on air water concentration deficit and on leaf area index. The Lohammar equation showed good agreement between estimated and measured canopy resistances over the whole course of leaf development. The stand was well-coupled to the atmosphere only for very small values of leaf area indices, less than one, and it was practically de-coupled for leaf area indices above two. From the point of view of factors controlling evaporation, this type of stand acts as a traditional forest at the beginning and end of the season and as an agricultural crop in the middle of the season.     

The impact of new land surface physics on the GCM simulation of climate and climate sensitivity

 Recent improvements to the Hadley Centre climate model include the introduction of a new land surface scheme called “MOSES” (Met Office Surface Exchange Scheme). MOSES is built on the previous scheme, but incorporates in addition an interactive plant photosynthesis and conductance module, and a new soil thermodynamics scheme which simulates the freezing and melting of soil water, and takes account of the dependence of soil thermal characteristics on the frozen and unfrozen components. The impact of these new features is demonstrated by comparing 1×CO₂ and 2×CO₂ climate simulations carried out using the old (UKMO) and new (MOSES) land surface schemes. MOSES is found to improve the simulation of current climate. Soil water freezing tends to warm the high-latitude land in the northern Hemisphere during autumn and winter, whilst the increased soil water availability in MOSES alleviates a spurious summer drying in the mid-latitudes. The interactive canopy conductance responds directly to CO₂, supressing transpiration as the concentration increases and producing a significant enhancement of the warming due to the radiative effects of CO₂ alone. Received: 16 March 1998 / Accepted: 4 August 1998     

A coupled model on land-atmosphere interactions — Simulating the characteristics of the PBL over a heterogeneous surface

An attempt is made to construct a model, coupling land surface and atmospheric processes in the planetary boundary layer (PBL). A grassland strip in a semi-desert (hereinafter called desert) is presupposed, so as to simulate the case of heterogeneous vegetation cover.Modeling results indicate that every term in the equation of the surface energy balance changes as the air flows over the grassland. The striking contrast of water and energy conditions between the grassland and the desert means that the air over the grassland is cooler and wetter than that over the desert. Consequently, in the heating and dynamic forcing of the air by the underlying surface, heterogeneities arise and are then transferred upward by the turbulent motions. Horizontal differences thus develop in the PBL, resulting in a local circulation. Meanwhile, the horizontal differences affect the free atmosphere through vertical motion at the top of the PBL.List of symbols d ₁,d ₂,d ₃ depths of surface, middle and lower layers of soil - T _c ,T ₁,T ₂,T ₃ temperatures of canopy, surface, middle and lower layers of soil - R _nc net radiation of canopy layer - _c shielding factor of vegetation - Ew, Etc evaporation from wet fraction of foliage and transpiration from dry fraction of foliage - Et ₁,Et ₂ transpiration of foliage water absorbed by the root in the upper and lower soil, respectively - H _c sensible heat of canopy - P _c ,D _c precipitation rate and drainage of canopy - C _s ,C _c ,C _w heat capacity of soil, canopy and water - _w , _s density of water and air near the surface - D hydraulic permeability of soil - _s saturated value of the ratio of volumetric soil moisture - S _g , _g solar radiation and surface reflection - H _g ,R _{L g} turbulent heat flux and long wave radiation of surface - P _g ,E _g precipitation rate and evaporation of soil surface - K _s soil thermal diffusivity - K ^(m),K ^(H),K ^(q) eddy coefficients of momentum, heat and moisture - u, v, w components of wind speed in three directions - air potential temperature - e turbulent kinetic energy - p atmospheric pressure - C _p specific heat of air under constant pressure - R _d gas constant - u _* friction velocity - _* feature temperature - h height of the PBL - f Coriolis parameter - L ₀ Monin-Obukhov length - latent heat of vaporization - q specific humidity - M _c ,M _cm interception water storage of canopy and its maximum - ₀ Exner number of largescale background field - perturbation Exner number - u _g ,v _g components of the geostrophic wind speed Sponsored by the National Natural Science Foundation of China.     

Heat transfer to and from vegetated surfaces: An analytical method for the bulk exchange coefficients

The semi-analytical model outlined in a previous study (Massman, 1987) to describe momentum exchange between the atmosphere and vegetated surfaces is extended to include the exchange of heat. The methods employed are based on one-dimensional turbulent diffusivities and use analytical solutions to the steady-state diffusion equation. The model is used to assess the influence that (1) the canopy foliage structure and density, (2) the wind profile structure within the canopy, and (3) the shelter factor can have upon the inverse surface Stanton number (kB ^–1) as well as to explore the consequences of introducing a scalar displacement height which can be different from the momentum displacement height. In general, the triangular foliage area density function gives results which agree more closely to observations than that for constant foliage area density. The intended application of this work and its predecessor (Massman, 1987) is for parameterizing the bulk aerodynamic resistances for heat and momentum exchange for use within large-scale models of plant-atmosphere exchanges.     

Heat and mass transfer to and from surfaces with dense vegetation or similar permeable roughness

A differential equation is obtained to describe the concentration of passive admixtures (water vapor, sensible heat, pollutants, CO₂, etc.) of turbulent flow inside a dense and uniform vegetational canopy. The profiles of eddy diffusivity, wind speed and shear stress are assumed to be exponential decay functions of depth below the top of the canopy. This equation is solved for the case of a vegetation with constant concentration of the admixture at the foliage surfaces. The solution is used to formulate bulk mass or heat transfer coefficients, which can be applied to practical problems involving surfaces covered with a vegetation or with similar porous or fibrous roughness elements. The results are shown to be consistent with experimental data presented by Chamberlain (1966), Garratt and Hicks (1973) and Garratt (1978). Calculations with the model illustrate that, as compared to its behavior over surfaces with bluff roughness elements, ln(z ₀/z _0c ) (wherez ₀ is the momentum roughness andz _0c , the scalar roughness) for permeable roughness elements is relatively insensitive tou _* and practically independent ofz ₀.     

Atmospheric Stability Effects on Wind Fields and Scalar Mixing Within and Just Above a Subalpine Forest in Sloping Terrain

Air temperature T _a , specific humidity q, CO₂ mole fraction χ _c , and three-dimensional winds were measured in mountainous terrain from five tall towers within a 1 km region encompassing a wide range of canopy densities. The measurements were sorted by a bulk Richardson number Ri _b . For stable conditions, we found vertical scalar differences developed over a “transition” region between 0.05 < Ri _b < 0.5. For strongly stable conditions (Ri _b > 1), the vertical scalar differences reached a maximum and remained fairly constant with increasing stability. The relationships q and χ _c have with Ri _b are explained by considering their sources and sinks. For winds, the strong momentum absorption in the upper canopy allows the canopy sublayer to be influenced by pressure gradient forces and terrain effects that lead to complex subcanopy flow patterns. At the dense-canopy sites, soil respiration coupled with wind-sheltering resulted in CO₂ near the ground being 5–7 μmol mol⁻¹ larger than aloft, even with strong above-canopy winds (near-neutral conditions). We found Ri _b -binning to be a useful tool for evaluating vertical scalar mixing; however, additional information (e.g., pressure gradients, detailed vegetation/topography, etc.) is needed to fully explain the subcanopy wind patterns. Implications of our results for CO₂ advection over heterogenous, complex terrain are discussed.     

Numerical Investigations of Mean Winds Within Canopies of Regularly Arrayed Cubical Buildings Under Neutral Stability Conditions

Recently, several attempts have been made to model the wind velocity in an urban canopy in order to accurately predict the mixing and transport of momentum, heat, and pollutants within and above the canopy on an urban scale. For this purpose, unverified assumptions made by Macdonald (Boundary-Layer Meteorol 97:25–45, 2000) to develop a model for the profile of the mean wind velocity within an urban canopy have been used. In the present study, in order to provide foundations for improving the urban canopy models, the properties of the spatially-averaged mean quantities used to make these assumptions have been investigated by performing large-eddy simulations (LES) of the airflow around square and staggered arrays of cubical blocks with the following plan area densities: λ _p = 0.05, 0.11, 0.16, 0.20, 0.25, and 0.33. The LES results confirm that the discrepancy between the spatial average of wind velocity and Macdonald’s five-point average of wind velocity can be large in both types of arrays for large λ _p . It is also confirmed that Prandtl’s mixing length varies significantly with height within the canopy, contrary to Macdonald’s assumption for both types of arrays and for both small and large λ _p . On the other hand, in accordance with Macdonald’s assumption, the sectional drag coefficient is found to be almost constant with height except in the case of staggered arrays with high λ _p .