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).     

“R1-R4” AND “BOISED” SEDIMENT PREDICTION MODEL TESTS USING FOREST ROADS IN GRANITICS1

ABSTRACT: Erosion and sedimentation data from research watersheds in the Silver Creek Study Area in central Idaho were used to test the prediction of logging road erosion using the R1-R4 sediment yield model, and sediment delivery using the “BOISED” sediment yield prediction model. Three small watersheds were instrumented and monitored such that erosion from newly constructed roads and sediment delivery to the mouths of the watersheds could be measured for four years following road construction. The errors for annual surface erosion predictions for the two standard road tests ranged from +31.2 t/ha/yr (+15 percent) to -30.3 t/ha/yr (-63 percent) with an average of zero t/ha/yr and a standard deviation of the differences of 18.7 t/ha/yr. The annual prediction errors for the three watershed scale tests had a greater range from -40.8 t/ha/yr (-70 percent) to +65.3 t/ha/yr (+38 percent) with a mean of -1.9 t/ha/yr and a standard deviation of the differences of 25.2 t/ha/yr. Sediment yields predicted by BOISED (watershed scale tests) were consistently greater (average of 2.5 times) than measured sediment yields. Hillslope sediment delivery coefficients in BOISED appear to be overly conservative to account for average site conditions and road locations, and thus over-predict sediment delivery. Mass erosion predictions from BOISED appear to predict volume well (465 tonnes actual versus 710 tonnes predicted, or a 35 percent difference) over 15 to 20 years, however mass wasting is more episodic than the model predicts.     

HYDROLOGIC PARAMETERIZATION OF WATERSHEDS FOR RUNOFF PREDICTION USING SWAT1

ABSTRACT: The use of continuous time, distributed parameter hydrologic models like SWAT (Soil and Water Assessment Tool) has opened several opportunities to improve watershed modeling accuracy. However, it has also placed a heavy burden on users with respect to the amount of work involved in parameterizing the watershed in general and in adequately representing the spatial variability of the watershed in particular. Recent developments in Geographical Information Systems (GIS) have alleviated some of the difficulties associated with managing spatial data. However, the user must still choose among various parameterization approaches that are available within the model. This paper describes the important parameterization issues involved when modeling watershed hydrology for runoff prediction using SWAT with emphasis on how to improve model performance without resorting to tedious and arbitrary parameter by parameter calibration. Synthetic and actual watersheds in Indiana and Mississippi were used to illustrate the sensitivity of runoff prediction to spatial variability, watershed decomposition, and spatial and temporal adjustment of curve numbers and return flow contribution. SWAT was also used to predict stream runoff from actual watersheds in Indiana that have extensive subsurface drainage. The results of this study provide useful information for improving SWAT performance in terms of stream runoff prediction in a manner that is particularly useful for modeling ungaged watersheds wherein observed data for calibration is not available.     

Riparian buffers and potentially unstable ground

The spatial coincidence between riparian buffers of various widths and extents and potentially unstable ground was quantified using a physically based model for shallow landslide initiation and GIS for two watersheds on the Olympic Peninsula, Washington, USA. The proportion of the potentially unstable ground in each watershed within riparian buffers is a function of both buffer width and the extent of the stream channel network being buffered. While current buffers required by Washington State cover less than 5% of the potentially unstable ground, buffering all stream channels in these watersheds with 100-m buffers covered 75%–90% of the potentially unstable areas. Our analyses further show that: (1) riparian buffers are not efficient mechanisms for protecting potentially unstable ground, and (2) identifying potentially unstable ground using a physically based model should prove more effective for designing methods to reduce shallow landsliding hazards than relying on extensive buffer zones along stream channels.     

REELFOOT LAKE SEDIMENTATION RATES AND SOURCES1

ABSTRACT: The three basins of Reelfoot Lake, which is located in northwestern Tennessee, were investigated using the Cs-137 tracer technique to determine rates of sediment deposition and to estimate the time before the basins will fill with sediment. Blue Basin, the largest of the three basins with 2922 ha, had an average annual sedimentation rate of 0.9 cm/yr from 1984 to 1984. The basin will become too shallow for most boating and recreational activities in about 200 years. Buck Basin, the central basin with 774 ha, had an average annual sedimentation rate of 1.1 cm/yr and will become too shallow for most recreational uses in about 100 years. Upper Blue Basin, the most upstream and smallest basin with 439 ha, had an average annual sedimentation rate of 1.7 cm/yr and will become too shallow for most recreational uses in about 60 years. Two important sources of sediment to Reelfoot Lake are erosion from a large number of soybean fields and channelization of many of the streams that flow into the lake. Changes in land management that would reduce erosion could increase the time the lake would remain usable for recreational activities.     

Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds1

Abstract: The watershed scale Soil and Water Assessment Tool (SWAT) model divides watersheds into smaller subwatersheds for simulation of rainfall‐runoff and sediment loading at the field level and routing through stream networks. Typically, the SWAT model first needs to be calibrated and validated for accurate estimation through adjustment of sensitive input parameters (i.e., Curve Number values, USLE P, slope and slope‐length, and so on). However, in some instances, SWAT‐simulated results are greatly affected by the watershed delineation and Digital Elevation Models (DEM) cell size. In this study, the SWAT ArcView GIS Patch II was developed for steep sloping watersheds, and its performance was evaluated for various threshold values and DEM cell size scenarios when delineating subwatersheds using the SWAT model. The SWAT ArcView GIS Patch II was developed using the ArcView GIS Avenue program and Spatial Analyst libraries. The SWAT ArcView GIS Patch II improves upon the SWAT ArcView GIS Patch I because it reflects the topographic factor in calculating the field slope‐length of Hydrologic Response Units in the SWAT model. The simulated sediment value for 321 subwatersheds (watershed delineation threshold value of 25 ha) is greater than that for 43 subwatersheds (watershed delineation threshold value of 200 ha) by 201% without applying the SWAT ArcView GIS Patch II. However, when the SWAT ArcView GIS Patch II was applied, the difference in simulated sediment yield decreases for the same scenario (i.e., difference in simulated sediment with 321 subwatersheds and 43 subwatersheds) was 12%. The simulated sediment value for DEM cell size of 50 m is greater than that for DEM cell size of 10 m by 19.8% without the SWAT ArcView GIS Patch II. However, the difference becomes smaller (3.4% difference) between 50 and 10 m with the SWAT ArcView GIS Patch II for the DEM scenarios. As shown in this study, the SWAT ArcView GIS Patch II can reduce differences in simulated sediment values for various watershed delineation and DEM cell size scenarios. Without the SWAT ArcView GIS Patch II, variations in the SWAT‐simulated results using various watershed delineation and DEM cell size scenarios could be greater than those from input parameter calibration. Thus, the results obtained in this study show that the SWAT ArcView GIS Patch II should be used when simulating hydrology and sediment yield for steep sloping watersheds (especially if average slope of the subwatershed is >25%) for more accurate simulation of hydrology and sediment using the SWAT model. The SWAT ArcView GIS Patch II is available at http://www.EnvSys.co.kr/~swat for free download.     

SEDIMENT DEPOSITION ESTIMATION FROM CESIUM-137 MEASUREMENTS1

ABSTRACT The movement of fallout ¹³⁷Cs carried by soil particles was studied as an indicator of erosion and sedimentation in the Allerton watersheds and 4-H Memorial Lake located near Monticello, Illinois. Sediment deposition was greater in the waterway draining from watershed IB than in the waterway from watershed IA. At the average rate of 2.3 cm/yr of sediment deposition in the lake (from 1954 to 1979), there will be a loss of over 2 meters of water depth in the next century. However, there appears to be a decreasing rate of sediment deposition in the 4-H Memorial Lake as a result of improved conservation practices on the watersheds and the increased effectiveness of vegetated waterways and buffers for retaining sediment.     

RESERVOIR SEDIMENTATION RATES LINKED TO LONG-TERM CHANGES IN AGRICULTURAL LANI) USE1

ABSTRACT: Long-term land use and reservoir sedimentation were quantified and linked in a small agricultural reservoir-watershed system without having historical data. Land use was determined from a time sequence of aerial photographs, and reservoir sedimentation was determined from cores with ¹³⁷Cs dating techniques. They were linked by relating sediment deposition to potential sediment production which was determined by the Universal Soil Loss Equation and by SCS estimates for gullied land. Sediment cores were collected from Tecumseh Lake, a 55-ha reservoir with a 1,189-ha agricultural watershed, constructed in 1934 in central Oklahoma. Reservoir sediment deposition decreased from an average of 5,933 Mg/yr from 1934 to 1954, to 3,179 Mg/yr from 1954 to 1962, and finally to 1,017 Mg/yr from 1962 to 1987. Potential sediment production decreased from an average of 29,892 to 11,122 and then to 3,589 Mg/yr for the same time periods as above, respectively. Reductions in deposition and sediment production corresponded to reductions in cultivated and abandoned cropland which became perennial pasture. Together, cultivated and abandoned cropland accounted for 59 percent of the watershed in 1937, 24 percent in 1954, and 10 percent in 1962. Roadway erosion, stream bank erosion, stored stream channel sediment, and long-term precipitation were considered, but none seemed to play a significant role in changing sediment deposition rates. Instead, the dominant factor was the conversion of fields to perennial pastures. The effect of conservation measures on reservoir sedimentation can now be quantified for many reservoirs where historical data is not available.     

Test of APEX for nine forested watersheds in East Texas

Hydrologic/water quality models are increasingly used to explore management and policy alternatives for managing water quality and quantity from intensive silvicultural practices with best management practices (BMPs) in forested watersheds due to the limited number of and cost of conducting watershed monitoring. The Agricultural Policy/Environmental eXtender (APEX) model was field-tested using 6 yr of data for flow, sediment, nutrient, and herbicide losses collected from nine small (2.58 to 2.74 ha) forested watersheds located in southwest Cherokee County in East Texas. Simulated annual average stream flow for each of the nine watersheds was within +/- 7% of the corresponding observed values; simulated annual average sediment losses were within +/- 8% of measured values for eight out of nine watersheds. Nash-Sutcliffe efficiency (EF) values ranged from 0.68 to 0.94 based on annual stream flow comparison and from 0.60 to 0.99 based on annual sediment comparison. Similar to what was observed, simulated flow, sediment, organic N, and P were significantly increased on clear-cut watersheds compared with the control watersheds. APEX reasonably simulated herbicide losses, with an EF of 0.73 and R(2) of 0.74 for imazapyr, and EF of 0.65 and R(2) of 0.68 for hexazinone based on annual values. Overall, the results show that APEX was able to predict the effects of silvicultural practices with BMPs on water quantity and quality and that the model is a useful tool for simulating a variety of responses to forest conditions.     

THE USE OF BINARY OPTIMIZATION AND HYDROLOGIC MODELS TO FORM RIPARIAN BUFFERS1

ABSTRACT: Riparian buffers are considered important management options for protecting water quality. Land costs and buffer performance, which are functions of local environmental characteristics, are likely to be key attributes in the selection process, especially when budgets are limited. In this article we demonstrate how a framework involving hydrologic models and binary optimization can be used to find the optimal buffer subject to a budget constraint. Two hydrologic models, SWAT and REMM, were used to predict the loads from different source areas with and without riparian buffers. These loads provided inputs for a binary optimization model to select the most cost efficient parcels to form a riparian buffer. This methodology was applied in a watershed in Delaware County, New York. The models were parameterized using readily available digital databases and were later compared against observed flow and water quality data available for the site. As a result of the application of this method, the marginal utility of incremental increases in buffer widths along the stream channel and the set of parcels to form the best affordable riparian buffer were obtained.     

Assessing Sediment Yield for Selected Watersheds in the Laurentian Great Lakes Basin Under Future Agricultural Scenarios

In the Laurentian Great Lakes Basin (GLB), corn acreage has been expanding since 2005 in response to high demand for corn as an ethanol feedstock. This study integrated remote sensing-derived products and the Soil and Water Assessment Tool (SWAT) within a geographic information system (GIS) modeling environment to assess the impacts of cropland change on sediment yield within four selected watersheds in the GLB. The SWAT models were calibrated during a 6 year period (2000–2005), and predicted stream flows were validated. The R ² values were 0.76, 0.80, 0.72, and 0.81 for the St. Joseph River, the St. Mary River, the Peshtigo River, and the Cattaraugus Creek watersheds, respectively. The corresponding E (Nash and Sutcliffe model efficiency coefficient) values ranged from 0.24 to 0.79. The average annual sediment yields (tons/ha/year) ranged from 0.12 to 4.44 for the baseline (2000 to 2008) condition. Sediment yields were predicted to increase for possible future cropland change scenarios. The first scenario was to convert all “other” agricultural row crop types (i.e., sorghum) to corn fields and switch the current/baseline crop rotation into continuous corn. The average annual sediment yields increased 7–42 % for different watersheds. The second scenario was to further expand the corn planting to hay/pasture fields. The average annual sediment yields increased 33–127 % compared with baseline conditions.     

Temporal and Spatial Relationships Between Watershed Land Use and Salt Marsh Disturbance in a Pacific Estuary

Historical and recent remote sensing data can be used to address temporal and spatial relationships between upland land cover and downstream vegetation response at the watershed scale. This is demonstrated for sub-watersheds draining into Elkhorn Slough, California, where salt marsh habitat has diminished because of the formation of sediment fans that support woody riparian vegetation. Multiple regression models were used to examine which land cover variables and physical properties of the watershed most influenced sediment fan size within 23 sub-watersheds (1.4 ha to 200 ha). Model explanatory power increased (adjusted R(2) = 0.94 vs. 0.75) among large sub-watersheds (>10 ha) and historical watershed variables, such as average farmland slope, flowpath slope, and flowpath distance between farmland and marsh, were significant. It was also possible to explain the increase in riparian vegetation by historical watershed variables for the larger sub-watersheds. Sub-watershed area is the overriding physical characteristic influencing the extent of sedimentation in a salt marsh, while percent cover of agricultural land use is the most influential land cover variable. The results also reveal that salt marsh recovery depends on relative cover of different land use classes in the watershed, with greater chances of recovery associated with less intensive agriculture. This research reveals a potential delay between watershed impacts and wetland response that can be best revealed when conducting multi-temporal analyses on larger watersheds.     

Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed1

Cho, Jaepil, Richard R. Lowrance, David D. Bosch, Timothy C. Strickland, Younggu Her, and George Vellidis, 2010. Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed. Journal of the American Water Resources Association (JAWRA) 46(3):586-602. DOI: 10.1111/j.1752-1688.2010.00436.x Abstract: The Soil and Water Assessment Tool (SWAT) does not fully simulate riparian buffers, but has a simple filter function that is responsive to filter strip width (FILTERW). The objectives of this study were to (1) evaluate SWAT hydrology and water quality response to changes in watershed subdivision levels and different FILTERW configurations and (2) provide guidance for selecting appropriate watershed subdivision for model runs that include the riparian buffer feature through the FILTERW parameter. Watershed subdivision level is controlled by the critical source area (CSA) which defines the minimum drainage area required to form the origin of a stream. SWAT was calibrated on a 15.7 km² subdrainage within the Little River Experimental Watershed, Georgia. The calibrated parameter set was applied to 32 watershed configurations consisting of four FILTERW representations for each of eight CSA levels. Streamflow predictions were stable regardless of watershed subdivision and FILTERW configuration. Predicted sediment and nutrient loads from upland areas decreased as CSA increased when spatial variations of riparian buffers are considered. Sediment and nutrient yield at the watershed outlet was responsive to different combinations of CSA and FILTERW depending on selected in-stream processes. CSA ranges which provide stable sediment and nutrient yields at the watershed outlet was suggested for avoiding significant modifications in selected parameter set.     

Reduction of stream sediment concentration by a riparian buffer: filtering of road runoff in disturbed headwater basins of montane mainland southeast Asia

We determined the extent that a riparian buffer reduces stream suspended sediment concentrations by filtering road runoff during 18 rain events in a 2.5-ha, multi-use watershed in northern Thailand. The dominant buffer species was the perennial sedge Fimbristylis aphylla Zoll. ex Steud. (Cyperaceae). We monitored stream sediment concentration for situations where road runoff either flowed into the riparian buffer or was diverted directly into the stream (buffer and no buffer scenarios). These data were used to develop the following relationships between instantaneous stream sediment concentration (Ci) and discharge (Qi): Ci= 28.329Qi(0.851) (buffer scenario) and Ci= 22.265Qi(1.579) (no buffer scenario). Using these functions to calculate total event suspended concentrations, we determined that the buffer reduced suspended sediment concentration by 34 to 87%, for the range of events monitored. Removal of sediment from runoff generated on a 2.4-m-wide, 165-m-long unpaved road section was achieved principally via ponding, which reduced the transport capacity as flow entered the relatively flat, saturated buffer. Sediment deposition occurred primarily within the first 10 m of the buffer. Some sediment was also deposited on the fillslope leading to the buffer. Maximum road sediment concentration during the largest buffer event approached 100,000 mg L(-1). Meanwhile, the corresponding maximum stream suspended sediment concentration was <4000 mg L(-1). In contrast, maximum stream concentrations when flow bypassed the buffer during smaller events were commonly 4000 to 7000 mg L(-1). Naturally occurring buffers represent an economical means of mitigating road-related impacts in upland basins in Southeast Asia, particularly if combined with measures limiting sediment and runoff production on contributing road sections.     

EXAMINING LAND USE INFLUENCES ON STREAM HABITATS AND MACROINVERTEBRATES: A GIS APPROACH1

ABSTRACT: Geographic Information Systems (GIS) were used to assess the relationships between land use patterns and the physical habitat and macroinvertebrate fauna of streams within similar sized watersheds. Eleven second or third order watersheds ranging from highly urbanized to heavily forested were selected along Lake Superior's North Shore. Land use patterns within the watersheds were quantified using readily available digital land use/land cover information, with a minimum mapping resolution of 16 ha. Physical habitat features, describing substrate characteristics and stream morphology, were characterized at sample points within each stream. Principle component and correlation analyses were used to identify relationships between macroinvertebrates and stream physical habitat, and between habitat and land use patterns. Substrate characteristics and presence of coarse woody debris were found to have the strongest correlations with macreinvertebrate assemblage richness and composition. Agricultural and urban land use was correlated with substrate characteristics. Algal abundance, associated with macroinvertebrate compositional differences, was correlated with housing density and non-forest land covers. The use of readily available spatial data, even at this relatively coarse scale, provides a means to detect the primary relationships between land use and stream habitat quality; finer-resolution GIS databases are needed to assess more subtle influences, such as those due to riparian conditions.     

Landscape Planning for Agricultural Nonpoint Source Pollution Reduction III: Assessing Phosphorus and Sediment Reduction Potential

Riparian buffers have the potential to improve stream water quality in agricultural landscapes. This potential may vary in response to landscape characteristics such as soils, topography, land use, and human activities, including legacies of historical land management. We built a predictive model to estimate the sediment and phosphorus load reduction that should be achievable following the implementation of riparian buffers; then we estimated load reduction potential for a set of 1598 watersheds (average 54 km²) in Wisconsin. Our results indicate that land cover is generally the most important driver of constituent loads in Wisconsin streams, but its influence varies among pollutants and according to the scale at which it is measured. Physiographic (drainage density) variation also influenced sediment and phosphorus loads. The effect of historical land use on present-day channel erosion and variation in soil texture are the most important sources of phosphorus and sediment that riparian buffers cannot attenuate. However, in most watersheds, a large proportion (approximately 70%) of these pollutants can be eliminated from streams with buffers. Cumulative frequency distributions of load reduction potential indicate that targeting pollution reduction in the highest 10% of Wisconsin watersheds would reduce total phosphorus and sediment loads in the entire state by approximately 20%. These results support our approach of geographically targeting nonpoint source pollution reduction at multiple scales, including the watershed scale.     

INTEGRATED ECHO SOUNDER,GPS, AND GIS FOR RESERVOIR SEDIMENTATION STUDIES: EXAMPLES FROM TWO ARKANSAS LAKES1

ABSTRACT: Bathymetric and sedimentation surveys were conducted using a dual frequency (28/200 kHz) echo sounder system in two reservoirs (Lee Creek Reservoir and Lake Shepherd Springs) in the Ozark Plateau of northwestern Arkansas. Echo sounder survey data were merged within geographic information system (GIS) software to provide detailed visualization and analyses of current depths, pre‐impoundment topography, distribution, thickness, and volume estimates of lacustrine sediment, time averaged sediment accumulation rates, long term average annual sediment flux, and water storage capacity. Calculated long term average sediment accumulation rates were used to model sediment infilling and projected lifetimes of each reservoir. Results from echo sounder surveys and GIS analyses suggest that the Lee Creek Reservoir has a projected lifetime of approximately 500 years compared to a projected lifetime for Lake Shepherd Springs of approximately 3,000 years. Estimated differences in projected lifetimes of these reservoirs reflected differences in initial reservoir volume and long term average annual sediment flux from the respective watersheds related to watershed area, physiography, land cover, and land use. The universal soil loss equation (USLE) model generated sediment fluxes an order of magnitude larger from the watersheds of both reservoirs compared to the geophysical data estimates. This study demonstrated the utility of merging geophysical survey (echo sounder) data within a GIS as an aid to understanding patterns of reservoir sedimentation. These data and analyses also provide a baseline relevant to understanding sedimentation processes and are necessary for development of long term management plans for these reservoirs and their watersheds.     

ROLE OF WATERSHED SUBDIVISION ON MODELING THE EFFECTIVENESS OF BEST MANAGEMENT PRACTICES WITH SWAT1

Distributed parameter watershed models are often used for evaluating the effectiveness of various best management practices (BMPs). Streamflow, sediment, and nutrient yield predictions of a watershed model can be affected by spatial resolution as dictated by watershed subdivision. The objectives of this paper are to show that evaluation of BMPs using a model is strongly linked to the level of watershed subdivision; to suggest a methodology for identifying an appropriate subdivision level; and to examine the efficacy of different BMPs at field and watershed scales. In this study, the Soil and Water Assessment Tool (SWAT) model was calibrated and validated for streamflow, sediment, and nutrient yields at the outlet of the Dreisbach (623 ha) and Smith Fry (730 ha) watersheds in Maumee River Basin, Indiana. Grassed waterways, grade stabilization structures, field borders, and parallel terraces are the BMPs that were installed in the study area in the 1970s. Sediment and nutrient outputs from the calibrated model were compared at various watershed subdivision levels, both with and without implementation of these BMPs. Results for the study watersheds indicated that evaluation of the impacts of these BMPs on sediment and nutrient yields was very sensitive to the level of subdivision that was implemented in SWAT. An optimal watershed subdivision level for representation of the BMPs was identified through numerical simulations. For the study watersheds, it would appear that the average subwatershed area corresponding to approximately 4 percent of total watershed area is needed to represent the influence of these BMPs when using the SWAT model.     

Nitrate and sediment fluxes from a California rangeland watershed

Long-term water quality records for assessing natural variability, impact of management, and that guide regulatory processes to safeguard water resources are rare for California oak woodland rangelands. This study presents a 20-yr record (1981-2000) of nitrate-nitrogen (NO(3)-N) and suspended sediment export from a typical, grazed oak woodland watershed (103 ha) in the northern Sierra Nevada foothills of California. Mean annual precipitation over the 20-yr period was 734 mm yr(-1) (range 366-1205 mm yr(-1)). Mean annual stream flow was 353 mm y(-1) (range 87-848 mm yr(-1)). Average annual stream flow was 48.1 +/- 16% of precipitation. Mean annual NO(3)-N export was 1.6 kg ha(-1) yr(-1) (range 0.18-3.6 kg ha(-1) yr(-1)). Annual NO(3)-N export significantly (P < 0.05) increased with increasing annual stream flow and precipitation. Mean daily NO(3)-N export was 0.004 kg ha(-1) d(-1) (range 10(-5) to 0.55 kg ha(-1) d(-1)). Mean annual suspended sediment export was 198 kg ha(-1) yr(-1) (range 23-479 kg ha(-1) yr(-1)). There was a positive relationship (P < 0.05) between annual suspended sediment export, annual stream flow and precipitation. Mean daily suspended sediment export was 0.54 kg ha(-1) d(-1) (range 10(-4) to 155 kg ha(-1) d(-1)). Virtually no sediment was exported during the dry season. The large variation in daily and annual fluxes highlights the necessity of using long-term records to establish quantitative water quality targets for rangelands and demonstrates the difficulty of designing a water quality monitoring program for these ecosystems.     

Testing of the Modified Streambank Erosion and Instream Phosphorus Routines for the SWAT Model

In some watersheds, streambanks are a source of two major pollutants, phosphorus (P) and sediment. P originating from both uplands and streambanks can be transported and stored indefinitely on floodplains, streambanks, and in closed depressions near the stream. The objectives of this study were to (1) test the modified streambank erosion and instream P routines for the Soil and Water Assessment Tool (SWAT) model in the Barren Fork Creek watershed in northeast Oklahoma, (2) predict P in the watershed with and without streambank‐derived P, and (3) determine the significance of streambank erosion P relative to overland P sources. Measured streambank and channel parameters were incorporated into a flow‐calibrated SWAT model and used to estimate streambank erosion and P for the Barren Fork Creek using modified streambank erosion and instream P routines. The predicted reach‐weighted streambank erosion was 40 kg/m vs. the measured 42 kg/m. Streambank erosion contributed 47% of the total P to the Barren Fork Creek and improved P predictions compared to observed data, especially during the high‐flow events. Of the total P entering the stream system, approximately 65% was removed via the watershed outlet and 35% was stored in the floodplain and stream system. This study successfully applied the SWAT model's modified streambank erosion and instream P routines and demonstrated that streambank‐derived P can improve P modeling at the watershed scale. 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.