Abstract
Accessing the deep terrestrial subsurface (greater than 1 km below the surface) presents significant practical challenges, leaving these ecosystems largely uncharacterized despite their extensive presence beneath Earth's landmasses. In this study, we introduce the BedrettoLab Deep Life Observatory (DELOS), a new underground laboratory to study the biogeochemical diversity of groundwater in a high-altitude Alpine catchment tens of meters to 1.6 km underground. Biogeochemical monitoring of DELOS over spatial and temporal scales highlight three dominant ecotypes throughout DELOS: (1) Shallow groundwater with low electrical conductivity enriched in Leptospirillia; (2) High-inflow fault zones enriched in ultra-small bacteria and archaea; (3) Bicarbonate-enriched waters that are enriched in Candidatus Kryptonia and Spirochaetota. Despite a consistent lithology throughout DELOS, groundwater from fractures that are spatially near each other are not always represented by the same ecotype and can be more similar to groundwater emitted from fractures thousands of meters away. Despite this heterogeneity, the biological and hydrochemical compositions of the groundwater of individual fractures remained relatively stable throughout the course of a 1-year monitoring period. An exception to this trend occurred after a series of seismic events near one groundwater-bearing fracture. Here, the microbial community and hydrochemical composition of the groundwater changed after the seismic events but returned to the site's “baseline” composition within 3 weeks. Taken together, these findings provide new insights into the spatial and temporal heterogeneity of deep subsurface ecosystems and the subsurface connectivity of an Alpine subsurface environment.
1 Introduction
Since the earliest reports of subsurface life in coal beds and oil fields in the 1920s (Bastin et al., ; Ginsburg-Karagitscheva, ; Lipman, ), the field of subsurface geomicrobiology has grown significantly (for a review, see Onstott, ) and it is now estimated that more bacteria and archaea reside in the continental subsurface than any other microbial habitat (for a review, see Magnabosco et al., ). All three domains of life and even animals have been found to be capable of surviving hundreds of meters to kilometers underground (Borgonie et al., ) and, in many cases, survive exclusively through the utilization of geologically-generated energy sources such as H2 and CO2 (for a review, see Magnabosco et al., ). Due to the practical challenges of accessing and studying deep (>500 m) subsurface ecosystems, a variety of subsurface experiments and monitoring programs have been pursued in underground laboratories such as the Äspö Hard Rock Laboratory (Äspö HRL, Sweden) (e.g. Kotelnikova and Pedersen, ; Pedersen, ; Kyle et al., ; Westmeijer et al., 2022), Mizunami Underground Research Laboratory (Japan) (e.g. Fukuda et al., ; Ino et al., ), Grimsel Test Site (Switzerland) (e.g. Konno et al., ), Mont Terri Underground Rock Laboratory (Switzerland) (e.g. Stroes-Gascoyne et al., 2007; Poulain et al., ; Bagnoud et al., ) and Sanford Underground Research Facility (USA) (e.g. Osburn et al., ; Casar et al., ; Momper et al., ). These research programs have revealed diverse communities of subsurface microorganisms capable of performing a wide range of biogeochemical functions such as methane, nitrogen, sulfur, and iron cycling (e.g. Pedersen, ; Fukuda et al., ; Konno et al., ; Ino et al., ; Wu et al., 2017; Lopez-Fernandez et al., , ). This study describes the BedrettoLab Deep Life Observatory (DELOS), a new granite-hosted underground laboratory for subsurface microbiology investigations operated by ETH Zürich and located within the Bedretto Tunnel (Switzerland).
The Bedretto Tunnel was originally constructed as an adit for Furka Base Railway Tunnel from 1971 to 1982 and intersects 3 geological units in the Gotthard massif (e.g. Hafner, ; Keller and Schneider, ; Rast et al., ): (1) the Tremola Series from tunnel meter (TM) 0 to TM-434; (2) the Prato Series (TM-434 to TM-1138); and (3) the Rotondo granite (TM-1138 to TM-5218). The metasedimentary Tremola Series is composed of mica-gneiss, amphibolites, schists, calc-silicate rocks, and quartzites (Hafner, ; Keller and Schneider, ) that change between rock types over meter scales (Labhart, ). The Prato Series also exhibits compositional heterogeneity of mica-gneisses, amphibolites, and schists over centimeter to meter scales (Labhart, ). A sharp contact between the Prato Series and Rotondo granite occurs at TM-1138 and large amounts of deformation in the granite is observed between TM-1,138 and TM-1,143 (Rast et al., ). The Rotondo granite is dissected by brittle faults (Lützenkirchen, ) and ductile shear zones (Rast et al., ) and is further divided into a biotite-rich porphyritic granite (RG2) in the region of TM-2,800 to TM-3,440 that is surrounded by an equigranular granite (RG1) composed of quartz, alkali feldspar, plagioclase, and biotite (Rast et al., ). 206Pb/238U Zircon dating of the Rotondo granite indicates intrusion ages of 285 to 319 Ma for RG1 and 280 to 335 Ma RG2. Groundwater residence time estimates based on tritium for natural fractures spanning approximately TM-3,800 to TM-5,218 indicate that glacially-fed fractures experience rapid recharge (1 to 1.5 years) while non-glacially-fed fractures exhibit a much more complex recharge history and contain a sub-modern component (Ofterdinger, ). It is worth noting that these hydrochemical measurements are limited to the mineralogically consistent last ~1,500 m of the Bedretto Tunnel due to a tunnel collapse that occurred while the tunnel was left un-maintained from 1982 to 2015 and the residence time of groundwater from TM 0 − 3, 800 is currently unknown.
In 2018, ETH Zürich opened the “Bedretto Underground Laboratory for Geosciences and Geoenergies” (BULGG) between the Bedretto Tunnel's TM-2,000 and TM-2,100 (Ma et al., ). The first deep (197–303 m in length) boreholes of the BedrettoLab were drilled in 2019 and ahead of the “Bedretto Reservoir Project” (BRP; May 2021 to July 2023) which executed a series of hydraulic stimulation experiments designed to aid in the development of engineered geothermal systems (Gischig et al., ; Ma et al., ). In preparation for the BRP hydraulic stimulations, the BedrettoLab Deep Life Observatory (DELOS) was initiated to study the subsurface microorganisms residing within the Gotthard Massif. DELOS was named as an homage to the ancient Athenian island that prohibited death, due to the slow turnover times observed in the continental subsurface (Onstott et al., ). Since the majority of the Bedretto Tunnel is not cemented, DELOS is not restricted to the experimental boreholes of BULGG. Instead, inflows from the densely fractured Rotondo granite and gneissic sections of the 5.2 km Bedretto Tunnel can be easily sampled and provide an in-depth view of fractured rock ecosystems from tens to ~1,600 meters below land surface (mbls). This study describes the microbiological and hydrochemical changes throughout the natural fractures of the Bedretto Tunnel over a 1-year period from 2020 to 2021 and captures a period ahead of, during, and after the first hydraulic stimulation experiments within BULGG.
2 Materials and methods
2.1 Area of Study
Surface water from above the Bedretto Tunnel (Figure 1) and groundwater from the BedrettoLab Deep Life Observatory (DELOS) and Grimsel Test site (GTS, https://grimsel.com/) were collected for this study. DELOS is accessed through the Ronco Portal of the Bedretto Tunnel (46.497518°N, 8.494992°E) in the Gotthard Massif (Switzerland, Figure 1) and is part of the Riale di Ronco Alpine catchment. The far end of the Bedretto Tunnel (Furka Portal) intersects the 15 km long Furka Base Tunnel. Only a few sections across the tunnel have been shotcreted or supported by steel arches (structural support is mainly required near the Ronco Portal), providing a unique opportunity to sample groundwater directly from fractures spanning over 5 km (Figure 1, Supplementary Data 1). The GTS is an underground research laboratory located approximately 15 km northwest of DELOS and situated within the granitic rock of the Aarmassif (Keusen et al., ). The GTS has been in operation since 1984 and has supported a variety of field-scale experiments related to radioactive waste disposal and enhanced geothermal systems. The similar geology and proximity of GTS to DELOS provided motivation to collect, analyze, and compare the microbial communities and water composition of GTS boreholes “Pinkel”, “ISC Injection 2” and “SB80001”.
Figure 1
2.2 Field Measurements And Sampling Strategy
DELOS groundwater samples were obtained from a total of 57 locations. 54 sampling sites were related to natural water inflows along the Bedretto tunnel and 3 sites were connected to tunnel inflows in smaller, uncased boreholes (Ø 10 cm, length unknown) in the tunnel sidewalls (Tunnel Meter (TM)-901, TM-755A and TM-2794, where TM refers to the distance in meters from the Val Bedretto tunnel entrance). Approximately monthly hydrochemical sampling from August 2020 to December 2021 and biological sampling from November 2020 to December 2021 are reported in this study. Additional samples were collected from three boreholes in the Grimsel Test Site in April 2021 and a variety of surface water sources in the DELOS region during the summer of 2021 (Figure 1). Temperature, pH, and electrical conductivity (EC) were directly measured on site using an electronic field probe (EXTECH, multi-parameter PCSTestr 35) that was 3-point calibrated (pH: 4.01, 7.00 and 10.00; EC: 84 μS cm−1, 147 μS cm−1 and 1413 μS cm−1). The redox potential (ORP) was measured with a KBM-100 (Voltcraft) field meter. The total alkalinity (carbonate alkalinity) of samples collected from November 2020 to March 2021 was determined through acid titration (Hach®, AL-AP Test) and flow rates (Q) were measured using a bucket-and-stopwatch approach.
2.3 Hydrochemical analyses
For major ion analysis, 2 × 5 mL filtered water samples were collected in the field, transported on ice and stored at 4°C before cation (Mg2+, Ca2+, K+, Na+, , Li+, Sr2+) and anion (, Cl−,F−, , , , Br−) analysis via ion chromatography (IC; Dionex DX-120 IC System) at the Institute of Geological Sciences ETH Zürich. Prior to ion analysis, the EC and pH of each sample was measured using the inoLab EC-pH measuring instrument from WTW. Samples that exceeded the IC working range of 0-500 μS cm−1 were diluted 10× to ensure ion concentrations were within the ranges of the standards. Water isotope samples were collected and transported to Zürich in 100 mL polyethylene bottles with airtight caps (Semadeni, CH). In Zürich, water isotope samples were filtered using a 0.45 μm filter, sealed in 3 × 2 mL glass autosampler vials and stored at 4°C until water isotope analysis was performed using a Picarro L2130-i analyzer at the Geological Institute of ETH Zürich. Raw measurements were post-processed using the Picarro software ChemCorrect which calculates the isotopic composition as δ18O or δ2H-value with respect to the Vienna Standard Mean Ocean Water (VSMOW). The local meteoric water line of DELOS was estimated from regional isotope data obtained on waterisotopes.org and derived from (Bowen and Revenaugh, ).
2.4 DNA extraction, 16S rRNA gene amplification, quantification, and sequencing
Biomass was collected on site by filtering fracture or surface water through 0.22 μm Sterivex filters (SVGPL10RC, Millipore) with a peristaltic pump. Prior to filtration, new sterile tubing was flushed with the fracture water of interest for 3–5 minutes. After filtration, Sterivex filters were stored in sterile, 50 mL falcon tubes, transported on ice to ETH Zürich the same day, and frozen at -80°C until DNA extraction. DNA extractions of samples and extraction blanks were performed following the manufacturer's standard protocol of the DNeasy PowerWater Kit (14900-100-NF, Qiagen). DNA concentrations were quantified by the Qubit 4 Fluorometer (ThermoFisher) and stored in a -80°C freezer for use in downstream applications.
All DNA samples were amplified in triplicate using the Platinum II Hot-Start PCR Master Mix Kit (14000013, Invitrogen) and following a previously defined procedure (McNichol et al., ). Briefly, a single-step PCR was performed to amplify a conserved part of the V4-V5 region of the 16S SSU rRNA gene from bacteria and archaea. The region was targeted using Golay-barcoded 515F forward primers (AATGATACGGCGACCACCGAGATCTACACGCT - XXXXXXXXXXXX - TATGGTAATTGT - GTGYCAGCMGCCGCGGTAA, with the X's denoting the variable Golay barcode location) and a 926R PCR reverse primer (CAAGCAGAAGACGGCATACGAGATAGTCAGCCAGGG - CCGYCAATTYMTTTRAGTTT) (Walters et al., 2016). After PCR, triplicates were pooled together and samples were cleaned using Agencourt AMPure XP beads (Beckman Coulter) at a 0.7× ratio and DNA concentrations were measured using the Qubit 4 Fluorometer (ThermoFisher). Using the Qubit measured amplicon concentrations, samples were normalized and pooled. The pooled library was cleaned with AMPure XP beads at a 0.7× ratio and the final, cleaned library was quantified using both the Qubit 4 and 4150 Tapestation (Agilent). The amplicon library was sequenced at the ETH Zürich Genetic Diversity Center on an Illumina MiSeq using a 500-cycle paired end V2 reagent kit (MS-102-2003, Illumina). PhiX was spiked in at a concentration of 8%. The final library loading concentration was 4 pM and custom sequencing primers were used: 515Fseq (TATGGTAATTGTGTGYCAGCMGCCGCGGTAA), 926Rseq (CGGCATACGAGATAGTCAGCCAGGGCCGYCAATTYMTTT RAGTTT), and a modified Illumina index sequence primer (AATGATACGGCGACCACCGAGATCTACACGCT).
16S rRNA gene amplicon copy numbers were quantified using a quantitative PCR (qPCR) with the Lightcycler 480 (Roche). Briefly, a solution containing the V4-V5 forward primer 515F (GTGNCAGCMGCCGCGGTAA, final concentration 0.5 μM), the V4-V5 reverse primer 926R (CCGYCAATTYMTTTRAGTTT, 0.5 μM), 5 μL SsoFast EvaGreen Supermix with Low ROX 2x (1725210, BioRad), and 1 μL Bovine Serum Albumin (10 mg mL−1) were combined with 2 μL of template DNA were prepared in triplicate and loaded onto a 96-well plate for qPCR. 16S rRNA gene standards ranging from 2.088 × 106 to 20.88 16S rRNA gene amplicon copies per μL were prepared using a consensus 16S synthetic construct and serial dilution. The qPCR program included an initialization of 95°C for 300 seconds, 40 cycles of denaturation (10 seconds at 95°C), annealing (30 seconds at 50°C) and elongation (15 seconds at 68°C) with a melting curve of 15 seconds at 95°C and 60 seconds at 55°C and one acquisition per °C from 55°C to 95°C. The fit (R2 = 0.99) and melt curve of the standards are shown in Supplementary Figures 1, 2.
2.5 Biodiversity analyses
The quality of 16S rRNA gene amplicon sequences was assessed using FastQC (Babraham Bioinformatics Group). Based on the sequence quality, the first 5 nucleotides of the forward and reverse reads were trimmed using the FASTX-trimmer, http://hannonlab.cshl.edu/fastx_toolkit/license.html. Additionally, nucleotides beyond position 230 were removed from the reverse read by the FASTX-trimmer. After trimming, Trimmomatic (Bolger et al., ) was used to quality filter the reads (SLIDINGWINDOW:100:28) and remove adapters. Paired end reads were joined using fastq-join (-p 3 -m 20; v-1.04.807), https://github.com/ExpressionAnalysis/ea-utils/blob/wiki/FastqJoin.md. Amplicon sequencing variants (ASVs) were then predicted from joined reads using the “pseudopool“ option of dada2 (Callahan et al., ) and annotated using dada2 and the Silva nr99 v138.2 dataset (Quast et al., ; Yilmaz et al., 2014). After annotation, all ASVs found in the control samples were removed from the ASV data table and omitted from all downstream analyses. Biodiversity and statistical analyses were performed in R and python with ASV data. All code is available via https://github.com/GeobiologyLab/DELOS-2021-time-series.
3 Results
3.1 Seasonal and spatial hydrochemical variability throughout DELOS
DELOS groundwaters exhibit annual mean temperatures from 5°C to 19.5°C and a pH ranging from 8 to 9.6 (Supplementary Data 2), depending on location. Groundwater temperature generally correlates with overburden with the exception of high inflow zones where the largest faults occur around tunnel meter (TM) 1306, TM-2848 and TM-4166 (Figure 2). An additional zone of high inflow is observed in a landslide scarp section between TM 300 − 400, consistent with earlier hydraulic studies (Masset and Loew, ; Halter et al., ). Nitrate was detected in 22 of the 57 DELOS locations with elevated concentrations observed within the high inflow zones surrounding TM-300, TM-1306, TM-2848, TM-3750, TM-3884 and TM-4166 and a maximum concentration of 3.5 ppm at TM-216 (Figure 2). The average, per-fracture electrical conductivity (EC) in DELOS spans 43.2 μS cm−1 (TM-1306) to 931.1 μS cm−1 (TM-901; Figure 2). Groundwaters exhibiting an EC in excess of 500 μS cm−1 are also enriched with millimolar levels of sulfate. In general, the EC measurements of all groundwater-filled fractures from tunnel meters 4, 200 − 5, 218 of DELOS exhibit very little variation throughout the 16 month hydrochemical monitoring period (mean ± standard deviation: 91 ± 12 μS cm−1). On the other hand, the EC throughout the first 2,000 DELOS tunnel meters is highly variable between fractures but relatively consistent within a fracture throughout the sampling period (Figure 2). Inter-fracture water isotope measurements also exhibit greater variability in the first 2,000 tunnel meters relative to the last 1,000 tunnel meters (Figure 2). It is worth noting that the last 1,000 tunnel meters of DELOS is likely influenced by the drainage effect of the ~15 km long Furka Base Tunnel (Amberg, ) and, therefore, the similarities in the hydrochemistry of these sites may be anthropogenically driven.
Figure 2
All DELOS groundwater and DELOS surface waters follow the local meteoric water line (LMWL, Figure 3). With the exception of TM-1494 whose water isotope values are more similar to groundwater from DELOS fractures with over 400 m greater overburden, DELOS groundwater exhibits isotopic depletion as overburden increases (Figure 3). In general, the per-fracture DELOS groundwater isotope values were relatively stable throughout the monitoring period and did not reveal strong seasonal pattern (Figure 2). Previous water isotope measurements from the Grimsel Test Site (Schneeberger et al., 2017) also follow the DELOS LMWL and are most similar to the shallowest DELOS groundwater sites (Figure 3). It is worth noting that the nearby (< 15 km away), granite hosted Grimsel Test Site (GTS) is located approximately 1,730 m above sea level (masl) (Blechschmidt and Vomvoris, ) and the GTS water isotope samples are derived from sampling locations with an overburden of 420 − 520 m (Schneeberger et al., 2017). The surface elevation of DELOS groundwater samples with similar isotopic composition as the GTS samples is approximately 1,760 masl (overburden of 280 m).
Figure 3
3.2 DELOS microbial diversity and abundance
14,508 unique amplicon sequencing variants (ASVs) spanning 64 phyla were identified throughout DELOS with the highest richness (2,944 ASVs; Shannon Diversity (H′): 5.8) observed in the November 2021 TM-1306 sample and the lowest richness (267 ASVs; H′ = 4.0) observed in the December 2020 sampling of TM-4166. DELOS microbial communities are dominated (85 to >99%) by bacteria. The most abundant bacterial phyla observed within DELOS are related to Candidatus Kryptonia (up to 30% at TM-4652), Nitrospirota (up to 55% at TM-444), Pseudomonadota (previously known as Proteobacteria; up to 70% at TM-2647) and Spirochaetota (up to 30% at TM-4652) (Supplementary Data 3). A large (>50%) population of chemolithotrophic Leptospirillia from the phylum Nitrospirota make the microbial community of the shallow, gneiss-hosted sample from TM-444 an outlier relative to other samples (Figure 4). Rather than Leptospirillia, other sites exhibiting large populations of Nitrospirota are represented by members of the nitrogen and sulfur cycling class Thermodesulfovibrionia (up to 30% at TM-1494, Figure 5). The DPANN archaeal superphylum dominates DELOS archaeal communities and members of this group are most abundant (up to 15% of the total microbial community) in samples from TM-1306, TM-2848 and TM-4447. Ultra-small bacteria within the Candidate Phylum Radiation (Patescibacteria) and bacterial class Omnitrophia within the Verrucomicrobiota phylum are also relatively abundant throughout DELOS, exhibiting their highest relative abundances at TM-4166 (up to 15%) and TM-4447 (up to 23%), respectively.
Figure 4
Figure 5
Fractures within the gneissic section of DELOS generally exhibit higher concentrations of 16S rRNA gene amplicon copies mL−1 (Q25%=7.6 × 103, Q75%=3.3 × 104) relative to granitic samples (Q25%=1.7 × 103, Q75%=6.4 × 103) (Supplementary Figure 3). An exception to this trend is observed at the granitic TM-2467 fracture (8.2 ± 2.1 × 104 16S rRNA gene amplicon copies mL−1). Here, a visible biofilm surrounds the fracture of TM-2647 (Supplementary Data 1) and may have contributed to the elevated copy numbers in the samples. Large populations of sulfur oxidizing Pseudomonadota related to Sulfurifustis (35%) and Thiobacillus (10%) that, in the case of Thiobacillus, have been linked to biofilm formation in the granite-hosted Äspö Hard Rock Laboratory (Lopez-Fernandez et al., ), further point to potential biofilm contamination within the TM-2647 sample. On the other hand, the fracture at TM-4652 intersects a quartz vein (Supplementary Data 1) and consistently exhibited the lowest concentrations (931 ± 378) of 16S rRNA gene amplicon copies mL−1 with a community largely composed of Candidatus Kryptonia and Spirochaetota.
Canonical correspondence analysis of 16S rRNA gene amplicon sequencing variants (ASVs) and environmental data from DELOS groundwater samples collected in November and December 2020 reveals that DELOS groundwater microbial communities are distributed along a δ18O, depth and temperature gradient that correlates to tunnel meter (Figure 4). As with the DELOS hydrochemical measurements, samples from the last 1,000 m exhibit more similar microbial communities than samples collected throughout the first 3,500 m of the tunnel (Figure 4). An exception to the otherwise tight clustering of samples from the last 1,000 m is a pH 9.0 sample with relatively high abundance (25%) of an uncultivated lineage (“Dadabacteria” also known as Candidatus Desulfobacterota) from TM-5132. Further separation of DELOS ecosystem types (ecotypes) are seen in the sulfate- and Thermodesulfovibrionia-enriched TM-1494 sample and high flow rate and DPANN/Patescibacteria-enriched samples (TM-1306, TM-2848, TM-4166). Together, these findings point to three major ecotypes in the DELOS groundwater system: (1) Shallow groundwater with low electrical conductivity enriched in Leptospirillia; (2) High-inflow fault zones enriched in ultra-small bacteria and archaea; (3) Bicarbonate-enriched waters predominantly represented throughout the last 1,000 tunnel meters and that are enriched in Candidatus Kryptonia and Spirochaetota.
Following the DELOS-wide survey in November and December of 2020, five fractures were selected for seasonal monitoring based on hydrochemical and sampling considerations. 16S rRNA gene analysis reveals that community membership is generally stable and that inter-group differences were significantly greater than intra-group differences (Figure 5, Supplementary Figure 4). TM-4652 exhibits the most phylum-level changes during the 12-month monitoring period, revealing a slow replacement of relatively abundant Pseudomonadota and Patescibacteria with Spirochaetota and Ignavibacteriota. A similar change in the abundance of these phyla is observed at TM-1494; however, the most dramatic change in taxonomic composition within this site occurs when the dominant population of Thermodesulfovibrionia decreases from 32% to 20% relative abundance from 5 − 21 May 2021. This decrease co-occurs with an increase in the relative abundance of Ignavibacteriota and two of unclassified bacterial ASVs that are also observed in TM-1494 throughout the various sampling time points and in TM-4652. This change in community composition follows the first hydraulic stimulation experiment in the Bedretto Reservoir Project's stimulation boreholes (ST, located near TM-2050 of BULGG) on 3–5 May 2021, seismic events near the TM-1494 fracture possibly related to the hydraulic stimulation on 6 May 2021, and seismic events near the TM-1494 fracture associated with heavy rainfall on 10 May 2021 (personal communication, BULGG). An increase in electrical conductivity that correlates to an increase in all major ion concentrations except magnesium is additionally observed at TM-1494 between 5 and 21 May 2021. A second hydraulic stimulation experiment that targeted a different depth range (interval) of the aforementioned experimental ST borehole was conducted from 22 − 24 May 2021. This experiment demonstrated a connection between the injection and production wells of the stimulation (Arnet, ) but did not result in seismic activity near the TM-1494 fault. The 1 June 2021 sampling of TM-1494 after the 22 − 24 May 2021 stimulation revealed that the fracture's hydrochemical and microbial community composition had recovered to its “baseline” state, indicating that the changes observed on 21 May were relatively short-lived (Figure 5). Interestingly, an increase in electrical conductivity along with an increase in temperature following the 22 − 24 May hydraulic stimulation can be seen in TM-1306 and may indicate that the two Bedretto Reservoir Project stimulation experiments affected two different regions of the Bedretto Reservoir (TM 2000-2100); however, a microbial sample was not collected from TM-1306 during this time period (Figure 5). Together, these observations indicate that both experimentally and naturally induced seismicity contribute to the temporal biodiversity and chemical patterns observed throughout DELOS.
3.3 DELOS groundwater and surface microbial communities exhibit minimal overlap
Three groundwater samples from the granite-hosted Grimsel Test Site (GTS, ~15 km from DELOS) and nine samples from DELOS surface waters (Figure 1) were collected for hydrochemical and microbial community analysis in spring and summer 2021, respectively. These samples represent water from two GTS experimental boreholes (“Pinkel” and “ISC Injection 2”), a GTS monitoring borehole (“SB80001”), the Ticino River, a high-altitude freshwater spring, four glacial lakes, and three samples (one above ice, two below ice) from the Witenwasseren glacier (Supplementary Data 2). Surface samples exhibited a lower mean temperature (6.8 ± 3.5°C) and pH (7.4 ± 0.6) than DELOS and GTS (Table 1). The electrical conductivity of the Ticino River is elevated (EC, 409 μS cm−1) relative to other surface sites, resulting in the high variance in EC across surface sites (Table 1). Groundwater collected from the GTS boreholes exhibit mean temperature, EC and 16S rRNA gene amplicon copies mL−1 measurements between the DELOS surface and DELOS groundwater samples (Table 1).
Table 1
| DELOS groundwater | Grimsel Test Site | DELOS surface | |
|---|---|---|---|
| Temperature (°C) | 14.1 ± 3.5 | 12.9 ± 0.4 | 6.8 ± 3.5 |
| pH | 9.0 ± 0.9 | 10.3 ± 0.4 | 7.4 ± 0.6 |
| Electrical Conductivity (μS cm−1) | 252 ± 250 | 108 ± 14 | 72 ± 144 |
| 16S rRNA copies (mL−1, Q25%) | 2.0 × 103 | 4.5 × 103 | 1.0 × 106 |
| 16S rRNA copies (mL−1, Q75%) | 8.2 × 103 | 5.0 × 105 | 4.5 × 106 |
The mean and standard deviation of temperature, pH, electrical conductivity and 16S rRNA gene amplicon copies per mL collected from DELOS groundwater, Grimsel Test Site groundwater and surface locations surrounding DELOS (DELOS surface).
Q25% = 25% Quartile. Q75% = 75% Quartile.
The GTS groundwater samples captured three distinct microbial communities: (1) a Desulfitobacteriia-dominated (phylum Bacillota) community within the Grimsel In-Situ Stimulation and Circulation Experiment experimental borehole (Doetsch et al., , “ISC Injection 2”) (2) a Pseudomonadota-dominated community composed primarily of Hydrogenophaga and Rhodocyclaceae in a radionuclide experimental borehole (Quinto et al., , “Pinkel”) and (3) a Pseudomonadota- and Thermodesulfovibrionia-dominated community in the monitoring borehole (“SB80001”). Two Silva-defined bacterial phyla (Caldisericota and WS1) were identified in GTS at < 0.2% relative abundance but not DELOS and 29 phyla are found in DELOS but not GTS. On the other hand, 8 Silva-defined bacterial phyla (Caldisericota, Campylobacterota, Fusobacteriota, Rhodothermota, Synergistota, LCP-89, TX1A-33, and WS4) were identified in the surface samples surrounding DELOS but not DELOS, while 21 bacterial phyla were found in DELOS but not the surface water samples. At the phylum level, surface microbial communities were more similar to each other than DELOS groundwater communities (Figure 6) and there is very little ASV overlap between DELOS fractures, the DELOS surface, and nearby Grimsel Test Site (Figure 7). Furthermore, there is only a small tendency for nearby fractures to have more similar microbial communities (Sørensen Distance-Decay R2 = 0.2) with four of the five DELOS time series fractures exhibiting a greater similarity in community composition with GTS samples than the surrounding surface waters (Figure 8).
Figure 6
Figure 7
Figure 8
4 Discussion
The BedrettoLab Deep Life Observatory (DELOS) provides a unique opportunity to study how hydrogeology affects subsurface microbiology. Unlike the majority of underground labs which make use of vertically excavated mines or deep boreholes (e.g. Hallbeck and Pedersen, ; Ino et al., ; Osburn et al., ), DELOS captures a 1, 600 meter depth gradient via a horizontally contiguous 5-km long tunnel through 2 main lithologies (gneiss-schist and granite) in the Alps. As with previous subsurface studies (e.g. Hallbeck and Pedersen, ; Ino et al., ; Osburn et al., ), the different geologies (gneiss-schist from TM-0 to TM-1138 and granite from TM-1138 to TM-5218) host hydrochemically and microbially distinct waters; however, within the same lithology, there are also high amounts of hydrochemical and biological variability correlated to depth, recharge, electrical conductivity and flow rate (Figures 2, 4, 6). Calcium- and sulfate-rich waters appear in the high EC sites (Figure 2) while the back half (>TM-2500) of the tunnel tends to have a larger sodium component. High inflow zones of the tunnel are depleted in sulfate and enriched in bicarbonate and nitrate (Figure 2, Supplementary Data 2). These bicarbonate-type waters are most similar to groundwater found in the granite-hosted Korea Atomic Energy Research Institute Underground Research Tunnel (KAERI URT; Ryu et al., ), Gotthard Rail Base Tunnel (< 15 km from DELOS; Bucher and Stober, ), and Grimsel Test Site (Konno et al., ). On the other hand, water from the Mizunami Underground Research Laboratory (Japan; Iwatsuki et al., ) and Chalk River Laboratory (Canada; Beaton et al., ) uniquely and respectively provide sodium chloride- and sodium bicarbonate-type end members among the granite underground laboratories reported in Figure 9.
Figure 9
Overall, DELOS exhibits a higher degree of hydrochemical diversity than reports of other granite-hosted underground labs (Figure 9) and, within DELOS, this hydrochemical heterogeneity is especially apparent within two granitic inflow zones (TM-1306 and TM-1494) separated by less than 200 m laterally and a difference in overburden of less than 50 m. Water from TM-1494 has a high sulfate concentration and a water isotope signature that is more similar to water encountered in deeper parts of DELOS. On the other hand, the water isotope signature and microbial community of TM-1306 is more similar to shallower parts of the tunnel and the Ticino River (Figure 3, Supplementary Figure 5). As stable water isotopes of oxygen and hydrogen are frequently used to trace groundwater and estimate recharge (Clark and Fritz, ; Blasch and Bryson, ), the hydrochemical and biological variation within the proximal but isotopically distinct fracture waters supports previous observations of groundwater origin and recharge shaping subsurface ecosystems (e.g. Ben Maamar et al., ; Hubalek et al., ; Yan et al., 2021). On the other hand, isotopically similar water from TM-1494 and TM-2848 (Figure 3) exhibit extremely different hydrochemical and microbial community compositions (Figure 5) while isotopically distinct water from TM-2848 and TM-4166 exhibit similar hydrochemistry and microbial community compositions (Figure 4, Supplementary Data 2, 3). The similarities between TM-2848 and TM-4166 appear to be correlated to the relatively high flow rates observed at these sites (Figures 2, 4) which also emerge near zones containing biotite-rich lamprophyres (Figure 2). Previous investigations into the microbial diversity of the metasedimentary Deep Mine Microbial Observatory (USA) have shown that mineralogy plays an important role in shaping subsurface hydrochemistry and microbial communities (Casar et al., , ). A local enrichment of vanadium in the TM-2848 and TM-4166 groundwater and lamprophyres relative to other locations in DELOS (Supplementary Table 1) is suggestive of a similar connection between mineralogy, hydrochemistry and microbial diversity. As most groundwater vanadium studies surround sites with vanadium concentrations several orders of magnitude higher than DELOS (e.g. Xiao et al., 2017; Zhang et al., 2019, 2020a), additional work is needed to discern whether there is a significant relationship between vanadium and microbial diversity in DELOS.
Fracture flow rate, a proxy for hydraulic conductivity, may play an additional and understudied role in shaping subsurface ecosystems. DELOS fractures TM-1306, TM-2848 and TM-4166 occur within the 3 regions of highest inflow and all show an increased proportion of ultra-small bacteria and archaea relative to other sampling sites throughout DELOS. Previous studies have suggested that smaller cells are transported more easily than larger cells due to mechanical filtration (e.g. Gannon et al., ; Abu-Ashour et al., ); however, regions of DELOS with a lower flow rate and, thus a lower permeability and higher potential for filtering (Bense et al., ), do not exhibit an ultra-small microbial enrichment. Mechanical filtering is, therefore, not sufficient to explain the distribution of ultra-small microorganisms in this region and other factors such as groundwater turbulence may play a role. In shallow groundwater, an increase in the detachment of ultra-small bacteria from soils into groundwater during high infiltration periods has been observed (Herrmann et al., , ). If similar processes control the abundance patterns of ultra-small microorganisms in DELOS, the increased turbulence associated within these highly permeable regions may similarly detach populations of ultra-small microorganisms from shallower biofilms and/or provide a medium to more easily carry small particles further.
Seasonal observation of DELOS fractures additionally indicates that groundwater mixing and/or alterations in flow paths can occur in DELOS (Figure 5). Research on the “critical zone” has found that changes in land cover and/or land use can change the biological and chemical properties of subsurface ecosystems and that extreme events such as heavy rainfall and earthquakes dominate the transfer of energy and matter from the surface to shallow (< 100 m below land surface) subsurface environments (Küsel et al., ; Skelton et al., 2019; Hosono et al., ). Seasonal mixing of more oxidizing shallow groundwater with deeper, more reducing fluids have also been observed and shown to introduce new substrates that stimulate the expansion of microorganisms into higher-energy niches (e.g. Magnabosco et al., ; Bochet et al., ; Zhang et al., 2023) and human activities, such as tunnel construction and fluid injection, have been shown to alter groundwater flow paths and water chemistry (e.g. Banwart et al., ; Dong et al., ; Daly et al., ; Zhang et al., 2020b, 2022). Throughout the monitoring period, very minor seasonal changes were observed throughout DELOS (Figure 2); however, a measurable change in electrical conductivity and microbial community composition was observed in TM-1494 following a hydraulic stimulation experiment and rainfall-associated seismic events near the main fault sampled at this location. Hydrochemical measurements taken throughout the tunnel during this period (April–June 2021) suggest that these events may have altered the groundwater flow paths (Figure 2; Arnet, ). In particular, the increase in electrical conductivity observed at TM-1306 and TM-1494 during this period may be stem from an increase in the contribution of higher EC groundwater from the “badboy” fault (TM-1993, EC = 707 ± 64 μS cm−1) to these fractures. Although biological samples from TM-1993 are not available for comparison, a related experimental strategy, known as a “flow test” involving a 10 month injection of fluids into a deep borehole of the metasedimentary Sanford Underground Research Facility, showed that groundwater can be displaced during injection and result in the replacement of pre-existing communities with new microorganisms via advection (Zhang et al., 2022). Together, these observations suggest that the hydraulic stimulation and/or seismic events caused temporary re-organization of TM-1494's groundwater flow paths, recovering to “baseline” conditions within a matter of weeks.
Although further monitoring is needed to fully describe the connectivity of DELOS and the effect hydraulic stimulation experiments have on the hydrochemical and microbial community diversity patterns, correlated changes between monitored fractures provide new insights into the hydro-biogeological subsurface connectivity of a mountain block in a typical alpine watershed. Recently, a study of the 100,000 km2 Death Valley Regional Flow System (DVRFS, USA) identified shared taxa in hydrogeologically connected groundwater separated by over 100 km (Merino et al., ). At a more local scale, a small number of “core” species (97% OTUs) were reported to increase in relative abundance (from 6% up to 47%) as water flowed downstream of the preferential recharge area in the approximately 5 km Hainich Critical Zone Observatory (HCZO, Germany; Yan et al., 2020). Although the observations and relative proportions of core taxa in the DVRFS, HCZO and DELOS are not directly comparable due to differences in filtration volume and sequencing depth, only 2 ASVs were identified in all DELOS sampling locations from the November and December 2020 time points. These ASVs account for 0.9 ± 0.7% of each samples community with the highest relative abundance of these ASVs (2.5%) in TM-1494. Interestingly, despite exhibiting extremely different hydrochemical compositions, DELOS groundwater sampling location TM-1494 and GTS borehole SB80001 share the most ASVs. These ASVs are “rare”, low abundance taxa whose presence in both systems may indicate a similar source of water and microorganisms into these systems. Importantly, neither DELOS TM-1494 or GTS borehole SB80001 have been utilized for hydraulic stimulation experiments; however, these kinds of experiments have been performed in other boreholes within DELOS and GTS (Gischig et al., ; Ma et al., ) and, therefore, the sampling locations may experience similar changes in hydrochemical and microbiological compositions as a result of the injection procedures associated with hydraulic stimulation (e.g. Zhang et al., 2022). Ultimately, the subsurface fractures that exhibit the highest degrees of microbial community overlap tend to correlate with both hydrochemical and hydrogeologic features in DELOS, providing a new insights into the subsurface connectivity of the Riale di Ronco Alpine catchment.
5 Conclusions
Situated in the Alpine catchment of Riale di Ronco, Switzerland, the BedrettoLab Deep Life Observatory (DELOS) provides a new window into the continental subsurface biosphere by capturing an approximately 1, 500 m depth gradient over 5.2 km, largely within the same host rock. Despite the proximity and geologic similarity of subsurface sampling locations within DELOS, a wide range of hydrochemical and microbial community compositions are observed. Despite this variability, broad ecotypes can still be defined. These ecotypes are largely differentiated by the concentration of sulfate, bicarbonate and ultra-small bacteria and archaea. Additional biogeochemical monitoring of DELOS over a one year period further demonstrates the utility of coupling microbial diversity and hydrochemical studies to help inform hydraulic connectivity and groundwater flow. In particular, hydrochemical responses in TM-1306 and TM-1494 following seismicity suggest that groundwater from these locations experience mixing with warmer, sulfate-rich water but recover to their baseline hydrochemical and biological1 compositions within three weeks (Figure 2). Taken together, the spatial and temporal heterogeneity of DELOS provide a first look at the hydrochemical and biological diversity within DELOS and provide new insights in the groundwater history, local mineralogy, and hydraulic conductivity of a subsurface Alpine environment.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, PRJNA1181539. The data and code is available via https://github.com/GeobiologyLab/DELOS-2021-time-series.
Author contributions
AA: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing. MA: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing. NG: Formal analysis, Investigation, Methodology, Writing – review & editing. AC: Formal Analysis, Investigation, Visualization, Writing – review & editing. PR: Formal analysis, Investigation, Writing – review & editing. HT: Formal analysis, Writing – review & editing. QW: Investigation, Writing – review & editing. EZ: Investigation, Writing – review & editing. MH: Investigation, Writing – review & editing. BB: Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Writing – review & editing. CM: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Funding for this project was provided by ETH start-up funding to CM and the Swiss National Science Foundation (Project 215626). Funding for NGD was provided by the Werner Siemens-Stiftung (project 2019-FE-313).
Acknowledgments
We would like to thank Fanny Leuenberger, Simon Loew, Niklaus Kressig, Negar Haghipour, Corey Archer and Thomas Blattmann for assistance with the hydrochemical analyses presented in this study, Erminio Molteni, Bekir Yuece, Lena Bakker, Jasmine Berg and the Bedretto Team for assistance in sampling surface and subsurface locations, Markus Rast for the helpful discussion and insights surrounding the geology of DELOS, and Ingo Blechschmidt and Raphael Schneeberger for access to the Grimsel Test Site.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1522714/full#supplementary-material
Supplementary Data 1.Catalog of DELOS sampling sites.
Supplementary Data 2.Environmental and sequencing metadata.
Supplementary Data 3.16S rRNA gene Amplicon Sequencing Variant and Annotation Table, available on https://github.com/GeobiologyLab/DELOS-2021-time-series.
Supplementary Table 1.A selection of DELOS granitic groundwater trace metal concentrations.
Supplementary Figure 1.The standard curve produced for the qPCR analysis of this study.
Supplementary Figure 2.The melt curve of the qPCR analysis of this study.
Supplementary Figure 3.DELOS 16S rRNA gene amplicon copies per mL.
Supplementary Figure 4.Analysis of Similarities (ANOSIM) of DELOS groundwater 1-year time series samples.
Supplementary Figure 5.Extended version of the water isotope plot of Figure 3.
Footnotes
1.^Microbial community data for TM-1306 during the hydraulic stimulation was not collected.
References
1
Abu-AshourJ.JoyD. M.LeeH.WhiteleyH. R.ZelinS. (1994). Transport of microorganisms through soil. Water, Air, Soil Pollut. 75, 141–158. 10.1007/BF01100406
2
AmbergR. (1983). Design and construction of the Furka Base Tunnel. Rock Mechanics and Rock Eng. 16, 215–231. 10.1007/BF01042358
3
ArnetM. (2021). Deep alpine Fluids: Origin, Pathways and Dynamic Remobilisation in Response to Hydraulic Stimulations at the Bedretto Underground Laboratory for Geoenergies (BULGG) (Master's thesis). Zürich: ETH Zurich.
4
BagnoudA.ChoureyK.HettichR. L.De BruijnI.AnderssonA. F.LeupinO. X.et al. (2016). Reconstructing a hydrogen-driven microbial metabolic network in Opalinus clay rock. Nat. Commun. 7:12770. 10.1038/ncomms12770
5
BanwartS.TullborgE.-L.PedersenK.GustafssonE.LaaksoharjuM.NilssonA.-C.et al. (1996). Organic carbon oxidation induced by large-scale shallow water intrusion into a vertical fracture zone at the Äspö Hard Rock Laboratory (Sweden). J. Contamin. Hydrol. 21, 115–125. 10.1016/0169-7722(95)00037-2
6
BastinE. S.GreerF. E.MerrittC.MoultonG. (1926). The presence of sulphate reducing bacteria in oil field waters. Science63, 21–24. 10.1126/science.63.1618.21
7
BeatonE.StevensonB. S.King-SharpK. J.StampsB. W.NunnH. S.StuartM. (2016). Local and regional diversity reveals dispersal limitation and drift as drivers for groundwater bacterial communities from a fractured granite formation. Front. Microbiol. 7:222054. 10.3389/fmicb.2016.01933
8
Ben MaamarS.AquilinaL.QuaiserA.PauwelsH.Michon-CoudouelS.Vergnaud-AyraudV.et al. (2015). Groundwater isolation governs chemistry and microbial community structure along hydrologic flowpaths. Front. Microbiol. 6:1457. 10.3389/fmicb.2015.01457
9
BenseV.GleesonT.LovelessS.BourO.ScibekJ. (2013). Fault zone hydrogeology. Earth-Sci. Rev. 127:171–192. 10.1016/j.earscirev.2013.09.008
10
BlaschK. W.BrysonJ. R. (2007). Distinguishing sources of ground water recharge by using δ2H and δ18O. Groundwater45, 294–308. 10.1111/j.1745-6584.2006.00289.x
11
BlechschmidtI.VomvorisS. (2015). “The fourth decade of underground research at the Grimsel Test Site-what we have learned and where we go from here,” in Proc. of WM (Phoenix), 16–19.
12
BochetO.BethencourtL.DufresneA.FarasinJ.PédrotM.LabasqueT.et al. (2020). Iron-oxidizer hotspots formed by intermittent oxic-anoxic fluid mixing in fractured rocks. Nature Geosci. 13, 149–155. 10.1038/s41561-019-0509-1
13
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics30, 2114–2120. 10.1093/bioinformatics/btu170
14
BorgonieG.Linage-AlvarezB.OjoA. O.MundleS. O.FreeseL. B.Van RooyenC.et al. (2015). Eukaryotic opportunists dominate the deep-subsurface biosphere in South Africa. Nat. Commun. 6:8952. 10.1038/ncomms9952
15
BowenG. J.RevenaughJ. (2003). Interpolating the isotopic composition of modern meteoric precipitation. Water Resour. Res. 39:10. 10.1029/2003WR002086
16
BucherK.StoberI. (2010). Fluids in the upper continental crust. Geofluids10, 241–253. 10.1002/9781444394900.ch17
17
CallahanB. J.McMurdieP. J.RosenM. J.HanA. W.JohnsonA. J. A.HolmesS. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods13, 581–583. 10.1038/nmeth.3869
18
CasarC. P.KrugerB. R.FlynnT. M.MastersonA. L.MomperL. M.OsburnM. R. (2020). Mineral-hosted biofilm communities in the continental deep subsurface, Deep Mine Microbial Observatory, SD, USA. Geobiology18, 508–522. 10.1111/gbi.12391
19
CasarC. P.KrugerB. R.OsburnM. R. (2021). Rock-hosted subsurface biofilms: mineral selectivity drives hotspots for intraterrestrial life. Front. Microbiol. 12:658988. 10.3389/fmicb.2021.658988
20
ClarkI. D.FritzP. (2013). Environmental Isotopes in Hydrogeology. Boca Raton, FL: CRC Press.
21
DalyR. A.BortonM. A.WilkinsM. J.HoytD. W.KountzD. J.WolfeR. A.et al. (2016). Microbial metabolisms in a 2.5-km-deep ecosystem created by hydraulic fracturing in shales. Nat. Microbiol. 1, 1–9. 10.1038/nmicrobiol.2016.146
22
DoetschJ.GischigV.KrietschH.VilligerL.AmannF.DutlerN.et al. (2018). “Grimsel isc experiment description,” in Technical Report (Zürich: ETH Zurich).
23
DongY.KumarC. G.ChiaN.KimP.-J.MillerP. A.PriceN. D.et al. (2014). Halomonas sulfidaeris-dominated microbial community inhabits a 1.8 km-deep subsurface Cambrian Sandstone reservoir. Environm. Microbiol. 16, 1695–1708. 10.1111/1462-2920.12325
24
FukudaA.HagiwaraH.IshimuraT.KoudukaM.IokaS.AmanoY.et al. (2010). Geomicrobiological properties of ultra-deep granitic groundwater from the Mizunami Underground Research Laboratory (MIU), central Japan. Microb. Ecol. 60, 214–225. 10.1007/s00248-010-9683-9
25
GannonJ.MingelgrinU.AlexanderM.WagenetR. (1991). Bacterial transport through homogeneous soil. Soil Biol. Biochem. 23, 1155–1160. 10.1016/0038-0717(91)90028-I
26
Ginsburg-KaragitschevaT. L. (1926). “Microbiological research in the sulphurous and salty waters of Apsheron,” in Azerbajdzanskoe Neftjanoe Khozjajstvo.
27
GischigV. S.GiardiniD.AmannF.HertrichM.KrietschH.LoewS.et al. (2020). Hydraulic stimulation and fluid circulation experiments in underground laboratories: stepping up the scale toward engineered geothermal systems. Geomech. Energy Environm. 24:100175. 10.1016/j.gete.2019.100175
28
HafnerS. (1958). Petrographie des südwestlichen Gotthardmassivs zwischen St. Gotthardpass und Nufenenpass (PhD thesis). Zurich: ETH Zurich.
29
HallbeckL.PedersenK. (2008). Characterization of microbial processes in deep aquifers of the Fennoscandian Shield. Appl. Geochem. 23, 1796–1819. 10.1016/j.apgeochem.2008.02.012
30
HalterT.Gholizadeh DoonechalyN.NotzonA.RybachL.HertrichM.GiardiniD. (2024). Exploring the feasibility of energy extraction from the Bedretto Tunnel in Switzerland. Energies17:3669. 10.3390/en17153669
31
HerrmannM.LehmannK.TotscheK. U.KüselK. (2023). Seepage-mediated export of bacteria from soil is taxon-specific and driven by seasonal infiltration regimes. Soil Biol. Biochem. 187:109192. 10.1016/j.soilbio.2023.109192
32
HerrmannM.WegnerC.-E.TaubertM.GeesinkP.LehmannK.YanL.et al. (2019). Predominance of Cand. Patescibacteria in groundwater is caused by their preferential mobilization from soils and flourishing under oligotrophic conditions. Front. Microbiol. 10:1407. 10.3389/fmicb.2019.01407
33
HosonoT.YamadaC.MangaM.WangC.-Y.TanimizuM. (2020). Stable isotopes show that earthquakes enhance permeability and release water from mountains. Nature Commun. 11, 1–9. 10.1038/s41467-020-16604-y
34
HubalekV.WuX.EilerA.BuckM.HeimC.DopsonM.et al. (2016). Connectivity to the surface determines diversity patterns in subsurface aquifers of the Fennoscandian Shield. ISME J. 10, 2447–2458. 10.1038/ismej.2016.36
35
InoK.KonnoU.KoudukaM.HirotaA.TogoY. S.FukudaA.et al. (2016). Deep microbial life in high-quality granitic groundwater from geochemically and geographically distinct underground boreholes. Environm. Microbiol. Reports8, 285–294. 10.1111/1758-2229.12379
36
IwatsukiT.FurueR.MieH.IokaS.MizunoT. (2005). Hydrochemical baseline condition of groundwater at the Mizunami Underground Research Laboratory (MIU). Appl. Geochem. 20, 2283–2302. 10.1016/j.apgeochem.2005.09.002
37
KellerF.SchneiderT. (1982). Geologie und geotechnik. Schweizer Ingenieur und Architekt100, 512–520.
38
KeusenH.GanguinJ.SchulerP.BulettiM. (1989). “Grimsel Test Site: geology,” in Technical Report, Nationale Genossenschaft fuer die Lagerung Radioaktiver Abfaelle (NAGRA). Available at: https://www.osti.gov/etdeweb/biblio/6997735
39
KonnoU.KoudukaM.KomatsuD. D.IshiiK.FukudaA.TsunogaiU.et al. (2013). Novel microbial populations in deep granitic groundwater from Grimsel Test Site, Switzerland. Microb. Ecol. 65, 626–637. 10.1007/s00248-013-0184-5
40
KotelnikovaS.PedersenK. (1998). Distribution and activity of methanogens and homoacetogens in deep granitic aquifers at Äspö Hard Rock Laboratory, Sweden. FEMS Microbiology Ecology26:121–134. 10.1016/S0168-6496(98)00028-2
41
KüselK.TotscheK. U.TrumboreS. E.LehmannR.SteinhäuserC.HerrmannM. (2016). How deep can surface signals be traced in the critical zone? Merging biodiversity with biogeochemistry research in a central German Muschelkalk landscape. Front. Earth Sci. 4:32. 10.3389/feart.2016.00032
42
KyleJ. E.EydalH. S.FerrisF. G.PedersenK. (2008). Viruses in granitic groundwater from 69 to 450 m depth of the Äspö Hard Rock Laboratory, Sweden. ISME J. 2, 571–574. 10.1038/ismej.2008.18
43
LabhartT. (2005). Erläuterungen zum Geologischen Atlas der Schweiz 1: 25000, Val Bedretto, Atlasblatt 68. Bern-Ittigen: Bundesamt für Wasser und Geologie.
44
LipmanC. B. (1928). The discovery of living micro-organisms in ancient rocks. Science68, 272–273. 10.1126/science.68.1760.272.b
45
Lopez-FernandezM.SimoneD.WuX.SolerL.NilssonE.HolmfeldtK.et al. (2018). Metatranscriptomes reveal that all three domains of life are active but are dominated by bacteria in the Fennoscandian crystalline granitic continental deep biosphere. MBio9, 10–1128. 10.1128/mBio.01792-18
46
Lopez-FernandezM.WestmeijerG.TurnerS.BromanE.StåhleM.BertilssonS.et al. (2023). Thiobacillus as a key player for biofilm formation in oligotrophic groundwaters of the Fennoscandian Shield. NPJ Biofilms Microb. 9:41. 10.1038/s41522-023-00408-1
47
LützenkirchenV. H. (2002). Structural Geology and Hydrogeology of Brittle Fault Zones in the Central and Eastern Gotthard Massif, Switzerland (PhD thesis). Zurich: ETH Zurich.
48
MaX.HertrichM.AmannF.BrökerK.Gholizadeh DoonechalyN.GischigV.et al. (2022). Multi-disciplinary characterizations of the BedrettoLab-a new underground geoscience research facility. Solid Earth13, 301–322. 10.5194/se-13-301-2022
49
MagnaboscoC.BiddleJ.CockellC.JungbluthS.TwingK. (2019). “Biogeography, ecology, and evolution of deep life,” in Deep Carbon (Cambridge: Cambridge University Press), 524–555.
50
MagnaboscoC.LinL.-H.DongH.BombergM.GhiorseW.Stan-LotterH.et al. (2018a). The biomass and biodiversity of the continental subsurface. Nat. Geosci. 11, 707–717. 10.1038/s41561-018-0221-6
51
MagnaboscoC.TimmersP. H.LauM. C.BorgonieG.Linage-AlvarezB.KuloyoO.et al. (2018b). Fluctuations in populations of subsurface methane oxidizers in coordination with changes in electron acceptor availability. FEMS Microbiol. Ecol. 94:fiy089. 10.1093/femsec/fiy089
52
MassetO.LoewS. (2010). Hydraulic conductivity distribution in crystalline rocks, derived from inflows to tunnels and galleries in the Central Alps, Switzerland. Hydrogeol. J. 18, 863–891. 10.1007/s10040-009-0569-1
53
McNicholJ.BerubeP. M.BillerS. J.FuhrmanJ. A. (2021). Evaluating and improving small subunit rRNA PCR primer coverage for bacteria, archaea, and eukaryotes using metagenomes from global ocean surveys. MSystems6, 10–1128. 10.1128/mSystems.00565-21
54
MerinoN.JacksonT. R.CampbellJ. H.KerstingA. B.SackettJ.FisherJ. C.et al. (2022). Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow. Sci. Total Environm. 842:156768. 10.1016/j.scitotenv.2022.156768
55
MomperL.CasarC. P.OsburnM. R. (2023). A metagenomic view of novel microbial and metabolic diversity found within the deep terrestrial biosphere at DeMMO: a microbial observatory in South Dakota, USA. Environm. Microbiol. 25, 3719–3737. 10.1111/1462-2920.16543
56
OfterdingerU. S. (2001). Ground Water Flow Systems in the Rotondo granite, Central Alps (Switzerland) (PhD thesis). Zurich: ETH Zurich.
57
OnstottT. C. (2016). Deep Life: The Hunt for the Hidden Biology of Earth, Mars, and Beyond.
58
OnstottT. C.MagnaboscoC.AubreyA.BurtonA.DworkinJ.ElsilaJ.et al. (2014). Does aspartic acid racemization constrain the depth limit of the subsurface biosphere?Geobiology12, 1–19. 10.1111/gbi.12069
59
OsburnM. R.KrugerB.MastersonA. L.CasarC. P.AmendJ. P. (2019). Establishment of the Deep Mine Microbial Observatory (DeMMO), South Dakota, USA, a geochemically stable portal into the deep subsurface. Front. Earth Sci. 7:196. 10.3389/feart.2019.00196
60
PedersenK. (1997). Microbial life in deep granitic rock. FEMS Microbiol. Rev. 20, 399–414. 10.1016/S0168-6445(97)00022-3
61
PoulainS.SergeantC.SimonoffM.Le MarrecC.AltmannS. (2008). Microbial investigations in Opalinus Clay, an argillaceous formation under evaluation as a potential host rock for a radioactive waste repository. Geomicrobiol. J. 25, 240–249. 10.1080/01490450802153314
62
QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. 10.1093/nar/gks1219
63
QuintoF.BlechschmidtI.Garcia PerezC.GeckeisH.GeyerF.GolserR.et al. (2017). Multiactinide analysis with accelerator mass spectrometry for ultratrace determination in small samples: application to an in situ radionuclide tracer test within the colloid formation and migration experiment at the Grimsel Test Site (Switzerland). Analyt. Chem. 89, 7182–7189. 10.1021/acs.analchem.7b01359
64
RastM.GalliA.RuhJ. B.GuillongM.MadonnaC. (2022). Geology along the Bedretto Tunnel: kinematic and geochronological constraints on the evolution of the Gotthard Massif (Central Alps). Swiss J. Geosci. 115:8. 10.1186/s00015-022-00409-w
65
RyuJ.-H.KohY.-K.ParkS.-W.KimG.-Y.ChoiJ.-W. (2012). Geochemical characterization of deep groundwater in KURT using geochemical modeling. J. Environm. Eng. 138, 351–359. 10.1061/(ASCE)EE.1943-7870.0000492
66
SchneebergerR.MäderU. K.WaberH. N. (2017). Hydrochemical and isotopic (δ2H, δ18O, 3H) characterization of fracture water in crystalline rock (Grimsel, Switzerland). Procedia Earth Planet. Sci. 17, 738–741. 10.1016/j.proeps.2016.12.187
67
SkeltonA.Liljedahl-ClaessonL.WästebyN.AndrénM.StockmannG.SturkellE.et al. (2019). Hydrochemical changes before and after earthquakes based on long-term measurements of multiple parameters at two sites in northern Iceland–a review. J. Geophys. Res. 124, 2702–2720. 10.1029/2018JB016757
68
Stroes-GascoyneS.SchippersA.SchwynB.PoulainS.SergeantC.SimonoffM.et al. (2007). Microbial community analysis of Opalinus clay drill core samples from the Mont Terri underground research laboratory, Switzerland. Geomicrobiol. J. 24, 1–17. 10.1080/01490450601134275
69
WaltersW.HydeE. R.Berg-LyonsD.AckermannG.HumphreyG.ParadaA.et al. (2016). Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. MSystems1, 10–1128. 10.1128/mSystems.00009-15
70
WestmeijerG.MehrshadM.TurnerS.AlakangasL.SachpazidouV.BunseC.et al. (2022). Connectivity of Fennoscandian Shield terrestrial deep biosphere microbiomes with surface communities. Commun. Biol. 5:37. 10.1038/s42003-021-02980-8
71
WuX.PedersenK.EdlundJ.ErikssonL.ÅströmM.AnderssonA. F.et al. (2017). Potential for hydrogen-oxidizing chemolithoautotrophic and diazotrophic populations to initiate biofilm formation in oligotrophic, deep terrestrial subsurface waters. Microbiome5, 1–13. 10.1186/s40168-017-0253-y
72
XiaoX.-Y.WangM.-W.ZhuH.-W.GuoZ.-H.HanX.-Q.ZengP. (2017). Response of soil microbial activities and microbial community structure to vanadium stress. Ecotoxicol. Environm. Safety142, 200–206. 10.1016/j.ecoenv.2017.03.047
73
YanL.HermansS. M.TotscheK. U.LehmannR.HerrmannM.KüselK. (2021). Groundwater bacterial communities evolve over time in response to recharge. Water Res. 201:117290. 10.1016/j.watres.2021.117290
74
YanL.HerrmannM.KampeB.LehmannR.TotscheK. U.KüselK. (2020). Environmental selection shapes the formation of near-surface groundwater microbiomes. Water Res. 170:115341. 10.1016/j.watres.2019.115341
75
YilmazP.ParfreyL. W.YarzaP.GerkenJ.PruesseE.QuastC.et al. (2014). The SILVA and “all-species living tree project (LTP)” taxonomic frameworks. Nucleic Acids Res. 42, D643–D648. 10.1093/nar/gkt1209
76
ZhangB.WangS.DiaoM.FuJ.XieM.ShiJ.et al. (2019). Microbial community responses to vanadium distributions in mining geological environments and bioremediation assessment. J. Geophys. Res.: Biogeosci. 124, 601–615. 10.1029/2018JG004670
77
ZhangH.ZhangB.WangS.ChenJ.JiangB.XingY. (2020a). Spatiotemporal vanadium distribution in soils with microbial community dynamics at vanadium smelting site. Environm. Pollut. 265:114782. 10.1016/j.envpol.2020.114782
78
ZhangX.GaoX.LiC.LuoW.WangY.LuoX. (2023). Seasonal dynamics and interaction of shallow groundwater geochemical properties and microbial community patterns. Chem. Geol. 638:121703. 10.1016/j.chemgeo.2023.121703
79
ZhangY.DekasA. E.HawkinsA. J.ParadaA. E.GorbatenkoO.LiK.et al. (2020b). Microbial community composition in deep-subsurface reservoir fluids reveals natural interwell connectivity. Water Resour. Res. 56:e2019WR025916. 10.1029/2019WR025916
80
ZhangY.HorneR. N.HawkinsA. J.PrimoJ. C.GorbatenkoO.DekasA. E. (2022). Geological activity shapes the microbiome in deep-subsurface aquifers by advection. Proc. Nat. Acad. Sci. 119:e2113985119. 10.1073/pnas.2113985119
Summary
Keywords
subsurface microbiology, underground laboratory, groundwater, microbial ecology, deep life observatory
Citation
Acciardo AS, Arnet M, Gholizadeh Doonechaly N, Ceccato A, Rodriguez P, Tran HNH, Wenning Q, Zimmerman E, Hertrich M, Brixel B and Magnabosco C (2025) Spatial and temporal groundwater biogeochemical variability help inform subsurface connectivity within a high-altitude Alpine catchment (Riale di Ronco, Switzerland). Front. Microbiol. 16:1522714. doi: 10.3389/fmicb.2025.1522714
Received
04 November 2024
Accepted
28 January 2025
Published
25 February 2025
Volume
16 - 2025
Edited by
William J. Brazelton, The University of Utah, United States
Reviewed by
Magdalena R. Osburn, Northwestern University, United States
George Westmeijer, Umeå University, Sweden
Daniel Gonzalez-Duque, The Ohio State University, United States
Updates
Copyright
© 2025 Acciardo, Arnet, Gholizadeh Doonechaly, Ceccato, Rodriguez, Tran, Wenning, Zimmerman, Hertrich, Brixel and Magnabosco.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Cara Magnabosco cara.magnabosco@eaps.ethz.ch
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.