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.     

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.     

Using Spatially Detailed Water‐Quality Data and Solute‐Transport Modeling to Support Total Maximum Daily Load Development1

Walton‐Day, Katherine, Robert L. Runkel, and Briant A. Kimball, 2012. Using Spatially Detailed Water‐Quality Data and Solute‐Transport Modeling to Support Total Maximum Daily Load Development. Journal of the American Water Resources Association (JAWRA) 48(5): 949‐969. DOI: 10.1111/j.1752‐1688.2012.00662.x Abstract: Spatially detailed mass‐loading studies and solute‐transport modeling using OTIS (One‐dimensional Transport with Inflow and Storage) demonstrate how natural attenuation and loading from distinct and diffuse sources control stream water quality and affect load reductions predicted in total maximum daily loads (TMDLs). Mass‐loading data collected during low‐flow from Cement Creek (a low‐pH, metal‐rich stream because of natural and mining sources, and subject to TMDL requirements) were used to calibrate OTIS and showed spatially variable effects of natural attenuation (instream reactions) and loading from diffuse (groundwater) and distinct sources. OTIS simulations of the possible effects of TMDL‐recommended remediation of mine sites showed less improvement to dissolved zinc load and concentration (14% decrease) than did the TMDL (53‐63% decrease). The TMDL (1) assumed conservative transport, (2) accounted for loads removed by remediation by subtracting them from total load at the stream mouth, and (3) did not include diffuse‐source loads. In OTIS, loads were reduced near their source; the resulting concentration was decreased by natural attenuation and increased by diffuse‐source loads during downstream transport. Thus, by not including natural attenuation and loading from diffuse sources, the TMDL overestimated remediation effects at low flow. Use of the techniques presented herein could improve TMDLs by incorporating these processes during TMDL development.     

A Tool for Automated Load Duration Curve Creation1

Abstract: A recent study by the Texas Bacteria Total Maximum Daily Load (TMDL) Task Force has recommended the use of load duration curves as a primary tool in calculating bacterial TMDLs. This method is attractive because it effectively integrates flow regimes into TMDL analyses, clearly communicates data through a method that is understandable to the general public, and has been successfully applied in TMDL studies in other states. To ease the creation of load duration curves, an automated load duration curve creation tool called LDCurve has been created within a Microsoft Excel framework. Web services and a webscraper are used to retrieve U.S. Geological Survey streamflow data and Texas Commission on Environmental Quality water quality data. Data are imported to the spreadsheet, combined to create flow and load duration curves, and plotted. Final steps result in a preliminary estimate of the overall load reductions needed to meet water quality standards in the modeled segment. LDCurve is currently only applicable in the state of Texas, but may be updated to model water quality throughout the nation using analogous web services from the EPA STORET database. By using automated data retrievals and computations, the LDCurve tool reduces the amount of time required to create curves and calculate load reductions to a matter of minutes. LDCurve and all supporting materials are available online for free download at: http://tools.crwr.utexas.edu/LDCurve/ .     

Assessment of a turfgrass sod best management practice on water quality in a suburban watershed

The disposal of manure on agricultural land has caused water quality concerns in many rural watersheds, sometimes requiring state environmental agencies to conduct total maximum daily load (TMDL) assessments of stream nutrients, such as nitrogen (N) and phosphorus (P). A best management practice (BMP) has been developed in response to a TMDL that mandates a 50% reduction of annual P load to the North Bosque River (NBR) in central Texas. This BMP exports composted dairy manure P through turfgrass sod from the NBR watershed to urban watersheds. The manure-grown sod releases P slowly and would not require additional P fertilizer for up to 20 years in the receiving watershed. This would eliminate P application to the sod and improve the water quality of urban streams. The soil and water assessment tool (SWAT) was used to model a typical suburban watershed that would receive the sod grown with composted dairy manure to assess water quality changes due to this BMP. The SWAT model was calibrated to simulate historical flow and estimated sediment and nutrient loading to Mary's Creek near Fort Worth, Texas. The total P stream loading to Mary's Creek was lower when manure-grown sod was transplanted instead of sod grown with inorganic fertilizers. Flow, sediment and total N yield were the same for both cases at the watershed outlet. The SWAT simulations indicated that the turfgrass BMP can be used effectively to import manure P into an urban watershed and reduce in-stream P levels when compared to sod grown with inorganic fertilizers.     

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.     

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.     

A Benefit‐Cost Analysis of Total Maximum Daily Load Implementation1

Abstract: Total Maximum Daily Load (TMDL) implementation generates benefits and costs from water quality improvements, which are rarely quantified. This analysis examines a TMDL written to address bacteria and aquatic‐life‐use impairments on Abrams and Opequon Creeks in Virginia. Benefits were estimated using a contingent valuation survey of local residents. Costs were based on the number and type of best management practices (BMPs) necessary to achieve TMDL pollution reduction goals. BMPs were quantified using watershed‐scale water quality simulation models (Generalized Watershed Loading Function and Hydrological Simulation Program‐FORTRAN). Based on our projections, the costs to achieve TMDL induced pollution reduction goals outweigh the estimated benefits. Benefit‐cost ratios ranged between 0.1 and 0.3.     

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.     

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.     

Development of the Spatial Rainfall Generator (SRGEN) for the Agricultural Policy/Environmental Extender Model

Accurate spatial representation of climatic patterns is often a challenge in modeling biophysical processes at the watershed scale, especially where the representation of a spatial gradient in rainfall is not sufficiently captured by the number of weather stations. The spatial rainfall generator (SRGEN) is developed as an extension of the “weather generator” (WXGEN), a component of the Agricultural Policy/Environmental eXtender (APEX) model. SRGEN generates spatially distributed daily rainfall using monthly weather statistics available at multiple locations in a watershed. The spatial rainfall generator as incorporated in APEX is tested on the Cowhouse watershed (1,178 km²) in central Texas. The watershed presented a significant spatial rainfall gradient of 2.9 mm/km in the lateral (north‐south) directions based on four rainfall gages. A comparative analysis between SRGEN and WXGEN indicates that SRGEN performs well (PBIAS = 2.40%). Good results were obtained from APEX for streamflow (NSE = 0.99, PBIAS = 8.34%) and NO₃‐N and soluble P loads (PBIAS ≈ 6.00% for each, respectively). However, APEX underpredicted sediment yield and organic N and P loads (PBIAS: 24.75‐27.90%) with SRGEN, although its uncertainty in output was lower than WXGEN results (PBIAS: ?13.02 to ?46.13%). The overall improvement achieved in rainfall generation by SRGEN is demonstrated to be effective in the improving model performance on flow and water quality output.     

LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART II: MODEL APPLICATION1

ABSTRACT: This paper describes the application of a river basin scale hydrologic model (described in Part I) to Richland and Chambers Creeks watershed (RC watershed) in upper Trinity River basin in Texas. The inputs to the model were accumulated from hydro-graphic and geographic databases and maps using a raster-based GIS. Available weather data from 12 weather stations in and around the watershed and stream flow data from two USGS stream gauge station for the period 1965 to 1984 were used in the flow calibration and validation. Sediment calibration was carried out for the period 1988 through 1994 using the 1994 sediment survey data from the Richland-Chambers lake. Sediment validation was conducted on a subwatershed (Mill Creek watershed) situated on Chambers Creek of the RC watershed. The model was evaluated by well established statistical and visual methods and was found to explain at least 84 percent and 65 percent of the variability in the observed stream flow data for the calibration and validation periods, respectively. In addition, the model predicted the accumulated sediment load within 2 percent and 9 percent from the observed data for the RC watershed and Mill Creek watershed, respectively.     

INTEGRATING DIFFUSEINONPOINT POLLUTION CONTROL AND WATER BODY RESTORATION INTO WATERSHED MANAGEMENT1

ABSTRACT: The objective of water quality/watershed management is attainment of water quality goals specified by the Clean Water Act. The Total Maximal Daily Load (TMDL) planning process is a tool to set up watershed management. However, TMDL methodologies and concepts have several problems, including determination of Loading Capacity for only low flow critical periods that preclude consideration of wet weather sources in water quality management. Research is needed to develop watershed pollutant loading and receiving waters Loading Capacity models that will link wet and dry weather pollution loads to the probability of the exceedence of water quality standards. The long term impact of traditional Best Management Practices as well as ponds and wetlands, must be reassessed to consider long term accumulation of conservative toxic compounds. Socioeconomic research should focus on providing information on economic and social feasibility of implementation of additional controls in water quality limited watersheds.     

MODELING ACID MINE DRAINAGE ON A WATERSHED SCALE FOR TMDL CALCULATIONS1

ABSTRACT: The Cheat River of West Virginia is impaired by acid mine drainage (AMD). Fifty‐five of its river segments were placed on the 303(d) list, which required calculations of total maximum daily load (TMDL) to meet the water quality criteria for pH, total iron, aluminum, manganese, and zinc. An existing watershed model was enhanced to simulate AMD as nonpoint source load. The model divided a watershed into a network of catchments and river segments. Each catchment was divided into soil layers, which could contain pyrite, calcite and other minerals. A kinetic expression was used to simulate pyrite oxidation as a function of oxygen in the soil voids. Oxygen in the soil voids was consumed by pyrite oxidation and replenished by earth breathing. The by‐products of pyrite oxidation were calculated according to its mass action equations. Chemical equilibrium was used to account for the speciation of ferrous and ferric irons and precipitation of metal hydroxides. Simulated hydrology and water quality were compared to available data. The USEPA used the calibrated model to calculate the TMDLs in the Cheat River Watershed.     

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.     

Stakeholder Analysis of a Collaborative Watershed Management Process: A Florida Case Study1

Borisova, Tatiana, Laila Racevskis, and Jennison Kipp, 2012. Stakeholder Analysis of a Collaborative Watershed Management Process: A Florida Case Study. Journal of the American Water Resources Association (JAWRA) 48(2): 277‐296. DOI: 10.1111/j.1752-1688.2011.00615.x Abstract: This study focuses on a Florida watershed where development of a total maximum daily load (TMDL) and its implementation plan resulted in conflicts among stakeholders. The overall goal is to build a better understanding of stakeholder perceptions of water quality problems, water policy processes and decisions, and water management plan development in a region where these issues have become contentious. Findings are based on a stakeholder analysis using qualitative data collected through focus groups with agricultural producers, local governments, and environmental groups, and supplemented with additional qualitative data on the watershed management process. Stakeholder conflicts in this case study are associated with perceived flaws in the structural and procedural characteristics of the stakeholder involvement process: (1) suboptimal watershed stakeholder representation on the TMDL executive committee, (2) an inappropriate voting procedure for making TMDL decisions, (3) limitations in information sharing between regulatory agencies and watershed stakeholders, and (4) stakeholders’ doubts about whether tradeoffs associated with achieving the water quality targets were assessed adequately throughout the TMDL planning and implementation process. This study contributes to the literature on collaborative watershed management by analyzing stakeholder involvement given Florida’s unique institutional settings, where implementation of TMDL pollution abatement is mandatory.     

STORMFLOW SIMULATION USING A GEOGRAPHICAL INFORMATION SYSTEM WITH A DISTRIBUTED APPROACH1

ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.     

Mass Load-Based Pollution Management of the Han River and Its Tributaries, Korea

Spatio-temporal variations of biochemical oxygen demand (BOD) and total coliform (TC) in the Han River, Korea, were investigated in terms of concentration-based and mass loading-based approaches. Considering the river water quality criteria regulated by the Ministry of Environment in Korea, the tributaries linked to the mainstream of the Han River were found to be highly contaminated with respect to both BOD and TC and, in fact, most of the tributaries exceeded the maximum water quality criteria. To evaluate the pollution impact of tributaries on the mainstream, the monthly water quality monitoring data for six years (from 1995 to 2000) were collected from the Han River basin, and statistically analyzed using Pearson’s correlation coefficient. The results revealed that mass loading-based approach was superior to the concentration-based approach for effective Han River watershed management. Overall results supported that the mass loading-based approach associated with total maximum daily loads (TMDL) management would be a useful and suitable protocol in watershed management for improving the water quality of the Han River and protecting public health. Therefore, this study supporting TMDL management can be applicable to a wide array of contaminants and watershed settings in Korea.     

Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho1

Stratton, Benjamin T., Venakataramana Sridhar, Molly M. Gribb, James P. McNamara, and Balaji Narasimhan, 2009. Modeling the Spatially Varying Water Balance Processes in a Semiarid Mountainous Watershed of Idaho. Journal of the American Water Resources Association (JAWRA) 45(6):1390‐1408. Abstract: The distributed Soil Water Assessment Tool (SWAT) hydrologic model was applied to a research watershed, the Dry Creek Experimental Watershed, near Boise Idaho to investigate its water balance components both temporally and spatially. Calibrating and validating SWAT is necessary to enable our understanding of the water balance components in this semiarid watershed. Daily streamflow data from four streamflow gages were used for calibration and validation of the model. Monthly estimates of streamflow during the calibration phase by SWAT produced satisfactory results with a Nash Sutcliffe coefficient of model efficiency 0.79. Since it is a continuous simulation model, as opposed to an event‐based model, it demonstrated the limited ability in capturing both streamflow and soil moisture for selected rain‐on‐snow (ROS) events during the validation period between 2005 and 2007. Especially, soil moisture was generally underestimated compared with observations from two monitoring pits. However, our implementation of SWAT showed that seasonal and annual water balance partitioning of precipitation into evapotranspiration, streamflow, soil moisture, and drainage was not only possible but closely followed the trends of a typical semiarid watershed in the intermountain west. This study highlights the necessity for better techniques to precisely identify and drive the model with commonly observed climatic inversion‐related snowmelt or ROS weather events. Estimation of key parameters pertaining to soil (e.g., available water content and saturated hydraulic conductivity), snow (e.g., lapse rates, melting), and vegetation (e.g., leaf area index and maximum canopy index) using additional field observations in the watershed is critical for better prediction.