Cost-effective sampling network design for contaminant plume monitoring under general hydrogeological conditions

A new simulation-optimization methodology is developed for cost-effective sampling network design associated with long-term monitoring of large-scale contaminant plumes. The new methodology is similar in concept to the one presented by Reed et al. (Reed, P.M., Minsker, B.S., Valocchi, A.J., 2000a. Cost-effective long-term groundwater monitoring design using a genetic algorithm and global mass interpolation. Water Resour. Res. 36 (12), 3731-3741) in that an optimization model based on a genetic algorithm is coupled with a flow and transport simulator and a global mass estimator to search for optimal sampling strategies. However, this study introduces the first and second moments of a three-dimensional contaminant plume as new constraints in the optimization formulation, and demonstrates the proposed methodology through a real-world application. The new moment constraints significantly increase the accuracy of the plume interpolated from the sampled data relative to the plume simulated by the transport model. The plume interpolation approaches employed in this study are ordinary kriging (OK) and inverse distance weighting (IDW). The proposed methodology is applied to the monitoring of plume evolution during a pump-and-treat operation at a large field site. It is shown that potential cost savings up to 65.6% may be achieved without any significant loss of accuracy in mass and moment estimations. The IDW-based interpolation method is computationally more efficient than the OK-based method and results in more potential cost savings. However, the OK-based method leads to more accurate mass and moment estimations. A comparison of the sampling designs obtained with and without the moment constraints points to their importance in ensuring a robust long-term monitoring design that is both cost-effective and accurate in mass and moment estimations. Additional analysis demonstrates the sensitivity of the optimal sampling design to the various coefficients included in the objective function of the optimization model.     

Application of a novel modeling tool with multistressor functionality to support management of organic contaminants in the Baltic Sea

Organic contaminants constitute one of many stressors that deteriorate the ecological status of the Baltic Sea. When managing environmental problems in this marine environment, it may be necessary to consider the interactions between various stressors to ensure that averting one problem does not exacerbate another. A novel modeling tool, BALTSEM-POP, is presented here that simulates interactions between climate forcing, hydrodynamic conditions, and water exchange, biogeochemical cycling, and organic contaminant transport and fate in the Baltic Sea. We discuss opportunities to use the model to support different aspects of chemicals management. We exemplify these opportunities with a case study where two emission-reduction strategies for a chemical used in personal care products (decamethylcyclopentasiloxane) are evaluated, and where the confounding influence of future climate change and eutrophication on the impact of the emission-reduction strategies are assessed.     

Migration of Cs through quartz sand: experimental results and modeling approaches

We present results from experiments on the migration of ¹³⁷Cs through columns containing quartz sand. Times for ¹³⁷Cs movement through these columns and the quantity of ¹³⁷Cs adsorbed by the sand decreased as the ionic strength of the pore water increased from 0.002 to 0.1 m. The breakthrough curves were characterized by a slow approach towards steady-state concentrations as well as by long tails, indicating that ¹³⁷Cs adsorption to the sand grains was, at least in part, controlled by rate-limited reactions. Various formulations for solute mass transfer were tested for their ability to fit the experimental breakthrough curves. Based on a statistical analysis, a nonlinear, two-site model was identified as the most appropriate for describing the suite of experimental data. Variation in the model parameter that describes the rate of ¹³⁷Cs adsorption to the sand showed no consistent pattern with changes in ionic strength. In contrast, model parameters describing the sorption capacity of the sand grains and the fraction of kinetic sorption sites on the sand decreased with increasing ionic strength. The parameter describing the rate of ¹³⁷Cs desorption varied directly with changes in ionic strength.     

Model-based experimental design for assessing effects of mixtures of chemicals

We exposed flour beetles (Tribolium castaneum) to a mixture of four poly aromatic hydrocarbons (PAHs). The experimental setup was chosen such that the emphasis was on assessing partial effects. We interpreted the effects of the mixture by a process-based model, with a threshold concentration for effects on survival. The behavior of the threshold concentration was one of the key features of this research. We showed that the threshold concentration is shared by toxicants with the same mode of action, which gives a mechanistic explanation for the observation that toxic effects in mixtures may occur in concentration ranges where the individual components do not show effects. Our approach gives reliable predictions of partial effects on survival and allows for a reduction of experimental effort in assessing effects of mixtures, extrapolations to other mixtures, other points in time, or in a wider perspective to other organisms.     

Cost-effective treatment of sludge conditioning using supernatant fluid polyelectrolyte

Environmental Science and Pollution Research - Treatment of activated sludge using pure polyelectrolyte and its supernatant fluids for sludge dewatering was examined. The supernatant was used for...     

Cost-effective reduction of NOx emissions from electricity generation.

This paper analyzes the benefits and costs of policies to reduce NOx emissions from electricity generation in the United States. Because emissions of NO contribute to the high concentration of atmospheric ozone in the eastern states associated with health hazards, the U.S. Environmental Protection Agency (EPA) has called on eastern states to formulate state implementation plans (SIPs) for reducing NOx emissions. Our analysis considers three NOx reduction scenarios: a summer seasonal cap in the eastern states covered by EPA's NOx SIP Call, an annual cap in the same SIP Call region, and a national annual cap. All scenarios allow for emissions trading. Although EPA's current policy is to implement a seasonal cap in the SIP Call region, this analysis indicates that an annual cap in the SIP Call region would yield about $400 million more in net benefits (benefits less costs) than would a seasonal policy, based on particulate-related health effects only. An annual cap in the SIP Call region is also the policy that is most likely to achieve benefits in excess of costs. Consideration of omissions from this accounting, including the potential benefits from reductions in ozone concentrations, strengthens the finding that an annual program offers greater net benefits than does a seasonal program.     

Impact of injection system design on ISCO performance with permanganate--mathematical modeling results

In situ chemical oxidation (ISCO) using permanganate (MnO(4)(-)) can be a very effective technique for remediation of soil and groundwater contaminated with chlorinated solvents. However, many ISCO projects are less effective than desired because of poor delivery of the chemical reagents to the treatment zone. In this work, the numerical model RT3D was modified and applied to evaluate the effect of aquifer characteristics and injection system design on contact and treatment efficiency. MnO(4)(-) consumption was simulated assuming the natural oxidant demand (NOD) is composed of a fraction that reacts instantaneously and a fraction that slowly reacts following a 2nd order relationship where NOD consumption rate increases with increasing MnO(4)(-) concentration. MnO(4)(-) consumption by the contaminant was simulated as an instantaneous reaction. Simulation results indicate that the mass of permanganate and volume of water injected has the greatest impact on aquifer contact efficiency and contaminant treatment efficiency. Several small injection events are not expected to increase contact efficiency compared to a single large injection event, and can increase the amount of un-reacted MnO(4)(-) released down-gradient. High groundwater flow velocities can increase the fraction of aquifer contacted. Initial contaminant concentration and contaminant retardation factor have only a minor impact on volume contact efficiency. Aquifer heterogeneity can have both positive and negative impacts on remediation system performance, depending on the injection system design.     

Adsorption study for the removal of a basic dye: experimental and modeling

An effective adsorbent is developed from saw dust and its various adsorption characteristics are studied for removing a basic dye (crystal violet) from its aqueous solution. Equilibrium data are fitted to various adsorption isotherms. It is seen that about 341 mg of crystal violet can be removed using 1g of the adsorbent at 298 K. Kinetic study is also carried out to observe the effects of various process parameters viz. particle size of the adsorbent, initial concentration of the dye, temperature and adsorbent amount. A generalized two-resistance mass transfer model, which includes a film mass transfer coefficient (k(f)) and an internal effective diffusivity (Dp), is used to interpret the adsorption kinetic data. The model parameters (k(f) and Dp) are estimated by fitting the experimental data to the model. The evaluated parameters are used to predict the concentration profiles at various other operating conditions. The average deviation of the predicted values lies within 10% in all the cases. Sensitivity analysis is performed to observe the sensitivity of the model to the variations in the model parameters.     

Spatial modeling for air pollution monitoring network design: example of residential woodsmoke

The purpose of this paper is to demonstrate how to develop an air pollution monitoring network to characterize small-area spatial contrasts in ambient air pollution concentrations. Using residential woodburning emissions as our case study, this paper reports on the first three stages of a four-stage protocol to measure, estimate, and validate ambient residential woodsmoke emissions in Vancouver, British Columbia. The first step is to develop an initial winter nighttime woodsmoke emissions surface using inverse-distance weighting of emissions information from consumer woodburning surveys and property assessment data. Second, fireplace density and a compound topographic index based on hydrological flow regimes are used to enhance the emissions surface. Third, the spatial variation of the surface is used in a location-allocation algorithm to design a network of samplers for the woodsmoke tracer compound levoglucosan and fine particulate matter. Measurements at these network sites are then used in the fourth stage of the protocol (not presented here): a mobile sampling campaign aimed at developing a high-resolution surface of woodsmoke concentrations for exposure assignment in health effects studies. Overall the results show that relatively simple data inputs and spatial analysis can be effective in capturing the spatial variability of ambient air pollution emissions and concentrations.     

Optimal design of a helical coil support for dewars in fuel cell applications

Environmental Science and Pollution Research - Fuel cells are gaining popularity because of their efficient energy production without causing environmental pollution. Recently, DRDO has developed a...     

Biogeochemical influence on transport of chromium in manganese sediments: experimental and modeling approaches

Hexavalent chromium (Cr(VI)) was reduced to immobile and nontoxic Cr(III) by a dissimilatory metal reducing bacteria, Shewanella alga Simidu (BrY-MT) ATCC 55627. A series of kinetic batch and dynamic column experiments were conducted to provide an understanding of Cr(VI) reduction by the facultative anaerobe BrY-MT. Reduction of Cr(VI) was rapid (within 1 h) in columns packed with quartz sand and bacteria, whereas Cr(VI) reduction by BrY-MT was delayed (57 h) in the presence of beta-MnO2-coated sand. A mathematical model was developed and evaluated against data obtained from column experiments. The model takes into account (1) advective-dispersive transport of Cr(III), Cr(VI), lactate, and protein (mobile and immobile bacteria); (2) first-order kinetic adsorption of Cr(III) and lactate; (3) conversion of solid phase beta-MnO2 to solid phase MnOOH due to oxidation of Cr(III); (4) dual-Monod kinetics, where Cr(VI) is the electron acceptor and lactate is the electron donor. The breakthrough data for Cr(III), Cr(VI), lactate, and protein (mobile and immobile bacteria) were fitted simultaneously. The breakthrough data are well described by the mathematical model that considers the above processes. This result demonstrates the ability of the coupled hydrobiogeochemical model to simulate chromium transport in complex reactive systems.     

Field verification of sound attenuation modeling and air emission testing in support of missile motor disposal activities

The U.S. Department of Defense-approved activities conducted at the Utah Test and Training Range (UTTR) include both operational readiness test firing of intercontinental ballistic missile (ICBM) motors, as well as the destruction of obsolete or otherwise unusable ICBM motors through open burn/open detonation (OB/OD). Within the Utah Division of Air Quality, these activities have been identified as having the potential to generate unacceptable noise levels, as well as significant amounts of volatile organic compounds (VOCs). Hill Air Force Base, UT, has completed a series of field tests at the UTTR in which sound-monitoring surveillance of OB/OD activities was conducted to validate the Sound Intensity Prediction System (SIPS) model. Using results generated by the SIPS model to support the decision to detonate, the UTTR successfully disposed of missile motors having an aggregate net explosive weight (NEW) of 81,374 lb without generating adverse noise levels within populated areas. In conjunction with collecting noise-monitoring data, air emissions were collected to support the development of air emission factors for both static missile motor firings and OB/OD activities. Through the installation of 15 ground-based air samplers, the generation of combustion-fixed gases, VOCs, and chlorides was monitored during the 81,374-lb NEW detonation event. Comparison of field measurements to predictions generated from the US Navy energetic combustion pollutant formation model, POLU4WN, indicated that, as the detonation fire ball expanded, organic compounds, as well as CO, continued to oxidize as the combustion gases mixed with ambient air. VOC analysis of air samplers confirmed the presence of chloromethane, vinyl chloride, benzene, toluene, and 2-methyl-1-propene. Qualitative chloride analysis indicated that gaseous HCl was generated at low concentrations, if at all.     

Investigation of various physicochemical and environmental parameter influence on pesticide sorption to ditch bed substratum by means of experimental design

Diffuse pollution by pesticide applied in rural catchments may contribute to alter water quality. Besides actions relative to the way the substances are introduced into the environment, it is also possible to limit the contamination by interfering on their transfer pathways from fields to the main river network. Especially, interface areas such as buffer strips or small ditches may play a major part in pesticide diffuse pollution decrease. In ditches a great variety of materials may act as sorbents for organic contaminants: grass, leaves, wood debris or sediments. In this study, laboratory experiments were designed to determine sorption characteristics for three herbicides with different physicochemical properties on sediment and leaves in decay commonly found in agricultural ditches. Sorption capacities were assessed for the herbicides isoproturon, diuron and diflufenican.

Experimental design was carried out to investigate the effects of five parameters on herbicide sorption on sediment and dead leaves. These parameters have been chosen according to parallel field experiment needs. Thus, the influence of initial sorbent moisture, herbicide form, i.e. active substance or commercial formulation, water quality (tap or natural ditch water), bromide ions (used as conservative tracers) and solid/liquid ratio have been tested. Within the parameters investigated, pesticide formulation and solid/liquid ratio were the most important parameters affecting pesticide sorption on both ditch materials.     

Removal of carbon disulfide via coupled reactions on a bi-functional catalyst: experimental and modeling results

A mathematical model describing the rate of carbon disulfide (CS₂) removal due to coupled reactions has been developed. Kinetic studies were carried out in a fixed bed reactor under atmospheric pressure and a range of temperatures (85–125 °C). The effects of flow rate, CS₂ inlet concentration, temperature and relative humidity were analyzed. A kinetic model based on axial dispersion, external and internal mass-transfer resistances, as well as effects of S deposition on the inner-face of the catalyst was in agreement with the CS₂ experimental breakthrough curves. The mathematical model can be used for process design and scale up of similar systems.     

Simulation modeling for nitrogen removal and experimental estimation of mass fractions of microbial groups in single-sludge system

Nitrification-denitrification in a single-sludge nitrogen removal system (SSNRS; with a sufficient carbon source for denitrification) was performed. With an increase in the mixed liquor recycle ratio (R(m)) from 1 to 2, the total nitrogen (TN) removal efficiency at a lower volumetric loading rate (VLR=0.21 NH(4)(+)-N m(-3) d(-1)) increased, but the TN removal efficiency at a higher VLR (0.35 kg NH(4)(+)-N m(-3) d(-1)) decreased. A kinetic model that accounts for the mass fractions of Nitrosomonas, Nitrobacter, nitrate reducer and nitrite reducer (f(n1), f(n2), f(dn1), and f(dn2)) in the SSNRS and an experimental approach for the estimation of the mass fractions of nitrogen-related microbial groups are also proposed. The estimated f(dn1) plus f(dn2) (0.65-0.83) was significantly larger than the f(n1) plus f(n2) (0.28-0.32); the f(n1) (0.21-0.26) was larger than the f(n2) (0.05-0.07); and the f(dn1) (0.32-0.45) varied slightly with the f(dn2) (0.33-0.38). At the lower VLR, the f(dn1) plus f(dn2) increased with increasing R(m); however at the higher VLR, the f(dn1) plus f(dn2) did not increase with increasing R(m). By using the kinetic model, the calculated residual NH(4)(+)-N and NO(2)(-)-N in the anoxic reactor and NO(2)(-)-N and NO(3)(-)-N in the aerobic reactor were in fairly good agreement with the experimental data; the calculated NO(3)(-)-N in the anoxic reactor was over-estimated and the calculated NH(4)(+)-N in the aerobic reactor was under-estimated.     

Using SFA indicators to support environmental policy

In order to improve the link from Substance Flow Analysis (SFA) studies to environmental policy, a translation is made from the SFA overview of flows and stocks into a limited set of indicators. This set is designed to evaluate a region’s substance regime with regard to environmental quality and sustainable development, including problem shifting in time and space.     

Bayesian hierarchical modeling of personal exposure to particulate matter

In the US EPA's 1998 Baltimore Epidemiology-Exposure Panel Study, a group of 16 residents of a single building retirement community wore personal monitors recording personal fine particulate air pollution concentrations (PM_2.5) for 27 days, while other monitors recorded concurrent apartment, central indoor, outdoor and ambient site PM_2.5 concentrations. Using the Baltimore panel study data, we develop a Bayesian hierarchical model to characterize the relationship between personal exposure and concentrations of PM_2.5 indoors and outdoors. Personal exposure is expressed as a linear combination of time spent in microenvironments and associated microenvironmental concentrations. The model incorporates all available monitoring data and accounts for missing data and sources of uncertainty such as measurement error and individual differences in exposure. We discuss the implications of using personal versus ambient PM_2.5 measurements in characterization of personal exposure to PM_2.5.     

Application of dispersion modeling to indoor gas release scenarios

Many complex models are available to study the dispersion of contaminants or ventilation effectiveness in indoor spaces. Because of the computationally complex numerical schemes employed, most of these models require mainframe computers or workstations. However, simple design tools or guidelines are needed, in addition to complicated models. A dispersion model based on the basic governing equations was developed and uses an analytical solution. Because the concentration is expressed by an analytical solution, the grid size and time steps are user definable. A computer program was used to obtain numerical results and to obtain release history from a thermodynamic source model. The model can be used to estimate three-dimensional spatial and temporal variations in concentrations resulting from transient gas releases in an enclosure. The model was used to study a gas release scenario from a pressurized cylinder into a large ventilated building, in this case, a transit parking and fueling facility.     

An experimental design to probe the interactions of dissolved organic matter and xenobiotics: bioavailability of pyrene and 2,2',5,5'-tetrachlorobiphenyl to Daphnia magna

Experiments were conducted to probe the interactions between natural dissolved organic matter (DOM) and two xenobiotics, and to determine how DOM influences their bioavailability. The experimental set-up, using dialysis bags, was designed to expose test organisms to the same constant concentration of free dissolved chemical, while increasing the concentration of the bound-to-DOM fraction. Daphnia magna S. were exposed to pyrene or 2,2',5,5'-tetrachlorobiphenyl in the presence of 0, 1, 2, 5, 10 or 20 mg L(-1) of a reference riverine humic acid (Suwannee River Humic Acid). The physico-chemical parameters were well constrained in the microcosm, demonstrating its potential usefulness. However bioaccumulation by D. magna showed important variability between replicate treatments, sufficient to mask any trends as a function of DOM concentration. The organic-carbon-normalised partition coefficients (K(OC)) ranged from 52000 to 92000 L kg(-1) for pyrene and from 8200 to 89000 L kg(-1) for 2,2',5,5'-tetrachlorobiphenyl, with a marked "concentration effect" for the latter compound.     

Inverse modeling to estimate LNAPL plume release timing

Methodologies are presented for dating releases of light nonaqueous phase liquids (LNAPLs) using an inverse modeling approach with simple analytical models. Models for LNAPL plume migration are presented to predict LNAPL plume velocity in the unsaturated and saturated zones as a function of basic soil and fluid properties. A relative mobility factor is introduced for LNAPL movement at the water table that depends primarily on the van Genuchten n parameter (related to the breadth of the soil pore size distribution) and the magnitude of water table fluctuations. Estimated LNAPL plume velocities compare reasonably with more rigorous numerical models, which may be used in cases where data availability warrant the greater effort entailed.Two methods of estimating release timing and its uncertainty are investigated. A direct estimation method is described that determines travel time for a single observed travel distance based on estimated soil and fluid properties. Release date uncertainty may be determined using the first order (FO) or Monte Carlo (MC) methods. The second method for estimating release date involves nonlinear parameter estimation utilizing distance vs. time measurements and other data.A case study is presented for a field site where independent estimates of release timing were obtained from a numerical modeling analysis. Release timing estimates based on direct inversion of the analytical timing model agree well with the numerical analysis. Results for a second field site indicate that release date confidence limits estimated by the FO method, assuming log-normally distributed travel times, are close to values determined by the MC method, which makes no assumption regarding the form of the travel time probability distribution.Results for a hypothetical problem indicate that LNAPL velocity and travel time may be accurately estimated if sufficient data on travel distance vs. time are available. Incorporating prior information on relevant soil and fluid properties into the objective function reduces the uncertainty in release date if prior estimates are accurate. However, biased prior estimates may lead to over- or underestimation of release date uncertainty. Simultaneous estimation of soil and fluid properties and release date is possible if prior information is available to condition the parameter estimates.