Steady nonuniform shallow flow within emergent vegetation
- Author(s): Wang, WJ
- Huai, WX
- Thompson, S
- Katul, GG
- et al.
Published Web Locationhttps://doi.org/10.1002/2015WR017658
© 2015. American Geophysical Union. All Rights Reserved. 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.