STOCHASTIC WATER QUALITY ANALYSIS USING RELIABILITY METHOD1

ABSTRACT: This study developed a QUAL2E‐Reliability Analysis (QUAL2E‐RA) model for the stochastic water quality analysis of the downstream reach of the main Han River in Korea. The proposed model is based on the QUAL2E model and incorporates the Advanced First‐Order Second‐Moment (AFOSM) and Mean‐Value First‐Order Second‐Moment (MFOSM) methods. After the hydraulic characteristics from standard step method are identified, the optimal reaction coefficients are then estimated using the Broyden‐Fletcher‐Goldfarb‐Shanno (BFGS) method. Considering variations in river discharges, pollutant loads from tributaries, and reaction coefficients, the violation probabilities of existing water quality standards at several locations in the river were computed from the AFOSM and MFOSM methods, and the results were compared with those from the Monte Carlo method. The statistics of the three uncertainty analysis methods show that the outputs from the AFOSM and MFOSM methods are similar to those from the Monte Carlo method. From a practical model selection perspective, the MFOSM method is more attractive in terms of its computational simplicity and execution time.     

Development of a water quality loading index based on water quality modeling

Water quality modeling is an ideal tool for simulating physical, chemical, and biological changes in aquatic systems. It has been utilized in a number of GIS-based water quality management and analysis applications. However, there is considerable need for a decision-making process to translate the modeling result into an understandable form and thereby help users to make relevant judgments and decisions. This paper introduces a water quality index termed QUAL2E water quality loading index (QWQLI). This new WQI is based on water quality modeling by QUAL2E, which is a popular steady-state model for the water quality of rivers and streams. An experiment applying the index to the Sapgyo River in Korea was implemented. Unlike other WQIs, the proposed index is specifically used for simulated water quality using QUAL2E to mainly reflect pollutant loading levels. Based on the index, an iterative modeling-judgment process was designed to make decisions to decrease input pollutants from pollutant sources. Furthermore, an indexing and decision analysis can be performed in a GIS framework, which can provide various spatial analyses. This can facilitate the decision-making process under various scenarios considering spatial variability. The result shows that the index can evaluate and classify the simulation results using QUAL2E and that it can effectively identify the elements that should be improved in the decision-making process. In addition, the results imply that further study should be carried out to automate algorithms and subsidiary programs supporting the decision-making process.     

COMPUTING THE RISKS ASSOCIATED WITH WASTELOAD) ALLOCATION MODELING1

ABSTRACT: The risks associated with a traditional wasteload allocation (WLA) analysis were quantified with data from a recent study of the Upper Trinity River (Texas). Risk is define here as the probability of failing to meet an established in-stream water quality standard. The QUAL-TX dissolved oxygen (DO) water quality model was modified to a Monte Carlo framework. Flow augmentation coding was also modified to allow an exact match to be computed between the predicted and an established DO concentration standard, thereby providing an avenue for linking input parameter uncertainty to the assignment of a wasteload permit (allowable mass loading rate). Monte Carlo simulation techniques were employed to propagate input parameter uncertainties, typically encountered during WLA analysis, to the computed effluent five-day carbonaceous biochemical oxygen demand requirements for a single major wastewater treatment plant (WWTP). The risk of failing to meet an established in-stream DO criterion may be as high as 96 percent. The uncertainty associated with estimation of the future total Kjeldahl nitrogen concentration for a single tributary was found to have the greatest impact on the determination of allowable WWTP loadings.     

Model‐Based Nitrogen and Phosphorus (Nutrient) Criteria for Large Temperate Rivers: 1. Model Development and Application

An initial inquiry into model‐based numeric nitrogen and phosphorus (nutrient) criteria for large rivers is presented. Field data collection and associated modeling were conducted on a segment of the lower Yellowstone River in the northwestern United States to assess the feasibility of deriving numeric nutrient criteria using mechanistic water‐quality models. The steady‐state one‐dimensional model QUAL2K and a transect‐based companion model AT2K were calibrated and confirmed against low‐flow conditions at a time when river loadings, water column chemistry, and diurnal indicators were approximately steady state. Predictive simulation was then implemented via nutrient perturbation to evaluate the steady‐state and diurnal response of the river to incremental nutrient additions. In this first part of a two‐part series, we detail our modeling approach, model selection, calibration and confirmation, sensitivity analysis, model outcomes, and associated uncertainty. In the second part (Suplee et al., 2015) we describe the criteria development process using the tools described herein. Both articles provide a fundamental understanding of the process required to develop site‐specific numeric nutrient criteria using models in applied regulatory settings.     

Dissolved phosphorus retention and release from a coastal plain in-stream wetland

Dissolved phosphorus (DP) can be released from wetlands as a result of flooding or shifts in water column concentrations. Our objectives were to determine the long-term (1460 d) DP retention and release characteristics of an in-stream wetland, and to evaluate how these characteristics respond to flooding, draining, and changes in DP concentrations. The studied in-stream wetland drains an agriculturally intensive subwatershed in the North Carolina Coastal Plain region. The wetland's DP retention and release characteristics were evaluated by measuring inflow and outflow DP concentrations, DP mass balance, and DP movement across the sediment-water column interface. Phosphorus sorption isotherms were measured to determine the sediment's equilibria P concentration (EPCo), and passive samplers were used to measure sediment pore water DP concentrations. Initially, the in-stream wetland was undersized (0.31 ha) and released 1.5 kg of DP. Increasing the in-stream wetland area to 0.67 ha by flooding resulted in more DP retention (28 kg) and low outflow DP concentrations. Draining the in-stream wetland from 0.67 to 0.33 ha caused the release of stored DP (12.1 kg). Shifts both in sediment pore water DP concentrations and sediment EPCo values corroborate the release of stored DP. Reflooding the wetland from 0.33 to 0.85 ha caused additional release of stored DP into the outflowing stream (10.9 kg). We conclude that for a time period, this in-stream wetland did provide DP retention. During other time periods, DP was released due to changes in wetland area, rainfall, and DP concentrations.     

MODELING THE BIG BLACK RIVER: A COMPARISON OF WATER QUALITY MODELS1

ABSTRACT: The Mississippi Department of Environmental Quality uses the Steady Riverine Environmental Assessment Model (STREAM) to establish effluent limitations. While the U.S. Environmental Protection Agency has approved of its use, questions arise regarding the model's simplicity. The objective of this research was to compare STREAM with the more commonly utilized Enhanced Stream Water Quality Model (QUAL2E). The comparison involved a statistical evaluation procedure based on sensitivity analyses, input probability distribution functions, and Monte Carlo simulation with site‐specific data from a 46‐mile (74‐km) reach of the Big Black River in central Mississippi. Site specific probability distribution functions were derived from measured rates of reaeration, sediment oxygen demand, photosynthesis, and respiration. Both STREAM and QUAL2E reasonably predicted daily average dissolved oxygen (DO) based on a comparison of output probability distributions with observed DO. Observed DO was consistently within 90 percent confidence intervals of model predictions. The STREAM approach generally overpredicted while QUAL2E generally matched observed DO. Using the more commonly assumed lognormal distribution as opposed to a Weibull distribution for two of the sensitive input parameters resulted in minimal differences in the statistical evaluations. The QUAL2E approach had distinct advantages over STREAM in simulating the growth cycle of algae.     

Modeling Water Quantity and Sulfate Concentrations in the Devils Lake Watershed Using Coupled SWAT and CE‐QUAL‐W2

Devils Lake is an endorheic lake in the Red River of the North basin in northeastern North Dakota. During the last two decades, the lake water level has risen by nearly 10 m, causing floods that have cost more than 1 billion USD in mitigation measures. Another increase of approximately 1.5 m in the lake water level would cause spillage into the Sheyenne River. To alleviate this potentially catastrophic spillage, two artificial outlets were constructed. However, the artificial drainage of water into the Sheyenne River raises water quality concerns because the Devils Lake water contains significantly higher concentrations of dissolved solids, particularly sulfate. In this study, the Soil and Water Assessment Tool (SWAT) was coupled with the CE‐QUAL‐W2 model to simulate both water balance and sulfate concentrations in the lake. The SWAT model performed well in simulating daily flow in tributaries with E_NS > 0.5 and |PBIAS| < 25%, and reproduced the lake water level with a root mean square error of 0.35 m for the study period from 1995 to 2014. The water temperature and sulfate concentrations simulated by CE‐QUAL‐W2 for the lake are in general agreement with the field observations. The model results show that the operation of the two outlets since August 2005 has lowered the lake level by 0.70 m. Furthermore, the models show pumping water from the two outlets raises sulfate concentrations in the Sheyenne River from ~100 to >500 mg/L. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.     

Water quality modelling of the river Yamuna (India) using QUAL2E-UNCAS

This paper describes the utility of QUAL2E as a modelling package in the evaluation of a water quality improvement programme. In this study, QUAL2E was applied to determine the pollution loads in the river Yamuna during its course through the national capital territory of Delhi, India. The study aimed at examining the influence of different scenarios on river water quality. Four different pollution scenarios were analysed besides the 'business as usual' situation. The study revealed that it was necessary to treat the discharge from drains to the river Yamuna and diversion of a substantial load to the Agra canal for further treatment was also essential. It was also established through this study that maintaining a flow rate of more than 10 m(3)/s in the river could also help preserve the river's water quality. To clearly display the model outcomes and demarcate polluted zones in the river stretch, model results were interfaced with a Geographical Information System (GIS) to produce cartographic outputs. In addition, uncertainty analysis in the form of first-order error analysis and Monte Carlo analysis was performed, to realise the effect of each model parameter on DO and BOD predictions. The uncertainty analysis gave satisfactory results on simulated data.     

Benthic Algal (Periphyton) Growth Rates in Response to Nitrogen and Phosphorus: Parameter Estimation for Water Quality Models

Nitrogen (N) and phosphorus (P) are significant pollutants that can stimulate nuisance blooms of algae. Water quality models (e.g., Water Quality Simulation Program, CE‐QUAL‐R1, CE‐QUAL‐ICM, QUAL2k) are valuable and widely used management tools for algal accrual due to excess nutrients in the presence of other limiting factors. These models utilize the Monod and Droop equations to associate algal growth rate with dissolved nutrient concentration and intracellular nutrient content. Having accurate parameter values is essential to model performance; however, published values for model parameterization are limited, particularly for benthic (periphyton) algae. We conducted a 10‐day mesocosm experiment and measured diatom‐dominated periphyton biomass accrual through time as chlorophyll a (chl a) and ash‐free dry mass (AFDM) in response to additions of N (range 5–11,995 µg nitrate as nitrogen [NO₃‐N]/L) and P (range 0.89–59.51 µg soluble reactive phosphorus/L). Resulting half‐saturation coefficients and growth rates are similar to other published values, but minimum nutrient quotas are higher than those previously reported. Saturation concentration for N ranged from 150 to 2,450 µg NO₃‐N/L based on chl a and from 8.5 to 60 µg NO₃‐N/L when based on AFDM. Similarly, the saturation concentration for P ranged from 12 to 29 µg‐P/L based on chl a, and from 2.5 to 6.1 µg‐P/L based on AFDM. These saturation concentrations provide an upper limit for streams where diatom growth can be expected to respond to nutrient levels and a benchmark for reducing nutrient concentrations to a point where benthic algal growth will be limited.     

MODELING THE IMPACT OF SUBSURFACE NUTRIENT FLUX ON WATER QUALITY IN THE LOWER TRUCKEE RWER,NEVADA1

ABSTRACT: The impacts of regional groundwater quality and local agricultural activities on in-stream water quality in the Lower Truckee River, Nevada, were assessed through a detailed program of monitoring and computer simulation. An agricultural diversion and return-flow were monitored in great detail to determine mass loading rates of nutrients from agriculture in the area. Once characterized, the cumulative impacts of agricultural diversions and return-flows were evaluated using the Water Quality Assessment Program (WASP) to model nitrogen, phosphorus, periphyton, and dissolved oxygen. Monitoring showed that a significant proportion of the water diverted for agricultural purposes returned to the river as surface point return-flow (estimated at 13.9 percent $ 0.1 percent), and as groundwater diffuse return-flow (estimated at 27 percent $ 6 percent). Modeling efforts demonstrated the significant effect of assumed regional groundwater quality (nitrate) upon predicted periphyton growth and associated diel fluctuations of dissolved oxygen.     

INTEGRATION OF WATER QUALITY MODELING,REMOTE SENSING,AND GIS1

ABSTRACT: Water quality modeling has been developed for more than three quarters of a century, but is limited to the study of trends instead of making accurate short-term forecasts. A major barrier to water quality forecasting is the lack of an efficient system for water quality monitoring. Traditional water quality sampling is time-consuming, expensive, and can only be taken for small sizes. Remote sensing provides a new technique to monitor water quality repetitively for a large area. The objective of this research is to use remotely sensed data in a water quality model - QUAL2E - in a case study of the Te-Chi Reservoir in Taiwan. The water quality variables developed from the simulations are displayed in map form. The developed forecasting system is designed to predict water quality variables using remote sensing data as an input to initialize and update water quality conditions.     

Controls on Nutrients Across a Prairie Stream Watershed: Land Use and Riparian Cover Effects

Nutrient inputs generally are increased by human-induced land use changes and can lead to eutrophication and impairment of surface waters. Understanding the scale at which land use influences nutrient loading is necessary for the development of management practices and policies that improve water quality. The authors assessed the relationships between land use and stream nutrients in a prairie watershed dominated by intermittent stream flow in the first-order higher elevation reaches. Total nitrogen, nitrate, and phosphorus concentrations were greater in tributaries occupying the lower portions of the watershed, closely mirroring the increased density of row crop agriculture from headwaters to lower-elevation alluvial areas. Land cover classified at three spatial scales in each sub-basin above sampling sites (riparian in the entire catchment, catchment land cover, and riparian across the 2 km upstream) was highly correlated with variation in both total nitrogen (r² = 53%, 52%, and 49%, respectively) and nitrate (r² = 69%, 65%, and 56%, respectively) concentrations among sites. However, phosphorus concentrations were not significantly associated with riparian or catchment land cover classes at any spatial scale. Separating land use from riparian cover in the entire watershed was difficult, but riparian cover was most closely correlated with in-stream nutrient concentrations. By controlling for land cover, a significant correlation of riparian cover for the 2 km above the sampling site with in-stream nutrient concentrations could be established. Surprisingly, land use in the entire watershed, including small intermittent streams, had a large influence on average downstream water quality although the headwater streams were not flowing for a substantial portion of the year. This suggests that nutrient criteria may not be met only by managing permanently flowing streams.     

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.     

RIVER PHOSPHORUS DYNAMICS AND RESERVOIR EUTROPHICATION POTENTIAL1

ABSTRACT: Sediments from the Pompton and Passaic Rivers at Two Bridges were analyzed for potentially available phosphorus fractions and total phosphorus (TP). Water samples from the same sites were analyzed for dissolved phosphorus, TP, suspended solids (SS), and volatile SS. Significant negative correlations between river TP concentrations and flow were observed. However, storm flows resulted in increases in TP and SS concentrations and flux (loadings). Most of the increase in river P loading at high flow was in the dissolved fraction, suggesting that the sediments may be a large source of dissolved P. Concentrations of potentially available P in the sediments ranged from 140 to 1310 times the TP concentration in the overlying water. According to a modified Vollenweider model, current P concentrations in the Pompton and Passaic Rivers will result in excessive P loading in the Wanaque Reservoir if even small volumes of river water are pumped to the reservoir through the recently completed Wanaque South pipeline. Reductions in sewage treatment plant effluent P concentrations alone will not produce sufficient decreases in river phosphorus concentrations to avoid this predicted overloading and eutrophication.     

WATERSHED AND LAKE WATER QUALITY ASSESSMENT: AN INTEGRATED MODELING APPROACH1

Watershed planning groups and action agencies seek to understand how lake water quality responds to changes in watershed management. This study developed and demonstrated the applicability of an integrated modeling approach for providing this information. An integrated model linking watershed conditions to water-quality of the receiving lake incorporated the following components: (1) an event-based AGNPS model to estimate watershed pollutant losses; (2) annualization of AGNPS results to produce annual lake pollutant loadings; (3) a base flow separation package, SAM, to estimate base flow; (4) estimates of nutrients in base flow and point sources; and (5) linkage of watershed loadings directly to EUTROMOD lake water quality algorithms. Results are presented for Melvern Lake, a 28-km² multipurpose reservoir with a 900-km² agricultural watershed in east central Kansas. Reasonable estimates of current lake quality were attained using an average phosphorus availability factor of 31 percent to calibrate model results to measured in-lake phosphorus. Comparison of a range of possible scenarios, including all cropland changed to no-till (best case) and all CRP and good-condition grasslands changed to cropland (worst case), indicated only a (4 percent change for in-lake phosphorus and a (2 percent change for chlorophyll a. These results indicated that this watershed is not sensitive to projected changes in land use and management.     

Assessing pollution prevention program by QUAL2E simulation analysis for the Kao-Ping River Basin, Taiwan

Wise and sustainable uses of water resources are essential for an effective river-basin-management planning. Recent management strategy further addresses the fact that quantity and quality of water are closely interrelated, and both must be considered simultaneously for all water resources and water quality management practices. The aim of this paper is to explore the impacts of water resources redistribution and pollution prevention actions between and within river basins simultaneously in South Taiwan. Much emphasis will be placed on assessing the impacts of water transfer over natural boundary to satisfy the needs of industrial development in the Tseng-Wen River system and its resultant influence on the water quality in the downstream area of the Kao-Ping River system where the pollution prevention program is to be implemented. The Kao-Ping River was further characterized hydraulically and environmentally, based on a full investigation of discharges and withdrawals in the river reaches. QUAL2E was successfully calibrated and validated using data collected between 1998 and 1999, and the model was capable of predicting the concentrations of biochemical oxygen demand, dissolved oxygen, total phosphate-phosphorus, and ammonia-nitrogen (NH3-N) for the entire river system. With the aid of QUAL2E simulation model, it shows eliminating the pig farming activities and constructing the sewer systems in the upstream area of Kao-Ping River cannot guarantee the full compliance with water quality standards in the downstream area and water transfer in the upstream area further increases negative impacts on the water quality in the wet season. The predicted situation of water quality in the dry season may even present worse condition. Additional water pollution control policy, such as the use of economic instruments, for controlling and reducing the waste-load of biochemical oxygen demand and ammonia-nitrogen is needed in the Kao-Ping River system in the long run.     

NUTRIENT LOADING ASSESSMENT IN THE ILLINOIS RIVER USING A SYNTHETIC APPROACH1

ABSTRACT: A synthetic relationship is developed between nutrient concentrations and discharge rates at two river gauging sites in the Illinois River Basin. Analysis is performed on data collected by the U.S. Geological Survey (USGS) on nutrients in 1990 through 1997 and 1999 and on discharge rates in 1988 through 1997 and 1999. The Illinois River Basin is in western Arkansas and northeastern Oklahoma and is designated as an Oklahoma Scenic River. Consistently high nutrient concentrations in the river and receiving water bodies conflict with recreational water use, leading to intense stakeholder debate on how best to manage water quality. Results show that the majority of annual phosphorus (P) loading is transported by direct runoff, with high concentrations transported by high discharge rates and low concentrations by low discharge rates. A synthetic relationship is derived and used to generate daily phosphorus concentrations, laying the foundation for analysis of annual loading and evaluation of alternative management practices. Total nitrogen (N) concentration does not have as clear a relationship with discharge. Using a simple regression relationship, annual P loadings are estimated as having a root mean squared error (RMSE) of 39.8 t/yr and 31.9 t/yr and mean absolute percentage errors of 19 percent and 28 percent at Watts and Tahlequah, respectively. P is the limiting nutrient over the full range of discharges. Given that the majority of P is derived from Arkansas, management practices that control P would have the most benefit if applied on the Arkansas side of the border.     

A PHOSPHORUS MANAGEMENT LP MODEL CASE STUDY1

ABSTRACT A linear programming model to assess the cost-effectiveness of appropriate point and nonpoint phosphorus control measures was constructed for Carnegie Lake, an eutrophic lake located in Mercer County, New Jersey. The resultant model was tested for present and future conditions. Feasible solutions were obtained only after significant relaxation of receiving water quality standards. The high levels of point source treatment required to meet the in-stream water quality standards and the mesotrophic loading conditions suggest that a source control solution may not be feasible.     

WASTE LOADINGS FROM A FRESHWATER ATLANTIC SALMON FARM IN SCOTLAND1

ABSTRACT: Recent studies suggest that waste generation from the freshwater phase of Atlantic salmon (Salmo salar L.) production varies considerably on an annual basis. A fish farm on the West Coast of Scotland was visited regularly during a two-year period to determine inflow and outflow water quality. Waste output budgets of suspended solids (SS), biochemical oxygen demand (BOD), total nitrogen (TN), total phosphorus (TP), total ammonia nitrogen (TAN = NH₃+NH₄⁺), dissolved reactive phosphorus (DRP) and total phosphorus (TP) were produced. The annual waste loadings obtained were 71 kg TN t fish^?1 yr^?1 (one year of data only), 10.9–11.1 kg TP t fish^?1 yr^?1, 1.2–2.1 kg DRP t fish^?1 yr^?1, 422–485 kg BOD₅ t fish^?1 yr^?1, 327–337 kg SS t fish^?1 yr^?1, and 30–35 kg TAN-N t fish^?1 yr^?1. Simple linear regression models relating waste parameter production to water temperature and feeding regime were developed. When compared to existing data for other salmonid production systems, greater ranges of daily waste loadings were observed. Wide variations in concentrations of these parameters during a daily cycle were also observed, suggesting that mass balance estimates of waste production will provide more robust estimates of waste output than frequent monitoring of outflow water quality.     

Dissolved phosphorus export from an animal waste impacted in-stream wetland: response to tropical storm and hurricane disturbance

The ability of wetlands to retain P makes them an important landscape feature that buffers P movement. However, their P retention ability can be compromised through hydrologic disturbances caused by hurricanes and tropical storms (TS). This study had three objectives: (i) to determine the effects of hurricanes and TS on dissolved phosphorus (DP) concentrations and loads discharged from a Coastal Plain in-stream wetland (ISW); (ii) to evaluate shifts in P storage pools that would reflect P accretion/removal patterns; and (iii) to determine if relationships exist between storm characteristics with releases of DP and water volume. From January 1996 to October 1999, the ISW's outflow DP concentrations and flow volumes (Q) were measured and they were used to calculate DP mass export loads. In addition, the sediment total phosphorus (TP) concentrations were measured, and both the water column and sediment pore water DP concentrations were examined using passive samplers. In several instances, TS facilitated greater DP releases than a single hurricane event. The largest release of DP occurred in 1999 after Hurricanes Dennis, Floyd, and Irene. The large differences in DP exports among the storms were explained by Q variations. Storm activity also caused changes in sediment pore water DP and sediment TP concentrations. This study revealed that some TS events caused higher DP releases than a single hurricane; however, multiple hurricanes delivering heavy precipitation totals significantly increased DP export.