in situ observation of helium and argon release during
fluid-pressure- triggered rock deformation

Clément Roques, Ulrich W Weber, Bernard Brixel, Hannes Krietsch, Nathan
Dutler, Matthias Brennwald, Linus Villiger, Joseph Doetsch, Mohammadreza
Jalali, Valentin Gischig, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Clément Roques, Ulrich W Weber, Bernard Brixel, Hannes Krietsch, Nathan Dutler, et al. . in situ observation of helium and argon release during fluid-pressure- triggered rock deformation. Scientific Reports, 2020, 10 (1), pp.Art. n°6949. 10.1038/s41598-020-63458-x . insu-02563199

HAL Id: insu-02563199

https://insu.hal.science/insu-02563199

Submitted on 5 May 2020

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# In situ observation of helium and argon release during fluid-pressure-triggered rock deformation

Clément Roques<fnref n="1" />,²✉, Ulrich W. Weber³,⁴, Bernard Brixel<fnref n="1" />, Hannes Krietsch<fnref n="1" />, Nathan Dutler⁵, Matthias S. Brennwald³, Linus Villiger<fnref n="1" />, Joseph Doetsch<fnref n="1" />, Mohammadreza Jalali⁶, Valentin Gischig<fnref n="1" />,⁷, Florian Amann⁶, Benoît Valley⁵, Maria Klepikova⁸ & Rolf Kipfer³,⁹,¹⁰

Temporal changes in groundwater chemistry can reveal information about the evolution of flow path connectivity during crustal deformation. Here, we report transient helium and argon concentration anomalies monitored during a series of hydraulic reservoir stimulation experiments measured with an in situ gas equilibrium membrane inlet mass spectrometer. Geodetic and seismic analyses revealed that the applied stimulation treatments led to the formation of new fractures (hydraulic fracturing) and the reactivation of natural fractures (hydraulic shearing), both of which remobilized (He, Ar)-enriched fluids trapped in the rock mass. Our results demonstrate that integrating geochemical information with geodetic and seismic data provides critical insights to understanding dynamic changes in fracture network connectivity during reservoir stimulation. The results of this study also shed light on the linkages between fluid migration, rock deformation and seismicity at the decimeter scale.

Changes in the crustal stress state caused by natural or human-induced subsurface fluid overpressures can lead to brittle rock mass damage, from the formation of grain-scale microcracks to the failure of kilometric-scale faults<fnref n="1" />. Characterizing the timing and spatial extent of crustal deformation is critical for both industrial applications, including geothermal, oil and gas production, and improving our mechanistic understanding of earthquakes. From this perspective, seismic and geodetic monitoring systems generally provide the vast majority of data collected. Several authors have also proposed that geochemical anomalies may be used as proxies for rock deformation²⁻⁵.

Fluids in the Earth's crust display highly variable geochemical traits due to the spatial and temporal variability in fluid recharge composition, fluid-rock interactions and residence times. Previous studies have shown extremely variable distributions of residence times, ranging from decades⁶ to several hundred or even millions of years for fluids stored in low-permeability rocks⁷,⁸. In this context, the fluid composition may evolve from diluted waters with modern signatures in recharge areas to saline fluids in the deeper crust, where the fluid composition tends to equilibrate with the host rock mineralogy through dissolution/precipitation processes. By analyzing specific dissolved chemical tracers, one can reconstruct the origin of fluids and gain insights into their recharge, percolation and storage conditions⁹,¹⁰.

From this perspective, noble gases such as radon (Rn), helium (He) and argon (Ar) have received particular attention in recent decades because of their widespread occurrence in the subsurface and their low chemical reactivity, proving to be ideal tracers to track the origin of fluids, analyze fluid-rock interactions and determine residence times⁷,⁸,¹¹⁻¹⁴. Concentrations of dissolved noble gases in subsurface fluids and their isotopic compositions are driven by two main mechanisms: 1) water/air partitioning during recharge, which is sensitive to atmospheric

<footnote n="1">¹ETH Zürich, Department of Earth Sciences, Sonneggstrasse 5, 8092, Zürich, Switzerland. ²University Rennes 1, Géosciences Rennes, UMR 6118, Av. du Général Leclerc, 35042, Rennes, France. ³Eawag - Swiss Federal Institute for Aquatic Science and Technology, Department of Water Resources and Drinking Water, Ueberlandstrasse 133, 8600, Dübendorf, Switzerland. ⁴University of Oslo, Department of Geosciences, Sem Sælands vei 1, 0371, Oslo, Norway. ⁵University of Neuchâtel, Center for Hydrogeology and Geothermics, Rue Emile-Argand 11, 2000, Neuchâtel, Switzerland. ⁶RWTH Aachen, Department of Engineering Geology and Hydrogeology, Lochnerstrasse 4-20, 52064, Aachen, Germany. ⁷CSD INGENIEURE AG, Hessstrasse 27D, 3097, Liebefeld, Switzerland. ⁸University of Lausanne, Applied and Environmental Geophysics group, Institute of Earth Sciences, Lausanne, 1015, Switzerland. ⁹ETH Zürich, Department of Environmental System Science, Universtaetstrasse 16, 8092, Zürich, Switzerland. ¹⁰ETH Zürich, Department of Earth Sciences, Institute of Geochemistry and Petrology, Sonneggstrasse 5, 8092, Zürich, Switzerland. ✉e-mail: clement.roques@univ-rennes1.fr</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Overview and location of the experimental test volume for the *In situ* Stimulation and Circulation (ISC) project: (a) geological map of Switzerland (modified from the website https://map.geo.admin.ch/, Federal Office of Topography), (b) tunnel network and location of the investigated volume, (c) geological model displaying the main shear zones (modified after Krietsch *et al.*, 2018²³), along with the main boreholes used for injection (INJ1 and INJ2) and monitoring (FBS), the locations of stimulation (HF#) and strain monitoring (Fiber-Bragg-Grating, FBG) intervals, and the location of the *in situ* gas equilibrium membrane inlet mass spectrometer (GE-MIMS[24]) used for monitoring dissolved noble gas concentrations. For a more complete description of the geology of the investigated volume and the monitoring system deployed, readers are referred to Krietsch *et al.*, 2018 and Amann *et al.*, 2018²²,²³.</figcaption></figure>

pressure, temperature and salinity, and 2) radiogenic production in the crust through radioactive alpha-decay of heavy elements (mainly uranium, thorium or potassium) combined with the accumulation of enriched fluids from deeper reservoirs (e.g., the upper mantle)⁹.

Chemical anomalies created by earthquakes have been explained by the mixing of fluids from different sources (with different chemical compositions) inhibited by changes in flow pathway geometry and hydraulic properties during stress build-up and failure¹⁵,¹⁶. Laboratory studies have also observed the emergence of concentration anomalies in noble gases during rock deformation experiments¹⁷⁻²⁰. The authors argued that the progressive increase in microcrack surfaces during increasing stress allows the remobilization of accumulated radiogenic gases and even prior macroscopic failure. While those previous studies suggested that changes in dissolved noble gas concentration could be used as a proxy to track *in situ* stress modifications, field evidence and a description of the mechanisms responsible for fluid remobilization and mixing at the scale of discrete fracture systems are still missing to our knowledge. Our aim in this study is to bridge this knowledge gap through a joint analysis of hydraulic, seismic, deformation and geochemical data gathered during an *in situ* stimulation experiment.

## Experiment description

The study was carried out as part of a series of six controlled hydraulic fracturing (HF) experiments²¹ that took place in May 2017 at the Grimsel Test Site, Switzerland (*In situ* Stimulation and Circulation project²², www.grimsel.com, Fig. 1). The site consists of a network of tunnels, operated by the Swiss National Cooperative for the Disposal of Radioactive Waste (NAGRA), and these tunnels are located up to 480 m below ground surface in the central Aar granite and granodiorite formations²³. The stimulated rock volume comprises a relatively intact granodiorite intersected by six subvertical shear zones that can be grouped into 2 sets: ductile shear zones S1 (red in Fig. 1c) and brittle-ductile shear zones S3 (green in Fig. 1c), which form the main transmissive fractured zones²³.

In this study, we refer to a 'fracture' as a single discontinuity that can be identified on optical or acoustic logs (resolution of ~1 mm). A 'shear zone' is described as a zone characterized by strong strain concentration that may contain several fractures visible in the optical or acoustic logs. Finally, 'intact rock' is understood to be a section of the rock mass with the absence of fractures detectable from the optical and acoustic logs.

From May 16th to 19th, six high-pressure fluid injections (resulting pressure below 10 MPa during the fracture propagation phase) were performed in two injection boreholes (INJ1-2 in Fig. 1c, labeled HF#) at different depth intervals by isolating 1 m intact rock with double packer systems. A dense network of boreholes was equipped with monitoring systems aimed at capturing hydraulic, strain and seismic responses during stimulation (see Methods). As a complement, a portable gas equilibrium membrane inlet mass spectrometer (GE-MIMS²⁴) was installed at a tunnel inflow, where groundwater naturally seeps from an extensional fracture located between the
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Temporal evolution of (a) injection pressure drawn as a continuous black line together with the magnitude of seismic events, which are represented by circle symbols. The diameter of the circles are scaled with the magnitude and the colors differentiate the responses from the different HF experiments; (b) normalized strain $\epsilon^* = \epsilon/\text{std}(\epsilon)$ [-] for two selected strain sensors located in both S1 and S3 shear zones. Strain has been normalized by its standard deviation to better appreciate the shapes of both time series. Positive values indicate compressional regimes, while negative values indicate extensional regimes; (c) fracture discharge rates monitored at the GE-MIMS station. The horizontal dashed line represents the average discharge rate prior to stimulation; (d) relative changes in He/N₂ from background [%] and (e) relative changes in Ar/N₂ from background [%]. The background average ratio in He/N₂ and Ar/N₂ are shown as blue squares. Background error bars are 2σ and represent the typical analytical error. Time series of raw data are shown as blue triangles, and low-pass filtered ratios are represented by magenta lines. The different stimulation experiments are labeled with HF# at the top of the figure. Colored triangles in subpanel c identify the peaks in He/N₂ concentrations following each HF# experiment, which are further analyzed in Fig. 3.</figcaption></figure>

two targeted S3 shear zones. The GE-MIMS allowed monitoring of near-real-time He, N₂, and Ar concentrations dissolved in the pore water with a sampling interval of 8 minutes (see Methods). We anticipated potential effects of the low flowrates involved at the seeping fracture that could lead to incomplete gas/water equilibration at the extraction membrane of the GE-MIMS, by normalizing the He and Ar concentrations by the concentrations of dissolved N₂ (only of atmospheric origin in the case of this experiment). This ensures that the observed He/N₂ and Ar/N₂ ratios only reflect the variations of dissolved He and Ar in the fluid (see Methods).

## Results

**Strain, seismic and hydraulic responses induced by high-pressure fluid injection.** The evolution of strain monitored in the S1 and S3 shear zones during the different high-pressure injections showed that both extension and compression deformation were involved. (Figure 2b, positive values for the compressional regime and negative values for the extensional regime). Deformation patterns show variability across the different HF experiments with amplitudes ranging from -386.3 µε extension to +28.5 µε compression (Table 1SM in Supplementary Material). As a general interpretation of the main processes involved, we argue that while the extensional regime is dominant at the locus of injection due to local increasing pore pressure, the nearby structures experience a compressive response to balance the stress perturbation. As the injection progresses,
<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Ar/N₂ as a function of He/N₂. The mixing line (short dashed line) from atmospheric (black reversed triangle) to (He, Ar)-enriched fluid end-members (olive star) was collected in 1982 during the realization of exploration boreholes and reported in Frick et al.<sup>30,31</sup>. The confidence area for binary mixing is drawn as a gray triangle accounting for the variability in the Ar and N₂ concentrations measured in the samples collected in 1982. The background average ratio and standard deviations are shown as blue squares and error bars. (b) is a zoom of the concentration range of interest (corresponding to the dashed rectangle in panel a) with peak concentrations presented with different symbols and colors for the different HF# experiments (refer to Fig. 2d for the locations of peak concentrations identified in the time series of He/N₂). Concentrations measured after the stimulation experiment are shown as a wine-colored circle. Error bars are 2σ.</figcaption></figure>

the stress perturbation diffuses away from the stimulation interval, enabling fracture initiation and reactivation. After a peak in deformation is monitored almost simultaneously with the maximum injection pressure, the strain decreases back to values close to the initial background. This recovery is mostly controlled by the combined effect of fluid pressure diffusion within the fracture network and elastic relaxation of the rock. Although the strain signal appears to be mainly dominated by reversible processes, the residual strain observed in the time series indicates that irreversible deformation due to imperfect normal fracture closing is involved<sup>21</sup>.

A total of 774 seismic events were located during the whole experiment, with magnitudes ranging from −6.2 to −3.1 (Fig. 2a). Seismic events strongly correlate with injection activities, and even though similar fluid volumes were injected during all experiments following similar injection protocols, strong variability in seismic characteristics was observed across a single HF experiment (Table 1SM in Supplementary Material) - i.e., located seismic events, maximum induced magnitudes, seismically activated areas (determined with planar seismic cluster geometries) and Gutenberg-Richter *b*-values, as well as seismogenic indices, $\Sigma$ (determined for experiments exhibiting more than 50 induced events). The Gutenberg-Richter *b*-values and the seismogenic indices, $\Sigma$, are of particular interest when comparing the intensity of seismic responses across different tests<sup>25</sup>. Note that a low *b*-value (i.e., comparably larger magnitude events with respect to small magnitude events) and a high seismogenic index (i.e., comparably increased seismogenic reaction to a unit volume of fluid) represents a comparable, more intensive seismic response<sup>26</sup> (Fig. 1SM in Supplementary Material).

High-pressure injection activities promoted increases in discharge rates monitored at the GE-MIMS station during HF2, 5, 6 and 8 (Fig. 2c). These fast reactions at the tunnel outlet reveal efficient connectivity and high hydraulic diffusivities (on the order of 0.1–1 m/s) of the stimulated fracture network.

The main differences in seismic, deformation and hydraulic responses among the HF experiments allow for the definition of two main clusters, which are primarily constrained by their location with respect to the S3 shear zones: (I) **HF3, 5 and 8** were all carried out south of the S3 shear zones. All of these experiments showed high injection pressures, fast pressure decay after injections and relatively low fluid recovery. This finding suggests that the stimulated fractures progressively connected to the preexisting transmissive fracture system, allowing the injected fluid to propagate efficiently within the network<sup>21</sup>. HF3 and HF8 induced 70 and 183 micro-seismic events, respectively, for which deformation patterns were mainly dominated by the extensional regime along the S3 shear zones. While HF3 showed a dispersed seismic cloud over a large rock volume, seismic events located during HF8 were organized along a well-identified fracture plane (estimated area of 235 m², see Table 1SM in Supplementary Material) propagating toward the S3 shear zone. The seismic response of experiment HF3 (b = 1.55, $\Sigma$: −4.8) increased compared to experiment HF8 (b = 2.66, $\Sigma$: −9.0), and the largest magnitude event of the experiment series was monitored during HF3. HF5 showed the lowest number of seismic events (11), with a smaller seismically activated area (<10 m²) and a compressional regime along S3<sup>21,26</sup>. (II) **HF2 and 6** were performed within the S1 shear zones – i.e., north of the S3 shear zones. Pressure diffused at a slower rate within the reservoir, and fluid recovery was higher after the shut-in phase, suggesting lower connectivity to the preexisting transmissive network. Extensional deformation along S1 was dominant during HF2, while HF6 was mainly linked to compressional deformation of S1 and the extension of S3<sup>21</sup>. HF6 did not induce a significant
number of microseismic events, which was most likely because the experiment was conducted in an interval that contained preexisting S1 structures that were hydraulically stimulated in a previous experiment series²⁶. HF2 promoted the highest number of induced seismic events, an activated surface area of approximately 95 m² and the largest seismic response of the HF experiment series²⁶ (b = 1.35, Σ: -4.0, also see Table 1SM and Fig. 1SM in Supplementary Material).

**Geochemical tracing of involved fluid origins before stimulation.** Pore water sampled at the monitoring station shows low mineralization with an electrical conductivity of 300 µS/cm, a pH of 9 and a temperature of 12 °C (Table 2SM in Supplementary Material). N₂ and Ar concentrations of 1.4·10⁻² ccSTP/g and 3.7·10⁻⁴ ccSTP/g, respectively, indicate that both elements are oversaturated by approximately 25% with respect to atmospheric-equilibrated water. This excess of N₂ and Ar is typically assigned to an “excess-air” phenomenon, which implies oversaturation in gas species due to the partial dissolution of entrapped air bubbles during recharge²⁷,²⁸. Ar shows 2.5% more oversaturation than N₂, while He is also significantly enriched with respect to atmospheric-equilibrated water of approximately 800% (3.4·10⁻⁷ ccSTP/g, Table 2SM in Supplementary Material). This oversaturation in He, combined with the slight enrichment in Ar over N₂ of ~2.5%, shows the accumulation of radiogenic-produced noble gases, which is confirmed by a low ³He/⁴He ratio of the excess He (taking a value of 2.2·10⁻⁷, which is 6.1 times lower than the atmospheric ratio²⁹ of 1.384·10⁻⁶; please refer to Table 2SM in the Supplementary Material and Methods).

Even larger excesses in He, Ar and N₂ (i.e., concentrations exceeding the expected atmospheric equilibrium concentrations) were also determined in pore waters collected during early sampling campaigns on exploration boreholes in 1982³⁰,³¹ (Fig. 3a and Table 2SM in Supplementary Material). These differences indicate that the initial pore water sampled in 1982 could have slightly different meteoric origins with higher soil temperatures at recharge, different residence times controlled by local geological heterogeneity, or that the present-day concentrations may be diluted by modern fluids injected during recent activities performed in the tunnel, i.e., extensive drilling, hydraulic testing, and in situ experiments. Nevertheless, the relationship between He/N₂ and Ar/N₂ (Fig. 3a) is remarkably aligned along a mixing line between the non-atmospheric fluid enriched in He and Ar that was collected in 1982 and modern air-equilibrated meteoric water. Note that the background concentration (Fig. 3a) plots in the mixing triangle (in gray on Fig. 3) defined by the variability in Ar/N₂ ratios measured in the pore water sampled in 1982. This confirms that binary mixing between the two geochemical end-members explains the fluid composition seep- ing from the fracture where the GE-MIMS was installed. It also allows us to estimate that the (He, Ar)-enriched end-member contributes approximately 25% to the background fluid composition.

**Geochemical anomalies involved during reservoir stimulation.** The first two injections, HF1 and HF3, performed on May 16th (Fig. 2) caused a significant increase in He/N₂ ratios of up to 17% relative to the background values, while Ar/N₂ ratios remained close to background values (Fig. 2d,e). After this sudden increase, the He/N₂ ratios tend to continuously decrease. This trend was interrupted by enhanced ratios in response to the following stimulation (HF2, 5, 6 and 8). Experiments HF3 and 8, which were carried out south of the transmissive shear zones and the closest to the GE-MIMS monitoring station (24 and 17.4 m geometrical distance from locus of injection, respectively, Table 1SM in Supplementary Material), were associated with an abrupt increase in the He/N₂ anomalies. HF5 fostered a much smaller anomaly, while HF2 and HF6 were hardly detectable against the given experimental uncertainty of the dissolved gas analysis. In response to HF8, an ~2% increase in the Ar/N₂ ratio was observed.

Enhanced He/N₂ ratios correlate with the occurrence of microseismic events (Fig. 2a) and increase in discharge rates (Fig. 2c), indicating that He-enriched fluids are remobilized within the fracture network during the injection cycles. The majority of normalized peak concentrations lines up subparallel to a binary mixing line of the two respective end-members and within the confidence area (Fig. 3b), suggesting that no other type of fluid is involved. Mixing fractions of the He/N₂ and Ar/N₂-enriched fluid ranged from 20 to 45%, whereby maximum values were linked to stimulation experiments HF3 (33%), HF5 (33%) and HF8 (45%), while HF2 and HF6 showed lower values of 28% and 21%, respectively. It is remarkable that these two sets of experiments identified based on the differences in geochemical anomalies match the two clusters that have been previously identified as having similar hydromechanical responses. Experiments HF3, 5 and 8, which were performed south of the S3 shear zone, have enabled the development of new fractures that effectively connected the preexisting transmissive fracture network and were followed by higher He/N₂ anomalies. However, experiments HF2 and HF6, which were undertaken north of the S3 shear zone, induced poor fracture connection with the transmissive S3 shear zone and were followed by lower He/N₂ anomalies.

The day after the last experiment, He/N₂ and Ar/N₂ reached a lower level than before the experiments, with a contribution of approximately 11.5% from the (He, Ar)-enriched fluid. Figure 3b also shows that the concentrations plot outside the confidence area for binary mixing. This difference may result from the modifications in the fracture network connectivity due to the stimulation, which tends to increase the contribution of a different source of water that was not previously identified (i.e., meteoric water recharged under different temperature and pressure conditions).

## Discussion

**He and Ar concentration anomalies provide insights into new flow pathways formed in response to reservoir stimulation.** Our assessment demonstrates the ability of GE-MIMS technology to monitor fluid remobilization in response to *in situ* rock deformation. The temporal scales of the processes involved not only require a high-frequency monitoring system²⁴,³² but also a multi-tracer approach, i.e., He, Ar and N₂ in the case of the GTS. Complementing the findings of previous experiments performed in the laboratory¹⁷⁻²⁰, as
well as those advanced from large-scale observations in the context of earthquakes⁵, we interpret these anomalies as the release of radiogenic He (and Ar) stored in the rock mass that could be transported within the newly created and pre-existing fracture system. While the identification of the source of the mobilized fluids is still speculative, two complementing hypotheses on its origin are proposed: (1) radiogenic He (and Ar) accumulated in the rock matrix is released by the formation of new fracture areas; and/or (2) the shearing of pre-existing fractures induces the remobilization of stagnant fluids, with higher residence times enriched in radiogenic He and Ar, trapped in lower transmissivity zones of the fracture network. In both scenarios, the improved connectivity resulting from the stimulation, combined with the high pore pressure gradients involved, allows the transport and dispersion of those fluids with their respective noble gas signature away from their original location.

The magnitude and temporal scaling of the produced concentration anomalies depend not only on concentration contrasts between the involved fluids (geochemical end members) but also on flow connectivity and transport properties of the fracture system through which the fluids move. The time lag between the high-pressure injection and the peak in noble gas concentrations may provide estimates of the average fluid velocities involved if the original location of the activated fluids is known, i.e., assuming that most of the (He, Ar)-enriched fluid originates from the vicinity of the injection point where the fractures are created or reactivated. In the case of this experiment, advective velocities of $10^{-3}$ to $10^{-2}$ m/s can explain the observed time lags. The rate at which concentration anomalies dissipate within the fracture network may provide information on how much the fluid is dispersed within the reservoir, with the dispersion length scaling with the square root of the residence time³³. We observed a He/N₂ restitution rate during HF3 that suggested higher dispersion than for the other experiments, with an anomaly that persisted during the whole week of the experiment. HF3 also showed a dispersed seismic response, suggesting the activation of multiple fractures in the surrounding volume, enhancing the dispersion of (He, Ar)-enriched fluids. In contrast, the abrupt and sharp He/N₂ anomaly linked to HF8 relaxes to the corresponding background values within 10 hrs. Since HF8 was performed closest to the GE-MIMs and induced the creation of a single fracture plane that rapidly connected to the transmissive S3 shear zone²¹,²⁶, the dispersion of (He, Ar)-enriched fluids is expected to be lower than in the case of HF3.

## Conclusion

In this paper, we provide *in situ* evidence of He and Ar anomalies occurring during the deformation of a crystalline fractured system at the decameter scale. We observed clear differences in the magnitude and temporal dynamics of these anomalies for different clusters of stimulation experiments that displayed differences in their hydro-seismo-mechanical responses. Complementing previous hypotheses drawn from laboratory experiments, we speculate that the remobilization of (He, Ar)-enriched fluids is principally controlled by (1) the extent of the deformed rock mass and the area of the ruptured fractures, which may control the amount of (He, Ar)-enriched fluids released; (2) the type of deformation, which may lead to different concentration anomalies. The generation of a single fracture plane in intact rocks leads to a stronger (He, Ar) signal, typically experiment HF8, than in the case of the reactivation of multiple fractures that may disperse the signal, typically experiment HF3; (3) the efficiency of the newly generated fracture of connecting with the preexisting transmissive fracture system, allowing its transport within the reservoir; and (4) the distance of the source of (He, Ar)-enriched fluids from the monitoring point.

Overall, high-frequency noble gas measurements in fractured rocks may help to improve the conceptualization of the nature and evolution of flow paths in response to stimulation operations, and complement the information obtained from hydraulic, deformation and seismic monitoring systems. Specifically, the analysis of the magnitude and temporal scaling of hydrogeochemical changes, such as the described noble gas anomalies, may inform about the evolution of the hydraulic connectivity and the transport properties of the stimulated fracture network, which is information that is not accessible through classic seismic and geodetic methods.

From similar perspectives, valuable information could be gained through high-frequency noble gas monitoring in both CO₂ injection and oil and gas extraction projects, where undesired gas migration pathways toward shallower aquifers could be identified during operation³⁴,³⁵. Monitoring noble gas concentrations may also help understand the impact of natural seismic events on regional aquifer connectivity, providing critical information to link seismicity to fluid migration in the Earth's crust.

## Methods

**Injection protocol.** Six high-pressure fluid injections were performed in two injection boreholes (INJ1-2 in Fig. 1, labeled HF#) by isolating 1 m of intact rock with double packer systems. The targeted injection intervals had a constant base length of 1 m. Prior to each HF test, the interval integrity and proper packer sealing were tested with a short pulse test. HF tests started with a flow rate-controlled (5 l/min) injection cycle that lasted for 30 seconds to initiate fracturing (i.e., formation breakdown), followed by shut-in and bleed-off phases. Subsequently, the fracture was propagated with flow rate-controlled injections (up to 100 l/min). A total of 1,000 liters of fluid was injected during each test over all cycles. Each injection cycle was followed by a shut-in phase. During HF5, HF6 and HF8, water and a rheology modifier (mixture with xanthan gum polymer) were used to achieve a viscosity-dominated propagation regime²¹. The latter consists of a xanthan-salt-water mixture providing 35 times higher viscosity than water. For the tests undertaken with xanthan-salt-water, a supplementary flushing cycle was performed to flush as much xanthan as possible out of the fracture. The last cycle was a pressure-controlled step test to evaluate the post-stimulation injectivity of the created hydraulic fracture and to estimate the stress acting normal to the hydraulic fracture (jacking pressure). Note that HF1 is not considered in this study because it suffered from logistical issues that prevent any reliable monitoring and interpretation.

**Deformation monitoring.** To obtain quantitative information about the deformation field inside the test volume, three boreholes (referred to as FBS boreholes in Fig. 1) were drilled, and each borehole was equipped
with 20 Fiber-Bragg-Grating strain sensors. The sensors provide one-dimensional borehole-parallel information about local deformations. Prior to installation, fractures were mapped along the borehole from optical and acoustic televiewer logs³⁶. Care was taken to ensure that the strain sensors covered the different geological features of interest, i.e., including preexisting fractures, shear zones and intact rock sections. The sensors have a length of 1 m, accuracy of 1 microstrain, and resolution of 0.1 microstrain and were measured at a sampling rate of 1000 Hz. The strain signals described in this paper (from borehole FBS1 at 22.35 m and 42.2 m depths) each cover one of the main shear zones (S3.1 and S1.3) targeted by the stimulations, providing information about the extent of deformation and its temporal variability.

**Seismic monitoring.** Seismicity was monitored during the injection²⁶, shut-in and bleed-off phases of the six HF experiments using 26 highly sensitive piezoelectric acoustic emission sensors. Eight of the acoustic emission sensors were deployed in four monitoring boreholes in proximity to the injection intervals. Seismicity was recorded continuously at a sampling rate of 200 kHz. Seismic event detection was performed using a standard recursive short-time-average/long-time-average coincidence trigger, whereby picking was carried out manually. Seismic event location was performed in an anisotropic, homogeneous velocity model. An amplitude-based magnitude $M_A$ was estimated for located seismic events and corrected for angle-dependent sensitivity variations and coupling quality of acoustic emission sensors. The amplitude-based magnitudes $M_A$ were finally adjusted to absolute magnitudes $M_W$ determined for a subset of larger magnitude seismic events. The absolute magnitudes $M_W$ were estimated from acoustic emission sensors installed in nearby tunnels that were collocated with calibrated accelerometers. The collocation of acoustic emission sensors and accelerometers allowed for a cross-calibration of the acoustic emission sensors. Finally, to determine absolute magnitudes $M_W$, the theoretical displacement source spectrum introduced by Boatwright (1978)³⁷ was fitted to the recorded displacement spectrum on the cross calibrated acoustic emission sensors²⁶.

**Continuous noble gas monitoring.** A portable gas equilibrium membrane inlet mass spectrometer (GE-MIMS, Gasometrix GmbH) was installed at a nearby natural seeping zone related to the transmissive S3 shear zones (Fig. 1c). The instrument allows for near-real-time on-site environmental quantification of partial pressures and hence the concentrations of gases dissolved in water²⁴. The fluid was collected from the whole area of the seeping shear zone. The fracture was sealed with a plastic foil and waterproof mortar, forcing the fluid to drain by gravity to the GE-MIMS. The sealing was designed to avoid any secondary gas exchange with the atmospheric gases. Background monitoring with GE-MIMS revealed anoxic conditions, confirming that the sampling system was airtight. Then, the collected water was transferred to GE-MIMS membrane contactors (two parallel Micro Modules 3 M Liqui-Cel MM-0.5–1×1 Series), which are specifically designed for low flow rates (from 5–1000 ml/min). Water temperature and total gas pressure in the membranes’ gas head space are determined by the respective sensors. The small head space further maintains a high temporal resolution at the respective water flow rate. He, Ne (in dry gas), Ar, Kr, $N_2$, and $O_2$ were analyzed every 8 min with an analytical uncertainty of 1–3%. Gas standards from an air-filled Plastigas bag were analyzed every 40 min to calibrate the sample data. Natural background gas concentrations were measured for two days prior to the start of stimulation activities, from May 12th to May 14th (Fig. 2).

The flow rate of the seeping pore water was measured at the outflow of the membrane contactors using a tipping bucket flow gauge connected to a Siemens micrologger starting on May 17 at 6 am before the initiation of HF2. The monitored seeping fracture provided flow rates ranging from 20 to 160 mL/min. These flow rates are lower than the recommended value for air-equilibrated water to maintain the mass balance of the gases in the membrane contactor at gas/water²⁴, implying that the partial pressures of the different gas species in the head-space could be controlled by diffusion-limited gas/water exchange across the membranes themselves. This would lead to incomplete gas/water equilibration, especially during periods of low water flow. To avoid any potential effects arising from incomplete gas/water equilibrium within the membrane contactor, we normalized the noble gas concentrations by the measured concentrations of dissolved $N_2$, which is only of atmospheric origin in our case study. With ratios of aqueous diffusion coefficients of 3.34 (He/$N_2$) and 1.06 (Ar/$N_2$), diffusive fractionation would result in an excessive He/$N_2$ ratio compared with Ar/$N_2$ if the mass balance of the gases in the membrane contactors would be limited by diffusion across the membrane. However, the reverse was observed in the dataset, which indicates that the observed He/$N_2$ and Ar/$N_2$ ratios are not fractionated by non-equilibrium effects in the membrane contactors, most likely due to the higher gas concentrations in the fluid compared to air-equilibrated water. Therefore, the observed He/$N_2$ and Ar/$N_2$ ratios accurately reflect the ratios of dissolved He, Ar and $N_2$ in the fluid.

The concentrations of He and Ne, as well as the ³He/⁴He ratios of individual water samples obtained in standard copper tubes, were analyzed at the noble gas laboratory at ETHZ according to standard analytical protocols. For technical details, the reader is referred to the work of Beyerle *et al*.³⁸

## Data availability

The Grimsel ISC Experiment Description is available at https://doi.org/10.3929/ethz-b-000310581. Data used in this publication are available from the ETH Zürich research collection. Partial pressures and discharge rates measured at the GE-MIMS station can be found at https://doi.org/10.3929/ethz-b-000365199. The seismic dataset, as well as hydraulic data of the Grimsel ISC hydraulic fracturing experiments, are available at https://doi.org/10.3929/ethz-b-000276170 and the hydromechanical dataset is available at https://doi.org/10.3929/ethz-b-000328270.

Received: 11 December 2019; Accepted: 31 March 2020;
Published online: 24 April 2020
## References

1. Zoback, M. D. Reservoir Geomechanics, https://doi.org/10.1017/CBO9780511586477 (Cambridge University Press, 2007).

2. Claesson, L. et al. Hydrogeochemical changes before and after a major earthquake. *Geology* **32**, 641 (2004).

3. Skelton, A. et al. Changes in groundwater chemistry before two consecutive earthquakes in Iceland. *Nat. Geosci.* **7**, 752–756 (2014).

4. Wakita, H., Igarashi, G. & Notsu, K. An anomalous radon decrease in groundwater prior to an M6.0 earthquake: A possible precursor? *Geophys. Res. Lett.* **18**, 629–632 (1991).

5. Sano, Y. et al. Groundwater helium anomaly reflects strain change during the 2016 Kumamoto earthquake in Southwest Japan. *Sci. Rep.* **6**, 1–7 (2016).

6. Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E. & Cardenas, M. B. The global volume and distribution of modern groundwater. *Nat. Geosci.* **9**, 1–15 (2015).

7. Holland, G. et al. Deep fracture fluids isolated in the crust since the Precambrian era. *Nature* **497**, 357–360 (2013).

8. Aquilina, L. et al. Impact of climate changes during the last 5 million years on groundwater in basement aquifers. *Sci. Rep.* **5**, 14132 (2015).

9. Kipfer, R., Aeschbach-Hertig, W., Peeters, F. & Stute, M. Noble gases in lakes and ground waters. *Noble gases geochemistry Cosmochem. Rev. Mineral. Geochemistry* **47**, 615–700 (2002).

10. *Environmental Tracers in Subsurface Hydrology*. Cook, Peter G. Herczeg, Andrew L., https://doi.org/10.1007/978-1-4615-4557-6 (Springer US, 2000).

11. Lowenstern, J. B., Evans, W. C., Bergfeld, D. & Hunt, A. G. Prodigious degassing of a billion years of accumulated radiogenic helium at Yellowstone. *Nature* **506**, 355–358 (2014).

12. Kennedy, B. M. & van Soest, M. C. Flow of Mantle Fluids Through the Ductile Lower Crust: Helium Isotope Trends. *Science* (80–). **318**, 1433–1436 (2007).

13. Barry, P. H. et al. Determining fluid migration and isolation times in multiphase crustal domains using noble gases. *Geology* **45**, 775–778 (2017).

14. *The Noble Gases as Geochemical Tracers.*, https://doi.org/10.1007/978-3-642-28836-4 (Springer Berlin Heidelberg, 2013).

15. Skelton, A. et al. Hydrochemical changes before and after earthquakes based on long term measurements of multiple parameters at 2 sites in northern Iceland - a review. *J. Geophys. Res. Solid Earth* 1–19, https://doi.org/10.1029/2018JB016757 (2019).

16. Manga, M. et al. Changes in permeability caused by transient stresses: Field observations, experiments, and mechanisms. *Rev. Geophys.* **50** (2012).

17. Bauer, S. J., Gardner, W. P. & Lee, H. Release of radiogenic noble gases as a new signal of rock deformation. *Geophys. Res. Lett.* **43**(10), 688–10,694 (2016).

18. Bauer, S. J., Gardner, W. P. & Heath, J. E. Helium release during shale deformation: Experimental validation. *Geochemistry, Geophys. Geosystems* **17**, 2612–2622 (2016).

19. Gardner, W. P., Bauer, S. J., Kuhlman, K. L. & Heath, J. E. Modeling Dynamic Helium Release as a Tracer of Rock Deformation. *J. Geophys. Res. Solid Earth* 1–15, https://doi.org/10.1002/2017JB014376 (2017).

20. Honda, M., Kurita, K., Hamano, Y. & Ozima, M. Experimental studies of He and Ar degassing during rock fracturing. *Earth Planet. Sci. Lett.* **59**, 429–436 (1982).

21. Dutler, N. et al. Hydraulic fracture propagation in a heterogeneous stress field in a crystalline rock mass. *Solid Earth Discuss.* 1–41, https://doi.org/10.5194/se-2019-111 (2019).

22. Amann, F. et al. The seismo-hydromechanical behavior during deep geothermal reservoir stimulations: open questions tackled in a decameter-scale *in situ* stimulation experiment. *Solid Earth* **9**, 115–137 (2018).

23. Krietsch, H. et al. Comprehensive geological data of a fractured crystalline rock mass analog for hydraulic stimulation experiments. *Sci. Data* **5**, 1–12 (2018).

24. Brennwald, M. S., Schmidt, M., Oser, J. & Kipfer, R. A Portable and Autonomous Mass Spectrometric System for On-Site Environmental Gas Analysis. *Environ. Sci. Technol.* **50**, 13455–13463 (2016).

25. Shapiro, S. A., Rothert, E., Rath, V. & Rindschwentner, J. Characterization of fluid transport properties of reservoirs using induced microseismicity. *Geophysics* **67**, 212–220 (2002).

26. Villiger, L. et al. Influence of reservoir geology on seismic response during decameter scale hydraulic stimulations in crystalline rock. *Solid Earth Discuss.* **2019**, 1–46 (2019).

27. Heaton, T. H. E. & Vogel, J. C. "Excess air" in groundwater. *J. Hydrol.* **50**, 201–216 (1981).

28. Holocher, J. et al. Experimental investigations on the formation of excess air in quasi-saturated porous media. *Geochim. Cosmochim. Acta* **66**, 4103–4117 (2002).

29. Clarke, W., Top, Z. & Eismont, W. Helium Isotope and Tritium Content of Lakes and Uranium Exploration in the NEA/IAEA Test area. in *Uranium Exploration in Athabasca Basin, Saskatchewan, Canada* (ed. Cameron, E. M.) 139–146 (1983).

30. Frick, U. et al. *Grimsel Test Site - The Radionuclide Migration Experiment - Overview of Investigations 1985-1990*. PSI Bericht Nr. 120 (Nagra NTB 91-04). (1992).

31. Baertschi, P., Antonsen, O., Meier, W. R. & Keil, R. Untersuchungen an Wasserproben der Sondierbohrunoen SB1-SB4 im Hauptstollen des Pumpkraftwerks Grimsel/Oberaar (Nagra Auftrag), EIR Internal Technical Report TM-44-82-13. (1982).

32. Chatton, E. et al. Field continuous measurement of dissolved gases with a CF-MIMS: Applications to the physics and biogeochemistry of groundwater flow. *Environ. Sci. Technol.* https://doi.org/10.1021/acs.est.6b03706, (2016).

33. Fetter, C. W. *Applied Hydrogeology*. (Pretince Hall, 2001).

34. Jackson, R. B. et al. Increased stray gas abundance in a subset of drinking water wells near Marcellus shale gas extraction. *Proc. Natl. Acad. Sci. USA* **110**, 11250–11255 (2013).

35. White, J. A. et al. Geomechanical behavior of the reservoir and caprock system at the in Salah CO₂ storage project. *Proc. Natl. Acad. Sci. USA* **111**, 8747–8752 (2014).

36. Krietsch, H. et al. Comprehensive geological dataset describing a crystalline rock mass for hydraulic stimulation experiments. *Sci. Data* **5**, 180269 (2018).

37. Boatwright, J. Detailed spectral analysis of two small New York state earthquakes. *Bull. Seismol. Soc. Am.* **68**, 1117–1131 (1978).

38. Beyerle, U. R. S. & Du, C. A Mass Spectrometric System for the Analysis of Noble Gases and Tritium from Water Samples. *Environ. Sci. Technol.* **34**, 2042–2050 (2000).

## Acknowledgements

The ISC is a project of the Deep Underground Laboratory at ETH Zurich, established by the Swiss Competence Center for Energy Research - Supply of Electricity (SCCER-SoE) with the support of Innosuisse. Funding for the ISC project was provided by the ETH Foundation with grants from Shell and EWZ and by the Swiss Federal Office of Energy. Nathan Dutler and the hydraulic fracture tests of the ISC are supported by SNF Grant No. 200021_165677. We are indebted to Nagra for hosting the ISC experiment in their GTS facility and to the Nagra technical staff, especially René Dorrer, for onsite support. We are grateful to two anonymous reviewers whose comments have greatly improved the content of this paper.
## Author contributions

C.R. and R.K. initiated this research; J.D., V.G., F.A., and B.V. initiated the *in situ* stimulation and circulation (ISC) experiment; C.R., R.K., U.W. and M.B. developed the geochemical monitoring setup; U.W. and C.R. installed the GE-MIMS and collected samples; J.D., V.G., F.A., M.K. and B.V. designed the hydraulic, deformation and seismic monitoring setup; N.D., M.J., J.D., V.G., F.A., and B.V. designed and conducted the hydraulic fracturing experiment; C.R., U.W., M.B. and R.K. processed and interpreted the geochemical data; H.K. and N.D. processed and interpreted the deformation data; L.V. processed and interpreted the seismic data; N.D., B.B., M.J. and M.K. processed and interpreted the hydraulic data; C.R., N.D., U.W., B.B., H.K. and J.D. interpreted the overall dataset; H.K. created Figure 1; C.R. and U.W. created Figures 2 and 3; L.V. created Figure 1SM; C.R., U.W., B.B., N.D., H.K. and L.V. wrote the initial manuscript; and all coauthors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-63458-x.

**Correspondence** and requests for materials should be addressed to C.R.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or

format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020