Studies of microbial sulfate reduction (MSR) within the hot biosphere are dominated by investigations of deep-sea hydrothermal vents and other modern marine ecosystems1,2,5. This is despite abundant drilling into deep, hot environments in continental settings related to industrial activities including hydrocarbon production and increasingly geological CO2 sequestration, hydrogen storage, and disposal of waste fluids. Understanding the potential geochemical and biological reactions related to the injection of surface-derived fluids into extreme environments (e.g., > 80°C6) assumed to be biologically inactive is critical.
Within modern marine environments the estimates for the upper temperature limit for MSR has increased from 110°C as measured at deep-sea hydrothermal vents1, potentially up to 120°C based on laboratory incubation experiments2,7,8. In these hot marine environments, production of readily degradable organic substrates (e.g., acetate) permits the development of highly active populations of sulfate reducing microorganisms with rapid doubling times of hours as opposed to weeks for other sulfate reducing microorganisms2. To date, reports of biological activity at such temperatures are almost exclusively associated with marine environments1,2,5. In comparison, biological activity in the deep biosphere, hosted in sedimentary basins, is thought to cease at temperatures in excess of 80°C based on observations of hydrocarbon biodegradation in petroleum systems6.
In the deep biosphere, drilling into continental environments where temperatures exceed 80°C provide exceptional opportunities to study biogeochemical processes in non-hydrothermal extreme environments. Advances in directional drilling and enhanced hydrocarbon recovery technology has led to exploitation of hydrocarbons from low-permeability rocks such as shales and siltstones via hydraulic fracturing (HF)9. This process enables sampling of the otherwise inaccessible terrestrial deep biosphere via continuously produced fluids permitting time-resolved sampling of biogeochemical reactions in the deep biosphere.
HF of a hydrocarbon reservoir involves the high-pressure injection of water and sand to create fracture networks thereby forming fluid exchange pathways while generating space for active microbial communities. Along with water and sand, different chemicals (Table SI-1) are co-injected during HF to promote efficient recovery of hydrocarbons (e.g., oxidizing breakers such as persulfates10), and to stabilize the fracture networks (e.g., sand, guar gum). In general, great efforts are made to limit biological reactions as previous research has shown that microbial metabolisms and growth negatively impacts production of hydrocarbons due to production of H2S11 which is toxic and corrosive even at low concentrations. Injected fluids are treated with biocides, yet some microbial groups including resilient endospore-forming lineages12 are consistently reported in produced fluids13. Despite originating from entirely different geological formations and generally at cooler temperatures, microbial communities in HF produced fluids are remarkably consistent11,13,14, leading to the interpretation that microorganisms are introduced via surface-derived populations from nearby water sources followed by environmental selection in the extreme subsurface environment during and subsequent to HF14.
Montney shale formation
Here, we use aqueous and isotopic geochemistry combined with microbial DNA sequencing of HF fluids to evaluate biogeochemical responses to HF within a hot (> 90°C) and deep (> 3 km) environment. We analyzed produced fluids from a hydraulically fractured horizontal well at 3,411 meters depth in Canada’s Montney Formation, an Early Triassic siltstone formation15(Fig. 1). Injected water was sourced from local rivers characterized by low total dissolved solids (TDS < 1,000 mg/L) that mixes upon injection with high salinity (TDS > 170,000 mg/L) formation water in the target reservoir. The measured temperature within the reservoir is 91.6°C and is likely a conservatively low value given that cooling is often associated with the drilling process. Other local Montney reservoirs temperature is typically ≥ 95°C (Fig. SI-1). Injection of HF fluids initially reduces the temperature of fracture-filling fluids as cooler surface water (12°C) interacts with the hot reservoir (91.6°C). Following injection, the well was “shut-in” for 30 days during which time fluid temperature increases, reaching equilibrium on the order of weeks16 with simulations indicating < 2°C difference between geothermal and fluid temperatures17 during the production phase.
Water samples were collected at high density sampling over the first week of production (samples 20 to 900 m3), and then monthly to capture the initial conditions immediately following the onset of production, and monthly changes over the course of one year of production. This enabled the investigation of initial evolution of produced water geochemistry and microbiology related to reactions of highly reactive oxygen-bearing compounds within a reducing environment, the degradation of HF chemicals, and the introduction of microbial populations within the injected fluids. Monthly sampling revealed biological activation of the subsurface sulfur cycle as redox conditions shifted back to a highly reducing environment over time. Previous work and industry experience indicates that the geochemistry of produced water reflects the composition of the original formation water following a year of production18,19.
Oxygenating an anoxic environment
Introducing surface-derived fluids into the subsurface exposes a hypersaline anoxic environment to low salinity, oxygen-rich fluids. Analysis of major anions and cations revealed that many species behave conservatively (e.g., Na, Cl, Li etc.; Fig. SI-2) displaying binary mixing of HF fluids and formation waters as noted by others19. However, non-conservative species including iron (Fe), sulfate, and sulfide reveal water-rock interactions and biological metabolism (i.e., MSR) that are consistent with the initiation of a biologically active subsurface sulfur cycle. Sulfate concentrations in produced water quickly increased to ~ 500 mg/L (Fig. 2a) and remained relatively constant throughout the sampling period (Fig. 3a). Core analysis of reservoir rocks from this region shows evidence of abundant anhydrite (CaSO4) within pore-space (Fig. SI-3) which provides the sulfate source that is buffering dissolved sulfate concentrations in produced fluids.
Injection of oxygen-bearing compounds like persulfates and dissolved O2 results in the oxidation of iron sulfides (i.e., pyrite occurs at up to 3.7 weight %) in the fractured host rock releasing Fe2+ ions. This is shown by the sharp increase in dissolved Fe in initial samples (Fig. 2a). Rapid iron sulfide oxidation through either microbial (Fig. 2) or abiotic reactions produces sulfate causing a corresponding shift in δ18OSO4 from + 12.8 to + 15.3‰, reflecting dissolved oxygen being close to the atmospheric oxygen isotope ratio (δ18Oatm = + 24.2‰) 20. Rapid consumption of reactive oxygen corresponded with δ18OSO4 decreasing to baseline values (Fig. 2a). Previous studies suggest that removal of oxygen via sulfide oxidation occurs rapidly at elevated temperatures consistent with our observations21. The trajectory of the slope of δ18OSO4 plotted vs. δ34SSO4 of initial samples (i.e., 20–100 m3 Fig. 2a) (Fig. 4) is analogous to those in other iron-rich environments where high rates of sulfate recycling22 result in considerable changes in δ18OSO4 coupled with minimal δ34SSO4 variability. Any H2S produced during this period is rapidly re-oxidized or sequestered via reactions with Fe2+ which is consistent with the presence of sulfate reducing and sulfide oxidizing microbial populations in the initial water samples (Fig. 2b).
Microbial life in a mixed redox environment
Initial microbial populations in produced water (samples 20 to 900 m3, Fig. 2b) reflect mixed redox conditions where contemporaneous sulfate reduction and oxidation occur. Relatives of the sulfate reducing, oxygen tolerant bacterium Desulfosporosinus fructosivorans (Fig. 2b) were observed in these initial samples, followed by elevated populations of putative sulfur oxidizing populations in subsequent samples (Fig. 2b). This suggests that biological activation of the subsurface sulfur cycle commences rapidly following injection of surface-derived fluids during HF. Desulfosporosinus spp. harbour multiple oxygen detoxification systems allowing these bacteria to persist in aerobic environments23 together with an ability to form endospores likely conferring resistance to biocides in what could be a surface-derived dormant “seed bank” population24,25 transported into the subsurface via injection fluids where the cells are activated by new environmental conditions.
Halotolerant species are enriched following injection especially Halanaerobium spp. (Fig. SI-4), which are routinely detected in shale reservoir produced waters 14,26,27 where they are suspected to metabolise guar gum and/or produce sulfides via thiosulfate reduction26. Thiosulfate was not detected in produced waters, suggesting organotrophic metabolism by Halanaerobium led to enrichment of these populations. Relatives of the halotolerant Pseudomonas songnenensis (Fig. 2b)28 that is capable of aerobic hydrocarbon degradation29 and Methylophaga spp. (Fig. SI-4) observed at hydrothermal vents and oil spill sites30 capable of degrading higher molecular weight hydrocarbons31 were also observed.
This mixed anoxic-hypoxic phase may be required for initiating subsurface biological sulfur cycling, given that many sulfate reducing microorganisms do not directly metabolize n-alkanes and instead require metabolic residues such as organic acids and alcohols generated as by-products of organotrophic metabolism including the aerobic degradation of hydrocarbons32. Acetate and propionate in produced fluids (Fig. SI-5) may result from the fermentation/degradation of organic compounds used in HF additives33 and/or in situ hydrocarbons or via biodegradation or thermal cracking of kerogen at temperatures between 80 and 200°C34. Increasing acetate concentrations over the sampling period suggest a fermentative source. Fermentation of carbohydrates by Halanaerobium spp. may produce acetate and propionate that subsequently supports the growth of sulfur utilizing microbial communities35. This hints at the importance of an aerobic phase to provide the necessary reductants for future anaerobic populations to take hold.
Re-establishment of anoxic conditions
Over time as HF chemicals are degraded and specific microbial populations become enriched conditions evolve from a mixed redox environment to an anoxic one. Over a 12-month period, we observed large changes to microbial populations and their impacts on fluids as revealed by geochemical and isotopic analyses. After the first three months iron concentrations decreased to below the detection limit coinciding with the onset of measurable H2S in the produced water (Fig. 3a). Dissolved iron is removed via sulfide precipitation reactions, signaling a change in redox state in the subsurface.
The isotopic composition of sulfate during the initial sampling period (20–900 m3 Fig. 2a) and those at the beginning of monthly sampling (Fig. 3a) resemble δ34SSO4 and δ18OSO4 values of in situ anhydrite minerals within the Montney Formation36. This coupled with relatively constant sulfate concentrations in produced waters through time indicates that anhydrite dissolution is the principal source of sulfate. Over the course of the monthly sampling δ34SSO4 and δ18OSO4 values of dissolved sulfate increased by 7.8‰ and 2.1‰, respectively (Fig. 3a) coincident with the onset of H2S production. This shift in δ34SSO4 and δ18OSO4 values is indicative of kinetic S and O isotope fractionation via semi-closed Rayleigh distillation during MSR. δ34SSO4 and δ18OSO4 values increase during MSR as enzymatic processes within the cell preferentially metabolise bonds with lighter sulfur (32S) and oxygen (16O) isotopes resulting in a residual pool of sulfate becoming progressively enriched in 34S and 18O37,38. It has been postulated that the magnitude of this isotope fractionation may be a function of microbial community and metabolites39 and the availability of sulfate 40,41. Additionally, laboratory and field studies indicate that higher MSR rates are associated with decreased sulfur37,42,43 and oxygen44 isotope fractionation. Moreover, lower oxygen isotope fractionations compared to sulfur isotope fractionations have been associated with higher sulfate reduction rates 44,45. The ratio of sulfur to oxygen isotope fractionation in remaining sulfate of produced waters is remarkably low (0.358) corresponding to a rapid sulfate reduction rate of approximately 3.1·10− 5 mol·cm− 3·year− 1 44 (Fig. 4) implying minimal sulfate recycling via disproportionation. A similar sulfur-oxygen isotope fractionation relationship and sulfate reduction rate was observed at an oil seep in the Gulf of Mexico42 suggesting that oil seeps may represent good analogues for processes occurring in the subsurface during HF operations. Indeed, many of the microbes (e.g., Pseudomonas songnenensis, Methylophaga thalassica) observed in this study are similar to those recovered from oil seeps and oil degradation sites31.
Sulfide δ34S values ranged from 17.1‰ to 19.7‰ equating to an average sulfur isotope enrichment factor of − 16.5‰ (δ34SSO4 – δ34SH2S). Sulfur isotope enrichment factors associated with MSR in natural populations typically range from − 16 to − 43‰46, but can be as large as − 72‰45. Smaller S isotope fractionation is associated with higher sulfate reduction rates coupled with ample sulfate supply46. The relatively small S isotope fractionation between sulfate and sulfide in this system implies high sulfate reduction rate, supporting the δ34SSO4–δ18OSO4 relationship described above. Desulfosporosinus spp. are the dominant sulfate reducing microorganism within most of the monthly samples (Fig. 3b), however species diversity increases toward later samples as the system evolves towards more anoxic conditions. Despite relatively low abundances compared to other bacteria, Desulfosporosinus is capable of high sulfate reduction rate47 making it a keystone species having a disproportionally large effect compared with its abundance48, here estimated to be 2% or less.
In concert with increasing δ34SSO4 we observed an increase in the diversity and relative abundance of sulfate reducing microorganism and higher H2S concentrations in produced gas (Fig. 3a,b). The abundance of sulfate reducing microorganisms increases until 7 months post-injection, coincident with peak H2S concentrations in produced gas. Following this peak, the abundance of sulfate reducing microorganisms and H2S concentrations decreased potentially signalling the end of nutrient availability or the system poisons itself. Acetate, propionate, and sulfate concentrations remain constant (Fig. SI-5) potentially favouring the system poisoning hypothesis. Together, this data suggests that rapid sulfate reduction via MSR persists in these environments despite ‘hot’ (> 90°C) ambient temperatures in excess of previously assumed MSR thresholds.