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.     

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.     

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.     

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.     

ESTIMATING NUTRIENT AND SEDIMENT THRESHOLD CRITERIA FOR BIOLOGICAL IMPAIRMENT IN PENNSYLVANIA WATERSHEDS1

ABSTRACT: This study employs a simple nonlinear statistical approach to establish nitrogen, phosphorus, and sediment concentration and unit area load thresholds to aid in the evaluation of aquatic biological health of watersheds within the state of Pennsylvania. Flow, nitrogen and phosphorus species, sediment, basin area, land cover, and biological assessment data were assembled for 29 Pennsylvania watersheds. For each watershed, rating curves depicting flow versus load relationships were developed using the U.S. Environmental Protection Agency's (USEPA's) storage and retrieval database (STORET) flow and concentration data, then applied to daily flow data obtained from U.S. Geological Survey (USGS) daily flow gauging stations to estimate daily load between 1989 and 1999. The load estimates and concentration data were then sorted into six sets of data: mean annual unit area nitrogen, phosphorus, and sediment loads; and average nitrogen, phosphorus, and sediment concentrations. Results of Mann‐Whitney tests conducted on each of the six datasets indicate that there is a statistically significant difference between the concentrations and unit area loads of nitrogen, phosphorus, and sediment in impaired and unimpaired watersheds. Concentration thresholds, calculated as the midpoint between the impaired and unimpaired watersheds’ 95 percent confidence interval for the median, were estimated to be 2.01 mg/L, 0.07 mg/L, and 197.27 mg/L for nitrogen, phosphorus, and sediment, respectively. Annual unit area load thresholds were estimated to be equal to 8.64 kg/ha, 0.30 kg/ha, and 785.29 kg/ha, respectively, for nitrogen, phosphorus, and sediment species.     

VELOCITY AND CONCENTRATION DISTRIBUTIONS OF SEDIMENT‐LADEN OPEN CHANNEL FLOW1

ABSTRACT: A mathematical model to predict both velocity and concentration distributions for sediment‐laden open channel flow is developed. Velocity profiles are derived by theoretical analysis and numerical method. Logarithmic law and semi‐empirical wake function concept are not adopted. An empirical equation for the ratio of sediment exchange and fluid diffusion coefficients is considered to solve the diffusion equation for suspended‐sediment concentration profiles. Four sets of experimental data from previous researchers are compared to numerical calculation. In the engineering applications, velocity and concentration profiles of sediment‐laden flow can be predicted simultaneously by the present model with the measured velocity and sediment‐concentration at reference level.     

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.     

Integrating channel design and assessment methods based on sediment transport capacity in gravel bed streams

Natural channel design (NCD) and analytical channel design (ACD) are two competing approaches to stable channel design that share fundamental similarities in accounting for sediment transport processes with designs based on hybrid fluvial geomorphology and hydraulic engineering methods. In this paper, we highlight the linkage between ACD's capacity/supply ratio (CSR) and NCD's sediment capacity models (FLOWSED/POWERSED), illustrating how ACD and NCD have reached a point of convergent evolution within the stream restoration toolbox. We modified an existing CSR analytical spreadsheet tool which enabled us to predict relative channel stability using both conventional bed load transport equations and regional sediment regression curves. The stable channel design solutions based on measured data most closely matched the Parker (ACD) and/or Pagosa good/fair (NCD) relationships, which also showed the greatest CSR sensitivity in response to channel alterations. We found that CSR differences among the transport relationships became more extreme the further the design width deviated from the supply reach, suggesting that a stable upstream supply reach may serve as the best design analog. With this paper, we take a step toward resolving lingering controversy in the field of stream restoration, advancing the science and practice by reconciling key differences between ACD and NCD in the context of reach scale morphodynamics.     

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.     

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.     

Long‐Term Trends of Nutrients and Sediment from the Nontidal Chesapeake Watershed: An Assessment of Progress by River and Season

To assess historical loads of nitrogen (N), phosphorus (P), and suspended sediment (SS) from the nontidal Chesapeake Bay watershed (NTCBW), we analyzed decadal seasonal trends of flow‐normalized loads at the fall‐line of nine major rivers that account for >90% of NTCBW flow. Evaluations of loads by season revealed N, P, and SS load magnitudes have been highest in January‐March and lowest in July‐September, but the temporal trends have followed similar decadal‐scale patterns in all seasons, with notable exceptions. Generally, total N (TN) load has dropped since the late 1980s, but particulate nutrients and SS have risen since the mid‐1990s. The majority of these rises were from Susquehanna River and relate to diminished net trapping at the Conowingo Reservoir. Substantial rises in SS were also observed, however, in other rivers. Moreover, the summed rise in particulate P load from other rivers is of similar magnitude as from Susquehanna. Dissolved nutrient loads have dropped in the upland (Piedmont and above) rivers, but risen in two small rivers in the Coastal Plain affected by lagged groundwater input. In addition, analysis of fractional contributions revealed consistent N trends across the upland watersheds. Finally, total N:total P ratios have declined in most rivers, suggesting the potential for changes in nutrient limitation. Overall, this integrated study of historical data highlights the value of maintaining long‐term monitoring at multiple watershed locations.     

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.     

Empirical Estimation of Stream Discharge Using Channel Geometry in Low‐Gradient,Sand‐Bed Streams of the Southeastern Plains

Manning's equation is used widely to predict stream discharge (Q) from hydraulic variables when logistics constrain empirical measurements of in‐bank flow events. Uncertainty in Manning's roughness (n_M) is the major source of error in natural channels, and sand‐bed streams pose difficulties because flow resistance is affected by flow‐dependent bed configuration. Our study was designed to develop and validate models for estimating Q from channel geometry easily derived from cross‐sectional surveys and available GIS data. A database was compiled consisting of 484 Q measurements from 75 sand‐bed streams in Alabama, Georgia, South Carolina, North Carolina (Southeastern Plains), and Florida (Southern Coastal Plain), with six New Zealand streams included to develop statistical models to predict Q from hydraulic variables. Model error characteristics were estimated with leave‐one‐site‐out jackknifing. Independent data of 317 Q measurements from 55 Southeastern Plains streams indicated the model (Q = A_cR_H^0.6906S^0.1216; where A_c is the channel area, R_H is the hydraulic radius, and S is the bed slope) best predicted Q, based on Akaike's information criterion and root mean square error. Models also were developed from smaller Q range subsets to explore if subsets increased predictive ability, but error fit statistics suggested that these were not reasonable alternatives to the above equation. Thus, we recommend the above equation for predicting in‐bank Q of unbraided, sandy streams of the Southeastern Plains.     

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.     

Critical Evaluation of How the Rosgen Classification and Associated “Natural Channel Design” Methods Fail to Integrate and Quantify Fluvial Processes and Channel Response1

Abstract: Over the past 10 years the Rosgen classification system and its associated methods of “natural channel design” have become synonymous to some with the term “stream restoration” and the science of fluvial geomorphology. Since the mid 1990s, this classification approach has become widely adopted by governmental agencies, particularly those funding restoration projects. The purposes of this article are to present a critical review, highlight inconsistencies and identify technical problems of Rosgen’s “natural channel design” approach to stream restoration. This paper’s primary thesis is that alluvial streams are open systems that adjust to altered inputs of energy and materials, and that a form‐based system largely ignores this critical component. Problems with the use of the classification are encountered with identifying bankfull dimensions, particularly in incising channels and with the mixing of bed and bank sediment into a single population. Its use for engineering design and restoration may be flawed by ignoring some processes governed by force and resistance, and the imbalance between sediment supply and transporting power in unstable systems. An example of how C5 channels composed of different bank sediments adjust differently and to different equilibrium morphologies in response to an identical disturbance is shown. This contradicts the fundamental underpinning of “natural channel design” and the “reference‐reach approach.” The Rosgen classification is probably best applied as a communication tool to describe channel form but, in combination with “natural channel design” techniques, are not diagnostic of how to mitigate channel instability or predict equilibrium morphologies. For this, physically based, mechanistic approaches that rely on quantifying the driving and resisting forces that control active processes and ultimate channel morphology are better suited as the physics of erosion, transport, and deposition are the same regardless of the hydro‐physiographic province or stream type because of the uniformity of physical laws.     

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.     

River Channel Geometry and Rating Curve Estimation Using Height above the Nearest Drainage

River channel geometry is an important input to hydraulic and hydrologic models. Traditional approaches to quantify river geometry have involved surveyed river cross sections, which cannot be extended to ungaged basins. In this paper, we describe a method for developing a synthetic rating curve to relate flow to water level in a stream reach based on reach‐averaged channel geometry properties developed using the Height above Nearest Drainage (HAND) method. HAND uses a digital elevation model (DEM) of the terrain and computes the elevation difference between each land surface cell and the stream bed cell to which it drains. Taking increments in water level in the stream, HAND defines the inundation zone and a water depth grid within this zone, and the channel characteristics are defined from this water depth grid. We apply our method to the Blanco River (Texas) and the Tar River (North Carolina) using 10‐m terrain data from the United States Geological Survey (USGS) 3D Elevation Program (3DEP) dataset. We evaluate the method's performance by comparing the reach‐average stage‐river geometry relationships and rating curves to those from calibrated Hydrologic Engineering Center's River Analysis System (HEC‐RAS) models and USGS gage observations. The results demonstrate that after some adjustment, the river geometry information and rating curves derived from HAND using national‐coverage datasets are comparable to those obtained from hydraulic models or gage measurements. We evaluate the inundation extent and show our approach is able to capture the majority of the Federal Emergency Management Agency (FEMA) 100‐year floodplain.     

Sediment Source Fingerprinting: Transforming From a Research Tool to a Management Tool1

Abstract: Information on the nature and relative contribution of different watershed sediment sources is recognized as a key requirement in the design and implementation of targeted management strategies for sediment control. A direct method of assessing sediment sources in a watershed that has attracted attention in recent years is sediment fingerprinting. The aim of this article is to describe the development of sediment fingerprinting as a research tool and to consider how the method might be transformed from a research tool to a management tool within a regulatory framework, with special reference to the United States total maximum daily load (TMDL) program. When compared with the current source assessment tools in developing sediment TMDLs, sediment fingerprinting offers considerable improvement as a tool for quantifying sources of sediment in terms of source type (e.g., channel vs. hillslope) as well as spatial location (subwatershed). While developing a conceptual framework for sediment TMDLs, we recognize sediment fingerprinting along with sediment budgeting and modeling as valuable tools in the TMDL process for developing justifiable sediment TMDLs. The discussions presented in this article may be considered as a first step toward streamlining the sediment fingerprinting approach for its wider application in a regulatory framework.     

Impact of Climate and Land Use Change on Streamflow and Sediment Yield in a Snow‐Dominated Semiarid Mountainous Watershed

This study investigates the impact of climate and land use change on the magnitude and timing of streamflow and sediment yield in a snow‐dominated mountainous watershed in Salt Lake County, Utah using a scenario approach and the Hydrological Simulation Program — FORTRAN model for the 2040s (year 2035–2044) and 2090s (year 2085–2094). The climate scenarios were statistically and dynamically downscaled from global climate models. Land use and land cover (LULC) changes were estimated in two ways — from a regional planning scenario and from a deterministic model. Results indicate the mean daily streamflow in the Jordan River watershed will increase by an amount ranging from 11.2% to 14.5% in the 2040s and from 6.8% to 15.3% in the 2090s. The respective increases in sediment load in the 2040s and 2090s is projected to be 6.7% and 39.7% in the canyons and about 7.4% to 14.2% in the Jordan valley. The historical 50th percentile timing of streamflow and sediment load is projected to be shifted earlier by three to four weeks by mid‐century and four to eight weeks by late‐century. The projected streamflow and sediment load results establish a nonlinear relationship with each other and are highly sensitive to projected climate change. The predicted changes in streamflow and sediment yield will have implications for water supply, flood control and stormwater management.     

A MONTE CARLO TEST OF LOAD CALCULATION METHODS,LAKE TAHOE BASIN,CALIFORNIA‐NEVADA1

ABSTRACT: The sampling of streams and estimation of total loads of nitrogen, phosphorus, and suspended sediment play an important role in efforts to control the eutrophication of Lake Tahoe. We used a Monte Carlo procedure to test the precision and bias of four methods of calculating total constituent loads for nitrate‐nitrogen, soluble reactive phosphorus, particulate phosphorus, total phosphorus, and suspended sediment in one major tributary of the lake. The methods tested were two forms of the Beale's Ratio Estimator, the Period Weighted Sample, and the Rating Curve. Intensive sampling in 1985 (a dry year) and 1986 (a wet year) provided a basis for estimating loads by the “worked record” method for comparison with estimates based on resampling actual data at the lower intensity that characterizes the present monitoring program. The results show that: (1) the Period Weighted Sample method was superior to the other methods for all constituents for 1985; and (2) for total phosphorus, particulate phosphorus, and suspended sediment, the Rating Curve gave the best results in 1986. Modification of the present sampling program and load calculation methods may be necessary to improve the precision and reduce the bias of estimates of total phosphorus loads in basin streams.