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 THE HYDROCHEMISTRY OF THE CANNONSVILLE WATERSHED WITH GENERALIZED WATERSHED LOADING FUNCTIONS (GWLF)1

ABSTRACT: A previous modeling study used the Generalized Watershed Loading Functions (GWLF) model to simulate stream‐flow, and nutrient and sediment loads to Cannonsville Reservoir from the West Branch Delaware River (WBDR). We made several model revisions, calibrated key parameters, and tested the original GWLF model and a revised GWLF model using more recent data. Model revisions included: addition of unsaturated leakage between unsaturated and saturated subsurface reservoirs; revised timing of sediment export; inclusion of urban sediments and dissolved nutrients; tracking of particulate nutrients from point sources; and revised timing of septic system loads. The revision of sediment yield timing resulted in significant improvements in monthly sediment and particulate phosphorus predictions as compared to the original model. Addition of unsaturated leakage improved hydrologic predictions during low flow months. The other model changes improve realism without adding significant model complexity or data requirements. Goodness of fit of revised model predictions versus stream measurements, as measured by the Nash‐Sutcliff coefficient of model efficiency, exceeded 0.8 for streamflow‐0.7 for sediment yield and dissolved nitrogen (N) and 0.6 for particulate and dissolved phosphorus (P). The revised GWLF model, with limited calibration, provides reasonable estimates of monthly streamflow, and nutrient and sediment loads in the Cannonsville watershed.     

EVALUATION OF RUNOFF,EROSION, AND PHOSPHORUS MODELING SYSTEM—SIMPLE1

ABSTRACT: The purpose of this study was to evaluate the performance of Spatially Integrated Models for Phosphorus Loading and Erosion (SIMPLE) in predicting runoff volume, sediment loss, and phosphorus loading from two watersheds. The modeling system was applied to the 334 ha QOD subwatershed, part of the Owl Run watershed, located in Fauquier County, Virginia, and to the 2240 ha watershed, Battle Branch, located in Delaware County, Oklahoma. Simulation runs were conducted at cell and field scales, and simulation results were compared with observed data. Runoff volume and dissolved phosphorus loading were measured at the Battle Branch watershed. Runoff volume, sediment yield, and total phosphorus loading were measured at the QOD site. SIMPLE tended to underestimate runoff volumes during the dormant period, from November to March. The comparison between observed and predicted dissolved phosphorus showed better correlation than for observed and predicted total phosphorus loading. Cell level simulations provided similar estimates of runoff volume and phosphorus loading when compared to field level simulations for both watersheds. However, observed sediment yields better compared with the values predicted from the cell level simulation when compared to field level simulation. Finally, results of model evaluation indicated that SIMPLE's predictive ability is acceptable for screening applications but not for site-specific quantitative predictions.     

Integrated Modular Modeling of Water and Nutrients From Point and Nonpoint Sources in the Patuxent River Watershed1

Abstract: We present a simple modular landscape simulation model that is based on a watershed modeling framework in which different sets of processes occurring in a watershed can be simulated separately with different models. The model consists of three loosely coupled submodels: a rainfall‐runoff model (TOPMODEL) for runoff generation in a subwatershed, a nutrient model for estimation of nutrients from nonpoint sources in a subwatershed, and a stream network model for integration of point and nonpoint sources in the routing process. The model performance was evaluated using monitoring data in the watershed of the Patuxent River, a tributary to the Chesapeake Bay in Maryland, from July 1997 through August 1999. Despite its simplicity, the landscape model predictions of streamflow, and sediment and nutrient loads were as good as or better than those of the Hydrological Simulation Program‐Fortran model, one of the most widely used comprehensive watershed models. The landscape model was applied to predict discharges of water, sediment, silicate, organic carbon, nitrate, ammonium, organic nitrogen, total nitrogen, organic phosphorus, phosphate, and total phosphorus from the Patuxent watershed to its estuary. The predicted annual water discharge to the estuary was very close to the measured annual total in terms of percent errors for both years of the study period (≤2%). The model predictions for loads of nutrients were also good (20‐30%) or very good (<20%) with exceptions of sediment (40%), phosphate (36%), and organic carbon (53%) for Year 1.     

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.     

Relationship of Water Quality and Pollutant Loads to Land Uses in adjoining watersheds1

ABSTRACT: The Grand and Saugeen Rivers in southern Ontario were chosen for study as pilot watersheds under the Pollution From Land Use Activities Reference Group (PLUARG) study. The pilot watersheds have adjacent headwater areas and are physically similar in geology, physiography, and climate. Significant differences in water quality between the watersheds at their outlets are attributed to land use and population differences. The major pollutant sources in the two pilot watersheds were identified as trace elements from urban runoff and point source discharges; phosphorus from agricultural and urban runoff and private waste disposal; chloride from transportation corridors; and sediment and nitrogen from agricultural runoff. Yields at the watershed outlets were similar for suspended sediment and two to three times as high in the Grand River for phosphours, nitrogen, chloride, and lead. The higher phosphorus and nitrogen levels were attributed to larger point source inputs and the higher proportion of agricultural activity, comprising 75 percent of the Grand River basin compared to 64 percent in the Saugeen River basin. Similarly, the higher chloride and lead levels were attributed to an order of magnitude larger population and three times as much urban land in the Grand River basin compared to the Saugeen River basin.     

WATERSHED WEIGHTING OF EXPORT COEFFICIENTS TO MAP CRITICAL PHOSPHOROUS LOADING AREAS1

ABSTRACT: The Export Coefficient model (ECM) is capable of generating reasonable estimates of annual phosphorous loading simply from a watershed's land cover data and export coefficient values (ECVs). In its current form, the ECM assumes that ECVs are homogeneous within each land cover type, yet basic nutrient runoff and hydrological theory suggests that runoff rates have spatial patterns controlled by loading and filtering along the flow paths from the upslope contributing area and downslope dispersal area. Using a geographic information system (GIS) raster, or pixel, modeling format, these contributing area and dispersal area (CADA) controls were derived from the perspective of each individual watershed pixel to weight the otherwise homogeneous ECVs for phosphorous. Although the CADA‐ECM predicts export coefficient spatial variation for a single land use type, the lumped basin load is unaffected by weighting. After CADA weighting, a map of the new ECVs addressed the three fundamental criteria for targeting critical pollutant loading areas: (1) the presence of the pollutant, (2) the likelihood for runoff to carry the pollutant offsite, and (3) the likelihood that buffers will trap nutrients prior to their runoff into the receiving water body. These spatially distributed maps of the most important pollutant management areas were used within New York's West Branch Delaware River watershed to demonstrate how the CADA‐ECM could be applied in targeting phosphorous critical loading areas.     

Quantitative Assessment of Agricultural Runoff and Soil Erosion Using Mathematical Modeling: Applications in the Mediterranean Region

Three mathematical models, the runoff curve number equation, the universal soil loss equation, and the mass response functions, were evaluated for predicting nonpoint source nutrient loading from agricultural watersheds of the Mediterranean region. These methodologies were applied to a catchment, the gulf of Gera Basin, that is a typical terrestrial ecosystem of the islands of the Aegean archipelago. The calibration of the model parameters was based on data from experimental plots from which edge-of-field losses of sediment, water runoff, and nutrients were measured. Special emphasis was given to the transport of dissolved and solid-phase nutrients from their sources in the farmers' fields to the outlet of the watershed in order to estimate respective attenuation rates. It was found that nonpoint nutrient loading due to surface losses was high during winter, the contribution being between 50% and 80% of the total annual nutrient losses from the terrestrial ecosystem. The good fit between simulated and experimental data supports the view that these modeling procedures should be considered as reliable and effective methodological tools in Mediterranean areas for evaluating potential control measures, such as management practices for soil and water conservation and changes in land uses, aimed at diminishing soil loss and nutrient delivery to surface waters. Furthermore, the modifications of the general mathematical formulations and the experimental values of the model parameters provided by the study can be used in further application of these methodologies in watersheds with similar characteristics.     

TRANSPORT AND DISTRIBUTION OF NUTRIENTS IN THE LOXAHATCHEE RIVER ESTUARY,SOUTHEASTERN FLORIDA, 1979–811

ABSTRACT: Concentrations of total nitrogen, total phosphorus, and total organic carbon in the Loxahatchee River estuary decreased with increasing salinity in a manner indicating that mixing and dilution of freshwater by seawater was the primary process controlling the down-stream concentrations of nutrients. Most of the nutrients in the surface freshwater inflows entered the estuary from five major tributaries; however, about 10 percent of the total nitrogen and 32 percent of the total phosphorus were from urban stormwater runoff. The input of nutrients was highly seasonal and storm related. During a 61-day period of above average rainfall that included Tropical Storm Dennis, the major tributaries discharged 2.7 metric tons of total phosphorus, 75 metric tons of total nitrogen, and 1,000 metric tons of organic carbon to the estuary. This period accounted for more than half of the total nutrient load from the major tributaries during the 1981 water year (October 1, 1980, through September 30, 1981). Inorganic phosphorus and nitrogen increased relative to total phosphorus and nitrogen during storm runoff. Nutrient yield from the basin, expressed as grams per square meter of basin area, was relatively low. However, because the basin area (544 square kilometers) is large compared with the volume of the estuary, the basin might be expected to contribute significantly to estuarine enrichment were it not for tidal flushing. Approximately 60 percent of the total volume of the estuary is flushed on each tide. Because the estuary is well flushed, it probably has a large tolerance for nutrient loading.     

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.     

Development of the Chesapeake Bay Watershed Total Maximum Daily Load Allocation

Nutrient load allocations and subsequent reductions in total nitrogen and phosphorus have been applied in the Chesapeake watershed since 1992 to reduce hypoxia and to restore living resources. In 2010, sediment allocations were established to augment nutrient allocations supporting the submerged aquatic vegetation resource. From the initial introduction of nutrient allocations in 1992 to the present, the allocations have become more completely applied to all areas and loads in the watershed and have also become more rigorously assessed and tracked. The latest 2010 application of nutrient and sediment allocations were made as part of the Chesapeake Bay total maximum daily load and covered all six states of the Chesapeake watershed. A quantitative allocation process was developed that applied principles of equity and efficiency in the watershed, while achieving all tidal water quality standards through an assessment of equitable levels of effort in reducing nutrients and sediments. The level of effort was determined through application of two key watershed scenarios: one where no action was taken in nutrient control and one where maximum nutrient control efforts were applied. Once the level of effort was determined for different jurisdictions, the overall load reduction was set watershed‐wide to achieve dissolved oxygen water quality standards. Further adjustments were made to the allocation to achieve the James River chlorophyll‐a standard.     

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.     

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.     

INTEGRATION OF GEOGRAPHIC INFORMATION SYSTEMS AND A COMPUTER MODEL TO EVALUATE IMPACTS OF AGRICULTURAL RUNOFF ON WATER QUALITY1

ABSTRACT: This study integrates an Agricultural Non-Point Source Pollution Model (AGNPS), the Geographic Resource Analysis Support System (GRASS) (U.S. Army Corps of Engineers, 1987), and GRASS WATERWORKS (a hydrologic modeling tool box being developed at the Michigan State University Center for Remote Sensing) to evaluate the impact of agricultural runoff on water quality in the Cass River, a subwatershed of Saginaw Bay. AGNPS is used to estimate the amounts, origin, and distribution of sediment, nitrogen (N), and phosphorus (P) in the watershed. GRASS and GRASS WATERWORKS are used to generate parameters needed for AGNPS from digital maps, which include soil association, land use, watershed boundaries, water features, and digital elevation. Outputs of the model include spatially distributed estimates of volume and peak runoff, overland and channel erosion, sediment yields, and concentrations of nitrogen and phosphorus. Management scenarios are explored in the AGNPS model to minimize sedimentation and nutrient loading. Scenarios evaluated include variations in crop cover, tillage methods, and other agricultural management practices. In addition, areas vulnerable to erosion are identified for best management practices.     

Eutrophication of Buttermilk Bay,a cape cod coastal embayment: Concentrations of nutrients and watershed nutrient budgets

Nutrient concentrations in Buttermilk Bay, a coastal embayment on the northern end of Buzzards Bay, MA, are higher in the nearshore where salinities are lower. This pattern suggests that freshwater sources may contribute significantly to nutrient inputs into Buttermilk Bay. To evaluate the relative importance of the various sources we estimated inputs of nutrients by each major source into the watershed and into the bay itself. Septic systems contributed about 40% of the nitrogen and phosphorus entering the watershed, with precipitation and fertilizer use adding the remainder. Groundwater transported over 85% of the nitrogen and 75% of the phosphorus entering the bay. Most nutrients entering the watershed failed to reach the bay; uptake by forests, soils, denitrification, and adsorption intercepted two-thirds of the nitrogen and nine-tenths of the phosphorus that entered the watershed. The nutrients that did reach the bay most likely originated from subsoil injections into groundwater by septic tanks, plus some leaching of fertilizers.Buttermilk Bay water has relatively low nutrient concentrations, probably because of uptake of nutrients by macrophytes and because of relatively rapid tidal flushing. Annual budgets of nutrients entering the watershed showed a low nitrogen-to-phosphorus ratio of 6, but passage of nutrients through the watershed raised N/P to 23, probably because of adsorption of PO₄ during transit. The N/P ratio of water that leaves the watershed and presumably enters the bay is probably high enough to maintain active growth of nitrogenlimited coastal producers. There is a seasonal shift in N/P in the water column of Buttermilk Bay. N/P exceeded the 161 Redfield ratio during midwinter; the remainder of the year N/P fell below 161. This suggests that annual budgets do not provide sufficiently detailed data with which to interpret nutrient-limitation of producers. Further, some idea of water turnover is also needed to evaluate impact of loading rates. Urbanization of watersheds seems to increase loadings to nearshore environments, and to shift the nutrient loadings delivered to coastal waters to relatively high N-to-P ratios, potentially stimulating growth of nitrogen-limited primary producers.     

NUTRIENT LOADING ASSESSMENT IN THE ILLINOIS RIVER USING A SYNTHETIC APPROACH1

ABSTRACT: A synthetic relationship is developed between nutrient concentrations and discharge rates at two river gauging sites in the Illinois River Basin. Analysis is performed on data collected by the U.S. Geological Survey (USGS) on nutrients in 1990 through 1997 and 1999 and on discharge rates in 1988 through 1997 and 1999. The Illinois River Basin is in western Arkansas and northeastern Oklahoma and is designated as an Oklahoma Scenic River. Consistently high nutrient concentrations in the river and receiving water bodies conflict with recreational water use, leading to intense stakeholder debate on how best to manage water quality. Results show that the majority of annual phosphorus (P) loading is transported by direct runoff, with high concentrations transported by high discharge rates and low concentrations by low discharge rates. A synthetic relationship is derived and used to generate daily phosphorus concentrations, laying the foundation for analysis of annual loading and evaluation of alternative management practices. Total nitrogen (N) concentration does not have as clear a relationship with discharge. Using a simple regression relationship, annual P loadings are estimated as having a root mean squared error (RMSE) of 39.8 t/yr and 31.9 t/yr and mean absolute percentage errors of 19 percent and 28 percent at Watts and Tahlequah, respectively. P is the limiting nutrient over the full range of discharges. Given that the majority of P is derived from Arkansas, management practices that control P would have the most benefit if applied on the Arkansas side of the border.     

Watershed scale assessment of nitrogen and phosphorus loadings in the Indian River Lagoon basin,Florida

There is a growing evidence that the ecological and biological integrity of the lagoon has declined during the last 50 years, probably due to the decline in water quality. Establishment of a watershed scale seagrass-based nutrient load assessment is the major aim of water quality management in the Indian River Lagoon (IRL). Best estimate loadings incorporate wet and dry deposition, surface water, groundwater, sediment nutrient flux, and point source effluent discharge data. On the average, the IRL is receiving annual external loadings of 832, 645 and 94,476kg of total nitrogen (TN) and total phosphorus (TP), respectively, from stormwater discharges and agricultural runoff. The average internal cycling of TN and TP from sediment deposits in the IRL was about 42,640kg TN and 1050kg TPyr(-1). Indirect evidence suggests that atmospheric deposition has played a role in the ongoing nutrient enrichment in the IRL. The estimated total atmospheric deposition of TN and TP was about 32,940 and 824kgyr(-1), while groundwater contribution was about 84,920 and 24,275kgyr(-1), respectively, to the surface waters of the IRL. The estimated annual contribution of point effluent discharge was about 60,408kg TN and 7248kg TP. In total, the IRL basin is receiving an annual loading of about 1,053,553kg TN and 127,873kg TP. With these results, it is clear that the current rate of nutrient loadings is causing a shift in the primary producers of the IRL from macrophyte to phytoplankton- or algal-based system. The goal is to reverse that shift, to attain and maintain a macrophyte-based estuarine system in the IRL.     

Intra-storm study of solute chemical composition of overland flow water in two agricultural fields

Few studies have assessed the transport of dissolved nutrients at the field scale under natural rainfall conditions. Hysteresis between dissolved nutrients and discharge behavior can complicate such assessments and this effect has only been examined qualitatively. In this study, we investigated factors contributing to short-term variations of dissolved cation (Ca, Mg, Na, and K) and anion (soluble reactive phosphorus [SRP], NO3, and SO4) concentrations in runoff water and developed a quantitative method to study their hysteretic behavior. Within-storm variations of dissolved nutrient concentrations were determined in two agricultural fields during four natural rainfall events along with discharge, sediment, antecedent soil water conditions, and nutrient contents. For each event, nutrient loads were plotted against discharge during the rising and falling limb of the runoff hydrograph. The resulting hysteresis curves were characterized by an index H, which is the ratio between the integrated areas under the rising and falling curves of the hydrograph. Results showed that nutrient concentrations increased with time during each event. Counterclockwise (H < 1) hysteresis, occurring when the falling limb had larger loads, was found when soils were initially dry whereas clockwise hysteresis (H > 1) was associated with prior wet soil conditions. Two hypotheses are suggested to explain these variations. First, suspended sediments could have acted as a sink for dissolved nutrients and the sensitivity of nutrients to hydrological conditions was determined by their preferential sorption on these sediments. Second, movement of nutrients into runoff occurred more readily as soils became wetter during an event.     

AN EXAMINATION OF LAND USE - NUTRIENT EXPORT RELATIONSHIPS1

ABSTRACT: Lake and watershed management strategies and recent environmental legislation dictate that nonpoht nutrient sources associated with storm water runoff must be assessed. Accordingly, a nutrient flu assessment for phosphorus and nitrogen is conducted through an extensive literature review of nutrient export studies. These studies are reevaluated. The nutrient export coefficients are screened according to sampling design criteria and compiled according to land use. The ecological mechanisms within each land use influencing the magnitude of nutrient flux are also discussed     

WATER QUALITY CHARACTERISTICS OF STORM WATER FROM MAJOR LAND USES IN SOUTH FLORIDA1

ABSTRACT: Starting in 1998, a study was conducted to characterize storm water quality from predominant land use types in a coastal watershed along the south central coast of Florida, namely citrus, pasture, urban, natural wetland, row crop, dairy, and golf courses. Sixty‐three sampling sites were located at strategic points on drainage conveyances for each of seven specific land use areas. Runoff samples were collected following storm events that met defined rainfall criteria for a period of 30 months. Nitrogen (N), phosphorus (P), heavy metals, pesticides, and other water quality parameters were determined, and the results were analyzed to compare and characterize land uses as relative sources for these constituents in runoff. Results showed that runoff from most land use types had low dissolved oxygen concentration and that sediment and nutrient concentrations were closely related to land use, particularly to the amount of fertilizer applied in each land use. Among the eight heavy metals tested, copper was the most frequently detected and was mostly associated with runoff from citrus and golf course land uses. High levels of arsenic were also detected in golf course runoff. The most frequently detected pesticide was simazine from citrus. The information and methodologies presented may facilitate pollution source characterization and ecological restoration efforts.