ENVIRONMENTAL AND ECONOMIC DAMAGE CAUSED BY SEDIMENT FROM AGRICULTURAL NONPOINT SOURCES1

Environmental and economic damages caused by agricultural nonpoint source inputs of sediment and associated pollutants are examined. Widespread water quality problems are identified in lakes, rivers, and estuaries in agricultural areas, and billions of dollars of on-site and offsite costs result from this eroded soil every year. Some water bodies have been irretrievably damaged, and expensive rehabilitation programs are needed to remedy in-place water pollution problems if Clean Water Act goals are to be achieved. Unless effective abatement and rehabilitation programs are established, billions of dollars of benefits to future generations will not be realized as more waters become irretrievably damaged, and billions more will continue to be spent by government to treat symptoms of these sediment-related problems.     

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.     

AN ALGORITHM FOR ESTIMATING SURFACE SUSPENDED SEDIMENT CONCENTRATIONS WITH LANDSAT MSS DIGITAL DATA1

ABSTRACT: Algorithms for Landsat MSS digital data are needed to reduce the necessity of calibrating each Landsat scene if these data are to be useful in monitoring programs for surface suspended sediments. In this study digital data were extracted from 16 Landsat Multispectral Scanner (MSS) scenes collected between March 1987 and August 1988 over Enid Reservoir in North Central Mississippi. These data were converted to radiance and reflectance data for comparison with field measurements of surface suspended sediment concentrations. Concentrations ranged from 2 to 168 mg/1 during the study with only four greater than 100 mg/l. Linear and polynomial regression analyses were used to relate the surface suspended sediment concentrations with radiance and reflectance. Reflectance in MSS band 2 (0.6 to 0.7 μm) and MSS band 3 (0.7 to 0.8 μm) were best related to the surface suspended sediment concentrations with coefficients of determination accounting for 71 percent and 68 percent of the variation in the data, respectively. Regressions with radiance data accounted for 36 percent (band 2) or less of the variation. Logarithmic transformations of either reflectance or sediment concentrations increase the coefficients of determination for MSS band 2 reflectance data to 81 percent. Regressions between the ratio of MSS band 1 to MSS band 2 reflectances and concentrations also accounted for 80 percent of the variation. An equation Log_e SS (mg/l) = 9.21R½+ 2.71R½2 + 8.45, where S is surface suspended sediment concentrations and R¹/₂ is the ratio of MSS band 1 to MSS band 2 reflectances, provided the best fit to the data with a coefficient of determination of 0.82. This equation is essentially the same as an algorithm proposed by Topliss et at. (1990), for estimating surface suspended sediment concentrations in Canadian coastal waters. These equations for Enid Reservoir and Canadian waters suggest that it may be possible to develop an algorithm for widespread use for estimating surface suspended sediments.     

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.     

EFFECTS OF WATERSHED REPRESENTATION ON RUNOFF AND SEDIMENT YIELD MODELING1

ABSTRACT: This paper evaluates the effects of watershed geometric representation (i.e., plane and channel representation) on runoff and sediment yield simulations in a semiarid rangeland watershed. A process based, spatially distributed runoff erosion model (KINEROS2) was used to explore four spatial representations of a 4.4 ha experimental watershed. The most complex representation included all 96 channel elements identifiable in the field. The least complex representation contained only five channel elements. It was concluded that oversimplified watershed representations greatly influence runoff and sediment yield simulations by inducing excessive infiltration on hillslopes and distorting runoff patterns and sediment fluxes. Runoff and sediment yield decrease systematically with decreasing complexity in watershed representation. However, less complex representations had less impact on runoff and sediment‐yield simulations for small rainfall events. This study concludes that the selection of the appropriate level of watershed representation can have important theoretical and practical implications on runoff and sediment yield modeling in semiarid environments.     

A DYNAMIC SYSTEM MODEL FOR SIMULATING SEDIMENT DISCHARGE1

ABSTRACT: A dynamic sediment discharge model was developed and proposed for the simulation of watershed systems. It war developed from an expansion of splash and flow erosion relationships under steady state conditions. It was described as a general erosion model that can be reduced to forms comparable to many conceptual soil erosion and sediment yield models. The model incorporates eight parameters such as rainfall intensities, runoff rates, and previous sediment discharges. The model was tested with two small watersheds with simulation results which were very satisfactory compared to the data.     

SETTING AN EFFECTIVE TMDL: SEDIMENT LOADING AND EFFECTS OF SUSPENDED SEDIMENT ON FISR1

ABSTRACT: The Agricultural Drainage and Pesticide Transport model was used to examine the relationship between fish and suspended sediment in the context of a proposed total maximum daily load (TMDL) in two agricultural watersheds in Minnesota. During a 50‐year simulation, Wells Creek, a third‐order cold water stream, had an estimated 1,164 events (i.e., one or more consecutive days of estimated sediment loading) and the Chippewa River, a fourth‐order warm water stream, had 906 events of measurable suspended sediment. Sublethal thresholds were exceeded for 970 events and lethal levels for 194 events for brown trout in Wells Creek, whereas adult nonsalmonids would have experienced sublethal levels for 923 events and lethal levels for 241 events. Sublethal levels were exceeded for 756 events and lethal thresholds were exceeded for 150 events in the Chippewa River. Nonsalmonids would have experienced 15 events of mortality between 0 and 20 percent in Wells Creek. In the Chippewa River, there were 35 events of mortality between 0 and 20 percent and one event in which mortality could have exceeded 20 percent. The Minnesota Pollution Control Agency has proposed listing stream reaches as being impaired for turbidity at 25 NTU, which is approximately 46 mg suspended sediment/1. We estimated that 46 mg/1 would be exceeded approximately 30 days in a year (d/yr) in both systems. A TMDL of 46 mg SS/1 may be too high to ensure that stream fishes are not negatively affected by suspended sediment. We recommend that an indicator incorporating the duration of exposure be applied.     

A MULTI-TIER APPROACH FOR AGRICULTURAL WATERSHED MANAGEMENT1

A multi-tier approach for agricultural watershed management has been proposed. The approach involves identification of a watershed management issue/problem, selection or development of simple conceptual model suitable for the exploration of the issue/problem identified and appropriate to the database available, and application of the model the address the identified issue/problem. The procedure is repeated by increasing the complexity in the conceptual model until the identified issue/problem has been addressed satisfactorily. An application of the procedure to an example watershed in southern Ontario conditions is shown. The application example has revealed that for identification of temporal pattern of runoff and sediment loads a simple conceptual model is adequate. For identification of spatial location of the sediment source areas and for the development of a monitoring program for the evaluation of remedial strategies a more complex distributed agricultural watershed model is necessary.     

SPNM,A MODEL FOR PREDICTING SEDIMENT,PHOSPHORUS, AND NITROGEN YIELDS FROM AGRICULTURAL BASINS1

ABSTRACT: A model called SPNM from the words “sediment-phosphorus-nitrogen model” was developed for simulating agricultural contributions to water pollution. SPNM is designed to predict sediment, P, and N yields for individual storms on small basins and to route these yields through streams and valleys of large basins. Users need no computer programming experience because the model is a problem-oriented computer language. SPNM is useful in planning water resources projects and in research. Tests of the model on a watershed provided realistic results.     

SEDIMENT RESUSPENSION AND DRAWDOWN IN A WATER SUPPLY RESERVOIR1

ABSTRACT: The magnitudes and patterns of sediment resuspension are assessed in Cannonsville Reservoir, New York, to quantify and characterize this internal source of sediment. The assessment is based on analyses of sediment trap collections from 10 sites over the spring to fall interval of two years. Temporal and spatial patterns in sediment deposition are demonstrated to be driven by resuspension/redeposition processes. Sediment that had been resuspended and redeposited represented 80 to 96 percent, on average, of the depositing solids collected along the main axis of the lake. About 90 percent of the redeposited sediment was inorganic. Increased resuspension caused by drawdown of the reservoir surface and fall turnover resulted in 10 to 50‐fold increases in deposition rates compared to levels observed when the reservoir was full and strongly thermally stratified. Elevated levels of redeposition from resuspension in the reservoir have been driven by both higher water column concentrations of suspended solids and settling velocities. Recurring longitudinal and lateral gradients in resuspension are delineated, establishing that resuspended solids are transported from the riverine to the lacustrine zone and from near‐shore to pelagic areas. Resuspension is demonstrated to cause increases in inanimate particle (tripton) concentrations. Higher tripton levels have been observed in years with greater drawdown. Water quality impacts of the resuspension phenomenon are considered.     

THE USE OF MULTIPLE REGRESSION MODELS IN PREDICTING SEDIMENT YIELD1

ABSTRACT. Four commonly used models for predicting sediment yield are analyzed and compared using previously published data. Three of these models involve logarithmic transformations. Some of the problems involved in transforming data are discussed in the context of logarithmic transformations. These problems are illustrated using the results of standard regression analyses and economic loss function analyses. For the data analyzed, the linear model is preferable to each of the logarithmic models on the basis of each analysis, and the usual multiple objective nature of the model choice problem is thus modified. The extent to which these results can be generalized is discussed in the context of model choice.     

BANK STABILITY ANALYSIS FOR PREDICTING REACH SCALE LAND LOSS AND SEDIMENT YIELD1

ABSTRACT: When extensive lengths of river become stabilized (e.g., by channel incision), riverbank erosion can result in considerable riparian land loss and the delivery of large volumes of sediment downstream. The ability to predict the stability and failure geometry of eroding riverbanks is therefore an important prerequisite in estimating the rate of bank erosion and sediment yield associated with bank erosion. In this paper, a new stability analysis for layered river banks is introduced. The new analysis differs from many previous analyses in that it takes into consideration the effects of positive pore water pressure in the saturated portion, and negative pore water pressure in the unsaturated portion, of the bank as well as the influence of hydrostatic confining pressure due to the water level in the river. In addition, the failure plane is not constrained to pass through the toe of the bank and the bank profile geometry is not restricted to an idealized special case. The predictive ability of the new bank stability analysis is assessed using data from two field sites. Subsequently, a methodology for applying the bank stability analysis at the scale of the river reach (0.1 to 10 km) is discussed. This method involves the use of empirical models of bed level adjustment to estimate the magnitude of incision at specific locations along the reach, with these estimates used to drive the stability analysis. Application of the new method is demonstrated with an example.     

OPTIMIZATION PROCEDURE FOR COST EFFECTIVE BMP PLACEMENT AT A WATERSHED SCALE1

ABSTRACT: A combinatorial optimization procedure for best management practice (BMP) placement at the watershed level facilitates selection of cost effective BMP scenarios to control non point source (NFS) pollution. A genetic algorithm (GA) was selected from among several optimization heuristics. The GA combines an optimization component written in the C++ language with spatially variable NFS pollution prediction and economic analysis components written within the Arc View geographic information system. The procedure is modular in design, allowing for component modifications while maintaining the basic conceptual framework. An objective function was developed to lexicographically optimize pollution reduction followed by cost increase. Scenario cost effectiveness is then calculated for scenario comparisons. The NPS pollutant fitness score allows for evaluation of multiple pollutants, based on prioritization of each pollutant. The economic component considers farm level public and private costs, cost distribution, and land area requirements. Development of a sediment transport function, used with the Universal Soil Loss Equation, allows the optimization procedure to run within a reasonable timeframe. The procedure identifies multiple near optimal solutions, providing an indication of which fields have a more critical impact on overall cost effectiveness and flexibility in the final solution selected for implementation. The procedure was demonstrated for a 1,014‐ha watershed in the Ridge and Valley physiographic region of Virginia.     

A COMPARISON OF SURFACE RUNOFF AND SEDIMENT YIELDS FROM LOW- AND HIGH-SEVERITY SITE PREPARATION BURNS1

ABSTRACI:. Slash burning is a common site preparation technique used after timber harvest throughout the Southeastern United States. Little quantitative information exists on the hydrologic response to burn severity. This study compared the effects of low-severity and high-severity burns on runoff and sediment yields during rainfall simulation and during natural rainfall in the Southern Appalachian Mountains. Fire severity was largely determined by moisture conditions of the forest floor prior to ignition. Runoff and sediment yield variability was high between plots within the same treatment area due to differences in forest floor characteristics and infiltration rates. Conditions of high-severity resulted when burning was conducted with relatively dry fuels. Sediment yields were 40-times greater for the high-severity treatment areas than the low-severity treatment areas.     

MODELING METALS TRANSPORT AND SEDIMENT/WATER INTERACTIONS IN A MINING IMPACTED MOUNTAIN STREAM1

ABSTRACT: The U.S. Environmental Protection Agency (USEPA) Water Quality Analysis Simulation Program (WASP5) was used to model the transport and sediment/water interactions of metals under low flow, steady state conditions in Tenmile Creek, a mountain stream supplying drinking water to the City of Helena, Montana, impacted by numerous abandoned hard rock mines. The model was calibrated for base flow using data collected by USEPA and validated using data from the U.S. Geological Survey (USGS) for higher flows. It was used to assess metals loadings and losses, exceedances of Montana State water quality standards, metals interactions in stream water and bed sediment, uncertainty in fate and transport processes and model parameters, and effectiveness of remedial alternatives that include leaving contaminated sediment in the stream. Results indicated that during base flow, adits and point sources contribute significant metals loadings to the stream, but that shallow ground water and bed sediment also contribute metals in some key locations. Losses from the water column occur in some areas, primarily due to adsorption and precipitation onto bed sediments. Some uncertainty exists in the metal partition coefficients associated with sediment, significance of precipitation reactions, and in the specific locations of unidentified sources and losses of metals. Standards exceedances are widespread throughout the stream, but the model showed that remediation of point sources and mine waste near water courses can help improve water quality. Model results also indicate, however, that alteration of the water supply scheme and increasing base flow will probably be required to meet all water quality standards.     

A DECISION MODEL TO PREDICT SEDIMENT YIELD FROM FOREST PRACTICES1

ABSTRACT: In order to choose among “best management practices,” forest managers need to predict sediment yield to perennial streams following various forest land operations. The “universal soil loss equation” (USLE) is not directly applicable to forest operations because of the heterogenous soil surface conditions left by harvesting, site preparation, and planting. A sediment hazard index is proposed, to be based on the amount of exposed mineral soil and its proximity to streams. The model offered includes rainfall erosivity, soil erodibility and average land slope, together with the index W. A paired watershed experiment in the central Georgia Piedmont was used to estimate parameters in the model. The experimental basin (80 acres) was clearcut, drum roller chopped twice, and planted by machine. The standard error of estimate of sediment yield was computed to be about 50/lbs/ac per sampling period (four months). Use of William's erogivity index (storm flow times peak flow) reduced the standard error to 33/lbs/ac.     

DEVELOPMENT OF A GIS‐BASED FLOOD INFORMATION SYSTEM FOR FLOODPLAIN MODELING AND DAMAGE CALCULATION1

ABSTRACT: This work presents a flexible system called GIS‐based Flood Information System (GFIS) for floodplain modeling, flood damages calculation, and flood information support. It includes two major components, namely floodplain modeling and custom designed modules. Model parameters and input data are gathered, reviewed, and compiled using custom designed modules. Through these modules, it is possible for GFIS to control the process of flood‐plain modeling, presentation of simulation results, and calculation of flood damages. Empirical stage‐damage curves are used to calculate the flood damages. These curves were generated from stage‐damage surveys of anthropogenic structures, crops, etc., in the coastal region of a frequently flooded area in Chia‐I County, Taiwan. The average annual flood damages are calculated with exceedance probability and flood damages for the designed rainfalls of 2, 5, 10, 25, 50, 100, and 200 year recurrence intervals with a duration of 24 hours. The average annual flood depth in this study area can also be calculated using the same method. The primary advantages of GFIS are its ability to accurately predict the locations of flood area, depth, and duration; calculate flood damages in the floodplain; and compare the reduction of flood damages for flood mitigation plans.     

AN EFFICIENT RESPONSE MATRIX METHOD FOR COUPLING A GROUNDWATER SIMULATOR AND A REGIONAL AGRICULTURAL MANAGEMENT MODEL1

ABSTRACT: By extending the concept of response matrix to consider “active” and “passive” effects, an efficient response matrix method is developed for coupling a groundwater simulator and a regional agricultural management model The method eliminates the need to store all of the recovery information from preceding time periods. Active effects are those which occur during the actual application of a pumping or recharge stress while passive effects represent the recovery of water levels from an initial departure from steady-state conditions at the beginning of a time step. Derivation of the required matrices and a numerical example are presented for the Salinas Valley groundwater basin in California.     

A COMPARISON OF SINGLE‐CELL AND MULTICELL CULVERTS FOR STREAM CROSSINGS1

ABSTRACT: Single‐barrel culverts are a common means of roadway crossings for smaller streams. While this culvert design provides an economical solution for a crossing, the adverse effects of conveying the stream through a single opening can be far reaching. The single‐barrel culvert is typically sized for a design storm much greater than the channel forming discharge. This oversizing causes an interruption of the normal flow patterns and sediment transport for the system. Shallow depths at low flow in the pipe and perching at the outlet can impede fish passage. Multicell culverts (where the main culvert at the channel invert is sized for bankfull discharge, and additional pipes are placed at the floodplain elevation to convey overbank flow up to the design discharge) have been recommended as a best management practice to minimize erosion and improve fish passage. This flume study scaled a prototype single‐barrel culvert to both a single‐cell model, and a multicell design to compare outlet scour and flow depths within the culvert. The results provide designers and planners with evidence of the benefits of multicell culverts to justify the higher cost of installation compared to single‐barrel culverts.