APPLICATION OF A SEDIMENT ROUTING MODEL TO A MIDDLE EUROPEAN WATERSHED1

ABSTRACT: A simulation model that computes sediment yield due to sheet and rill erosion at the outlet of a large watershed requires daily precipitation and the soil, topographic, and vegetative characteristics of the watershed. An important problem, particularly in a large watershed, is the transport of sediment produced in the sub-watersheds to the outlet of the whole watershed. This problem is approached mathematically by a sediment routing model that is used as a component of the total model.     

MODELING UPLAND EROSION1

ABSTRACT: We have developed a computer model of soil loss on an upland watershed from the continuity considerations for sediment transport and from equations describing rill and interrill erosion. The model is based on dividing the upland area into a grid containing rill and interill zones, on the Universal Soil Loss Equation (USLE), and on equations describing detachment and transport capacity of rill flow. The USLE estimates the sediment load from the contributing areas. The location and amount of total erosion and deposition are determined by comparing the transport and detachment capacity of rill flow for specific storms. The model considers the mechanics of erosion process and can serve as basis for reservoir and channel design and land use planning.     

Modeling Pre‐ and Post‐Dam Removal Sediment Dynamics: The Kalamazoo River,Michigan1

Abstract:  The state of Michigan is interested in removing two low‐head dams in an 8.8 km reach of the Kalamazoo River between Plainwell and Otsego, Michigan, while minimizing impacts locally and to downstream reaches. The study was designed to evaluate the erosion, transport, and deposition of sediments over a 37.3‐year period using the channel evolution model CONCEPTS for three simulation scenarios: Dams In (DI), Dams Out (DO), and Design (D). The total mass of sediment emanating from the channel boundary, for the DI case, shows net deposition of 4,100 T/y for the study reach, with net transport (suspended and bed load) of 10,500 T/y passing the downstream boundary. For the DO case, net erosion is 19,200 T/y with net transport of 30,100 T/y (187% increase) passing the downstream boundary. For the D case, net deposition is 2,570 T/y (37% decrease) with transport of 14,200 T/y (35% increase) passing the downstream boundary. The most significant findings were: (1) removal of the low‐head dams will cause significant erosion of sediments stored behind the dams and increased sediment loads passing the downstream boundary and (2) sediment loads for the proposed channel design are similar to existing conditions and offer reduced fine‐sediment loadings.     

Identification of Sediment Sources in Forested Watersheds With Surface Coal Mining Disturbance Using Carbon and Nitrogen Isotopes1

Abstract: Sediments and soils were analyzed using stable carbon and nitrogen isotope ratio mass spectrometry and carbon and nitrogen elemental analyses to evaluate the their ability to indicate land‐use and land management disturbance and pinpoint loading from sediment transport sources in forested watersheds disturbed by surface coal mining. Samples of transported sediment particulate organic matter were collected from four watersheds in the Southern Appalachian forest region of southeastern Kentucky. The four watersheds had different surface coal mining history that were classified as undisturbed, active mining, and reclaimed conditions. Soil samples were analyzed including reclaimed grassland soils, undisturbed forest soils, geogenic organic matter associated with coal fragments in mining spoil, and soil organic matter from un‐mined grassland soils. Statistically significant differences were found for all biogeochemical signatures when comparing transported sediments from undisturbed watersheds and surface coal mining disturbed watersheds, and the results were attributed to differences in erosion sources and the presence of geogenic organic matter. Sediment transport sources in the surface coal mining watersheds were analyzed using Monte Carlo mass balance un‐mixing and it was found that: δ¹⁵N showed the ability to differentiate streambank erosion and surface soil erosion; and δ¹³C showed the ability to differentiate soil organic matter and geogenic organic matter. Results from the analyses suggest that streambank erosion downstream of surface coal mining sites is an especially significant source of sediment in coal mining disturbed watersheds. Further, the results suggest that the sediment transport processes governing streambank erosion loads are taking longer to reach geomorphologic equilibrium in the watershed as compared with the surface erosion processes. The dual‐isotope technique provides a useful method for further investigation of the impact of surface coal mining in the uplands of the watershed upon the geomorphologic state of the channel and the source of organic matter in aquatic systems impacted by surface coal mining.     

THE TWO‐DIMENSIONAL UPLAND EROSION MODEL CASC2D‐SED1

ABSTRACT: The two‐dimensional soil erosion model CASC2D‐SED simulates the dynamics of upland erosion during single rainstorms. The model is based on the raster‐based surface runoff calculations from CASC2D. Rainfall precipitation is distributed in time and space. Infiltration is calculated from the Green‐Ampt equations. Surface runoff is calculated from the diffusive wave approximation to the Saint‐Venant equations in two‐dimensions. Watershed data bases in raster Geographical Information System (GIS) provide information on the soil type, size fractions, soil erodibility, cropping management, and conservation practice factors for soil erosion calculations. Upland sediment transport is calculated for the size fractions (sand, silt, and clay), and the model displays the sediment flux, the amount of suspended sediment, and the net erosion and deposition using color graphics. The model has been tested on Goodwin Creek, Mississippi. The peak discharge and time to peak are within ± 20 percent and sediment transport rates within ?50 percent to 200 percent.     

Impacts of Land‐Cover Change on Suspended Sediment Transport in Two Agricultural Watersheds1

Schilling, Keith E., Thomas M. Isenhart, Jason A. Palmer, Calvin F. Wolter, and Jean Spooner, 2011. Impacts of Land‐Cover Change on Suspended Sediment Transport in Two Agricultural Watersheds. Journal of the American Water Resources Association (JAWRA) 47(4):672‐686. DOI: 10.1111/j.1752‐1688.2011.00533.x Abstract: Suspended sediment is a major water quality problem, yet few monitoring studies have been of sufficient scale and duration to assess the effectiveness of land‐use change or conservation practice implementation at a watershed scale. Daily discharge and suspended sediment export from two 5,000‐ha watersheds in central Iowa were monitored over a 10‐year period (water years 1996‐2005). In Walnut Creek watershed, a large portion of land was converted from row crop to native prairie, whereas in Squaw Creek land use remained predominantly row crop agriculture. Suspended sediment loads were similar in both watersheds, exhibiting flashy behavior typical of incised channels. Modeling suggested that expected total soil erosion in Walnut Creek should have been reduced 46% relative to Squaw Creek due to changes in land use, yet measured suspended sediment loads showed no significant differences. Stream mapping indicated that Walnut Creek had three times more eroding streambank lengths than did Squaw Creek suggesting that streambank erosion dominated sediment sources in Walnut Creek and sheet and rill sources dominated sediment sources in Squaw Creek. Our results demonstrate that an accounting of all sources of sediment erosion and delivery is needed to characterize sediment reductions in watershed projects combined with long‐term, intensive monitoring and modeling to account for possible lag times in the manifestation of the benefits of conservation practices on water quality.     

A Sediment Budget for an Urbanizing Watershed, 1951‐1996, Montgomery County, Maryland,U.S.A.1

Abstract: Despite widespread interest, few sediment budgets are available to document patterns of erosion and sedimentation in developing watersheds. We assess the sediment budget for the Good Hope Tributary, a small watershed (4.05 km²) in Montgomery County, Maryland, from 1951‐1996. Lacking monitoring data spanning the period of interest, we rely on a variety of indirect and stratigraphic methods. Using regression equations relating sediment yield to construction, we estimated an upland sediment production of 5,700 m³ between 1951 and 1996. Regression equations indicate that channel cross‐sectional area is correlated with the extent of development; these relationships, when combined with historical land use data, suggest that upland sediment yield was augmented by 6,400 m³ produced by enlargement of first‐order and second‐order stream channels. We used dendrochronology to estimate that 4,000 m³ of sediment was stored on the floodplain from 1951‐1996. The sediment yield from the watershed, obtained by summing upstream contributions, totals 8,100 m³ of sediment, or 135 tons/km²/year. These results indicate that upland erosion, channel enlargement, and floodplain storage are all significant components of the sediment budget of our study area, and all three are approximately equal in magnitude. Erosion of “legacy” floodplain sediments originally deposited during poor agricultural practices of the 19th and early 20th Centuries has likely contributed between 0 and 20% of the total sediment yield, indicating that these remobilized deposits are not a dominant component of the sediment yield of our study area.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

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.     

Quantifying Reductions of Mass‐Failure Frequency and Sediment Loadings From Streambanks Using Toe Protection and Other Means: Lake Tahoe,United States1

Abstract: Streambank erosion by mass‐failure processes represents an important form of channel adjustment and a significant source of sediment in disturbed streams. Mass failures regularly occur by a combination of hydraulic processes that undercut bank toes and geotechnical processes that cause bank collapse by gravity. Little if any quantitative information is available on the effectiveness of bank treatments on reducing erosion. To evaluate potential reduction in sediment loadings emanating from streambanks, the hydraulic and geotechnical processes responsible for mass failure were simulated under existing and mitigated conditions using a Bank‐Stability and Toe‐Erosion Model (BSTEM). Two critical erosion sites were selected from each of the three watersheds known to contribute the greatest amounts of fine sediment by streambank processes in the Lake Tahoe Basin. A typical high‐flow annual hydrograph was selected for analysis. Bank‐material strength data were collected for each layer as were species‐specific root‐reinforcement values. The effects of the first flow event on bank‐toe erosion were simulated using an excess shear‐stress approach. The resulting geometry was then exported into the bank‐stability submodel to test for the relative stability of the bank under peak flow and drawdown conditions. In this way, BSTEM was used iteratively for all flow events for both existing and mitigated conditions. On average, 13.6% of the material was eroded by hydraulic shear, the remainder by mass failures, which occurred about five times over the simulation period. Simulations with 1.0 m‐high rock‐toe protection showed a dramatic reduction in streambank erosion (69‐100%). Failure frequency for the simulation period was reduced in most cases to a single episode. Thus, an almost 90% reduction in streambank loadings was achieved by virtually eliminating the erosion of only 14% of the material that was entrained by hydraulic forces. Consequently, simulations show average load reductions of about an order of magnitude. Results stress the critical importance of protecting the bank toe‐region from steepening by hydraulic forces that would otherwise entrain previously failed and in situ bank materials, thereby allowing the upper bank to flatten (by failure) to a stable slope.     

Contribution of In‐Channel Processes to Sediment Yield of an Urbanizing Watershed1

Abstract: A study was conducted between September 2003 and September 2006 to obtain baseline sediment inventories and monitor sediment transport and storage along a 3.7 km length of the channel of Valley Creek within Valley Forge National Historical Park, Pennsylvania. Valley Creek is a tributary of the Schuylkill River and drains an urbanizing 60.6 km² watershed that currently has 18% impervious land cover. Numerous field methods were employed to measure the suspended sediment yield, longitudinal profile, cross‐sections, banklines, and particle size distribution of the streambed. Suspended sediment yield for the watershed was measured at a USGS gage located just upstream of the park boundary between July 2004 and July 2005, the period corresponding to field surveys of bank erosion and channel change. The estimated suspended sediment yield of 95.7 t/km²/year is representative of a year with unusually high discharge, including a storm event that produced a peak of 78 m³/s, the second highest discharge on record for the USGS gage. Based on the median annual streamflow for the 24 years of record at the USGS gage from 1983 to 2006, the median annual sediment yield is estimated to be closer to 34 t/km²/year, considerably lower than median and mean values for other sites within the region. The mass of silt, clay, and fine sand derived from bank erosion along the 3.7 km study reach during the field survey period accounts for an estimated 2,340 t, equivalent to about 43% of the suspended sediment load. The mass of fine sediment stored in the bed along the study reach was estimated at 1,500 t, with about 330 t of net erosion during the study period. Although bank erosion appears to be a potentially dominant source of sediment by comparison with annual suspended sediment load, bed sediment storage and potential for remobilization is of the same order of magnitude as the mass of sediment derived from bank erosion.     

A GIS‐ASSISTED DISTRIBUTED WATERSHED MODEL FOR SIMULATING FLOODING AND INUNDATION1

ABSTRACT: A distributed watershed model combining kinematic wave routing, 1‐D dynamic channel‐flow routing, and 2‐D diffusive overland‐flow routing has been developed to simulate flooding and inundation levels of large watersheds. The study watershed was linked to a GIS database and was divided into an upstream mountainous area and a downstream alluvial plain. A kinematic wave routing was adopted at the mountainous area to compute the discharge flowing into the alluvial plain. A 1‐D dynamic channel routing solving the St. Venant equations by the Preissmann method was performed for the main channel of the alluvial plain, whereas a 2‐D overland‐flow routing solving the diffusion wave equation with the Alternating Direction Explicit scheme was used for floodplains. The above two routings were connected by weir‐link discharge formula. The parameters in the model were calibrated and independently verified by single‐event storms. An example application of flooding/inundation analysis was conducted for the Taichung station and the Woozi depot (Taiwan High Speed Rail). Suggested inundation‐proofing measures ‐ including raising ground surface elevation of the station and depot and building a waterproofing exterior wall and their combination ‐ were investigated. It was concluded that building the waterproofing exterior wall had a strong tendency to decrease peak inundation depth.     

Using Fly Ash as a Marker to Quantify Culturally‐Accelerated Sediment Accumulation in Playa Wetlands

Wetlands in the Rainwater Basin in Nebraska are vulnerable to sediment accumulation from the surrounding watershed. Sediment accumulation has a negative impact on wetland quality by decreasing the depth and volume of water stored, and the plant community species composition and density growing in the wetland. The objective of this study was to determine the amount of sediment that has accumulated in five selected wetlands in the Rainwater Basin in Nebraska. Soil cores were taken at five or six locations along transects across each wetland. This study used the fly ash, which is generated by coal‐burning locomotives that were present generally in the late 1800s and early 1900s, as a marker to quantify the sediment deposition rates. The cores were divided into 5 cm sections and the soils were analyzed using a fly ash extraction and identification technique. Results indicate that the average depth of sediment ranged from 23.00 to 38.00 cm. The annual average depth of sediment accumulation ranged from 0.18 cm/yr to 0.29 cm/yr. The annual sediment accumulation rate from both wind erosion and water erosion in these five sampling wetlands was between 1.946 and 3.225 kg/m²/yr. The results of this research can be used to develop restoration plans for wetlands. The fly ash testing technology can also be applied to other areas with the railroads across the United States.     

Implications of Conceptual Channel Representation on SWAT Streamflow and Sediment Modeling

Hydrologic modeling outputs are influenced by how a watershed system is represented. Channel routing is a typical example of the mathematical conceptualization of watershed landscape and processes in hydrologic modeling. We investigated the sensitivity of accuracy, equifinality, and uncertainty of Soil and Water Assessment Tool (SWAT) modeling to channel dimensions to demonstrate how a conceptual representation of a watershed system affects streamflow and sediment modeling. Results showed the amount of uncertainty and equifinality strongly responded to channel dimensions. On the other hand, the model performance did not significantly vary with the changes in the channel representation due to the degree of freedom allowed by the conceptual nature of hydrologic modeling in the parameter calibration. Such findings demonstrated good modeling performance statistics do not necessarily mean small output uncertainty, and partial improvements in the watershed representation may neither increase modeling accuracy nor reduce uncertainty. We also showed the equifinality and uncertainty of hydrologic modeling are case‐dependent rather than specific to models or regions, suggesting great caution should be used when attempting to transfer uncertainty analysis results to other modeling studies, especially for ungauged watersheds. Editor's note: This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Sediment Availability and Erosion Rates on Off‐Highway Vehicle Trails in the Ouachita Mountains,USA

Factors influencing sediment availability are assessed and erosion rates are quantified for an off‐highway vehicle (OHV) trail system in the Ouachita Mountains of Arkansas. As of May 2012, the Wolf Pen Gap trail system included 77.0 km of "trails" which consist of county roads; open and closed Forest Service roads; and open and closed OHV trails. For a given trail length, the sediment volume available to be eroded is determined by bare trail width and sediment depth. Four condition types are defined that group trail sections based on statistically different trail widths or depths. Trail construction method appears to influence sediment availability differences more than erosion potential (as indexed by trail slope gradient and length). The range for annual trail erosion rates is estimated as 75 and 210 tonne/ha/yr. The high and low rates are obtained using two independent methods. The 210 tonne/ha/yr rate is computed from mean sediment capture at 30 sediment traps installed for 0.5–1.0 year. The 75 tonne/ha/yr rate is computed assuming all available trail sediment measured in a one‐time sampling is eroded over the next year. We argue in support of this assumption and suggest both rate values may be conservative. Trail erosion rates and sediment trap observations indicate frequent trap cleanout will be needed to continue sediment capture from All Terrain Vehicle trails.     

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.     

Evaluation of SWAT Impoundment Modeling Methods in Water and Sediment Simulations

Worldwide studies show 80%–90% of all sediments eroded from watersheds is trapped within river networks such as reservoirs, ponds, and wetlands. To represent the impact of impoundments on sediment routing in watershed modeling, Soil and Water Assessment Tool (SWAT) developers recommend to model reservoirs, ponds, and wetlands using impoundment tools (ITs). This study evaluates performance of SWAT ITs in the modeling of a small, agricultural watershed dominated by lakes and wetlands. The study demonstrates how to incorporate impoundments into the SWAT model, and discusses and evaluates involved parameters. The study then recommends an appropriate calibration sequence, i.e., landscape parameters calibration, followed by pond/wetlands calibration, then channel parameter calibrations, and lastly, reservoir parameter calibration. Results of this study demonstrate not following SWAT recommendation regarding modeling water land use as an impoundment depreciates SWAT performance, and may lead to misplaced calibration efforts and model over‐calibration. Further, the chosen method to model impoundments’ outflow significantly impacts sediment loads in the watershed, while streamflow simulation is not very sensitive. This study also allowed calculation of mass accumulation rates in modeled impoundments where the annual mass accumulation rate in wetlands (2.3 T/ha/yr) was 39% higher than mass accumulation rate in reservoirs (1.4 T/ha/yr).     

Sediment Source Assessment in a Lowland Watershed Using Nitrogen Stable Isotopes1

Fox, James F., Charles M. Davis, and Darren K. Martin, 2010. Sediment Source Assessment in a Lowland Watershed Using Nitrogen Stable Isotopes. Journal of the American Water Resources Association (JAWRA) 46(6):1192–1204. DOI: 10.1111/j.1752-1688.2010.00485.x Abstract: Sediment sources and transported sediments were sampled in a lowland watershed with pronounced fine sediment storage in the streambed. Sediments were analyzed for carbon and nitrogen content and stable nitrogen isotopic composition. Analysis of the data shows that temporarily stored streambed sediments dominate the sediment load during moderate- and low-flow hydrologic events. Modeling of sediment transport and nitrogen elemental and isotopic mass balance was performed for the watershed for a 12-month time period using a continuous, conceptual-based model. The model results show that during moderate- and low-flow hydrologic events, the streambed is slowly downcutting. During very high-flow hydrologic events, deposition is pronounced in the streambed and sediment is replenished to the bed. Nitrogen model results show that elemental and isotopic nitrogen of streambed sediments vary substantially over the simulation period. In this manner, the streambed in a lowland watershed functions as a temporary storage zone that, in turn, can impact the nitrogen elemental and isotopic signature of sediments. The variation could significantly impact estimates of sediment provenance using nitrogen tracer-based methods. Future work should consider both hydrologic and biogeochemical control on the nitrogen isotopic signature of sediments in small lowland watersheds and streams where a significant portion of deposited fines are temporarily stored.     

Monitored and Modeled Correlations of Sediment and Nutrients with Chesapeake Bay Water Clarity

This article analyzes the correlations of the observed and modeled light attenuation coefficient, Kd, with in situ total suspended solids (TSS) and chlorophyll‐a concentrations in Chesapeake Bay (CB) tidal waters, and with sediment and nutrient loads from the Chesapeake watershed. Light attenuation is closely related to in situ TSS and chlorophyll‐a concentrations, however, the strength of the correlation differs among the CB segments. There are distinct differences between saline and tidal fresh segments in the main Bay, but less distinction among saline and tidal fresh segments in the tidal tributaries. The correlation between Kd with sediment and nutrient loads is complicated by the lag times of TSS and the chlorophyll‐a responses to reductions in nutrient and sediment loads from the watershed, and also due to the diverse load sources. Three sets of model sensitivity scenarios were performed with: (1) differential sediment and nutrient loads; (2) selective sediment source types; and (3) geographically isolated inputs. The model results yield similar findings as those based on observed data and provide information regarding the effect of sediment on specific water bodies. Based on the model results a method was developed to determine sediment and nutrient load reductions needed to achieve the water clarity standards of the CB segments.