Development and Application of the 2010 Chesapeake Bay Watershed Total Maximum Daily Load Model

The Phase 5.3 Watershed Model simulates the Chesapeake watershed land use, river flows, and the associated transport and fate of nutrient and sediment loads to the Chesapeake Bay. The Phase 5.3 Model is the most recent of a series of increasingly refined versions of a model that have been operational for more than two decades. The Phase 5.3 Model, in conjunction with models of the Chesapeake airshed and estuary, provides estimates of management actions needed to protect water quality, achieve Chesapeake water quality standards, and restore living resources. The Phase 5.3 Watershed Model tracks nutrient and sediment load estimates of the entire 166,000 km² watershed, including loads from all six watershed states. The creation of software systems, input datasets, and calibration methods were important aspects of the model development process. A community model approach was taken with model development and application, and the model was developed by a broad coalition of model practitioners including environmental engineers, scientists, and environmental managers. Among the users of the Phase 5.3 Model are the Chesapeake watershed states and local governments, consultants, river basin commissions, and universities. Development and application of the model are described, as well as key scenarios ranging from high nutrient and sediment load conditions if no management actions were taken in the watershed, to low load estimates of an all‐forested condition.     

Computing Atmospheric Nutrient Loads to the Chesapeake Bay Watershed and Tidal Waters

Application of integrated Chesapeake Bay models of the airshed, watershed, and estuary support air and water nitrogen controls in the Chesapeake. The models include an airshed model of the Mid‐Atlantic region which tracks the estimated atmospheric deposition loads of nitrogen to the watershed, tidal Bay, and adjacent coastal ocean. The three integrated models allow tracking of the transport and fate of nitrogen air emissions, including deposition in the Chesapeake watershed, the subsequent uptake, transformation, and transport to Bay tidal waters, and their ultimate influence on Chesapeake water quality. This article describes the development of the airshed model, its application to scenarios supporting the Chesapeake Total Maximum Daily Load (TMDL), and key findings from the scenarios. Key findings are that the atmospheric deposition loads are among the largest input loads of nitrogen in the watershed, and that the indirect nitrogen deposition loads to the watershed, which are subsequently delivered to the Bay are larger than the direct loads of atmospheric nitrogen deposition to Chesapeake tidal waters. Atmospheric deposition loads of nitrogen deposited in coastal waters, which are exchanged with the Chesapeake, are also estimated. About half the atmospheric deposition loads of nitrogen originate from outside the Chesapeake watershed. For the first time in a TMDL, the loads of atmospheric nitrogen deposition are an explicit part of the TMDL load reductions.     

Response of Nitrogen Loading to the Chesapeake Bay to Source Reduction and Land Use Change Scenarios: A SPARROW‐Informed Analysis

In response to concerns regarding the health of streams and receiving waters, the United States Environmental Protection Agency established a total maximum daily load for nitrogen in the Chesapeake Bay watershed for which practices must be in place by 2025 resulting in an expected 25% reduction in load from 2009 levels. The response of total nitrogen (TN) loads delivered to the Bay to nine source reduction and land use change scenarios was estimated using a Spatially Referenced Regression on Watershed Attributes model. The largest predicted reduction in TN load delivered to the Bay was associated with a scenario in which the mass of TN as fertilizer applied to agricultural lands was decreased. A 25% decrease in the mass of TN applied as fertilizer resulted in a predicted reduction in TN loading to the Bay of 11.3%, which was 2.5–5 times greater than the reductions predicted by other scenarios. Eliminating fertilizer application to all agricultural land in the watershed resulted in a predicted reduction in TN load to the Bay of 45%. It was estimated that an approximate 25% reduction in TN loading to the Bay could be achieved by eliminating fertilizer applied to the 7% of subwatersheds contributing the greatest fertilizer‐sourced TN loads to the Bay. These results indicate that management strategies aimed at decreasing loading from a small number of subwatersheds may be effective for reducing TN loads to the Bay, and similar analyses are possible in other watersheds.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

Stochastic Watershed Water Quality Simulation for TMDL Development – A Case Study in the Newport Bay Watershed1

Abstract: Systematic consideration of uncertainty in data, model structure, and other factors is generally unaddressed in most Total Maximum Daily Load (TMDL) calculations. Our previous studies developed the Management Objectives Constrained Analysis of Uncertainty (MOCAU) approach as an uncertainty analysis technique specifically for watershed water quality models, based on a synthetic case. In this study, we applied MOCAU to analyze diazinon loading in the Newport Bay watershed (Southern California). The study objectives included (1) demonstrating the value of performing stochastic simulation and uncertainty analysis for TMDL development, using MOCAU as the technique and (2) evaluating the existing diazinon TMDL and generating insights for the development of scientifically sound TMDLs, considering uncertainty. The Watershed Analysis Risk Management Framework model was used as an example of a complex watershed model. The study revealed the importance and feasibility of conducting stochastic watershed water quality simulation for TMDL development. The critical role of management objectives in a systematic uncertainty assessment was well demonstrated. The results of this study are intuitive to TMDL calculation, model structure improvement and sampling strategy design.     

Scale Effects on Spatially Varying Relationships Between Urban Landscape Patterns and Water Quality

Scientific interpretation of the relationships between urban landscape patterns and water quality is important for sustainable urban planning and watershed environmental protection. This study applied the ordinary least squares regression model and the geographically weighted regression model to examine the spatially varying relationships between 12 explanatory variables (including three topographical factors, four land use parameters, and five landscape metrics) and 15 water quality indicators in watersheds of Yundang Lake, Maluan Bay, and Xinglin Bay with varying levels of urbanization in Xiamen City, China. A local and global investigation was carried out at the watershed-level, with 50 and 200 m riparian buffer scales. This study found that topographical features and landscape metrics are the dominant factors of water quality, while land uses are too weak to be considered as a strong influential factor on water quality. Such statistical results may be related with the characteristics of land use compositions in our study area. Water quality variations in the 50 m buffer were dominated by topographical variables. The impact of landscape metrics on water quality gradually strengthen with expanding buffer zones. The strongest relationships are obtained in entire watersheds, rather than in 50 and 200 m buffer zones. Spatially varying relationships and effective buffer zones were verified in this study. Spatially varying relationships between explanatory variables and water quality parameters are more diversified and complex in less urbanized areas than in highly urbanized areas. This study hypothesizes that all these varying relationships may be attributed to the heterogeneity of landscape patterns in different urban regions. Adjustment of landscape patterns in an entire watershed should be the key measure to successfully improving urban lake water quality.     

Monitored and Modeled Correlations of Sediment and Nutrients with Chesapeake Bay Water Clarity

This article analyzes the correlations of the observed and modeled light attenuation coefficient, Kd, with in situ total suspended solids (TSS) and chlorophyll‐a concentrations in Chesapeake Bay (CB) tidal waters, and with sediment and nutrient loads from the Chesapeake watershed. Light attenuation is closely related to in situ TSS and chlorophyll‐a concentrations, however, the strength of the correlation differs among the CB segments. There are distinct differences between saline and tidal fresh segments in the main Bay, but less distinction among saline and tidal fresh segments in the tidal tributaries. The correlation between Kd with sediment and nutrient loads is complicated by the lag times of TSS and the chlorophyll‐a responses to reductions in nutrient and sediment loads from the watershed, and also due to the diverse load sources. Three sets of model sensitivity scenarios were performed with: (1) differential sediment and nutrient loads; (2) selective sediment source types; and (3) geographically isolated inputs. The model results yield similar findings as those based on observed data and provide information regarding the effect of sediment on specific water bodies. Based on the model results a method was developed to determine sediment and nutrient load reductions needed to achieve the water clarity standards of the CB segments.     

Watershed associations in West Virginia: their impact on environmental protection

Grassroots watershed associations have formed as an avenue to facilitate public involvement in protecting watersheds. Growth in the number of watershed associations has created a need to evaluate whether their existence enhances environmental protection by local communities. In this research, environmental protection was measured by the number of watershed protection activities conducted and amount of funding directed towards protecting surface waters by non-profit, volunteer organizations. Regression models were utilized to examine the influence of population demographic and watershed characteristics variables to explain the formation of watershed associations along with watershed protection activities and funding. Regression results showed that the presence of watershed associations had positive impacts on both the number of activities and financial resources acquired. The results confirm that watershed associations enhance the ability of communities to obtain funding and to implement watershed protection actions. Thus, watershed associations serve as an effective institution for promoting environmental protection.     

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.     

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.     

Using Ecological Risk Assessment to Identify the Major Anthropogenic Stressor in the Waquoit Bay Watershed,Cape Cod,Massachusetts

The Waquoit Bay Watershed ecological risk assessment was performed by an interdisciplinary and interagency workgroup. This paper focuses on the steps taken to formulate the analysis plan for this watershed assessment. The workgroup initially conducted a series of meetings with the general public and local and state managers to determine environmental management objectives for the watershed. The workgroup then decided that more information was needed on the impacts of six stressors: nutrient enrichment, physical alteration of habitat, altered freshwater flow, toxic chemicals, pathogens, and fisheries harvesting. Assessment endpoints were selected to establish the link between environmental management objectives, impacts of stressors, and scientifically measurable endpoints. The following assessment endpoints were selected: estuarine eelgrass cover, scallop abundance, finfish diversity and abundance, wetland bird distribution and abundance, piping plover distribution and abundance, tissue contaminant levels, and brook trout distribution and abundance in streams. A conceptual model was developed to show the pathways between human activities, stressors, and ecological effects. The workgroup analyzed comparative risks, by first ranking stressors in terms of their potential risk to biotic resources in the watershed. Then stressors were evaluated by considering the components of stressors (e.g., the stressor chemical pollution included both heavy metals and chlorinated solvents components) in terms of intensity and extensiveness. The workgroup identified nutrient enrichment as the major stressor. Nutrient enrichment comprised both phosphorus enrichment in freshwater ponds and nitrogen enrichment within estuaries. Because phosphorus impacts were being analyzed and mitigated by the Air Force Center for Environmental Excellence, this assessment focused on nitrogen. The process followed to identify the predominant stressor and focus the analyses on nitrogen impacts on eelgrass and scallops will serve as an example of how to increase the use of the findings of a watershed assessment in decision making.     

Hydrologic Landscape Classification to Estimate Bristol Bay,Alaska Watershed Hydrology

Hydrologic landscapes (HLs) have proven to be a useful tool for broad scale assessment and classification of landscapes across the United States as they help organize larger geographical areas into areas of similar hydrologic characteristics. We developed a HL classification for the Bristol Bay watershed of southwest Alaska that incorporates indices of annual climate and seasonality, terrain, geology, and the influences of large lakes and glaciers. A HL classification is particularly useful in this large watershed because of its hydrologic and landscape variability, important salmon fishery, variety of environmental and potential anthropogenic stressors, and lack of widespread hydrologic data. Following creation of Bristol Bay basin‐wide HL classes, we compared the HL distributions within watersheds grouped by two calculated runoff parameters derived from available long‐term streamflow records and found HL distributions within these groups provided predictive insight on hydrologic behavior. Using these developed runoff groups, we estimated expected hydrologic behavior in watersheds across the larger Bristol Bay watershed that lacked gauged streamflow records. The HL approach provides a scientific basis for estimating the first‐order hydrologic behavior of watersheds and landscapes that lack detailed hydrologic information.     

SWATmodel: A Multi‐Operating System,Multi‐Platform SWAT Model Package in R

The Soil and Water Assessment Tool (SWAT) model (Arnold et al., 1998) is a popular watershed management tool. Currently, the SWAT model, actively supported by the U.S. Department of Agriculture and Texas A&M, operates only on Microsoft^® Windows, which hinders modelers that use other operating systems (OS). This technical note introduces the Comprehensive R Archive Network (CRAN) distributed “SWATmodel” package which allows SWAT 2005 and 2012 to be widely distributed and run as a linear model‐like function on multiple OS and processor platforms. This allows researchers anywhere in the world using virtually any OS to run SWAT. In addition to simplifying the use of SWAT across computational platforms, the SWATmodel package allows SWAT modelers to utilize the analytical capabilities, statistical libraries, modeling tools, and programming flexibility inherent to R. The software allows watershed modelers to develop a simple hydrological watershed model conceptualization of the SWAT model and to obtain a first approximation of the minimum expected results a more complicated model should deliver. As a proof of concept, we test the SWAT model by initializing and calibrating 314 U.S. Geological Survey stream gages in the Chesapeake Bay watershed and present the results.     

Using Multiple Watershed Models to Predict Water,Nitrogen, and Phosphorus Discharges to the Patuxent Estuary1

Boomer, Kathleen M.B., Donald E. Weller, Thomas E. Jordan, Lewis Linker, Zhi‐Jun Liu, James Reilly, Gary Shenk, and Alexey A. Voinov, 2012. Using Multiple Watershed Models to Predict Water, Nitrogen, and Phosphorus Discharges to the Patuxent Estuary. Journal of the American Water Resources Association (JAWRA) 1‐25. DOI: 10.1111/j.1752‐1688.2012.00689.x Abstract: We analyzed an ensemble of watershed models that predict flow, nitrogen, and phosphorus discharges. The models differed in scope and complexity and used different input data, but all had been applied to evaluate human impacts on discharges to the Patuxent River or to the Chesapeake Bay. We compared predictions to observations of average annual, annual time series, and monthly discharge leaving three basins. No model consistently matched observed discharges better than the others, and predictions differed as much as 150% for every basin. Models that agreed best with the observations in one basin often were among the worst models for another material or basin. Combining model predictions into a model average improved overall reliability in matching observations, and the range of predictions helped describe uncertainty. The model average was not the closest to the observed discharge for every material, basin, and time frame, but the model average had the highest Nash–Sutcliffe performance across all combinations. Consistently poor performance in predicting phosphorus loads suggests that none of the models capture major controls. Differences among model predictions came from differences in model structures, input data, and the time period considered, and also to errors in the observed discharge. Ensemble watershed modeling helped identify research needs and quantify the uncertainties that should be considered when using the models in management decisions.     

Watershed conservation and preservation: Environmental engagement as helping behavior

When will a community's residents take action against urban sprawl that threatens the watershed where they live? Drawing on theoretical and empirical studies of helping behavior, we predicted that individuals will be most likely to respond to environmental challenges when they are aware of the environmental threat, believe the danger posed by the threat to be great, and feel responsible for addressing environmental problems. We tested this awareness–appraisal–responsibility (AAR) model by surveying watershed residents’ awareness of watershed features, appraisal of watershed quality, sense of responsibility for protecting the watershed, and behavioral and contribution intentions. Structural equation modeling supported the model by confirming that resident's awareness and appraisal of their local watershed are related to their pro-environment behavioral and contribution intentions, but also their perception of responsibility for protecting the watershed. Mediational analyses confirmed that these relationships are likely sequentially ordered, with awareness leading to appraisal to responsibility and then behavioral and contribution intentions. The discussion considers the implications of these findings for interventions designed to increase environmental engagement.     

Modeling atmospheric nitrogen deposition and transport in the Chesapeake Bay watershed

Atmospheric deposition of nitrate nitrogen and ammonium nitrogen has been identified as a major factor in the decline of water quality in the Chesapeake Bay. Reports have indicated that atmospheric deposition may account for 25 to 80% of the total nitrogen load entering the bay. However, uncertainties exist regarding the accuracy of the atmospheric deposition inputs, nitrogen retention coefficients, and in-stream nutrient uptake rates used in these studies. This project was designed to reassess the potential inputs of atmospheric nitrogen deposition to the bay through the use of a high-resolution wet deposition model, improved wet and dry deposition and nutrient retention estimates, existing soils and land use data, and geographic information systems software. Model results indicate that the methods used in previous studies may overestimate the contribution of atmospheric nitrate and ammonium deposition to the Chesapeake Bay watershed (CBW). Wet and dry atmospheric nitrate and ammonium nitrogen deposition estimates to the CBW ranged from 52.7 to 141.9 and 41.9 to 60.1 million kg/yr, respectively, between 1984 and 1996. Dry and total atmospheric deposition loads to the watershed are substantially less than previous estimates. Estimates of the percent contribution of atmospherically deposited nitrogen to the Chesapeake Bay represent between 20 and 32% of the total nitrate and ammonium nitrogen load to the watershed from all nitrogen sources. While these estimates are lower than many other published estimates, regression analysis of model parameters, nitrogen retention coefficients, output, and measured in-stream nitrogen loads indicate that the calculated nitrogen loads may still be too high.     

Improvement in HSPF's Low‐Flow Predictions by Implementation of a Power Law Groundwater Storage‐Discharge Relationship

We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km² drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.     

面向管理决策的标准化流域水环境模型评估验证技术框架研究

不规范的流域水环境模型应用增加了决策风险。从过程管理的角度来看,我国尚未针对流域水环境模型的评估与验证建立标准化的技术流程,模型标准化应用水平较低。在总结已有研究成果以及先进管理经验的基础上,本文构建了标准化的流域水环境模型评估验证技术框架,提出了对应用于流域水环境管理决策的模型开展评估验证的基本原则、工作流程和技术要求,并通过案例研究验证了技术框架的可行性。技术框架引入了结构合理性评估、参数识别与灵敏度分析、模拟效果评估、不确定性分析等模型评估验证的关键技术,结合不同的模型类型、决策功能等特征给出了原则性的技术要求和应用建议。研究成果充分考虑了我国的环境管理需求,与现阶段环境模拟技术要求、环境监测能力和数据条件相适应,在理论探讨和技术实现层面具备明确的可行性,将促进我国流域水环境模型的规范化、标准化和本地化应用。     

Urbanization and the Loss of Resource Lands in the Chesapeake Bay Watershed

We made use of land cover maps, and land use change associated with urbanization, to provide estimates of the loss of natural resource lands (forest, agriculture, and wetland areas) across the 168,000 km² Chesapeake Bay watershed. We conducted extensive accuracy assessments of the satellite-derived maps, most of which were produced by us using widely available multitemporal Landsat imagery. The change in urbanization was derived from impervious surface area maps (the built environment) for 1990 and 2000, from which we estimated the loss of resource lands that occurred during this decade. Within the watershed, we observed a 61% increase in developed land (from 5,177 to 8,363 km²). Most of this new development (64%) occurred on agricultural and grasslands, whereas 33% occurred on forested land. Some smaller municipalities lost as much as 17% of their forest lands and 36% of their agricultural lands to development, although in the outlying counties losses ranged from 0% to 1.4% for forests and 0% to 2.6% for agriculture. Fast-growing urban areas surrounded by forested land experienced the most loss of forest to impervious surfaces. These estimates could be used for the monitoring of the impacts of development across the Chesapeake Bay watershed, and the approach has utility for other regions nationwide. In turn, the results and the approach can help jurisdictions set goals for resource land protection and acquisition that are consistent with regional restoration goals.     

NONPOINT SOURCE DISCHARGES OF NUTRIENTS FROM PIEDMONT WATERSHEDS OF CHESAPEAKE BAY1

ABSTRACT: We measured annual discharges of water, sediments, and nutrients from 10 watersheds with differing proportions of agricultural lands in the Piedmont physiographic province of the Chesapeake Bay drainage. Flow-weighted mean concentrations of total N, nitrate, and dissolved silicate in watershed discharges were correlated with the proportion of cropland in the watershed. In contrast, concentrations of P species did not correlate with cropland. Organic P and C correlated with the concentration of suspended particles, which differed among watersheds. Thus, the ratio of N:P:Si in discharges differed greatly among watersheds, potentially affecting N, P or Si limitation of phytoplankton growth in the receiving waters. Simple regression models of N discharge versus the percentage of cropland suggest that croplands discharge 29–42 kg N ha^-1 yr^-1 and other lands discharge 1.2–5.8 kg N ha^-1 yr^-1. We estimated net anthropogenic input of N to croplands and other lands using county level data on agriculture and N deposition from the atmosphere. For most of the study watersheds, N discharge amounted to less than half of the net anthropogenic N.