Estimating the Instantaneous Drag–Wind Relationship for a Horizontally Homogeneous Canopy

Abstract

The mean drag–wind relationship is usually investigated assuming that field data are representative of spatially-averaged metrics of statistically stationary flow within and above a horizontally homogeneous canopy. Even if these conditions are satisfied, large-eddy simulation (LES) data suggest two major issues in the analysis of observational data. Firstly, the streamwise mean pressure gradient is usually neglected in the analysis of data from terrestrial canopies, which compromises the estimates of mean canopy drag and provides misleading information for the dependence of local mean drag coefficients on local velocity scales. Secondly, no standard approach has been proposed to investigate the instantaneous drag–wind relationship, a critical component of canopy representation in LES. Here, a practical approach is proposed to fit the streamwise mean pressure gradient using observed profiles of the mean vertical momentum flux within the canopy. Inclusion of the fitted mean pressure gradient enables reliable estimates of the mean drag–wind relationship. LES data show that a local mean drag coefficient that characterizes the relationship between mean canopy drag and the velocity scale associated with total kinetic energy can be used to identify the dependence of the local instantaneous drag coefficient on instantaneous velocity. Iterative approaches are proposed to fit specific models of velocity-dependent instantaneous drag coefficients that represent the effects of viscous drag and the reconfiguration of flexible canopy elements. LES data are used to verify the assumptions and algorithms employed by these new approaches. The relationship between mean canopy drag and mean velocity, which is needed in models based on the Reynolds-averaged Navier-Stokes equations, is parametrized to account for both the dependence on velocity and the contribution from velocity variances. Finally, velocity-dependent drag coefficients lead to significant variations of the calculated displacement height and roughness length with wind speed.