Philip M. Jardine

Significant Association between Sulfate-Reducing Bacteria and Uranium-Reducing Microbial Communities as Revealed by a Combined Massively Parallel Sequencing-Indicator Species Approach

ABSTRACT Massively parallel sequencing has provided a more affordable and high-throughput method to study microbial communities, although it has mostly been used in an exploratory fashion. We combined pyrosequencing with a strict indicator species statistical analysis to test if bacteria specifically responded to ethanol injection that successfully promoted dissimilatory uranium(VI) reduction in the subsurface of a uranium contamination plume at the Oak Ridge Field Research Center in Tennessee. Remediation was achieved with a hydraulic flow control consisting of an inner loop, where ethanol was injected, and an outer loop for flow-field protection. This strategy reduced uranium concentrations in groundwater to levels below 0.126 μM and created geochemical gradients in electron donors from the inner-loop injection well toward the outer loop and downgradient flow path. Our analysis with 15 sediment samples from the entire test area found significant indicator species that showed a high degree of adaptation to the three different hydrochemical-created conditions. Castellaniella and Rhodanobacter characterized areas with low pH, heavy metals, and low bioactivity, while sulfate-, Fe(III)-, and U(VI)-reducing bacteria (Desulfovibrio, Anaeromyxobacter, and Desulfosporosinus) were indicators of areas where U(VI) reduction occurred. The abundance of these bacteria, as well as the Fe(III) and U(VI) reducer Geobacter, correlated with the hydraulic connectivity to the substrate injection site, suggesting that the selected populations were a direct response to electron donor addition by the groundwater flow path. A false-discovery-rate approach was implemented to discard false-positive results by chance, given the large amount of data compared.

In Situ Bioremediation of Uranium with Emulsified Vegetable Oil as the Electron Donor

A field test with a one-time emulsified vegetable oil (EVO) injection was conducted to assess the capacity of EVO to sustain uranium bioreduction in a high-permeability gravel layer with groundwater concentrations of (mM) U, 0.0055; Ca, 2.98; NO3(-), 0.11; HCO3(-), 5.07; and SO4(2-), 1.23. Comparison of bromide and EVO migration and distribution indicated that a majority of the injected EVO was retained in the subsurface from the injection wells to 50 m downgradient. Nitrate, uranium, and sulfate were sequentially removed from the groundwater within 1-2 weeks, accompanied by an increase in acetate, Mn, Fe, and methane concentrations. Due to the slow release and degradation of EVO with time, reducing conditions were sustained for approximately one year, and daily U discharge to a creek, located approximately 50 m from the injection wells, decreased by 80% within 100 days. Total U discharge was reduced by 50% over the one-year period. Reduction of U(VI) to U(IV) was confirmed by synchrotron analysis of recovered aquifer solids. Oxidants (e.g., dissolved oxygen, nitrate) flowing in from upgradient appeared to reoxidize and remobilize uranium after the EVO was exhausted as evidenced by a transient increase of U concentration above ambient values. Occasional (e.g., annual) EVO injection into a permeable Ca and bicarbonate-containing aquifer can sustain uranium bioreduction/immobilization and decrease U migration/discharge.

Dynamics of microbial community composition and function during in situ bioremediation of a uranium-contaminated aquifer.

ABSTRACT A pilot-scale system was established to examine the feasibility of in situ U(VI) immobilization at a highly contaminated aquifer (U.S. DOE Integrated Field Research Challenge site, Oak Ridge, TN). Ethanol was injected intermittently as an electron donor to stimulate microbial U(VI) reduction, and U(VI) concentrations fell to below the Environmental Protection Agency drinking water standard (0.03 mg liter−1). Microbial communities from three monitoring wells were examined during active U(VI) reduction and maintenance phases with GeoChip, a high-density, comprehensive functional gene array. The overall microbial community structure exhibited a considerable shift over the remediation phases examined. GeoChip-based analysis revealed that Fe(III)-reducing bacterial (FeRB), nitrate-reducing bacterial (NRB), and sulfate-reducing bacterial (SRB) functional populations reached their highest levels during the active U(VI) reduction phase (days 137 to 370), in which denitrification and Fe(III) and sulfate reduction occurred sequentially. A gradual decrease in these functional populations occurred when reduction reactions stabilized, suggesting that these functional populations could play an important role in both active U(VI) reduction and maintenance of the stability of reduced U(IV). These results suggest that addition of electron donors stimulated the microbial community to create biogeochemical conditions favorable to U(VI) reduction and prevent the reduced U(IV) from reoxidation and that functional FeRB, SRB, and NRB populations within this system played key roles in this process.

CXTFIT/Excel-A modular adaptable code for parameter estimation, sensitivity analysis and uncertainty analysis for laboratory or field tracer experiments

We implemented the widely used CXTFIT code in Excel to provide flexibility and added sensitivity and uncertainty analysis functions to improve transport parameter estimation and to facilitate model discrimination for multi-tracer experiments on structured soils. Analytical solutions for one-dimensional equilibrium and nonequilibrium convection dispersion equations were coded as VBA functions so that they could be used as ordinary math functions in Excel for forward predictions. Macros with user-friendly interfaces were developed for optimization, sensitivity analysis, uncertainty analysis, error propagation, response surface calculation, and Monte Carlo analysis. As a result, any parameter with transformations (e.g., dimensionless, log-transformed, species-dependent reactions, etc.) could be estimated with uncertainty and sensitivity quantification for multiple tracer data at multiple locations and times. Prior information and observation errors could be incorporated into the weighted nonlinear least squares method with a penalty function. Users are able to change selected parameter values and view the results via embedded graphics, resulting in a flexible tool applicable to modeling transport processes and to teaching students about parameter estimation. The code was verified by comparing to a number of benchmarks with CXTFIT 2.0. It was applied to improve parameter estimation for four typical tracer experiment data sets in the literature using multi-model evaluation and comparison. Additional examples were included to illustrate the flexibilities and advantages of CXTFIT/Excel. The VBA macros were designed for general purpose and could be used for any parameter estimation/model calibration when the forward solution is implemented in Excel. A step-by-step tutorial, example Excel files and the code are provided as supplemental material.

Application of a hybrid MPI/OpenMP approach for parallel groundwater model calibration using multi-core computers

Calibration of groundwater models involves hundreds to thousands of forward solutions, each of which may solve many transient coupled nonlinear partial differential equations, resulting in a computationally intensive problem. We describe a hybrid MPI/OpenMP approach to exploit two levels of parallelisms in software and hardware to reduce calibration time on multi-core computers. HydroGeoChem 5.0 (HGC5) is parallelized using OpenMP for direct solutions for a reactive transport model application, and a field-scale coupled flow and transport model application. In the reactive transport model, a single parallelizable loop is identified to account for over 97% of the total computational time using GPROF. Addition of a few lines of OpenMP compiler directives to the loop yields a speedup of about 10 on a 16-core compute node. For the field-scale model, parallelizable loops in 14 of 174 HGC5 subroutines that require 99% of the execution time are identified. As these loops are parallelized incrementally, the scalability is found to be limited by a loop where Cray PAT detects over 90% cache missing rates. With this loop rewritten, similar speedup as the first application is achieved. The OpenMP-parallelized code can be run efficiently on multiple workstations in a network or multiple compute nodes on a cluster as slaves using parallel PEST to speedup model calibration. To run calibration on clusters as a single task, the Levenberg-Marquardt algorithm is added to HGC5 with the Jacobian calculation and lambda search parallelized using MPI. With this hybrid approach, 100-200 compute cores are used to reduce the calibration time from weeks to a few hours for these two applications. This approach is applicable to most of the existing groundwater model codes for many applications.

Simulation of carbon cycling, including dissolved organic carbon transport, in forest soil locally enriched with 14C

The DyDOC model was used to simulate the soil carbon cycle of a deciduous forest at the Oak Ridge Reservation (Tennessee, USA). The model application relied on extensive data from the Enriched Background Isotope Study (EBIS), which exploited a short-term local atmospheric enrichment of radiocarbon to establish a large-scale manipulation experiment with different inputs of 14C from both above-ground and below-ground litter. The model was first fitted to hydrological data, then observed pools and fluxes of carbon and 14C data were used to fit parameters describing metabolic transformations of soil organic matter (SOM) components and the transport and sorption of dissolved organic matter (DOM). This produced a detailed quantitative description of soil C cycling in the three horizons (O, A, B) of the soil profile. According to the parameterised model, SOM turnover within the thin O-horizon rapidly produces DOM (46xa0gCxa0m−2xa0a−1), which is predominantly hydrophobic. This DOM is nearly all adsorbed in the A- and B-horizons, and while most is mineralised relatively quickly, 11xa0gCxa0m−2xa0a−1 undergoes a “maturing” reaction, producing mineral-associated stable SOM pools with mean residence times of 100–200xa0years. Only a small flux (~1xa0gCxa0m−2xa0a−1) of hydrophilic DOM leaves the B-horizon. The SOM not associated with mineral matter is assumed to be derived from root litter, and turns over quite quickly (mean residence time 20–30xa0years). Although DyDOC was successfully fitted to C pools, annual fluxes and 14C data, it accounted less well for short-term variations in DOC concentrations.

Influence of soil geochemical and physical properties on chromium(VI) sorption and bioaccessibility.

The Department of Defense (DoD) is faced with the daunting task of possible remediation of numerous soil-Cr(VI) contaminated sites throughout the continental U.S. The primary risk driver at these sites is hand-to-mouth ingestion of contaminated soil by children. In the following study we investigate the impact of soil geochemical and physical properties on the sorption and bioaccessibility of Cr(VI) in a vast array of soils relevant to neighboring DoD sites. For the 35 soils used in this study, A-horizon soils typically sorbed significantly more Cr(VI) relative to B-horizon soils. Multiple linear regression analysis suggested that Cr(VI) sorption increased with increasing soil total organic C (TOC) and decreasing soil pH. The bioaccessibility of total Cr (CrT) and Cr(VI) on the soils decreased with increasing soil TOC content. As the soil TOC content approached 0.4%, the bioaccessibility of soil bound Cr systematically decreased from approximately 65 to 10%. As the soil TOC content increased from 0.4 to 4%, the bioaccessibility of Cr(VI) and CrT remained relatively constant at approximately 4% and 10%, respectively. X-ray absorption near edge structure (XANES) spectroscopy suggested that Cr(VI) reduction to Cr(III) was prevalent and that the redox transformation of Cr(VI) increased with increasing soil TOC. XANES confirmed that nearly all bioaccessible soil Cr was the Cr(VI) moiety. Multiple linear regression analysis suggested that the bioaccessibility of Cr(VI) and its reduced counterpart Cr(III), decreased with increasing soil TOC and increasing soil pH. This is consistent with the observation that the reduction reaction and formation of Cr(III) increased with increasing soil TOC and that Cr(III) was significantly less bioaccessible relative to Cr(VI). The model was found to adequately describe CrT bioaccessibility in soils from DoD facilities where Cr(VI) contaminated sites were present. The results of this study illustrate the importance of soil properties on Cr(VI) sorption and bioassessability and help define what soil types have the greatest risk associated with Cr(VI) exposure.

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.

Firing range soils yield a diverse array of fungal isolates capable of organic acid production and Pb-mineral solubilization

ABSTRACT Anthropogenic sources of lead contamination in soils include mining and smelting activities, effluents and wastes, agricultural pesticides, domestic garbage dumps, and shooting ranges. While Pb is typically considered relatively insoluble in the soil environment, some fungi may potentially contribute to mobilization of heavy metal cations by means of secretion of low-molecular-weight organic acids (LMWOAs). We sought to better understand the potential for metal mobilization within an indigenous fungal community at an abandoned shooting range in Oak Ridge, TN, where soil Pb contamination levels ranged from 24 to >2,700 mg Pb kg dry soil−1. We utilized culture-based assays to determine organic acid secretion and Pb-carbonate dissolution of a diverse collection of soil fungal isolates derived from the site and verified isolate distribution patterns within the community by 28S rRNA gene analysis of whole soils. The fungal isolates examined included both ascomycetes and basidiomycetes that excreted high levels (up to 27 mM) of a mixture of LMWOAs, including oxalic and citric acids, and several isolates demonstrated a marked ability to dissolve Pb-carbonate at high concentrations up to 10.5 g liter−1 (18.5 mM) in laboratory assays. Fungi within the indigenous community of these highly Pb-contaminated soils are capable of LMWOA secretion at levels greater than those of well-studied model organisms, such as Aspergillus niger. Additionally, these organisms were found in high relative abundance (>1%) in some of the most heavily contaminated soils. Our data highlight the need to understand more about autochthonous fungal communities at Pb-contaminated sites and how they may impact Pb biogeochemistry, solubility, and bioavailability, thus consequently potentially impacting human and ecosystem health.

Microbial communities biostimulated by ethanol during uranium (VI) bioremediation in contaminated sediment as shown by stable isotope probing

Stable isotope probing (SIP) was used to identify microbes stimulated by ethanol addition in microcosms containing two sediments collected from the bioremediation test zone at the US Department of Energy Oak Ridge site, TN, USA. One sample was highly bioreduced with ethanol while another was less reduced. Microcosms with the respective sediments were amended with 13C labeled ethanol and incubated for 7 days for SIP. Ethanol was rapidly converted to acetate within 24 h accompanied with the reduction of nitrate and sulfate. The accumulation of acetate persisted beyond the 7 d period. Aqueous U did not decline in the microcosm with the reduced sediment due to desorption of U but continuously declined in the less reduced sample. Microbial growth and concomitant 13C-DNA production was detected when ethanol was exhausted and abundant acetate had accumulated in both microcosms. This coincided with U(VI) reduction in the less reduced sample. 13C originating from ethanol was ultimately utilized for growth, either directly or indirectly, by the dominant microbial community members within 7 days of incubation. The microbial community was comprised predominantly of known denitrifiers, sulfate-reducing bacteria and iron (III) reducing bacteria including Desulfovibrio, Sphingomonas, Ferribacterium, Rhodanobacter, Geothrix, Thiobacillus and others, including the known U(VI)-reducing bacteria Acidovorax, Anaeromyxobacter, Desulfovibrio, Geobacter and Desulfosporosinus. The findings suggest that ethanol biostimulates the U(VI)-reducing microbial community by first serving as an electron donor for nitrate, sulfate, iron (III) and U(VI) reduction, and acetate which then functions as electron donor for U(VI) reduction and carbon source for microbial growth.