Development of bacteria and benthic total maximum daily loads: a case study, Linville Creek, Virginia

Two total maximum daily load (TMDL) studies were performed for Linville Creek in Rockingham County, Virginia, to address bacterial and benthic impairments. The TMDL program is an integrated watershed management approach required by the Clean Water Act. This paper describes the procedures used by the Center for TMDL and Watershed Studies at Virginia Tech to develop the Linville Creek TMDLs and discusses the key lessons learned from and the ramifications of the procedures used in these and other similar TMDL studies. The bacterial impairment TMDL was developed using the Hydrological Simulation Program-Fortran (HSPF). Fecal coliform loads were estimated through an intensive source characterization process. The benthic impairment TMDL was developed using the Generalized Watershed Loading Function (GWLF) model and the reference watershed approach. The bacterial TMDL allocation scenario requires a 100% reduction in cattle manure direct-deposits to the stream, a 96% reduction in nonpoint-source loadings to the land surface, and a 95% reduction in wildlife direct-deposits to the stream. Sediment was identified as the primary benthic stressor. The TMDL allocation scenario for the benthic impairment requires an overall reduction of 12.3% of the existing sediment loads. Despite the many drawbacks associated with using watershed-scale models like HSPF and GWLF to develop TMDLs, the detailed watershed and pollutant-source characterization required to use these and similar models creates information that stakeholders need to select appropriate corrective measures to address the cause of the water quality impairment when implementing the TMDL.     

TMDL Balance: A Model for Coastal Water Pollutant Loadings

Bacterial contamination accounts for more than 60% of the impairments included on the 2008 Texas 303(d) List. Many of these bacterial impairments are along the Texas Gulf Coast because coastal waters often are regulated for oyster harvesting, which have strict water quality standards. Under the Clean Water Act, each one of these impaired waterbodies requires a total maximum daily load (TMDL) study to be performed. A recent, statewide study recommended the development and application of simple modeling approaches to address the majority of Texas's bacteria TMDLs, including “… simple load duration curve, GIS [geographic information systems], and/or mass balance models.” We developed the TMDL Balance model in response to this recommendation. TMDL Balance is a steady state, mass balance, GIS‐based model for simulating pollutant loads and concentrations in coastal systems. The model uses plug‐flow reactor and continuously‐stirred tank reactor equations to route spatially distributed point and nonpoint source loads through a watershed via overland flow, non‐tidal flow, and tidal flow, decaying the loads via first‐order kinetics. In this paper, we explain the development of the watershed loading portion of the TMDL Balance model, demonstrating the methodology through a case study: computing bacterial loads in the Copano Bay watershed of southeast Texas. The application highlights an example of distributing bacterial sources spatially based on land use data.     

Linear Superposition in Total Maximum Daily Loads Applications: The Steady‐State Response Matrix Revisited1

Abstract: The steady‐state response matrix has historically proved a valuable tool in computing the water quality response to loadings and in providing insight into the relative impact of individual loading sources. The insight obtained may be is particularly useful in modern applications of increasingly complex water quality models to problems involving multiple point and nonpoint sources, such as in the assessment of total maximum daily loads (TMDLs). Where appropriate and the underlying equations linear, the steady‐state response matrix can be used to synthesize the results of more complicated models and present them in a way easily understood by policy makers. A straightforward method is presented for generating the response matrix using complex models, and example applications discussed. Example applications include a simple demonstration; incorporation of the method into the Mississippi Department of Environmental Quality’s STREAM model used in TMDL development; a TMDL modeling study of the Grand Calumet River and Indiana Harbor Canal, Indiana, using CE‐QUAL‐ICM; and a TMDL modeling study of the Big Sunflower River, Mississippi, using the Water Analysis Simulation Program model.     

Water Quality Model Uncertainty Analysis of a Point‐Point Source Phosphorus Trading Program1

Kardos, Josef S. and Christopher C. Obropta, 2011. Water Quality Model Uncertainty Analysis of a Point‐Point Source Phosphorus Trading Program. Journal of the American Water Resources Association (JAWRA) 47(6):1317–1337. DOI: 10.1111/j.1752‐1688.2011.00591.x Abstract: Water quality modeling is a major source of scientific uncertainty in the Total Maximum Daily Load (TMDL) process. The effects of these uncertainties extend to water quality trading programs designed to implement TMDLs. This study examines the effects of water quality model uncertainty on a nutrient trading program. The study builds on previous work to design a phosphorus trading program for the Nontidal Passaic River Basin in New Jersey that would implement the watershed TMDL for total phosphorus (TP). The study identified how water quality model uncertainty affects outcomes of potential trades of TP between wastewater treatment plants. The uncertainty analysis found no evidence to suggest that the outcome of trades between wastewater treatment plants, as compared with command and control regulation, will significantly increase uncertainty in the attainment of dissolved oxygen surface water quality standards, site‐specific chlorophyll a criteria, and reduction targets for diverted TP load at potential hot spots in the watershed. Each simulated trading scenario demonstrated parity with or improvement from the command and control approach at the TMDL critical locations, and low risk of hot spots elsewhere.     

ENHANCING FECAL COLIFORM TOTAL MAXIMUM DAILY LOAD MODELS THROUGH BACTERIAL SOURCE TRACKING1

ABSTRACT: Surface water impairment by fecal coliform bacteria is a water quality issue of national scope and importance. In Virginia, more than 400 stream and river segments are on the Commonwealth's 2002 303(d) list because of fecal coliform impairment. Total maximum daily loads (TMDLs) will be developed for most of these listed streams and rivers. Information regarding the major fecal coliform sources that impair surface water quality would enhance the development of effective watershed models and improve TMDLs. Bacterial source tracking (BST) is a recently developed technology for identifying the sources of fecal coliform bacteria and it may be helpful in generating improved TMDLs. Bacterial source tracking was performed, watershed models were developed, and TMDLs were prepared for three streams (Accotink Creek, Christians Creek, and Blacks Run) on Virginia's 303(d) list of impaired waters. Quality assurance of the BST work suggests that these data adequately describe the bacteria sources that are impairing these streams. Initial comparison of simulated bacterial sources with the observed BST data indicated that the fecal coliform sources were represented inaccurately in the initial model simulation. Revised model simulations (based on BST data) appeared to provide a better representation of the sources of fecal coliform bacteria in these three streams. The coupled approach of incorporating BST data into the fecal coliform transport model appears to reduce model uncertainty and should result in an improved TMDL.     

Concentration‐Duration‐Frequency Curves for Stream Turbidity: Possibilities for Assessing Biological Impairment1

Abstract: Siltation and subsequent biological impairment is a national problem prompting state regulatory agencies to develop sediment total maximum daily loads (TMDL) for many streams. To support TMDL targets for reduced sediment yield in disturbed watersheds, a critical need exists for stream assessments to identify threshold concentrations of suspended sediment that impact aquatic biota. Because of the episodic nature of stream sediment transport, thresholds should not only be a function of sediment concentration, but also of duration and dose frequency. Water quality sondes can collect voluminous amounts of turbidity data, a surrogate for suspended sediment, at intervals that can be used to characterize concentration, duration, and frequency of elevated turbidity events. To characterize turbidity sonde data in an ecologically relevant manner, a methodology for concentration‐duration‐frequency (CDF) curves was developed using turbidity doses that relate to different levels of biological impairment. To illustrate this methodology, turbidity CDF curves were generated for two sites on Little Pigeon River in the Great Smoky Mountains National Park, Tennessee, using over 30,000 sonde data measurements per site for a one‐year period. Utilizing a Poisson arrival approach, turbidity spikes were analyzed stochastically by observing the frequency and duration of recorded events over a turbidity level that relates to a biological dose response. An exponential equation was used to fit duration and frequency of a specified turbidity level to generate concentric‐shaped CDF curves, where at specific turbidities longer durations occurred less frequently and conversely shorter durations occurred more frequently. The significance of the equation fit to the data was accomplished with a Kolmogorov‐Smirnov goodness‐of‐fit test. Our findings showed that the CDF curves derived by an exponential function performed reasonable well, with most curves significant at a 95% confidence level. These CDF curves were then used to demonstrate how they could be used to assess biological impairment, and identify future research needs for improved development of sediment TMDLs.     

Magnitude,Frequency, and Duration Relations for Suspended Sediment in Stable (“Reference”) Southeastern Streams1

Abstract: Sediment is listed as one of the leading causes of water‐quality impairments in surface waters of the United States (U.S.). A water body becomes listed by a State, Territory or Tribe if its designated use is not being attained (i.e., impaired). In many cases, the prescribed designated use is aquatic health or habitat, indicating that total maximum daily loads (TMDL) targets for sediment should be functionally related to this use. TMDL targets for sediment transport have been developed for many ecoregions over the past several years using suspended‐sediment yield as a metric. Target values were based on data from “reference” streams or reaches, defined as those exhibiting geomorphic characteristics of equilibrium. This approach has proved useful to some states developing TMDLs for suspended sediment, although one cannot conclude that if a stream exceeds the target range, the aquatic ecosystem will be adversely impacted. To address this problem, historical flow‐transport and sediment‐transport data from hundreds of sites in the Southeastern U.S. were re‐examined to develop parameters (metrics) such as frequency and duration of sediment concentrations. Sites determined as geomorphically stable from field evaluations and from analysis of gauging‐station records were sorted by ecoregion. Mean‐daily flow data obtained from the U.S. Geological Survey were applied to sediment‐transport rating relations to determine suspended‐sediment load for each day of record. The frequency and duration that a given concentration was equaled or exceeded were then calculated to produce a frequency distribution for each site. “Reference” distributions were created using the stable sites in each ecoregion by averaging all of the distributions at specified exceedance intervals. As with the “reference” suspended‐sediment yields, there is a broad range of frequency and duration distributions that reflects the hydrologic and sediment‐transport regimes of the ecoregions. Ecoregions such as the Mississippi Valley Loess Plains (#74) maintain high suspended sediment concentrations for extended periods, whereas coastal plain ecoregions (#63 and 75) show much lower concentrations.     

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.     

Sediment Source Fingerprinting: Transforming From a Research Tool to a Management Tool1

Abstract: Information on the nature and relative contribution of different watershed sediment sources is recognized as a key requirement in the design and implementation of targeted management strategies for sediment control. A direct method of assessing sediment sources in a watershed that has attracted attention in recent years is sediment fingerprinting. The aim of this article is to describe the development of sediment fingerprinting as a research tool and to consider how the method might be transformed from a research tool to a management tool within a regulatory framework, with special reference to the United States total maximum daily load (TMDL) program. When compared with the current source assessment tools in developing sediment TMDLs, sediment fingerprinting offers considerable improvement as a tool for quantifying sources of sediment in terms of source type (e.g., channel vs. hillslope) as well as spatial location (subwatershed). While developing a conceptual framework for sediment TMDLs, we recognize sediment fingerprinting along with sediment budgeting and modeling as valuable tools in the TMDL process for developing justifiable sediment TMDLs. The discussions presented in this article may be considered as a first step toward streamlining the sediment fingerprinting approach for its wider application in a regulatory framework.     

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.     

Modeling Effects of Natural Flow Restoration on Metals Fate and Transport in a Mountain Stream Impacted by Mine Waste1

Abstract: The effects of natural flow restoration on metals fate and transport in the Upper Tenmile Creek Watershed, Montana, were modeled using the Water Quality Analysis Simulation Program developed by the U.S. Environmental Protection Agency (USEPA). This 50‐km² watershed has over 150 historic abandoned mines, including mine waste rock and tailings, as well as adits discharging acid mine drainage, and is the primary drinking water supply for the City of Helena. Water supply diversions almost completely dewater some stream reaches during summer low flows, but the city is considering a new drinking water source and restoration of natural flows in Tenmile Creek as part of acid mine drainage remediation and broader aquatic habitat restoration. One dimensional steady‐state simulation of total recoverable cadmium, copper, lead, and zinc in the mainstem was performed, and the model was calibrated to June 2000 base‐flow data. Representative low‐flows in August and high‐flow snowmelt conditions in June were modeled using mean monthly natural flow estimates from the U.S. Geological Survey and representative USEPA metals concentrations data. The modeling showed that total recoverable metals concentrations, and especially loads, can vary significantly among input locations and over time in the watershed. Some data gaps limit evaluation of variability and increase uncertainty in several locations. Model results indicated, however, that natural low‐ and high‐flow restoration by itself can reduce some metals concentrations in the mainstem compared to June 2000 values, which were influenced by significant water diversion. Some values (such as Zn) may still exceed standards during natural August low flow due to the remaining high concentrations and loads in the primary inputs to the mainstem. Others (such as Cu) can increase during high flow due to remaining mine waste sources and loading of particulate Cu associated with erosion and transport of solids. Greater than 50% reduction in concentrations and loads from some of the main tributaries may be necessary to meet all standards, especially for potential particulate loads with higher flows in June.     

Surface‐Water Nutrient Conditions and Sources in the United States Pacific Northwest1

Wise, Daniel R. and Henry M. Johnson, 2011. Surface‐Water Nutrient Conditions and Sources in the United States Pacific Northwest. Journal of the American Water Resources Association (JAWRA) 47(5):1110‐1135. DOI: 10.1111/j.1752‐1688.2011.00580.x Abstract: The SPAtially Referenced Regressions On Watershed attributes (SPARROW) model was used to perform an assessment of surface‐water nutrient conditions and to identify important nutrient sources in watersheds of the Pacific Northwest region of the United States (U.S.) for the year 2002. Our models included variables representing nutrient sources as well as landscape characteristics that affect nutrient delivery to streams. Annual nutrient yields were higher in watersheds on the wetter, west side of the Cascade Range compared to watersheds on the drier, east side. High nutrient enrichment (relative to the U.S. Environmental Protection Agency’s recommended nutrient criteria) was estimated in watersheds throughout the region. Forest land was generally the largest source of total nitrogen stream load and geologic material was generally the largest source of total phosphorus stream load generated within the 12,039 modeled watersheds. These results reflected the prevalence of these two natural sources and the low input from other nutrient sources across the region. However, the combined input from agriculture, point sources, and developed land, rather than natural nutrient sources, was responsible for most of the nutrient load discharged from many of the largest watersheds. Our results provided an understanding of the regional patterns in surface‐water nutrient conditions and should be useful to environmental managers in future water‐quality planning efforts.     

INCORPORATING NATURAL VARIABILITY,UNCERTAINTY, AND RISK INTO WATER QUALITY EVALUATIONS USING DURATION CURVES1

ABSTRACT: Quantifying natural variability, uncertainty, and risk with minimal data is one of the greatest challenges facing those engaged in water quality evaluations, such as development of total maximum daily loads (TMDL), because of regulatory, natural, and analytical constraints. Quantification of uncertainty and variability in natural systems is illustrated using duration curves (DCs), plots that illustrate the percent of time that a particular flow rate (FDC), concentration (CDC), or load rate (LDC; “TMDL”) is exceeded, and are constructed using simple derived distributions. Duration curves require different construction methods and interpretations, depending on whether there is a statistically significant correlation between concentration (C) and flow (Q), and on the sign of the C‐Q regression slope (positive or negative). Flow DCs computed from annual runoff data vary compared with an FDC developed using all data. Percent exceedance for DCs can correspond to risk; however, DCs are not composed of independent quantities. Confidence intervals of data about a regression line can be used to develop confidence limits for the CDC and LDC. An alternate expression to a fixed TMDL is suggested as the risk of a load rate being exceeded and lying between confidence limits. Averages over partial ranges of DCs are also suggested as an alternative expression of TMDLs. DCs can be used to quantify watershed response in terms of changes in exceedances, concentrations, and load rates after implementation of best management practices.     

ESTIMATING TMDL BACKGROUND SUSPENDED SEDIMENT LOADING TO GREAT LAKES TRIBUTARIES FROM EXISTING DATA1

ABSTRACT: The total maximum daily load (TMDL) for suspended sediment is the maximum quantity of suspended sediment that can enter a waterway without affecting the beneficial uses of the waterway. It is calculated as the sum of permissible allotments of point sources of suspended sediment, permissible allotments of nonpoint sources of suspended sediment, background (natural) loading of suspended sediment, and a margin of safety. The goal of this project was to develop methods for estimating background levels of sediment loads in tributaries of the Great Lakes. Such quantification is key to determining permissible TMDL in waters that do not meet water quality standards under the Clean Water Act of 1972. Suspended sediment loading for 46 rivers was estimated from data collected at U.S. Geological Survey (USGS) gages. Land use and physiographic attributes were estimated for these gaged basins with a geographic information system (GIS). Basin attributes and sediment yield data are the basis for examining two approaches to estimating background suspended sediment loads. One method, based upon envelope curves of sediment yield and drainage area, will be shown to have considerable merit. A second method, based upon correlation of sediment yield to various basin attributes such as drainage area and land use, will be shown to be fraught with difficulties.     

Dry Weather Water Quality Loadings in Arid,Urban Watersheds of the Los Angeles Basin,California, USA1

Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.     

Source and Delivery of Nutrients to Receiving Waters in the Northeastern and Mid‐Atlantic Regions of the United States1

Moore, Richard B., Craig M. Johnston, Richard A. Smith, and Bryan Milstead, 2011. Source and Delivery of Nutrients to Receiving Waters in the Northeastern and Mid‐Atlantic Regions of the United States. Journal of the American Water Resources Association (JAWRA) 47(5):965‐990. DOI: 10.1111/j.1752‐1688.2011.00582.x Abstract: This study investigates nutrient sources and transport to receiving waters, in order to provide spatially detailed information to aid water‐resources managers concerned with eutrophication and nutrient management strategies. SPAtially Referenced Regressions On Watershed attributes (SPARROW) nutrient models were developed for the Northeastern and Mid‐Atlantic (NE US) regions of the United States to represent source conditions for the year 2002. The model developed to examine the source and delivery of nitrogen to the estuaries of nine large rivers along the NE US Seaboard indicated that agricultural sources contribute the largest percentage (37%) of the total nitrogen load delivered to the estuaries. Point sources account for 28% while atmospheric deposition accounts for 20%. A second SPARROW model was used to examine the sources and delivery of phosphorus to lakes and reservoirs throughout the NE US. The greatest attenuation of phosphorus occurred in lakes that were large relative to the size of their watershed. Model results show that, within the NE US, aquatic decay of nutrients is quite limited on an annual basis and that we especially cannot rely on natural attenuation to remove nutrients within the larger rivers nor within lakes with large watersheds relative to the size of the lake.     

Improving the TMDL Process Using Watershed Risk Assessment Principles

Ecological risk assessment (ERA) evaluates potential causal relationships between multiple sources and stressors and impacts on valued ecosystem components. ERAs applied at the watershed scale have many similarities to the place-based analyses that are undertaken to develop Total Maximum Daily Loads (TMDLs), in which linkages are established between stressors, sources, and water quality standards, including support of designated uses. TMDLs focus on achieving water quality standards associated with attainment of designated uses. In attempting to attain the water quality standard, many TMDLs focus on the stressor of concern rather than the ecological endpoint or indicators of the designated use that the standard is meant to protect. A watershed ecological risk assessment (WERA), at least in theory, examines effects of most likely stressors, as well as their probable sources in the watershed, to prioritize management options that will most likely result in meeting environmental goals or uses. Useful WERA principles that can be applied to TMDL development include: development and use of comprehensive conceptual models in the Problem Identification step of TMDLs; use of a transparent process for selecting Numeric Targets for TMDLs based on assessment endpoints derived from the management goal or designated use under consideration; analysis of co-occurring stressors likely to cause beneficial use impairment based on the conceptual model; use of explicit uncertainty analyses in the Linkage Analysis step of TMDL development; and frequent stakeholder interactions throughout the process. WERA principles are currently most applicable to those TMDLs in which there is no numeric standard and, therefore, indicators and targets need to be developed, such as many nutrient or sediment TMDLs. WERA methods can also be useful in determining TMDL targets in situations where simply targeting the water quality standard may re-attain the numeric criterion but not the broader designated use. Better incorporation of problem formulation principles from WERA into the TMDL development process would be helpful in improving the scientific rigor of TMDLs.     

GROUPWISE MODELING STUDY OF BACTERIALLY IMPAIRED WATERSHEDS IN TEXAS: CLUSTERING ANALYSIS1

ABSTRACT: Under the Clean Water Act (CWA) program, the Texas Commission on Environmental Quality (TCEQ) listed 110 stream segments in the year 2000 with pathogenic bacteria impairment. A study was conducted to evaluate the probable sources of pollution and characterize the watersheds associated with these impaired water bodies. The primary aim of the study was to group the water bodies into clusters having similar watershed characteristics and to examine the possibility of studying them as a group by choosing models for total maximum daily load (TMDL) development based on their characteristics. This approach will help to identify possible sources and determine appropriate models and hence reduce the number of required TMDL studies. This in turn will help in reducing the effort required to restore the health of the impaired water bodies in Texas. The main characteristics considered for the classification of water bodies were land use distribution within the watershed, density of stream network, average distance of land of a particular use to the closest stream, household population, density of on‐site sewage facilities (OSSFs), bacterial loading from different types of farm animals and wildlife, and average climatic conditions. The climatic data and observed instream fecal coliform bacteria concentrations were analyzed to evaluate seasonal variability of instream water quality. The grouping of water bodies was carried out using the multivariate statistical techniques of factor analysis/principal component analysis, cluster analysis, and discriminant analysis. The multivariate statistical analysis resulted in six clusters of water bodies. The main factors that differentiated the clusters were found to be bacterial contribution from farm animals and wildlife, density of OSSFs, density of households connected to public sewers, and land use distribution.     

Total Maximum Daily Load Criteria Assessment Using Monitoring and Modeling Data

Applications of Total Maximum Daily Load (TMDL) criteria for complex estuarine systems like Chesapeake Bay have been limited by difficulties in estimating precisely how changes in input loads will impact ambient water quality. A method to deal with this limitation combines the strengths of the Chesapeake Bay's Water Quality Sediment Transport Model (WQSTM), which simulates load response, and the Chesapeake Bay Program's robust historical monitoring dataset. The method uses linear regression to apply simulated relative load responses to historical observations of water quality at a given location and time. Steps to optimize the application of regression analysis were to: (1) determine the best temporal and spatial scale for applying the WQSTM scenarios, (2) determine whether the WQSTM method remained valid with significant perturbation from calibration conditions, and (3) evaluate the need for log transformation of both dissolved oxygen (DO) and chlorophyll a (CHL) datasets. The final method used simple linear regression at the single month, single WQSTM grid cell scale to quantify changes in DO and CHL resulting from simulated load reduction scenarios. The resulting linear equations were applied to historical monitoring data to produce a set of “scenario‐modified” DO or CHL concentration estimates. The utility of the regression method was validated by its ability to estimate progressively increasing attainment in support of the 2010 Chesapeake Bay TMDL.     

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.