Underground Rock Laboratory Research Articles

Retention of 226Ra in the sandy Opalinus Clay facies from the Mont Terri rock laboratory, Switzerland

The safety assessment of deep geological disposal of nuclear waste requires a sound understanding of radionuclide retention in the repository host rocks. Here, we investigate the retention behaviour of the safety relevant radionuclide 226Ra in the heterogeneous sandy facies of the Opalinus Clay (OPA-SF) at the Mont Terri underground rock laboratory. Various rock samples were selected from a drill core (BAD-1) located in the lower sandy facies, which were representative for its textural and mineralogical heterogeneity. The samples were carefully characterized with respect to mineralogy and chemistry using X-Ray diffraction and electron microscopy. The amounts of quartz, carbonates and clay minerals in the samples ranged between 58 wt.% – 65 wt.%, 13 wt.% – 35 wt.% and 5 wt.% – 25 wt.%, respectively, underpinning the mineralogical heterogeneity. Rare (Ba,Sr)-sulphate precipitates were identified in the clay matrix by electron microscopy. The uptake behaviour of Ra was evaluated using batch sorption experiments. Corresponding to the mineralogical variability of the rock samples, the measured distribution ratios Rd for the 226Ra retention varied between 20 L kg−1 and 260 L kg−1 after 120 days. The experimental results and thermodynamic modelling approaches suggest that various processes contribute to the uptake of 226Ra on different time scales: fast cation exchange and surface complexation reactions by clay minerals and the (slower) formation of solid solutions during recrystallisation of (Ba,Sr)-sulphates as well as retention by carbonate minerals; the latter mechanism still requiring further investigation.

Cation exchange parameters for Opalinus Clay and its confining units

The knowledge of cation exchange parameters is important to understand and model the porewater chemistry in clayrocks. The Opalinus Clay (OPA) is a well-studied formation in view of its importance as host rock for radioactive waste repositories, but uncertainties exist regarding its cation exchange parameters, in particular its selectivity coefficients (Kc). In fact, porewater chemistry models for OPA have so far been based on generic selectivity coefficients.Thanks to a recent extensive siting programme, involving eight vertical deep boreholes in northern Switzerland, a large number of porewater and cation exchange data from the OPA and its clay-rich confining units could be analysed. Thus, from the combination of porewater compositions (obtained from high-pressure squeezing and advective displacement), and exchangeable cation data (obtained from Ni ethylenediamine extraction), reliable selectivity coefficients could be derived. The values for KcNa/K, KcNa/Ca, KcNa/Mg are consistent with those obtained in parallel from the Mont Terri underground rock laboratory although the scatter is somewhat larger in the latter dataset. The derived selectivity coefficients indicate no dependency on ionic strength or in-situ temperature. Minor differences in Kc values between OPA and confining units can be attributed to differences in clay mineralogy. The validity of a two-site cation exchange model including illite and smectite for estimating exchangeable Na, Ca and Mg from measured porewater compositions could be shown, while this simple model underpredicts the share of exchangeable K. The exchangeable cation extraction procedure induces the dissolution of Sr-bearing minerals whose amount is in the same range as exchangeable Sr, thus not enabling the derivation of selectivity coefficients for Sr.Finally, derived selectivity coefficients confirm results of previous studies, thus providing confidence in the applied experimental and modelling protocols on the one hand and in the outcome of those studies on the other hand.

Gas equilibrium membrane inlet mass spectrometry (GE-MIMS) for water at high pressure

Abstract. Gas species are widely used as natural or artificial tracers to study fluid dynamics in environmental and geological systems. The recently developed gas equilibrium membrane inlet mass spectrometry (GE-MIMS) method is most useful for accurate and autonomous on-site quantification of dissolved gases in aquatic systems. GE-MIMS works by pumping water through a gas equilibrator module containing a gas headspace, which is separated from the water by a gas-permeable membrane. The partial pressures of the gas species in the headspace equilibrate with the gas concentrations in the water according to Henry's Law and are quantified with a mass spectrometer optimized for low gas consumption (miniRUEDI or similar). However, the fragile membrane structures of the commonly used equilibrator modules break down at water pressures ≳3 bar. These modules are therefore not suitable for use in deep geological systems or other environments with high water pressures. To this end, the SysMoG® MD membrane module (Solexperts AG, Switzerland; “SOMM”) was developed to withstand water pressures of up to 100 bar. Compared to the conventionally used GE-MIMS equilibrator modules, the mechanically robust construction of the SOMM module entails slow and potentially incomplete gas–water equilibration. We tested the gas equilibration efficiency of the SOMM and developed an adapted protocol that allows correct operation of the SOMM for GE-MIMS analysis at high water pressures. This adapted SOMM GE-MIMS technique exhibits a very low gas consumption from the SOMM to maintain the gas–water equilibrium according to Henry's Law and provides the same analytical accuracy and precision as the conventional GE-MIMS technique. The analytical potential of the adapted SOMM GE-MIMS technique was demonstrated in a high-pressure fluid migration experiment in an underground rock laboratory. The new technique overcomes the pressure limitations of conventional gas equilibrators and thereby opens new opportunities for efficient and autonomous on-site quantification of dissolved gases in high-pressure environments, such as in research and monitoring of underground storage of CO2 and waste deposits or in the exploration of natural resources.

Characterization of heterogeneities in the sandy facies of the Opalinus Clay (Mont Terri underground rock laboratory, Switzerland)

SUMMARY This contribution is presenting a multidisciplinary investigation of heterogeneities in a clay rock formation, based on seismic tomography, logging and core analysis, as a reconnaissance study for a diffusion experiment. Diffusion experiments in clay rock formations provide crucial experimental data on diffusive transport of radionuclides (RN) in extremely low hydraulic conductivity media. Previous diffusion experiments, conducted, for example, in the Mont Terri underground rock laboratory within the relatively homogeneous shaly facies of Opalinus Clay, and modelling studies of these experiments have demonstrated that the clay rock could sufficiently well be described as a homogeneous anisotropic medium. For other lithofacies, characterized by larger heterogeneity, such simplification may be unsuitable, and the description of heterogeneity over a range of scales will be important. The sandy facies of the Opalinus Clay exhibits a significantly more pronounced heterogeneity compared to the shaly facies, and a combined characterization and RN diffusion study has been initiated to investigate various approaches of heterogeneity characterization and subsequent diffusion in a heterogeneous environment. As an initial step, two inclined exploratory boreholes have been drilled to access the margins of the experiment location. These boreholes have been used to acquire a cross-hole tomographic seismic data set. Optical, natural gamma and backscattering logging were applied and rock cores were analysed. The integrated results of these investigations allowed the identification of an anomalous brighter layer within the investigated area of the sandy facies of approximately 1 m thickness and with its upper bound at roughly 10 m depth within the inclined exploratory wells. Mineralogical analyses revealed only slight variations throughout the rock cores and indicated that the anomalous layer exhibited a slightly higher quartz content, and locally significantly higher calcite contents, accompanied by a lower content of clay minerals. The anomalous layer was characterized by reduced natural gamma emissions, due to the lower clay content, and increased neutron backscattering likely indicating an increased porosity. Seismic P-wave velocities, derived from anisotropic tomography, exhibited a maximal gradient near the top of this layer. The transition from the overlaying darker rock matrix into this layer has been identified as an appropriate location for the setup of a tracer diffusion experiment in a heterogeneous environment.

Multi-disciplinary characterization of shrinkage effects in Opalinus clay

Multi-disciplinary characterization of shrinkage effects in Opalinus clay

Carbon steel corrosion in a bentonite buffer: A comparison between in situ exposure and lab based experiments

AbstractIn this paper, carbon steel corrosion rates from experiments performed in situ (performed in an underground rock laboratory) and ex situ (performed in a conventional laboratory) test methods in anaerobic saturated bentonite are compared. The results indicate that the long‐term corrosion rate follows a power law decay curve, with a higher initial rate and greater rate of decay at higher temperature. In compacted bentonite blocks, varying the density has no significant effect on the corrosion rate measured during in situ testing, whereas granular bentonite of an equal dry density leads to a higher corrosion rate. Precorrosion of test specimens in an aerated solution at room temperature has a negligible effect on the subsequent anaerobic corrosion rate when tested ex situ for durations up to 11.8 years. The main difference between in situ and ex situ tests is the formation of silicon‐rich corrosion products in the ex situ tests but not in the in situ tests. Despite these differences, the corrosion rates in both tests exhibit the same general temporal evolution and similar magnitudes, suggesting that the main steel degradation processes are maintained in both configurations.

Performance of high-resolution impact and vibration sources for underground seismic exploration of clay formations at the scale of a rock laboratory

SUMMARY Low permeability, high retention capacity and self-sealing ability are advantageous characteristics that are attributed to argillaceous rocks. In contrast, other properties of clay, such as internal heterogeneities, strong attenuation and anisotropic behaviour, pose major challenges for underground exploration techniques. Although with regard to the underground storage of nuclear waste, the seismic exploration in the underground itself is of great importance to fill the gap between surface and borehole investigations. Furthermore, to prevent destruction of the host rock during exploration this demands low to non-invasive techniques. To approach these issues, a seismic survey was carried out in the Mont Terri Underground Rock Laboratory (Switzerland) using a gallery-based acquisition with an operating range up to several decametres. The seismic campaign included three-component borehole sensors and two different seismic source types (pneumatic impact and magnetostrictive vibroseis source). An executed source comparison analysed the characteristics of the different source types, for example frequency or amplitude behaviour of the generated wavefield, to assess their performance under similar conditions at the meso scale and to reveal their strengths and weaknesses in clay. Based on these findings, we performed traveltime and reflection analyses that demonstrate their potential to characterize clay formations and to map internal structures. The highest seismic velocities are found in the carbonate-rich sandy facies (vPmax = 4000 m s−1, vSmax = 2050 m s−1), slower velocities are found in the sandy facies (vPmax = 3720 m s−1, vSmax = 1840 m s−1) and the slowest velocities are found in the shaly facies (vPmax = 3220 m s−1, vSmax = 1480 m s−1). The seismic velocity anisotropy is larger within the shaly facies (AvP = 23 per cent, AvS = 32 per cent) compared to the sandy facies (AvP = 9 per cent, AvS = 12 per cent) and it is more pronounced for S-waves than P-waves. Thus, non-invasive meso-scale seismic techniques are suited to characterize the Opalinus Clay in great detail.

Fracture growth leading to mechanical spalling around deposition boreholes of an underground nuclear waste repository

This study presents a three-dimensional numerical analysis of multiple fracture growth leading to spalling around nuclear waste deposition boreholes. Mechanical spalling due to stress amplification after drilling is simulated using a finite element-based fracture growth simulator. Fractures initiate in tension based on a damage criterion and grow by evaluating stress intensity factors at each fracture tip. Tip propagation is multi-modal, resulting in final fracture patterns that are representative of both tensile and shear failure. Their geometries are represented by smooth parametric surfaces, which evolve during growth using lofting. The corresponding surface and volumetric meshes are updated at every growth step to accommodate the evolving fracture geometries. The numerical model is validated by comparing simulated fracture patterns against those observed in the AECL Underground Rock Laboratory Mine-By Experiment. It is subsequently calibrated to simulate fracture initiation and growth around boreholes drilled in the Forsmark granodiorite, subjected to a far-field anisotropic triaxial stress that corresponds to the in situ stress model from the Swedish Forsmark site. The deposition tunnel is implicitly simulated by attaching the deposition borehole to a free domain boundary. Several geomechanical cases are investigated, in which fracture growth is numerically evaluated as a function of in situ stress state, tunnel orientation, borehole geometry, total number of boreholes and borehole spacing. Numerical results show that spalling occurs in all cases, given the underground conditions at Forsmark, with borehole geometry, spacing and stresses affecting the extent of fracture nucleation and growth patterns. The uncertainty in underground stress conditions is evaluated through varying stress magnitudes and orientations relative to the tunnel floor. Whereas tunnel orientation influences the relative locations where fractures initiate with respect to the tunnel floor, fracture growth and its final extent is determined by the relative magnitudes of the in situ stresses. Higher stress differential causes higher spalling depths, but in all cases, failure extent is localised to the borehole vicinity, not exceeding one borehole radius. The cylindrical borehole is modified at the top to provide an access ramp for the spent fuel canisters and fracture growth around several deposition boreholes is simulated for borehole tops having cylindrical, conical, and wedge shapes. The enlargement of the borehole top induces higher stress concentrations at the borehole–tunnel junction, increasing the severity of spalling closer to the tunnel floor. Massive failure occurs when a multiple borehole model is considered and the inter-borehole distance is small enough that adjacent “spalled” areas interact. At Forsmark, through-going fractures are predicted to develop when the borehole spacing is less than 4 m. The effect of spalling on the structural integrity of the deposition boreholes is illustrated for each test case and quantified in terms of maximum spalling depth, spalling width and total fractured surface area.

Hydrogen gas transfer between a borehole and claystone: experiment and geochemical model

The Hydrogen Transfer experiment, located at the Mont Terri underground rock laboratory in Switzerland, is an in situ experimental study of the interactions and transport of hydrogen injected into a borehole installed within Opalinus Clay, a claystone formation. A Python-based model was developed to analyse and model the experimental data, for diffusion of dissolved gases and solutes in claystone pore water, for thermodynamic modelling of gas–water–solid phase equilibria in the injection interval and for simulations of chemical equilibria and reaction kinetics in claystone and injection interval. The developed model reproduces the temporal evolution of gas pressure, composition and solute concentrations measured in situ with a minimum set of adjustable parameters. The effective diffusion coefficients for dissolved gases in Opalinus Clay derived by the modelling of experimental data were found to be very close to values measured in other experimental studies. It was discovered that an accurate description of the temporal variations in hydrogen injection and temporal changes in the inflow of formation water is essential for modelling of microbially mediated hydrogen consumption in the injection interval.

Evolution of excavation damaged zones in Opalinus Clay shale inferred from seismic investigations

We present the results of new and repeat gallery- and borehole-scale seismic experiments carried out in two excavation damaged zones (EDZs) of 10 (gallery 08) and 20 (gallery 98) years of age at the Mont Terri Underground Rock Laboratory in St. Ursanne, Switzerland. The initial geometry and rock mass properties of the EDZs were studied during and shortly after excavation—allowing for a direct comparison of measurements collected in 2018 to those completed in the past.In gallery 08 parts of the EDZ remain unchanged over ten years, whilst other zones show improvement (i.e., self-sealing) or considerable degradation. In addition, velocity tomograms from 2018 measure extensive areas of p-wave velocity reduction in comparison to the surrounding rock mass at depths shallower than 3.5 m. Isolated extensile fractures are still found in boreholes up to a maximum radial extent of 4.2 m. Data from borehole interval velocity measurements and cross-hole seismic surveys from this gallery suggest the damaged zone still extends up to 1.9 m and 3.1 m radial depth.Refraction seismic results in gallery 98, in contrast, suggest a less pervasive, shallower (<1.0 m) EDZ with a lower frequency of extensile fractures. The smaller circumference of this gallery can in part explain this; however, when compared to various measurements collected during and after excavation, some reduction in EDZ extent is inferable. Both borehole interval velocity measurements and seismic refraction data record p-wave velocity values higher on average than those in gallery 08 and without obvious areas of velocity reduction. Most interval velocity measurements likewise do not measure zones of decreased velocity (e.g., corresponding to rock mass damage) at this location. This suggests more self-sealing has occurred over the twenty-year post-excavation period in gallery 98 in comparison to the ten-year period in gallery 08. As repository safety depends upon the integrity of the rock mass surrounding an excavation, these results are relevant for mid-term assessments of safety (i.e., especially prior to repository closure).

Seismic anisotropy of Opalinus Clay: tomographic investigations using the infrastructure of an underground rock laboratory (URL)

Seismic anisotropy and attenuation make claystone formations difficult to characterize. On the other hand, in many geotechnical environments, precise knowledge of structure and elastic properties of clay formations is needed. In crystalline and rock salt underground structures, high-resolution seismic tomography and reflection imaging have proven a useful tool for structural and mechanical characterization at the scale of underground infrastructure (several deca- to hundreds of meters). This study investigates the applicability of seismic tomography for the characterization of claystone formations from an underground rock laboratory under challenging on-site conditions including anisotropy, strong attenuation and restricted acquisition geometry. The seismic tomographic survey was part of a pilot experiment in the Opalinus Clay of the Mont Terri Rock Laboratory, using 3-component geophones and rock anchors, which are installed 2 m within the rock on two levels, thus suppressing effects caused by the excavation damage zone. As a source, a pneumatic impact source was used. The survey covers two different facies types (shaly and carbonate-rich sandy), for which the elliptical anisotropy is calculated for direct ray paths by fitting an ellipse to the separated data for each facies. The tomographic inversion was done with a code providing a good grid control and enabling to take the seismic anisotropy into account. A-priori anisotropy can be attributed to the grid points, taking various facies types or other heterogeneities into account. Tomographic results, compared to computations using an isotropic velocity model, show that results are significantly enhanced by considering the anisotropy and demonstrate the ability of the approach to characterize heterogeneities of geological structures between the galleries of the rock laboratory.

Unloading induced absolute negative pore pressures in a low permeable clay shale

The stability analysis of underground excavations constructed in low-permeable rocks requires a comprehensive understanding of the short- and long-term pore pressure response around the excavation. This paper presents the results of a novel in-situ experiment, which was conducted in Opalinus Clay at the Mont Terri Underground Rock Laboratory, Switzerland. The primary aim of the project was to understand unloading-induced pore pressure changes in Opalinus Clay with a focus on the possibility of the development of absolute negative pore pressures. In this experiment, the pore pressure response was initiated by overcoring a small diameter pilot borehole equipped with a pressure transducer and a full-range tensiometer. The results of the experiment showed that the stress relief was accompanied by a noticeable pore pressure drop in the core. An absolute negative pore pressure of −755 kPa was measured. To replicate the pore pressure evolution during the in-situ experiment, a 3D hydro-mechanical-coupled finite element model was developed. It is shown that the magnitude of the pore pressure drop cannot be explained by the classical poroelasticity theory for fully saturated Opalinus Clay samples. Good agreement between the experimental data and numerical results can only be achieved when unloading induced desaturation is considered. Rapid unloading may cause pore water cavitation and desaturation, which significantly increases the pore fluid compressibility. Our considerations based on the nucleation theory suggest that during unloading the required conditions for pore water cavitation in pores larger than a critical value were met.

Mechanical and hydraulic properties of the excavation damaged zone (EDZ) in the Opalinus Clay of the Mont Terri rock laboratory, Switzerland

Abstract. Construction of cavities in the subsurface is always accompanied by excavation damage. Especially in the context of deep geological nuclear waste disposal, the evolving excavation damaged zone (EDZ) in the near field of emplacement tunnels is of utmost importance concerning safety aspects. As the EDZ differs from the intact host rock due to enhanced hydraulic transmissivity and altered geomechanical behavior, reasonable and location-dependent input data on hydraulic and mechanical properties are crucial. Thus, in this study, a hydromechanical characterization of an EDZ in the Mont Terri underground rock laboratory, Switzerland, was performed using three different handheld devices: (1) air permeameter, (2) microscopic camera and (3) needle penetrometer. The discrete fracture network (DFN), consisting of artificially induced unloading joints and reactivated natural discontinuities, was investigated by a portable air permeameter and combined microscopic imaging with automatic evaluation. Geomechanical and geophysical characterization of the claystone was conducted based on needle penetrometer testing at the exposed rock surface. Within the EDZ, permeable fractures with a mean hydraulic aperture of 84 ± 23 µm are present. Under open conditions, self-sealing of fractures is suppressed, and cyclic long-term fracture aperture oscillations in combination with closure resulting from convergence processes is observed. Based on measured needle penetration indices, a uniaxial compressive strength of 30 ± 13 MPa (normal to bedding) and 18 ± 8 MPa (parallel to bedding) was determined. Enhanced strength and stiffness are directly related to near-surface desaturation of the claystone and a sharp decrease in water content from 6.6 wt % to 3.7 wt %. The presented methodological approach is particularly suitable for time-dependent monitoring of EDZs since measurements are nondestructive and do not change the actual state of the rock mass. This allows for a spatially resolved investigation of hydraulic and mechanical fracture apertures, fracture surface roughness, and physico-mechanical rock parameters and their intra-facies variability.

Reactive Transport Simulation of Low-pH Cement Interacting with Opalinus Clay Using a Dual Porosity Electrostatic Model

Strong chemical gradients between clay and concrete porewater lead to diffusive transport across the interface and subsequent mineral reactions in both materials. These reactions may influence clay properties such as swelling behaviour, permeability or radionuclide retention, which are relevant for the safety of a radioactive waste repository. Different cement types lead to different interactions with Opalinus Clay (OPA), which must be understood to choose the most suitable material. The consideration of anion-depleted porosity due to electrostatic repulsion in clay modelling substantially influences overall diffusive transport and pore clogging at interfaces. The identical dual porosity model approach previously used to predict interaction between Portland cement and OPA is now applied to low-alkali cement—OPA interaction. The predictions are compared with corresponding samples from the cement-clay interaction (CI) experiment in the Mont Terri underground rock laboratory (Switzerland). Predicted decalcification of the cement at the interface (depletion of C–S–H and absence of ettringite within 1 mm from the interface), the Mg enrichment in clay and cement close to the interface (neoformation of up to 17 vol% Mg hydroxides in concrete, and up to 6 vol% in OPA within 0.6 mm at the interface), and the slightly increased S content in the cement 3–4 mm away from the interface qualitatively match the sample characterisation. Simulations of Portland cement—OPA interaction indicate a weaker chemical disturbance over a larger distance compared with low-pH cement—OPA. In the latter case, local changes in porosity are stronger and lead to predicted pore clogging.

Fault sealing and caprock integrity for CO&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt; storage: an in situ injection experiment

Abstract. The success of geological carbon storage depends on the assurance of permanent containment for injected carbon dioxide (CO2) in the storage formation at depth. One of the critical elements of the safekeeping of CO2 is the sealing capacity of the caprock overlying the storage formation despite faults and/or fractures, which may occur in it. In this work, we present an ongoing injection experiment performed in a fault hosted in clay at the Mont Terri underground rock laboratory (NW Switzerland). The experiment aims to improve our understanding of the main physical and chemical mechanisms controlling (i) the migration of CO2 through a fault damage zone, (ii) the interaction of the CO2 with the neighboring intact rock, and (iii) the impact of the injection on the transmissivity in the fault. To this end, we inject CO2-saturated saline water in the top of a 3 m thick fault in the Opalinus Clay, a clay formation that is a good analog of common caprock for CO2 storage at depth. The mobility of the CO2 within the fault is studied at the decameter scale by using a comprehensive monitoring system. Our experiment aims to close the knowledge gap between laboratory and reservoir scales. Therefore, an important aspect of the experiment is the decameter scale and the prolonged duration of observations over many months. We collect observations and data from a wide range of monitoring systems, such as a seismic network, pressure temperature and electrical conductivity sensors, fiber optics, extensometers, and an in situ mass spectrometer for dissolved gas monitoring. The observations are complemented by laboratory data on collected fluids and rock samples. Here we show the details of the experimental concept and installed instrumentation, as well as the first results of the preliminary characterization. An analysis of borehole logging allows for identifying potential hydraulic transmissive structures within the fault zone. A preliminary analysis of the injection tests helped estimate the transmissivity of such structures within the fault zone and the pressure required to mechanically open such features. The preliminary tests did not record any induced microseismic events. Active seismic tomography enabled sharp imaging the fault zone.

Assessing the hydraulic and mechanical impacts of heat generating radioactive waste at the whole repository scale

Radioactive waste disposal facilities for heat generating waste are designed to maintain a compressive stress regime in the host rock. Field-scale tests have been undertaken to study how heating of the host rock can reduce the effective stress on the rock. However, it is not possible to carry out experiments at the scale of a whole repository (km-scale), so variation of host rock properties that might be encountered at this scale is not captured by the experiments.Previous numerical modelling of field scale experiments has demonstrated that the physical processes leading to changes in the stress regime are well understood. Here those models are applied at the km-scale, considering approaches to representing a whole repository and the effect of variability in the properties of the host rock. A dataset from the Callovo-Oxfordian Claystone at the Meuse/Heute-Marne Underground Rock Laboratory is used to characterise the natural variability of a potential host rock at the scale of a disposal facility.The modelling demonstrates that for understanding the generation of thermal stresses, considerable learning can be gained from models that employ symmetry conditions and represent a small part of the repository, hence reducing the size of the computational problem. The host rock at the mid-point between waste cells remained in vertical compression for all the parameter combinations used, which builds confidence that horizontal fracturing of the host rock between waste cells due to thermal pressurisation is unlikely for the disposal concept discussed here. Ground surface uplift was also considered and is in the range of 9–13 cm for disposal system studied.Spatial variability in the properties of the host rock has the potential to lead to significant variability in temperature, pressure, effective stress and displacement around the facility. In particular, thermal conductivity, permeability, Young’s modulus and coefficient of thermal expansion need to be well characterised.

Empirical Investigations of the Instrument Response for Distributed Acoustic Sensing (DAS) across 17 Octaves

ABSTRACT With the potential of high temporal and spatial sampling and the capability of utilizing existing fiber-optic infrastructure, distributed acoustic sensing (DAS) is in the process of revolutionizing geophysical ground-motion measurements, especially in remote and urban areas, where conventional seismic networks may be difficult to deploy. Yet, for DAS to become an established method, we must ensure that accurate amplitude and phase information can be obtained. Furthermore, as DAS is spreading into many different application domains, we need to understand the extent to which the instrument response depends on the local environmental properties. Based on recent DAS response research, we present a general workflow to empirically quantify the quality of DAS measurements based on the transfer function between true ground motion and observed DAS waveforms. With a variety of DAS data and reference measurements, we adapt existing instrument-response workflows typically in the frequency band from 0.01 to 10 Hz to different experiments, with signal frequencies ranging from 1/3000 to 60 Hz. These experiments include earthquake recordings in an underground rock laboratory, hydraulic injection experiments in granite, active seismics in agricultural soil, and icequake recordings in snow on a glacier. The results show that the average standard deviations of both amplitude and phase responses within the analyzed frequency ranges are in the order of 4 dB and 0.167π radians, respectively, among all experiments. Possible explanations for variations in the instrument responses include the violation of the assumption of constant phase velocities within the workflow due to dispersion and incorrect ground-motion observations from reference measurements. The results encourage further integration of DAS-based strain measurements into methods that exploit complete waveforms and not merely travel times, such as full-waveform inversion. Ultimately, our developments are intended to provide a quantitative assessment of site- and frequency-dependent DAS data that may help establish best practices for upcoming DAS surveys.

Mechanical behavior of sandy facies of Opalinus Clay under different load conditions

Abstract The mechanical behavior of sandy facies of Opalinus Clay at the Mont Terri underground rock laboratory (URL) in Switzerland was investigated with drained and undrained triaxial compression and extension, cyclic compression, and creep tests. Samples were taken from boreholes drilled parallel to bedding. Most of the samples were reconditioned to minimize sampling effects of desaturation and micro-cracking. The compression was accomplished by increasing axial stress at constant radial stress. The extension was carried out by increasing radial stress at constant axial stress. Moreover, extension was also achieved by simultaneously increasing radial stress and decreasing axial stress under constant mean stress. The test results showed elastoplastic stress–strain behavior with volumetric compaction until onset of dilatancy at high deviatoric stresses above 80%–90% of the peak failure strength. The strength is dependent upon load path and mean stress. The strength under triaxial compression is higher than that under extension. The respective strength increases with increasing mean stress. Desaturation enhances the stiffness and strength of the claystone. The deformation and strength of the claystone are time-dependent. Under constant deviatoric stress, the claystone crept continuously with time, which can be characterized by a transient phase and a following stationary phase, and even a tertiary phase at high deviatoric stresses to rupture.

Middle‐Atmosphere Dynamics Observed With a Portable Muon Detector

AbstractIn the past years, large particle physics experiments have shown that muon rate variations detected in underground laboratories are sensitive to regional, middle‐atmosphere temperature variations. Potential applications include tracking short‐term atmosphere dynamics, such as Sudden Stratospheric Warmings. We report here that such sensitivity is not only limited to large surface detectors under high‐opacity conditions. We use a portable muon detector conceived for muon tomography for geophysical applications, and we study muon rate variations observed over 1 year of measurements at the Mont Terri Underground Rock Laboratory, Switzerland (opacity of ~700 meter water equivalent). We observe a direct correlation between middle‐atmosphere seasonal temperature variations and muon rate. Muon rate variations are also sensitive to the abnormal atmosphere heating in January–February 2017, associated to a Sudden Stratospheric Warming. Estimates of the effective temperature coefficient for our particular case agree with theoretical models and with those calculated from large neutrino experiments under comparable conditions. Thus, portable muon detectors may be useful to (1) study seasonal and short‐term middle‐atmosphere dynamics, especially in locations where data are lacking such as midlatitudes, and (2) improve the calibration of the effective temperature coefficient for different opacity conditions. Furthermore, we highlight the importance of assessing the impact of temperature on muon rate variations when considering geophysical applications. Depending on latitude and opacity conditions, this effect may be large enough to hide subsurface density variations due to changes in groundwater content and should therefore be removed from the time series.

Evolution of 10-20 year old Excavati-on Damaged Zones at the Mont Terri Underground Rock Laborato-ry and implications for self-sealing

Evolution of 10-20 year old Excavati-on Damaged Zones at the Mont Terri Underground Rock Laborato-ry and implications for self-sealing