Prediction of Total Sediment Load in Sand‐Bed Rivers in Korea Using Lateral Distribution Method

A new method for numerically predicting the total sediment load in a river is proposed. The method can be used to predict the total sediment load with information on channel geometry and slope, flow, and bed materials. The conventional method uses a 1D approach that assumes the channel has a wide rectangular shape. However, the proposed method computes depth‐averaged velocity over the width and predicts the total sediment load based on the flow computations. The new method, therefore, is expected to predict better if the flow changes significantly in the lateral direction. The proposed method was applied to three large sand‐bed rivers in Korea, where information is available regarding suspended sediment. Five formulas were tested of use in making total sediment load computations, namely Engelund‐Hansen's, Ackers‐White's, Yang's, Brownlie's, and Karim's formulas. The predicted total sediment loads are compared not only with measured data but also with results calculated using the 1D approach. Discrepancy ratios between the predicted and measured total sediment loads are given and the results are discussed.     

HYSTAR Sediment Model: Distributed Two‐Dimensional Simulation of Watershed Erosion and Sediment Transport Using Time‐Area Routing

An erosion and sediment transport component incorporated in the HYdrology Simulation using Time‐ARea method (HYSTAR) upland watershed model provides grid‐based prediction of erosion, transport and deposition of sediment in a dynamic, continuous, and fully distributed framework. The model represents the spatiotemporally varied flow in sediment transport simulation by coupling the time‐area routing method and sediment transport capacity approach within a grid‐based spatial data model. This avoids the common, and simplistic, approach of using the Universal Soil Loss Equation (USLE) to estimate erosion rates with a delivery ratio to relate gross soil erosion to sediment yield of a watershed, while enabling us to simulate two‐dimensional sediment transport processes without the complexity of numerical solution of the partial differential governing equations. In using the time‐area method for routing sediment, the model offers a novel alternative to watershed‐scale sediment transport simulation that provides detailed spatial representation. In predicting four‐year sediment hydrographs of a watershed in Virginia, the model provided good performance with R² of 0.82 and 0.78 and relative error of ?35% and 11% using the Yalin and Yang's sediment transport capacity equations, respectively. Prediction of spatiotemporal variation in sediment transport processes was evaluated using maps of sediment transport rates, concentrations, and erosion and deposition mass, which compare well with expected behavior of flow hydraulics and sediment transport processes.     

Comparison of Unsteady and Quasi‐Unsteady Flow Models in Simulating Sediment Transport in an Ephemeral Arizona Stream1

Hummel, Ryan, Jennifer G. Duan, and Shiyan Zhang, 2012. Comparison of Unsteady and Quasi‐Unsteady Flow Models in Simulating Sediment Transport in an Ephemeral Arizona Stream. Journal of the American Water Resources Association (JAWRA) 48(5): 987‐998. DOI: 10.1111/j.1752‐1688.2012.00663.x Abstract: Hydrodynamic and sediment transport models are useful engineering tools for predicting unsteady flood flow and sediment transport. Many models such as HEC‐RAS, HEC‐6, and IALLUVIAL apply quasi‐unsteady flow model, whereas others apply the unsteady flow model. It remains unknown if a quasi‐unsteady flow model is sufficiently accurate for simulating sediment transport in rapidly varied unsteady flood events, especially in ephemeral rivers in arid and semiarid regions. This study compared the quasi‐unsteady HEC‐RAS 4.1 model with one‐dimensional (1D) Finite Volume Method (FVM) based model in simulating flood flow and sediment transport in the Pantano Wash, a dryland river in the state of Arizona. The objective is to determine which sediment transport method is appropriate in predicting bed elevation changes in an ephemeral stream, Pantano Wash, and if an unsteady model is more accurate than a quasi‐unsteady flow model in predicting sediment transport. Results showed that the quasi‐unsteady HEC‐RAS model and the 1D FVM yielded similar results of bed degradation and aggradation for this dryland stream, although the FVM model predicted better flood hydrographs. Among the seven sediment transport formulas embedded in HEC‐RAS, Yang’s and Engelund‐Hansen’s equations gave the best matches with the field measurements for this particular case study.     

EFFECT OF FLOW VELOCITY ON SEDIMENT OXYGEN DEMAND: COMPARISON OF THEORY AND EXPERIMENTS1

ABSTRACT: Theoretical equations that establish the relationship between sediment oxygen demand (SOD) in a lake and the flow velocity and dissolved oxygen concentration in the bulk water already exist. These theoretical equations for oxygen consumption in the sediment express biological consumption with Michaelis-Menten kinetics, and chemical consumption by a first order reaction. Data from laboratory experiments that were conducted to validate the theoretical equations also exist. These experiments were performed in a laboratory channel with well defined flow characteristics for three types of sediments. Herein, the theoretical equations are used to model the experimental data for the three types of sediments. The values used for the parameters in the theoretical equations are determined by iteration until a best fit is obtained for the relationship of SOD to flow velocity from both the theoretical model and experimental data. The goodness of fit is measured by the standard error of prediction and the regression coefficient.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.     

INVESTIGATION OF TURBIDITY MAXIMUM IN A MESOTIDAL ESTUARY,TAIWAN1

ABSTRACT: To comprehend the distributions of salinity, temperature, and suspended sediment in the Danshuei River estuary in Taiwan, monthly field surveys were conducted in 2003. These included several high and low slackwater surveys and intensive surveys. The results show that the Danshuei River estuary is predominately a partially mixed estuary. The highest concentration of suspended sediment is typically observed at the Chung‐Hsin Bridge, the most upstream sampling station. The suspended sediment concentration exhibits a general decreasing trend in the downstream direction. It may be concluded that the sediments mostly come from the upstream reach. A locally high concentration of suspended sediment is found at the Kuan‐Du station because of the local deep channel bathymetry and two‐layered estuarine circulation. A vertical two‐dimensional hydrodynamic and sediment transport model is applied to investigate the tidally averaged salinity distribution, residual circulation, and suspended sediment concentration. The modeling results reveal that, under the Q₇₅ flow condition (i.e., low flow), a turbidity maximum occurs at the Kuan‐Du station due to the strong estuarine circulation. The model simulation with a much higher river flow condition results in a weaker residual circulation and weaker turbidity maximum.     

Streamflow,Sediment Transport,and Geomorphic Change during the 2011 Flood on the Missouri River Near Bismarck–Mandan,ND

Geomorphic change from extreme events in large managed rivers has implications for river management. A steady‐state, quasi‐three‐dimensional hydrodynamic model was applied to a 29‐km reach of the Missouri River using 2011 flood data. Model results for an extreme flow (500‐year recurrence interval [RI]) and an elevated managed flow (75‐year RI) were used to assess sediment mobility through examination of the spatial distribution of boundary or bed shear stress (τ_b) and longitudinal patterns of average τ_b, velocity, and kurtosis of τ_b. Kurtosis of τ_b was used as an indicator of planform channel complexity and can be applied to other river systems. From differences in longitudinal patterns of sediment mobility for the two flows we can infer: (1) under extreme flow, the channel behaves as a single‐thread channel controlled primarily by flow, which enhances the meander pattern; (2) under elevated managed flows, the channel behaves as multithread channel controlled by the interaction of flow with bed and channel topography, resulting in a more complex channel; and (3) for both flows, the model reach lacks a consistent pattern of deposition or erosion, which indicates migration of areas of erosion and deposition within the reach. Despite caveats and limitations, the analysis provides useful information about geomorphic change under extreme flow and potential implications for river management. Although a 500‐year RI is rare, extreme hydrologic events such as this are predicted to increase in frequency.     

Impact of Temporal Resolution of Flow‐Duration Curve on Sediment Load Estimation1

Chen, Li, Rina Schumer, Anna Knust, and William Forsee, 2011. Impact of Temporal Resolution of Flow‐Duration Curve on Sediment Load Estimation. Journal of the American Water Resources Association (JAWRA) 48(1): 145‐155. DOI: 10.1111/j.1752‐1688.2011.00602.x Abstract: Estimates of a channel’s annual sediment transport capacity typically incorporate annualized flow‐duration curves. Average daily flow data, commonly used to develop flow‐duration curves, may not adequately describe sediment‐transporting flows in arid and semiarid ephemeral streams. In this study, we examined impacts of varied temporal resolution flow data on annual sediment load estimation. We derived flow‐duration curves for eight sites in the Southwestern United States based on both 15‐min and daily‐averaged flow data. We then estimated sediment loads for both flow‐duration curves using the Sediment Impact Analysis Method, implemented in HEC‐RAS. When average daily flow is used to generate flow‐duration curves, sediment load estimation is lower by up to an order of magnitude. This trend is generally unaffected by uncertainty associated with sediment particle size or hydraulic roughness. The ratio of sediment loads estimated by 15‐min versus daily‐averaged flow‐duration curves is strongly correlated with channel slope, being greater on steep‐slope channels. Sediment loads estimated by the two types of flow‐duration curves are closely correlated, suggesting possible relationships for improving predictions when high‐temporal resolution data are unavailable. Results also suggest that the largest flow contributes significantly to total sediment load, and thus will greatly impact ephemeral stream geomorphology in arid and semiarid regions.     

Flow,Organic, and Inorganic Sediment Yields from a Channelized Watershed in the South Carolina Lower Coastal Plain

Many small streams in coastal watersheds in the southeastern United States are modified for agricultural, residential, and commercial development. In the South Carolina Lower Coastal Plain, low‐relief topography and a shallow water table make stream channelization ubiquitous. To quantify the impacts of urbanization and stream channelization, we measured flow and sediment from an urbanizing watershed and a small forested watershed. Flow and sediment export rates were used to infer specific yields from forested and nonforested regions of the urbanizing watershed. Study objectives were to: (1) quantify the range of runoff‐to‐rainfall ratios; (2) quantify the range of specific sediment yields; (3) characterize the quantity and quality of particulate matter exported; and (4) estimate sediment yield attributable to agriculture, development, and channelization activities in the urbanizing watershed. Our results showed that the urban watershed exported over five times more sediment per unit area compared with the forested watershed. Sediment concentration was related to flow flashiness in the urban watershed and to flow magnitude in the forested watershed. Sediments from the forested watershed were dominated by organic matter, whereas mineral matter dominated sediment from the urban stream. Our results indicated that a significant shift in sediment quality and quantity are likely to occur as forested watersheds are transformed by urbanization in coastal South Carolina.     

A Depth‐Proportional Intake Device for Automatic Water Samplers1

Abstract: This paper describes the construction and testing of a device for pumping water samplers that collects suspended sediment samples by moving the intake vertically to keep it at the same proportion of flow depth. The device uses a simple sprocket mechanism that can be mounted vertically on the downstream side of culverts and bridge pilings to protect against damage from floating debris during storms. Suspended sediment samples collected from an urban stream with the depth‐proportional device were compared with manual samples taken with a depth‐integrated sampler. Scatter in the relationship between pumped and manual samples (R² = 0.76) are probably explained by horizontal variability in concentrations, poor mixing associated with lateral sediment inputs from construction site erosion, the downstream orientation of the intake, and the failure of the concentration at 60% of the flow depth to match the average vertical concentration.     

Modeling of Vertical Solute Dispersion in a Sediment Bed Enhanced by Wave‐Induced Interstitial Flow1

Abstract: Mass (solute) transport in a stream or lake sediment bed has a significant effect on chemical mass balances and microbial activities in the sediment. A “1D vertical dispersion model” is a useful tool to analyze or model solute transfer between river or lake water and a sediment bed. Under a motionless water column, solute transfer into and within the sediment bed is by molecular diffusion. However, surface waves or bed forms create periodic pressure waves along the sediment/water interface, which in turn induce flows in the pores of the sediment bed. The enhancement of solute transport by these interstitial periodic flows in the pores has been incorporated in a 1D depth‐dependent “enhanced dispersion coefficient (D_E).” Typically, D_E diminishes exponentially with depth in the sediment bed. Relationships have been developed to estimate D_E as a function of the characteristics of sediment (particle size, hydraulic conductivity, and porosity) and pressure waves (wave length and height). In this paper, we outline and illustrate the calculation of D_E as well as the penetration depth (d_p) of the flow effect. Sample applications to illustrate the computational procedure are provided for dissolved oxygen transfer into a stream gravel bed and release of phosphorus from a lake bed. The sensitivity of the results to input parameter values is illustrated, and compared with the errors obtained when interstitial flow is ignored. Maximum values of D_E near the sediment surface can be on the order of 1 cm²/s in a stream gravel bed with standing waves, and 0.001 cm²/s in a fine sand lake bed under progressive surface waves, much larger than molecular diffusion coefficients.     

DISCHARGE AND SUSPENDED SEDIMENT PATTERNS OF AN INTERMITTENT COLD DESERT STREAM1

Sage Creek in south‐central Wyoming is listed as impaired by the U.S. Environmental Protection Agency (USEPA) due to its sediment contribution to the North Platte River. Despite the magnitude of sediment impacts on streams, little research has been conducted to characterize patterns of sediment transport or to model suspended sediment concentration in many arid western U.S. streams. This study examined the relationship between stream discharge and suspended sediment concentration near the Sage Creek and North Platte River confluence from 1998 through 2003. The objectives were to determine patterns of stream discharge and suspended sediment concentration, produce a sediment prediction model, and compare sediment concentrations for the six‐year period. Stream discharge and suspended sediment transport responded rapidly to convective storms and spring runoff events. During the study period, events exceeding 0.23 m³/s accounted for 92 percent of the sediment load, which is believed to originate from erodible headwater uplands. Further analysis of these data indicates that time series modeling is superior to simple linear regression in predicting sediment concentration. Significant increases in suspended sediment concentration occurred in all years except 2003. This analysis suggests that a six‐year monitoring record was insufficient to factor out impacts from climate, geology, and historical sediment storage.     

HYDRAULIC GEOMETRY AND STREAM CHANNEL BEHAVIOR: A UNCERTAIN LINK1

ABSTRACT: The at-a-station hydraulic geometry of stream channels can serve as a predictor of alluvial stream channel behavior. This geometry is the empirical relations describing changes in water surface width, mean depth, and mean velocity with changing discharge. The exponent values are correlated with channel morphology and behavior such as scour and fill, flow resistance, bank resistance, and competence. Channel behavior and morphology are apparently related, but some causes for effects are uncertain. Several studies, using empirical and theoretical bases, are reviewed here to illustrate the relation between hydraulic geometry and channel behavior, but the relations are not always consistent. Hydraulic geometry variables are easy to measure and readily available, but they do not always reflect what may be more important ones such as turbulence, the velocity distribution profile, and distribution and cohesion of sediment particles. This paper illustrates some of these problems, provides some solutions, and addresses need for more work to better predict stream channel behavior from hydraulic geometry     

Concentration‐Duration‐Frequency Curves for Stream Turbidity: Possibilities for Assessing Biological Impairment1

Abstract: Siltation and subsequent biological impairment is a national problem prompting state regulatory agencies to develop sediment total maximum daily loads (TMDL) for many streams. To support TMDL targets for reduced sediment yield in disturbed watersheds, a critical need exists for stream assessments to identify threshold concentrations of suspended sediment that impact aquatic biota. Because of the episodic nature of stream sediment transport, thresholds should not only be a function of sediment concentration, but also of duration and dose frequency. Water quality sondes can collect voluminous amounts of turbidity data, a surrogate for suspended sediment, at intervals that can be used to characterize concentration, duration, and frequency of elevated turbidity events. To characterize turbidity sonde data in an ecologically relevant manner, a methodology for concentration‐duration‐frequency (CDF) curves was developed using turbidity doses that relate to different levels of biological impairment. To illustrate this methodology, turbidity CDF curves were generated for two sites on Little Pigeon River in the Great Smoky Mountains National Park, Tennessee, using over 30,000 sonde data measurements per site for a one‐year period. Utilizing a Poisson arrival approach, turbidity spikes were analyzed stochastically by observing the frequency and duration of recorded events over a turbidity level that relates to a biological dose response. An exponential equation was used to fit duration and frequency of a specified turbidity level to generate concentric‐shaped CDF curves, where at specific turbidities longer durations occurred less frequently and conversely shorter durations occurred more frequently. The significance of the equation fit to the data was accomplished with a Kolmogorov‐Smirnov goodness‐of‐fit test. Our findings showed that the CDF curves derived by an exponential function performed reasonable well, with most curves significant at a 95% confidence level. These CDF curves were then used to demonstrate how they could be used to assess biological impairment, and identify future research needs for improved development of sediment TMDLs.     

DELIVERY OF SUSPENDED SEDIMENT AND POLLUTANTS FROM NONPOINT SOURCES DURING OVERLAND FLOW1

ABSTRACT: Delivery of sediment and particulate pollutants from diffuse sources is shown to be related to the loss of sediment carrying energy of runoff during the overland flow phase. The loss is caused by the termination of rainfall and by reduction of flow energy during the recession phase of the overland flow hydrograph. It has been demonstrated both by theoretical analyses and experimental measurements that the saturated sediment concentration in overland flow is a function of rainfall erosivity and the runoff flow rate. The hypotheses were verified by field measurements from a small homogeneous watershed.     

Magnitude,Frequency, and Duration Relations for Suspended Sediment in Stable (“Reference”) Southeastern Streams1

Abstract: Sediment is listed as one of the leading causes of water‐quality impairments in surface waters of the United States (U.S.). A water body becomes listed by a State, Territory or Tribe if its designated use is not being attained (i.e., impaired). In many cases, the prescribed designated use is aquatic health or habitat, indicating that total maximum daily loads (TMDL) targets for sediment should be functionally related to this use. TMDL targets for sediment transport have been developed for many ecoregions over the past several years using suspended‐sediment yield as a metric. Target values were based on data from “reference” streams or reaches, defined as those exhibiting geomorphic characteristics of equilibrium. This approach has proved useful to some states developing TMDLs for suspended sediment, although one cannot conclude that if a stream exceeds the target range, the aquatic ecosystem will be adversely impacted. To address this problem, historical flow‐transport and sediment‐transport data from hundreds of sites in the Southeastern U.S. were re‐examined to develop parameters (metrics) such as frequency and duration of sediment concentrations. Sites determined as geomorphically stable from field evaluations and from analysis of gauging‐station records were sorted by ecoregion. Mean‐daily flow data obtained from the U.S. Geological Survey were applied to sediment‐transport rating relations to determine suspended‐sediment load for each day of record. The frequency and duration that a given concentration was equaled or exceeded were then calculated to produce a frequency distribution for each site. “Reference” distributions were created using the stable sites in each ecoregion by averaging all of the distributions at specified exceedance intervals. As with the “reference” suspended‐sediment yields, there is a broad range of frequency and duration distributions that reflects the hydrologic and sediment‐transport regimes of the ecoregions. Ecoregions such as the Mississippi Valley Loess Plains (#74) maintain high suspended sediment concentrations for extended periods, whereas coastal plain ecoregions (#63 and 75) show much lower concentrations.     

Fine‐Sediment Loadings to Lake Tahoe1

Abstract: Over the past 35 years, a trend of decreasing water clarity has been documented in Lake Tahoe, attributable in part to the delivery of fine‐grained sediments emanating from upland and channel sources. The overall objective of the research reported here was to determine the amount of fine sediment delivered to Lake Tahoe from each of the 63 contributing watersheds. The research described in this report used combinations of field‐based observations of channel and bank stability with measured and simulated data on fine‐sediment loadings to estimate fine‐sediment loadings from unmonitored basins throughout the Lake Tahoe Basin. Loadings were expressed in the conventional format of mass per unit time but also in the number of particles finer than 20 μm, the latter for future use in a lake‐clarity model. The greatest contributors of fine sediment happened to be those with measured data, not requiring extrapolation. In descending order, they are as follows: Upper Truckee River [1,010 tonnes per year (T/year)], Blackwood Creek (846 T/year), Trout Creek (462 T/year), and Ward Creek (412 T/year). Summing estimated values from the contributing watersheds provided an average, annual estimate of fine‐sediment (<0.063 mm) loadings to the lake of 5,206 T/year. A total of 7.79E + 19 particles in the 5‐20 μm fraction were calculated to enter Lake Tahoe in an average year with the Upper Truckee River accounting for almost 25% of the total. Contributions from Blackwood, Ward, Trout, and Third creeks account for another 23% of these very fine particles. Thus, these five streams making up about 40% of the basin area, account for almost 50% of all fine‐sediment loadings to the lake. Contribution of fine sediment from streambank erosion were estimated by developing empirical relations between measured or simulated bank‐erosion rates with a field‐based measure of the extent of bank instability along given streams. An average, annual fine‐sediment loading from streambank erosion of 1,305 T/year was calculated. This represents about 25% of the average, annual fine‐sediment load delivered to the lake from all sources. The two largest contributors, the Upper Truckee River (639 T/year) and Blackwood Creek (431 T/year), account for slightly more than 80% of all fines emanating from streambanks, representing about 20% of the fine sediment delivered to Lake Tahoe from all sources. Extrapolations of fine‐sediment loadings to the unmonitored watersheds are based on documented empirical relations, yet contain a significant amount of uncertainty. Except for those values derived directly from measured data, reported results should be considered as estimates.     

THE LODGING VELOCITY FOR EMERGENT AQUATIC PLANTS IN OPEN CHANNELS1

ABSTRACT: The growth of aquatic plants in open‐channels has many adverse environmental effects including, but not limited to, impeding the transport of water, hindering navigation, increasing flood elevations, increasing sediment deposition, and degrading water quality. Existing control strategies include physical removal and chemical treatment. Physical removal is only a temporary solution and chemical treatment is unacceptable if the water will be consumed by humans. The hydrodynamic method can wash out the encroached aquatic plants by keeping flow velocity higher than the critical velocity required to bend and rupture (lodge) their stems. This approach is a promising, physically‐based, efficient, economic, and permanent solution for this problem. However, the success of this approach requires the accurate prediction of the critical lodging velocity. This paper presents an analytic study of the lodging velocity for the submerged portion of aquatic plants of narrow leaved emergent stems that are wider at bottom than the top. Based on the principles of engineering materials and the theory of turbulent flow, a semi‐empirical formula is derived for the prediction of the critical lodging velocity. It indicates that the lodging of aquatic plants is controlled not only by flow conditions but also the geometric and mechanical characteristics of the plants. These analytic results provide a satisfactory explanation of the lodging phenomena observed in the field and are verified by the available experimental data.     

Road Erosion and Delivery Index (READI): A Model for Evaluating Unpaved Road Erosion and Stream Sediment Delivery

The Road Erosion and Delivery Index (READI) is a new geographic information system–based model to assess erosion and delivery of water and sediment from unpaved road networks to streams. READI quantifies the effectiveness of existing road surfacing and drain placements in reducing road sediment delivery and guides upgrades to optimize future reductions. Roads are draped on a digital elevation model and parsed into hydrologically distinct segments. Segments are further divided by engineered drainage structures. For each segment, a kinematic wave approximation generates runoff hydrographs for specified storms, with discharge directly to streams at road–stream crossings and onto overland‐flow plumes at other discharge points. Plumes are attenuated by soil infiltration, which limits their length, with delivery occurring if plumes intersect streams. Sediment production and sediment delivery can be calculated as a relative dimensionless index. READI predicts only a small proportion of new drains and new surfacing results in the majority of sediment delivery reductions. The model illustrates how the spatial relationships between road and stream networks, controlled by topography and network geometries, influence patterns of road–stream connectivity. READI was applied in seven northern California basins. The model was also applied in a recent burn area to examine how reduced hillslope infiltration can result in increased hydrologic connectivity and sediment delivery.