Chemical composition of arsenic-based acid mine drainage in the downstream of a gold mine: Fuzzy regression and clustering analysis.

This study employs fuzzy regression and fuzzy multivariate clustering techniques to analyze arsenic-polluted water samples originating from acid rock drainage in waste rock dumps. The research focuses on understanding the complex relationships between variables associated with arsenic contamination, such as water arsenic concentration, pH levels, and soil characteristics. To this end, fuzzy regression models were developed to estimate the relationships between water arsenic concentration and independent variables, thus, incorporating the inherent uncertainties into the analysis. Furthermore, multivariate fuzzy k-means clustering analysis facilitated the identification of fuzzy-based clusters within the dataset, providing insights into spatial patterns and potential sources of arsenic pollution. The pairwise comparisons indicated the strongest correlation of 0.62 between soil total arsenic and pH, while the weakest correlation of 0.13 was observed between soil-soluble arsenic and soil iron, providing valuable insights into their relationships and impact on water arsenic levels. The associated uncertainties in the relationships among the variables were determined based on the degree of belongingness of each data point to various fuzzy sets. Three distinct clusters emerged from the analysis: Cluster 1 comprised Points 5, 6, and 7; Cluster 2 included Points 1, 2, 3, 4, 8, and 9; and Cluster 3 consisted of Points 10, 11, 12, and 13. The findings enhance our understanding of the factors influencing arsenic contamination to provide an effective mitigation strategy in acid rock drainage scenarios. This research also demonstrates the applicability and effectiveness of fuzzy regression and fuzzy multivariate clustering in the analysis of arsenic-polluted water samples.

Interpreting water demands of forests and grasslands within a new Budyko formulation of evapotranspiration using percolation theory.

The relationship between carbon cycle and water demand is key to understanding global climate change, vegetation productivity, and predicting the future of water resources. The water balance, which enumerates the relative fractions of precipitation P that run off, Q, or are returned to the atmosphere through evapotranspiration, ET, links drawdown of atmospheric carbon with the water cycle through plant transpiration. Our theoretical description based on percolation theory proposes that dominant ecosystems tend to maximize drawdown of atmospheric carbon in the process of growth and reproduction, thus providing a link between carbon and water cycles. In this framework, the only parameter is the fractal dimensionality df of the root system. Values of df appear to relate to the relative roles of nutrient and water accessibility. Larger values of df lead to higher ET values. Known ranges of grassland root fractal dimensions predict reasonably the range of ET(P) in such ecosystems as a function of aridity index. Forests with shallower root systems, should be characterized by a smaller df and, therefore, ET that is a smaller fraction of P. The prediction of ET/P using the 3D percolation value of df matches rather closely results deemed typical for forests based on a phenomenology already in common use. We test predictions of Q with P against data and data summaries for sclerophyll forests in southeastern Australia and the southeastern USA. Applying PET data from a nearby site constrains the data from the USA to lie between our ET predictions for 2D and 3D root systems. For the Australian site, equating cited "losses" with PET underpredicts ET. This discrepancy is mostly removed by referring to mapped values of PET in that region. Missing in both cases is local PET variability, more important for reducing data scatter in southeastern Australia, due to the greater relief.

Mathematical analysis of laboratory microbial experiments demonstrating deterministic chaotic dynamics.

Presented is a system of four ordinary differential equations and a mathematical analysis of microbiological experiments in a four-component chemostat-nutrient n, rods r, cocci c, and predators p. The analysis is consistent with the conclusion that previous experiments produced features of deterministic chaotic and classical dynamics depending on dilution rate. The surrogate model incorporates as much experimental detail as possible, but necessarily contains unmeasured parameters. The objective is to understand better the differences between model simulations and experimental results in complex microbial populations. The key methodology for simulation of chaotic dynamics, consistent with the measured dilution rate and microbial volume averages, was to cause the preference of p for r vs. c to vary with the r and c concentrations, to make r more competitive for nutrient than c, and to recycle some dying p biomass, leading to a modified version of the Monod kinetics model. Our mathematical model demonstrated that the occurrence of chaotic dynamics requires a predator, p, preference for r versus c to increase significantly with increases in r and c populations. Also included is a discussion of several generalizations of the existing model and a possible involvement of the minimum energy dissipation principle. This principle appears fundamental to thermodynamic systems including living systems. Several new experiments are suggested.

Measurement of Volatile Compounds for Real-Time Analysis of Soil Microbial Metabolic Response to Simulated Snowmelt.

Snowmelt dynamics are a significant determinant of microbial metabolism in soil and regulate global biogeochemical cycles of carbon and nutrients by creating seasonal variations in soil redox and nutrient pools. With an increasing concern that climate change accelerates both snowmelt timing and rate, obtaining an accurate characterization of microbial response to snowmelt is important for understanding biogeochemical cycles intertwined with soil. However, observing microbial metabolism and its dynamics non-destructively remains a major challenge for systems such as soil. Microbial volatile compounds (mVCs) emitted from soil represent information-dense signatures and when assayed non-destructively using state-of-the-art instrumentation such as Proton Transfer Reaction-Time of Flight-Mass Spectrometry (PTR-TOF-MS) provide time resolved insights into the metabolism of active microbiomes. In this study, we used PTR-TOF-MS to investigate the metabolic trajectory of microbiomes from a subalpine forest soil, and their response to a simulated wet-up event akin to snowmelt. Using an information theory approach based on the partitioning of mutual information, we identified mVC metabolite pairs with robust interactions, including those that were non-linear and with time lags. The biological context for these mVC interactions was evaluated by projecting the connections onto the Kyoto Encyclopedia of Genes and Genomes (KEGG) network of known metabolic pathways. Simulated snowmelt resulted in a rapid increase in the production of trimethylamine (TMA) suggesting that anaerobic degradation of quaternary amine osmo/cryoprotectants, such as glycine betaine, may be important contributors to this resource pulse. Unique and synergistic connections between intermediates of methylotrophic pathways such as dimethylamine, formaldehyde and methanol were observed upon wet-up and indicate that the initial pulse of TMA was likely transformed into these intermediates by methylotrophs. Increases in ammonia oxidation signatures (transformation of hydroxylamine to nitrite) were observed in parallel, and while the relative role of nitrifiers or methylotrophs cannot be confirmed, the inferred connection to TMA oxidation suggests either a direct or indirect coupling between these processes. Overall, it appears that such mVC time-series from PTR-TOF-MS combined with causal inference represents an attractive approach to non-destructively observe soil microbial metabolism and its response to environmental perturbation.

Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest.

Deep-water access is arguably the most effective, but under-studied, mechanism that plants employ to survive during drought. Vulnerability to embolism and hydraulic safety margins can predict mortality risk at given levels of dehydration, but deep-water access may delay plant dehydration. Here, we tested the role of deep-water access in enabling survival within a diverse tropical forest community in Panama using a novel data-model approach. We inversely estimated the effective rooting depth (ERD, as the average depth of water extraction), for 29 canopy species by linking diameter growth dynamics (1990-2015) to vapor pressure deficit, water potentials in the whole-soil column, and leaf hydraulic vulnerability curves. We validated ERD estimates against existing isotopic data of potential water-access depths. Across species, deeper ERD was associated with higher maximum stem hydraulic conductivity, greater vulnerability to xylem embolism, narrower safety margins, and lower mortality rates during extreme droughts over 35 years (1981-2015) among evergreen species. Species exposure to water stress declined with deeper ERD indicating that trees compensate for water stress-related mortality risk through deep-water access. The role of deep-water access in mitigating mortality of hydraulically-vulnerable trees has important implications for our predictive understanding of forest dynamics under current and future climates.

Predicting Water Cycle Characteristics from Percolation Theory and Observational Data.

The fate of water and water-soluble toxic wastes in the subsurface is of high importance for many scientific and practical applications. Although solute transport is proportional to water flow rates, theoretical and experimental studies show that heavy-tailed (power-law) solute transport distribution can cause chemical transport retardation, prolonging clean-up time-scales greatly. However, no consensus exists as to the physical basis of such transport laws. In percolation theory, the scaling behavior of such transport rarely relates to specific medium characteristics, but strongly to the dimensionality of the connectivity of the flow paths (for example, two- or three-dimensional, as in fractured-porous media or heterogeneous sediments), as well as to the saturation characteristics (i.e., wetting, drying, and entrapped air). In accordance with the proposed relevance of percolation models of solute transport to environmental clean-up, these predictions also prove relevant to transport-limited chemical weathering and soil formation, where the heavy-tailed distributions slow chemical weathering over time. The predictions of percolation theory have been tested in laboratory and field experiments on reactive solute transport, chemical weathering, and soil formation and found accurate. Recently, this theoretical framework has also been applied to the water partitioning at the Earth's surface between evapotranspiration, ET, and run-off, Q, known as the water balance. A well-known phenomenological model by Budyko addressed the relationship between the ratio of the actual evapotranspiration (ET) and precipitation, ET/P, versus the aridity index, ET0/P, with P being the precipitation and ET0 being the potential evapotranspiration. Existing work was able to predict the global fractions of P represented by Q and ET through an optimization of plant productivity, in which downward water fluxes affect soil depth, and upward fluxes plant growth. In the present work, based likewise on the concepts of percolation theory, we extend Budyko's model, and address the partitioning of run-off Q into its surface and subsurface components, as well as the contribution of interception to ET. Using various published data sources on the magnitudes of interception and information regarding the partitioning of Q, we address the variability in ET resulting from these processes. The global success of this prediction demonstrated here provides additional support for the universal applicability of percolation theory for solute transport as well as guidance in predicting the component of subsurface run-off, important for predicting natural flow rates through contaminated aquifers.

Precipitation mediates sap flux sensitivity to evaporative demand in the neotropics.

Transpiration in humid tropical forests modulates the global water cycle and is a key driver of climate regulation. Yet, our understanding of how tropical trees regulate sap flux in response to climate variability remains elusive. With a progressively warming climate, atmospheric evaporative demand [i.e., vapor pressure deficit (VPD)] will be increasingly important for plant functioning, becoming the major control of plant water use in the twenty-first century. Using measurements in 34 tree species at seven sites across a precipitation gradient in the neotropics, we determined how the maximum sap flux velocity (vmax) and the VPD threshold at which vmax is reached (VPDmax) vary with precipitation regime [mean annual precipitation (MAP); seasonal drought intensity (PDRY)] and two functional traits related to foliar and wood economics spectra [leaf mass per area (LMA); wood specific gravity (WSG)]. We show that, even though vmax is highly variable within sites, it follows a negative trend in response to increasing MAP and PDRY across sites. LMA and WSG exerted little effect on vmax and VPDmax, suggesting that these widely used functional traits provide limited explanatory power of dynamic plant responses to environmental variation within hyper-diverse forests. This study demonstrates that long-term precipitation plays an important role in the sap flux response of humid tropical forests to VPD. Our findings suggest that under higher evaporative demand, trees growing in wetter environments in humid tropical regions may be subjected to reduced water exchange with the atmosphere relative to trees growing in drier climates.

Modeling the impact of liana infestation on the demography and carbon cycle of tropical forests.

There is mounting empirical evidence that lianas affect the carbon cycle of tropical forests. However, no single vegetation model takes into account this growth form, although such efforts could greatly improve the predictions of carbon dynamics in tropical forests. In this study, we incorporated a novel mechanistic representation of lianas in a dynamic global vegetation model (the Ecosystem Demography Model). We developed a liana-specific plant functional type and mechanisms representing liana-tree interactions (such as light competition, liana-specific allometries, and attachment to host trees) and parameterized them according to a comprehensive literature meta-analysis. We tested the model for an old-growth forest (Paracou, French Guiana) and a secondary forest (Gigante Peninsula, Panama). The resulting model simulations captured many features of the two forests characterized by different levels of liana infestation as revealed by a systematic comparison of the model outputs with empirical data, including local census data from forest inventories, eddy flux tower data, and terrestrial laser scanner-derived forest vertical structure. The inclusion of lianas in the simulations reduced the secondary forest net productivity by up to 0.46 tC  ha-1  year-1 , which corresponds to a limited relative reduction of 2.6% in comparison with a reference simulation without lianas. However, this resulted in significantly reduced accumulated above-ground biomass after 70 years of regrowth by up to 20 tC /ha (19% of the reference simulation). Ultimately, the simulated negative impact of lianas on the total biomass was almost completely cancelled out when the forest reached an old-growth successional stage. Our findings suggest that lianas negatively influence the forest potential carbon sink strength, especially for young, disturbed, liana-rich sites. In light of the critical role that lianas play in the profound changes currently experienced by tropical forests, this new model provides a robust numerical tool to forecast the impact of lianas on tropical forest carbon sinks.

Use of carbon stable isotopes to monitor biostimulation and electron donor fate in chromium-contaminated groundwater.

Hexavalent chromium Cr(VI) is a common inorganic contaminant in industrial areas and represents a serious threat to human health due its toxicity. Here we report experimental results from a field-scale investigation of Cr(VI) bio-immobilization at Hanford 100H reservation, a U.S Department of Energy facility (Washington State, USA). Microbial Cr(VI) reduction was stimulated via injection of a13C-labeled sodium lactate solution into the high-permeability aquifer consisting of gravel and coarse sand sediments. Concentrations and carbon isotope ratios of metabolites, including dissolved inorganic carbon and total organic carbon, and compound-specific analysis of acetate and propionate, together with phospholipid fatty acids (biomass) have been analyzed to help provide an understanding of the predominant redox processes accompanying Cr(VI) reduction. Results of our study indicate that the injection of an electron donor caused a sharp decrease of Cr(VI) concentration from ∼32 to ∼10 nM. Cr(VI) reduction was associated with a decrease in the concentration of carboxylic acids, such as lactate (∼6 mM to undetectable), propionate (∼9 mM to undetectable), and acetate (∼6 mM to undetectable), as well as dissolved inorganic carbon (30-10 mM C). Carbon isotope data indicate carbon transfers from the original substrate to organic byproducts and mineralized carbon. Concentrations of metabolites and stable isotope data as well as carbon isotope mass balance calculations were used to monitor biologically mediated reduction of Cr(VI).

Climate change impact on residual contaminants under sustainable remediation.

This study investigates the potential impact of climate change on residual contaminants in vadose zones and groundwater. We assume that the effect of climate changes can be represented by perturbations in the natural recharge through the aquifer system. We perform numerical modeling of unsaturated/saturated flow and transport and consider different performance metrics: contaminant concentrations at observation wells and contaminant export at the site's boundary. We evaluate the effect of increasing and decreasing recharge as well as the impact of potential failure of surface capping structures employed to immobilize vadose zone contaminants. Our approach is demonstrated in a real case study by simulating transport of non-reactive radioactive tritium at the U.S. Department of Energy's Savannah River Site. Results show that recharge changes significantly affect well concentrations: after an initial slight dilution we identify a significant concentration increase at different observation wells some years after the recharge increase and/or the cap failure, as a consequence of contaminants' mobilization. This effect is generally emphasized and occurs earlier as the recharge increases. Under decreased aquifers' recharge the concentration could slightly increase for some years, due to a decrease of dilution, depending on the magnitude of the negative recharge shift. We identify trigger levels of recharge above which the concentration/export breakthrough curves and the time of exceedance of the Maximum Contaminant Level for tritium are remarkably affected. Moreover, we observe that the contaminant export at the control plane, identified as the risk pathway to the downgradient population, may only be minimally affected by shifts in the natural recharge regime, except for some extreme cases. We conclude that more frequent sampling and in-situ monitoring near the source zone should be adopted to better explain concentrations' anomalies under changing climatic conditions. Moreover, the maintenance of the cap is critical not only to sequester residual contaminants in the vadose zone, but also to reduce the uncertainty associated with future precipitation changes. Finally, realistic flow and transport simulations achieved through proper calibration processes, rather than conservative modeling, should be adopted to identify non-trivial trade-offs which enable better allocation of resources towards reducing uncertainty in decision making.

In Situ Monitoring of Groundwater Contamination Using the Kalman Filter.

This study presents a Kalman filter-based framework to establish a real-time in situ monitoring system for groundwater contamination based on in situ measurable water quality variables, such as specific conductance (SC) and pH. First, this framework uses principal component analysis (PCA) to identify correlations between the contaminant concentrations of interest and in situ measurable variables. It then applies the Kalman filter to estimate contaminant concentrations continuously and in real-time by coupling data-driven concentration-decay models with the previously identified data correlations. We demonstrate our approach with historical groundwater data from the Savannah River Site F-Area: We use SC and pH data to estimate tritium and uranium concentrations over time. Results show that the developed method can estimate these contaminant concentrations based on in situ measurable variables. The estimates remain reliable with less frequent or no direct measurements of the contaminant concentrations, while capturing the dynamics of short- and long-term contaminant concentration changes. In addition, we show that data mining, such as PCA, is useful to understand correlations in groundwater data and to design long-term monitoring systems. The developed in situ monitoring methodology is expected to improve long-term groundwater monitoring by continuously confirming the contaminant plume's stability and by providing an early warning system for unexpected changes in the plume's migration.

Variation in hydroclimate sustains tropical forest biomass and promotes functional diversity.

The fate of tropical forests under climate change is unclear as a result, in part, of the uncertainty in projected changes in precipitation and in the ability of vegetation models to capture the effects of drought-induced mortality on aboveground biomass (AGB). We evaluated the ability of a terrestrial biosphere model with demography and hydrodynamics (Ecosystem Demography, ED2-hydro) to simulate AGB and mortality of four tropical tree plant functional types (PFTs) that operate along light- and water-use axes. Model predictions were compared with observations of canopy trees at Barro Colorado Island (BCI), Panama. We then assessed the implications of eight hypothetical precipitation scenarios, including increased annual precipitation, reduced inter-annual variation, El Niño-related droughts and drier wet or dry seasons, on AGB and functional diversity of the model forest. When forced with observed meteorology, ED2-hydro predictions capture multiple BCI benchmarks. ED2-hydro predicts that AGB will be sustained under lower rainfall via shifts in the functional composition of the forest, except under the drier dry-season scenario. These results support the hypothesis that inter-annual variation in mean and seasonal precipitation promotes the coexistence of functionally diverse PFTs because of the relative differences in mortality rates. If the hydroclimate becomes chronically drier or wetter, functional evenness related to drought tolerance may decline.

Using strontium isotopes to evaluate the spatial variation of groundwater recharge.

Recharge of alluvial aquifers is a key component in understanding the interaction between floodplain vadose zone biogeochemistry and groundwater quality. The Rifle Site (a former U-mill tailings site) adjacent to the Colorado River is a well-established field laboratory that has been used for over a decade for the study of biogeochemical processes in the vadose zone and aquifer. This site is considered an exemplar of both a riparian floodplain in a semiarid region and a post-remediation U-tailings site. In this paper we present Sr isotopic data for groundwater and vadose zone porewater samples collected in May and July 2013 to build a mixing model for the fractional contribution of vadose zone porewater (i.e. recharge) to the aquifer and its variation across the site. The vadose zone porewater contribution to the aquifer ranged systematically from 0% to 38% and appears to be controlled largely by the microtopography of the site. The area-weighted average contribution across the site was 8% corresponding to a net recharge of 7.5 cm. Given a groundwater transport time across the site of ~1.5 to 3 years, this translates to a recharge rate between 5 and 2.5 cm/yr, and with the average precipitation to the site implies a loss from the vadose zone due to evapotranspiration of 83% to 92%, both ranges are in good agreement with previously published results by independent methods. A uranium isotopic (234U/238U activity ratios) mixing model for groundwater and surface water samples indicates that a ditch across the site is hydraulically connected to the aquifer and locally significantly affects groundwater. Groundwater samples with high U concentrations attributed to natural bio-reduced zones have 234U/238U activity ratios near 1, suggesting that the U currently being released to the aquifer originated from the former U-mill tailings.

Identifying key controls on the behavior of an acidic-U(VI) plume in the Savannah River Site using reactive transport modeling.

Acidic low-level waste radioactive waste solutions were discharged to three unlined seepage basins at the F-Area of the Department of Energy (DOE) Savannah River Site (SRS), South Carolina, USA, from 1955 through 1989. Despite many years of active remediation, the groundwater remains acidic and contaminated with significant levels of U(VI) and other radionuclides. Monitored Natural Attenuation (MNA) is a desired closure strategy for the site, based on the premise that regional flow of clean background groundwater will eventually neutralize the groundwater acidity, immobilizing U(VI) through adsorption. An in situ treatment system is currently in place to accelerate this in the downgradient portion of the plume and similar measures could be taken upgradient if necessary. Understanding the long-term pH and U(VI) adsorption behavior at the site is critical to assess feasibility of MNA along with the in-situ remediation treatments. This paper presents a reactive transport (RT) model and uncertainty quantification (UQ) analyses to explore key controls on the U(VI)-plume evolution and long-term mobility at this site. Two-dimensional numerical RT simulations are run including the saturated and unsaturated (vadose) zones, U(VI) and H(+) adsorption (surface complexation) onto sediments, dissolution and precipitation of Al and Fe minerals, and key hydrodynamic processes are considered. UQ techniques are applied using a new open-source tool that is part of the developing ASCEM reactive transport modeling and analysis framework to: (1) identify the complex physical and geochemical processes that control the U(VI) plume migration in the pH range where the plume is highly mobile, (2) evaluate those physical and geochemical parameters that are most controlling, and (3) predict the future plume evolution constrained by historical, chemical and hydrological data. The RT simulation results show a good agreement with the observed historical pH and concentrations of U(VI), nitrates and Al concentrations at multiple locations. Mineral dissolution and precipitation combined with adsorption reactions on goethite and kaolinite (the main minerals present with quartz) could buffer pH at the site for long periods of time. UQ analysis using the Morris one-at-a-time (OAT) method indicates that the model/parameter is most sensitive to the pH of the waste solution, discharge rates, and the reactive surface area available for adsorption. However, as a key finding, UQ analysis also indicates that this model (and parameters) sensitivity evolves in space and time, and its understanding could be crucial to assess the temporal efficiency of a remediation strategy in contaminated sites. Results also indicate that residual U(VI) and H(+) adsorbed in the vadose zone, as well as aquifer permeability, could have a significant impact on the acidic plume long-term mobility.

Microbial community succession during lactate amendment and electron acceptor limitation reveals a predominance of metal-reducing Pelosinus spp.

The determination of the success of in situ bioremediation strategies is complex. By using controlled laboratory conditions, the influence of individual variables, such as U(VI), Cr(VI), and electron donors and acceptors on community structure, dynamics, and the metal-reducing potential can be studied. Triplicate anaerobic, continuous-flow reactors were inoculated with Cr(VI)-contaminated groundwater from the Hanford, WA, 100-H area, amended with lactate, and incubated for 95 days to obtain stable, enriched communities. The reactors were kept anaerobic with N(2) gas (9 ml/min) flushing the headspace and were fed a defined medium amended with 30 mM lactate and 0.05 mM sulfate with a 48-h generation time. The resultant diversity decreased from 63 genera within 12 phyla to 11 bacterial genera (from 3 phyla) and 2 archaeal genera (from 1 phylum). Final communities were dominated by Pelosinus spp. and to a lesser degree, Acetobacterium spp., with low levels of other organisms, including methanogens. Four new strains of Pelosinus were isolated, with 3 strains being capable of Cr(VI) reduction while one also reduced U(VI). Under limited sulfate, it appeared that the sulfate reducers, including Desulfovibrio spp., were outcompeted. These results suggest that during times of electron acceptor limitation in situ, organisms such as Pelosinus spp. may outcompete the more-well-studied organisms while maintaining overall metal reduction rates and extents. Finally, lab-scale simulations can test new strategies on a smaller scale while facilitating community member isolation, so that a deeper understanding of community metabolism can be revealed.

Microbial community response to addition of polylactate compounds to stimulate hexavalent chromium reduction in groundwater.

To evaluate the efficacy of bioimmobilization of Cr(VI) in groundwater at the Department of Energy Hanford site, we conducted a series of microcosm experiments using a range of commercial electron donors with varying degrees of lactate polymerization (polylactate). These experiments were conducted using Hanford Formation sediments (coarse sand and gravel) immersed in Hanford groundwater, which were amended with Cr(VI) and several types of lactate-based electron donors (Hydrogen Release Compound, HRC; primer-HRC, pHRC; extended release HRC) and the polylactate-cysteine form (Metal Remediation Compound, MRC). The results showed that polylactate compounds stimulated an increase in bacterial biomass and activity to a greater extent than sodium lactate when applied at equivalent carbon concentrations. At the same time, concentrations of headspace hydrogen and methane increased and correlated with changes in the microbial community structure. Enrichment of Pseudomonas spp. occurred with all lactate additions, and enrichment of sulfate-reducing Desulfosporosinus spp. occurred with almost complete sulfate reduction. The results of these experiments demonstrate that amendment with the pHRC and MRC forms result in effective removal of Cr(VI) from solution most likely by both direct (enzymatic) and indirect (microbially generated reductant) mechanisms.

In situ long-term reductive bioimmobilization of Cr(VI) in groundwater using hydrogen release compound.

The results of a field experiment designed to test the effectiveness of a novel approach for long-term, in situ bioimmobilization of toxic and soluble Cr(VI) in groundwater using a hydrogen release compound (HRC)--a slow release glycerol polylactate--are described. The field experiment was conducted at the Hanford Site (Washington), a U.S. Department of Energy nuclear production facility, using a combination of hydrogeological, geophysical, geochemical, and microbiological measurements and analyses of water samples and sediments. The results of this experiment show that a single HRC injection into groundwater stimulates an increase in biomass, a depletion of terminal electron acceptors O2, NO3-, and SO4(2-), and an increase in Fe2+, resulting in a significant decrease in soluble Cr(VI). The Cr(VI) concentration has remained below the background concentration in the downgradient pumping/ monitoring well, and below the detection limit in the injection well for more than 3 years after the HRC injection. The degree of sustainability of Cr(VI) reductive bioimmobilization under different redox conditions at this and other contaminated sites is currently under study.

Geophysical monitoring of hydrological and biogeochemical transformations associated with Cr(VI) bioremediation.

Understanding how hydrological and biogeochemical properties change over space and time in response to remedial treatments is hindered by our ability to monitor these processes with sufficient resolution and over field relevant scales. Here, we explored the use of geophysical approaches for monitoring the spatiotemporal distribution of hydrological and biogeochemical transformations associated with a Cr(VI) bioremediation experiment performed at Hanford, WA. We first integrated hydrological wellbore and geophysical tomographic data sets to estimate hydrological zonation at the study site. Using results from laboratory biogeophysical experiments and constraints provided by field geochemical data sets, we then interpreted time-lapse seismic and radar tomographic data sets, collected during thirteen acquisition campaigns over a three year experimental period, in terms of hydrological and biogeochemical transformations. The geophysical monitoring data sets were used to infer: the spatial distribution of injected electron donor; the evolution of gas bubbles; variations in total dissolved solids (nitrate and sulfate) as a function of pumping activity; the formation of precipitates and dissolution of calcites; and concomitant changes in porosity. Although qualitative in nature, the integrated interpretation illustrates how geophysical techniques have the potential to provide a wealth of information about coupled hydrobiogeochemical responses to remedial treatments in high spatial resolution and in a minimally invasive manner. Particularly novel aspects of our study include the use of multiple lines of evidence to constrain the interpretation of a long-term, field-scale geophysical monitoring data set and the interpretation of the transformations as a function of hydrological heterogeneity and pumping activity.