Pool Formation in Boulder-Bed Streams: Implications From 1-D and 2-D Numerical Modeling

In mountain rivers of Southern California, boulder-large roughness elements strongly influence flow hydraulics and pool formation and maintenance. In these systems, boulders appear to control the stream morphology by converging flow and producing deep pools during channel forming discharges. Our research goal is to develop quantitative relationships between boulder roughness elements, temporal patterns of scour and fill, and geomorphic processes that are important in producing pool habitat. The longitudinal distribution of shear stress, unit stream power and velocity were estimated along a 48 m reach on Rattlesnake Creek, using the HEC-RAS v 3.0 and River 2-D numerical models. The reach has an average slope of 0.02 and consists of a pool-riffle sequence with a large boulder constriction directly above the pool. Model runs were performed for a range of stream discharges to test if scour and fill thresholds for pool and riffle environments could be identified. Results from the HEC-RAS simulations identified that thresholds in shear stress, unit stream power and mean velocity occur above a discharge of 5.0 cms. Results from the one-dimensional analysis suggest that the reversal in competency is likely due to changes in cross-sectional width at varying flows. River 2-D predictions indicated that strong transverse velocity gradients were present through the pool at higher modeled discharges. At a flow of 0.5 cms (roughly 1/10th bankfull discharge), velocities are estimated at 0.6 m/s and 1.3 m/s for the pool and riffle, respectively. During discharges of 5.15 cms (approximate bankfull discharge), the maximum velocity in the pool center increased to nearly 3.0 m/s, while the maximum velocity over the riffle is estimated at approximately 2.5 cms. These results are consistent with those predicted by HEC-RAS, though the reversal appears to be limited to a narrow jet that occurs through the pool head and pool center. Model predictions suggest that the velocity reversal is produced by a boulder-bedrock constriction that rapidly decreases the channel width above the pool by roughly 25 percent. The width constriction creates highly turbulent flow capable of scouring bed material through the pool. The high velocity core that is produced through the pool center appears to be enhanced by the formation of a large eddy directly below the boulder. Values of unit stream power and shear stress indicate that the pool exit is an area of deposition of bed material due to a decrease in tractive force. The presence of a strong transverse velocity gradient suggests that only a portion of the flow is responsible for scouring bed material. After we eliminate the dead water zone, the lowest five percent of the velocity range, patterns of effective width between pools and riffles begin to emerge. The ratio of flow width between adjacent pools and riffles is one measure of flow convergence. At a discharge of 0.5 cms, the ratio of effective width between pools and riffles is roughly 1:1, implying that there is uniform flow with little flow convergence. At a discharge of 5.15 cms the width ratio between the pool and riffle is about 1:3, demonstrating the strong convergent flow patterns at the pool head. The observed effective width relationship suggests that when considering restoration designs, boulders should be placed in areas that replicate natural convergence and divergence patterns in order to maximize pool area and depth.