Richard R. McDonald

Modeling the Effect of Flow and Sediment Transport on White Sturgeon Spawning Habitat in the Kootenai River, Idaho

Kootenai River white sturgeon spawn in an 18-km reach of the Kootenai River, Id. Since completion of Libby Dam upstream from the spawning reach in 1972, 1974 is the only year with documented significant recruitment of juvenile fish. Where successful in other rivers, white sturgeon spawn over clean coarse material of gravel size or larger. The channel substrate in the current (2008) 18-km spawning reach is composed primarily of sand and some buried gravel; within a few kilometers upstream there is an extended reach of clean gravel, cobble, and bedrock. We used a quasi-three-dimensional flow and sediment-transport model along with the locations of collected sturgeon eggs as a proxy for spawning location from 1994 to 2002 to gain insight into spawning-habitat selection in a reach which is currently unsuitable due to the lack of coarse substrate. Spatial correlations between spawning locations and simulated velocity and depth indicate fish select regions of higher velocity and greater depth within any river cross section to spawn. These regions of high velocity and depth occur in the same locations regardless of the discharge magnitude as modeled over a range of pre- and postdam flow conditions. A flow and sediment-transport simulation shows high discharge, and relatively long-duration flow associated with predam flow events is sufficient to scour the fine sediment overburden, periodically exposing existing lenses of gravel and cobble as lag deposits in the current spawning reach. This is corroborated by video observations of bed surface material following a significant flood event in 2006, which show gravel and cobble present in many locations in the current spawning reach. Thus, both modeling and observations suggest that the relative rarity of extremely high flows in the current regulated flow regime is at least partly responsible for the lack of successful spawning; in the predam flow regime, frequent high flows removed the fine sediment overburden, unveiling coarse material and providing suitable substrate in the current spawning reach.

Nonperiodic eddy pulsations

Recirculating flow in lateral separation eddies is typically weaker than main stem flow and provides an effective environment for trapping sediment. Observations of recirculating flow and sedimentary structures demonstrate that eddies pulsate in size and in flow velocity even when main stem flow is steady. Time series measurements of flow velocity and location of the reattachment point indicate that these pulsations are nonperiodic. Nonperiodic flow in the lee of a channel margin constriction is grossly different from the periodic flow in the lee of a cylinder that is isolated in a flow. Our experiments demonstrate that placing a flow-parallel plate adjacent to a cylinder is sufficient to cause the leeside flow to change from a periodic sequence of vortices to a nonperiodically pulsating lateral separation eddy, even if flow conditions are otherwise unchanged. Two processes cause the leeside flow to become nonperiodic when the plate is added. First, vortices that are shed from the cylinder deform and become irregular as they impact the plate or interfere with remnants of other vortices near the reattachment point. Second, these deformed vortices and other flow structures are recirculated in the lateral separation eddy, thereby influencing the future state (pressure and momentum distribution) of the recirculating flow. The vortex deformation process was confirmed experimentally by documenting spatial differences in leeside flow; vortex shedding that is evident near the separation point is undetectable near the reattachment point. Nonlinear forecasting techniques were used in an attempt to distinguish among several possible kinds of nonperiodic flows. The computational techniques were unable to demonstrate that any of the nonperiodic flows result from low-dimensional nonlinear processes.

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Estimating floodwater depths from flood inundation maps and topography

Remote sensing analysis is routinely used to map flood inundation during flooding events or retrospectively for planning and research activities. Quantification of the depth of floodwater is important for emergency response, relief operations, damage assessment etc. The Floodwater Depth Estimation Tool (FwDET) calculates water depth based on topographic analysis using standard GIS tools within a Python script. FwDET’s low input requirements (DEM and inundation polygon) and high computational efficiency lend it as a useful tool for emergency response and large-scale applications. Operational use of FwDET is described herein as part of emergency response activation of the Global Flood Partnership (GFP) during the 2017 USA Hurricane Season and May 2018 flooding in Sri Lanka. Use of FwDET during Hurricanes Harvey (Texas and Louisiana), Irma (Florida) and Maria (Puerto Rico) demonstrated its utility by producing large-scale water depth products at near-real-time at relatively high spatial resolution. Despite FwDET’s success, limitations of the tool stemmed from bureaucratic disallowance of non-governmental remote sensing products by U.S. federal emergency response agencies, misclassified remotely sensed floodwaters and challenges obtaining global high resolution DEMs specifically for the aforementioned Sri Lankian flooding. While global-scale DEM products at 30m resolution are freely available, these datasets are of integer precision and thus have limited vertical resolution. This limitation is significant primarily in flat (e.g. coastal) locations and flooded domains comprised of relatively small patches of water.