U(VI) bioreduction with emulsified vegetable oil as the electron donor--microcosm tests and model development.

We conducted microcosm tests and biogeochemical modeling to study U(VI) reduction in contaminated sediments amended with emulsified vegetable oil (EVO). Indigenous microorganisms in the sediments degraded EVO and stimulated Fe(III), U(VI), and sulfate reduction, and methanogenesis. Acetate concentration peaked in 100-120 days in the EVO microcosms versus 10-20 days in the oleate microcosms, suggesting that triglyceride hydrolysis was a rate-limiting step in EVO degradation and subsequent reactions. Acetate persisted 50 days longer in oleate- and EVO- than in ethanol-amended microcosms, indicating that acetate-utilizing methanogenesis was slower in the oleate and EVO than ethanol microcosms. We developed a comprehensive biogeochemical model to couple EVO hydrolysis, production, and oxidation of long-chain fatty acids (LCFA), glycerol, acetate, and hydrogen, reduction of Fe(III), U(VI) and sulfate, and methanogenesis with growth and decay of multiple functional microbial groups. By estimating EVO, LCFA, and glycerol degradation rate coefficients, and introducing a 100 day lag time for acetoclastic methanogenesis for oleate and EVO microcosms, the model approximately matched observed sulfate, U(VI), and acetate concentrations. Our results confirmed that EVO could stimulate U(VI) bioreduction in sediments and the slow EVO hydrolysis and acetate-utilizing methanogens growth could contribute to longer term bioreduction than simple substrates (e.g., ethanol, acetate, etc.) in the subsurface.

U(VI) bioreduction with emulsified vegetable oil as the electron donor--model application to a field test.

We amended a shallow fast-flowing uranium (U) contaminated aquifer with emulsified vegetable oil (EVO) and subsequently monitored the biogeochemical responses for over a year. Using a biogeochemical model developed in a companion article (Tang et al., Environ. Sci. Technol.2013, doi: 10.1021/es304641b) based on microcosm tests, we simulated geochemical and microbial dynamics in the field test during and after the 2-h EVO injection. When the lab-determined parameters were applied in the field-scale simulation, the estimated rate coefficient for EVO hydrolysis in the field was about 1 order of magnitude greater than that in the microcosms. Model results suggested that precipitation of long-chain fatty acids, produced from EVO hydrolysis, with Ca in the aquifer created a secondary long-term electron donor source. The model predicted substantial accumulation of denitrifying and sulfate-reducing bacteria, and U(IV) precipitates. The accumulation was greatest near the injection wells and along the lateral boundaries of the treatment zone where electron donors mixed with electron acceptors in the groundwater. While electron acceptors such as sulfate were generally considered to compete with U(VI) for electrons, this work highlighted their role in providing electron acceptors for microorganisms to degrade complex substrates thereby enhancing U(VI) reduction and immobilization.

Improved Method for Estimating Reaction Rates During Push-Pull Tests.

The breakthrough curve obtained from a single-well push-pull test can be adjusted to account for dilution of the injection fluid in the aquifer fluid. The dilution-adjusted breakthrough curve can be analyzed to estimate the reaction rate of a solute. The conventional dilution-adjusted method assumes that the ratios of the concentrations of the nonreactive and reactive solutes in the injection fluid vs. the aquifer fluid are equal. If this assumption is invalid, the conventional method will generate inaccurate breakthrough curves and may lead to erroneous conclusions regarding the reactivity of a solute. In this study, a new method that generates a dilution-adjusted breakthrough curve was theoretically developed to account for any possible combination of nonreactive and reactive solute concentrations in the injection and aquifer fluids. The newly developed method was applied to a field-based data set and was shown to generate more accurate dilution-adjusted breakthrough curves. The improved dilution-adjusted method presented here is simple, makes no assumptions regarding the concentrations of the nonreactive and reactive solutes in the injection and aquifer fluids, and easily allows for estimating reaction rates during push-pull tests.

Dissolved plume attenuation with DNAPL source remediation, aqueous decay and volatilization--analytical solution, model calibration and prediction uncertainty.

A vertically-integrated analytical model for dissolved phase transport is described that considers a time-dependent DNAPL source based on the upscaled dissolution kinetics model of Parker and Park with extensions to consider time-dependent source zone biodecay, partial source mass reduction, and remediation-enhanced source dissolution kinetics. The model also considers spatial variability in aqueous plume decay, which is treated as the sum of aqueous biodecay and volatilization due to diffusive transport and barometric pumping through the unsaturated zone. The model is implemented in Excel/VBA coupled with (1) an inverse solution that utilizes prior information on model parameters and their uncertainty to condition the solution, and (2) an error analysis module that computes parameter covariances and total prediction uncertainty due to regression error and parameter uncertainty. A hypothetical case study is presented to evaluate the feasibility of calibrating the model from limited noisy field data. The results indicate that prediction uncertainty increases significantly over time following calibration, primarily due to propagation of parameter uncertainty. However, differences between the predicted performance of source zone partial mass reduction and the known true performance were reasonably small. Furthermore, a clear difference is observed between the predicted performance for the remedial action scenario versus that for a no-action scenario, which is consistent with the true system behavior. The results suggest that the model formulation can be effectively utilized to assess monitored natural attenuation and source remediation options if careful attention is given to model calibration and prediction uncertainty issues.

A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions.

This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing N(E) equilibrium reactions and a set of reactive transport equations of M-N(E) kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.

Estimating radiation dose from time to emesis and lymphocyte depletion.

Lymphocyte depletion kinetics and time to emesis have previously been shown to correlate with radiation dose. However, the method for estimating dose from lymphocyte counts was cumbersome, and a tabulation of estimated dose vs. time to emesis published by the International Atomic Energy Agency did not agree well with a regression of data from many (>100) radiation accident cases. The time-to-emesis data have been reanalyzed, and the new regression corroborates the previously published table. Also, dose estimation from post-exposure lymphocyte counts has been simplified and no longer requires serial calculations. Instead, dose can be estimated by a simple table lookup, given the ratio of two lymphocyte counts and the time between blood samples.

An upscaled approach for transport in media with extended tailing due to back-diffusion using analytical and numerical solutions of the advection dispersion equation.

The mono-continuum advection-dispersion equation (mADE) is commonly regarded as unsuitable for application to media that exhibit rapid breakthrough and extended tailing associated with diffusion between high and low permeability regions. This paper demonstrates that the mADE can be successfully used to model such conditions if certain issues are addressed. First, since hydrodynamic dispersion, unlike molecular diffusion, cannot occur upstream of the contaminant source, models must be formulated to prevent "back-dispersion." Second, large variations in aquifer permeability will result in differences between volume-weighted average concentration (resident concentration) and flow-weighted average concentration (flux concentration). Water samples taken from wells may be regarded as flux concentrations, while soil samples may be analyzed to determine resident concentrations. While the mADE is usually derived in terms of resident concentration, it is known that a mADE of the same mathematical form may be written in terms of flux concentration. However, when solving the latter, the mathematical transformation of a flux boundary condition applied to the resident mADE becomes a concentration type boundary condition for the flux mADE. Initial conditions must also be consistent with the form of the mADE that is to be solved. Thus, careful attention must be given to the type of concentration data that is available, whether resident or flux concentrations are to be simulated, and to boundary and initial conditions. We present 3-D analytical solutions for resident and flux concentrations, discuss methods of solving numerical models to obtain resident and flux concentrations, and compare results for hypothetical problems. We also present an upscaling method for computing "effective" dispersivities and other mADE model parameters in terms of physically meaningful parameters in a diffusion-limited mobile-immobile model. Application of the latter to previously published studies of systems that exhibit early breakthrough and extended tailing shows that the upscaled mADE model is able to describe the observed behavior with reasonable accuracy given only known physical parameters for the systems without any model calibration.

Fate and transport of uranium (VI) in weathered saprolite.

Batch and column experiments were conducted to investigate sorption and transport of uranium (U) in the presence of saprolite derived from interbedded shale, limestone, and sandstone sequences. Sorption kinetics were measured at two initial concentrations (C0; 1, 10 µM) and three soil:solution ratios (Rs/w; 0.005, 0.25, 2 kg/L) at pH 4.5 (pH of the saprolite). The rate of U loss from solution (µmole/L/h) increased with increasing Rs/w. Uranium sorption exhibited a fast phase with 80% sorption in the first eight hours for all C0 and Rs/w values and a slow phase during which the reaction slowly approached (pseudo)equilibrium over the next seven days. The pH-dependency of U sorption was apparent in pH sorption edges. U(VI) sorption increased over the pH range 4-6, then decreased sharply at pH > 7.5. U(VI) sorption edges were well described by a surface complexation model using calibrated parameters and the reaction network proposed by Waite et al. (1994). Sorption isotherms measured using the same Rs/w and pH values showed a solids concentration effect where U(VI) sorption capacity and affinity decreased with increasing solids concentration. This effect may have been due to either particle aggregation or competition between U(VI) and exchangeable cations for sorption sites. The surface complexation model with calibrated parameters was able to predict the general sorption behavior relatively well, but failed to reproduce solid concentration effects, implying the importance of appropriate design if batch experiments are to be utilized for dynamic systems. Transport of U(VI) through the packed column was significantly retarded. Transport simulations were conducted using the reactive transport model HydroGeoChem (HGC) v5.0 that incorporated the surface complexation reaction network used to model the batch data. Model parameters reported by Waite et al. (1994) provided a better prediction of U transport than optimized parameters derived from our sorption edges. The results presented in this study highlight the challenges in defining appropriate conditions for batch-type experiments used to extrapolate parameters for transport models, and also underline a gap in our ability to transfer batch results to transport simulations.

A spreadsheet program for two-well tracer test data analysis.

Two-well tracer tests are often conducted to investigate subsurface solute transport in the field. Analyzing breakthrough curves in extraction and monitoring wells using numerical methods is nontrivial due to highly nonuniform flow conditions. We extended approximate analytical solutions for the advection-dispersion equation for an injection-extraction well doublet in a homogeneous confined aquifer under steady-state flow conditions for equal injection and extraction rates with no transverse dispersion and negligible ambient flow, and implemented the solutions in Microsoft Excel using Visual Basic for Application (VBA). Functions were implemented to calculate concentrations in extraction and monitoring wells at any location due to a step or pulse injection. Type curves for a step injection were compared with those calculated by numerically integrating the solution for a pulse injection. The results from the two approaches are similar when the dispersivity is small. As the dispersivity increases, the latter was found to be more accurate but requires more computing time. The code was verified by comparing the results with published-type curves and applied to analyze data from the literature. The method can be used as a first approximation for two-well tracer test design and data analysis, and to check accuracy of numerical solutions. The code and example files are publicly available.

Modeling uranium transport in acidic contaminated groundwater with base addition.

This study investigates reactive transport modeling in a column of uranium(VI)-contaminated sediments with base additions in the circulating influent. The groundwater and sediment exhibit oxic conditions with low pH, high concentrations of NO(3)(-), SO(4)(2-), U and various metal cations. Preliminary batch experiments indicate that additions of strong base induce rapid immobilization of U for this material. In the column experiment that is the focus of the present study, effluent groundwater was titrated with NaOH solution in an inflow reservoir before reinjection to gradually increase the solution pH in the column. An equilibrium hydrolysis, precipitation and ion exchange reaction model developed through simulation of the preliminary batch titration experiments predicted faster reduction of aqueous Al than observed in the column experiment. The model was therefore modified to consider reaction kinetics for the precipitation and dissolution processes which are the major mechanism for Al immobilization. The combined kinetic and equilibrium reaction model adequately described variations in pH, aqueous concentrations of metal cations (Al, Ca, Mg, Sr, Mn, Ni, Co), sulfate and U(VI). The experimental and modeling results indicate that U(VI) can be effectively sequestered with controlled base addition due to sorption by slowly precipitated Al with pH-dependent surface charge. The model may prove useful to predict field-scale U(VI) sequestration and remediation effectiveness.

Prediction of uranium and technetium sorption during titration of contaminated acidic groundwater.

This study investigates uranium and technetium sorption onto aluminum and iron hydroxides during titration of acidic groundwater. The contaminated groundwater exhibits oxic conditions with high concentrations of NO(3)(-), SO(4)(2-), U, Tc, and various metal cations. More than 90% of U and Tc was removed from the aqueous phase as Al and Fe precipitated above pH 5.5, but was partially resolublized at higher pH values. An equilibrium hydrolysis and precipitation reaction model adequately described variations in aqueous concentrations of metal cations. An anion exchange reaction model was incorporated to simulate sulfate, U and Tc sorption onto variably charged (pH-dependent) Al and Fe hydroxides. Modeling results indicate that competitive sorption/desorption on mixed mineral phases needs to be considered to adequately predict U and Tc mobility. The model could be useful for future studies of the speciation of U, Tc and co-existing ions during pre- and post-groundwater treatment practices.

Kinetic analysis and modeling of oleate and ethanol stimulated uranium (VI) bio-reduction in contaminated sediments under sulfate reduction conditions.

Microcosm tests with uranium contaminated sediments were performed to explore the feasibility of using oleate as a slow-release electron donor for U(VI) reduction in comparison to ethanol. Oleate degradation proceeded more slowly than ethanol with acetate produced as an intermediate for both electron donors under a range of initial sulfate concentrations. A kinetic microbial reduction model was developed and implemented to describe and compare the reduction of sulfate and U(VI) with oleate or ethanol. The reaction path model considers detailed oleate/ethanol degradation and the production and consumption of intermediates, acetate and hydrogen. Although significant assumptions are made, the model tracked the major trend of sulfate and U(VI) reduction and describes the successive production and consumption of acetate, concurrent with microbial reduction of aqueous sulfate and U(VI) species. The model results imply that the overall rate of U(VI) bioreduction is influenced by both the degradation rate of organic substrates and consumption rate of intermediate products.

Geochemical modeling of reactions and partitioning of trace metals and radionuclides during titration of contaminated acidic sediments.

Many geochemical reactions that control aqueous metal concentrations are directly affected by solution pH. However, changes in solution pH are strongly buffered by various aqueous phase and solid phase precipitation/dissolution and adsorption/desorption reactions. The ability to predict acid-base behavior of the soil-solution system is thus critical to predict metal transport under variable pH conditions. This studywas undertaken to develop a practical generic geochemical modeling approach to predict aqueous and solid phase concentrations of metals and anions during conditions of acid or base additions. The method of Spalding and Spalding was utilized to model soil buffer capacity and pH-dependent cation exchange capacity by treating aquifer solids as a polyprotic acid. To simulate the dynamic and pH-dependent anion exchange capacity, the aquifer solids were simultaneously treated as a polyprotic base controlled by mineral precipitation/ dissolution reactions. An equilibrium reaction model that describes aqueous complexation, precipitation, sorption and soil buffering with pH-dependent ion exchange was developed using HydroGeoChem v5.0 (HGC5). Comparison of model results with experimental titration data of pH, Al, Ca, Mg, Sr, Mn, Ni, Co, and SO4(2-) for contaminated sediments indicated close agreement suggesting that the model could potentially be used to predictthe acid-base behavior of the sediment-solution system under variable pH conditions.