Development of the Chesapeake Bay Watershed Total Maximum Daily Load Allocation

Nutrient load allocations and subsequent reductions in total nitrogen and phosphorus have been applied in the Chesapeake watershed since 1992 to reduce hypoxia and to restore living resources. In 2010, sediment allocations were established to augment nutrient allocations supporting the submerged aquatic vegetation resource. From the initial introduction of nutrient allocations in 1992 to the present, the allocations have become more completely applied to all areas and loads in the watershed and have also become more rigorously assessed and tracked. The latest 2010 application of nutrient and sediment allocations were made as part of the Chesapeake Bay total maximum daily load and covered all six states of the Chesapeake watershed. A quantitative allocation process was developed that applied principles of equity and efficiency in the watershed, while achieving all tidal water quality standards through an assessment of equitable levels of effort in reducing nutrients and sediments. The level of effort was determined through application of two key watershed scenarios: one where no action was taken in nutrient control and one where maximum nutrient control efforts were applied. Once the level of effort was determined for different jurisdictions, the overall load reduction was set watershed‐wide to achieve dissolved oxygen water quality standards. Further adjustments were made to the allocation to achieve the James River chlorophyll‐a standard.     

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.     

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.     

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.     

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.     

ESTIMATING NUTRIENT AND SEDIMENT THRESHOLD CRITERIA FOR BIOLOGICAL IMPAIRMENT IN PENNSYLVANIA WATERSHEDS1

ABSTRACT: This study employs a simple nonlinear statistical approach to establish nitrogen, phosphorus, and sediment concentration and unit area load thresholds to aid in the evaluation of aquatic biological health of watersheds within the state of Pennsylvania. Flow, nitrogen and phosphorus species, sediment, basin area, land cover, and biological assessment data were assembled for 29 Pennsylvania watersheds. For each watershed, rating curves depicting flow versus load relationships were developed using the U.S. Environmental Protection Agency's (USEPA's) storage and retrieval database (STORET) flow and concentration data, then applied to daily flow data obtained from U.S. Geological Survey (USGS) daily flow gauging stations to estimate daily load between 1989 and 1999. The load estimates and concentration data were then sorted into six sets of data: mean annual unit area nitrogen, phosphorus, and sediment loads; and average nitrogen, phosphorus, and sediment concentrations. Results of Mann‐Whitney tests conducted on each of the six datasets indicate that there is a statistically significant difference between the concentrations and unit area loads of nitrogen, phosphorus, and sediment in impaired and unimpaired watersheds. Concentration thresholds, calculated as the midpoint between the impaired and unimpaired watersheds’ 95 percent confidence interval for the median, were estimated to be 2.01 mg/L, 0.07 mg/L, and 197.27 mg/L for nitrogen, phosphorus, and sediment, respectively. Annual unit area load thresholds were estimated to be equal to 8.64 kg/ha, 0.30 kg/ha, and 785.29 kg/ha, respectively, for nitrogen, phosphorus, and sediment species.     

Model Calculations of Total Maximum Daily Loads of Mercury for Drainage Lakes1

Abstract: The Watershed Analysis Risk Management Framework watershed model was enhanced to simulate the transport and fate of mercury and to calculate the fish mercury concentrations (FMC) attained by fish through the food web. The model was applied to Western Lake Superior Basin of Minnesota, which has many peat lands and lakes. Topographic, land use, and soil data were used to set up the model. Meteorology and precipitation chemistry data from nearby monitoring stations were compiled to drive the model. Simulated flow and mercury concentrations for several stream stations were comparable to available data. The model was used to perform mercury total maximum daily load calculations for two contrasting drainage lakes (Wild Rice Lake and Whiteface Reservoir). The model results for wet deposition, dry deposition, evasion, watershed yield, and soil sequestration of mercury were comparable with available actual data. The model predicted lake ice cover from November to April and weak stratification in summer, typical of shallow lakes in cold regions. The simulated sulfate decrease and methylmercury increase near the lake bottom in late summer are caused by sulfate reduction and mercury methylation that occur in the surficial sediment. Simulated FMC were within the range of observed values and the R² of correlation between the simulated and observed FMC was 0.77. Under the 1989‐2004 base condition, the average simulated FMC of four‐year‐old walleye was 0.31 μg/g for Whiteface Reservoir and 0.15 μg/g for Wild Rice Lake. The FMC criterion in Minnesota is 0.2 μg/g. Wild Rice Lake already meets this criterion without any load reduction. The model showed that a 65% reduction in atmospheric mercury deposition will not, by itself, allow Whiteface Reservoir to meet the criterion in 15 years. Additional best management practices will be needed to reduce 50% of the watershed input.     

Sources of Suspended-Sediment Flux in Streams of the Chesapeake Bay Watershed: A Regional Application of the SPARROW Model1

Brakebill, John W., Scott W. Ator, and Gregory E. Schwarz, 2010. Sources of Suspended-Sediment Flux in Streams of the Chesapeake Bay Watershed: A Regional Application of the SPARROW Model. Journal of the American Water Resources Association (JAWRA) 46(4): 757-776. DOI: 10.1111/j.1752-1688.2010.00450.x Abstract: We describe the sources and transport of fluvial suspended sediment in nontidal streams of the Chesapeake Bay watershed and vicinity. We applied SPAtially Referenced Regressions on Watershed attributes, which spatially correlates estimated mean annual flux of suspended sediment in nontidal streams with sources of suspended sediment and transport factors. According to our model, urban development generates on average the greatest amount of suspended sediment per unit area (3,928 Mg/km²/year), although agriculture is much more widespread and is the greatest overall source of suspended sediment (57 Mg/km²/year). Factors affecting sediment transport from uplands to streams include mean basin slope, reservoirs, physiography, and soil permeability. On average, 59% of upland suspended sediment generated is temporarily stored along large rivers draining the Coastal Plain or in reservoirs throughout the watershed. Applying erosion and sediment controls from agriculture and urban development in areas of the northern Piedmont close to the upper Bay, where the combined effects of watershed characteristics on sediment transport have the greatest influence may be most helpful in mitigating sedimentation in the bay and its tributaries. Stream restoration efforts addressing floodplain and bank stabilization and incision may be more effective in smaller, headwater streams outside of the Coastal Plain.     

Interactions of warming and altered nutrient load timing on the phenology of oxygen dynamics in Chesapeake Bay

The effects of nutrient loading on estuaries are well studied, given the multitude of negative water quality and ecosystem effects that have been attributed to excess nitrogen and phosphorus. A current gap in this knowledge involves the sensitivity of seasonal cycles of estuarine biogeochemical processes to direct (warming) and indirect influences (nutrient load timing) of climate change. We used a coupled hydrologic–biogeochemical model to investigate changes in the phenology of hypoxia and related biogeochemical processes in Chesapeake Bay under three different hydrologic regimes. Shifts to earlier nutrient load timing during idealized simulations reduced the overall annual hypoxic volume, resulting from discernable, but relatively small reductions in phytoplankton biomass and both sediment and water-column respiration. Simulated increases in water temperature caused an increase in spring/early summer hypoxic volume associated with elevated respiration rates, but an associated exhaustion of organic matter in the early summer caused a decrease in late summer/fall hypoxic volume due to lowered respiration. Warming effects on hypoxia were larger than nutrient timing effects in scenarios where warming was restricted to spring and when it was applied to all months of the year. These idealized simulations begin the process of understanding the potential impacts of future climatic changes in the seasonal timing of key biogeochemical processes associated with eutrophication.     

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.     

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.     

Nutrient and Sediment Production,Watershed Characteristics,and Land Use in the Tahoe Basin,California‐Nevada1

Abstract: In efforts to control the degradation of water quality in Lake Tahoe, public agencies have monitored surface water discharge and concentrations of nitrogen, phosphorus, and suspended sediment in two separate sampling programs. The first program focuses on 20 watersheds varying in size from 162 to 14,000 ha, with continuous stream gaging and periodic sampling; the second focuses on small urbanized catchments, with automated sampling during runoff events. Using data from both programs, we addressed the questions (1) what are the fluxes and concentrations of nitrogen and phosphorus entering the lake from surface runoff; (2) how do the fluxes and concentrations vary in space and time; and (3) how are they related to land use and watershed characteristics? To answer these questions, we calculated discharge‐weighted average concentrations and annual fluxes and used multiple regression to relate those variable to a suite of GIS‐derived explanatory variables. The final selected regression models explain 47‐62% of the variance in constituent concentrations in the stormwater monitoring catchments, and 45‐72% of the variance in mean annual yields in the larger watersheds. The results emphasize the importance of impervious surface and residential density as factors in water quality degradation, and well‐developed soil as a factor in water quality maintenance.     

Hydrologic and Phosphorus Export Behavior of Small Streams in Commercial Poultry‐Pasture Watersheds1

Romeis, J. Joshua, C. Rhett Jackson, L. Mark Risse, Andrew N. Sharpley, and David E. Radcliffe, 2011. Hydrologic and Phosphorus Export Behavior of Small Streams in Commercial Poultry‐Pasture Watersheds. Journal of the American Water Resources Association (JAWRA) 1‐19. DOI: 10.1111/j.1752‐1688.2011.00521.x Abstract: Few watershed‐scale studies have evaluated phosphorus export in streamflow from commercial poultry‐pasture operations. Continuous streamflow and mixed‐frequency water quality datasets were collected from nine commercial poultry‐pasture (AG) and three forested (FORS) headwater streams (2.4‐44 ha) in the upper Etowah River basin of Georgia to estimate total P (TP) loads and examine variability of hydrologic response and water quality of storm and nonstorm‐flow regimes. Data collection duration ranged from 18 to 22 months, and approximately 1,600 water quality samples were collected. Significant (p < 0.1) inverse relationships were detected between peak flow response variables and both drainage area and fraction of forest cover. Order‐of‐magnitude differences in TP and dissolved reactive P (DRP) concentration were observed between AG and FORS sites and among AG sites. TP yields of FORS sites ranged from 0.01 to 0.1 kg P/ha. Yields of AG sites ranged from 0.031 to 3.17 kg P/ha (median = 0.354 kg P/ha). With 95% confidence intervals, AG yields ranged from 0.025 to 13.1 kg P/ha. These small‐watershed‐scale yields were similar to field‐scale yields measured in other studies in other regions. TP yields were significantly related to area‐weighted Mehlich‐1 soil test P concentrations (p = 0.0073) and base‐flow water sample P concentrations (p ≤ 0.0005). Water quality sampling during base‐flow conditions may be a useful screening tool for P risk‐based management programs.     

Maximizing Returns from Payments for Water‐Based Ecosystem Services: Incorporating Externality Effects of Land Management

Given the expansion of payments for water‐based ecosystem services (PWES) worldwide, two relevant issues are as follows: (1) determination of efficient allocations of payments among land managers, and (2) how this might change when paying one manager to implement a best management practice (BMP) to enhance an ecosystem service impacts the cost‐effectiveness of BMPs considered by other land managers not currently involved in PWES. Such externalities may be negative if diminishing returns dominate, or positive if mechanisms such as “social diffusion” dominate. We analyze how a planner should optimally allocate payments, depending on whether the expected externalities are negligible, negative, or positive. We employ (1) an optimal control model to gain insights on the problem’s dynamics, and (2) stochastic dynamic programming to determine optimal funding strategies using a specific application. The study contributes to the literature by identifying dynamically optimal PWES payment patterns, and illustrates how they should change when one accounts for externalities induced by the program. Because such impacts have not been addressed previously in a rigorous way, this treatment provides useful value added for PWES design and implementation.