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.     

Nutrient Inputs to the Laurentian Great Lakes by Source and Watershed Estimated Using SPARROW Watershed Models1

Robertson, Dale M. and David A. Saad, 2011. Nutrient Inputs to the Laurentian Great Lakes by Source and Watershed Estimated Using SPARROW Watershed Models. Journal of the American Water Resources Association (JAWRA) 47(5):1011‐1033. DOI: 10.1111/j.1752‐1688.2011.00574.x Abstract: Nutrient input to the Laurentian Great Lakes continues to cause problems with eutrophication. To reduce the extent and severity of these problems, target nutrient loads were established and Total Maximum Daily Loads are being developed for many tributaries. Without detailed loading information it is difficult to determine if the targets are being met and how to prioritize rehabilitation efforts. To help address these issues, SPAtially Referenced Regressions On Watershed attributes (SPARROW) models were developed for estimating loads and sources of phosphorus (P) and nitrogen (N) from the United States (U.S.) portion of the Great Lakes, Upper Mississippi, Ohio, and Red River Basins. Results indicated that recent U.S. loadings to Lakes Michigan and Ontario are similar to those in the 1980s, whereas loadings to Lakes Superior, Huron, and Erie decreased. Highest loads were from tributaries with the largest watersheds, whereas highest yields were from areas with intense agriculture and large point sources of nutrients. Tributaries were ranked based on their relative loads and yields to each lake. Input from agricultural areas was a significant source of nutrients, contributing ～33‐44% of the P and ～33‐58% of the N, except for areas around Superior with little agriculture. Point sources were also significant, contributing ～14‐44% of the P and 13‐34% of the N. Watersheds around Lake Erie contributed nutrients at the highest rate (similar to intensively farmed areas in the Midwest) because they have the largest nutrient inputs and highest delivery ratio.     

FTABLE Generation Method Effects on Instream Fecal Bacteria Concentrations Simulated With HSPF1

Abstract: Computer simulation models are used extensively for the development of total maximum daily loads (TMDLs). Specifically, the Hydrological Simulation Program‐FORTRAN (HSPF) is used in Virginia for the development of TMDLs for bacteria impairments. HSPF estimates discharge from a reach using function tables (FTABLES). The FTABLE relates stream stage, surface area, and volume to discharge from a reach. In this study, five FTABLE estimation methods were assessed by comparing their effect on various simulation outputs. Four “field‐based” methods used detailed cross‐sectional data collected via site surveys. A fifth “digital‐based” method used digital elevation data in combination with the Natural Resources Conservation Service Regional Hydraulic Geometry Curves. Sets of FTABLEs created using each method were used in simulations of instream bacteria concentration for a Virginia watershed. Several statistics relating to instream bacteria including long‐term average concentration, die‐off, and the violation rate of Virginia’s bacteria criterion were compared. The pair‐wise Student’s t‐test was used for the comparison. The HSPF simulations that used FTABLES estimated from digitally based data consistently produced significantly higher long‐term average instream fecal bacteria concentrations, significantly lower instream fecal bacteria die‐off, which is related to differences in residence time in the streams, and significantly higher water quality criterion violation rates.     

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

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

EXAMINING LAND USE INFLUENCES ON STREAM HABITATS AND MACROINVERTEBRATES: A GIS APPROACH1

ABSTRACT: Geographic Information Systems (GIS) were used to assess the relationships between land use patterns and the physical habitat and macroinvertebrate fauna of streams within similar sized watersheds. Eleven second or third order watersheds ranging from highly urbanized to heavily forested were selected along Lake Superior's North Shore. Land use patterns within the watersheds were quantified using readily available digital land use/land cover information, with a minimum mapping resolution of 16 ha. Physical habitat features, describing substrate characteristics and stream morphology, were characterized at sample points within each stream. Principle component and correlation analyses were used to identify relationships between macroinvertebrates and stream physical habitat, and between habitat and land use patterns. Substrate characteristics and presence of coarse woody debris were found to have the strongest correlations with macreinvertebrate assemblage richness and composition. Agricultural and urban land use was correlated with substrate characteristics. Algal abundance, associated with macroinvertebrate compositional differences, was correlated with housing density and non-forest land covers. The use of readily available spatial data, even at this relatively coarse scale, provides a means to detect the primary relationships between land use and stream habitat quality; finer-resolution GIS databases are needed to assess more subtle influences, such as those due to riparian conditions.     

“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.     

Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed1

Cho, Jaepil, Richard R. Lowrance, David D. Bosch, Timothy C. Strickland, Younggu Her, and George Vellidis, 2010. Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed. Journal of the American Water Resources Association (JAWRA) 46(3):586-602. DOI: 10.1111/j.1752-1688.2010.00436.x Abstract: The Soil and Water Assessment Tool (SWAT) does not fully simulate riparian buffers, but has a simple filter function that is responsive to filter strip width (FILTERW). The objectives of this study were to (1) evaluate SWAT hydrology and water quality response to changes in watershed subdivision levels and different FILTERW configurations and (2) provide guidance for selecting appropriate watershed subdivision for model runs that include the riparian buffer feature through the FILTERW parameter. Watershed subdivision level is controlled by the critical source area (CSA) which defines the minimum drainage area required to form the origin of a stream. SWAT was calibrated on a 15.7 km² subdrainage within the Little River Experimental Watershed, Georgia. The calibrated parameter set was applied to 32 watershed configurations consisting of four FILTERW representations for each of eight CSA levels. Streamflow predictions were stable regardless of watershed subdivision and FILTERW configuration. Predicted sediment and nutrient loads from upland areas decreased as CSA increased when spatial variations of riparian buffers are considered. Sediment and nutrient yield at the watershed outlet was responsive to different combinations of CSA and FILTERW depending on selected in-stream processes. CSA ranges which provide stable sediment and nutrient yields at the watershed outlet was suggested for avoiding significant modifications in selected parameter set.