“Nutrient Inputs to the Laurentian Great Lakes by Source and Watershed Estimated Using SPARROW Watershed Models” by Dale M. Robertson and David A. Saad2

Richards, R. Peter, Ibrahim Alameddine, J. David Allan, David B. Baker, Nathan S. Bosch, Remegio Confesor, Joseph V. DePinto, David M. Dolan, Jeffrey M. Reutter, and Donald Scavia, 2012. Discussion –“Nutrient Inputs to the Laurentian Great Lakes by Source and Watershed Estimated Using SPARROW Watershed Models” by Dale M. Robertson and David A. Saad. Journal of the American Water Resources Association (JAWRA) 1‐10. DOI: 10.1111/jawr.12006 Abstract: Results from the Upper Midwest Major River Basin (MRB3) SPARROW model and underlying Fluxmaster load estimates were compared with detailed data available in the Lake Erie and Ohio River watersheds. Fluxmaster and SPARROW estimates of tributary loads tend to be biased low for total phosphorus and high for total nitrogen. These and other limitations of the application led to an overestimation of the relative contribution of point sources vs. nonpoint sources of phosphorus to eutrophication conditions in Lake Erie, when compared with direct estimates for data‐rich Ohio tributaries. These limitations include the use of a decade‐old reference point (2002), lack of modeling of dissolved phosphorus, lack of inclusion of inputs from the Canadian Lake Erie watersheds and from Lake Huron, and the choice to summarize results for the entire United States Lake Erie watershed, as opposed to the key Western and Central Basin watersheds that drive Lake Erie’s eutrophication processes. Although the MRB3 SPARROW model helps to meet a critical need by modeling unmonitored watersheds and ranking rivers by their estimated relative contributions, we recommend caution in use of the MRB3 SPARRROW model for Lake Erie management, and argue that the management of agricultural nonpoint sources should continue to be the primary focus for the Western and Central Basins of Lake Erie.     

Factors Affecting Stream Nutrient Loads: A Synthesis of Regional SPARROW Model Results for the Continental United States1

Preston, Stephen D., Richard B. Alexander, Gregory E. Schwarz, and Charles G. Crawford, 2011. Factors Affecting Stream Nutrient Loads: A Synthesis of Regional SPARROW Model Results for the Continental United States. Journal of the American Water Resources Association (JAWRA) 47(5):891‐915. DOI: 10.1111/j.1752‐1688.2011.00577.x Abstract: We compared the results of 12 recently calibrated regional SPARROW (SPAtially Referenced Regressions On Watershed attributes) models covering most of the continental United States to evaluate the consistency and regional differences in factors affecting stream nutrient loads. The models – 6 for total nitrogen and 6 for total phosphorus – all provide similar levels of prediction accuracy, but those for major river basins in the eastern half of the country were somewhat more accurate. The models simulate long‐term mean annual stream nutrient loads as a function of a wide range of known sources and climatic (precipitation, temperature), landscape (e.g., soils, geology), and aquatic factors affecting nutrient fate and transport. The results confirm the dominant effects of urban and agricultural sources on stream nutrient loads nationally and regionally, but reveal considerable spatial variability in the specific types of sources that control water quality. These include regional differences in the relative importance of different types of urban (municipal and industrial point vs. diffuse urban runoff) and agriculture (crop cultivation vs. animal waste) sources, as well as the effects of atmospheric deposition, mining, and background (e.g., soil phosphorus) sources on stream nutrients. Overall, we found that the SPARROW model results provide a consistent set of information for identifying the major sources and environmental factors affecting nutrient fate and transport in United States watersheds at regional and subregional scales.     

SPARROW Modeling of Nitrogen Sources and Transport in Rivers and Streams of California and Adjacent States,U.S.

The SPARROW (SPAtially Referenced Regressions On Watershed attributes) model was used to evaluate the spatial distribution of total nitrogen (TN) sources, loads, watershed yields, and factors affecting transport and decay in the stream network of California and portions of adjacent states for the year 2002. The two major TN sources to local catchments on a mass basis were fertilizers and manure (51.7%) and wastewater discharge (15.9%). Other sources contributed < 12%. Fertilizer use is widespread in the Central Valley region of California, and also important in several other regions because of the diversity of California agriculture. Precipitation, sand content of surficial soils, wetlands, and tile drains were important for TN movement to stream reaches. Median streamflow in the study area is about 0.04 m³/s. Aquatic losses of nitrogen were found to be most important in intermittent and small to medium sized streams (0.2‐14 m³/s), while larger streams showed less loss, and therefore are important for TN transport. Nitrogen loss in reservoirs was found to be insignificant, possibly because most of the larger ones are located upstream of nitrogen sources. The model was used to show loadings, sources, and tributary inputs to several major rivers. The information provided by the SPARROW model is useful for determining both the major sources contributing nitrogen to streams and the specific tributaries that transport the load.     

A Multi‐Agency Nutrient Dataset Used to Estimate Loads,Improve Monitoring Design,and Calibrate Regional Nutrient SPARROW Models1

Saad, David A., Gregory E. Schwarz, Dale M. Robertson, and Nathaniel L. Booth, 2011. A Multi‐Agency Nutrient Dataset Used to Estimate Loads, Improve Monitoring Design, and Calibrate Regional Nutrient SPARROW Models. Journal of the American Water Resources Association (JAWRA) 47(5):933‐949. DOI: 10.1111/j.1752‐1688. 2011.00575.x Abstract: Stream‐loading information was compiled from federal, state, and local agencies, and selected universities as part of an effort to develop regional SPAtially Referenced Regressions On Watershed attributes (SPARROW) models to help describe the distribution, sources, and transport of nutrients in streams throughout much of the United States. After screening, 2,739 sites, sampled by 73 agencies, were identified as having suitable data for calculating long‐term mean annual nutrient loads required for SPARROW model calibration. These sites had a wide range in nutrient concentrations, loads, and yields, and environmental characteristics in their basins. An analysis of the accuracy in load estimates relative to site attributes indicated that accuracy in loads improve with increases in the number of observations, the proportion of uncensored data, and the variability in flow on observation days, whereas accuracy declines with increases in the root mean square error of the water‐quality model, the flow‐bias ratio, the number of days between samples, the variability in daily streamflow for the prediction period, and if the load estimate has been detrended. Based on compiled data, all areas of the country had recent declines in the number of sites with sufficient water‐quality data to compute accurate annual loads and support regional modeling analyses. These declines were caused by decreases in the number of sites being sampled and data not being entered in readily accessible databases.     

Estimating Potential E. coli Sources in a Watershed Using Spatially Explicit Modeling Techniques1

Riebschleager, K.J., R. Karthikeyan, R. Srinivasan, and K. McKee, 2012. Estimating Potential E. coli Sources in a Watershed Using Spatially Explicit Modeling Techniques. Journal of the American Water Resources Association (JAWRA) 48(4): 745‐761. DOI: 10.1111/j.1752‐1688.2012.00649.x Abstract: The Spatially Explicit Load Enrichment Calculation Tool (SELECT) was automated to characterize waste and the associated pathogens from various sources within a mixed land use watershed. Potential Escherichia coli loads in Lake Granbury watershed were estimated using spatially variable governing factors, such as land use, soil condition, and distance to streams. A new approach for characterizing E. coli loads resulting from malfunctioning on‐site wastewater treatment systems (OWTSs) was incorporated into SELECT along with the Pollutant Connectivity Factor (PCF) module. The PCF component was applied to identify areas contributing E. coli loads during runoff events by incorporating the influence of potential E. coli loading, runoff potential, and travel distance to waterbodies. Simulation results indicated livestock and wildlife are potential E. coli contributing sources in the watershed. The areas in which these sources are potentially contributing are not currently monitored for E. coli. The bacterial water quality violations seen around Lake Granbury are most likely the result of malfunctioning OWTSs and pet wastes. SELECT results demonstrate the need to evaluate each contributing source separately to effectively allocate site specific best management practices (BMPs) utilizing stakeholder inputs. It also serves as a powerful screening tool for determining areas where detailed investigation is merited.     

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

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

Nutrient Sources and Transport in the Missouri River Basin,with Emphasis on the Effects of Irrigation and Reservoirs1

Brown, Juliane B., Lori A. Sprague, and Jean A. Dupree, 2011. Nutrient Sources and Transport in the Missouri River Basin, With Emphasis on the Effects of Irrigation and Reservoirs. Journal of the American Water Resources Association (JAWRA) 47(5):1034‐1060. DOI: 10.1111/j.1752‐1688.2011.00584.x Abstract: SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were used to relate instream nutrient loads to sources and factors influencing the transport of nutrients in the Missouri River Basin. Agricultural inputs from fertilizer and manure were the largest nutrient sources throughout a large part of the basin, although atmospheric and urban inputs were important sources in some areas. Sediment mobilized from stream channels was a source of phosphorus in medium and larger streams. Irrigation on agricultural land was estimated to decrease the nitrogen load reaching the Mississippi River by as much as 17%, likely as a result of increased anoxia and denitrification in the soil zone. Approximately 16% of the nitrogen load and 33% of the phosphorus load that would have otherwise reached the Mississippi River was retained in reservoirs and lakes throughout the basin. Nearly half of the total attenuation occurred in the eight largest water bodies. Unlike the other major tributary basins, nearly the entire instream nutrient load leaving the outlet of the Platte and Kansas River subbasins reached the Mississippi River. Most of the larger reservoirs and lakes in the Platte River subbasin are upstream of the major sources, whereas in the Kansas River subbasin, most of the source inputs are in the southeast part of the subbasin where characteristics of the area and proximity to the Missouri River facilitate delivery of nutrients to the Mississippi River.     

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.     

DETERMINING NUTRIENT EXPORT COEFFICIENTS AND SOURCE LOADING UNCERTAINTY USING IN‐STREAM MONITORING DATA1

ABSTRACT: Over a three‐year period, flow and nutrients were monitored at 13 sites in the upper North Bosque River watershed in Texas. Drainage areas above sampling sites differed in percent of dairy waste application fields, forage fields, wood/range, and urban land area. A multiple regression approach was used to develop total phosphorus (TP) and total nitrogen (TN) export coefficients for the major land uses in these heterogeneous drainage areas. The largest export coefficients were associated with dairy waste application fields followed by urban, forage fields, and wood/range. An empirical model was then established to assess nutrient contribution by major sources using developed export coefficients and point source loadings from municipal wastewater treatment. This model was verified by comparison of estimated loadings to measured in‐stream data. Monte Carlo simulation techniques were applied to provide an uncertainty analysis for nutrient loads by source, based on the variance associated with each export coefficient. The largest sources of nutrients contributing to the upper North Bosque River were associated with dairy waste application fields and forage fields, while the greatest relative uncertainty in source contribution was associated with loadings from urban and wood/range land uses.     

The Great Lakes Hydrography Dataset: Consistent,Binational Watersheds for the Laurentian Great Lakes Basin

Ecosystem‐based management of the Laurentian Great Lakes, which spans both the United States and Canada, is hampered by the lack of consistent binational watersheds for the entire Basin. Using comparable data sources and consistent methods, we developed spatially equivalent watershed boundaries for the binational extent of the Basin to create the Great Lakes Hydrography Dataset (GLHD). The GLHD consists of 5,589 watersheds for the entire Basin, covering a total area of approximately 547,967 km², or about twice the 247,003 km² surface water area of the Great Lakes. The GLHD improves upon existing watershed efforts by delineating watersheds for the entire Basin using consistent methods; enhancing the precision of watershed delineation using recently developed flow direction grids that have been hydrologically enforced and vetted by provincial and federal water resource agencies; and increasing the accuracy of watershed boundaries by enforcing embayments, delineating watersheds on islands, and delineating watersheds for all tributaries draining to connecting channels. In addition, the GLHD is packaged in a publically available geodatabase that includes synthetic stream networks, reach catchments, watershed boundaries, a broad set of attribute data for each tributary, and metadata documenting methodology. The GLHD provides a common set of watersheds and associated hydrography data for the Basin that will enhance binational efforts to protect and restore the Great Lakes.     

The Lake Erie Agricultural Systems for Environmental Quality project: an introduction

In the last part of the twentieth century, recognition became widespread of the important effect of agricultural runoff on the health of aquatic ecosystems in the Lake Erie basin and elsewhere. Because of the efforts to remediate Lake Erie, the "dead lake" among the Laurentian Great Lakes, a number of research and demonstration projects were undertaken in the Lake Erie basin to evaluate and foster adoption of conservation tillage and other farming techniques that would reduce runoff while maintaining productivity. In addition, intensive water quality studies of long duration were begun on major tributaries to Lake Erie during this time. The Lake Erie Agricultural Systems for Environmental Quality (LEASEQ) project examined governmental programs, changes in agriculture, and changes in water and soil quality during the period 1975-1995, and sought to evaluate the linkages among these factors. The study area is characterized by extensive agricultural land use of soils developed from glacial materials deposited on Paleozoic sedimentary bedrock, mostly limestone. Tile drainage is extensive, particularly in slow-draining clay-rich lacustrine soils in the lower reaches of the watersheds. This paper introduces the study area, its geology, geography, soils, and agricultural history. In addition, we provide an overview of the LEASEQ concept and introduce the 11 other papers in this series, which provide a detailed exposition of the results of our studies.     

Sources and Delivery of Nutrients to the Northwestern Gulf of Mexico from Streams in the South‐Central United States1

Rebich, Richard A., Natalie A. Houston, Scott V. Mize, Daniel K. Pearson, Patricia B. Ging, and C. Evan Hornig, 2011. Sources and Delivery of Nutrients to the Northwestern Gulf of Mexico From Streams in the South‐Central United States. Journal of the American Water Resources Association (JAWRA) 47(5):1061‐1086. DOI: 10.1111/j.1752‐1688.2011.00583.x Abstract: SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were developed to estimate nutrient inputs [total nitrogen (TN) and total phosphorus (TP)] to the northwestern part of the Gulf of Mexico from streams in the South‐Central United States (U.S.). This area included drainages of the Lower Mississippi, Arkansas‐White‐Red, and Texas‐Gulf hydrologic regions. The models were standardized to reflect nutrient sources and stream conditions during 2002. Model predictions of nutrient loads (mass per time) and yields (mass per area per time) generally were greatest in streams in the eastern part of the region and along reaches near the Texas and Louisiana shoreline. The Mississippi River and Atchafalaya River watersheds, which drain nearly two‐thirds of the conterminous U.S., delivered the largest nutrient loads to the Gulf of Mexico, as expected. However, the three largest delivered TN yields were from the Trinity River/Galveston Bay, Calcasieu River, and Aransas River watersheds, while the three largest delivered TP yields were from the Calcasieu River, Mermentau River, and Trinity River/Galveston Bay watersheds. Model output indicated that the three largest sources of nitrogen from the region were atmospheric deposition (42%), commercial fertilizer (20%), and livestock manure (unconfined, 17%). The three largest sources of phosphorus were commercial fertilizer (28%), urban runoff (23%), and livestock manure (confined and unconfined, 23%).     

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.     

Long‐Term Trends of Nutrients and Sediment from the Nontidal Chesapeake Watershed: An Assessment of Progress by River and Season

To assess historical loads of nitrogen (N), phosphorus (P), and suspended sediment (SS) from the nontidal Chesapeake Bay watershed (NTCBW), we analyzed decadal seasonal trends of flow‐normalized loads at the fall‐line of nine major rivers that account for >90% of NTCBW flow. Evaluations of loads by season revealed N, P, and SS load magnitudes have been highest in January‐March and lowest in July‐September, but the temporal trends have followed similar decadal‐scale patterns in all seasons, with notable exceptions. Generally, total N (TN) load has dropped since the late 1980s, but particulate nutrients and SS have risen since the mid‐1990s. The majority of these rises were from Susquehanna River and relate to diminished net trapping at the Conowingo Reservoir. Substantial rises in SS were also observed, however, in other rivers. Moreover, the summed rise in particulate P load from other rivers is of similar magnitude as from Susquehanna. Dissolved nutrient loads have dropped in the upland (Piedmont and above) rivers, but risen in two small rivers in the Coastal Plain affected by lagged groundwater input. In addition, analysis of fractional contributions revealed consistent N trends across the upland watersheds. Finally, total N:total P ratios have declined in most rivers, suggesting the potential for changes in nutrient limitation. Overall, this integrated study of historical data highlights the value of maintaining long‐term monitoring at multiple watershed locations.     

Sources and Transport of Phosphorus to Rivers in California and Adjacent States,U.S., as Determined by SPARROW Modeling

The SPARROW (SPAtially Referenced Regression on Watershed attributes) model was used to simulate annual phosphorus loads and concentrations in unmonitored stream reaches in California, U.S., and portions of Nevada and Oregon. The model was calibrated using de‐trended streamflow and phosphorus concentration data at 80 locations. The model explained 91% of the variability in loads and 51% of the variability in yields for a base year of 2002. Point sources, geological background, and cultivated land were significant sources. Variables used to explain delivery of phosphorus from land to water were precipitation and soil clay content. Aquatic loss of phosphorus was significant in streams of all sizes, with the greatest decay predicted in small‐ and intermediate‐sized streams. Geological sources, including volcanic rocks and shales, were the principal control on concentrations and loads in many regions. Some localized formations such as the Monterey shale of southern California are important sources of phosphorus and may contribute to elevated stream concentrations. Many of the larger point source facilities were located in downstream areas, near the ocean, and do not affect inland streams except for a few locations. Large areas of cultivated land result in phosphorus load increases, but do not necessarily increase the loads above those of geological background in some cases because of local hydrology, which limits the potential of phosphorus transport from land to streams.     

Assessing Risk in Operational Decisions Using Great Lakes Probabilistic Water Level Forecasts

/ A method adapted from the National Weather Service's Extended Streamflow Prediction technique is applied retrospectively to three Great Lakes case studies to show how risk assessment using probabilistic monthly water level forecasts could have contributed to the decision-mak-ing process. The first case study examines the 1985 International Joint Commission (IJC) decision to store water in Lake Superior to reduce high levels on the downstream lakes. Probabilistic forecasts are generated for Lake Superior and Lakes Michigan-Huron and used with riparian inundation value functions to assess the relative impacts of the IJC's decision on riparian interests for both lakes. The second case study evaluates the risk of flooding at Milwaukee, Wisconsin, and the need to implement flood-control projects if Lake Michigan levels were to continue to rise above the October 1986 record. The third case study quantifies the risks of impaired municipal water works operation during the 1964-1965 period of extreme low water levels on Lakes Huron, St. Clair, Erie, and Ontario. Further refinements and other potential applications of the probabilistic forecast technique are discussed.KEY WORDS: Great Lakes; Water levels; Forecasting; Risk; Decision making     

Water Quality Functions of Riparian Forest Buffers in Chesapeake Bay Watersheds

/ Maryland, Virginia, and Pennsylvania, USA, have agreed to reduce nutrient loadings to Chesapeake Bay by 40% by the year 2000. This requires control of nonpoint sources of nutrients, much of which comes from agriculture. Riparian forest buffer systems (RFBS) provide effective control of nonpoint source (NPS) pollution in some types of agricultural watersheds. Control of NPS pollution is dependent on the type of pollutant and the hydrologic connection between pollution sources, the RFBS, and the stream. Water quality improvements are most likely in areas of where most of the excess precipitation moves across, in, or near the root zone of the RFBS. In areas such as the Inner Coastal Plain and Piedmont watersheds with thin soils, RFBS should retain 50%-90% of the total loading of nitrate in shallow groundwater, sediment in surface runoff, and total N in both surface runoff and groundwater. Retention of phosphorus is generally much less. In regions with deeper soils and/or greater regional groundwater recharge (such as parts of the Piedmont and the Valley and Ridge), RFBS water quality improvements are probably much less. The expected levels of pollutant control by RFBS are identified for each of nine physiographic provinces of the Chesapeake Bay Watershed. Issues related to of establishment, sustainability, and management are also discussed.KEY WORDS: Riparian forest buffers; Chesapeake Bay; Nonpoint source pollution; Nitrogen; Phosphorus; Sediment