The higher-order scalar concentration fluctuation properties are examined in the context of Monin–Obukhov similarity theory for a variety of greenhouse gases that have distinct and separate source/sink locations along an otherwise ideal micrometeorological field site. Air temperature and concentrations of water vapour, carbon dioxide and methane were measured at high frequency (10 Hz) above a flat and extensive peat-land soil in the San Joaquin–Sacramento Delta (California, USA) area, subjected to year-round grazing by beef cattle. Because of the heterogeneous distribution of the sources and sinks of CO2 and especially CH4 emitted by cattle, the scaling behaviour of the higher-order statistical properties diverged from predictions based on a balance between their production and dissipation rate terms, which can obtained for temperature and H2O during stationary conditions. We identify and label these departures as ‘exogenous’ because they depend on heterogeneities and non-stationarities induced by boundary conditions on the flow. Spectral analysis revealed that the exogenous effects show their signatures in regions with frequencies lower than those associated with scalar vertical transport by turbulence, though the two regions may partially overlap in some cases. Cospectra of vertical fluxes appear less influenced by these exogenous effects because of the modulating role of the vertical velocity at low frequencies. Finally, under certain conditions, the presence of such exogenous factors in higher-order scalar fluctuation statistics may be ‘fingerprinted’ by a large storage term in the mean scalar budget.

Abstract: The interpretation of tower‐based eddy‐covariance (EC) turbulent flux measurements above forests hinges on three key assumptions: (1) steadiness in the flow statistics, (2) planar homogeneity of scalar sources or sinks, and (3) planar homogeneity in the flow statistics. Large eddy simulations (LESs) were used to control the first two so as to explore the break‐down of the third for idealized and real gentle topography such as those encountered in Amazonia. The LES runs were conducted using uniformly distributed sources inside homogeneous forests covering complex terrain to link the spatial patterns of scalar turbulent fluxes to topographic features. Results showed strong modulation of the fluxes by flow features induced by topography, including large area with negative fluxes compensating “chimney” regions with fluxes almost an order of magnitude larger than the landscape flux. Significant spatial heterogeneity persisted up to at least two canopy heights, where most eddy‐covariance measurements are performed above tall forests. A heterogeneity index was introduced to characterize and contrast different scenarios, and a topography categorization was shown to have predictive capabilities in identifying regions of negative and enhanced fluxes.

Sweep and ejection events in turbulent boundary layer flows have been explored for half a century now to describe eddies impacting turbulent stresses. Yet, moving these studies from their current diagnostic phase to a prognostic form remains a formidable challenge. Here, a cumulant expansion is used to derive a link between the transport of shear stress and the balance of local sweep and ejection events. Cumulant expansion is further used to connect this transport to a metric of asymmetry in the streamwise velocity distribution. These relations are employed to develop two so-called structural models for predicting the turbulent stress transport, which is traditionally neglected in first-order closure of the shear stress budget. Several datasets collected in rough-wall conditions are used to show the importance of the transport term in the roughness sublayer and to demonstrate the predictive skill of the two structural models. The model parameters are invariant to the tested range of Reynolds number and surface roughness, indicating the structural similarity between the velocity asymmetry, sweep/ejection balance, and stress transport may be universal and independent of roughness. Finally, the implementation of the structural models for improved closure schemes of the shear stress budget in modeling applications and wall-modeling in large-eddy simulations are discussed.

Turbulent velocity, temperature, water vapor concentration, and other scalars were measured at the canopyatmosphere interface of a 13-14-m-tall uniform pine forest and a 33-m-tall nounuiform hardwood forest. These measurement were used to investigate whether the mixing layer (ML) analogy of Raupach et al. predicts eddy sizes and now characteristics responsible for much of the turbulent stresses and vertical scalar fluxes. For this purpose, wavelet spectra and cospectra were derived and analyzed. It was found that the MI. analogy predicts well vertical velocity variances and integral timescales. However, at low wavenumbers, inactive eddy motion signatures were present in horizontol velocity wavelet spectra, suggesting that MI. may not be suitable for scaling horizontal velocity perturbations. Momentum and scalar wavelet cospectra of turbulent stresses and scalar fluxes demonstrated that active eddy motion, which was shown by Raupach et al. to be the main energy contributor to vertical velocity (w) spectral energy (Em). is also the main scalar flux-transporting eddy motion. Predictions using ML of the peak E, frequency are in excellent agreement with measured waveled cospectral peaks of vertical fluxes (Kh = 1.5, where K is wavenumber and h is canopy height). Using Lorentz wavelet thresholding of vertical velocity time series, wavelet coefficients associated with active turbulence were identified. It was demonstrated that detection frequency of organized structures, as predicted from Lorentz wavelet filtering, relate to the arrival frequency /h and integral timescale, where is the mean horizontal velocity at height z = h. The newly proposed wavelet thresholding approach, which relies on a"global" wavelet threshold formulation for the energy in w, provides simultaneous energy-covariance-preserving characterization of "active" turbulence at the canopy-atmosphere interface.

Surface flow redistribution on flat ground from crusted bare soil to vegetated patches following intense rainfall events elevates plant available water above that provided by rainfall. The significance of this surface water redistribution to sustaining vegetation in arid and semiarid regions is undisputed. What is disputed is the quantity and spatial distribution of the redistributed water. In ecohydrological models, such nonuniform flows are described using the Saint-Venant equation (SVE) subject to a Manning roughness coefficient closure. To explore these assumptions in the most idealized setting, flume experiments were conducted using rigid cylinders representing rigid vegetation with varying density. Flow was induced along the streamwise x direction by adjusting the free water surface height H(x) between the upstream and downstream boundaries mimicking the nonuniformity encountered in nature. In natural settings, such H(x) variations arise due to contrasts in infiltration capacity and ponded depths during storms. The measured H(x) values in the flume were interpreted using the SVE augmented with progressively elaborate approximations to the roughness representation. The simplest approximation employs a friction factor derived from a drag coefficient (Cd) for isolated cylinders in a locally (but not globally) uniform flow and upscaled using the rod density that was varied across experiments. Comparison between measured and modeled H(x) suggested that such a "naive" approach overpredicts H(x). Blockage was then incorporated into the SVE model calculations but resulted in underestimation of H(x). Biases in modeled H(x) suggest that Cd must be varying in x beyond what a local or bulk Reynolds number predicts. Inferred Cd(x) from the flume experiments exhibited a near-parabolic shape most peaked in the densest canopy cases. The outcome of such Cd(x) variations is then summarized in a bulk resistance formulation that may be beneficial to modeling runon-runoff processes on shallow slopes using SVE.

Momentum and mass exchanges between the atmosphere and forests situated on complex terrain continue to draw significant research attention primarily because of their significance to a plethora of applications. In this paper, we investigated flows behavior on the leeward side of a two-dimensional forested ridge under neutrally stratified conditions using large-eddy simulations (LESs). The goal is to understand how variations in leaf area index (LAI), vertical canopy foliage distributions, and forest edge positions affect mean/turbulent flow statistics, momentum fluxes, and onset of recirculation patterns. Although pressure perturbations are dominated by the hill shape, it is demonstrated here that changes in canopy foliage distribution modulate intensities and patterns of the leeward adverse pressure gradients. Such changes in the adverse pressure gradients alter the mean velocity streamlines including the patterns and magnitudes of the leeward downward mean vertical velocity and the velocity variances and momentum flux in the wake region. While a downwind recirculation zone develops in all cases, the details regarding the incipient location and recirculation zone size vary including positions of the separation and reattachment points. Furthermore, changes in the strength and depth of the zone occur due to canopy-induced changes in adverse pressure gradients, advection, and canopy drag. Because the recirculation zone impacts the local mean advective terms in momentum and scalar exchanges, the simulations here indicate that canopy morphology-induced changes in the leeward flows have significant implications to both measurements and models of biosphere-atmosphere exchange over complex terrain.

Seed dispersal alters gene flow, reproduction, migration and ultimately spatial organization of dryland ecosystems. Because many seeds in drylands lack adaptations for long-distance dispersal, seed transport by secondary processes such as tumbling in the wind or mobilization in overland flow plays a dominant role in determining where seeds ultimately germinate. Here, recent developments in modeling runoff generation in spatially complex dryland ecosystems are reviewed with the aim of proposing improvements to mechanistic modeling of seed dispersal processes. The objective is to develop a physically-based yet operational framework for determining seed dispersal due to surface runoff, a process that has gained recent experimental attention. A Buoyant OBject Coupled Eulerian - Lagrangian Closure model (BOB-CELC) is proposed to represent seed movement in shallow surface flows. The BOB-CELC is then employed to investigate the sensitivity of seed transport to landscape and storm properties and to the spatial configuration of vegetation patches interspersed within bare earth. The potential to simplify seed transport outcomes by considering the limiting behavior of multiple runoff events is briefly considered, as is the potential for developing highly mechanistic, spatially explicit models that link seed transport, vegetation structure and water movement across multiple generations of dryland plants.

Droughts in a warming climate have become more common and more extreme, making understanding forest responses to water stress increasingly pressing. Analysis of water stress in trees has long focused on water potential in xylem and leaves, which influences stomatal closure and water flow through the soil-plant-atmosphere continuum. At the same time, changes of vegetation water content (VWC) are linked to a range of tree responses, including fluxes of water and carbon, mortality, flammability, and more. Unlike water potential, which requires demanding in situ measurements, VWC can be retrieved from remote sensing measurements, particularly at microwave frequencies using radar and radiometry. Here, we highlight key frontiers through which VWC has the potential to significantly increase our understanding of forest responses to water stress. To validate remote sensing observations of VWC at landscape scale and to better relate them to data assimilation model parameters, we introduce an ecosystem-scale analog of the pressure-volume curve, the non-linear relationship between average leaf or branch water potential and water content commonly used in plant hydraulics. The sources of variability in these ecosystem-scale pressure-volume curves and their relationship to forest response to water stress are discussed. We further show to what extent diel, seasonal, and decadal dynamics of VWC reflect variations in different processes relating the tree response to water stress. VWC can also be used for inferring belowground conditions-which are difficult to impossible to observe directly. Lastly, we discuss how a dedicated geostationary spaceborne observational system for VWC, when combined with existing datasets, can capture diel and seasonal water dynamics to advance the science and applications of global forest vulnerability to future droughts.