Search

We use multitemporal analyses based on Synthetic Aperture Radar differential interferometry (DInSAR) to study the slope adjacent to the large Punatsangchhu-I hydropower plant, a concrete gravity dam under construction in Bhutan since 2009. Several slope failures affected the site since 2013, probably as a consequence of toe undercutting of a previously unrecognised active landslide. Our results indicate that downslope displacement, likely related to the natural instability, was already visible in 2007 on various sectors of the entire valley flank. Moreover, the area with active displacements impinging on the dam site has continuously increased in size since 2007 and into 2018, even though stabilization measures have been implemented since 2013. Stabilisation measures currently only focus on a small portion of the slope, however, the unstable area is larger than previously evaluated. Highly damaged rock is present across many areas of the entire valley flank, indicating that the volumes involved may be orders of magnitude higher than the area on which stabilisation efforts have been concentrated after the 2013 failure. The results highlight that satellite-based DInSAR could be systematically used to support decision making processes in the different phases of a complex hydropower project, from the feasibility study, to the dam site selection and construction phase. [ABSTRACT FROM AUTHOR]

Water infiltration into fractures is ubiquitous in crustal rocks. However, little is known about how such a progressive wetting process affects fracture stiffness and seismic wave propagation, which are highly relevant for characterizing fracture systems in situ. We study the acousto‐mechanical behavior of a free‐standing fractured granite subjected to gradual water infiltration with a downward‐moving wetting front over 12 days. We observe significant differences (i.e., by an order of magnitude) in wave amplitudes across the fractured granite compared to an intact granite, with both cases showing a strong correlation between wave amplitudes and wetting front movement. Effects of water infiltration into the fracture and surrounding matrix on seismic attenuation are captured by a numerical model with parameters constrained by experimental data. Back‐calculated fracture stiffness decreases exponentially with the wetting front migration along the fracture. We propose that moisture‐induced matrix expansion around the fracture increases asperity mismatch, leading to reduced fracture stiffness. Plain Language Summary: In the shallow layers of the Earth, hydrological cycles such as snowmelt, fog, dew, and rain have been shown to change the moisture content of crustal rocks, which can alter the elastic properties of natural fractures and affect the propagation of seismic waves. Understanding how seismic waves propagate in the near‐surface environment is crucial for the assessment of earthquake hazards and the characterization of geologic heterogeneities. In this work, we perform well‐controlled laboratory experiments to study the acousto‐mechanical behavior of a single fracture in granitic rock subjected to progressive wetting over 12 days. We report that the fracture stiffness decreases exponentially as the wetting front advances along the fracture. Our research sheds light on an important question in fracture characterization: how elastic waves propagate across a fracture undergoing moisture‐induced expansion. Key Points: A laboratory study establishes a relationship among water imbibition, seismic attenuation, and stiffness evolution in a wetted fractureWave amplitudes across a fracture correlate strongly with the wetting front movement of infiltrated water within the fractureFracture stiffness exponentially decreases with the advance of the wetting front along the fracture [ABSTRACT FROM AUTHOR]

The paper represents a summary about our research projects conducted between 2003 and 2015 related to the mechanical behaviour of Opalinus Clay at Mont Terri. The research summarized covers a series of laboratory and field tests that address the brittle failure behaviour of Opalinus Clay, its undrained and effective strength, the dependency of petro-physical and mechanical properties on total suction, hydro-mechanically coupled phenomena and the development of a damage zone around excavations. On the laboratory scale, even simple laboratory tests are difficult to interpret and uncertainties remain regarding the representativeness of the results. We show that suction may develop rapidly after core extraction and substantially modifies the strength, stiffness, and petro-physical properties of Opalinus Clay. Consolidated undrained tests performed on fully saturated specimens revealed a relatively small true cohesion and confirmed the strong hydro-mechanically coupled behaviour of this material. Strong hydro-mechanically coupled processes may explain the stability of cores and tunnel excavations in the short term. Pore-pressure effects may cause effective stress states that favour stability in the short term but may cause longer-term deformations and damage as the pore-pressure dissipates. In-situ observations show that macroscopic fracturing is strongly influenced by bedding planes and faults planes. In tunnel sections where opening or shearing along bedding planes or faults planes is kinematically free, the induced fracture type is strongly dependent on the fault plane frequency and orientation. A transition from extensional macroscopic failure to shearing can be observed with increasing fault plane frequency. In zones around the excavation where bedding plane shearing/shearing along tectonic fault planes is kinematically restrained, primary extensional type fractures develop. In addition, heterogeneities such as single tectonic fault planes or fault zones substantially modify the stress redistribution and thus control zones around the excavation where new fractures may form. [ABSTRACT FROM AUTHOR]

Catastrophic landslides characterized by runaway slope failures remain difficult to predict. Here, we develop a physics‐based framework to prospectively assess slope failure potential. Our method builds upon the physics of extreme events in natural systems: the extremes so‐called "dragon‐kings" (e.g., slope tertiary creeps prior to failure) exhibit statistically different properties than other smaller‐sized events (e.g., slope secondary creeps). We develop statistical tools to detect the emergence of dragon‐kings during landslide evolution, with the secondary‐to‐tertiary creep transition quantitatively captured. We construct a phase diagram characterizing the detectability of dragon‐kings against "black‐swans" and informing on whether the slope evolves toward a catastrophic or slow landslide. We test our method on synthetic and real data sets, demonstrating how it might have been used to forecast three representative historical landslides. Our method can in principle considerably reduce the number of false alarms and identify with high confidence the presence of true hazards of catastrophic landslides. Plain Language Summary: Catastrophic slope failures that pose great threats to life and property remain difficult to predict due to the strong variability of slope behavior. As a result, only a limited number of large rock slope failures have been so far successfully forecasted with associated risks mitigated. Here, we propose a novel predictive framework to prospectively and quantitatively detect slope failure precursors with high confidence. Our research sheds light on one of the most challenging questions in landslide prediction: Would an active landslide slowly move or catastrophically fail in the future? Our method adds a new conceptual framework and operational methodology with a significant potential to support existing early warning systems and hence reduce landslide risks. Key Points: Tertiary creeps of catastrophic landslides accommodate dragon‐kings showing statistically different properties than secondary slope creepsA predictive framework is developed to forecast catastrophic landslides by detecting signatures typical of the emergence of dragon‐kingsA phase diagram characterizes the detectability of dragon‐kings against black‐swans and discriminates catastrophic and slow landslides [ABSTRACT FROM AUTHOR]

[Display omitted] • We explore THM coupled processes in a new TBM tunnel in the Opalinus Clay shale. • We present the evolution of in-situ state variables acquired over four years. • Radial total pressures (RTP) on the tunnel support are heterogeneous. • RTP reflect construction activities and long-term tunnel climate variations. • RTP reveal swelling, gap grout yielding, and are consistent with lab swelling tests. The old Belchen tunnel tubes in the Swiss Jura Mountains were excavated with drill-and-blast in swellable sedimentary rocks, i.e., anhydrite-rich marls (Gipskeuper) and Opalinus Clay shale (OPA). Already during construction in the 1960s both rock formations caused substantial damage to the tunnel support through high swelling pressures and heave, and in later years the tubes had to be refurbished again. Important maintenance and repair prompted the construction of a new, third Belchen tunnel tube (2016–2021) with a tunnel-boring machine (TBM). In this study we present in-situ datasets acquired to investigate the stress evolution and controlling mechanisms over more than four years at a monitoring section located in a strongly faulted OPA section of the new Belchen tunnel tube. The main datasets include time series of total radial pressure, radial strain, rock water content, rock and concrete temperatures, as well as details of the geological structures obtained from analyses of borehole logs and three-dimensional photogrammetric excavation face models. Finally, a series of idealised numerical simulations explore the impact of measured temperature variations on the measured total pressures, which confirm a strong temperature effect on radial pressures related to the setting of concrete and seasonal climatic variations. We find that in our monitoring section radial pressures on the tunnel support are very heterogeneous, i.e., they range between 0.5 MPa and 1.5 MPa, and still gently increasing 4 years after excavation. The measured pressures are 2–5 times greater than measured in the old Belchen tunnel tubes and similar in magnitude to swelling pressures obtained in laboratory tests. EDZ permeability measurements, water content evolution, and radial strain data from the tunnel invert suggest that swelling processes contribute to the long term radial pressure build-up. Thermo-elastic deformation and swelling might be superimposed by local reactivation of tectonic faults and gap grout cracking at crack-initiation stress levels. [ABSTRACT FROM AUTHOR]

Correction to: I Scientific Reports i https://doi.org/10.1038/s41598-021-89015-8, published online 06 May 2021 The original version of this Article contained an error in Table 1, where the formula given for the MCC score was incorrect. [Extracted from the article]

We present a combined experimental and numerical study of failure and damage in faulted Opalinus Clay shale around an opening at the Mont Terri Rock Laboratory, Switzerland. An experiment borehole with a diameter of 0.6 m and a length of 12.9 m intersecting a major fault zone at an acute angle of ∼ 40 ° was drilled, with subsequent overbreaks closely monitored. To investigate the underlying mechanisms that lead to the observed overbreaks, we developed a 3D geomechanics model to simulate the deformation and failure behaviour of faulted Opalinus Clay shale. To represent geological structures including the major fault zone and secondary fracture sets, a site-specific fracture network was constructed. Our model captured many important geomechanical properties and responses of the faulted Opalinus Clay shale, such as anisotropy of the shale matrix, deformation of the fault zone, dislocation of secondary fractures, and growth of new cracks as well as generation of overbreaks around the opening. We compared our simulation results with in-situ experimental observation of short-term overbreak patterns along the borehole and found that the overbreak occurrence is strongly controlled by both stress conditions and geological structures. Our results indicate that the damaged zone is characterised by an inner shell and an outer shell, where the former (dominated by extensive tensile/shear cracks) has a thickness of about one quarter of the borehole diameter and the latter (dominated by sparse shear cracks) extends by half to one borehole diameter into the surrounding rock. We further elucidated the impacts of the major fault zone and secondary fracture sets as well as the borehole–fault intersection angle on the deformation, damage, and failure characteristics. The findings and insights obtained in this study have important implications for the stability evaluation of waste emplacement drifts and long-term safety assessments of nuclear waste galleries in faulted argillaceous rocks. • 3D simulation captures overbreak and damage of Opalinus Clay around an opening. • Interplay of stress and heterogeneity drives overbreak formation and pattern. • Damaged zone of a two-shell structure rotates with the distance to a fault zone. • Acute angled faults exert strong controls on overbreak and damage generation. [ABSTRACT FROM AUTHOR]

Abstract: In this study, the rock mass response ahead of an advancing test tunnel in the Opalinus Clay at the Mont Terri Rock Laboratory (Switzerland) was investigated. Characterisation of the excavation-induced damage zone at Mont Terri is a challenging task due to the anisotropic and heterogeneous nature of the shale: pronounced bedding leads to intact rock anisotropy and prevalent small-scale tectonic shears lead to rock mass heterogeneity. Rock mass damage ahead of an experimental tunnel or niche was characterised through single-hole seismic wave velocity logging, borehole digital optical televiewer imaging, and geological drillcore mapping. Three-dimensional elastic stress analyses were completed and showed that rock mass degradation can be correlated to changes in the maximum to minimum principal stress ratio (i.e., spalling limit). Numerical results showed that close to the niche boundary, unloading lowers stress ratios, which correspond with decreasing seismic wave amplitudes and velocities; thus, indicating that strength degradation resulted from increasing crack-induced damage. Considerations of tectonic shears and distance from a previously stressed volume of rock were necessary in understanding both the damage state and extent ahead of the face. By integrating field and numerical data, the investigation showed that geological structures (i.e., bedding and bedding-parallel tectonic shears) were most influential near the entrance but played a lesser role as the niche deepened. Additionally, a portion of the niche is located in the perturbed zone of the intersecting Gallery04. [Copyright &y& Elsevier]

Abstract: The Opalinus Clay is currently under investigation as a potential host rock for deep geological disposal of nuclear waste at the Mont Terri Rock Laboratory in Switzerland. Bedding in the Opalinus Clay at Mont Terri is ubiquitous and highly persistent leading to mechanical transverse isotropy. Adding to the complexity at the Rock Laboratory is the frequent occurrence of small-scale tectonic shears. This paper explores the influences of millimetre-thick tectonic shears and bedding on the development of induced fractures mapped in the EZ-B field experiment at the research facility. A series of numerical analyses were carried out by increasing the geological complexity of the host rock and comparing the redistributed stress field with geological maps of the induced fractures. The analyses show that if tectonic shears are not kinematically constrained, mobilisation of the shears can play a key role in the development of the induced fracture network and therefore, be a primary factor in the development of the excavation damaged zone. This illustrates that under certain conditions rock mass heterogeneity (in this case, resulting from the tectonic shears) may dominate over rock matrix anisotropy (in this case, resulting from bedding) and must be considered when predicting the induced fracture network of the excavation damaged zone. [ABSTRACT FROM AUTHOR]

In this work we present a methodology based on remote sensing data to map and classify unstable slopes in high alpine areas and apply this methodology to the Himalaya of northwestern Bhutan. We leverage on radar images acquired by the ESA Envisat and the JAXA ALOS-1 satellites between 2007 and 2011. Synthetic aperture radar differential interferometry (DInSAR) is used to identify and quantify recent ground surface displacements associated with potentially unstable slopes. Optical satellite images and a high-resolution digital surface model is used for the geomorphological classification of individual landforms and deposits associated with displacements. The analysis of the geomorphological characteristics reveals that displacements are associated with: rock slides, mountain slope deformation and rock slope deformation (77.5%), soil creep (4%), soil slides (1.9%) and rock glaciers (12.3%); 4.3% of the detected displacements remain of unknown origin. The method applied here for the assessment of activity likelihood uses a very large number of interferograms obtained with images acquired by different satellites. The method, which is new and replicable, takes advantage of the combination of different acquisition geometries and then combines indices of vegetation cover, layover, velocity estimate and temporal sampling associated with each detected active area to assign a weight. The analysis shows that a cluster of unstable slopes exists in the northwest, potentially related to the presence of an important regional lineament and the distribution of permafrost. Comparatively little ongoing activity is detected in the south, but field validation documents substantial landslide activity in the recent geological past. The dataset here generated, containing 700 potentially unstable areas with an associated activity likelihood, is a novel dataset that makes an essential basis of future analysis to investigate rock slope instability distribution and predisposing factors. Moreover, this dataset represents a new important step towards future landslide hazard assessment in the Himalaya of Bhutan. • 693 previously unknown unstable slopes in NW Bhutan were found with DInSAR • Geomorphological classification revealed that almost 78% are rock slope instabilities • The spatial distribution of potentially unstable slopes in the region is clearly non uniform • Activity of unstable slopes is likely influenced by regional tectonic structures • Unstable slopes seem to be less abundant in the region in recent times than in the past [ABSTRACT FROM AUTHOR]

Joint set spacing is a fundamental parameter in the determination of rock mass quality and can be measured in situ by means of, e.g., scanline surveys and/or by remote sensing techniques, such as photogrammetric analyses. In many Alpine areas, rock mass outcrops are not easily accessible, and geomechanical propexrties can be measured only in a few unevenly distributed locations, which are often separated by large distances (in the order of hundreds of meters). Geostatistical techniques have been explored to achieve a reliable estimate of rock mass properties in unreachable zones. This work aims to estimate joint set spacing of outcropping rock masses and the associated uncertainty, using photogrammetric models and geostatistical modeling, in an area of about 30 km2, located in the Aar Massif of the Swiss Central Alps. Since the joint set spacing is strongly related to fracture genesis, joints were subdivided according to their type and age into three younger exfoliation joint sets and older tectonic joint sets. Each fracture set spacing was analyzed by variography and its spatial distribution was estimated using Sequential Gaussian Simulations. The present research proves that photogrammetric techniques combined with geostatistical modeling can be satisfactorily applied to develop predictive maps of joint spacing, providing that the geological processes governing the formation of joints are taken into account. Moreover, predictive maps must be associated with the corresponding uncertainty maps. Finally, in mountainous area, whenever a full three-dimensional approach is not feasible due to the absence of subsurface data, at least the ground altitude, besides latitude and longitude, has to be included in the geostatistical modeling. [ABSTRACT FROM AUTHOR]

The water adsorption into pore spaces in brittle rocks affects wave velocity and transmitted amplitude of elastic waves. Experimental and theoretical studies have been performed to characterize moisture-induced elastodynamic variations due to macroporous effects; however, little attention has been paid to the manner in which wetting of nanopores affects elastic wave transmission. In this work, we extend our understanding of moisture-induced elastic changes in a microcracked nanopore-dominated medium where 80 % of the surface area exhibits pore diameters (also include microcrack widths) below 10 nm. We study acousto-mechanical response resulting from a gradual wetting on a free-standing intact Herrnholz granite specimen over 98 h using time-lapse ultrasonic and digital imaging techniques. Linkages between ultrasonic attributes and adsorption-induced stress/strain are established during the approach of the wetting front. We find that Gassmann theory, previously validated in channel-like nanoporous media, does not work properly in predicting the P-wave velocity increase of microcracked nanopore-dominated media at ultrasonic frequency. However, squirt flow – a theory recognized to characterize wave velocity increase and attenuation in microcracked macropore-dominated media at the pore scale – also accounts for the observed increase of P-wave velocity in microcracked nanopore-dominated media. The transmitted amplitude changes in direct P waves are explained and predicted by the elastic wave propagation within P-wave first Fresnel zone and reflection/refraction on the wetting front. • Acousto-mechanical response of a free-standing intact specimen to gradual wetting. • Squirt flow dominates P-wave velocity increase in microcracked nanoporous media. • Amplitude changes are explained by elastic wave propagation around Fresnel zone. [ABSTRACT FROM AUTHOR]

The Gotthard Base Tunnel (GBT), constructed between 2000 and 2011, is a 57 km long and up to 2.5 km deep high-speed railway tunnel located in the Swiss Alps. Significant ground surface displacements reaching about 10 cm were observed during and after the tunnel construction. To gain a better understanding of the causal mechanism of such conspicuous ground deformations, we develop a three-dimensional (3D) fully-coupled hydro-mechanical model to simulate GBT-induced groundwater drainage, stress redistribution, rock mass consolidation, and fault zone deformation at the regional scale. First, we construct a geological model with topographical features, lithological units, and natural faults realistically represented. We constrain the material properties of fault zones and rock masses based on available extensive laboratory testing results and site characterisation datasets. We then simulate the tunnelling process over time, with the resulting coupled hydro-mechanical responses of faulted rock masses well captured and ground surface/subsurface displacements quantitatively analysed. The simulation results in general show a good agreement with the field monitoring data of ground surface displacement, subsurface tunnel settlement, and groundwater inflow into the tunnel. Our model indicates that ground surface displacements originate from GBT-induced water drainage and rock mass consolidation in the deep subsurface. Our results also show that the GBT construction could trigger faults to shear via drainage-induced pressure diffusion and poroelastic stressing. The research findings from our work have important implications for many groundwater drainage-related geoengineering activities such as underground excavation in alpine mountains and fluid withdrawal in subsurface reservoirs. [Display omitted] • 3D simulation captures tunnelling-induced groundwater flow and rock mass deformation. • The simulation explains ground subsidence in crystalline rock during deep tunnelling. • Drainage and consolidation in deep subsurface cause ground surface displacements. • Tunnelling-induced pressure diffusion and poroelastic stressing drive fault movement. [ABSTRACT FROM AUTHOR]

Valley glaciers have traditionally been expected to significantly influence the stability and movement rates of adjacent paraglacial landslides. However, detailed studies related to the mechanical and displacement interactions between glacier ice and unstable rock slopes are very rare. Here we present a detailed in-situ investigation of the spatial variations of the displacement field of the Great Aletsch Glacier in the surroundings of a large active instability, i.e., the Moosfluh Landslide, a Deep-Seated Gravitational Slope Deformation, with superimposed large (1–5 million m3) secondary rockslides formed during fall 2016. We performed repeat UAV-based photogrammetric surveys during ~3 days (74 h) and applied Digital Image Correlation techniques to record high-resolution surface displacement vector fields of the landslide, stable slopes and adjacent glacier. Our results show that the secondary rockslide adjacent to the glacier is composed of two parts of 1.5 and 2.8 million m3 volume respectively, showing significant differences in mean displacement velocities (0.4 and 0.9 m in 74 h respectively, excluding rapid movements from detached blocks). Both rockslide compartments induce clear deflections of the glacier flow field, moving with a maximum velocity of about 0.3 to 0.4 m in 74 h. This influence is highest near the ice-contact boundary and decreases within a distance of about 100–200 m from the rock slope instability. We investigate the viscous forces at the landslide/glacier contact using a straightforward analytical model for an incompressible rockslide block sliding along a planar, cohesionless surface into ice. These forces are then applied to a slope stability model based on the limit equilibrium concept, representing the real geometry at the interface boundary to quantitatively explore the true buttressing effects of valley glaciers on an already moving slope instability. We show that the viscous ice deformation plays an important role in mediating the displacement velocities of landslides in unstable conditions, while, on the other hand, the slope support from the valley glacier has very little influence on the stability of the investigated rockslides. As most valley glaciers are currently strongly retreating due to global warming, uncovering significant numbers of pre-LIA slope instabilities, this detailed investigation provides important hints on their potential displacement behavior. Understanding the factors controlling landslide velocity is of great importance for hazard analyses and early warning. Hence, this study has implications beyond academic interest, e.g., for the planning and operation of alpine infrastructure, such as cable cars or hydropower systems. • We monitor glacier displacement fields at contacts to stable slopes and an active rockslide. • Distortions of glacier displacements fields extend over 100-200 m. and correlate with variations in rockslide velocity. • A simplified analytical model with viscous ice forces can reproduce observed rockslide velocities. • Limit equilibrium analyses demonstrate negligible effect of glacial debuttressing on mature rockslide stability. [ABSTRACT FROM AUTHOR]

Slope deformation in mountainous terrain can be driven by different processes, the nature of which is either gravitational and irreversible or seasonal and reversible, the latter induced by permafrost variations or by hydromechanical or thermomechanical effects. The importance of identifying such deformations is not only related to the hazard they can pose, but also to the understanding of changes that permafrost or local hydrological conditions undergo. Here, we carry out a multi-temporal InSAR analysis over a mountainous area 8000 km2 large, straddling north-western Bhutan and southern Tibet. We propose a methodology to separate locally deforming areas from the effects of long-wavelength signals through the analysis of the spatio-temporal behaviour of 4-years long time series of surface displacements. We present the mapping of hundreds of small-scale features that appear to be actively deforming, as well as several examples of reversible deformation rarely detected at this scale in such a challenging and vast region. The analysis of the multi-annual trend of ground deformation shows a relatively small number of irreversible gravitational movements clearly related to rock slides which attests a low level of recent activity of large rock slope instabilities in the region. In the southernmost, lower elevation parts of the study area, we quantify reversible surface displacements with amplitudes ranging between 5 and 17 mm, and showing maximum displacements towards the satellite in summer months, thus compatible with hydro-mechanical effects related to groundwater table variations. In addition, the ground movement induced by the active layer's response to thawing and freezing over the gentler slopes and high-elevation permafrost regions of Bhutan and southern Tibet is on average around 10 mm, with maxima up to 28 mm. The localised displacements appear to be largely associated to braided stream plains, glacier outwash plains or low angle, fine sediment covered slopes. • Inventory of gravitational and reversible surface displacements for NW Bhutan • Classification of slope processes made on the basis of time series trends • Reversible, mostly vertical deformation, observed in the permafrost region • Reversible deformation with horizontal component observed at lower elevations • Rock glaciers appear to have lower creep rates than in the Alps. [ABSTRACT FROM AUTHOR]

Past exploration missions have revealed that the lunar topography is eroded through mass wasting processes such as rockfalls and other types of landslides, similar to Earth. We have analyzed an archive of more than 2 million high-resolution images using an AI and big data-driven approach and created the first global map of 136.610 lunar rockfall events. Using this map, we show that mass wasting is primarily driven by impacts and impact-induced fracture networks. We further identify a large number of currently unknown rockfall clusters, potentially revealing regions of recent seismic activity. Our observations show that the oldest, pre-Nectarian topography still hosts rockfalls, indicating that its erosion has been active throughout the late Copernican age and likely continues today. Our findings have important implications for the estimation of the Moon's erosional state and other airless bodies as well as for the understanding of the topographic evolution of planetary surfaces in general.

High-resolution 3D geological models are crucial for underground development projects and corresponding numerical simulations with applications in e.g., tunneling, hydrocarbon exploration, geothermal exploitation and mining. Most geological models are based on sparse geological data sampled pointwise or along lines (e.g., boreholes), leading to oversimplified model geometries. In the framework of a hydraulic stimulation experiment in crystalline rock at the Grimsel Test Site, we collected geological data in 15 boreholes using a variety of methods to characterize a decameter-scale rock volume. The experiment aims to identify and understand relevant thermo-hydro-mechanical-seismic coupled rock mass responses during high-pressure fluid injections. Prior to fluid injections, we characterized the rock mass using geological, hydraulic and geophysical prospecting. The combination of methods allowed for compilation of a deterministic 3D geological analog that includes five shear zones, fracture density information and fracture locations. The model may serve as a decameter-scale analog of crystalline basement rocks, which are often targeted for enhanced geothermal systems. In this contribution, we summarize the geological data and the resulting geological interpretation.