Modeling the Effects of Macrophytes on Hydrodynamics

A computer model was created as a scientific and management tool for understanding the effects of macrophytes on hydrodynamics and water quality. A model was required that could simulate macrophytes in a complex water body and could be coupled to a multicompartment water quality model of phytoplankton, dissolved oxygen, nutrients, pH, and organic matter. This would permit the investigation of water resource issues where macrophyte growth, phytoplankton growth, nutrient loadings, and flood control were all contributing factors. The model was added as a compartment to the U.S. Army Corps of Engineers two-dimensional, laterally averaged, dynamic water quality model, CE-QUAL-W2 (Corps of Engineers, water quality, width averaged, two dimensional) and applied to the Columbia Slough, Ore. Features of the macrophyte model include the capability to simulate multiple submerged macrophyte species; transport of nutrient fluxes between plant biomass and the water column and/or sediments; growth limitation due to nutrient, light and temperature; simulation of the spatial distribution of macrophytes vertically and horizontally; the modeling of light attenuation in the water column caused by macrophyte concentration; and the modeling of open channel flow with channel friction due to macrophytes. The macrophyte model was tested through mass balances and sensitivity analyses. The modeling of channel friction was evaluated by comparing predicted water levels with data from tests conducted in a laboratory flume. Use of the model in the Columbia Slough showed reasonable predictive capability regarding estimated biomass and water level dynamics.     

Integrated Hydrodynamic and Water Quality Modeling System to Support Nutrient Total Maximum Daily Load Development for Wissahickon Creek, Pennsylvania

This paper presents a hydrodynamic and water quality modeling system for Wissahickon Creek, Pa. Past data show that high nutrient levels in Wissahickon Creek were linked to large diurnal fluctuations in oxygen concentration, which combining with the deoxygenation effect of carbonaceous biological oxygen demand (CBOD) causes violations of dissolved oxygen (DO) standards. To obtain quantitative knowledge about the cause of the DO impairment, an integrated modeling system was developed based on a linked environmental fluid dynamics code (EFDC) and water quality simulation program for eutrophication (WASP/EUTRO5) modeling framework. The EFDC was used to simulate hydrodynamic and temperature in the stream, and the resulting flow information were incorporated into the WASP/EUTRO5 to simulate the fate and transport of nutrients, CBOD, algae, and DO. The standard WASP/EUTRO5 model was enhanced to include a periphyton dynamics module and a diurnal DO simulation module to better represent the prototype. The integrated modeling framework was applied to simulate the creek for a low flow period when monitoring data are available, and the results indicate that the model is a reasonable numerical representation of the prototype.     

Jackson River Modeling: 50‐Year Perspective

The development of water quality models, and also the nature of water quality impairment, is uniquely presented in the point source dissolved oxygen (DO) modeling completed in the Jackson River (Virginia) over the past 50?years. Various water quality modeling studies have been completed in the Jackson River over the years starting with the earliest of modeling frameworks, the Streeter–Phelps equation (1950s and 1960s); progressing to a biochemical oxygen demand–DO model (1970s and 1990s) including diurnal photosynthetic effects (DIURNAL); a Monte Carlo DO analysis using the DIURNAL model (1990s); to the most recent modeling that is currently developing a periphyton model to assess the impact of nutrient loadings on the periphyton community and ultimately DO levels (2000). These early modeling studies were completed by such modeling forefathers as Clarence J. Velz and Donald J. O'Connor, both completing their work at academic institutions (Manhattan College and the University of Michigan) and private consulting firms (Hydroscience and HydroQual, Inc.). Interesting to note is that Earle B. Phelps taught Clarence J. Velz, Donald J. O’Connor’s eventual professor at Manhattan College. Other work completed on the river by early environmental engineers included reaeration studies by Ernest C. Tsivoglou (1966) and the first activated sludge wastewater treatment design for a pulp and paper mill by Wesley Eckenfelder (1950s). The studies investigated: how to improve existing DO conditions in the river; the effects of color reductions on diurnal DO swings; proposed upstream flow regulation effects on water quality and river temperature; and the impact of instream oxygen addition.     

Influence of Physical Forcing on Bottom-Water Dissolved Oxygen within Caloosahatchee River Estuary, Florida

Environmental Fluid Dynamics Code, a numerical estuarine and coastal ocean circulation hydrodynamic and eutrophication model, was used to simulate the distributions of dissolved oxygen (DO), salinity, water temperature, and nutrients in the Caloosahatchee River Estuary. Modeled DO, salinity, and water temperature were in good agreement with field observational data from the Florida Department of Environmental Protection and South Florida Water Management District. Sensitivity analyses identified the effects of river discharge, atmospheric winds, and tidal forcing on the spatial and temporal distributions of DO. Simulation results indicated that vertical mixing due to wind forcing increased the bottom DO concentration. River discharge enhanced stratification in deep locations but propagated vertical mixing in the shallow upper estuary. Finally, tidal forcing heavily influenced bottom layer DO concentrations throughout the whole river estuary.     

In Situ Characterization of Resuspended-Sediment Oxygen Demand in Bubbly Creek, Chicago, Illinois

Sediment oxygen demand (SOD) can be a significant oxygen sink in various types of water bodies, particularly slow-moving waters with substantial organic sediment accumulation. In most settings in which SOD is a concern, the prevailing hydraulic conditions are such that the impact of sediment resuspension on SOD is not considered. However, in the case of Bubbly Creek in Chicago, the prevailing slack water conditions are interrupted by infrequent intervals of very high flow rates associated with pumped combined sewer overflow (CSO) during intense hydrologic events. These events can cause resuspension of the highly organic, nutrient-rich bottom sediments, resulting in precipitous drawdown of dissolved oxygen (DO) in the water column. To address this issue, a new in situ experimental apparatus designed to achieve high flow velocities was implemented to characterize SOD, both with and without sediment resuspension. In the case of resuspension, the suspended sediment concentration was analyzed as a function of bed shear stress, and a formulation was developed to characterize resuspended-sediment oxygen demand (SODR) as a function of suspended sediment concentration in a form similar to first-order biochemical oxygen demand (BOD) kinetics with the DO term in the form of Monod kinetics. The results obtained can be implemented into a model containing hydrodynamic, sediment transport, and water-quality components to yield oxygen demand varying in both space and time for specific flow events. The results are used to evaluate water quality improvement alternatives that take into account the impact of SOD under various flow conditions.     

Hydraulic Study of the Onset of Hypoxia in the Tone River Estuary

Field measurements were conducted to investigate the onset and growth of hypoxia in the Tone River Estuary. Vertical profiles of dissolved oxygen (DO), salinity, and temperature were measured along the deepest line of the estuary. Rates of oxygen consumption by water and sediment in a salt wedge were obtained using laboratory tests. The measurements showed that hypoxia frequently occurred in the front part of the salt wedge and expanded its area toward the river mouth during the summer when the river flow rate was small. The data also suggested that the onset of hypoxia was delayed by the estuarine circulation which supplied oxygen-rich seawater to the salt wedge. To simulate this phenomenon, a two-dimensional flow model in the vertical-longitudinal plane was constructed by transversely integrating the k–ε model equations. The results of model simulation for three months in the summer of 1997 closely matched the field data. The model simulation proved that DO degradation is highly correlated with the residence time of salt water in the estuary.     

Surface Analysis of Chesapeake Bay Water Quality Response to Different Nutrient and Sediment Loads

Based on a set of Chesapeake Bay Estuarine Model (CBEM) scenarios, a three-dimensional response surface of a water quality index, such as chlorophyll concentration, versus a pair of loading constituents, e.g., nitrogen and phosphorus, is constructed. The responses of water quality, such as dissolved oxygen, chlorophyll, and water clarity, to nitrogen, phosphorus, and sediment loads are analyzed. From the response surface, a water quality response is estimated under loading conditions beyond that of a limited set of scenarios. Response surfaces may be used to determine the possible universe of nutrient and sediment load reductions needed to obtain a particular water quality standard and to examine the tradeoffs among nutrient and sediment load reductions that achieve the same water quality objective.     

Long-Term Dynamic Modeling Approach to Quantifying Attached Algal Growth and Associated Impacts on Dissolved Oxygen in the Lower Truckee River, Nevada

Nutrient loads enter the lower Truckee River of western Nevada, affecting the growth of attached algae (periphyton) which causes depressed nighttime dissolved oxygen (DO) levels. The lower Truckee River is home to the endangered cui-ui and threatened Lahontan cut-throat trout, with DO standards being established to in part protect these species. Hydrodynamics, nutrient concentrations, periphyton biomass, and DO data spanning August 2000–December 2001 were used to calibrate and verify a modified version of the Water Quality Analysis Simulation Program Version 5 (WASP5). Under typical loading conditions the periphyton community is nitrogen limited, however nitrogen loading from an upstream wastewater treatment facility increased greatly during the analysis period due to approved site construction activities (discharge permit excursion) causing the periphyton community to temporarily become phosphorus limited. The developed modeling approach, with limited calibration, was able to accurately track dynamic system responses. Removing the impact of the noted discharge permit excursion resulted in a minimum computed DO value of 4.13?mg/L, occurring at the downstream end of the modeling domain on August 8, 2001. Additionally removing the impact of all nutrient loads from area agriculture resulted in a predicted minimum DO value of 4.54?mg/L, while also shifting its location significantly upstream and its timing to April 26, 2001. Meeting all prescribed DO standards required establishing a minimum in-stream flow value of 1.81?m3/s (64.0?ft3/s) downstream of Derby Dam.     

Evaluation of Nutrient Release from Sediments of Artificial Lake

This study examined nutrient fluxes between sediment and water with a laboratory-scale benthic chamber. This research targeted an artificial lake that had undergone water-quality problems. Two sites at Lake Asan, Site A in the vicinity of the dam and Site B at the inflow of the lake, were selected and characteristics of the sediments and their influence on the water quality of the lake were evaluated. Most of the inorganic phosphorus in the study area was in the form of apatite and nonapatite phosphorus (91.9% of Site A, 83.3% of Site B). Site B, with a fast-stream velocity, had larger particle size, smaller nutrient level, and smaller amounts of inorganic phosphorus than Site A. The benthic chamber experiment showed that overall fluxes of Site A were as follows: ammonia is 0.003??μmol?cm-2?day-1, nitrate is -0.067??μmol?cm-2?day-1, and phosphate is 1.049??nmol?cm-2?day-1. Site B showed an increase in phosphate concentration after a dissolved oxygen (DO) drop (<3??mg/l), which resulted in smaller negative nitrate fluxes (-0.043??μmol?cm-2?day-1), larger positive ammonia (0.137??μmol?cm-2?day-1), and larger phosphate fluxes (2.120??nmol?cm-2?day-1) than at Site A. The movement of nutrients at the sediment-water interface was more sensitive to environmental conditions such as DO than other factors, such as sediment characteristics and chemical forms of nutrients. On the basis of the fluxes obtained from the benthic chamber, positive values are estimated for both phosphorus and nitrogen release rates. This indicates that during the sampling period sediment acted as a source of nitrogen as well as phosphorus to the overlying waterbody.     

Dissolved Oxygen Demand at the Sediment-Water Interface of a Stream: Near-Bed Turbulence and Pore Water Flow Effects

A microbial dissolved oxygen (DO) uptake model was developed for a stream bed, including the effect of turbulence in the flow over the bed and pore water flow in the porous bed. The fine-grained sediment bed has hydraulic conductivities 0.01 ≤ k ≤ 1??cm/s, i.e., sediment particle diameter 0.006 ≤ ds ≤ 0.06??cm. The pore water flow is driven by pressure fluctuations at the sediment-water interface, mostly attributable to near-bed coherent motions in the turbulent boundary layer above the sediment bed. An effective mass transfer coefficient (De) coupled to a pore water flow model was used in the DO transport and DO uptake model. DO flux across the sediment-water interface and into the sediment, i.e., sedimentary oxygen demand (SOD), was related to hydraulic conductivity and microbial oxygen uptake rate in the sediment and shear velocity at the sediment-water interface. Simulated SOD values were validated against experimental data. For hydraulic conductivities of the sediment bed up to k ≈ 0.01??cm/s, the pore water flow effect on SOD was found negligible. Above this threshold, the effective mass (DO) transfer coefficient in the sediment bed (De) becomes larger as the hydraulic conductivity (k) becomes larger as the interstitial flow velocities increase; consequently, DO penetration depth increases with larger hydraulic conductivity of the sediment bed (k), and SOD increases as well. The enhancement of vertical DO transport into the sediment bed is strongest near the sediment-water interface, and rapidly diminishes with depth into the sediment layer. An increase in shear velocity at the sediment-water interface also enhances DO transfer. Shear velocity increases at the sediment-water interface will raise SOD regardless of the maximum oxidation rate if the hydraulic conductivity is above the threshold of k ≈ 1??cm/s. The relationship is nearly linear when U*<0.8??cm/s. At shear velocity U* = 1.6??cm/s, SOD for oxidation rates μ = 1000 and 2000??mg?l-1?d-1 are almost five times larger than those with no pore water flow. When pore water transport of DO is not limiting, SOD is a linear function of oxygen demand rate μ in the sediment when 0 ≤ μ ≤ 200??mg?l-1?d-1.     

Diffusional Mass Transfer at Sediment–Water Interface of Cylindrical Sediment Oxygen Demand Chamber

Diffusional mass transfer of dissolved substances across the sediment–water interface in coastal waters is an important factor for realistic determination of sediment oxygen demand (SOD) and nutrient recycle. The benthic diffusive boundary layer inside a cylindrical chamber commonly deployed for in situ measurements of sediment oxygen demand is studied. In a series of laboratory experiments, the SOD is measured with the chamber operated in both continuous flow and batch modes, and a microelectrode is employed to measure the near bed dissolved oxygen (DO) profile for different chamber flows and sediment types. The dependence of the diffusive boundary layer thickness and the sediment–water mass transfer coefficient on the hydraulic parameters are quantified. Using the derived mass transfer coefficient, it is shown that for a given sediment type, the SOD is a function of the bulk DO concentration and chamber flowrate. The theoretical predictions are validated by both laboratory and field SOD data.     

Development and Application of a Phosphorus Model for a Shallow Oxbow Lake

A three-dimensional numerical model was developed for simulating the phosphorus concentration in shallow lakes. In this model, the computational domain was divided into two parts: the water column and the bed sediment layer. The processes of mineralization, settling, adsorption, desorption, bed release (diffusion), growth, and death of phytoplankton were taken into account, and the concentration of organic phosphorus, phosphate, and related water quality constituents was simulated. The concentrations of adsorbed (particulate) and dissolved phosphate due to adsorption-desorption were calculated using two formulas derived based on the Langmuir equation. The release rate of phosphorus from the bed sediment layer was calculated by considering the effects of the concentration gradient across the water-sediment interface, pH, temperature, dissolved oxygen concentration, and flow conditions. The adsorption and desorption of phosphate from sediment particles, as well as its release from bed sediment, were verified using laboratory experimental data. The model was calibrated and applied to Deep Hollow Lake in the Mississippi alluvial plain. The simulated trends and magnitudes of phosphorus concentration were compared with field observations. The simulation results show that there are strong interactions between sediment-related processes and phosphorus concentration.     

Managing for Water Clarity in Chesapeake Bay

Diminished clarity has been listed as a water quality impairment in Chesapeake Bay. The CE-QUAL-ICM eutrophication model has been revised and recalibrated to provide management guidance in alleviating impaired clarity. The algorithms used to model light attenuation and suspended solids are presented herein. Computed and observed total and volatile solids and light attenuation are examined in several formats. A solids budget constructed for the bay identifies major solids sources as internal production, bank erosion, and watershed loading. Sensitivity to loading sources and a key management scenario are examined. Major but feasible reductions in solids and nutrient loads, coupled with reductions in bank erosion, are calculated to meet clarity goals at the 1-m depth in the main bay and major eastern embayments. Careful examination of model results at small scales is required to verify large-scale findings, however. We recommend major improvements in monitoring, computation of light attenuation, and in sediment transport modeling to improve the state of the art in modeling and management of water clarity.     

Development of Refined BOD and DO Models for Highly Polluted Kali River in India

Most commonly used river water quality models for biochemical oxygen demand (BOD) and dissolved oxygen (DO) simulations are mainly based on advection, decay, settling, and loading functions. Using these concepts, refined river water quality models for BOD and DO simulations are developed in the present work considering a large number of physically based parameters and input variables. The refined models developed can be transformed to some of the commonly used river water quality models, if physically based parameters and input variables are omitted or removed. To test the applicability of the refined models developed and commonly used models, a total of 732 water quality and flow data sets are collected during March 1999–February 2000 from 22 sampling stations of the River Kali in India. River Kali is a highly polluted river in India and receives continuous inflow of untreated point source pollution from municipal and industrial wastes and nonpoint source pollution from agricultural areas. Newton–Raphson technique is used to optimize the model parameters during calibration and the performance of different models are evaluated using error estimation, viz. standard error and mean multiplicative error, and correlation statistics (r2). The results indicate that the BOD–DO models proposed by Camp in 1963 provide better results in comparison to other commonly used models. Moreover, the refined models developed for BOD and DO simulations minimize error estimates and improve correlation between observed and computed BOD and DO values of River Kali.     

Modeling the Erosion of Mixtures of Organic and Granular In-Sewer Sediments

The release of fine-grained organic sediments from sediment deposits can have a detrimental impact on water quality in a number of situations. This paper examines the release of such sediment in the context of the erosion of mixed organic/granular sediment in-sewer deposits. In the European Union, sewer flow quality modeling software uses equations derived from uniform granular sediment studies. Actual sewer sediments are mixtures of organic and granular material and interactions between these fractions may account for the poor performance of current models. Laboratory experiments were carried out using surrogate sewer sediment mixtures. Impaction of the bed surface by saltating granular particles increased the erosion of fine-grained organic sediments. Changes in the composition of the bed surface over the duration of a test resulted in change in the availability of fine-grained sediment. A model that attempted to simulate these mechanisms, using an empirically based correction factor to account for the impaction mechanism and an active bed layer to account for changes in the bed surface composition was developed. The limited success of the simulations indicated that such simple modeling approaches may not be appropriate for organic/granular deposits in which grain sorting occurs.     

Oxygen Demand by a Sediment Bed of Finite Length

A model of sedimentary oxygen demand (SOD) for a sediment bed of finite length is presented. The responses of diffusive oxygen transfer in turbulent flow above the sediment surface and of microbial activity inside the sediment to a developing diffusive boundary layer are modeled numerically. The developing diffusive boundary layer above the sediment/water interface is modeled based on shear velocity and turbulent boundary layer concepts, and dissolved oxygen (DO) uptake inside the sediment is modeled as a function of the microbial growth rate. The model predicts that the diffusive boundary layer above the sediment/water interface thickens in flow direction, and that DO penetration depth into the sediment is practically constant over the length of the sediment bed. The effect of the developing diffusive boundary layer on SOD is minor, except at very low shear/flow velocities (shear velocity U*<0.01?cm/s) and/or high microbial density inside the sediment. The average SOD over the sediment bed therefore varies only slightly with its length. SOD varies somewhat in flow direction, i.e., SOD is largest near the leading edge (x = 0), decreases with distance, and finally, approaches a nearly constant value for fully developed boundary layer. Including microbial activity in the sediment makes the change of SOD in flow direction much smaller than is predicted by a pure vertical diffusive flux model. The diffusive boundary layer is nearly fully developed at a dimensionless distance x+ = 10,000, regardless of microbial activity inside the sediment. Longer sediment beds are required to eliminate the small leading edge effect on any measured average SOD value. SOD depends strongly on the diffusion coefficient of DO inside the sediment bed. This effect becomes more significant as shear/flow velocity is increased. Overall, SOD is found to be controlled principally by shear velocity of the water flowing above the sediment/water interface, microbial activity inside the sediment, and diffusion of DO inside the sediment. The length of the sediment bed is of lesser influence.     

Three-Dimensional Ecological-Eutrophication Model for Singapore

A 3D multilevel eutrophication model has been developed and applied to the tropical coastal waters of Singapore. The model considers the nitrogen, phosphorus, and carbon cycles; phytoplankton and zooplankton dynamics; and dissolved oxygen balance in the water column coupled to benthic processes and is driven by a 3D multilevel hydrodynamic model. The eutrophication model was calibrated with field data obtained from the Southwest Monsoon in August 1998 and subsequently validated using data from the Northeast Monsoon in December 1998. Predictions of nutrient and phytoplankton concentrations from the model simulation agreed well with the observed data. The model results show that the coastal waters of Singapore are relatively well-mixed and differences in concentration with depth for all state variables are generally <20%. In the Singapore Strait, horizontal spatial differences for all the variables are generally <40%. In the Southwest Monsoon, slightly larger values for phytoplankton and dissolved oxygen but smaller nutrient values are found in the eastern region of Singapore, compared to the west. The reverse trend is found in the Northeast Monsoon.     

南京生态科技岛海绵城市河道系统植物技术设施对水体重金属及营养盐净化能力评价

为探究海绵城市植物技术设施对重金属及营养盐净化能力,以江心洲南京生态科技岛河道系统为研究对象,对其上下游水体中重金属及营养盐含量进行分析,通过综合污染指数法对上下游重金属风险进行评估,采用冗余分析及Spearman相关系数探究水体环境因素对重金属含量的影响,利用河道区域四种植物技术设施(河岸缓冲带、植被过滤带、生态浮岛、台地式石笼护岸)探究不同植物组合对水体中污染物的净化能力。结果显示,江心洲河道上游重金属的枯水期、丰水期及平水期WQI值分别为1.85、1.74及2.90,分别对应为重金属轻度污染、轻度污染及中度污染。而河道下游重金属的枯水期、丰水期及平水期的WQI值分别为0.18、0.30及0.52,均未存在污染现象。河道上游水体中pH是影响重金属含量的最重要因素,pH与溶解氧(DO)、总氮(TN)、五日化学需氧量(COD₅)及总磷(TP)呈正相关。河道下游水体的pH也是影响水中重金属含量最重要的环境因素,其与溶解氧(DO)呈显著正相关。水体中营养盐净化能力大小为河岸缓冲带>植被过滤带>生态浮岛>台地式石笼护岸。相比其它植物组合,乔灌木群落栽植对水体中营养盐净化最具有潜力。     

Comparison of the Monod and Droop Methods for Dynamic Water Quality Simulations

The Monod method is widely used to model nutrient limitation and primary productivity in water bodies. It offers a straightforward approach to simulate the main processes governing eutrophication and it allows the proper representation of many aquatic systems. The Monod method is not able to represent the nutrient luxury uptake by algae, which consists of the excess nutrient uptake during times of high nutrient availability in the water column. The Droop method, which is also used to model nutrient limitation and primary productivity, takes into account the luxury uptake of nutrients. Because of the relative complexity of the Droop method, it has not been systematically adopted for the simulation of large stream networks. The Water Quality Analysis Simulation Program (WASP) version 7.1 was updated to include nutrient luxury uptake for periphyton growth. The objective of this paper is to present the new nutrient limitation processes simulated by WASP 7.1 and to compare the performance of the Droop and the Monod methods for a complex stream network where periphyton is the main organism responsible for primary productivity. Two applications of WASP 7.1 with the Droop and Monod methods were developed for the Raritan River Basin in New Jersey. Water quality parameters affecting the transport and fate of nutrients were calibrated based on observed data collected for the Raritan River total maximum daily load. The dissolved oxygen and nutrients simulated with WASP 7.1, obtained with the Droop and Monod methods, were compared at selected monitoring stations under different flows and nutrient availability conditions. The comparison of the WASP 7.1 applications showed the importance of using the Droop method when periphyton was the main organism responsible for primary productivity. The data simulated with the Droop method resulted in good agreement with the observed data for dissolved oxygen, ammonia-nitrogen, nitrate-nitrogen, and dissolved orthophosphate at the selected stations. The Monod method was not able to capture the diel dissolved oxygen variation when nutrients were scarce, and it resulted in unrealistic diel variations of nutrients at times of strong primary productivity at some locations.     

Eutrophication and Pathogen Abatement in the San Juan Bay Estuary

The San Juan Bay Estuary is a tropical lagoon system severely impacted by eutrophication and by elevated pathogen concentrations. This investigation examined alternatives for pollution abatement primarily through physical modifications to the system. The investigation included field surveys, computation of pollutant loads, and application of hydrodynamic and eutrophication models. A eutrophication model developed for temperate estuaries was successfully transferred to San Juan Bay Estuary. Results indicate two primary modifications reduce eutrophication. The first clears a constricted channel to the interior of the system and promotes flushing. The second fills deep holes in which anoxic conditions promote sediment nutrient release to the water column. Major reductions in pathogen concentration require regional controls on sources.