Richard T. Amos

Vadose zone attenuation of organic compounds at a crude oil spill site - interactions between biogeochemical reactions and multicomponent gas transport.

Contaminant attenuation processes in the vadose zone of a crude oil spill site near Bemidji, MN have been simulated with a reactive transport model that includes multicomponent gas transport, solute transport, and the most relevant biogeochemical reactions. Dissolution and volatilization of oil components, their aerobic and anaerobic degradation coupled with sequential electron acceptor consumption, ingress of atmospheric O(2), and the release of CH(4) and CO(2) from the smear zone generated by the floating oil were considered. The focus of the simulations was to assess the dynamics between biodegradation and gas transport processes in the vadose zone, to evaluate the rates and contributions of different electron accepting processes towards vadose zone natural attenuation, and to provide an estimate of the historical mass loss. Concentration distributions of reactive (O(2), CH(4), and CO(2)) and non-reactive (Ar and N(2)) gases served as key constraints for the model calibration. Simulation results confirm that as of 2007, the main degradation pathway can be attributed to methanogenic degradation of organic compounds in the smear zone and the vadose zone resulting in a contaminant plume dominated by high CH(4) concentrations. In accordance with field observations, zones of volatilization and CH(4) generation are correlated to slightly elevated total gas pressures and low partial pressures of N(2) and Ar, while zones of aerobic CH(4) oxidation are characterized by slightly reduced gas pressures and elevated concentrations of N(2) and Ar. Diffusion is the most significant transport mechanism for gases in the vadose zone; however, the simulations also indicate that, despite very small pressure gradients, advection contributes up to 15% towards the net flux of CH(4), and to a more limited extent to O(2) ingress. Model calibration strongly suggests that transfer of biogenically generated gases from the smear zone provides a major control on vadose zone gas distributions and vadose zone carbon balance. Overall, the model was successful in capturing the complex interactions between biogeochemical reactions and multicomponent gas transport processes. However, despite employing a process-based modeling approach, honoring observed parameter ranges, and generally obtaining good agreement between field observations and model simulations, accurate quantification of natural attenuation rates remains difficult. The modeling results are affected by uncertainties regarding gas phase saturations, tortuosities, and the magnitude of CH(4) and CO(2) flux from the smear zone. These findings highlight the need to better delineate gas fluxes at the model boundaries, which will help constrain contaminant degradation rates, and ultimately source zone longevity.

The structure of alkali germanophosphate glasses by Raman spectroscopy

Abstract A series of alkali germanophosphate glasses ((R2O)x(GeO2:P2O5)1−x where R=Na, K and Rb) with variable GeO2:P2O5 ratios (8:1, 6:1, 4:1 and 2:1) have been investigated using Raman spectroscopy. The glass network may be treated as being made up of separate germanate and phosphate components. Addition of alkali cations indicates that the alkalis preferentially modify the phosphate part of the network. Network depolymerization occurs by formation of Q2 and Q1 PO4 tetrahedra. At high Ge:P ratios the glasses exhibit a density anomaly. The anomaly is attributed to the formation of small three-membered GeO4 rings. However, alkali cation size and mass are contributing factors to the shape of the density anomaly in the K2O and Rb2O containing glasses. Depolymerization of the 2:1 alkali germanophosphate glasses is predominantly by formation of Q2 PO4 tetrahedra. At low Ge:P ratios (2:1) density trends are linear. No ring transition is observed in these glasses. Furthermore, the presence of [5]Ge cannot be ruled out, but if present, does not play a major role in generating a ‘germanate anomaly’.

Evidence for iron‐mediated anaerobic methane oxidation in a crude oil‐contaminated aquifer

In a methanogenic crude oil contaminated aquifer near Bemidji, Minnesota, the decrease in dissolved CH(4) concentrations along the groundwater flow path, along with the positive shift in δ(13) C(CH) (4) and negative shift in δ(13) C(DIC) , is indicative of microbially mediated CH(4) oxidation. Calculations of electron acceptor transport across the water table, through diffusion, recharge, and the entrapment and release of gas bubbles, suggest that these processes can account for at most 15% of the observed total reduced carbon oxidation, including CH(4) . In the anaerobic plume, the characteristic Fe(III)-reducing genus Geobacter was the most abundant of the microbial groups tested, and depletion of labile sediment iron is observed over time, confirming that reduced carbon oxidation coupled to iron reduction is an important process. Electron mass balance calculations suggest that organic carbon sources in the aquifer, BTEX and non-volatile dissolved organic carbon, are insufficient to account for the loss in sediment Fe(III), implying that CH(4) oxidation may also be related to Fe(III) reduction. The results support a hypothesis of Fe(III)-mediated CH(4) oxidation in the contaminated aquifer.

34S/32S fractionation during sulfate reduction in groundwater treatment systems: reactive transport modeling.

Isotope ratio measurements provide a tool for indicating the relative significance of biogeochemical reactions and for constraining estimates of the extent and rate of reactions in passive treatment systems. In this paper, the reactive transport model MIN3P is used to evaluate sulfur isotope fractionation in column experiments designed to simulate treatment of contaminated water by microbially mediated sulfate reduction occurring within organic carbon-based and iron and carbon-based permeable reactive barriers. A mass dependent fractionation model was used to determine reaction rates for 32S and 34S compounds during reduction, precipitation, and dissolution reactions and to track isotope-dependent mass transfer during SO4 removal. The δ34S values obtained from the MIN3P model were similar to those obtained from the Rayleigh equation, indicating that there was not a significant difference between the conceptual models. Differences between the MIN3P derived α value and the Rayleigh equation derived value were attributed to minor changes in the dissolution and precipitation rate of gypsum and mathematical differences in the fitting models. The results indicated that the prediction of δ34S was fairly insensitive to differences in the fractionation factor at the concentration ranges measured in the current study. However, more significant differences would be expected at low sulfate conditions.

Modeling gas formation and mineral precipitation in a granular iron column.

In granular iron permeable reactive barriers (PRBs), hydrogen gas formation, entrapment and release of gas bubbles, and secondary mineral precipitation have been known to affect the permeability and reactivity. The multicomponent reactive transport model MIN3P was enhanced to couple gas formation and release, secondary mineral precipitation, and the effects of these processes on hydraulic properties and iron reactivity. The enhanced model was applied to a granular iron column, which was studied for the treatment of trichloroethene (TCE) in the presence of dissolved CaCO(3). The simulation reasonably reproduced trends in gas formation, secondary mineral precipitation, permeability changes, and reactivity changes observed over time. The simulation showed that the accumulation of secondary minerals reduced the reactivity of the granular iron over time, which in turn decreased the rate of mineral accumulation, and also resulted in a gradual decrease in gas formation over time. This study provides a quantitative assessment of the evolving nature of geochemistry and permeability, resulting from coupled processes of gas formation and mineral precipitation, which leads to a better understanding of the processes controlling the granular iron reactivity, and represents an improved method for incorporating these factors into the design of granular iron PRBs.

Influence of uranyl speciation and iron oxides on uranium biogeochemical redox reactions

Uranium is a pollutant of concern to both human and ecosystem health. Uraniums redox state often dictates its partitioning between the aqueous- and solid-phases, and thus controls its dissolved concentration and, coupled with groundwater flow, its migration within the environment. In anaerobic environments, the more oxidized and mobile form of uranium (UO2 2+ and associated species) may be reduced, directly or indirectly, by microorganisms to U(IV) with subsequent precipitation of UO2. However, various factors within soils and sediments may limit biological reduction of U(VI), inclusive of alterations in U(VI) speciation and competitive electron acceptors. Here we elucidate the impact of U(VI) speciation on the extent and rate of reduction with specific emphasis on speciation changes induced by dissolved Ca, and we examine the impact of Fe(III) (hydr)oxides (ferrihydrite, goethite and hematite) varying in free energies of formation on U reduction. The amount of uranium removed from solution during 100 h of incubation with S. putrefaciens was 77% with no Ca or ferrihydrite present but only 24% (with ferrihydrite) and 14% (no ferrihydrite) were removed for systems with 0.8 mM Ca. Imparting an important criterion on uranium reduction, goethite and hematite decrease the dissolved concentration of calcium through adsorption and thus tend to diminish the effect of calcium on uranium reduction. Dissimilatory reduction of Fe(III) and U(VI) can proceed through different enzyme pathways, even within a single organism, thus providing a potential second means by which Fe(III) bearing minerals may impact U(VI) reduction. We quantify rate coefficients for simultaneous dissimilatory reduction of Fe(III) and U(VI) in systems varying in Ca concentration (0 to 0.8 mM), and using a mathematical construct implemented with the reactive transport code MIN3P, we reveal the predominant influence of uranyl speciation, specifically the formation of uranyl-calcium-carbonato complexes, and ferrihydrite on the rate and extent of uranium reduction in complex geochemical systems.

Reactive transport modeling of geochemical controls on secondary water quality impacts at a crude oil spill site near Bemidji, MN

Anaerobic biodegradation of organic amendments and contaminants in aquifers can trigger secondary water quality impacts that impair groundwater resources. Reactive transport models help elucidate how diverse geochemical reactions control the spatiotemporal evolution of these impacts. Using extensive monitoring data from a crude oil spill site near Bemidji, Minnesota (USA), we implemented a comprehensive model that simulates secondary plumes of depleted dissolved O2 and elevated concentrations of Mn2+, Fe2+, CH4, and Ca2+ over a two-dimensional cross section for 30 years following the spill. The model produces observed changes by representing multiple oil constituents and coupled carbonate and hydroxide chemistry. The model includes reactions with carbonates and Fe and Mn mineral phases, outgassing of CH4 and CO2 gas phases, and sorption of Fe, Mn, and H+. Model results demonstrate that most of the carbon loss from the oil (70%) occurs through direct outgassing from the oil source zone, greatly limiting the amount of CH4 cycled down-gradient. The vast majority of reduced Fe is strongly attenuated on sediments, with most (91%) in the sorbed form in the model. Ferrous carbonates constitute a small fraction of the reduced Fe in simulations, but may be important for furthering the reduction of ferric oxides. The combined effect of concomitant redox reactions, sorption, and dissolved CO2 inputs from source-zone degradation successfully reproduced observed pH. The model demonstrates that secondary water quality impacts may depend strongly on organic carbon properties, and impacts may decrease due to sorption and direct outgassing from the source zone.

Methane oxidation in a crude oil contaminated aquifer: Delineation of aerobic reactions at the plume fringes

High resolution direct-push profiling over short vertical distances was used to investigate CH(4) attenuation in a petroleum contaminated aquifer near Bemidji, Minnesota. The contaminant plume was delineated using dissolved gases, redox sensitive components, major ions, carbon isotope ratios in CH(4) and CO(2), and the presence of methanotrophic bacteria. Sharp redox gradients were observed near the water table. Shifts in δ(13)C(CH4) from an average of -57.6‰ (±1.7‰) in the methanogenic zone to -39.6‰ (±8.7‰) at 105m downgradient, strongly suggest CH(4) attenuation through microbially mediated degradation. In the downgradient zone the aerobic/anaerobic transition is up to 0.5m below the water table suggesting that transport of O(2) across the water table is leading to aerobic degradation of CH(4) at this interface. Dissolved N(2) concentrations that exceeded those expected for water in equilibrium with the atmosphere indicated bubble entrapment followed by preferential stripping of O(2) through aerobic degradation of CH(4) or other hydrocarbons. Multivariate and cluster analysis were used to distinguish between areas of significant bubble entrapment and areas where other processes such as the infiltration of O(2) rich recharge water were important O(2) transport mechanisms.

A mass balance approach to investigating geochemical controls on secondary water quality impacts at a crude oil spill site near Bemidji, MN

Secondary water quality impacts can result from a broad range of coupled reactions triggered by primary groundwater contaminants. Data from a crude-oil spill research site near Bemidji, MN provide an ideal test case for investigating the complex interactions controlling secondary impacts, including depleted dissolved oxygen and elevated organic carbon, inorganic carbon, CH4, Mn, Fe, and other dissolved ions. To better understand these secondary impacts, this study began with an extensive data compilation of various data types, comprising aqueous, sediment, gas, and oil phases, covering a 260m cross-sectional domain over 30years. Mass balance calculations are used to quantify pathways that control secondary components, by using the data to constrain the sources and sinks for the important redox processes. The results show that oil constituents other than BTEX (benzene, toluene, ethylbenzene, o-, m- and p-xylenes), including n-alkanes and other aromatic compounds, play significant roles in plume evolution and secondary water quality impacts. The analysis underscores previous results on the importance of non-aqueous phases. Over 99.9% of the Fe(2+) plume is attenuated by immobilization on sediments as Fe(II) and 85-95% of the carbon biodegradation products are outgassed. Gaps identified in carbon and Fe mass balances and in pH buffering mechanisms are used to formulate a new conceptual model. This new model includes direct out-gassing of CH4 and CO2 from organic carbon biodegradation, dissolution of directly produced CO2, and sorption with H(+) exchange to improve pH buffering. The identification of these mechanisms extends understanding of natural attenuation of potential secondary impacts at enhanced reductive dechlorination sites, particularly for reduced Fe plumes, produced CH4, and pH perturbations.

Investigating dominant processes in ZVI permeable reactive barriers using reactive transport modeling.

The reactive and hydraulic efficacy of zero valent iron permeable reactive barriers (ZVI PRBs) is strongly affected by geochemical composition of the groundwater treated. An enhanced version of the geochemical simulation code MIN3P was applied to simulate dominating processes in chlorinated hydrocarbons (CHCs) treating ZVI PRBs including geochemical dependency of ZVI reactivity, gas phase formation and a basic formulation of degassing. Results of target oriented column experiments with distinct chemical conditions (carbonate, calcium, sulfate, CHCs) were simulated to parameterize the model. The simulations demonstrate the initial enhancement of anaerobic iron corrosion due to carbonate and long term inhibition by precipitates (chukanovite, siderite, iron sulfide). Calcium was shown to enhance long term corrosion due to competition for carbonate between siderite, chukanovite, and aragonite, with less inhibition of iron corrosion by the needle like aragonite crystals. Application of the parameterized model to a field site (Bernau, Germany) demonstrated that temporarily enhanced groundwater carbonate concentrations caused an increase in gas phase formation due to the acceleration of anaerobic iron corrosion.