GENERALIZED WATERSHED LOADING FUNCTIONS FOR STREAM FLOW NUTRIENTS1

ABSTRACT: Loading functions are proposed as a general model for estimating monthly nitrogen and phosphorus fluxes in stream flow. The functions have a simple mathematical structure, describe a wide range of rural and urban nonpoint sources, and couple surface runoff and ground water discharge. Rural runoff loads are computed from daily runoff and erosion and monthly sediment yield calculations. Urban runoff loads are based on daily nutrient accumulation rates and exponential wash off functions. Ground water discharge is determined by lumped parameter unsaturated and saturated zone soil moisture balances. Default values for model chemical parameters were estimated from literature values. Validation studies over a three-year period for an 850 km² watershed showed that the loading functions explained at least 90 percent of the observed monthly variation in dissolved and total nitrogen and phosphorus fluxes in stream flow. Errors in model predictions of mean monthly fluxes were: dissolved phosphorus - 4 percent; total phosphorus - 2 percent; dissolved nitrogen - 18 percent; and total nitrogen - 28 percent. These results were obtained without model calibration.     

MODELING OF NONPOINT SOURCE POLLUTION OF NITROGEN AT THE WATERSHED SCALE1

ABSTRACT: The Watershed Nutrient Transport and Transformation (NTT-Watershed) model is a physically based, energy-driven, multiple land use, distributed model that is capable of simulating water and nutrient transport in a watershed. The topographic features and subsurface properties of the watershed are refined into uniform, homogeneous square grids. The vertical discretization includes vegetation, overland flow, soil water redistribution and groundwater zones. The chemical submodel simulates the nitrogen dynamics in terrestrial and aquatic systems. Three chemical state variables are considered (NO₃^-_-, NH₄⁺, and Org-N). The NTT-Watershed model was used to simulate the fate and transport of nitrogen in the Muddy Brook watershed in Connecticut. The model was shown to be capable of capturing the hydrologic and portions of the nitrogen dynamics in the watershed. Watershed planners could use this model in developing strategies of best management practices that could result in maximizing the reductions of nitrogen export from a watershed.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

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.     

HYDROLOGIC CALIBRATION OF THE SWAT MODEL IN A WATERSHED CONTAINING FRAGIPAN SOILS1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) model, designed for use on rural ungaged basins and incorporating a GRASS GIS interface, was used to model the hydrologic response of the Ariel Creek watershed of northeastern Pennsylvania. Model evaluation of daily flow prior to calibration revealed a deviation of runoff volumes (D_v) of 68.3 percent and a Nash-Sutcliffe coefficient of-0.03. Model performance was affected by unusually large observed snowmelt events and the inability of the model to accurately simulate baseflow, which was influenced by the presence of fragipans. Seventy-five percent of the soils in the watershed contain fragipans. Model calibration yielded a D_v of 39.9 percent and a Nash-Sutcliffe coefficient of 0.04, when compared on a daily basis. Monthly comparisons yielded a Nash-Sutcliffe coefficient of 0.14. Snowmelt events in the springs of 1993 and 1994, which were unusually severe, were not adequately simulated. Neglecting these severe events, which produced the largest and third largest measured flows for the period of record, a D_v of 4.1 percent and Nash-Sutcliffe coefficient of 0.20 were calculated on a daily comparison, while on a monthly basis the Nash-Sutciffe coefficient was 0.55. These results suggest that the SWAT model is better suited to longer period simulations of hydrologic yields. Baseflow volumes were accurately simulated after calibration (D_v= -0.2 percent). Refinements made to the algorithms controlling subsurface hydrology and snowmelt, to better represent the presence of fragipans and snowmelt events, would likely improve model performance.     

EXPORT COEFFICIENT MODELING TO ASSESS PHOSPHORUS LOADING IN AN URBAN WATERSHED1

ABSTRACT: An export coefficient modeling approach was used to assess the influence of land use on phosphorus loading to a Southern Ontario stream. A model was constructed for the 1995–1996 water year and calibrated within ± 3 percent of the observed mean concentration of total phosphorus. It was found that runoff from urban areas contributed most to the loading of phosphorus to the stream. When the model was assessed by running it for the 1977–1978 water year, using water quality and land use data collected independently, agreement within ± 7 percent was obtained. The model was then used to forecast the impact of future urban development proposed for the watershed, in terms of phosphorus loading, and to evaluate the reduction in loading resulting from several urban best management practices (BMP). It was determined that phosphorus removal will have to be applied to all the urban runoff from the watershed to appreciably reduce stream phosphorus concentration. Of the BMP designs assessed, an infiltration pond system resulted in the greatest phosphorus load reduction, 50 percent from the 1995–1996 baseline.     

EVALUATION OF Ann AGNPS NITROGEN LOADING IN AN AGRICULTURAL WATERSHED1

ABSTRACT: Pollutants entering a water system can be very destructive to the health of that system. Best Management Practices (BMPs) are used to reduce these pollutants, but understanding the most effective practices is very difficult. Watershed models are an effective tool to aid in the decision‐making process of selecting the BMPs that are most effective in reducing the pollutant loading and are also the most cost effective. The Annualized Agricultural Nonpoint Source Pollution model (AnnAGNPS 2.0) is a technological tool that can be used to estimate watershed response to agricultural management practices. The main purpose of this paper is to test the performance of AnnAGNPS 2.0 on nitrogen loading using comparisons with measurements from the Deep Hollow watershed of the Mississippi Delta Management Systems Evaluation Area (MDMSEA) project. Previous work has demonstrated the capability of the model to simulate runoff and sediment. From sensitivity analyses in this study, initial nitrogen concentration in the soil and crop nitrogen uptake had the most impact on the nitrogen loadings. AnnAGNPS simulations of monthly nitrogen loadings are poor. However, statistical test (t‐test) showed that the predicted nitrogen loading is not significantly different from observed nitrogen loading at the 95 percent level of confidence.     

WATERSHED AND LAKE WATER QUALITY ASSESSMENT: AN INTEGRATED MODELING APPROACH1

Watershed planning groups and action agencies seek to understand how lake water quality responds to changes in watershed management. This study developed and demonstrated the applicability of an integrated modeling approach for providing this information. An integrated model linking watershed conditions to water-quality of the receiving lake incorporated the following components: (1) an event-based AGNPS model to estimate watershed pollutant losses; (2) annualization of AGNPS results to produce annual lake pollutant loadings; (3) a base flow separation package, SAM, to estimate base flow; (4) estimates of nutrients in base flow and point sources; and (5) linkage of watershed loadings directly to EUTROMOD lake water quality algorithms. Results are presented for Melvern Lake, a 28-km² multipurpose reservoir with a 900-km² agricultural watershed in east central Kansas. Reasonable estimates of current lake quality were attained using an average phosphorus availability factor of 31 percent to calibrate model results to measured in-lake phosphorus. Comparison of a range of possible scenarios, including all cropland changed to no-till (best case) and all CRP and good-condition grasslands changed to cropland (worst case), indicated only a (4 percent change for in-lake phosphorus and a (2 percent change for chlorophyll a. These results indicated that this watershed is not sensitive to projected changes in land use and management.     

A FRAMEWORK FOR PHOSPHORUS TRANSPORT MODELING IN THE LAKE OKEECHOBEE WATERSHED1

ABSTRACT: A modeling framework was developed to determine phosphorus loadings to Lake Okeechobee from watersheds located north of the lake. This framework consists of the land-based model CREAMS-WT, the in-stream transport model QUAL2E, and an interface procedure to format the land-based model output for use by the in-stream model. QUAL2E hydraulics and water quality routines were modified to account for flow routing and phosphorus retention in both wetlands and stream channels. Phosphorus loadings obtained from previous applications of CREAMS-WT were used by QUAL2E, and calibration and verification showed that QUAL2E accurately simulated seasonal and annual phosphorus loadings from a watershed. Sensitivity and uncertainty analyses indicated that the accuracy of monthly loadings can be improved by using better estimates of in-stream phosphorus decay rates, ground water phosphorus concentrations, and runoff phosphorus concentrations as input to QUAL2E.     

GIS‐BASED WATER QUALITY MODELING IN THE SANDUSKY WATERSHED,OHIO, USA1

ABSTRACT: This study focused on the Sandusky Watershed (SW) in Ohio, located within the Great Lakes Basin, with emphasis on two of its subwatersheds, namely Honey Creek (HC) and Rock Creek (RC). The goal was to assess the capabilities of the Soil and Water Assessment Tool (SWAT) to simulate suspended sediment (SS), phosphorus (P) and nitrogen (N) yield in the SW that contribute major sediment and nutrient loads into Lake Erie. The model was calibrated using water flow and water quality parameters for water years 1998 to 1999 and validated model simulations covering the period of water years 2000 to 2001 for monthly conditions. The validation of SS showed correlation coefficients of 0.29 (SW), 0.75 (HC) and 0.69 (RC). Correlation coefficients for P were 0.68 (SW), 0.78 (HC) and 0.37 (RC); for N0₂‐N 0.84 (HC) and 0.38 (RC); for N0₃‐N 0.27 (HC) and 0.76 (RC); for NH₃‐N 0.57 (SW), 0.49 (HC), and 0.13 (RC). In addition, mean errors, root mean square errors, Nash‐Sutcliffe coefficients, and graphs were used to compare simulated to measured data. Simulation success was variable with poor and good simulations, but in most cases, simulated water quality values followed the trend of measured data except for extreme (or intense) rainfall/runoff events. Reviews of 17 applications indicated that the SWAT is suitable for long term continuous simulations but not for storm events. A spatially distributed modeling approach generated maps showing the spatial distribution of SS, P, and N for each simulation element across the Sandusky Watershed.     

MODELING THE BUILDUP AND WASHOFF OF POLLUTANTS ON URBAN WATERSHEDS1

ABSTRACT: A model for urban stormwater quality was developed in this study. The basis for the model is the process by which pollutants build up on the watershed surface. For the wet climate of the study site, it was assumed that there exists an interval of time over which the pollutant buildup equals the pollutant washoff (no accumulation of pollutant). The buildup model was represented by a linear function of the antecedent dry time. The buildup function was then linked with a pollutant washoff model represented by a power function of the storm runoff volume. Various time intervals for no net accumulation were tested to calibrate the model. The model was calibrated to observed data for two small urban basins in Baton Rouge, Louisiana, and model results were used to analyze the behavior of phosphorus concentrations in storm runoff from these basins over a long period of time.     

DETERMINATION OF BEST TIMING FOR POULTRY WASTE DISPOSAL: A MODELING APPROACH1

ABSTRACT: Confined production of poultry results in significant volumes of waste material which are typically disposed of by land application. Concerns over the potential environmental impacts of poultry waste disposal have resulted in ongoing efforts to develop management practices which maintain high quality of water downstream of disposal areas. The timing of application to minimize waste constituent losses is a management practice with the potential to ensure high quality of streams, rivers, and lakes downstream of receiving areas. This paper describes the development and application of a method to identify which time of year is best, from the standpoint of surface water quality, for land application of poultry waste. The procedure consists of using a mathematical simulation model to estimate average nitrogen and phosphorus losses resulting from different application timings, and then identifying the timings which minimize losses of these nutrients. The procedure was applied to three locations in Arkansas, and three different criteria for optimality of application timing were investigated. One criterion was oriented strictly to water quality, one was oriented only to crop production, and the last was a combination. The criteria resulted in different windows of time being identified as optimal. Optimal windows also varied with location of the receiving area. The results indicate that it is possible to land-apply poultry waste at times which both minimize nutrient losses and maximize crop yield.     

IMPLICATIONS OF ON‐GROUND NITROGEN LOADING AND SOIL TRANSFORMATIONS ON GROUND WATER QUALITY MANAGEMENT1

ABSTRACT: This paper presents a modeling approach based on a geographic information system (GIS) to estimate the variability of on‐ground nitrogen loading and the corresponding nitrate leaching to ground water. The methodology integrates all point and nonpoint sources of nitrogen, the national land cover database, soil nitrogen transformations, and the uncertainty of key soil and land use‐related parameters to predict the nitrate mass leaching to ground water. The analysis considered 21 different land use classes with information derived from nitrogen sources such as fertilizer and dairy manure applications, dairy lagoons, septic systems, and dry and wet depositions. Simulations were performed at a temporal resolution of one month to capture seasonal trends. The model was applied to a large aquifer of 376 square miles in Washington State that serves more than 100,000 residents with drinking water. The results showed that dairy manure is the main source of nitrogen in the area followed by fertilizers. It was also seen that nitrate leaching is controlled by the recharge rate, and there can be a substantial buildup of soil nitrogen over long periods of time. Uncertainty analysis showed that denitrification rate is the most influential parameter on nitrate leaching. The results showed that combining management alternatives is a successful strategy, especially with the use of nitrification inhibitors. Also, change in the land use pattern has a noticeable impact on nitrate leaching.     

THE EFFECTS OF CLIMATE CHANGE ON STREAM FLOW AND NUTRIENT LOADING'

ABSTRACT: This study assesses the potential impact of climate change on stream flow and nutrient loading in six watersheds of the Susquehanna River Basin using the Generalized Watershed Loading Function (GWLF). The model was used to simulate changes in stream flow and nutrient loads under a transient climate change scenario for each watershed. Under an assumption of no change in land cover and land management, the model was used to predict monthly changes in stream flow and nutrient loads for future climate conditions. Mean annual stream flow and nutrient loads increased for most watersheds, but decreased in one watershed that was intensively cultivated. Nutrient loading slightly decreased in April and late summer for several watersheds as a result of early snowmelt and increasing evapotranspiration. Spatial and temporal variability of stream flow and nutrient loads under the transient climate scenario indicates that different approaches for future water resource management may be useful.     

LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART II: MODEL APPLICATION1

ABSTRACT: This paper describes the application of a river basin scale hydrologic model (described in Part I) to Richland and Chambers Creeks watershed (RC watershed) in upper Trinity River basin in Texas. The inputs to the model were accumulated from hydro-graphic and geographic databases and maps using a raster-based GIS. Available weather data from 12 weather stations in and around the watershed and stream flow data from two USGS stream gauge station for the period 1965 to 1984 were used in the flow calibration and validation. Sediment calibration was carried out for the period 1988 through 1994 using the 1994 sediment survey data from the Richland-Chambers lake. Sediment validation was conducted on a subwatershed (Mill Creek watershed) situated on Chambers Creek of the RC watershed. The model was evaluated by well established statistical and visual methods and was found to explain at least 84 percent and 65 percent of the variability in the observed stream flow data for the calibration and validation periods, respectively. In addition, the model predicted the accumulated sediment load within 2 percent and 9 percent from the observed data for the RC watershed and Mill Creek watershed, respectively.     

PHYSICAL DETERMINANTS OF ECONOMIC VALUE OF RIPARIAN BUFFERS IN AN AGRICULTURAL WATERSHED1

ABSTRACT: Economic values of riparian buffers in a watershed are evaluated by the changes in the net economic return for farming with and without riparian buffers when achieving the same water quality objectives. The simulated water quality impacts of alternative farming systems using SWAT and experimental data for riparian buffers are used in a mathematical optimization model to estimate net economic return for farming subject to a water quality objective. Physical characteristics such as stream length, channel slope, average land slope, cropland percentage and several soil attributes are identified in the watershed using ARC/INFO GIS. A regression model is then used to evaluate the impacts of these physical characteristics on the estimated economic values of buffers. The study is conducted in Goodwater Creek watershed, Missouri. The results show the estimated economic value of buffers is significantly affected by some soil properties, stream length, and cropland percentage in watershed and can be used to improve the effectiveness of riparian buffers at watershed and regional levels.     

PRIORITIZING NONPOINT SOURCE PHOSPHORUS LOADING USING A GRASS-MODELING SYSTEM1

ABSTRACT: A UNIX-based windows application was developed to integrate a phosphorus transport model with the Geographic Resources Analysis Support System (GRASS). The system prioritizes potential phosphorus loading from fields or cells in a watershed and can evaluate the effects of alternate management practices on phosphorus yield. The model simulates phosphorus loading by using a daily mass balance on a unit-area basis and incorporates the effects from rainfall, topography, soil properties, fertilizer and animal waste application, and management. Model predictions include dissolved and sediment-bound phosphorus yield, runoff volume, and sediment yield. Within the integrated GRASS-modeling system, the user can obtain model input data, execute the model using various options, rank model input and out. put data, and display them as GRASS-based maps or data tables. All functions are menu driven, developed using C language and X-window tools to run on a SUN workstation platform. The system provides a powerful and efficient tool for prioritizing phosphorus loading from nonpoint sources.     

IDENTIFICATION OF CRITICAL NONPOINT POLLUTION SOURCE AREAS USING GEOGRAPHIC INFORMATION SYSTEMS AND WATER QUALITY MODELING1

n integrated approach coupling water quality computer simulation modeling with a geographic information system (GIS) was used to delineate critical areas of nonpoint source (NPS) pollution at the watershed level. Two simplified pollutant export models were integrated with the Virginia Geographic Information System (VirGIS) to estimate soil erosion, sediment yield, and phosphorus (P) loading from the Nomini Creek watershed located in Westmoreland County, Virginia. On the basis of selected criteria for soil erosion rate, sediment yield, and P loading, model outputs were used to identily watershed areas which exhibit three categories (low, medium, high) of non-point source pollution potentials. The percentage of the watershed area in each category, and the land area with critical pollution problems were also identified. For the 1505-ha Nomini Creek watershed, about 15, 16, and 21 percent of the watershed area were delineated as sources of critical soil erosion, sediment, and phosphorus pollution problems, respectively. In general, the study demonstrated the usefulness of integrating GIS with simulation modeling for nonpoint source pollution control and planning. Such techniques can facilitate making priorities and targeting nonpoint source pollution control programs.     

MODELING THE DISTRIBUTION OF DIFFUSE NITROGEN SOURCES AND SINKS IN THE NEUSE RIVER BASIN OF NORTH CAROLINA,USA1

This study quantified nonpoint source nitrogen (NPS‐N) sources and sinks across the 14,582 km² Neuse River Basin (NRB) located in North Carolina, to provide tabular data summaries and graphic overlay products to support the development of management approaches to best achieve established N reduction goals. First, a remote sensor derived, land cover classification was performed to support modeling needs. Modeling efforts included the development of a mass balance model to quantify potential N sources and sinks, followed by a precipitation event driven hydrologic model to effectively transport excess N across the landscape to individual stream reaches to support subsequent labeling of transported N values corresponding to source origin. Results indicated that agricultural land contributed 55 percent of the total annual NPS‐N loadings, followed by forested land at 23 percent (background), and urban areas at 21 percent. Average annual N source contributions were quantified for agricultural (1.4 kg/ha), urban (1.2 kg/ha), and forested cover types (0.5 kg/ha). Nonpoint source‐N contributions were greatest during the winter (40 percent), followed by spring (32 percent), summer (28 percent), and fall (0.3 percent). Seasonal total N loadings shifted from urban dominated and forest dominated sources during the winter, to agricultural sources in the spring and summer. A quantitative assessment of the significant NRB land use activities indicated that high (greater than 70 percent impervious) and medium (greater than 35 percent impervious) density urban development were the greatest contributors of NPS‐N on a unit area basis (1.9 and 1.6 kg/ha/yr, respectively), followed by row crops and pasture/hay cover types (1.4 kg/ha/yr).     

URBAN INFLUENCE ON PHOSPHORUS AND SEDIMENT LOADING OF WEST POINT LAKE,GEORGIA1

ABSTRACT: West Point Lake is a 10,360 ha mainstream impoundment of the Chattahoochee River located 95 kilometers downstream of Atlanta, Georgia. Origins and magnitude of external total phosphorus (TP) and total suspended solids (TSS) loads from the West Point Lake basin were estimated over a one-year period. Partitioning the drainage basin allowed the sources of these loads to be determined. The upper subbasin area, from Franklin, Georgia, to the headwaters of the Chattahoochee River, contributed 96 percent of the discharge and 97 percent of the TP and TSS loads into West Point Lake. The lower subbasin area, from Franklin to West Point Lake dam, only contributed 3 percent of the TP and TSS loads. Ninety-one percent and 87 percent of the TP and TSS loads, respectively, from the upper subbasin originated from the Atlanta area. Point sources discharged 70 percent and 3 percent of the upper subbasin TP and TSS loads, respectively. A large portion (66 percent) of the TP from the upper subbasin was in the bioavailable form.