Deep Subsurface Research Articles

Artificial subsurface lithoautotrophic microbial ecosystems and gas storage in deep subsurface.

Over the next few years, it is planned to convert all or part of the underground gas storage (UGS) facilities used for natural gas (salt caverns, depleted hydrocarbon reservoirs, and deep aquifers) into underground dihydrogen (H2) storage reservoirs. These deep environments host microbial communities, some of which are hydrogenotrophic (sulfate reducers, acetogens, and methanogens). The current state of microbiological knowledge is thus presented for the three types of UGS facilities. In the mid-1990s, the concept of anaerobic subsurface lithoautotrophic microbial ecosystems, or SLiMEs, emerged. It is expected that the large-scale injection of H2 into subsurface environments will generate new microbial ecosystems called artificial SLiMEs, which could persist over time. These artificial SLiMEs could lead to H2 loss, an intense methanogenic activity, a degradation of gas quality and a risk to installations through sulfide production. However, recent studies on salt caverns and deep aquifers suggest that hydrogenotrophic microbial activity also leads to alkalinization (up to pH 10), which can constrain hydrogenotrophy. Therefore, studying and understanding these artificial SLiMEs is both a necessity for the development of the H2 industry and presents an opportunity for ecologists to monitor the evolution of deep environments in real time.

Spatial and temporal dynamics of groundwater biogeochemistry in the deep subsurface of a high-energy beach

Intertidal sandy beach systems are considered complex biogeochemical reactors. At beach sites that are subject to high tidal and wave energy, seawater circulation can reach tens of meters deep into the subsurface and changing environmental conditions are assumed to lead to dynamic groundwater flow paths, saltwater-freshwater mixing zones, and a spatio-temporally variable groundwater biogeochemistry. Previous studies mainly focused on the upper meters of subterranean estuaries (STE), while the deep subsurface remained a black box. This study presents spatial (cross-shore) and temporal (~six-weekly, over 1.5 years) dynamics of the groundwater biogeochemistry that were observed down to 24 m below the ground surface (mbgs) of a sandy high-energy beach on Spiekeroog Island (Germany).In addition to redox conditions along a cross-shore transect ranging from oxic to Fe oxide reducing/slightly sulfidic, we found a previously unknown, distinct vertical redox zonation as well. Temporal variations of the biogeochemistry within low salinity groundwater at the most landward station close to the dune base were mainly driven by storm flood related seawater infiltration. Around the high water line, the extent of the upper saline plume (USP) varied over time. Furthermore, temporal dynamics of the O2 saturation at 6 mbgs indicated a seasonally shifting depth of the oxycline at this location. In the lower intertidal zone, groundwater solute concentrations displayed a temporally variable zone of deep freshwater discharge.Regarding the impact of the deep STE on the groundwater biogeochemistry of the discharge zone, our data revealed that nutrient, Mn, and Fe release along the deep flow paths through the USP towards the discharge zone was limited, likely due decreasing availability of labile organic matter and subsequent slowing down of metabolic processes with depth. High concentrations of metabolites in the upper ~2 mbgs of the discharge zone were, therefore, rather attributed to the incorporation of labile organic matter during continuous and storm flood related sediment relocation and/or the contribution of older waters, e.g., the subtidal saltwater wedge.

BROADBAND SEISMIC NODE BASED ON LOW-NOISE HIGH-SENSITIVITY MOLECULAR-ELECTRONIC SENSORS

The developed wideband compact digital seismic module designed for use in ground nodal seismic systems with deployment assistance from unmanned aerial vehicle (UAV) is discussed. The operation frequency band is 1–300 Hz, which is by far wider than that of the seismic nodes based on using standard (10Hz) electrodynamic geophones. Such improvement was achieved through the use of highly sensitive seismic sensors with molecular electronic transfer (MET) that allow for covering a low-frequency part of the spectrum that is needed for broadband processing and receiving information on deep subsurface formation. Low self-noise of the MET seismic sensors used for this seismic station has been proven in special tests. According to the test results, measurements of the self-noise of the developed sensor turned out to be 27-30 dB lower than the self-noise of similar sensors with additional mass. Besides of the 3 mutually orthogonal wideband MET sensors the developed module consists of a digital electronics block enabling autonomous registration and featuring wireless data transmission capability; a battery compartment providing autonomous operation of the seismic module; a device for attaching the seismic module to UAVs and lowering it to the ground. The developed device ensures data time synchronization with satellites, records data at frequencies up to 1000 samples per second with 24 bits resolution and stores the recorded information on the internal 32 Gb storage. A distinctive feature of the presented seismic node is the implementation of quality control of seismic data using a drone. This function is implemented on a Raspberry Pi single-board computer installed on board of the UAV, which has the ability to retrieve portion of information from the node for data quality control and initiate calibration procedure.

Envelope normalized reflection waveform inversion

AbstractThe reflection waveform inversion has the capability to reconstruct the background velocity model using only the reflection data by employing a migration/demigration process. Utilizing the waveform discrepancy to update the background velocity model, the conventional reflection waveform inversion method heavily relies on the true‐amplitude migration/demigration technique to reproduce the primary amplitude information from the observed reflections. We can reproduce the amplitude of observed reflections by performing least‐squares reverse time migration to estimate the reflectivity in each iteration. However, this strategy is quite time‐consuming. To avoid the need for the true‐amplitude migration/demigration or least‐squares reverse time migration, we develop an amplitude‐independent reflection waveform inversion method that uses an envelope‐normalized objective function. The envelope‐normalized waveform difference can extract the phase residuals accurately as a function of time. Compared with the global energy–normalized misfit, our proposed envelope‐normalized objective function is essentially a phase‐matched measurement. At the same time, due to the amplitude independence of our proposed objective function, the subsequent weak reflections contribute with a similar weight to the total value of the misfit as the strong early reflections do. This makes it possible to recover the deep subsurface velocity. Synthetic data of the Sigsbee model and marine streamer field data applications validate that our amplitude‐independent reflection waveform inversion method can further improve the resolution and accuracy by aligning the reflection events of synthetic and observed data phase to phase without the need to perform true‐amplitude migration/demigration or least‐squares reverse time migration as in conventional reflection waveform inversion.

13C-enriched carbonate precipitates reveal intense methanogenic oil degradation in the upper Wuerhe Formation, Northwest China

Abstract Microbial methanogenesis from crude oil is an important source of CH4 for gas reservoirs and atmospheric greenhouse gases. However, its petrological records have not been found in natural environments, and the geological conditions under which it may occur remain unclear. Here, we provide the first petrological evidence of Fe(III)-mediated methanogenic hydrocarbon degradation in a deep subsurface oil reservoir in Northwest China. The major findings are as follows: (1) highly positive δ13C values (up to +16‰) of secondary calcite attributed to methanogenesis; (2) paragenetic relation of high-δ13C calcite to biodegraded hydrocarbon; and (3) remarkably high FeO contents (up to 8 wt%) and heavy δ56Fe ratios (up to +0.52‰) in calcite, indicative of microbial Fe(III) reduction. Our study shows that methanogenic hydrocarbon degradation can occur in Fe(III)-reducing environments. This process transformed hydrocarbons into CO2 and CH4, where the former mostly precipitated as Fe-rich calcite (the carbon sink), while the latter, representing an estimated ~1968 Tg, might have escaped into the overburden and atmosphere from the Permian reservoir during the Late Triassic to Early Jurassic, which may have acted as an important CH4 source in changing global climate in the geological past.

Metabolic adaptations underpin high productivity rates in relict subsurface water

Groundwater aquifers are ecological hotspots with diverse microbes essential for biogeochemical cycles. Their ecophysiology has seldom been studied on a basin scale. In particular, our knowledge of chemosynthesis in the deep aquifers where temperatures reach 60 °C, is limited. Here, we investigated the diversity, activity, and metabolic potential of microbial communities from nine wells reaching ancient groundwater beneath Israel’s Negev Desert, spanning two significant, deep (up to 1.5 km) aquifers, the Judea Group carbonate and Kurnub Group Nubian sandstone that contain fresh to brackish, hypoxic to anoxic water. We estimated chemosynthetic productivity rates ranging from 0.55 ± 0.06 to 0.82 ± 0.07 µg C L−1 d−1 (mean ± SD), suggesting that aquifer productivity may be underestimated. We showed that 60% of MAGs harbored genes for autotrophic pathways, mainly the Calvin–Benson–Bassham cycle and the Wood–Ljungdahl pathway, indicating a substantial chemosynthetic capacity within these microbial communities. We emphasize the potential metabolic versatility in the deep subsurface, enabling efficient carbon and energy use. This study set a precedent for global aquifer exploration, like the Nubian Sandstone Aquifer System in the Arabian and Western Deserts, and reconsiders their role as carbon sinks.

Methylomonas rivi sp. nov., Methylomonas rosea sp. nov., Methylomonas aurea sp. nov. and Methylomonas subterranea sp. nov., type I methane-oxidizing bacteria isolated from a freshwater creek and the deep terrestrial subsurface.

Four methane-oxidizing bacteria, designated as strains WSC-6T, WSC-7T, SURF-1T, and SURF-2T, were isolated from Saddle Mountain Creek in southwestern Oklahoma, USA, and the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The strains were Gram-negative, motile, short rods that possessed intracytoplasmic membranes characteristic of type I methanotrophs. All four strains were oxidase-negative and weakly catalase-positive. Colonies ranged from pale pink to orange in colour. Methane and methanol were the only compounds that could serve as carbon and energy sources for growth. Strains WSC-6T and WSC-7T grew optimally at lower temperatures (25 and 20 °C, respectively) compared to strains SURF-1T and SURF-2T (40 °C). Strains WSC-6T and SURF-2T were neutrophilic (optimal pH of 7.5 and 7.3, respectively), while strains WSC-7T and SURF-1T were slightly alkaliphilic, with an optimal pH of 8.8. The strains grew best in media amended with ≤0.5% NaCl. The major cellular fatty acids were C14 : 0, C16 : 1 ω8c, C16 : 1 ω7c, and C16 : 1 ω5c. The DNA G+C content ranged from 51.5 to 56.0 mol%. Phylogenetic analyses indicated that the strains belonged to the genus Methylomonas, with each exhibiting 98.6-99.6% 16S rRNA gene sequence similarity to closely related strains. Genome-wide estimates of relatedness (84.5-88.4% average nucleotide identity, 85.8-92.4% average amino acid identity and 27.4-35.0% digital DNA-DNA hybridization) fell below established thresholds for species delineation. Based on these combined results, we propose to classify these strains as representing novel species of the genus Methylomonas, for which the names Methylomonas rivi (type strain WSC-6T=ATCC TSD-251T=DSM 112293T), Methylomonas rosea (type strain WSC-7T=ATCC TSD-252T=DSM 112281T), Methylomonas aurea (type strain SURF-1T=ATCC TSD-253T=DSM 112282T), and Methylomonas subterranea (type strain SURF-2T=ATCC TSD-254T=DSM 112283T) are proposed.

Offset VSP acquired with a single-use fiber-optic probe: Gorgon CCS case study

There are many challenges associated with monitoring industrial subsurface CO2 storage. Utilizing a reliable strategy for monitoring and assurance is paramount. Numerous aspects and trade-offs need to be considered when designing a cost-effective monitoring solution that must satisfy regulatory and social licensing requirements and minimize operational complexity. Among various geophysical techniques that are included in the suite of carbon capture and storage (CCS) monitoring programs, seismic surveys provide reliable observations of the deep subsurface. However, active surface seismic techniques are expensive, resulting in significant time separation between surveys (several months or years). Downhole seismic approaches combined with fiber-optic sensing provide flexible options in designing on-demand monitoring strategies. Distributed acoustic sensing offers the opportunity to utilize injection wells for downhole time-lapse seismic imaging and microseismic detection, minimizing the impact on production operations. This paper presents the results of a multioffset vertical seismic profiling survey acquired in CO2 injection wells using bare fibers. The fibers were deployed inside the production tubing through FiberLine Intervention technology developed by Well-SENSE. The test utilized two injection wells for active seismic acquisition. The study took place on Barrow Island (Western Australia) as part of the Gorgon CCS project. The study demonstrates the technological advantages of fast deployment of the sensing fiber and the high quality of acquired seismic data. The paper highlights characteristic noise patterns in the recorded data and suggests approaches to attenuate their effects on seismic image quality. Finally, the impact of survey geometry and directional sensitivity of fiber on the recorded wavefield and the quality of seismic imaging is discussed.

Heavy Mineral and Zircon Age Constraints on Provenance of Cenozoic Sandstones in the Gulf of Mexico Subsurface

Combined heavy mineral analysis and detrital zircon geochronology have enabled us to track detritus supplied by the ancestral river systems draining the North American continent into the deep subsurface of the Gulf of Mexico, in both the coastal plain and the offshore deep water areas. During deposition of the Paleocene–Eocene Wilcox Group, sandstones in the western part of the area are interpreted as the products of the Rosita system derived via paleo-Rio Grande material, with a large component of sediment shed from the Western Cordillera. By contrast, samples from wells further east have high proportions of zircons derived from the Yavapai-Mazatzal Province and are attributed to the Rockdale system with sediment fed predominantly by the paleo-Colorado or paleo-Colorado-Brazos. There is evidence that sediment from the Rosita system occasionally extended into the central Gulf of Mexico, and, likewise, data indicate that the Rockdale system sporadically supplied sediment to the western part of the basin. During the Late Eocene of the central Gulf of Mexico (Yegua Formation) there was a distinct shift in provenance. The earlier Yegua sandstones have a large Grenville zircon component and are most likely to have had a paleo-Mississippi origin, whereas the later Yegua sandstones are dominated by zircons of Western Cordilleran origin, similar to Wilcox sandstones fed by the Rosita system via the paleo-Rio Grande. The switch from paleo-Mississippi to paleo-Rio Grande sourcing implies there was a major reorganisation of drainage patterns during the Late Eocene. Miocene sandstones in the deepwater Gulf of Mexico were principally sourced from the paleo-Mississippi, although the paleo-Red River is inferred to have contributed to the more westerly-located wells.

Non-thermodynamic factors affect competition between thermophilic chemolithoautotrophs from deep-sea hydrothermal vents.

The deep subsurface is the largest reservoir of microbial biomass on Earth and serves as an analog for life on the early Earth and extraterrestrial environments. Methanogenesis and sulfur reduction are among the more common chemolithoautotrophic metabolisms found in hot anoxic hydrothermal vent environments. Competition between H2-oxidizing sulfur reducers and methanogens is primarily driven by the thermodynamic favorability of redox reactions with the former outcompeting methanogens. This study demonstrated that competition between the hydrothermal vent chemolithoautotrophs Methanocaldococcus jannaschii, Methanothermococcus thermolithotrophicus, and Desulfurobacterium thermolithotrophum is also influenced by other overlapping factors such as staggered optimal growth temperatures, stochasticity, and hydrology. By modeling all aspects of microbial competition coupled with field data, a better understanding is gained on how methanogens can outcompete thiosulfate reducers in hot anoxic environments and how the deep subsurface contributes to biogeochemical cycling.

Microbial anabolic and catabolic utilization of hydrocarbons in deep subseafloor sediments of Guaymas Basin.

Guaymas Basin, located in the Gulf of California, is a hydrothermally active marginal basin. Due to steep geothermal gradients and localized heating by sill intrusions, microbial substrates like short-chain fatty acids and hydrocarbons are abiotically produced from sedimentary organic matter at comparatively shallow depths. We analyzed the effect of hydrocarbons on uptake of hydrocarbons by microorganisms via nano-scale secondary ion mass spectrometry (NanoSIMS) and microbial sulfate reduction rates (SRR), using samples from two drill sites sampled by IODP Expedition 385 (U1545C and U1546D). These sites are in close proximity of each other (ca. 1km) and have very similar sedimentology. Site U1546D experienced the intrusion of a sill that has since then thermally equilibrated with the surrounding sediment. Both sites currently have an identical geothermal gradient, despite their different thermal history. The localized heating by the sill led to thermal cracking of sedimentary organic matter and formation of potentially bioavailable organic substrates. There were low levels of hydrocarbon and nitrogen uptake in some samples from both sites, mostly in surficial samples. Hydrocarbon and methane additions stimulated SRR in near-seafloor samples from Site U1545C, while samples from Site U1546D reacted positively only on methane. Our data indicate the potential of microorganisms to metabolize hydrocarbons even in the deep subsurface of Guaymas Basin.

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

AbstractTerrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m2/year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m2/year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m2/year, compared to 0.12 ± 0.02 gC/m2/year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m2/year) than in the deep subsurface (0.36 ± 0.02 gC/m2/year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

AbstractThe hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three‐dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high‐frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia‐Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river‐stage fluctuations. Adding 6.64 × 106 m3 year−1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5‐m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

Impact of mineral reactions and surface complexation on the transport of dissolved species in a subterranean estuary: Application of a comprehensive reactive transport modeling approach

Subterranean estuaries (STE) are hotspots of biogeochemical reactions. Here, dissolved constituents in waters of terrestrial and marine origin are transformed before they discharge to the coastal oceans. The involved biogeochemical reactions are complex and non-linear, calling for the application of numerical reactive transport modeling (RTM) to improve the process understanding. The aim of this study was to assess the roles of organic matter degradation and coupled secondary mineral reactions for the fate of dissolved species in STEs of sandy beaches. A comprehensive RTM approach was applied for this purpose, accounting for the effects of ion activities, pH, pe, redox reactions, mineral equilibria (calcite, goethite, siderite, iron sulfide, hydroxyapatite and vivianite) as well as surface complexation. Results show that the STE biogeochemistry and associated species fluxes are very sensitive to the assumed reaction network. For example, inorganic carbon and pH were largely controlled by calcite and siderite dynamics, and dissolved Fe2+ and HS- were precipitated as goethite, siderite and/or iron sulfides. Moreover, PO43- concentrations were affected by both the formation of vivianite or hydroxyapatite as well as surface complexation. This work helped to establish the relative importance of some of the major biogeochemical processes in the STE. However, further field studies are needed to understand which processes play a role in real-world STEs, including an exploration of the deep subsurface of STEs. Such field-based observations will improve our conceptual process understanding, which is key to developing well-constrained RTMs.

Formation mechanism of carbonate cement in tight sandstone reservoirs in the depression zone of foreland basin and its impact on reservoir heterogeneity: The upper Triassic Xujiahe Formation in Western Sichuan foreland basin, China

The carbonate cements in tight sandstone reservoirs are diverse, complex, and important determinants of reservoir quality. This study provides a comprehensive analysis of the characteristics and formation mechanism of carbonate cement in the tight sandstone reservoirs of the Xu3 Member within the Xujiahe Formation in the depression zone of the Western Sichuan foreland basin. Core samples were examined on thin sections using cathodoluminescence and scanning electron microscopy, fluid inclusion analysis, in situ carbon and oxygen stable isotopes, and rare earth element analysis. The study reveals that the Xu3 Member primarily comprises sublitharenite and litharenite, which has three stages of carbonate cementation within the reservoirs. On average, the reservoirs exhibit a porosity of 3.82% and a permeability of 0.1mD. In the eodiagenetic stage, calcite and dolomite (δ13CPDB = 0‰ to −4.07‰; δ18OPDB = −5.64‰ to −12.69‰) formed, filling both intergranular pores and intragranular dissolution pores. This process is driven by the expulsion of organic acids through compaction from the clay and coal-bearing sections during shallow burial. The carbon sources originate from external organic acids and internal carbonate rock fragments, which combine with Ca2+ and Mg2+ ions released from the dissolution of carbonate rock fragments. During the mesodiagenetic stage A, the predominant formation of ferroan calcite and ankerite (δ13CPDB = −5.98‰ to −7.64‰; δ18OPDB = −14.47‰ to −19.87‰) occurred, filling intergranular residual pores and secondary intragranular dissolution pores. This phenomenon is caused by the expulsion of organic acids formed by the compaction of matured organic matter at greater burial depths. The carbon sources primarily originate from external organic acids and combine with Ca2+, Mg2+, and Fe2+ ions in the sandstone, alongside accompanying organic acids. During the mesodiagenetic stage B, coarse-macro crystalline calcites (δ13CPDB = −1.06‰ to −1.64‰; δ18OPDB = −17.95‰ to −20.36‰) with significant positive Eu anomaly-filled fractures were formed. These calcites are associated with hydrothermal fluids in the deep subsurface, precipitating directly in the fractures due to the presence of hot hydrothermal fluids rich in Ca2+ and CO32−. The Xu3 Member displays two sandstone superimposed types. In the mudstone overlain by sandstone type, carbonate cement is mainly present in the middle of the cycle, and reservoir quality is relatively high at both the bottom and top of the cycle. In the sandstone overlain by sandstone type, carbonate cement is predominantly found in the middle and top sections of the cycle, with higher-quality reservoirs primarily forming at the bottom of the cycle. It resulted in significant heterogeneity observed in the tight sandstone reservoirs of the depression zone in the Western Sichuan foreland basin.

Single-cell analysis reveals an active and heterotrophic microbiome in the Guaymas Basin deep subsurface with significant inorganic carbon fixation by heterotrophs.

The global subsurface is the largest reservoir of microbial life on the planet yet remains poorly characterized. The activity of life in this realm has implications for long-term elemental cycling, particularly of carbon, as well as how life survives in extreme environments. Here, we recovered cells from the deep subsurface of the Guaymas Basin and investigated the level and distribution of microbial activity, the physicochemical drivers of activity, and the relative significance of organic versus inorganic carbon to subsurface biomass. Using a sensitive single-cell assay, we found that the majority of cells are active, that activity is likely driven by the availability of energy, and that although heterotrophy is the dominant metabolism, both organic and inorganic carbon are used to generate biomass. Using a new approach, we quantified inorganic carbon assimilation by heterotrophs and highlighted the importance of this often-overlooked mode of carbon assimilation in the subsurface and beyond.

Illuminating the “Invisible”: Substantial Deep Respiration and Lateral Export of Dissolved Carbon From Beneath Soil

AbstractDissolved organic and inorganic carbon (DOC and DIC) influence water quality, ecosystem health, and carbon cycling. Dissolved carbon species are produced by biogeochemical reactions and laterally exported to streams via distinct shallow and deep subsurface flow paths. These processes are arduous to measure and challenge the quantification of global carbon cycles. Here we ask: when, where, and how much is dissolved carbon produced in and laterally exported from the subsurface to streams? We used a catchment‐scale reactive transport model, BioRT‐HBV, with hydrometeorology and stream carbon data to illuminate the “invisible” subsurface processes at Sleepers River, a carbonate‐based catchment in Vermont, United States. Results depict a conceptual model where DOC is produced mostly in shallow soils (3.7 ± 0.6 g/m2/yr) and in summer at peak root and microbial respiration. DOC is flushed from soils to the stream (1.0 ± 0.2 g/m2/yr) especially during snowmelt and storms. A large fraction of DOC (2.5 ± 0.2 g/m2/yr) percolates to the deeper subsurface, fueling deep respiration to generate DIC. DIC is exported predominantly from the deeper subsurface (7.1 ± 0.4 g/m2/yr, compared to 1.3 ± 0.3 g/m2/yr from shallow soils). Deep respiration reduces DOC and increases DIC concentrations at depth, leading to commonly observed DOC flushing (increasing concentrations with discharge) and DIC dilution patterns (decreasing concentrations with discharge). Surprisingly, respiration processes generate more DIC than weathering in this carbonate‐based catchment. These findings underscore the importance of vertical connectivity between the shallow and deep subsurface, highlighting the overlooked role of deep carbon processing and export.

Microbial ecology of the deep terrestrial subsurface.

The terrestrial subsurface hosts microbial communities that, collectively, are predicted to comprise as many microbial cells as global surface soils. Although initially thought to be associated with deposited organic matter, deep subsurface microbial communities are supported by chemolithoautotrophic primary production, with hydrogen serving as an important source of electrons. Despite recent progress, relatively little is known about the deep terrestrial subsurface compared to more commonly studied environments. Understanding the composition of deep terrestrial subsurface microbial communities and the factors that influence them is of importance because of human-associated activities including long-term storage of used nuclear fuel, carbon capture, and storage of hydrogen for use as an energy vector. In addition to identifying deep subsurface microorganisms, recent research focuses on identifying the roles of microorganisms in subsurface communities, as well as elucidating myriad interactions-syntrophic, episymbiotic, and viral-that occur among community members. In recent years, entirely new groups of microorganisms (i.e. candidate phyla radiation bacteria and Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoloarchaeota, Nanoarchaeota archaea) have been discovered in deep terrestrial subsurface environments, suggesting that much remains unknown about this biosphere. This review explores the historical context for deep terrestrial subsurface microbial ecology and highlights recent discoveries that shape current ecological understanding of this poorly explored microbial habitat. Additionally, we highlight the need for multifaceted experimental approaches to observe phenomena such as cryptic cycles, complex interactions, and episymbiosis, which may not be apparent when using single approaches in isolation, but are nonetheless critical to advancing our understanding of this deep biosphere.

Draft genome sequence of the deep-subsurface Ciceribacter sp. strain T2.26MG-112.2, a second Rhizobiaceae isolated from the Iberian Pyrite Belt at 492.6 mbs.

T2.26MG-112.2 is a Ciceribacter strain that has been isolated from the deep subsurface of the Iberian Pyrite Belt. We report its draft genome consisting of a chromosome of ≈4.9 Mb and a plasmid of 357 kb. The annotation reveals 4,824 coding sequences, 48 tRNA genes, and 1 rRNA operon.

Acceleration of Deep Subsurface Fluid Fluxes in the Anthropocene

AbstractThe Anthropocene has been framed around humanity's impact on atmospheric, biologic, and near‐surface processes, such as land use and vegetation change, greenhouse gas emissions, and the above‐ground hydrologic cycle. Groundwater extraction has lowered water tables in many key aquifers but comparatively little attention has been given to the impacts in the deeper subsurface. Here, we show that fluid fluxes from the extraction and injection of fluids associated with oil and gas production and inflow of water into mines likely exceed background flow rates in deep (>500 m) groundwater systems at a global scale. Projected carbon capture and sequestration (CCS), geothermal energy production, and lithium extraction to facilitate the energy transition will require fluid production rates exceeding current oil and co‐produced water extraction. Natural analogs and geochemical modeling indicate that subsurface fluid manipulation in the Anthropocene will likely appear in the rock record. The magnitude and importance of these changes are unclear, due to a lack of understanding of how deep subsurface hydrologic and geochemical cycles and associated microbial life interact with the rest of the Earth system.