Laserfiche WebLink
<br /> <br />472 <br /> <br />Topographic Response of a Oar <br /> <br />, <br /> <br />~ 1 [-->- -->-] <br />or = - Cb V . qs <br /> <br />(2) <br /> <br />where '1 is the bed elevation, q; is the volumetric vector sediment flux and ~ <br /> <br />=....65 is the percentage of the bed made up of sediment (Le., unity minus the <br />porosity). <br />Using computed sediment fluxes, the rate of erosion is easily calculated from <br />Equation (2). The erosion rate, in conjunction with an assumed time increment, <br />predicts the modification to the bed in that interval. Thus, this equation may be <br />employed in an iterative manner to predict temporal adjustment of the bed to flow <br />conditions, as discussed below. <br />The other way in which the channel may be modified is through changes in <br />stage. Although this effect does not modify the actual bed topography (except <br />through erosion and deposition, as discussed in the previous paragraph), it does have <br />a profound effect on the planfonn and depth of the channel. For example, the <br />centerline trace of a typical meandering stream is typically more sinuous at low flow <br />than at bankfull; the centerline is defined as the locus of points midway between the <br />two water edges. Similarly, the relative cross-sectional area at a given longitudinal <br />location in the reach may change dramatically from low to high flow. TheSe effects <br />may seem obvious or even trivial, but they play a major role in determining the <br />difference in patterns of erosion and deposition at different stages in a given reach. <br />The method used in developing low flow initial channels from bankfull <br />topographic surveys is relatively straightforward. First, a comprehensive <br />topographic map is constructed from surveyed points. In order to obtain <br />approximate channel alignment and topography below bankfull elevation, the <br />bankfuII channel is "sliced" along a surface parallel to the bankfull free surface. <br />This surface is positioned at a distance below the bankfull free surface equivalent to <br />the desired stage drop. Points below the plane are interpolated using a tensioned <br />spline and resampled in a regular grid. The potential error in this approach arises <br />from assuming that the pattern of superelevation in the free surface is relatively <br />invariant with stage change. This approximation is never precisely correct, but the <br />error incurred in the depths is typically very small. Furthermore, the size of the <br />error can be calculated using the flow model, and a new free surface can be <br />calculated iteratively to improve the estimated depth. For the Green River reach <br />desj:ribed herein, the error was typically less than two or three centimeters. This is <br />well ~ithin the uncertainty in the surveyed topography, and was assumed to be <br />negligible. <br />The technique described above also can be employed to predict variations in <br />channel form due to stage change while erosion and deposition are taking place on <br />the bed. In other words, the channel topography can be adjusted with each <br />increment of time to account for both erosion and deposition and stage changes in a <br />quasi~teady sense. This allows the investigation of the response of the channel <br />topography in cases where the change in stage occurs over time scales which are <br />short compared to the bed modification and in cases where discharge and stage vary <br />slowly compared to the bed evolution. <br /> <br />The Evolution Model <br /> <br />The coupling of the various components of the full evolution model is shown in <br />Figure 3. As shown, the inputs to the evolution model consist of the following: the <br />initial planform and topography of the channel, the sediment size making up the <br />bed, and the sta~e and discharge as a function of time. Using the specified channel <br />geometry and iDltial discharge, the flow model described herein can be employed to <br />calculate the spatial distribution of boundary shear stress. These values are inserted <br />