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.     

Evaluation of Three Watershed‐Scale Pesticide Environmental Transport and Fate Models1

Abstract: The U.S. Environmental Protection Agency (USEPA) Office of Pesticide Programs (OPP) has completed an evaluation of three watershed‐scale simulation models for potential use in Food Quality Protection Act pesticide drinking water exposure assessments. The evaluation may also guide OPP in identifying computer simulation tools that can be used in performing aquatic ecological exposure assessments. Models selected for evaluation were the Soil Water Assessment Tool (SWAT), the Nonpoint Source Model (NPSM), a modified version of the Hydrologic Simulation Program‐Fortran (HSPF), and the Pesticide Root Zone Model‐Riverine Water Quality (PRZM‐RIVWQ) model. Simulated concentrations of the pesticides atrazine, metolachlor, and trifluralin in surface water were compared with field data monitored in the Sugar Creek watershed of Indiana’s White River basin by the National Water Quality Assessment (NAWQA) program. The evaluation not only provided USEPA with experience in using watershed models for estimating pesticide concentration in flowing water but also led to the development of improved statistical techniques for assessing model accuracy. Further, it demonstrated the difficulty of representing spatially and temporally variable soil, weather, and pesticide applications with relatively infrequent, spatially fixed, point estimates. It also demonstrated the value of using monitoring and modeling as mutually supporting tools and pointed to the need to design monitoring programs that support modeling.     

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 Regional Modeling Framework of Phosphorus Sources and Transport in Streams of the Southeastern United States1

García, Ana María, Anne B. Hoos, and Silvia Terziotti, 2011. A Regional Modeling Framework of Phosphorus Sources and Transport in Streams of the Southeastern United States. Journal of the American Water Resources Association (JAWRA) 47(5):991‐1010. DOI: 10.1111/j.1752‐1688.2010.00517.x Abstract: We applied the SPARROW model to estimate phosphorus transport from catchments to stream reaches and subsequent delivery to major receiving water bodies in the Southeastern United States (U.S.). We show that six source variables and five land‐to‐water transport variables are significant (p < 0.05) in explaining 67% of the variability in long‐term log‐transformed mean annual phosphorus yields. Three land‐to‐water variables are a subset of landscape characteristics that have been used as transport factors in phosphorus indices developed by state agencies and are identified through experimental research as influencing land‐to‐water phosphorus transport at field and plot scales. Two land‐to‐water variables – soil organic matter and soil pH – are associated with phosphorus sorption, a significant finding given that most state‐developed phosphorus indices do not explicitly contain variables for sorption processes. Our findings for Southeastern U.S. streams emphasize the importance of accounting for phosphorus present in the soil profile to predict attainable instream water quality. Regional estimates of phosphorus associated with soil‐parent rock were highly significant in explaining instream phosphorus yield variability. Model predictions associate 31% of phosphorus delivered to receiving water bodies to geology and the highest total phosphorus yields in the Southeast were catchments with already high background levels that have been impacted by human activity.     

Investigating the Fate and Transport of Escherichia coli in the Charles River,Boston, Using High‐Resolution Observation and Modeling1

Abstract: The processes affecting the fate and transport of Escherichia coli in surface waters were investigated using high‐resolution observation and modeling. The concentration patterns in Boston’s Charles River were observed during four sampling events with a total of 757 samples, including two spatial surveys with two along‐river (1,500 m length) and three across‐river (600 m length) transects at approximately 25‐m intervals, and two temporal surveys at a fixed location (Community Boating) over seven days at hourly intervals. The data reveal significant spatial and temporal structure at scales not resolved by typical monitoring programs. A mechanistic, time‐variable, three‐dimensional coupled hydrodynamic and water quality model was developed using the ECOMSED and RCA modeling frameworks. The computational grid consists of 3,066 grid cells with average length dimension of 25 m. Forcing functions include upstream and downstream boundary conditions, Stony Brook, and Muddy River (major tributaries) combined sewer overflow (CSO) and non‐CSO discharge and wind. The model generally reproduces the observed spatial and temporal patterns. This includes the presence and absence of a plume in the study area under similar loading, but different hydrodynamic conditions caused by operation of the New Charles River Dam (downstream) and wind. The model also correctly predicts an episode of high concentrations at the time‐series station following seven days of no rainfall. The model has an overall root mean square error (RMSE) of 250 CFU/100 ml and an error rate (above or below the USEPA‐recommended single sample criteria value of 235 CFU/100 ml) of 9.4%. At the time series station, the model has an RMSE of 370 CFU/100 ml and an error rate of 15%.     

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.     

Dissolved Organic Carbon and Estrogen Transport in Surface Runoff From Agricultural Land Receiving Poultry Litter1

Dutta, Sudarshan, Shreeram Inamdar, Jerry Tso, Diana S. Aga, and J. Tom Sims, 2012. Dissolved Organic Carbon and Estrogen Transport in Surface Runoff from Agricultural Land Receiving Poultry Litter. Journal of the American Water Resources Association (JAWRA) 48(3): 558-569. DOI: 10.1111/j.1752-1688.2011.00634.x Abstract: Dissolved organic carbon (DOC) provides a reactive substrate for the transport of organic contaminants with runoff. Very few studies have investigated the export of DOC from agricultural land, especially those receiving manure applications. We investigated exports of DOC in surface runoff from agricultural fields receiving various treatments of poultry litter (raw vs. pelletized). In addition, we also investigated how estrogens in runoff were associated with DOC. Different forms of estrogens studied were: estrone, 17β-estradiol, estriol, and their conjugates. Experimental agricultural plots were 12 m × 5 m long and had reduced tillage and no-till management practices. The aromatic content of DOC was characterized using specific ultraviolet absorbance (SUVA). Flow-weighted concentrations of DOC and SUVA in surface runoff from plots with poultry litter were significantly (p ≤ 0.10) greater than the control (no litter) plots. Compared to pelletized poultry litter, reduced-tillage plots with raw litter yielded higher DOC concentrations and SUVA values. No significant differences (p ≥ 0.10) in DOC and SUVA were observed between litter treatments for plots with no-till. Total estrogen concentrations (including all forms) were positively and significantly (p ≤ 0.10) correlated with DOC. These results can help select and guide agricultural management practices that can reduce the exports of DOC and associated contaminant from agricultural land receiving manure applications.     

Release,Dispersion, and Resuspension of Escherichia coli From Direct Fecal Deposits Under Controlled Flows1

Abstract: Water‐quality standards have been placed on fecal indicator organisms such as Escherichia coli in an attempt to limit the concentrations in water bodies. Cattle can be a significant source of bacteria to water systems, particularly when they are allowed direct access to streams. A flume study was conducted to quantify the effect and understand the transport of E. coli from directly deposited cattle manure. Five steady‐state flows, ranging from 0.00683 to 0.0176 m³/s, were studied and loads from a single cowpie exceeded the U.S. Environmental Protection Agency’s recommended water‐quality standards (235 CFU/100 ml) at each flow over the hour study period. Average E. coli concentrations ranged from 10² to 10⁵ CFU/100 ml over the hour sampling period for all flows. High spatial variations in E. coli concentrations were often seen at each sampling time, with higher concentrations typically at the bottom of the flume. E. coli resuspension was initially greater at 0.5 min after deposition, for the lowest flow (10⁵ CFU/m²/s); however, resuspension rates became similar over time, on the order of 10³ CFU/m²/s. This study demonstrates that the concentrations of E. coli can vary over the water column, and therefore grab samples may inaccurately measure bacteria concentrations and loads in streams. In addition, resuspension rates were often high, so the incorporation of this process into water‐quality models is important for bacteria prediction.     

Population Dynamics of Escherichia coli in Surface Water1

Bucci, Vanni, Marin Vuli?, Xiaodan Ruan, and Ferdi L. Hellweger, 2011. Population Dynamics of Escherichia coli in Surface Water. Journal of the American Water Resources Association (JAWRA) 47(3):611‐619. DOI: 10.1111/j.1752‐1688.2011.00528.x Abstract: Present models of the fate of enteric bacteria in surface water usually assume simple first‐order decay. However, observed decay kinetics are typically biphasic and there is increasing evidence for growth in the ambient environment. These features are similar to the stationary phase of long‐term batch cultures, where the dynamics are the result of repeated cycles of population takeovers by Growth Advantage in Stationary Phase (GASP) mutants able to scavenge nutrients released from dead cells. We investigated the possible role of this mechanism (i.e., growth of an adapted subpopulation) in surface water using laboratory experiments with various strains of Escherichia coli in different types of water. Addition of ampicillin (which only kills dividing cells) to these cultures caused a decrease in the total cell density, demonstrating the presence of a growing subpopulation. Furthermore, long‐term survivors from surface water cultures outcompete naïve cells and take over the population in surface water. This growth advantage is heritable (i.e., aged strains were isolated from a single cell) and consistent with the GASP mechanism, although the responsible loci were not identified. Our results show that E. coli populations in surface water are dynamic and consist of dying naïve cells and growing surface water adapted cells. These results can explain the apparent biphasic decay pattern observed in laboratory experiments. This mechanism may also be important in the ambient environment.     

Hot Spots and Hot Moments in Riparian Zones: Potential for Improved Water Quality Management1

Vidon, Philippe, Craig Allan, Douglas Burns, Tim P. Duval, Noel Gurwick, Shreeram Inamdar, Richard Lowrance, Judy Okay, Durelle Scott, and Steve Sebestyen, 2010. Hot Spots and Hot Moments in Riparian Zones: Potential for Improved Water Quality Management. Journal of the American Water Resources Association (JAWRA) 46(2):278-298. DOI: 10.1111/j.1752-1688.2010.00420.x Abstract: Biogeochemical and hydrological processes in riparian zones regulate contaminant movement to receiving waters and often mitigate the impact of upland sources of contaminants on water quality. These heterogeneous processes have recently been conceptualized as “hot spots and moments” of retention, degradation, or production. Nevertheless, studies investigating the importance of hot phenomena (spots and moments) in riparian zones have thus far largely focused on nitrogen (N) despite compelling evidence that a variety of elements, chemicals, and particulate contaminant cycles are subject to the influence of both biogeochemical and transport hot spots and moments. In addition to N, this review summarizes current knowledge for phosphorus, organic matter, pesticides, and mercury across riparian zones, identifies variables controlling the occurrence and magnitude of hot phenomena in riparian zones for these contaminants, and discusses the implications for riparian zone management of recognizing the importance of hot phenomena in annual solute budgets at the watershed scale. Examples are presented to show that biogeochemical process-driven hot spots and moments occur along the stream/riparian zone/upland interface for a wide variety of constituents. A basic understanding of the possible co-occurrence of hot spots and moments for a variety of contaminants in riparian systems will increase our understanding of the influence of riparian zones on water quality and guide management strategies to enhance nutrient or pollutant removal at the landscape scale.     

Impact of Dredging on Phosphorus Transport in Agricultural Drainage Ditches of the Atlantic Coastal Plain1

Abstract: Drainage ditches can be a key conduit of phosphorus (P) between agricultural soils of the Atlantic Coastal Plain and local surface waters, including the Chesapeake Bay. This study sought to quantify the effect of a common ditch management practice, sediment dredging, on fate of P in drainage ditches. Sediments from two drainage ditches that had been monitored for seven years and had similar characteristics (flow, P loadings, sediment properties) were sampled (0‐5 cm) after one of the ditches had been dredged, which removed fine textured sediments (clay = 41%) with high organic matter content (85 g/kg) and exposed coarse textured sediments (clay = 15%) with low organic matter content (2.2 g/kg). Sediments were subjected to a three‐phase experiment (equilibrium, uptake, and release) in recirculating 10‐m‐long, 0.2‐m‐wide, and 5‐cm‐deep flumes to evaluate their role as sources and sinks of P. Under conditions of low initial P concentrations in flume water, sediments from the dredged ditch released 13 times less P to the water than did sediments from the ditch that had not been dredged, equivalent to 24 mg dissolved P. However, the sediments from the dredged ditch removed 19% less P (76 mg) from the flume water when it was spiked with dissolved P to approximate long‐term runoff concentrations. Irradiation of sediments to destroy microorganisms revealed that biological processes accounted for up to 30% of P uptake in the coarse textured sediments of the dredged ditch and 18% in the fine textured sediments of the undredged ditch. Results indicate that dredging of coastal plain drainage ditches can potentially impact the P buffering capacity of ditches draining agricultural soils with a high potential for P runoff.     

Observations and Modeling of Stream Plunging in an Urban Lake1

Owens, Emmet M., Steven W. Effler, Anthony R. Prestigiacomo, David A. Matthews, and Susan M. O’Donnell, 2012. Observations and Modeling of Stream Plunging in an Urban Lake. Journal of the American Water Resources Association (JAWRA) 48(4): 707‐721. DOI: 10.1111/j.1752‐1688.2012.00646.x Abstract: The plunging behavior of two tributaries in Onondaga Lake, New York, is quantified based on a program of monitoring, process studies, and modeling. The dynamics of buoyancy of the tributaries are resolved with hourly measurements of temperature (T), specific conductance (SC), and turbidity (Tn) at the mouths, and observations every 6 h in the lake. Negative buoyancy of the tributaries is found to diminish and change rapidly during runoff events compared to dry periods. In‐lake patterns of the transport of plunging inflow are resolved for dry weather conditions using a dye tracer, and following a runoff event through measurements of T, SC, and Tn. The hydrodynamic/transport model ELCOM (Estuary Lake and Coastal Ocean Model) is demonstrated to perform well in simulating these patterns. Analyses conducted with the model establish the importance of diurnal effects and in‐lake mixing mediated by wind, the need for temporally detailed measurements during runoff events, and modifications of the plunging behavior of the urban tributary as it passes through a harbor. The model provides critical information to support rehabilitation programs for the lake by quantifying the transport of the two largest tributaries, particularly the distribution of the loads between the upper waters and stratified layers. The model predicts that 10% of the urban tributary load enters the upper waters of the lake within 24 h for a dry weather period; this portion increases to 30% for a runoff event.     

The Use of Carbon and Nitrogen Isotopes to Study Watershed Erosion Processes1

Abstract:  Tracer studies are needed to better understand watershed soil erosion and calibrate watershed erosion models. For the first time, stable nitrogen and carbon isotopes (δ¹⁵N and δ¹³C) and the carbon to nitrogen atomic ratio (C/N) natural tracers are used to investigate temporal and spatial variability of erosion processes within a sub‐watershed. Temporal variability was assessed by comparing δ¹⁵N, δ¹³C, and C/N of eroded‐soils from a non‐equilibrium erosion event immediately following freezing and thawing of surface soils with two erosion events characterized by equilibrium conditions with erosion downcutting. Spatial variability was assessed for the equilibrium events by using the δ¹⁵N and δ¹³C signatures of eroded‐soils to measure the fraction of eroded‐soil derived from rill/interrill erosion on upland hillslopes as compared to headcut erosion on floodplains. In order to perform this study, a number of tasks were carried out including: (1) sampling source‐soils from upland hillslopes and floodplains, (2) sampling eroded‐soils with an in situ trap in the stream of the sub‐watershed, (3) isotopic and elemental analysis of the samples using isotope ratio mass spectrometry, (4) fractioning eroded‐soil to its upland rill/interrill and floodplain headcut end‐members using an unmixing model within a Bayesian Markov Chain Monte Carlo framework, and (5) evaluating tracer unmixing model results by comparison with process‐based erosion prediction models for rill/interrill and headcut erosion processes. Results showed that finer soil particles eroded during the non‐equilibrium event were enriched in δ¹⁵N and δ¹³C tracers and depleted in C/N tracer relative to coarser soil particles eroded during the equilibrium events. Correlation of tracer signature with soil particle size was explainable based on known biogeochemical processes. δ¹⁵N and δ¹³C were also able to distinguish between upland rill/interrill erosion and floodplain headcut erosion, which was due to different plant cover at the erosion sources. Results from the tracer unmixing model highlighted future needs for coupling rill/interrill and headcut erosion prediction models.     

Modeling of Sediment Yield From Anjeni-Gauged Watershed,Ethiopia Using SWAT Model1

Setegn, Shimelis G., Bijan Dargahi, Ragahavan Srinivasan, and Assefa M. Melesse, 2010. Modeling of Sediment Yield From Anjeni-Gauged Watershed, Ethiopia Using SWAT Model. Journal of the American Water Resources Association (JAWRA) 46(3):514-526. DOI: 10.1111/j.1752-1688.2010.00431.x Abstract: The Soil and Water Assessment Tool (SWAT) was tested for prediction of sediment yield in Anjeni-gauged watershed, Ethiopia. Soil erosion and land degradation is a major problem on the Ethiopian highlands. The objectives of this study were to evaluate the performance and applicability of SWAT model in predicting monthly sediment yield and assess the impacts of subbasin delineation and slope discretization on the prediction of sediment yield. Ten years monthly meteorological, flow and sediment data were used for model calibration and validation. The annual average measured sediment yield was 24.6 tonnes/ha. The annual average simulated sediment yield was 27.8 and 29.5 tones/ha for calibration and validation periods, respectively. The study found that the observed values showed good agreement with the simulated sediment yield with Nash-Sutcliffe efficiency (NSE) = 0.81, percent bias (PBIAS) = 28%, RMSE-observations standard deviation ratio (RSR) = 0.23, and coefficient of determination (R²) = 0.86 for calibration and NSE = 0.79, PBIAS = 30%, RSR = 0.29, and R² = 0.84 for validation periods. The model can be used for further analysis of different management scenarios that could help different stakeholders to plan and implement appropriate soil and water conservation strategies.     

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.     

Modeling of Vertical Solute Dispersion in a Sediment Bed Enhanced by Wave‐Induced Interstitial Flow1

Abstract: Mass (solute) transport in a stream or lake sediment bed has a significant effect on chemical mass balances and microbial activities in the sediment. A “1D vertical dispersion model” is a useful tool to analyze or model solute transfer between river or lake water and a sediment bed. Under a motionless water column, solute transfer into and within the sediment bed is by molecular diffusion. However, surface waves or bed forms create periodic pressure waves along the sediment/water interface, which in turn induce flows in the pores of the sediment bed. The enhancement of solute transport by these interstitial periodic flows in the pores has been incorporated in a 1D depth‐dependent “enhanced dispersion coefficient (D_E).” Typically, D_E diminishes exponentially with depth in the sediment bed. Relationships have been developed to estimate D_E as a function of the characteristics of sediment (particle size, hydraulic conductivity, and porosity) and pressure waves (wave length and height). In this paper, we outline and illustrate the calculation of D_E as well as the penetration depth (d_p) of the flow effect. Sample applications to illustrate the computational procedure are provided for dissolved oxygen transfer into a stream gravel bed and release of phosphorus from a lake bed. The sensitivity of the results to input parameter values is illustrated, and compared with the errors obtained when interstitial flow is ignored. Maximum values of D_E near the sediment surface can be on the order of 1 cm²/s in a stream gravel bed with standing waves, and 0.001 cm²/s in a fine sand lake bed under progressive surface waves, much larger than molecular diffusion coefficients.     

Modeling of Escherichia coli Fluxes on a Catchment and the Impact on Coastal Water and Shellfish Quality1

Bougeard, Morgane, Jean‐Claude Le Saux, Nicolas Pérenne, Claire Baffaut, Marc Robin, and Monique Pommepuy, 2011. Modeling of Escherichia coli Fluxes on a Catchment and the Impact on Coastal Water and Shellfish Quality. Journal of the American Water Resources Association (JAWRA) 1‐17. DOI: 10.1111/j.1752‐1688.2011.00520.x Abstract: The simulation of the impact of Escherichia coli loads from watersheds is of great interest for assessing estuarine water quality, especially in areas with shellfish aquaculture or bathing activities. For this purpose, this study investigates a model association based on the Soil and Water Assessment Tool (SWAT) coupled with a hydrodynamic model (MARS 2D; IFREMER). Application was performed on the catchment and estuary of Daoulas area (France). The daily E. coli fluxes simulated by SWAT are taken as an input in the MARS 2D model to calculate E. coli concentrations in estuarine water and shellfish. Model validation is based on comparison of frequencies: a strong relationship was found between calculated and measured E. coli concentrations for river quality (r² = 0.99) and shellfish quality (r² = 0.89). The important influence of agricultural practices and rainfall events on the rapid and large fluctuations in E. coli fluxes from the watershed (reaching three orders of magnitude in <24 hours) is one main result of the study. Response time in terms of seawater quality degradation ranges from one to two days after any important rainfall event (greater than 10 mm/day) and the time for estuary to recover good water quality also mainly depends on the duration of the rainfall. In the estuary, three effects (rainfall, tidal dilution, and manure spreading) have been identified as important influences.     

The Regionalization of National‐Scale SPARROW Models for Stream Nutrients1

Schwarz, Gregory E., Richard B. Alexander, Richard A. Smith, and Stephen D. Preston, 2011. The Regionalization of National‐Scale SPARROW Models for Stream Nutrients. Journal of the American Water Resources Association (JAWRA) 47(5):1151‐1172. DOI: 10.1111/j.1752‐1688.2011.00581.x Abstract: This analysis modifies the parsimonious specification of recently published total nitrogen (TN) and total phosphorus (TP) national‐scale SPAtially Referenced Regressions On Watershed attributes models to allow each model coefficient to vary geographically among three major river basins of the conterminous United States. Regionalization of the national models reduces the standard errors in the prediction of TN and TP loads, expressed as a percentage of the predicted load, by about 6 and 7%. We develop and apply a method for combining national‐scale and regional‐scale information to estimate a hybrid model that imposes cross‐region constraints that limit regional variation in model coefficients, effectively reducing the number of free model parameters as compared to a collection of independent regional models. The hybrid TN and TP regional models have improved model fit relative to the respective national models, reducing the standard error in the prediction of loads, expressed as a percentage of load, by about 5 and 4%. Only 19% of the TN hybrid model coefficients and just 2% of the TP hybrid model coefficients show evidence of substantial regional specificity (more than ±100% deviation from the national model estimate). The hybrid models have much greater precision in the estimated coefficients than do the unconstrained regional models, demonstrating the efficacy of pooling information across regions to improve regional models.     

Modeling Biochemical Oxygen Demand Through the Middle and Lower Savannah River1

Abstract: In order to improve modeling accuracy and general understanding of lotic biochemical oxygen demand (BOD), this study characterized river metabolism with the current Georgia Environmental Protection Division method for the middle and lower Savannah River basin (MLSRB) and several alternative methods developed with 120‐day, long‐term biochemical oxygen demand (LTBOD) data from the MLSRB. The data were a subset of a larger two‐year LTBOD study to characterize and understand BOD in the MLSRB, located approximately between Augusta, Georgia, and Savannah, Georgia, along the border of Georgia and South Carolina. The LTBOD data included total oxygen loss and nitrogen speciation for separately quantifying nitrification. Results support the following insights and opportunities for modeling methods: (1) it is important to modeling accuracy that residuals be checked for even dispersion to avoid areas of over‐ and underprediction; (2) modeling with bounded, yet unfixed, rates is a sufficiently simple alternative to fixed‐rate modeling that can eliminate the need for manual adjustments and provide additional system understanding to inform regulation; (3) if fixed rates modeling is desired, model quality for this system might be improved through revising the current low rate (along with the associated f‐ratio updates) from 0.02/day rate to 0.006/day and potentially adding a new rate at 1.0/day in some cases; and (4) the current 57/43 ratio of slow/fast BOD is reasonable based on the 52/45/3 slow/fast/faster BOD proportions of this study.     

The Role of Headwater Streams in Downstream Water Quality1

Abstract: Knowledge of headwater influences on the water‐quality and flow conditions of downstream waters is essential to water‐resource management at all governmental levels; this includes recent court decisions on the jurisdiction of the Federal Clean Water Act (CWA) over upland areas that contribute to larger downstream water bodies. We review current watershed research and use a water‐quality model to investigate headwater influences on downstream receiving waters. Our evaluations demonstrate the intrinsic connections of headwaters to landscape processes and downstream waters through their influence on the supply, transport, and fate of water and solutes in watersheds. Hydrological processes in headwater catchments control the recharge of subsurface water stores, flow paths, and residence times of water throughout landscapes. The dynamic coupling of hydrological and biogeochemical processes in upland streams further controls the chemical form, timing, and longitudinal distances of solute transport to downstream waters. We apply the spatially explicit, mass‐balance watershed model SPARROW to consider transport and transformations of water and nutrients throughout stream networks in the northeastern United States. We simulate fluxes of nitrogen, a primary nutrient that is a water‐quality concern for acidification of streams and lakes and eutrophication of coastal waters, and refine the model structure to include literature observations of nitrogen removal in streams and lakes. We quantify nitrogen transport from headwaters to downstream navigable waters, where headwaters are defined within the model as first‐order, perennial streams that include flow and nitrogen contributions from smaller, intermittent and ephemeral streams. We find that first‐order headwaters contribute approximately 70% of the mean‐annual water volume and 65% of the nitrogen flux in second‐order streams. Their contributions to mean water volume and nitrogen flux decline only marginally to about 55% and 40% in fourth‐ and higher‐order rivers that include navigable waters and their tributaries. These results underscore the profound influence that headwater areas have on shaping downstream water quantity and water quality. The results have relevance to water‐resource management and regulatory decisions and potentially broaden understanding of the spatial extent of Federal CWA jurisdiction in U.S. waters.