LOADING FUNCTIONS FOR PREDICTING NUTRIENT LOSSES FROM COMPLEX WATERSHEDS1

ABSTRACT: A loading function methodology is presented for predicting runoff, sediment, and nutrient losses from complex watersheds. Separate models are defined for cropland, forest, urban and barnyard sources, and procedures for estimating baseflow nutrients are provided. The loading functions are designed for use as a preliminary screening tool to isolate the major contributors in a watershed. Input data sources are readily available and the functions do not require costly calibrations. Data requirements include watershed land use and soil information, daily precipitation and temperature records and rainfall erosivities. Comparison of predicted and measured water, sediment, and nutrient runoff fluxes for the West Branch Deleware River in New York, indicated that runoff was underpredicted by about 14 percent while dissolved nutrients were within 30 percent of observed values. Sediment and solid-phase nutrients were overpredicted by about 50 percent. An annual nutrient budget for the West Branch Delaware River showed that cornland was the major source of sediment, solid phase nutrients, and total phosphorus. Waste water treatment plants and ground water discharge contributed the most dissolved phosphorus and dissolved nitrogen, respectively.     

Phosphorus loss and runoff characteristics in three adjacent agricultural watersheds with claypan soils

Effects of precipitation, runoff, and management on total phosphorus (TP) loss from three adjacent, row-cropped watersheds in the claypan region of northeastern Missouri were examined from 1991 to 1997 to understand factors affecting P loss in watersheds dominated by claypan soils. Runoff samples from each individual runoff event were analyzed for TP and sediment concentration. The annual TP loss ranged from 0.29 to 3.59 kg ha(-1) with a mean of 1.36 kg ha(-1) across all the watersheds during the study period. Significantly higher loss of TP from the watersheds was observed during the fallow period. Multiple small runoff events or several large runoff events contributed to loss of TP from the watersheds. Total P loss in 1993, a year with above-normal precipitation, accounted for 30% of the total TP loss observed over seven years. The five largest runoff events out of a total of 66 events observed over seven years accounted for 27% of the TP loss. The five largest sediment losses were responsible for 24% of the TP loss over seven years. Runoff volume and sediment loss explained 64 to 73% and 47 to 58% of the variation in TP loss on watersheds during the study. Flow duration and maximum flow accounted for 49 and 66% of TP loss, respectively. The results of this study suggest that management practices that reduce runoff volume, flow duration, maximum flow, and sediment loss, and that maintain a suitable vegetative cover throughout the year could lower P loss in claypan soils.     

Nonpoint Source Loading Rates From Selected Land Uses1

ABSTRACT: Loading rates derived from monitoring natural runoff from selected land uses are compared. Land uses selected for evaluation are construction sites, barnyards, and agriculture (dairying). Runoff volumes, sediment, and nutrient fractions were monitored and expressed as areal loadings for comparison purposes. Sediment yield and total phosphorus (total P) loss was directly proportional to runoff (m³/ha). In decreasing order, the loadings for sediment and total P were as follows: construction site > barnyard > general dairying. Runoff from the barnyard area was approximately 10 times higher in soluble phosphorus and ammonium nitrogen than the other land uses under investigation. Areal loss for nitrate nitrogen was highest from the construction site and was attributed to the higher volume of runoff per unit area. Results show that barnyards in a dairying watershed are potentially a major source of sediment and nutrients, especially those dissolved fractions which have the potential for immediate water quality impacts. Relative to general agricultural land, urban construction sites also appear to be a major source of sediment and nutrients. As with barnyard sites, however, the effect of such sites on water quality likely depends on proximity to surface water bodies and other watershed characteristics affecting delivery ratios of contaminants.     

Evaluating agricultural best management practices in tile-drained subwatersheds of the Mackinaw River, Illinois

Best management practices (BMPs) are widely promoted in agricultural watersheds as a means of improving water quality and ameliorating altered hydrology. We used a paired watershed approach to evaluate whether focused outreach could increase BMP implementation rates and whether BMPs could induce watershed-scale (4000 ha) changes in nutrients, suspended sediment concentrations, or hydrology in an agricultural watershed in central Illinois. Land use was >90% row crop agriculture with extensive subsurface tile drainage. Outreach successfully increased BMP implementation rates for grassed waterways, stream buffers, and strip-tillage within the treatment watershed, which are designed to reduce surface runoff and soil erosion. No significant changes in nitrate-nitrogen (NO-N), total phosphorus (TP), dissolved reactive phosphorus, total suspended sediment (TSS), or hydrology were observed after implementation of these BMPs over 7 yr of monitoring. Annual NO-N export (39-299 Mg) in the two watersheds was equally exported during baseflow and stormflow. Mean annual TP export was similar between the watersheds (3.8 Mg) and was greater for TSS in the treatment (1626 ± 497 Mg) than in the reference (940 ± 327 Mg) watershed. Export of TP and TSS was primarily due to stormflow (>85%). Results suggest that the BMPs established during this study were not adequate to override nutrient export from subsurface drainage tiles. Conservation planning in tile-drained agricultural watersheds will require a combination of surface-water BMPs and conservation practices that intercept and retain subsurface agricultural runoff. Our study emphasizes the need to measure conservation outcomes and not just implementation rates of conservation practices.     

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.     

Evaluation of phosphorus transport in surface runoff from packed soil boxes

Evaluation of phosphorus (P) management strategies to protect water quality has largely relied on research using simulated rainfall to generate runoff from either field plots or shallow boxes packed with soil. Runoff from unmanured, grassed field plots (1 m wide x 2 m long, 3-8% slope) and bare soil boxes (0.2 m wide and 1 m long, 3% slope) was compared using rainfall simulation (75 mm h(-1)) standardized by 30-min runoff duration (rainfall averaged 55 mm for field plots and 41 mm for packed boxes). Packed boxes had lower infiltration (1.2 cm) and greater runoff (2.9 cm) and erosion (542 kg ha(-1)) than field plots (3.7 cm infiltration; 1.8 cm runoff; 149 kg ha(-1) erosion), yielding greater total phosphorus (TP) losses in runoff. Despite these differences, regressions of dissolved reactive phosphorus (DRP) in runoff and Mehlich-3 soil P were consistent between field plots and packed boxes reflecting similar buffering by soils and sediments. A second experiment compared manured boxes of 5- and 25-cm depths to determine if variable hydrology based on box depth influenced P transport. Runoff properties did not differ significantly between box depths before or after broadcasting dairy, poultry, or swine manure (100 kg TP ha(-1)). Water-extractable phosphorus (WEP) from manures dominated runoff P, and translocation of manure P into soil was consistent between box types. This study reveals the practical, but limited, comparability of field plot and soil box data, highlighting soil and sediment buffering in unamended soils and manure WEP in amended soils as dominant controls of DRP transport.     

Agroforestry buffers for nonpoint source pollution reductions from agricultural watersheds

Despite increased attention and demand for the adoption of agroforestry practices throughout the world, rigorous long-term scientific studies confirming environmental benefits from the use of agroforestry practices are limited. The objective was to examine nonpoint-source pollution (NPSP) reduction as influenced by agroforestry buffers in watersheds under grazing and row crop management. The grazing study consists of six watersheds in the Central Mississippi Valley wooded slopes and the row crop study site consists of three watersheds in a paired watershed design in Central Claypan areas. Runoff water samples were analyzed for sediment, total nitrogen (TN), and total phosphorus (TP) for the 2004 to 2008 period. Results indicate that agroforestry and grass buffers on grazed and row crop management sites significantly reduce runoff, sediment, TN, and TP losses to streams. Buffers in association with grazing and row crop management reduced runoff by 49 and 19%, respectively, during the study period as compared with respective control treatments. Average sediment loss for grazing and row crop management systems was 13.8 and 17.9 kg ha yr, respectively. On average, grass and agroforestry buffers reduced sediment, TN, and TP losses by 32, 42, and 46% compared with the control treatments. Buffers were more effective in the grazing management practice than row crop management practice. These differences could in part be attributed to the differences in soils, management, and landscape features. Results from this study strongly indicate that agroforestry and grass buffers can be designed to improve water quality while minimizing the amount of land taken out of production.     

Rate of fall-applied liquid swine manure: effects on runoff transport of sediment and phosphorus

Reducing the delivery of phosphorus (P) from land-applied manure to surface water is a priority in many watersheds. Manure application rate can be controlled to manage the risk of water quality degradation. The objective of this study was to evaluate how application rate of liquid swine manure affects the transport of sediment and P in runoff. Liquid swine manure was land-applied and incorporated annually in the fall to runoff plots near Morris, Minnesota. Manure application rates were 0, 0.5, 1, and 2 times the rate recommended to supply P for a corn (Zea mays L.)-soybean [Glycine max (L.) Merr.] rotation. Runoff volume, sediment, and P transport from snowmelt and rainfall were monitored for 3 yr. When manure was applied at the highest rate, runoff volume and sediment loss were less than the control plots without manure. Reductions in runoff volume and soil loss were not observed for spring runoff when frozen soil conditions controlled infiltration rates. The reduced runoff and sediment loss from manure amended soils compensated for addition of P, resulting in similar runoff losses of total P among manure application rates. However, losses of dissolved P increased with increasing manure application rate for runoff during the spring thaw period. Evaluation of water quality risks from fall-applied manure should contrast the potential P losses in snowmelt runoff with the potential that incorporated manure may reduce runoff and soil loss during the summer.     

SAGEBRUSH RANGELAND HYDROLOGY AND EVALUATION OF THE SPUR HYDROLOGY MODEL1

ABSTRACT: An excellent hydrologic record on sagebrush range-land has been developed at the Reynolds Creek Experimental Watershed in southwestern Idaho. The objectives of this paper were two-fold: (1) to analyze and describe the hydrologic record (8–18 years) from four sagebrush watersheds (1–83 ha); and (2) to evaluate the hydrology component of SPUR, a comprehensive rangeland model. The watersheds represent a gradient in elevation (1180–1658 m) and precipitation (240–350 mm/yr). Runoff was a small fraction (> 2 percent) of the total water budget for all of the watersheds. It occurred very infrequently at the three lower elevation watersheds: Summit, Flats, and Nancy Gulch. At Lower Sheep, the highest elevation watershed, runoff occurred most years for a period of 1 to 17 weeks in the winter. Frozen soil combined with rainfall or snowmelt was associated with most of the runoff from Flats and Nancy Gulch. At Summit summertime thunderstorms produced all of the runoff. The average annual sediment yield from all of the watersheds was low (17–950 kg/ha). It was highest from Summit, which had well developed alluvial channels and very steep slopes. SPUR was able to simulate runoff with reasonable accuracy only at Summit, where frozen soils were not a factor. There was poor correlation between predicted and actual annual 8ediment loss. The model tended to overpredict evapotranspiration early in the growing season and underpredict it in the late summer.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

ABSTRACT: Significant land cover changes have occurred in the watersheds that contribute runoff to the upper San Pedro River in Sonora, Mexico, and southeast Arizona. These changes, observed using a series of remotely sensed images taken in the 1970s, 1980s, and 1990s, have been implicated in the alteration of the basin hydrologic response. The Cannonsville subwatershed, located in the Catskill/Delaware watershed complex that delivers water to New York City, provides a contrast in land cover change. In this region, the Cannonsville watershed condition has improved over a comparable time period. A landscape assessment tool using a geographic information system (GIS) has been developed that automates the parameterization of the Soil and Water Assessment Tool (SWAT) and KINEmatic Runoff and EROSion (KINEROS) hydrologic models. The Automated Geospatial Watershed Assessment (AGWA) tool was used to prepare parameter input files for the Upper San Pedro Basin, a subwatershed within the San Pedro undergoing significant changes, and the Cannonsville watershed using historical land cover data. Runoff and sediment yield were simulated using these models. In the Cannonsville watershed, land cover change had a beneficial impact on modeled watershed response due to the transition from agriculture to forest land cover. Simulation results for the San Pedro indicate that increasing urban and agricultural areas and the simultaneous invasion of woody plants and decline of grasslands resulted in increased annual and event runoff volumes, flashier flood response, and decreased water quality due to sediment loading. These results demonstrate the usefulness of integrating remote sensing and distributed hydrologic models through the use of GIS for assessing watershed condition and the relative impacts of land cover transitions on hydrologic response.     

Stormwater runoff and export changes with development in a traditional and low impact subdivision

Development continues at a rapid pace throughout the country. Runoff from the impervious surfaces in these watersheds continues to be a major cause of degradation to freshwater bodies and estuaries. Low impact development techniques have been recommended to reduce these impacts. In this study, stormwater runoff and pollutant concentrations were measured as development progressed in both a traditional development, and a development that used low impact development techniques. Increases in total impervious area in each watershed were also measured. Regression relationships were developed between total impervious area and stormwater runoff/pollutant export. Significant, logarithmic increases in stormwater runoff and nitrogen and phosphorus export were found as development occurred in the traditional subdivision. The increases in stormwater runoff and pollutant export were more than two orders of magnitude. TN and TP export after development was 10 and 1 kg ha(-1) yr(-1), respectively, which was consistent with export from other urban/developed areas. In contrast, stormwater runoff and pollutant export from the low impact subdivision remained unchanged from pre-development levels. TN and TP export from the low impact subdivision were consistent with export values from forested watersheds. The results of this study indicate that the use of low impact development techniques on a watershed scale can greatly reduce the impacts of development on local waterways.     

Calibration and Verification of SWMM for Low Impact Development

The Storm Water Management Model was used to simulate runoff and nutrient export from a low impact development (LID) watershed and a watershed using traditional runoff controls. Predictions were compared to observed values. Uncalibrated simulations underpredicted weekly runoff volume and average peak flow rates from the multiple subcatchment LID watershed by over 80%; the single subcatchment traditional watershed had better predictions. Saturated hydraulic conductivity, Manning's n for swales, and initial soil moisture deficit were sensitive parameters. After calibration, prediction of total weekly runoff volume for the LID and traditional watersheds improved to within 12 and 5% of observed values, respectively. For the validation period, predicted total weekly runoff volumes for the LID and traditional watersheds were within 6 and 2% of observed values, respectively. Water quality simulation was less successful, Nash–Sutcliffe coefficients >0.5 for both calibration and validation periods were only achieved for prediction of total nitrogen export from the LID watershed. Simulation of a 100‐year, 24‐h storm resulted in a runoff coefficient of 0.46 for the LID watershed and 0.59 for the traditional watershed. Results suggest either calibration is needed to improve predictions for LID watersheds or expanded look‐up tables for Green–Ampt infiltration parameter values that account for compaction of urban soil and antecedent conditions are needed.     

Manure composition and incorporation effects on phosphorus in runoff following corn biomass removal

Greater demand for corn ( L.) stover for bioenergy use may lead to increased corn production acreage with minimal surface residue cover, resulting in greater risk for soil erosion and phosphorus (P) losses in runoff. A rainfall simulation study was conducted to determine the effects of spring-applied dairy cow () manure (none, in-barn composted, and exterior walled-enclosure pit) with >200 g kg organic solids content following fall corn biomass removal with and without incorporation (chisel plow [CP] and no-till [NT]) on sediment and P in runoff. Runoff was collected from a 0.83-m area for 60 min following the onset of rainfall simulation (76 mm h), once in spring and once in fall. Runoff dissolved reactive P (DRP) and dissolved organic P (DOP) concentrations were positively correlated with manure P rate and were higher in NT compared with CP. Conversely, sediment and particulate P (PP) concentrations in runoff were inversely correlated with manure P rate (and manure solids) and were higher in CP compared with NT. Runoff volume where no manure was applied was higher in NT than in CP in spring but similar in fall. The addition of manure reduced runoff volumes by an average of 82% in NT and 42% in CP over spring and fall. Results from this study indicate that surface application of dairy manure with relatively high solids content may reduce sediment and PP losses in runoff without increasing the risk of increased DRP and DOP losses in the year of application where corn biomass is harvested.     

Modeling phosphorus transport in an agricultural watershed using the WEPP model

The Water Erosion Prediction Project (WEPP) model has been tested for its ability to predict soil erosion, runoff, and sediment delivery over a wide range of conditions and scales for both hillslopes and watersheds. Since its release in 1995, there has been considerable interest in adding a chemical transport element to it. Total phosphorus (TP) loss at the watershed outlet was simulated as the product of TP in the soil, amount of sediment at the watershed outlet, and an enrichment ratio (ER) factor. WEPP can be coupled with a simple algorithm to simulate phosphorus transport bound to sediment at the watershed outlet. The objective of this work was to incorporate and test the ability of WEPP in estimatingTP loss with sediment at the small watershed scale. Two approaches were examined. One approach (P-EER) estimated ER according to an empirical relationship; the other approach used the ER calculated by WEPP (P-WER).The data used for model performance test were obtained from two side-by-side watersheds monitored between 1976 and 1980. The watershed sizes were 5.05 and 6.37 ha, and each was in a corn (Zea mays L.)-soybean [Glycine max (L.) Merr.] rotation. Measured and simulated results were compared for the period April to October in each year. There was no statistical difference between the mean measured and simulated TP loss. The Nash-Sutcliffe coefficient was 0.80 and 0.78 for the P-EER and P-WER methods, respectively. It was critical for both methods that WEPP adequately represent the biggest sediment yield events because sediment is the main driver for TP loss so that the model can adequately simulate TP losses bound to sediment. The P-WER method is recommended because it does not require use of empirical parameters to estimate TP loss at the watershed outlet.     

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.     

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.     

EFFECTS OF PRECIPITATION AND LAND USE ON STORM RUNOFF1

ABSTRACT: Storm-runoff quantity and quality were studied in three watersheds located near St. Paul in Ramsey County, Minnesota, from April 15 through September 15 of 1984, 1985, and 1986 to qualitatively determine the effects of precipitation and selected land uses on storm runoff. In respect to precipitation effects, differences in stormrunoff quantity between years in an urban watershed that lacks wetlands appear to be related to the average storm size (amount of precipitation) during the study period of each year. In contrast, the differences in storm-runoff quantity from watersheds that contain wetlands appear to be related to total precipitation during study period of each year. In respect to land use, the differences in storm-runoff quantity appear to be related to the amounts of impervious and wetland area. The watershed that contains the largest amount of impervious area and smallest amount of wetland area has the largest amount of storm runoff. Differences in storm-runoff quality appear to be related to the amounts of wetland and lake area. The watershed that contains the largest amounts of wetland and lake area has the smallest storm-runoff loading of suspended solids, phosphorus, and nitrogen. The wetland and lake areas likely retain the loading and, subsequently, lower the amount of storm-runoff loading exported from a watershed.     

Runoff phosphorus losses from surface-applied biosolids

Runoff losses of dissolved and particulate phosphorus (P) may occur when rainfall interacts with manures and biosolids spread on the soil surface. This study compared P levels in runoff losses from soils amended with several P sources, including 10 different biosolids and dairy manure (untreated and treated with Fe or Al salts). Simulated rainfall (71 mm h(-1)) was applied until 30 min of runoff was collected from soil boxes (100 x 20 x 5 cm) to which the P sources were surfaced applied. Materials were applied to achieve a common plant available nitrogen (PAN) rate of 134 kg PAN ha(-1), resulting in total P loading rates from 122 (dairy manure) to 555 (Syracuse N-Viro biosolids) kg P ha(-1). Two biosolids produced via biological phosphorus removal (BPR) wastewater treatment resulted in the highest total dissolved phosphorus (13-21.5 mg TDP L(-1)) and total phosphorus (18-27.5 mg TP L(-1)) concentrations in runoff, followed by untreated dairy manure that had statistically (p = 0.05) higher TDP (8.5 mg L(-1)) and TP (10.9 mg L(-1)) than seven of the eight other biosolids. The TDP and TP in runoff from six biosolids did not differ significantly from unamended control (0.03 mg TDP L(-1); 0.95 mg TP L(-1)). Highest runoff TDP was associated with P sources low in Al and Fe. Amending dairy manure with Al and Fe salts at 1:1 metal-to-P molar ratio reduced runoff TP to control levels. Runoff TDP and TP were not positively correlated to TP application rate unless modified by a weighting factor reflecting the relative solubility of the P source. This suggests site assessment indices should account for the differential solubility of the applied P source to accurately predict the risk of P loss from the wide variety of biosolids materials routinely land applied.     

EVALUATION OF INFLOW TO MIRROR LAKE,NEW HAMPSHIRE1

ABSTRACT: Measured stream discharge plus calculated ground water discharge (total measured runoff) were compared with runoff calculated by the unit-runoff method for the two largest watersheds of Mirror Lake for 1981–1983. Runoff calculated by the unit-runoff method, using Hubbard Brook watershed 3 as the index watershed, was greater than the total measured runoff into Mirror Lake during periods of high flow and slightly less than the total measured runoff into Mirror Lake during periods of low flow. Annual calculated unit runoff was 17 to 37 percent greater than total measured runoff. Differences in monthly runoff are far greater, ranging from 0 to greater than 100 percent. For high flows the calculated unit runoff is about 2 times greater than total measured runoff. For low flows the northwest basin of Mirror Lake has the greatest ground water contribution compared to the other two basins. In contrast, Hubbard Brook watershed 3 has the least ground water contribution.     

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.