Tracing Molybdenum Attenuation in Mining Environments Using Molybdenum Stable Isotopes.

Molybdenum contamination is a concern in mining regions worldwide. Better understanding of processes controlling Mo mobility in mine wastes is critical for assessing potential impacts and developing water-quality management strategies associated with this element. Here, we used Mo stable isotope (δ98/95Mo) analyses to investigate geochemical controls on Mo mobility within a tailings management facility (TMF) featuring oxic and anoxic environments. These isotopic analyses were integrated with X-ray absorption spectroscopy, X-ray diffraction, Raman spectroscopy, transmission electron microscopy, and aqueous chemical data. Dissolved Mo concentrations were inversely correlated with δ98/95Mo values such that enrichment of heavy Mo isotopes in solution reflected attenuation processes. Inner-sphere complexation of Mo(VI) with ferrihydrite was the primary driver of Mo removal and was accompanied by a ca. 1‰ isotope fractionation. Limited Mo attenuation and isotope fractionation were observed in Fe(II)- and Mo-rich anoxic TMF seepage, while attenuation and isotope fractionation were greatest during discharge and oxidation of this seepage after discharge into a pond where Fe-(oxyhydr)oxide precipitation promoted Mo sorption. Overall, this study highlights the role of sorption onto Fe-(oxyhydr)oxides in attenuating Mo in oxic environments, a process which can be traced by Mo isotope analyses.

Two-dimensional modeling of CO2 mineral trapping through the oxalate­carbonate pathway: Influence of the root system model.

Two-dimensional reactive transport models, one with a simplified root system and the other accounting for dynamically evolving root architecture, were constructed to examine the influence of model complexity on capturing the effect of soil-root dynamics relating to the Oxalate Carbonate Pathway (OCP) of the Iroko tree over 170 years. Oxidation of oxalate from fallen tree tissue by soil bacteria enables local soil pH increase, leading to the sequestration of atmospheric carbon in carbonate minerals (calcite) in the shallow soil surrounding the tree. Simulations of both root models corroborate previous one-dimensional models of the OCP focused on Ca and C mass balance, where high weathering rates of Ca-containing silicate minerals in bedrock, along with contributions from groundwater, provided sufficient Ca for precipitation of observed quantities of calcite. Both simulations demonstrate the development of a distinct high pH zone where oxalate is oxidized, Ca accumulates, and calcite precipitates (OCP zone); and a low pH zone where roots collect Ca, later returned to the top soil as calcium oxalate (Total Root Extent/TRE zone) via litterfall. While the extent of OCP zone development near the ground surface was very similar between simulations, differences in localized root water uptake between the two approaches resulted in variation in water and solute transport and influenced the geometry of the OCP zone at depth, with implications for calcite precipitation in the soil. Trends in CO2 and O2 partial pressures in the OCP zone were mirrored in the TRE zone, suggesting linkage between the two zones with regard to gas transport. Near the end of the tree's lifespan, results indicate that soil permeability decreases due to calcite precipitation may limit O2 ingress and availability in the shallow soil, while trapping CO2 released from the oxidation of organics in the shallow soil, with implications for the long-term sustainability of the OCP itself.

What Can Groundwater Monitoring Tell Us About Gas Migration? A Numerical Modeling Study.

Groundwater monitoring to measure a variety of indicator parameters including dissolved gas concentrations, total dissolved gas pressure (TDGP), and redox indicators is commonly used to evaluate the impacts of gas migration (GM) from energy development in shallow aquifer systems. However, these parameters can be challenging to interpret due to complex free-phase gas source architecture, multicomponent partitioning, and biogeochemical reactions. A series of numerical simulations using a gas flow model and a reactive transport model were conducted to delineate the anticipated evolution of indicator parameters following GM in an aquifer under a variety of physical and biogeochemical conditions. The simulations illustrate how multicomponent mass transfer processes and biogeochemical reactions create unexpected spatial and temporal variations in several analytes. The results indicate that care must be taken when interpreting measured indicator parameters including dissolved hydrocarbon concentrations and TDGP, as the presence of dissolved gases in background groundwater and biogeochemical processes can cause potentially misleading conclusions about the impact of GM. Based on the consideration of multicomponent gas partitioning in this study, it is suggested that dissolved background gases such as N2 and Ar can provide valuable insights on the presence, longevity and fate of free-phase natural gas in aquifer systems. Overall, these results contribute to developing a better understanding of indicators for GM in groundwater, which will aid the planning of future monitoring networks and subsequent data interpretation.

Travel time-based modelling of nitrate reduction in a fractured limestone aquifer by pyrite and iron carbonates under pore size limitation.

We investigate denitrification in a ferric iron-containing fractured micritic limestone aquifer (Triassic Upper Muschelkalk) in south-west Germany by numerical simulations. Low porosity values (average value of 1%), partly small pore sizes of the rock matrix (~ 0.1 µm), and thus potential absence of microbial activity in the rock matrix suggest that denitrification is taking place solely in the fracture. A key question is whether the nitrate reduction derived from groundwater observations at 25 locations in the study area can be explained by a model that restricts microbial denitrification to the fractures. A travel time-based reactive transport model is developed to efficiently simulate long-term nitrate reduction on the catchment scale. The model employs a 2-D numerical reaction model describing the fracture-rock matrix system and parametric travel time distributions. The role of (i) biotic and abiotic iron oxidation, (ii) the type and amount of iron bearing minerals, and (iii) mass transfer between matrix and fracture are investigated. The simulations show that pyrite and siderite (used as surrogate for iron carbonates) together as a source of electron donors provide enough reduction potential to decrease the nitrate concentrations as observed in the field. This confirms the hypothesis that diffusion-controlled mass transfer of electron donors from the matrix to the fracture is sufficient to establish considerable denitrification in the fracture. Uncertainty in modelled concentrations is demonstrated as a result of both the geochemical aquifer properties and the unknown shape of travel time distributions.

The nature of gas production patterns associated with methanol degradation in natural aquifer sediments: A microcosm study.

With growing global use of methanol as a fuel additive and extensive use in other industrial processes, there is the potential for unintended release and spills into soils and aquifers. In these subsurface systems it is likely that methanol will be readily biodegraded; however, degradation may lead to the production of by-products, most importantly methane possibly resulting in explosion hazards and volatile fatty acids (VFAs) causing aesthetic issues for groundwater. In this study, the formation of these potentially harmful by-products due to methanol biodegradation was investigated in natural sand and silt sediments using microcosms inoculated with neat methanol (100%) ranging in concentration from 100 to 100,000 ppm. To assess the rate of degradation and by-product formation, water and headspace samples were collected and analyzed for methanol, volatile fatty acids (VFAs, including acetic, butyric, and propionic acid), cation (metal) concentrations (Al, Ca, Fe, K, Mg, Mn and Na), microbial community structure and activity, headspace pressure, gas composition (CH4, CO2, O2 and N2), and compound specific isotopes. Methanol was completely biodegraded in sand and silt up to concentrations of 1000 ppm and 10,000 ppm, respectively. Degradation was initially aerobic, consuming oxygen (O2) and producing carbon dioxide (CO2). When O2 was depleted, the microcosms became anaerobic and a lag in methanol degradation occurred (ranging from 41 to 87 days). Following this lag, methanol was preferentially degraded to acetate, coupled with CO2 reduction. Microcosms with high methanol concentrations (10,000 ppm) were driven further down the redox ladder and exhibited fermentation, leading to concurrent acetate and methane (CH4) generation. In all cases acetate was an intermediate product, further degraded to the final products of CH4 and CO2. Carbonates present in the microcosm sediments helped buffer VFA acidification and replenished CO2. Methane generation in the anaerobic microcosms was short-lived, but temporarily reached high rates up to 13 mg kg-1 day-1. Under the conditions of these experiments, methanol degradation occurred rapidly, after initial lag periods, which were a function of methanol concentration and sediment type. Our experiment also showed that methanol degradation and associated methane production can occur in a stepwise fashion.

CO2-efflux measurements for evaluating source zone natural attenuation rates in a petroleum hydrocarbon contaminated aquifer.

In order to gain regulatory approval for source zone natural attenuation (SZNA) at hydrocarbon-contaminated sites, knowledge regarding the extent of the contamination, its tendency to spread, and its longevity is required. However, reliable quantification of biodegradation rates, an important component of SZNA, remains a challenge. If the rate of CO(2) gas generation associated with contaminant degradation can be determined, it may be used as a proxy for the overall rate of subsurface biodegradation. Here, the CO(2)-efflux at the ground surface is measured using a dynamic closed chamber (DCC) method to evaluate whether this technique can be used to assess the areal extent of the contaminant source zone and the depth-integrated rate of contaminant mineralization. To this end, a field test was conducted at the Bemidji, MN, crude oil spill site. Results indicate that at the Bemidji site the CO(2)-efflux method is able to both delineate the source zone and distinguish between the rates of natural soil respiration and contaminant mineralization. The average CO(2)-efflux associated with contaminant degradation in the source zone is estimated at 2.6 µmol m(-2) s(-1), corresponding to a total petroleum hydrocarbon mineralization rate (expressed as C(10)H(22)) of 3.3 g m(-2) day(-1).

Transport implications resulting from internal redistribution of arsenic and iron within constructed soil aggregates.

Soils are an aggregate-based structured media that have a multitude of pore domains resulting in varying degrees of advective and diffusive solute and gas transport. Consequently, a spectrum of biogeochemical processes may function at the aggregate scale that collectively, and coupled with solute transport, determine element cycling in soils and sediments. To explore how the physical structure impacts biogeochemical processes influencing the fate and transport of As, we examined temporal changes in speciation and distribution of As and Fe within constructed aggregates through experimental measurement and reactive transport simulations. Spherical aggregates were made with As(V)-bearing ferrihydrite-coated sand inoculated with Shewanella sp. ANA-3; aerated solute flow around the aggregate was then induced. Despite the aerated aggregate exterior, where As(V) and ferrihydrite persist as the dominant species, anoxia develops within the aggregate interior. As a result, As and Fe redox gradients emerge, and the proportion of As(III) and magnetite increases toward the aggregate interior. Arsenic(III) and Fe(II) produced in the interior migrate toward the aggregated exterior and result in coaccumulation of As and Fe(III) proximal to preferential flow paths as a consequence of oxygenic precipitation. The oxidized rind of aggregates thus serves as a barrier to As release into advecting pore-water, but also leads to be a buildup of this hazardous element at preferential flow boundaries that could be released upon shifting geochemical conditions.

Quantifying natural source zone depletion at petroleum hydrocarbon contaminated sites: A comparison of 14C methods.

Surficial CO2 efflux surveys have been used to delineate hydrocarbon source zones in contaminated aquifers and provide estimates of hydrocarbon biodegradation rates. This approach requires distinguishing between CO2 derived from petroleum degradation and CO2 produced from natural soil respiration. To this end, radiocarbon has been used to differentiate between 14C-depleted CO2 from hydrocarbon degradation and 14C-enriched CO2 from natural soil respiration to effectively quantify the contribution of each source to total CO2 efflux, and by deduction natural source zone depletion (NSZD) rates. In this study, a systematic method comparison has been conducted to evaluate available approaches for collecting CO2 gas samples for radiocarbon analysis used to correct total CO2 efflux measurements for quantifying natural source zone depletion rates. Gas samples for radiocarbon analysis were sampled from (i) dynamic closed chambers (located at ground surface), (ii) static chambers (also at ground surface), (iii) shallow soil gas probes (0.3 m bgs), and (iv) soil gas monitoring wells (~0.6 m below ground surface) during a CO2 efflux survey conducted at the site of a historical pipeline rupture near Bemidji, MN. The mean fraction of radiocarbon (F14C) obtained from samples overlying the source zone were (i) 0.93 ± 0.01, (ii) 0.73 ± 0.03, (iii) 0.71 ± 0.04, and (iv) 0.41 ± 0.06, for the four methods respectively. These F14C values were used to apportion total CO2 efflux measurements into contributions of contaminant-derived CO2 efflux and natural soil respiration to evaluate natural source zone depletion processes. Results suggest that the method of radiocarbon sampling has a significant effect on the calculated fraction of the CO2 efflux originating from contaminant-related soil respiration, with contributions ranging between 27% and 59% of total soil respiration. Results indicate that radiocarbon sampled from static chambers and shallow soil gas probes methods offer the best compromise between CO2 sample yield and sample representativeness, providing the most reliable estimates of CO2 effluxes originating from contaminant degradation. However, the results also show that for this study, all methods agree within a factor of <2.3 regarding the inferred NSZD rates.

Impacts of water table fluctuations on actual and perceived natural source zone depletion rates.

A viable means of quantifying the rate of natural source zone depletion (NSZD) at hydrocarbon contaminated sites is by the measurement of carbon dioxide (CO2) and methane (CH4) effluxes at the surface. This methodology assumes that gas effluxes are reflective of actual contaminant degradation rates in the subsurface, which is only accurate for quasi-steady state conditions. However, in reality, subsurface systems are highly dynamic, often resulting in fluctuations of the water table. To quantify the effects of water table fluctuations on NSZD rates, a simulated biodiesel spill in a 400 cm long, 100 cm wide and 150 cm tall sandtank was subjected to lowering and raising the water table, while soil-gas chemistry and surface CO2 and CH4 effluxes were measured. Results show that water table fluctuations have both short-term (perceived) and long-term (actual) effects on NSZD rates, interpreted using surface efflux measurements. When the water table was lowered, surface effluxes immediately increased up to 3 and 344 times higher than baseline for CO2 and CH4 effluxes, respectively, due to the liberation of anaerobically produced gas accumulated below the water table. After this immediate release, the system then reached quasi-steady state conditions 1.4 to 1.6 times higher for CO2 than baseline conditions, attributed to increased aerobic degradation in the broadened and exposed smear zone. When the water table was raised, quasi-steady state CO2 and CH4 effluxes declined to values of 0.9 and 0.4 times baseline effluxes, respectively, implying that contaminant degradation rates were reduced due to submergence of the smear zone. The findings of this study show that the dynamic effects of water table fluctuations and redistribution of the contaminants affect surface effluxes as well as short-term (perceived) and long-term (actual) contaminant degradation rates. Therefore, water table fluctuations need to be considered when quantifying NSZD at contaminated sites using sparse temporal rate measurements to estimate NSZD rates for extended periods of time (e.g., annual rates).

Manganese valence in oxides formed from in situ chemical oxidation of TCE by KMnO4.

Batch and column experiments designed to simulate in situ chemical oxidation (ISCO) in a sand aquifer were conducted to create Mn-oxides (MnOx) by oxidation of trichloroethylene (TCE) with permanganate (MnO4-). Electron energy-loss spectroscopy (EELS) and X-ray photoelectron spectroscopy (XPS) were used to quantify Mn valence in the oxides. The valence of Mn in the MnOx generated in near-source ISCO conditions was 2.2 and 2.3 when formed at low (<3) and neutral (6-7) pH conditions, respectively. There is no significant difference between these values. Valence was found to be sensitive to the preparation method and to aging. When formed in the presence of excess MnO4-, or aged for 3 months, Mn valence ranged from 2.5 to 3.6. Aging in a lower pH environment inhibited Mn oxidation. The EELS and XPS methods provided similar results, but there was a slight bias to higher values for XPS. This work demonstrates that MnO2(s) may not be the main product of MnO4- reaction with chlorinated solvents as is commonly assumed and that the efficiency of ISCO treatment may be greater than previously known.

Scale dependence of effective geochemical rates in weathering mine waste rock.

Hydrogeochemical models for the prediction of drainage quality from full-scale mine waste-rock piles are often parameterized using data from small-scale laboratory or field experiments of short duration. Yet, many model parameters and processes (e.g., sulfide-oxidation rates) vary strongly with the spatiotemporal dimensions of the experiment: the "upscaling" of prediction models remains a critical challenge for mine-waste management worldwide. Here, we investigate scale dependence in laboratory and field experiments that spanned orders-of-magnitude in size (i.e. 2 kg to 100,000 kg) at the Antamina mine in Peru. Normalized drainage mass loading rates systematically decreased with increasing scale, irrespective of waste-rock type. A process-based reactive-transport model was used to simulate observed rates and reproduce the geochemical composition of drainage across scales. Long-term trends in drainage quality could be quantitatively reproduced when the model was parameterized with mostly scale- and experiment-specific measured bulk properties or literature values, leaving geochemical rate coefficients the sole calibrated model parameters. Analysis of these fitted parameters revealed that the scale dependence of geochemical rates was largely explained by reactive mineral surface area. This work demonstrates that practical drainage quality predictions for full-scale waste-rock piles can be established from readily available bulk parameters determined at multiple scales.

Barometric-pumping controls fugitive gas emissions from a vadose zone natural gas release.

Subsurface natural gas release from leaking oil and gas wells is a major environmental concern. Gas migration can cause aquifer contamination, explosive conditions in soil gas, and greenhouse gas emissions. Gas migration is controlled by complex interacting processes, thus constraining the distribution and magnitude of "fugitive gas" emissions remains a challenge. We simulated wellbore leakage in the vadose zone through a controlled release experiment and demonstrate that fugitive gas emissions can be directly influenced by barometric pressure changes. Decreases in barometric-pressure led to surface gas breakthroughs (>20-fold increase in <24 hours), even in the presence of low-permeability surficial soils. Current monitoring strategies do not consider the effect of barometric pressure changes on gas migration and may not provide adequate estimates of fugitive gas emissions. Frequent or continuous monitoring is needed to accurately detect and quantify fugitive gas emissions at oil and gas sites with a deep water table.

Long-term monitoring of waste-rock weathering at the Antamina mine, Peru.

The weathering of mine waste rock can cause release of metal-laden and acidic drainage that requires long-term and costly environmental management. To identify and quantify the geochemical processes and physical transport mechanisms controlling drainage quality, we monitored the weathering of five large-scale (20,000 t) instrumented waste-rock piles of variable and mixed-composition at the Antamina mine, Peru, in a decade-long monitoring program. Fine-grained, sulfidic waste rock with low-carbonate content exhibited high sulfide oxidation rates (>1 g S kg-1 waste rock yr-1) and within 7 years produced acidic (pH < 3) drainage with high Cu and Zn concentrations in the g L-1 range. In contrast, drainage from coarse, carbonate-rich waste rock remained neutral for >10 years and had significantly lower metal loads. Efficient metal retention (>99%) caused by sorption and secondary mineral formation of e.g., gypsum, Fe-(oxy)hydroxides, and Cu/Zn-hydroxysulfates enforced strong (temporary) controls on drainage quality. Furthermore, reactive waste-rock fractions, as small as 10% of total mass, dominated the overall drainage chemistry from the waste-rock piles through internal mixing. This study demonstrates that a reliable prediction of the timing and quality of waste-rock drainage on practice-relevant spatiotemporal scales requires a quantitative understanding of the prevailing in-situ porewater conditions, secondary mineralogy, and spatial distribution of reactive waste-rock fractions in composite piles.

Laboratory-scale experimental and modelling investigations of 222Rn profiles in chemically heterogeneous LNAPL contaminated vadose zones.

The potential of LNAPL delineation by 222Rn soil-gas monitoring in a chemically heterogeneous vadose zone was investigated in this study based on laboratory (batch and columns) experiments and numerical modelling. An enhanced version of the MIN3P reactive transport code was used to simulate Rn transport in both uncontaminated and NAPL-contaminated vadose zones and results were validated against analytical solutions and laboratory experiments. Results show that 222Rn activity profiles are mainly controlled by porous media 222Rn production, vadose zone fluid saturations and especially the type and distribution of NAPL in contaminated areas. The results also show that decreases in 222Rn activity and variations in activity gradients provide evidence for the presence and saturation of NAPL. This study demonstrates that LNAPL delineation via 222Rn gas surveys at contaminated sites works best, if gas measurements extend as deep as possible and include regions where 222Rn activity decreases due to elevated NAPL content. In addition, collection and analysis of depth-discrete gas samples allows the characterization of vertical NAPL distribution based on the 222Rn activity gradient. The determination of 222Rn production in the unsaturated zone, as well as water capillary pressure curves are of key importance in enabling the discrimination of an uncontaminated from a NAPL-contaminated area.

Mobilization of Metal(oid) Oxyanions through Circumneutral Mine Waste-Rock Drainage.

Most studies on the weathering of mine waste rock focus on the generation of acidic drainage with high metal concentrations, whereas metal(loid) release under neutral-rock drainage (NRD) conditions has received limited attention. Here, we present geochemical and mineralogical data from a long-term (>10 years) kinetic testing program with 50 waste-rock field barrels at the polymetallic Antamina mine in Peru. The weathering of most rock lithologies in the field experiments generated circumneutral to alkaline drainage (6 < pH < 9) but with concentrations of the oxyanion-forming metal(loid)s As, Mo, Se, and Sb in the mg/L range. The mobilization of As and Sb was particularly efficient from intrusive, marble and hornfels rocks that contained labile As- and Sb-sulfides, irrespective of bulk elemental content or waste-rock reactivity. High-alkalinity drainage from these materials sustained neutral-pH conditions that are unfavorable to oxyanion adsorption onto Fe-(oxyhydr)oxides and, therefore, enhanced As and Sb leaching. The release of Mo and Se from sulfidic skarn and intrusive waste rock was more proportional to elemental content but equally enhanced by pH-inhibited adsorption and negligible secondary mineral precipitation under NRD conditions. Our results demonstrate that oxyanion concentrations of environmental concern may be conveyed by neutral- to alkaline-pH waste-rock drainage and should be a focus of mine wastewater monitoring programs.

Hydro-biogeochemical impacts of fugitive methane on a shallow unconfined aquifer.

Oil and gas development can result in natural gas migration into shallow groundwater. Methane (CH4), the primary component of natural gas, can subsequently react with solutes and minerals in the aquifer to create byproducts that affect groundwater chemistry. Hydro-biogeochemical processes induced by fugitive gas from leaky oil and gas wells are currently not well understood. We monitored the hydro-biogeochemical responses of a controlled natural gas release into a well-studied Pleistocene beach sand aquifer (Canadian Forces Base Borden, Ontario, Canada). Groundwater samples were collected before, during, and up to 700 days after gas injection and analyzed for pH, major and minor ions, alkalinity, dissolved gases, stable carbon isotope ratios of CO2 and CH4, and microbial community composition. Gas injection resulted in a dispersed plume of free and dissolved phase natural gas, affecting groundwater chemistry in two distinct temporal phases. Initially (i.e. during and immediately after gas injection), pH declined and major ions and trace elements fluctuated; at times increasing above baseline concentrations. Changes in the short-term were due to invasion of deep groundwater with elevated total dissolved solids entrained with the upward migration of free phase gas and, reactions that were instigated through the introduction of constituents other than CH4 present in the injected gas (e.g. CO2). At later times, more pronounced aerobic and anaerobic CH4 oxidation led to subtle increases in major ions (e.g. Ca2+, H4SiO4) and trace elements (e.g. As, Cr). Microbial community profiling indicated a persistent perturbation to community composition with a conspicuous ingrowth of taxa implicated in aerobic CH4 oxidation as well anaerobic S, N and Fe species metabolism.

Modeling vadose zone processes during land application of food-processing waste water in California's Central Valley.

Land application of food-processing waste water occurs throughout California's Central Valley and may be degrading local ground water quality, primarily by increasing salinity and nitrogen levels. Natural attenuation is considered a treatment strategy for the waste, which often contains elevated levels of easily degradable organic carbon. Several key biogeochemical processes in the vadose zone alter the characteristics of the waste water before it reaches the ground water table, including microbial degradation, crop nutrient uptake, mineral precipitation, and ion exchange. This study used a process-based, multi-component reactive flow and transport model (MIN3P) to numerically simulate waste water migration in the vadose zone and to estimate its attenuation capacity. To address the high variability in site conditions and waste-stream characteristics, four food-processing industries were coupled with three site scenarios to simulate a range of land application outcomes. The simulations estimated that typically between 30 and 150% of the salt loading to the land surface reaches the ground water, resulting in dissolved solids concentrations up to sixteen times larger than the 500 mg L(-1) water quality objective. Site conditions, namely the ratio of hydraulic conductivity to the application rate, strongly influenced the amount of nitrate reaching the ground water, which ranged from zero to nine times the total loading applied. Rock-water interaction and nitrification explain salt and nitrate concentrations that exceed the levels present in the waste water. While source control remains the only method to prevent ground water degradation from saline wastes, proper site selection and waste application methods can reduce the risk of ground water degradation from nitrogen compounds.

High resolution spatial and temporal evolution of dissolved gases in groundwater during a controlled natural gas release experiment.

Fugitive gas comprised primarily of methane (CH4) with traces of ethane and propane (collectively termed C1-3) may negatively impact shallow groundwater when unintentionally released from oil and natural gas wells. Currently, knowledge of fugitive gas migration, subsurface source identification and oxidation potential in groundwater is limited. To advance understanding, a controlled release experiment was performed at the Borden Research Aquifer, Canada, whereby 51m3 of natural gas was injected into an unconfined sand aquifer over 72days with dissolved gases monitored over 323days. During active gas injection, a dispersed plume of dissolved C1-3 evolved in a depth discrete and spatially complex manner. Evolution of the dissolved gas plume was driven by free-phase gas migration controlled by small-scale sediment layering and anisotropy. Upon cessation of gas injection, C1-3 concentrations increased to the greatest levels observed, particularly at 2 and 6m depths, reaching up to 31.5, 1.5 and 0.1mg/L respectively before stabilizing and persisting. At no time did groundwater become fully saturated with natural gas at the scale of sampling undertaken. Throughout the experiment the isotopic composition of injected methane (δ13C of -42.2‰) and the wetness parameter (i.e. the ratio of C1 to C2+) constituted excellent tracers for the presence of fugitive gas at concentrations >2mg/L. At discrete times C1-3 concentrations varied by up to 4 orders of magnitude over 8m of aquifer thickness (e.g. from <0.01 to 30mg/L for CH4), while some groundwater samples lacked evidence of fugitive gas, despite being within 10m of the injection zone. Meanwhile, carbon isotope ratios of dissolved CH4 showed no evidence of oxidation. Our results show that while impacts to aquifers from a fugitive gas event are readily detectable at discrete depths, they are spatially and temporally variable and dissolved methane has propensity to persist.

Microbial and geochemical controls on waste rock weathering and drainage quality.

Bacteria can adversely affect the quality of drainage released from mine waste by catalyzing the oxidation of sulfide minerals and thereby accelerating the release of acidity and metals. However, the microbiological and geochemical controls on drainage quality from unsaturated and geochemically heterogeneous waste rock remain poorly understood. Here, we identified coexisting neutrophilic and acidophilic bacteria in different types of waste rock, indicating that robust endemic consortia are sustained within pore-scale microenvironments. Subsequently, natural weathering was simulated in laboratory column experiments with waste rock that contained either in-situ microbial consortia or suppressed populations with up to 1000 times smaller abundance and reduced phenotypic diversity after heating and drying. Drainage from waste rock with in-situ populations was up to two pH units lower and contained up to 16 times more sulfate and heavy metals compared to drainage from waste rock bearing treated populations, indicating significantly higher sulfide-oxidation rates. The drainage chemistry was further affected by sorption and formation of secondary-mineral phases (e.g., gypsum and hydroxy-carbonates). This study provides direct evidence for the existence of diverse microbial communities in waste rock and their important catalytic role on weathering rates, and illustrates the mutual controls of microbiology and geochemistry on waste-rock drainage quality.

Unintentional contaminant transfer from groundwater to the vadose zone during source zone remediation of volatile organic compounds.

Historical heavy use of chlorinated solvents in conjunction with improper disposal practices and accidental releases has resulted in widespread contamination of soils and groundwater in North America and worldwide. As a result, remediation of chlorinated solvents is required at many sites. For source zone treatment, common remediation strategies include in-situ chemical oxidation (ISCO) using potassium or sodium permanganate, and the enhancement of biodegradation by primary substrate addition. It is well known that these remediation methods tend to generate gas (carbon dioxide (CO2) in the case of ISCO using permanganate, CO2 and methane (CH4) in the case of bioremediation). Vigorous gas generation in the presence of chlorinated solvents, which are categorized as volatile organic contaminants (VOCs), may cause gas exsolution, ebullition and stripping of the contaminants from the treatment zone. This process may lead to unintentional 'compartment transfer', whereby VOCs are transported away from the contaminated zone into overlying clean sediments and into the vadose zone. To this extent, benchtop column experiments were conducted to quantify the effect of gas generation during remediation of the common chlorinated solvent trichloroethylene (TCE/C2Cl3H). Both ISCO and enhanced bioremediation were considered as treatment methods. Results show that gas exsolution and ebullition occurs for both remediation technologies. Facilitated by ebullition, TCE was transported from the source zone into overlying clean groundwater and was subsequently released into the column headspace. For the case of enhanced bioremediation, the intermediate degradation product vinyl chloride (VC) was also stripped from the treatment zone. The concentrations measured in the headspace of the columns (TCE ∼300ppm in the ISCO column, TCE ∼500ppm and VC ∼1380ppm in the bioremediation column) indicate that substantial transfer of VOCs to the vadose zone is possible. These findings provide direct evidence for the unintended spreading of contaminants as a result of remediation efforts, which can, under some circumstances, result in enhanced risks for soil vapour intrusion.