John D. Albertson

THE EFFECT OF VEGETATION DENSITY ON CANOPY SUB-LAYER TURBULENCE

The canonical form of atmospheric flows near theland surface, in the absence of a canopy, resembles a rough-wallboundary layer. However, in the presence of an extensive and densecanopy, the flow within and just above the foliage behaves as aperturbed mixing layer. To date, no analogous formulation existsfor intermediate canopy densities. Using detailed laser Dopplervelocity measurements conducted in an open channel over a widerange of canopy densities, a phenomenological model that describesthe structure of turbulence within the canopy sublayer (CSL) isdeveloped. The model decomposes the space within the CSL intothree distinct zones: the deep zone in which the flow field isshown to be dominated by vortices connected with vonKármán vortex streets, butperiodically interrupted by strong sweep events whose features areinfluenced by canopy density. The second zone, which is near thecanopy top, is a superposition of attached eddies andKelvin–Helmholtz waves produced by inflectional instability in themean longitudinal velocity profile. Here, the relative importanceof the mixing layer and attached eddies are shown to vary withcanopy density through a coefficient α. We show that therelative enhancement of turbulent diffusivity over its surface-layer value near the canopy top depends on the magnitude ofα. In the uppermost zone, the flow follows the classicalsurface-layer similarity theory. Finally, we demonstrate that thecombination of this newly proposed length scale and first-orderclosure models can accurately reproduce measured mean velocity andReynolds stresses for a wide range of roughness densities. Withrecent advancement in remote sensing of canopy morphology, thismodel offers a promising physically based approach to connect theland surface and the atmosphere without resorting to empiricalmomentum roughness lengths.

Surface length scales and shear stress: Implications for land-atmosphere interaction over complex terrain

A large eddy simulation (LES) code of the atmospheric boundary layer (ABL) has been developed and applied to study the effect of spatially variable surface properties on the areally averaged surface shear stress at the land-atmosphere interface. The LES code simulates the space and time evolution of the large-scale turbulent eddies and their transport effects in the ABL. We report here on simulations of flow over spatially variable roughness fields. The dynamics are simulated, and the resulting space-time fields are averaged to explore the effects of the surface variability length scales on the average surface shear stress, as used in large-scale models to estimate scalar fluxes, such as evaporation. We observe asymmetrical response of the smooth-to-rough and rough-to- smooth transitions, such that the effects of the transitions accumulate rather than cancel. It is shown that the presence of abrupt changes in surface roughness and the atmospheres response to these patches create a marked dependence of the statistical structure of surface shear stress on the length scale of the surface patches. An increase in regionally averaged surface stress for decreasing horizontal patch length scale is found.

Natural integration of scalar fluxes from complex terrain

Large eddy simulations of turbulent flow and transport in the atmospheric boundary layer were conducted over heterogeneous sources of heat and water vapor to identify the blending properties of the turbulent mixing in an unstably stratified boundary layer. The numerical simulations show that the concept of blending in the ABL is in fact a useful one, even under convective conditions, for a range of surface conditions. Since the transport eddies that are responsible for the blending have sizes that are constrained by the boundary layer depth, and since the vertical motion is so important under the unstable density stratification studied here, we see that a hen the length scales of the source variability on the land surface become significantly greater than the ABL depth the blending is lost. In this case the source fields remain relatively uncoupled by the important eddy motion. However, for smaller surface length scales, the dynamic eddy motion couples the surface patches. Hence, there is good reason that the land surface exchange phenomenon would not be scale invariant over the entire range of scales, Because of the active role of temperature the effects of inhomogeneous surface sources of sensible heat persist higher into the ABL than do the effects of surface sources from more passive scalars, such as water vapor. Moreover, the mean fields of potential temperature and specific humidity blend at much lower heights than do the vertical turbulent flux fields of these two scalars. We propose a useful measure of blending efficiency for simulation studies and show how this bridges from the dynamics responsible for the blending to the horizontal homogeneity of scalar flux fields at measurement heights in the ABL.

Determining land surface fractional cover from NDVI and rainfall time series for a savanna ecosystem

Savanna ecosystems are water limited and responsive to rainfall on short time scales, characteristics that can be exploited to estimate fractional cover of trees, grass, and bare soil over large-scale areas from synthesis of remote sensing and rainfall measurements. A method is presented to estimate fractional cover based upon the differing ways in which grasses and trees respond to rainfall, and implementation of this method is demonstrated along the Kalahari Transect (KT), an aridity gradient in southern Africa. Seasonally averaged normalized difference vegetation index (NDVI) and the sensitivity of the NDVI to interannual variations in wet season rainfall are used as state-space variables in a linear unmixing model. End members for this analysis were determined on the basis of best fit to the observed data. The realized end members were consistent with the qualitative characteristics of trees (high NDVI, low sensitivity of NDVI to interannual variations in rainfall), bare soil (low NDVI, low sensitivity), and the transient grass/ bare soil area (moderate NDVI, high sensitivity). Observed sensitivity of NDVI to rainfall was measured as the relationship between wet season NDVI and normalized rainfall over a 16-year period (1983–1998). The unmixing model yields a north-to-south decrease in tree fractional cover that corresponds to the decrease in mean wet season precipitation from 1600 to 300 mm along the KT (R 2 =.87). The fractional tree cover results compare favorably with available ground-based observations. The potential extent of grass cover is limited by the dominance of trees on the northern end of the transect, peaks at the location with approximately 450 mm of mean wet season rainfall, and is limited by rainfall in the arid southern portion of the transect. With mean NDVI for grass inferred from the data, yearly estimates of tree, grass, and bare soil fractional cover can be derived. These annual estimates, which are difficult to obtain from traditional unmixing procedures, are important parameters in fuel load and land–atmosphere exchange models. No calibration or training sets were required for this analysis, and this method has the additional capability to predict fractional-cover components for future rainfall scenarios. D 2002 Published by Elsevier Science Inc.

Multiscale analysis of vegetation surface fluxes: from seconds to years

Abstract The variability in land surface heat (H), water vapor (LE), and CO2 (or net ecosystem exchange, NEE) fluxes was investigated at scales ranging from fractions of seconds to years using eddy-covariance flux measurements above a pine forest. Because these fluxes significantly vary at all these time scales and because large gaps in the record are unavoidable in such experiments, standard Fourier expansion methods for computing the spectral and cospectral statistical properties were not possible. Instead, orthonormal wavelet transformations ( OWT ) are proposed and used. The OWT are ideal at resolving process variability with respect to both scale and time and are able to isolate and remove the effects of missing data (or gaps) from spectral and cospectral calculations. Using the OWT spectra, we demonstrated unique aspects in three appropriate ranges of time scales: turbulent time scales (fractions of seconds to minutes), meteorological time scales (hour to weeks), and seasonal to interannual time scales corresponding to climate and vegetation dynamics. We have shown that: (1) existing turbulence theories describe the short time scales well, (2) coupled physiological and transport models (e.g. CANVEG) reproduce the wavelet spectral characteristics of all three land surface fluxes for meteorological time scales, and (3) seasonal dynamics in vegetation physiology and structure inject strong correlations between land surface fluxes and forcing variables at monthly to seasonal time scales. The broad implications of this study center on the possibility of developing low-dimensional models of land surface water, energy, and carbon exchange. If the bulk of the flux variability is dominated by a narrow band or bands of modes, and these modes “resonate” with key state and forcing variables, then low-dimensional models may relate these forcing and state variables to NEE and LE.

Soil moisture controls on canopy‐scale water and carbon fluxes in an African savanna

[1] Tower-based measurements of mass and energy exchanges at the end of the growing season in central Botswana were used to evaluate functional relationships commonly applied to predict water and carbon fluxes between savanna landscapes and the atmosphere. Following a large rainfall event, daily evapotranspiration (ETdaily) exhibited an exponential decay consistent with a derived analytical expression based on critical and wilting-point soil moisture limits for savanna vegetation native to the study region. A piecewise linear soil moisture limitation function provided good estimates of ETdaily as a function of potential evapotranspiration and soil moisture (R 2 = 0.92). Comparison of a soil moisture mass balance with measured ETdaily indicated deeper root water uptake at a site with more woody vegetation compared with a grass-dominated site. Linear correlation (R 2 = 0.90) of daytime CO2 flux and evapotranspiration supported a constant water use efficiency to estimate carbon fluxes from water fluxes. Daytime and nighttime CO2 fluxes responded similarly to soil drying, enabling estimation of total daily CO2 flux from ETdaily. These experimental results support a simple model of savanna land-atmosphere exchange over interstorm periods. INDEX TERMS: 1818 Hydrology: Evapotranspiration; 3322 Meteorology and Atmospheric Dynamics: Land/atmosphere interactions; 1866 Hydrology: Soil moisture; 1878 Hydrology: Water/energy interactions; KEYWORDS: African savanna, evapotranspiration, land-atmosphere exchange, soil moisture, water and carbon flux, water limitation

Recession flow analysis for aquifer parameter determination

The recession flow analysis of Brutsaert and Nieber (1977) extended by Troch et al. (1993) to estimate aquifer parameters (saturated hydraulic conductivity and mean aquifer depth) is examined by means of a numerical model. It is found to be reliable for the estimation of the catchment-scale saturated hydraulic conductivity and mean aquifer depth. Increasing the complexity of the synthetic watershed had no impact on the accuracy of the estimated parameters.

SPATIAL VARIABILITY OF TURBULENT FLUXES IN THE ROUGHNESS SUBLAYER OF AN EVEN-AGED PINE FOREST

The spatial variability of turbulent flow statistics in the roughness sublayer (RSL) of a uniform even-aged 14 m (= h) tall loblolly pine forest was investigated experimentally. Using seven existing walkup towers at this stand, high frequency velocity, temperature, water vapour and carbon dioxide concentrations were measured at 15.5 m above the ground surface from October 6 to 10 in 1997. These seven towers were separated by at least 100m from each other. The objective of this study was to examine whether single tower turbulence statistics measurements represent the flow properties of RSL turbulence above a uniform even-aged managed loblolly pine forest as a best-case scenario for natural forested ecosystems. From the intensive space-time series measurements, it was demonstrated that standard deviations of longitudinal and vertical velocities (σu, σw) and temperature (σT) are more planar homogeneous than their vertical flux of momentum (u*2) and sensible heat (H) counterparts. Also, the measured H is more horizontally homogeneous when compared to fluxes of other scalar entities such as CO2 and water vapour. While the spatial variability in fluxes was significant (>15 %), this unique data set confirmed that single tower measurements represent the ‘canonical’ structure of single-point RSL turbulence statistics, especially flux-variance relationships. Implications to extending the ‘moving-equilibrium’ hypothesis for RSL flows are discussed. The spatial variability in all RSL flow variables was not constant in time and varied strongly with spatially averaged friction velocity u*, especially when u* was small. It is shown that flow properties derived from two-point temporal statistics such as correlation functions are more sensitive to local variability in leaf area density when compared to single point flow statistics. Specifically, that the local relationship between the reciprocal of the vertical velocity integral time scale (Iw) and the arrival frequency of organized structures (ū/h) predicted from a mixing-layer theory exhibited dependence on the local leaf area index. The broader implications of these findings to the measurement and modelling of RSL flows are also discussed.

A Note On The Contribution Of Dispersive Fluxes To Momentum Transfer Within Canopies

Dispersive flux terms are formed when the time-averaged meanmomentum equation is spatially averaged within the canopy volume.These fluxes represent a contribution to momentum transfer arisingfrom spatial correlations of the time-averaged velocity componentswithin a horizontal plane embedded in the canopy sublayer (CSL).Their relative importance to CSL momentum transfer is commonlyneglected in model calculations and in nearly all fieldmeasurement interpretations. Recent wind-tunnel studies suggestthat these fluxes may be important in the lower layers of thecanopy; however, no one study considered their importance acrossall regions of the canopy and for a wide range of canopy roughnessdensities. Using detailed laser Doppler anemometry measurementsconducted in a model canopy composed of cylinders within a largeflume, we demonstrate that the dispersive fluxes are onlysignificant (i.e., >10%) for sparse canopies. These fluxes arein the same direction as the turbulent flux in the lower layers ofthe canopy but in the opposite direction near the canopy top. Fordense canopies, we show that the dispersive fluxes are <5% atall heights. These results appear to be insensitive to theReynolds number (at high Reynolds numbers).

Regional scale evaporation and the atmospheric boundary layer

Note: 33(1): 99-124 Reference EFLUM-ARTICLE-1995-006doi:10.1029/94RG03112 Record created on 2005-09-08, modified on 2017-02-23