Optimal Wasteload Allocation Procedure for Achieving Dissolved Oxygen Water Quality Objectives. I: Sensitivity Analysis

A method is presented to compute sensitivities of in-stream dissolved oxygen (DO) with respect to perturbations in the load vector and the reaction coefficients that make up the eutrophication cycle. It is shown that the direct sensitivity method, i.e., the repetitive solution of the direct problem, produces the desired information, however at a large computational cost. The utilization of the adjoint sensitivity method proves to be a much more efficient way to compute these sensitivities as large subsets of the sensitivity information domain can be easily extracted with just a few runs. It is found that for the given problem setup, in-stream DO is most sensitive to ammonia, effluent DO, algae, and carbonaceous biochemical oxygen demand (CBOD) loads. Additionally, the computed sensitivities vary considerably in their general trend over the simulation time. Sensitivities also are computed with respect to the reaction coefficients (26 total) that govern the interdependency of all constituents. Sediment oxygen demand proves to be the coefficient with the highest influence that is three orders of magnitude higher than the next set of coefficients comprised of reaeration, CBOD degradation, nitrification, and denitrification. All other coefficients have a negligible influence on DO concentrations. Computations are carried out using a two-dimensional model formulation applied to a long rectangular channel with varying width and slope and periodic but unsteady flow conditions.     

Optimal Wasteload Allocation Procedure for Achieving Dissolved Oxygen Water Quality Objectives. II: Optimal Load Control

This paper presents the development of an efficient strategy for achieving in-stream dissolved oxygen (DO) water quality standards (WQSs) via optimized point-load control strategies using the adjoint method. To this end, a least-squares-type objective function is formulated that measures the difference between desired (WQSs) and current DO concentrations at strategically selected monitoring points in the domain. The goal is to minimize the difference between actual DO concentration and the WQS, hence allowing time-variant loadings to utilize the assimilative capacity of the receiving water body at an optimal level. Time-variant discharge rates for a number of discharge locations are considered as control parameters, while different zone-specific critical DO levels are imposed as constraints. The selection of the control is kept flexible and a number of different scenarios are tested. First, only carbonaceous biochemical oxygen demand is used, which allows for a reduction of the number of equations that need to be solved. In the other tests, all constituents are switched on and different variables at each load node are selected as a control by first varying the concentrations individually, and then linking them through control of the volumetric flow rate. Optimization is achieved using a conjugate gradient search method, for which the gradients are computed through the solution of both the direct and adjoint problems. It is shown that the large amount of gradient information (parameter space has a dimension of several thousands) can be computed very efficiently using the adjoint, and that optimized results are achieved after only a few iterations irrespective of the initial guess. Computations are carried out using both two-dimensional model formulation applied to a long rectangular channel with varying width and slope and a model for the upper Potomac River estuary.     

High-Frequency Diel Dissolved Oxygen Stream Data Modeled for Variable Temperature and Scale

Diel dissolved oxygen (DO) concentrations and temperature were sensed at high-frequency and modeled in an eastern Iowan stream, Clear Creek, in an agricultural setting. The magnitude of the diel changes in DO and temperature were largest at the upstream (headwater) station. Inclusion of temperature change factors increased the accuracy of modeling results and yielded estimates of the reaeration rate constant, primary production rate, and respiration rate. The DO modeling of the high-frequency measurements (15-min intervals) revealed a temperature-driven nonlinear reaeration process that led to increases in nighttime DO concentrations. The DO modeling results from three sensing stations in the watershed revealed decreasing trends in primary productivity, respiration, and the reaeration rate constant with increasing drainage area. Light extinction from suspended solids was the main factor limiting net primary production. As a result, the P/R ratio also decreased with increasing drainage area. High-frequency sensor data and DO modeling revealed the effects of temperature and watershed scale on the primary factors that dictate diel DO dynamics in a stream setting.     

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.     

Assessing the Impacts of Nutrient Load Uncertainties on Predicted Truckee River Water Quality

This study examined the effects of uncertain model boundary conditions on dissolved oxygen (DO) predictions for the lower Truckee River, Nevada using an augmented version of the EPA’s Water Quality Analysis Simulation Program Version 5 (WASP5) that included periphyton, or attached algae, in eutrophication kinetics. Uncertainty analyses were performed on selected organic nitrogen (ON) and carbonaceous biochemical oxygen demand boundary conditions using Monte Carlo techniques. The stochastic model was run using boundary concentrations assigned from observed probability distributions. Ranges of simulated values were used to construct confidence intervals, the magnitudes of which indicated the uncertainty associated with model predictions. Uncertainty in agricultural ditch return concentrations had minimal effects on in-stream model predictions, as predicted values of daily minimum and maximum DOs, daily average ON, and periphyton biomass all failed to show significant variability as a result of ditch concentration uncertainty. This result indicates that while ditch return nutrient loads are not trivial, their exact concentrations are not needed to make relatively accurate predictions of in-stream DO. However, uncertainty in the upstream ON boundary did result in significant uncertainty during summer months with regard to in-stream model predictions of ON, periphyton biomass, and DO. The model is clearly more sensitive to changes in this boundary than to changes in agricultural ditch concentrations.     

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.     

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.     

Parameter Sensitivity and Predictive Uncertainty in a New Water Quality Model, Q2

A new dynamic model of water quality, Q2, has recently been developed, capable of simulating large branched river systems. This paper describes the application of a generalized sensitivity analysis (GSA) to Q2 for single reaches of the River Thames in southern England. Focusing on the simulation of dissolved oxygen (DO) (since this may be regarded as a proxy for the overall health of a river); the GSA is used to identify key parameters controlling model behavior and provide a probabilistic procedure for model calibration. It is shown that, in the River Thames at least, it is more important to obtain high quality forcing functions than to obtain improved parameter estimates once approximate values have been estimated. Furthermore, there is a need to ensure reasonable simulation of a range of water quality determinands, since a focus only on DO increases predictive uncertainty in the DO simulations. The Q2 model has been applied here to the River Thames, but it has a broad utility for evaluating other systems in Europe and around the world.     

Incorporating Both Physical and Kinetic Limitations in Quantifying Dissolved Oxygen Flux to Aquatic Sediments

Traditionally, dissolved oxygen (DO) fluxes have been calculated using the thin-film theory with DO microstructure data in systems characterized by fine sediments and low velocities. However, recent experimental evidence of fluctuating DO concentrations near the sediment-water interface suggests that turbulence and coherent motions control the mass transfer, and the surface renewal theory gives a more mechanistic model for quantifying fluxes. Both models involve quantifying the mass transfer coefficient (k) and the relevant concentration difference (ΔC). This study compared several empirical models for quantifying k based on both thin-film and surface renewal theories, as well as presents a new method for quantifying ΔC (dynamic approach) that is consistent with the observed DO concentration fluctuations near the interface. Data were used from a series of flume experiments that includes both physical and kinetic uptake limitations of the flux. Results indicated that methods for quantifying k and ΔC using the surface renewal theory better estimated the DO flux across a range of fluid-flow conditions.     

Eutrophication Model for the Patuxent Estuary: Advances in Predictive Capabilities

A robust eutrophication and sediment diagenesis model has been developed for the Patuxent Estuary to study the impact of different nutrient loadings on phytoplankton biomass and dissolved oxygen (DO) levels. The modeling approach was to begin with an existing water quality model (CE-QUAL-W2) for the Patuxent Estuary (hereafter referred to as the Estuary). First, formulations for the water column kinetics were completely replaced with routines based on the WASP/EUTRO5 water quality model. Then, a sediment diagenesis component was added to simulate the accumulation and mineralization of organic matter in the sediment, the generation of sediment oxygen demand, and the flux of phosphate and ammonia from the sediment. Loadings from the tributaries for nutrients and flow were based on a combination of watershed modeling and sampling by scientists at the Smithsonian Environmental Research Center. The new model was able to reproduce the ambient water quality data from 1997 to 1999 by adequately simulating the high concentrations of phytoplankton and low DO levels in the Estuary. The model was then used to evaluate the response to various hypothetical nutrient loading scenarios. Model results show that phytoplankton growth in the upper Estuary is much more sensitive to nutrient loading from tributaries than in the lower estuary. Further, model results indicate that DO concentrations in the lower Estuary are largely influenced by levels of nutrients and organic carbon at the mouth of the Estuary.     

Clustered distribution of calcium sensitivities: an indication of hetero-tetrameric gating components in Ca2+-activated K+ channels reconstituted from avian nasal gland cells

Calcium-activated potassium channels (maxi K+ channels) isolated from avian nasal salt gland cells were reconstituted into lipid bilayers and characterized. The 266 pS channel is blocked discretely by charybdotoxin from the external solution at nanomolar concentrations and by Ba2+ from the cytosolic side at micromolar concentrations. Fast tetraethylammonium (TEA) block is seen as apparent reductions in amplitude of the unitary currents. From the extent of the reductions, TEA binding affinity was calculated to be 0.16 mM from the external solution and 37 mm from internal solution. The overall channel properties conform to those of maxi K+ channels in other epithelial tissues. The calcium sensitivity of the channel was found to be variable from channel to channel, extending over a wide range of concentrations from 1 to 1,000 microM. Examination of the pooled calcium titration curves, revealed that these curves are grouped into five clusters, and the probability distribution of the clusters matches a binomial distribution. The Hill coefficient derived from the titration curves varies from 1 to 5 and is linearly correlated to calcium binding with a slope of 1 per 10-fold change in Kd. Clustered titration curves with such a characteristic suggest that the gating components and the calcium binding sites of the maxi K+ channels in the avian nasal gland are hetero-tetrameric and may result from random mixing of two distinct subunits possessing high and low calcium sensitivities, respectively.     

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.     

Depth-Integrated Estimation of Dissolved Oxygen in a Lake

The majority of variable estimation studies in water resources investigate the temporal variation of the variable. In this study, we examined the depth-dependent estimation of a lake’s dissolved oxygen (DO) using two artificial neural network (ANN) methods: (1) the radial basis functions (RBFs) and the feed forward back-propagation (FFBP), and (2) the multilinear regression (MLR). We tested two different input layer configurations. In the first case, we employed all other available lake parameters—total dissolved solids (TDS), pH, conductivity, lake depth, and lake temperature—to estimate DO. In the second case, we considered only depth and temperature to estimate DO. The performance evaluation criteria of these two cases were close. ANN estimation performances were noticeably superior to those of MLR, as reflected in the performance evaluation criteria and DO lake depth plots. We saw that the spatial variation of the lake’s DO can be captured by ANNs satisfactorily, even if available measurements are quite limited.     

Density Intrusion and Variation in Dissolved Oxygen Concentrations in a Bay with a Sill at Its Mouth

This paper describes the influence of water outside a bay on water quality within the bay for a typical enclosed bay with a sill at its mouth. The influence of water outside Ohfunato Bay, Japan was expressed in a one-dimensional vertical dissolved oxygen (DO) model as the rate of seawater exchange caused by tides and the density difference between water inside and outside the bay. A nonhydrostatic, three-dimensional model was used to confirm the occurrence of density intrusions. A one-dimensional vertical DO model was used to examine long-term changes in DO concentrations in the bottom layer of the bay. Results indicated that seawater exchange caused by density differences had a significant effect on the formation and disappearance of anoxic water in the bay’s bottom layer.     

Safe hemoglobin or hematocrit levels in surgical patients

The terminology and fundamental aspects of the delivery, consumption, and deficits of oxygen are recalled. In chronic and acute, nonseptic states, red blood cell (RBC) transfusion is capable of increasing oxygen consumption (VO2). In acute septic states, the response of VO2 to RBC transfusion is variable and unpredictable, but attempts to increase oxygen delivery (DO2) should be made if the clinical picture raises the suspicion of a potentially lethal oxygen deficit. Therapeutic interventions raising the cardiac index to "supranormal" values in critically ill patients improve their chances of survival; and maintenance of hemoglobin or hematocrit values around 11 g/dl or 33%, respectively, is one part of such interventions. Opinions differ on the general tolerance of anemia, as witnessed by postulated "critical levels" of the hemoglobin concentration between approximately 11 and 4 to 5 g/dl or hematocrit values between 33% and 12% to 15%, respectively. The common denominator underlying these vastly different opinions is the variable behavior of several "non-Hb variables," which influence the venous oxygen tensions apart from the hemoglobin or hematocrit. Abnormalities of these non-Hb variables-typically encountered in the critically ill-increase the dependence of patients on hemoglobin or hematocrit levels that suffice to protect them against an oxygen deficit. For this reason, the "critical" hemoglobin or hematocrit is an individual value, and a generally valid "transfusion trigger" does not exist. Finally, the entity now known as silent myocardial ischemia (SMI) is a decisive factor for the tolerance of anemia. Solid clinical evidence is now available to support the concept that patients over age 40 should not, as an elective procedure, be subjected to levels < 10 g/dl or < 30%, respectively, without prior exclusion of SMI by appropriate investigations.     

Critical Analysis of Computing Equations for Determination of BOD by Dilution Methods

The dilution method has been accepted widely as the standard procedure for determining biochemical oxygen demand (BOD). A mass-balance approach was used to analyze four equations for computing BOD: two involving external seed correction methods (Standard Methods) and two involving internal correction methods (Sawyer and McCarty in 1978 and Kline and Gibbs in 1979). Five testing conditions were identified as being pertinent to the underlying assumptions: seed addition method, uniform BOD bottle volume, uniform initial dissolved oxygen (DO), seed uptake linearity, and data handling method. The results of this study showed that internal methods may provide more reproducible results, primarily due to the inherent correction for dilution water demand and better data handling (regression) technique to cope with the variation in BOD bottle volume and initial DO. However, internal methods may violate the requirements imposed by Standard Methods for minimum final DO (1?mg/L) for the high spike runs and minimum 5?day DO depletion (2?mg/L) for the low spike runs. Recommendations are also made to improve the reliability of these methods.     

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.     

高压下的碳——氧平衡

计算和实测了１８７３Ｋ，气相总压力０．１，０．２，０．３ＭＰａ时的碳－氧平衡曲线，并分析了两者偏差的原因，研究表明，随着气相总压的提高，碳－氧平衡曲线上移，这可能很好的解释密封转炉的炉腔高压下吹氧脱碳的强化。     

The kinetics of disperse order development in α-cual alloys

The kinetic parameters of the Mehl-Johnson-Avrami rate equation for disperse order (DO) formation in a-CuAl alloy were determined from nonisothermal experiments. Activation energy values of 145 ± 7/155 ± 10/178 ± 14 kJ/mol were obtained from rate dependent curves for alloys containing, respectively, 19, 13, and 6.5 pct Al. These values are consistent with those for self-diffusion. Further analysis allowed for the determination of the preexponential factorko and the exponentn (= 1.39 ± 0.2), the latter being found to be insensitive to alloy composition. This indicates a diffusion limited growth of particles from a nonnegligible size, suggesting that the preexisting short range order acts as nuclei for DO development. By employing the absolute reaction rate theory, a virtually constant activation free energy, ΔG* = 163 ± 15 kJ/mol, was measured for the three alloys investigated. These results indicate that the transformation proceeds more rapidly with increasing aluminum content of the alloy. Values for domain size and concentration, computed by the classical theory for nucleation and growth processes, are compared with those obtained experimentally. Also, an (α + DO) field is proposed in the Cu-Al equilibrium diagram.     

Genetic analysis of adult-onset cataract in a city-based ophthalmic hospital

When oxygen delivery (DO2) critically decreases, oxygen consumption (VO2) becomes supply dependent. We examined whether end-tidal PCO2 (PetCO2) would identify supply dependency during shock. Five dogs (Group I) underwent progressive hemorrhage to decrease DO2 until they could no longer maintain a stable blood pressure. Five additional animals (Group II) were bled until VO2 decreased to 70% of baseline, followed by resuscitation. The PetCO2 versus time inflection point was compared with the DO2 at onset of supply dependency (DO2crit). DO2crit for Groups I and II were 6.9 +/- .4 and 8.1 +/- 1.3, respectively (p = NS), and not statistically different from the DO2 values at which PetCO2 decreased (6.6 +/- .7 and 6.3 +/- .7 mL/kg per min, respectively). AT constant minute volume, PetCO2 effectively indicated the onset of supply dependency and rapidly increased during resuscitation, paralleling the changes in VO2 in this model of hemorrhagic shock.