Determining Stationary Periods across Multiple Sensors: An Application to Observed Canopy Turbulence Response to Atmospheric Stability

A recently proposed multisensor stationarity analysis technique (MSATv1) is improved to eliminate the initial interrogation of time-averaged wind directions, a redundant and potentially biasing procedure for a technique capable of detecting changes in mean wind directions. The new technique, MSATv2, satisfies two basic expectations that are not guaranteed in MSATv1: 1) a nonstationary event should not belong to any stationary interval identified with a given stringency, and 2) nonstationary events identified with an arbitrary stringency should continue to be identified as nonstationary with increasing stringency. These expectations are confirmed by applying MSATv2 to two long periods, during the defoliated phase of the Canopy Horizontal Array Turbulence Study (CHATS), whose durations are determined solely by data availability. MSATv2 successfully determines visually trivial and nontrivial nonstationary transitions, uncovering details of the time evolution of dynamic processes. MSATv2 yields ensemble-average estimates of mean wind speeds and directions with well-controlled and quantifiable uncertainties for atmospheric stability conditions ranging from near neutral to free convection. These results enable interrogation of the observed canopy turbulence response to atmospheric stability in isolation from contamination by spatial variation with position relative to canopy elements. MSATv2 results also reveal the connection between the presence of organized convective structures and variability in mean shear, showing the role of organized convective structures in the observed relationship between the bulk drag coefficient and atmospheric instability.

Revisiting Kelvin Helmholtz Instabilities and von Kármán Vortices in Canopy Turbulence

Studying turbulence in vegetation canopies is important in the context of a number of micrometeorological and hydrological applications. While recent focus has shifted more towards exploring different kinds of canopy heterogeneities, there are still gaps in the existing knowledge on the multiple types of dynamics involved in the case of horizontally homogeneous canopies. For example, experimental studies have indicated that turbulence in the canopy sublayer (CSL) can be divided into three regimes. In the deep-zone, the flow-field is dominated by von Kármán vortex streets and interrupted by strong sweep events. The second zone near the canopy top is dominated by attached eddies and Kelvin-Helmholtz waves associated with the velocity inflection point in the mean longitudinal velocity profile. Above the canopy, the flow resembles classical boundary layer flow. In this study, these different kinds of dynamics are studied together by means of a large eddy simulation (LES). The main theme of this work is to address the question whether the parametrization of the canopy by a distributed drag force in numerical simulations instead of placing real solid obstacles is consistent with the three layer conceptual model. Unique techniques such as measures from information theory and coupled oscillator analysis are used to extract the coherent structures associated with the two motions. It can be stated that a better understanding of the rich dynamics associated with the simplest case of canopy turbulence can lead to more efficient simulations and more importantly improve the interpretation of more complex scenarios.

Explaining the convector effect in canopy turbulence by means of large-eddy simulation

Semi-arid forests are found to sustain a massive sensible heat flux in spite of having a low surface to air temperature difference by lowering the aerodynamic resistance to heat transfer (rH) – a property called the canopy convector effect (CCE). In this work large-eddy simulations are used to demonstrate that the CCE appears more generally in canopy turbulence. It is indeed a generic feature of canopy turbulence: rH of a canopy is found to reduce with increasing unstable stratification, which effectively increases the aerodynamic roughness for the same physical roughness of the canopy. This relation offers a sufficient condition to construct a general description of the CCE. In addition, we review existing parameterizations for rH from the evapotranspiration literature and test to what extent they are able to capture the CCE, thereby exploring the possibility of an improved parameterization.