Toxicity Assessment of Groundwater Contaminated by Petroleum Hydrocarbons at a Well-Characterized, Aged, Crude Oil Release Site.

Management of petroleum-impacted waters by monitored natural attenuation requires an understanding of the toxicology of both the original compounds released and the transformation products formed during natural breakdown. Here, we report data from a groundwater plume consisting of a mixture of crude oil compounds and transformation products in an effort to bridge the gap between groundwater quality information and potential biological effects of human exposures. Groundwater samples were characterized for redox processes, concentrations of nonvolatile dissolved organic carbon (NVDOC) and total petroleum hydrocarbons in the diesel range, as well as for activation of human nuclear receptors (hNR) and toxicologically relevant transcriptional pathways. Results show upregulation of several biological pathways, including peroxisome proliferator-activated receptor gamma and alpha, estrogen receptor alpha, and pregnane X receptor (PXR) with higher levels of hNR activity observed in more contaminated samples. Our study of affected groundwater contaminated by a crude-oil release 39 years ago shows these types of waters may have the potential to cause adverse impacts on development, endocrine, and liver functioning in exposed populations. Additionally, positive trends in activation of some of the molecular targets (e.g., PXR) with increasing NVDOC concentrations (including polar transformation products) demonstrate the importance of improving our understanding of the toxicity associated with the unknown transformation products present in hydrocarbon-impacted waters. Our results begin to provide insight into the potential toxicity of petroleum-impacted waters, which is particularly timely given the ubiquitous nature of waters impacted by petroleum contamination not only recently but also in the past and the need to protect drinking-water quality.

Examining Natural Attenuation and Acute Toxicity of Petroleum-Derived Dissolved Organic Matter with Optical Spectroscopy.

Groundwater samples containing petroleum-derived dissolved organic matter (DOMHC) originating from the north oil body within the National Crude Oil Spill Fate and Natural Attenuation Research Site near Bemidji, MN, USA were analyzed by optical spectroscopic techniques (i.e., absorbance and fluorescence) to assess relationships that can be used to examine natural attenuation and toxicity of DOMHC in contaminated groundwater. A strong correlation between the concentration of dissolved organic carbon (DOC) and absorbance at 254 nm ( a254) along a transect of the DOMHC plume indicates that a254 can be used to quantitatively assess natural attenuation of DOMHC. Fluorescence components, identified by parallel factor (PARAFAC) analysis, show that the composition of the DOMHC beneath and adjacent to the oil body is dominated by aliphatic, low O/C compounds ("protein-like" fluorescence) and that the composition gradually evolves to aromatic, high O/C compounds ("humic-/fulvic-like" fluorescence) as a function of distance downgradient from the oil body. Finally, a direct, positive correlation between optical properties and Microtox acute toxicity assays demonstrates the utility of these combined techniques in assessing the spatial and temporal natural attenuation and toxicity of the DOMHC in petroleum-impacted groundwater systems.

Rates of As and Trace-Element Mobilization Caused by Fe Reduction in Mixed BTEX-Ethanol Experimental Plumes.

Biodegradation of organic matter, including petroleum-based fuels and biofuels, can create undesired secondary water-quality effects. Trace elements, especially arsenic (As), have strong adsorption affinities for Fe(III) (oxyhydr)-oxides and can be released to groundwater during Fe-reducing biodegradation. We investigated the mobilization of naturally occurring As, cobalt (Co), chromium (Cr), and nickel (Ni) from wetland sediments caused by the introduction of benzene, toluene, ethylbenzene, and xylenes (BTEX) and ethanol mixtures under iron- and nitrate-reducing conditions, using in situ push-pull tests. When BTEX alone was added, results showed simultaneous onset and similar rates of Fe reduction and As mobilization. In the presence of ethanol, the maximum rates of As release and Fe reduction were higher, the time to onset of reaction was decreased, and the rates occurred in multiple stages that reflected additional processes. The concentration of As increased from <1 µg/L to a maximum of 99 µg/L, exceeding the 10 µg/L limit for drinking water. Mobilization of Co, Cr, and Ni was observed in association with ethanol biodegradation but not with BTEX. These results demonstrate the potential for trace-element contamination of drinking water during biodegradation and highlight the importance of monitoring trace elements at natural and enhanced attenuation sites.

An integrated Markov chain Monte Carlo algorithm for upscaling hydrological and geochemical parameters from column to field scale.

Predicting and controlling the concentrations of redox-sensitive elements are primary concerns for environmental remediation of contaminated sites. These predictions are complicated by dynamic flow processes as hydrologic variability is a governing control on conservative and reactive chemical concentrations. Subsurface heterogeneity in the form of layers and lenses further complicates the flow dynamics of the system impacting chemical concentrations including redox-sensitive elements. In response to these complexities, this study investigates the role of heterogeneity and hydrologic processes in an effective parameter upscaling scheme from the column to the landfill scale. We used a Markov chain Monte Carlo (MCMC) algorithm to derive upscaling coefficients for hydrological and geochemical parameters, which were tested for variations across heterogeneous systems (layers and lenses) and interaction of flow processes based on the output uncertainty of dominant biogeochemical concentrations at the Norman Landfill site, a closed municipal landfill with prevalent organic and trace metal contamination. The results from MCMC analysis indicated that geochemical upscaling coefficients based on effective concentration ratios incorporating local heterogeneity across layered and lensed systems produced better estimates of redox-sensitive biogeochemistry at the field scale. MCMC analysis also suggested that inclusion of hydrological parameters in the upscaling scheme reduced the output uncertainty of effective mean geochemical concentrations by orders of magnitude at the Norman Landfill site. This was further confirmed by posterior density plots of the scaling coefficients that revealed unimodal characteristics when only geochemical processes were involved, but produced multimodal distributions when hydrological parameters were included. The multimodality again suggests the effect of heterogeneity and lithologic variability on the distribution of redox-sensitive elements at the Norman Landfill site.

Enhanced biogeochemical cycling and subsequent reduction of hydraulic conductivity associated with soil-layer interfaces in the vadose zone.

Biogeochemical dynamics in the vadose zone are poorly understood due to the transient nature of chemical and hydrologic conditions but are nonetheless critical to understanding chemical fate and transport. This study explored the effects of a soil layer on linked geochemical, hydrological, and microbiological processes. Three laboratory soil columns were constructed: a homogenized medium-grained sand, a homogenized organic-rich loam, and a sand-over-loam layered column. Upward and downward infiltration of water was evaluated during experiments to simulate rising water table and rainfall events, respectively. In situ collocated probes measured soil water content, matric potential, and Eh. Water samples collected from the same locations were analyzed for Br, Cl, NO, SO, NH, Fe, and total sulfide. Compared with homogeneous columns, the presence of a soil layer altered the biogeochemistry and water flow of the system considerably. Enhanced biogeochemical cycling was observed in the layered column over the texturally homogeneous soil columns. Enumerations of iron- and sulfate-reducing bacteria showed 1 to 2 orders of magnitude greater community numbers in the layered column. Mineral and soil aggregate composites were most abundant near the soil-layer interface, the presence of which likely contributed to an observed order-of-magnitude decrease in hydraulic conductivity. These findings show that quantifying coupled hydrologic-biogeochemical processes occurring at small-scale soil interfaces is critical to accurately describing and predicting chemical changes at the larger system scale. These findings also provide justification for considering soil layering in contaminant fate and transport models because of its potential to increase biodegradation or to slow the rate of transport of contaminants.

In situ rates of sulfate reduction in response to geochemical perturbations.

Rates of in situ microbial sulfate reduction in response to geochemical perturbations were determined using Native Organism Geochemical Experimentation Enclosures (NOGEEs), a new in situ technique developed to facilitate evaluation of controls on microbial reaction rates. NOGEEs function by first trapping a native microbial community in situ and then subjecting it to geochemical perturbations through the introduction of various test solutions. On three occasions, NOGEEs were used at the Norman Landfill research site in Norman, Oklahoma, to evaluate sulfate-reduction rates in wetland sediments impacted by landfill leachate. The initial experiment, in May 2007, consisted of five introductions of a sulfate test solution over 11 d. Each test stimulated sulfate reduction with rates increasing until an apparent maximum was achieved. Two subsequent experiments, conducted in October 2007 and February 2008, evaluated the effects of concentration on sulfate-reduction rates. Results from these experiments showed that faster sulfate-reduction rates were associated with increased sulfate concentrations. Understanding variability in sulfate-reduction rates in response to perturbations may be an important factor in predicting rates of natural attenuation and bioremediation of contaminants in systems not at biogeochemical equilibrium.

Inverse estimation of parameters for multidomain flow models in soil columns with different macropore densities.

Soil and crop management practices have been found to modify soil structure and alter macropore densities. An ability to accurately determine soil hydraulic parameters and their variation with changes in macropore density is crucial for assessing potential contamination from agricultural chemicals. This study investigates the consequences of using consistent matrix and macropore parameters in simulating preferential flow and bromide transport in soil columns with different macropore densities (no macropore, single macropore, and multiple macropores). As used herein, the term"macropore density" is intended to refer to the number of macropores per unit area. A comparison between continuum-scale models including single-porosity model (SPM), mobile-immobile model (MIM), and dual-permeability model (DPM) that employed these parameters is also conducted. Domain-specific parameters are obtained from inverse modeling of homogeneous (no macropore) and central macropore columns in a deterministic framework and are validated using forward modeling of both low-density (3 macropores) and high-density (19 macropores) multiple-macropore columns. Results indicate that these inversely modeled parameters are successful in describing preferential flow but not tracer transport in both multiple-macropore columns. We believe that lateral exchange between matrix and macropore domains needs better accounting to efficiently simulate preferential transport in the case of dense, closely spaced macropores. Increasing model complexity from SPM to MIM to DPM also improved predictions of preferential flow in the multiple-macropore columns but not in the single-macropore column. This suggests that the use of a more complex model with resolved domain-specific parameters is recommended with an increase in macropore density to generate forecasts with higher accuracy.

Determination of dominant biogeochemical processes in a contaminated aquifer-wetland system using multivariate statistical analysis.

Determining the processes governing aqueous biogeochemistry in a wetland hydrologically linked to an underlying contaminated aquifer is challenging due to the complex exchange between the systems and their distinct responses to changes in precipitation, recharge, and biological activities. To evaluate temporal and spatial processes in the wetland-aquifer system, water samples were collected using cm-scale multi-chambered passive diffusion samplers (peepers) to span the wetland-aquifer interface over a period of 3 yr. Samples were analyzed for major cations and anions, methane, and a suite of organic acids resulting in a large dataset of over 8000 points, which was evaluated using multivariate statistics. Principal component analysis (PCA) was chosen with the purpose of exploring the sources of variation in the dataset to expose related variables and provide insight into the biogeochemical processes that control the water chemistry of the system. Factor scores computed from PCA were mapped by date and depth. Patterns observed suggest that (i) fermentation is the process controlling the greatest variability in the dataset and it peaks in May; (ii) iron and sulfate reduction were the dominant terminal electron-accepting processes in the system and were associated with fermentation but had more complex seasonal variability than fermentation; (iii) methanogenesis was also important and associated with bacterial utilization of minerals as a source of electron acceptors (e.g., barite BaSO(4)); and (iv) seasonal hydrological patterns (wet and dry periods) control the availability of electron acceptors through the reoxidation of reduced iron-sulfur species enhancing iron and sulfate reduction.

Analysis of recharge-induced geochemical change in a contaminated aquifer.

Recharge events that deliver electron acceptors such as O2, NO3, SO4, and Fe3+ to anaerobic, contaminated aquifers are likely important for natural attenuation processes. However, the specific influence of recharge on (bio)geochemical processes in ground water systems is not well understood. The impact of a moderate-sized recharge event on ground water chemistry was evaluated at a shallow, sandy aquifer contaminated with waste fuels and chlorinated solvents. Multivariate statistical analyses coupled with three-dimensional visualization were used to analyze ground water chemistry data (including redox indicators, major ions, and physical parameters) to reveal associations between chemical parameters and to infer processes within the ground water plume. Factor analysis indicated that dominant chemical associations and their interpreted processes (anaerobic and aerobic microbial processes, mineral precipitation/dissolution, and temperature effects) did not change significantly after the spring recharge event of 2000. However, the relative importance of each of these processes within the plume changed. After the recharge event, the overall importance of aerobic processes increased from the fourth to the second most important factor, representing the variability within the data set. The anaerobic signatures became more complex, suggesting that zones with multiple terminal electron-accepting processes (TEAPs) likely occur in the same water mass. Three-dimensional visualization of well clusters showed that water samples with similar chemical associations occurred in distinct water masses within the aquifer. Water mass distinctions were not based on dominant TEAPs, suggesting that the recharge effects on TEAPs occurred primarily at the interface between infiltrating recharge water and the aquifer.

Spatial and temporal changes in microbial community structure associated with recharge-influenced chemical gradients in a contaminated aquifer.

In a contaminated water-table aquifer, we related microbial community structure on aquifer sediments to gradients in 24 geochemical and contaminant variables at five depths, under three recharge conditions. Community amplified ribsosomal DNA restriction analysis (ARDRA) using universal 16S rDNA primers and denaturing gradient gel electrophoresis (DGGE) using bacterial 16S rDNA primers indicated: (i). communities in the anoxic, contaminated central zone were similar regardless of recharge; (ii). after recharge, communities at greatest depth were similar to those in uncontaminated zones; and (iii). after extended lack of recharge, communities at upper and lower aquifer margins differed from communities at the same depths on other dates. General aquifer geochemistry was as important as contaminant or terminal electron accepting process (TEAP) chemistry in discriminant analysis of community groups. The Shannon index of diversity (H) and the evenness index (E), based on DGGE operational taxonomic units (OTUs), were statistically different across community groups and aquifer depths. Archaea or sulphate-reducing bacteria 16S rRNA abundance was not clearly correlated with TEAP chemistry indicative of methanogenesis or sulphate reduction. Eukarya rRNA abundance varied by depth and date from 0 to 13% of the microbial community. This contaminated aquifer is a dynamic ecosystem, with complex interactions between physical, chemical and biotic components, which should be considered in the interpretation of aquifer geochemistry and in the development of conceptual or predictive models for natural attenuation or remediation.

Evaluating behavior of oxygen, nitrate, and sulfate during recharge and quantifying reduction rates in a contaminated aquifer.

This study evaluates the biogeochemical changes that occur when recharge water comes in contact with a reduced aquifer. It specifically addresses (1) which reactions occur in situ, (2) the order in which these reactions will occur if terminal electron acceptors (TEAs) are introduced simultaneously, (3) the rates of these reactions, and (4) the roles of the aqueous and solid-phase portions of the aquifer. Recharge events of waters containing various combinations of O2, NO3, and SO4 were simulated at a shallow sandy aquifer contaminated with waste fuels and chlorinated solvents using modified push-pull tests to quantify rates. In situ rate constants for aerobic respiration (14.4 day(-1)), denitrification (5.04-7.44 day(-1)), and sulfate reduction (4.32-6.48 day(-1)) were estimated. Results show that when introduced together, NO3 and SO4 can be consumed simultaneously at similar rates. To distinguish the role of aqueous phase from that of the solid phase of the aquifer, groundwater was extracted, amended with NO3 and SO4, and monitored overtime. Results indicate that neither NO3 nor SO4 was reduced during the course of the aqueous-phase study, suggesting that NO3 and SO4 can behave conservatively in highly reduced water. It is clear that sediments and their associated microbial communities are important in driving redox reactions.