A Critical Review of State-of-the-Art and Emerging Approaches to Identify Fracking-Derived Gases and Associated Contaminants in Aquifers.

High-volume, hydraulic fracturing (HVHF) is widely applied for natural gas and oil production from shales, coals, or tight sandstone formations in the United States, Canada, and Australia, and is being widely considered by other countries with similar unconventional energy resources. Secure retention of fluids (natural gas, saline formation waters, oil, HVHF fluids) during and after well stimulation is important to prevent unintended environmental contamination, and release of greenhouse gases to the atmosphere. Here, we critically review state-of-the-art techniques and promising new approaches for identifying oil and gas production from unconventional reservoirs to resolve whether they are the source of fugitive methane and associated contaminants into shallow aquifers. We highlight future research needs and propose a phased program, from generic baseline to highly specific analyses, to inform HVHF and unconventional oil and gas production and impact assessment studies. These approaches may also be applied to broader subsurface exploration and development issues (e.g., groundwater resources), or new frontiers of low-carbon energy alternatives (e.g., subsurface H2 storage, nuclear waste isolation, geologic CO2 sequestration).

Estimates of water and solute release from a coal waste rock dump in the Elk Valley, British Columbia, Canada.

Long term (1999 to 2014) flow and water quality data from a rock drain located at the base of a coal waste rock dump constructed in the Elk Valley, British Columbia was used to characterize the release of three solutes (NO3-, Cl- and SO42-) from the dump and obtain whole dump estimates of net percolation (NP). The concentrations of dump derived solutes in the rock drain water were diluted by snowmelt waters from the adjacent natural watershed during the spring freshet and reached a maximum concentration during the winter baseflow period. Historical peak baseflow concentrations of conservative ions (NO3- and Cl-) increased until 2006/07 after which they decreased. This decrease was attributed to completion of the flushing of the first pore volume of water stored within the dump. The baseflow SO42- concentrations increased proportionally with NO3- and Cl- to 2007, but then continued to slowly increase as NO3- and Cl- concentrations decreased. This was attributed to ongoing production of SO42- due to oxidation of sulfide minerals within the dump. Based on partitioning of the annual volume of water discharged from the rock drain to waste rock effluent (NP) and water entering the rock drain laterally from the natural watershed, the mean NP values were estimated to be 446±50mm/a (area normalized net percolation/year) for the dump and 172±71mm/a for the natural watershed. The difference was attributed to greater rates of recharge in the dump from summer precipitation compared to the natural watershed where rainfall interception and enhanced evapotranspiration will increase water losses. These estimates included water moving through subsurface pathways. However, given the limitations in quantifying these flows the estimated NP rates for both the natural watershed and the waste rock dump are considered to be low, and could be much higher (e.g. ~450mm/a and ~800mm/a).

Mineralogical controls on aluminum and magnesium in uranium mill tailings: Key Lake, Saskatchewan, Canada.

The mineralogy and evolution of Al and Mg in U mill tailings are poorly understood. Elemental analyses (ICP-MS) of both solid and aqueous phases show that precipitation of large masses of secondary Al and Mg mineral phases occurs throughout the raffinate neutralization process (pH 1-11) at the Key Lake U mill, Saskatchewan, Canada. Data from a suite of analytical methods (ICP-MS, EMPA, laboratory- and synchrotron-based XRD, ATR-IR, Raman, TEM, EDX, ED) and equilibrium thermodynamic modeling showed that nanoparticle-sized, spongy, porous, Mg-Al hydrotalcite is the dominant mineralogical control on Al and Mg in the neutralized raffinate (pH ≥ 6.7). The presence of this secondary Mg-Al hydrotalcite in mineral samples of both fresh and 15-year-old tailings indicates that the Mg-Al hydrotalcite is geochemically stable, even after >16 years in the oxic tailings body. Data shows an association between the Mg-Al hydrotalcite and both As and Ni and point to this Mg-Al hydrotalcite exerting a mineralogical control on the solubility of these contaminants.

Controls and rates of acid production in commercial-scale sulfur blocks.

Acidic drainage (pH 0.4-1.0) from oxidizing elemental sulfur (S(0)) blocks is an environmental concern in regions where S(0) is stockpiled. In this study, the locations, controls, and rates of H(2)SO(4) production in commercial-scale S(0) blocks ( approximately 1-2 x 10(6) m(3)) in northern Alberta, Canada, were estimated. In situ modeling of O(2) concentrations ([O(2)]) suggest that 70 to >97% of the annual H(2)SO(4) production occurs in the upper 1 m of the blocks where temperatures increase to >15 degrees C during the summer. Laboratory experiments show that S(0) oxidation rates are sensitive to temperature (Q(10) = 4.3) and dependent on the activity of autotrophic S(0)-oxidizing microbes. The annual efflux of SO(4) in drainage water from a S(0) block (5.5 x 10(5) kg) was within the estimated range of SO(4) production within the block (2.7 x 10(5) to 1.2 x 10(6) kg), suggesting that H(2)SO(4) production and removal rates were approximately equal during the study period. The low mean relative humidity within the block (68%; SD = 17%; n = 21) was attributed to osmotic suction from elevated H(2)SO(4) concentrations and suggests a mean in situ pH of approximately -2.1. The low pH of drainage waters was attributed to the mixing of fresh infiltrating water and low-pH in situ water. Heat generation during S(0) oxidation was an important factor in maintaining elevated temperatures (mean, 11.1 degrees C) within the block. The implications of this research are relevant globally because construction methods and the physical properties of S(0) blocks are similar worldwide.

Microbial distribution and diversity in saturated, high pH, uranium mine tailings, Saskatchewan, Canada.

Microbiological analyses were conducted on core samples collected along a vertical profile (0-66 m below surface) from the tailings management facility (TMF) at the Rabbit Lake uranium mine in northern Saskatchewan, Canada. Bacterial numbers in the core materials were similar to surrounding soils and surface waters, regardless of the seemingly unfavorable pH (mean=9.9) and temperature (approximately 0 degrees C) in the TMF. The greatest number of viable cells (105 CFU/g) was detected at the interface between the tailings and overlying standing water, below which cell counts decreased rapidly with depth. Whole-community metabolic profiles for samples from the different depths grouped into 3 clusters; however, these groups could not be positively correlated with sampling depth, temperature, redox potential, pH, or ore-mill feed. Flow-cell studies demonstrated microbial communities in the tailings surface water could develop biofilms and maintain cell activity at both pH 10 and 7, and altering the pH between these 2 values had little effect on biofilm viability. These results demonstrate the resilience and adaptive nature of naturally occurring microbial communities and signify a potential role of microbial activity in the long-term geochemical evolution of the TMF.

Long-term reactive solute transport in an aquitard using a centrifuge model.

Characterizing and predicting reactive solute transport in low hydraulic conductivity (K) clay-rich media is challenging because the very long transport time for solutes renders conventional column tests impractical. In this study, a centrifugation technique was developed to assess the transport of a simple aqueous solution (NaCl) by accelerating flow by centrifugal force through low K (1.1 x 10(-11) m/s) core samples. Duplicate cores (52-mm length x 33-mm diameter) were centrifuged at 330 xg for 90 d to model the migration of saline pore water (0.5 M NaCl) under in situ conditions through an approximately 17-m-thick clay prototype over approximately 24,000 years. A PHREEQC one-dimensional reactive solute transport code simulated effluent breakthrough of the NaCl during centrifugation, with best-fit cation exchange coefficients similar to batch tests. The calibrated code was used to predict solute profile development over the long term in the prototype or simulated field-scale conditions. Chromatographic separation of solutes due to ion exchange was evident over several meters in the simulated prototype and the field profile. The applicability of centrifugation methods to predict transport of more complex suites of reactive solutes over the long term is yet to be verified.

Microbial response to environmental gradients in a ceramic-based diffusion system.

A solid, porous matrix was used to establish steady-state concentration profiles upon which microbial responses to concentration gradients of nutrients or antimicrobial agents could be quantified. This technique relies on the development of spatially defined concentration gradients across a ceramic plate resulting from the diffusion of solutes through the porous ceramic matrix. A two-dimensional, finite-element numerical transport model was used to predict the establishment of concentration profiles, after which concentration profiles of conservative tracers were quantified fluorometrically and chemically at the solid-liquid interface to verify the simulated profiles. Microbial growth responses to nutrient, hypochloride, and antimicrobial concentration gradients were then quantified using epifluorescent or scanning confocal laser microscopy. The observed microbial response verified the establishment and maintenance of stable concentration gradients along the solid-liquid interface. These results indicate the ceramic diffusion system has potential for the isolation of heterogeneous microbial communities as well as for testing the efficacy of antimicrobial agents. In addition, the durability of the solid matrix allowed long-term investigations, making this approach preferable to conventional gel-stabilized systems that are impeded by erosion as well as expansion or shrinkage of the gel.

High resolution pore water delta2H and delta18O measurements by H2O(liquid)-H2O(vapor) equilibration laser spectroscopy.

A new H2O(liquid)-H2O(vapor) pore water equilibration and laser spectroscopy method provides a fast way to obtain accurate high resolution deltaD and delta18O profiles from single core samples from saturated and unsaturated geologic media. The precision and accuracy of the H2O(liquid)-H2O(vapor) equilibration method was comparable to or better than conventional IRMS-based methods, and it can be conducted on geologic cores that contain volumetric water contents as low as 5%. Significant advantages of the H2O(liquid)-H2O(vapor) pore water equilibration method and laser isotopic analysis method include dual hydrogen- and oxygen-isotope assays on single small core samples, low consumable and instrumentation costs, and the potential for field-based hydrogeologic profiling. A single core is sufficient to obtain detailed vertical isotopic depth profiles in geologic, soil, and lacustrine pore water, dramatically reducing the cost of obtaining pore water by conventional wells or physical water extraction methods. In addition, other inherent problems like contamination of wells by leakage and drilling fluids can be eliminated.

High-precision laser spectroscopy D/H and 18O/16O measurements of microliter natural water samples.

Newly available gas analyzers based on off-axis integrated cavity output spectroscopy (OA-ICOS) lasers have been advocated as an alternative to conventional isotope-ratio mass spectrometers (IRMS) for the stable isotopic analysis of water samples. In the case of H2O, OA-ICOS is attractive because it has comparatively low capital and maintenance costs, the instrument is small and field laboratory portable, and provides simultaneous D/H and 16O/18O ratio measurements directly on H2O molecules with no conversion of H2O to H2, CO, or H2/CO2-water equilibration required. Here we present a detailed assessment of the performance of a liquid-water isotope analyzer, including instrument precision, estimates of sample memory and sample mass effects, and instrumental drift. We provide a recommended analysis procedure to achieve optimum results using OA-ICOS. Our results show that, by using a systematic sample analysis and data normalization procedure routine, measurement accuracies of +/-0.8 per thousand for deltaD and +/-0.1 per thousand delta18O are achievable on nanoliter water samples. This is equivalent or better than current IRMS-based methods and at a comparable sample throughput rate.

Dynamics and stable isotope composition of gaseous and dissolved oxygen.

The vadose zone and ground water environments are a sink for atmospheric O(2). The pathways and rates of O(2) consumption are primarily related to the availability and rate of oxidation of key reductants (e.g., organics, sulfides), through a combination of biological or abiotic reactions. The range in delta(18)O of O(2) in the subsurface is large, from +20 per thousand to +39 per thousand (Vienna Standard Mean Ocean Water) in the vadose zone and from +12 per thousand to +46 per thousand in ground water. The aggregated O(2) isotope fractionation by consumption (alpha(k)) was found to range from 0.970 to 1.300 and 0.980 to 1.030 in vadose zones and aquifers, respectively. These data suggest the delta(18)O patterns in both unsaturated zones and aquifers can be attributed to microbially mediated reactions (alpha(k)= range from 0.975 to 1.000), but there are apparently other inverse isotope fractionating processes (alpha(k) > 1.000). Circumstantial evidence suggested O(2) processed during the sulfide oxidation and precipitation of Fe-oxyhydroxides process (or other unidentified processes) could be the cause of the significant (18)O depletions. Overall, delta(18)O data from vadose zones and ground water revealed that isotope fractionation by consumption of gaseous and dissolved O(2) in the subsurface and ground water environments is more complicated than what has classically been attributed solely to geomicrobial respiration. Given the questions and inexplicable data arising from this study, further detailed research on O(2) consuming processes in the Earth's subsurface and ground water is warranted.

Characterization of aqueous lead removal by phosphatic clay: equilibrium and kinetic studies.

Immobilization of heavy metals from contaminated environments is an emerging field of interest from both resource conservation and environmental remediation points of view. This study investigated the feasibility of using phosphatic clay, a waste by-product of the phosphate mining industry, as an effective sorbent for Pb from aqueous effluents. The major parameters controlling aqueous Pb removal, viz. initial metal ion concentrations, solution pH, sorbent amounts, ionic strength and presence of both inorganic and organic ligands were evaluated using batch experiments. Results demonstrated that aqueous Pb removal efficiency of phosphatic clay is controlled mainly by dissolution of phosphatic clay associated fluoroapatite [Ca(10)(PO(4))(5)CaCO(3)(F,Cl,OH)(2)], followed by subsequent precipitation of geochemically stable pyromorphite [Pb(10)(PO(4))(6)(F,Cl,OH)(2)], which was confirmed by both X-ray diffraction (XRD) and scanning electron microscopic (SEM) analysis. Lead removal efficiency of phosphatic clay increased with increasing pH, sorbent amount and decreasing ionic strength. It also depends on the nature of complexing ligands. Formation of insoluble calcium oxalate and lead oxalate in the presence of oxalic acid explained high uptake of Pb by phosphatic clay from aqueous solution. However, Pb sorption kinetics onto phosphatic clay were biphasic, with initially fast reactions followed by slow and continuous Pb removal reactions. The slow reactions may include surface sorption, co-precipitation and diffusion. The exceptional capability of phosphatic clay to remove aqueous Pb demonstrated its potential as a cost effective way to remediate Pb-contaminated water, soils and sediments.

An assessment of a mesocosm approach to the study of microbial respiration in a sandy unsaturated zone.

Microbial respiration rates were determined through a 3.2 m thick, sandy unsaturated zone in a 2.4 m diameter x 4.6 m high mesocosm. The mesocosm was maintained under near constant temperature (18 degrees to 23 degrees C) and reached steady moisture content conditions after several hundred days. Soil-gas CO2 concentrations in the mesocosm ranged from 0.09% to 3.31% and increased with depth. Respiration rates within the mesocosm were quantified over a 342-day period using measured CO2 concentrations and a transient, one-dimensional finite-element model. Microbial respiration rates were 2 x 10(-1) micrograms C.g-1.d-1 throughout most of the system, but decreased to 10(-4) to 10(-3) micrograms C.g-1.d-1 within the capillary fringe. Microbial respiration rates were also determined in minicosms (500 g sample mass) over a range in temperatures (4 degrees to 30 degrees C) and volumetric moisture contents (0.044 to 0.37). The functional dependence of CO2 production on temperature and soil-moisture content was similar for the two scales of laboratory observation. Respiration rates in the minicosms, for temperatures and moisture contents in the mesocosm, were up to an order of magnitude greater than those determined for the mesocosm. The higher respiration rates in the minicosms, compared to the mesocosm, were attributed to greater disturbance of the samples and to shorter acclimation time in the minicosms. Extrapolating the laboratory respiration rates to field conditions yielded rates that were two to three orders of magnitude greater than rates previously determined in situ for C-horizon material. Results show that in situ microbial reaction rates determined using disturbed samples in minicosms and mesocosms yielded respiration rates that greatly exceeded field conditions. Mesocosms can, however, provide a useful environment for conducting process-related research in unsaturated environments.

Speciation of dissolved iron(III) and iron(II) in water by on-line coupling of flow injection separation and preconcentration with inductively coupled plasma mass spectrometry.

A method has been developed for the speciation of trace dissolved Fe(II) and Fe(II) in water by on-line coupling of flow injection separation and preconcentration with inductively coupled plasma mass spectrometry (ICPMS). Selective determination of Fe(III) in the presence of Fe(II) was made possible by on-line formation and sorption of the Fe(III)-pyrrolidinecarbodithioate (PDC) complex in a PTFE knotted reactor over a sample acidity range of 0.07-0.4 mol L(-1) HCl, elution with 1 mol L(-1) HNO3, and detection by ICPMS. Over a sample acidity range of 0.001-0.004 mol L(-1) HCl, the sum of Fe(III) and Fe(II), i.e., Fe(III + II), could be determined without the need for preoxidation of Fe(II) to Fe(III). The concentration of Fe(II) was obtained as the difference between those of Fe(III + II) and Fe(III). With a sample flow rate of 5 mL min(-1) and a 30-s preconcentration time, an enhancement factor of 12, a retention efficiency of 80%, and a detection limit (3s) of 0.08 microg L(-1) were obtained at a sampling frequency of 21 samples h(-1). The relative standard deviation (n = 11) was 2.9% at the 10 microg L(-1) Fe(III) level. Recoveries of spiked Fe(III) and Fe(II) in local tap water, river water, and groundwater samples ranged from 95% to 103%. The concentrations of Fe(III) and Fe(II) in synthetic aqueous mixtures obtained by the proposed method were in good agreement with the spiked values. The result for total iron concentration in the river water reference material SLRS-3 was in good agreement with the certified value. The method was successfully applied to the determination of trace dissolved Fe(III) and Fe(II) in local tap water, river water, and groundwater samples.

An automated technique for measuring deltaD and delta18O values of porewater by direct CO2 and H2 equilibration.

The stable-oxygen and -hydrogen isotopic values (deltaD, delta18O) of porewater in geologic media are commonly determined on water obtained by extraction techniques such as centrifugation, mechanical squeezing, vacuum heating and cryogenic microdistillation, and azeotropic distillation. Each of these techniques may cause isotopic fractionation as part the extraction process and each is laborious. Here we demonstrate a new approach to obtain automated, high-precision deltaD and delta18O measurements of porewater in geologic sediments by direct H2- and CO2-porewater equilibration using a modified commercial CO2-water equilibrator. This technique provides an important and cost-effective improvement over current extraction methods, because many samples can be rapidly analyzed with minimal handling, thereby reducing errors and potential for isotopic fractionation. The precision and accuracy of direct H2- and CO2-porewater equilibration is comparable to or better than current porewater extraction methods. Finally, the direct equilibration technique allows investigators to obtain high-resolution (cm scale) porewater deltaD and delta18O profiles using cores from individual boreholes, eliminating the need for costly piezometers or conventional porewater extractions.

Distribution and Biogeochemical Importance of Bacterial Populations in a Thick Clay-Rich Aquitard System.

To investigate the distribution of microbial biomass, populations and activities within a clay-rich, low hydraulic conductivity (10-11 to 10-12 m s-1) aquitard complex, cores were aseptically obtained from a series of overlying clayey deposits; a fractured till, unfractured till (20-30 ka BP), a disturbed interfacial zone (20-30 ka BP), and a Cretaceous clay aquitard (71-72 Ma BP). The results of confocal microscopy studies, culture methods, molecular approaches, and extractive fatty acid analyses all indicated low bacterial numbers that were non-homogeneously distributed within the sediments. Various primers for catabolic genes were used to amplify extracted DNA. Results indicated the presence of eubacterial 23S rDNA, and the narH gene for nitrate reductase and ribulose-1,5-biphosphate carboxylase (RuBP carboxylase). Although there was no evidence of limitation by electron acceptors or donors, sulfate-reducing bacteria were not detected below the fractured till zone, using PCR, enrichment, or culture techniques. Denaturing gradient gel electrophoresis (DGGE) analyses indicated differences in community composition and abundance between the various geologic units. Results of FAME analyses of sediments yielded detectable extractable fatty acids throughout the aquitard complex. Bacterial activities were demonstrated by measuring mineralization of (14C) glucose. Porewater chemistry and stable isotope data were in keeping with an environment in which extremely slow growing, low populations of bacteria exert little impact. The observations also support the contention that in low permeability sediments bacteria may survive for geologic time periods.

Determination of (ultra)trace amounts of arsenic(III) and arsenic(V) in water by inductively coupled plasma mass spectrometry coupled with flow injection on-line sorption preconcentration and separation in a knotted reactor.

A method has been developed for determination of (ultra)trace amounts of As(III) and As(V) in water by flow injection on-line sorption preconcentration and separation coupled with inductively coupled plasma mass spectrometry (ICPMS) using a knotted reactor (KR). The determination of As(III) was achieved by selective formation of the As(III)-pyrrolidine dithiocarbamate complex over a sample acidity range of 0.01-0.7 mol L-1 HNO3, its adsorption onto the inner walls of the KR made from 150-cm-long, 0.5-mm-i.d. PTFE tubing, elution with 1 mol L-1 HNO3, and detection by ICPMS. Total inorganic arsenic was determined after prereduction of As(V) to As(III) in a 1% (m/v) L-cysteine-0.03 mol L-1 HNO3 media. The concentration of As(V) was calculated by difference (the total inorganic arsenic and As(III)). Owing to the group-specific character of the chelating agent, and the use of an efficient rinsing step before elution, the interferences encountered in conventional ICPMS from common major matrix, alkali and alkaline earth metals, and chlorides were eliminated. The presence of organoarsenic species such as monomethylarsonate and dimethylarsinate in water samples had no effect on the results of As(III) and As(V). Thus, the method can be applied to the speciation analysis of inorganic arsenic at submicrogram per liter levels in aqueous solutions with high total content of dissolved solid and/or high content of chlorides. Using a preconcentration time of 60 s and a sample flow rate of 5 mL min-1, an enhancement factor of 22 was achieved in comparison with conventional ICPMS. The time required for a single determination was 200 s. The detection limits (3s) was evaluated to be 0.021 microgram L-1 for As(III) and 0.029 microgram L-1 for total inorganic arsenic. The precision for 14 replicate determinations of 1 microgram L-1 As(III) was 2.8% (RSD) with drift correction and 3.9% (RSD) without drift correction. The concentrations of As(III) and As(V) in synthetic mixtures obtained by the present method were in good agreement with expected values. Results obtained by the proposed method for total arsenic in a river water reference material agreed well with certified and recently reevaluated values. The method was also applied to the speciation analysis of inorganic arsenic in porewaters.

Development of steady-state diffusion gradients for the cultivation of degradative microbial consortia.

A diffusion gradient plate was constructed and evaluated for its potential use in the isolation of degradative microbial consortia from natural habitats. In this model, a steady-state concentration gradient of diclofop methyl, established by diffusion through an agarose gel, provided the carbon for microbial growth. Colonization of the gel surface was observed with epifluorescence and scanning confocal laser microscopy to determine microbial responses to the diclofop gradient. A detectable gradient developed over a narrow band (<10 mm). Consequently, quantitative analyses of the microbial response to the gradient were difficult to obtain. A two-dimensional, finite-element numerical transport model for advective-diffusive transport was used to simulate concentration and flux profiles in the physical model. The simulated profiles were correlated with the measured concentration gradient (R = 0.89) and the cell numbers on the gel surface (R = 0.85). The numerical model was subsequently used to redesign the physical model. The detectable concentration gradient in the modified physical model extended over the length of the gel (38 mm). The simulated profile again showed a good correlation with the measured profile (R = 0.96) and the microbial responses to the concentration gradient (R = 0.99). It was concluded that these gradients provide the steady-state environments needed to sustain steady-state consortia. They also provide a physical pathway for the development of degradative biofilms from low to high concentrations of toxicants and simulate conditions under which low concentrations of toxicant are supplied at a constant flux over long periods of time, such as the conditions that could occur in natural environments.