Vol.:(0123456789)1 3 Environmental Earth Sciences (2021) 80:509 https://doi.org/10.1007/s12665-021-09787-0 THEMATIC ISSUE GeomInt: geomechanical integrity of host and barrier rocks–experiments, models and analysis of discontinuities Olaf Kolditz1,7 · Thomas Fischer1 · Thomas Frühwirt2 · Uwe‑Jens Görke1 · Carolin Helbig1 · Heinz Konietzky2 · Jobst Maßmann3 · Mathias Nest4 · Daniel Pötschke2 · Karsten Rink1 · Amir Sattari5 · Patrick Schmidt6 · Holger Steeb6 · Frank Wuttke5 · Keita Yoshioka1 · Bernhard Vowinckel3,8 · Gesa Ziefle3 · Thomas Nagel2 Received: 19 June 2020 / Accepted: 13 July 2021 / Published online: 7 August 2021 © The Author(s) 2021 Abstract The present paper gives an overview of the GeomInt project “Geomechanical integrity of host and barrier rocks—experi- ment, modelling and analysis of discontinuities” which has been conducted from 2017–2020 within the framework of the “Geo:N Geosciences for Sustainability” program. The research concept of the collaborative project is briefly introduced followed by a summary of the most important outcomes. The research concept puts geological discontinuities into the centre of investigations—as these belong to the most interesting and critical elements for any subsurface utilisation. Thus, while research questions are specific, they bear relevance to a wide range of applications. The specific research is thus integrated into a generic concept in order to make the results more generally applicable and transferable. The generic part includes a variety of conceptual approaches and their numerical realisations for describing the evolution of discontinuities in the most important types of barrier rocks. An explicit validation concept for the generic framework was developed and realised by specific “model-experiment-exercises” (MEX) which combined experiments and models in a systematic way from the very beginning. 16 MEX have been developed which cover a wide range of fundamental fracturing mechanisms, i.e. swelling/ shrinkage, fluid percolation, and stress redistribution processes. The progress in model development is also demonstrated by field-scale applications, e.g. in the analysis and design of experiments in underground research laboratories in Opalinus Clay (URL Mont Terri, Switzerland) and salt rock (research mine Springen, Germany). Keywords GeomInt · Fracture flow · Fracture mechanics · Barrier integrity · Discontinuities · Open source · OpenGeoSys Introduction The use of the subsurface as a source of resources, a stor- age space and for installing underground municipal or traffic infrastructure has increased both in intensity and diversity in recent years Volchko et al. (2020). In addition to clas- sical anthropogenic interventions such as mining, oil and gas production or tunnel construction, other forms of under- ground use have come into the focus of economic, politi- cal and scientific activities, particularly in connection with the transformation of energy systems and a global trend to urbanization McCartney et al. (2016). These include, for example, the extraction of energy (e.g. geothermal energy) and energy sources using new technologies (e.g. unconven- tional gas extraction) Field et al. (2018) as well as the geo- logical short-term and long-term storage of energy carriers such as compressed air, hydrogen, and methane Kabuth et al. (2017); Bauer et al. (2012) and the safe storage of waste This article is a part of the Topical Collection in Environmental Earth Sciences on “Sustainable Utilization of Geosystems” guest edited by Ulf Hünken, Peter Dietrich and Olaf Kolditz. * Olaf Kolditz olaf.kolditz@ufz.de 1 Helmholtz Centre for Environmental Research - UFZ, Leipzig, Germany 2 Technische Universität Bergakademie Freiberg - TUBAF, Freiberg, Germany 3 Federal Institute for Geosciences and Natural Resources - BGR, Hanover, Germany 4 Institut für Gebirgsmechanik GmbH - IfG, Leipzig, Germany 5 Christian-Albrechts-Universität zu Kiel - CAU , Kiel, Germany 6 University of Stuttgart, Stuttgart, Germany 7 Technische Universität Dresden, Dresden, Germany 8 Technische Universität Braunschweig, Braunschweig, Germany http://crossmark.crossref.org/dialog/?doi=10.1007/s12665-021-09787-0&domain=pdf Environmental Earth Sciences (2021) 80:509 1 3 509 Page 2 of 20 generated during energy production or in other industrial sectors, e.g. carbon dioxide, radioactive waste, or chemo- toxic waste Hudson et al. (2001); Gens (2010); Birkholzer et al. (2019); Ringrose and Meckel (2019). The rocks and rock mass themselves are characterised by complex material behaviour. Irreversible deformation, rate-dependence, swelling and shrinkage effects occur under real loading conditions. Damping aspects, for example in the range of low-frequency seismic waves, damage and physical- chemical ageing (e.g. solution and/or precipitation reactions) play an important role in many geotechnical applications and have their roots in the multiphase nature of geomate- rials Buscarnera (2012); Steeb et al. (2014); Loret et al. (2002). One common theme in the deformation behaviour of geomaterials is localization leading to structural features such as shear bands, compaction bands, unloading fractures. The formation, presence and evolution of these features can quickly dominate the system evolution and is implicated in most geological processes such as tectonics, orogenesis, hydrothermal systems, etc. Hobbs et al. (1990); Jirásek (2007); Duretz et al. (2018); Döhmann et al. (2019); Parisio et al. (2019). Despite their importance, the topics of dam- age, crack formation and propagation as well as interface problems are aspects that require further clarification from a scientific perspective. Because of prevailing knowledge gaps in these areas, many such processes can currently not be adequately modelled with numerical simulation systems established in geotechnical practice and therefore represent topics in urgent need of research. The corresponding pro- cesses manifest themselves in diverse micro- and macro- mechanical phenomena in the considered materials and are here summarised under the umbrella term discontinuities (Fig. 1). Geological systems in their natural state feature a range of discontinuities which, from a modelling perspective, can be loosely categorised into weak discontinuities, such as material interfaces, and strong discontinuities, such as frac- tures Armero and Garikipati (1996); Hansbo and Hansbo (2004). These existing discontinuities determine the in-situ stress-state as well as hydrothermal flow fields encountered when attempting to establish a particular geotechnology at a given site. Such anthropogenic intervention locally disturbs the quasi-steady state of the system and shifts it further away from equilibrium, thus triggering a range of multi-physical processes. These thermo-hydro-mechanical- chemical-microbiological (THMCB) processes can in turn trigger the formation of new discontinuities in the form of stress fractures, hydraulic fractures, quenching fractures or chemical wormholes Regenauer-Lieb et al. (2013); Liu et al. (2019). From the perspectives of operational and environ- mental safety, several consequences are of interest, such as induced seismicity due to the formation of new fractures or Fig. 1 Graphical abstract of the GeomInt project: geological reservoir-barrier systems and geomechanical integrity. The formation of new dis- continuities and their connectivity to existing fracture networks can dominate system evolution and needs to be thoroughly understood Environmental Earth Sciences (2021) 80:509 1 3 Page 3 of 20 509 the reactivation of existing faults Guglielmi et al. (2015); Ucar et al. (2017); Chang and Segall (2016); Parisio et al. (2019). Very often the creation of new fluid pathways and of connectivity to existing fractures or fracture networks is of primary importance: it is desired in fossil or geothermal resource extraction and to be avoided in hazardous waste disposal applications Quintal et al. (2014); Armand et al. (2014); Lei et al. (2017). However, in all applications the formation of new discontinuities has to occur in a controlled manner within a data-poor environment. Therefore, a sound theoretical understanding of the evolution of discontinuities is required for virtually all geotechnologies (Fig. 1). The aim of the GeomInt project is to contribute to this knowledge base and its transfer to practical engineering. GeomInt is dealing with specific research questions con- cerning the evolution of discontinuities: – Swelling and shrinkage processes in clay rocks are related to the key question of barrier integrity under changing moisture conditions. To what degree can the low permeability of clay rock be maintained under desat- uration-induced shrinkage during the heating process and can it be fully re-established by swelling processes upon subsequent resaturation? (see work package WP1) – Fluid percolation processes have been studied in both clay and salt rocks. How does fluid percolation depend on the surrounding stress regimes? How can the related threshold be defined for both rock types? How can the evolution of the discontinuities be described effectively in the presence of initial stress effects? (see work package WP2) – Stress redistribution is a fundamental mechanism of frac- ture initiation in brittle rocks. How will cyclic loading of rock fractures influence the evolution of fracture perme- ability? (see work package WP3) but also more general questions concerning the transferabil- ity of experimental and numerical methods for various rock types: – How to determine anisotropic and lamination effects in different clay rocks, i.e. from sandy and clayey facies (see “Experimental platform”)? – How to compare and combine different conceptual (con- tinuum and discontinuum) and numerical approaches (lattice, finite and distinct elements) for accurate descrip- tion of discontinuity evolution in various brittle and duc- tile rock types (see “Numerical platform”)? – How to combine experimental and modelling approaches in an optimal way to build cross-scale approaches for the analyses of discontinuity propagation (see Model- Experiment (MEX) concept)? Methodology In order to address the above posed research questions, Geo- mInt was following an integrated approach. The understand- ing and quantification of interactions with transient geologi- cal rock properties (e.g. permeability), which determine the geomechanical integrity and tightness of geological reser- voir-barrier systems, are at the centre of this work. Included in the investigations are discontinuities conceptualised as volumetrically distributed damage types as often used for describing the damaged zone of solid rocks, or as discrete crack and fracture networks. Three physical effects leading to the emergence and development of specific discontinuities are considered as main research areas: Swelling and shrinkage processes, pres- sure-driven fluid percolation and stress redistribution—and are forming the structure of the GeomInt project (Fig. 2). WP1 is dealing with swelling and shrinkage processes in clay stone caused by saturation changes. These processes are of interest mainly for understanding the evolution of the rock properties around ventilated tunnels as well as under elevated temperature conditions as may occur around emplacements of canisters containing heat-emitting radio- active waste Gens et al. (2010); Bossart et al. (2017). One of the questions is to what degree the low permeability of clay rock can be maintained in lieu of desaturation-induced shrinkage processes and whether potential permeability increases are reversible upon resaturation of the clay rock. BGR obtained Opalinus Clay samples from the Mont Terri Underground Research Lab (URL) in Switzerland and pro- vide them for laboratory testing to the partners CAU and IfG in Kiel and Leipzig, respectively. Samples come from various clayey and sandy facies. Experimental work with specimens from the sandy facies is particularly challenging due to the significant heterogeneity of the material. WP1 is closely linked with the Mont Terri Project in cooperation with swisstopo Bossart and Milnes (2018). Fluid percolation processes are studied in both clay and salt rocks, i.e. brittle-to-ductile materials. Samples of Opa- linus Clay came from the Mont Terri URL (see above). Salt rock samples are mainly obtained from the Springen site in Thuringia, Germany. Experimental work concerning perco- lation is conducted in the IfG and CAU labs in Leipzig and Kiel. WP2 investigates the existence and nature of a perco- lation threshold for both rock types depending on hydro- mechanical (HM) processes (i.e. fluid pressure and mechani- cal stress field). The development of discontinuities due to fluids moving through a newly formed percolation network in previously tight rock is of particular relevance for applica- tions of fluid storage in the subsurface, thermally induced fluid expansion and migration around heat sources, as well as intended hydraulic fracturing of shale rocks Minkley et al. Environmental Earth Sciences (2021) 80:509 1 3 509 Page 4 of 20 (2015); Ghanbarzadeh et al. (2015); Harrington et al. (2017); Armand et al. (2017); Yoshioka et al. (2019). WP3 is investigating discontinuities formed by and sub- jected to stress redistribution in brittle materials. Stress redistribution constitutes a fundamental fracturing mecha- nism that is also closely linked to the processes studied in WPs 1 and 2, where the multiphysical couplings are stronger McDermott et al. (2018); Parisio et al. (2019). Granite rock samples are obtained from locations in the Ore Mountains, i.e. from Kirchberg and Freiberg (URL Reiche Zeche). Experimental investigations are conducted in the Freiberg (TU Freiberg) and Stuttgart labs (University of Stuttgart). Constant Normal Load (CNL) and Constant Normal Stiff- ness (CNS) experiments are conducted to study fluid flow in rough fractures under confining stresses. Rock samples from Freiberg will be also used within the “Crystalline Task” of the new DECOVALEX-2023 phase. For building strong links between the above described work packages, laboratory work (LAB), modelling (MOD) and in-situ experiments (URL for Underground Research Laboratories) are clustered (Fig. 2 upper right). Experimen- tal and modelling platforms have been established during project in order to foster collaboration between the rock type oriented groups (see sections “Experimental platform” and “Numerical platform”, respectively). The clear structure of the GeomInt project is part of its specific research method- ology. The process-oriented work packages are interlinked with synthesis activities such as data and model integration using virtual reality (VR) methods (Fig. 2 lower right). A corresponding pilot demonstrator is being implemented in the Mont Terri project Rink et al. (2012). GeomInt relies on a very close link between experimental and modelling works. To this purpose 16 Model-Experi- ment-Exercises (MEX) have been defined at the beginning of the project (see section “Model-Experiment-Exercises (MEX)”). Within these exercises various numerical methods have been applied to same experimental data sets in order to determine advantages / disadvantages of particular numeri- cal methods (see section “In-situ experiments (URL)”) and e.g. for experimental design of the CD-A experiment in Mont Terri (see section “Synthesis–Numerical methods capabilities”). To link the laboratory experiments to the in-situ experi- ments and related applications of current interest, samples for the lab work are taken from the URLs. Figure 3 illus- trates the geographical WP workflows from in-situ sam- pling to geomechanical laboratories and modelling. The main sources for rock samples are (i) Mont Terri for clay (ii) Springen for salt and (iii) Freiberg/Kirchberg for crystal- line rock specimens. Moreover, there is collaboration with other URLs (Bure, Grimsel) concerning experimental and modelling work—mainly for testing transferability of the methodology to other rock types (e.g. Callovo-Oxfordian Clay—COx). Fig. 2 Work package structure of GeomInt according to processes and rock types (left) and combining experimental as well as model- ling works in a synergistic workflow (upper right). Work packages are interlinked at methodological levels (lower left) and GeomInt is using Virtual Reality (VR) methods for data and model integration (lower right) Environmental Earth Sciences (2021) 80:509 1 3 Page 5 of 20 509 Experimental platform For the coordinated experimental analysis of the THM prop- erties and behaviour of clay, salt and crystalline rocks an experimental platform has been established within GeomInt, which comprises both laboratory (LAB) and in-situ test in Underground Research Laboratories (URL) (Fig. 2). After a brief introduction of the experimental platform of GeomInt (section “Laboratory experiments (LAB)”) one laboratory research result is highlighted (section “LAB highlight: Ani- sotropy of clay rocks”). Laboratory experiments (LAB) The laboratory facilities of the project partners were used for the characterisation of rock samples from clay, salt, and crystalline for the description of discontinuities under vari- ous loads (THM conditions) and at various spatial scales. Specimens have been gathered from various URLs such as Mt. Terri, Springen and from the Ore Mountains (see Fig. 3) to cover the required barrier rock types. Newly accessible samples have been tested, e.g. from core drilling campaigns in the sandy facies of the Opalinus Clay in the Mont Terri URL. The experimental work on the various rock types by the project partners has been distributed in principle as follows: – CAU 1: investigated Opalinus Clay from the URL Mont Terri (sandy and shaley facies). The rock specimens were sampled and provided by BGR2 – IfG3: contributed their long-term experience salt mechan- ics and investigated salt rock from the Springen salt mine. – TUBAF4: studied fracture behaviour of crystalline rock samples from the Kirchberg mine. – UoS5: used their µXRCT facility for the structural char- acterisation of crystalline rocks. The following types of laboratory experiments have been conducted with the aim of obtaining basic material param- eters on the one hand and of increasing the understanding of more complex multiphysical phenomena on the other: – three-point bending tests – Brazilian tests – conventional triaxial deformation tests – shrinkage and swelling behaviour of Opalinus Clay under free and in-situ stress conditions – fluid-driven percolation in clay and salt using either gas or brine – cyclic (mechanical, hydraulic) percolation tests to study healing and sealing effects – direct shear tests, mechanical and hydraulic fracture properties under normal and shear load conditions – cyclic tests to study the evolution of fracture properties Figure 4 illustrates an overview of the experimental facilities of the LAB platform. The results of the experimental investi- gations are published in a series of papers Rizvi et al. (2019); Vowinckel et al. (2020); Frühwirt et al. (2021); Sattari et al. (2020) and a book publication Kolditz et al. (2021). Here, a general selection is presented: – characterisation of rock properties, e.g. fracture tough- ness, splitting strength, mechanical and thermal proper- ties under various loading conditions (THM), – influence of embedded layering on material strength, deformation behaviour and fracture propagation – influence of anisotropy on shrinkage and swelling effects, – fluid-driven (gas and brine) percolation in clay and salt and its impact on permeability changes, Fig. 3 Geographical workflow including interlinked experimental and modelling works in GeomInt 1 Christian-Albrechts University of Kiel. 2 Federal Institute for Geosciences and Natural Resources. 3 Institut für Gebirgsmechanik GmbH. 4 Technische Universität Bergakademie Freiberg. 5 University of Stuttgart. Environmental Earth Sciences (2021) 80:509 1 3 509 Page 6 of 20 Fig. 4 Experimental facilities of the LAB platform for the mechanical characterisation barrier rocks (clay, crystalline, and salt rocks): (sec- ond row) THM Laboratory Tests - e.g. mechanical properties of all rock types; (third row) Shrinkage and Swelling Tests for clay rocks; Percolation Tests for salt rocks; (lower row) Stress Redistribution Tests: Dynamic big shear box device at TU Freiberg (left: Overview of apparatus; mid: Preparation of sample for shear test, right: Shear box with displacement transducers ready for testing) Environmental Earth Sciences (2021) 80:509 1 3 Page 7 of 20 509 – dependence of crack paths on principal stress directions in anisotropic media. LAB highlight: Anisotropy of clay rocks The published data for Opalinus Clay rock indicate a strong dependence of its thermo-hydro-mechanical properties on sample heterogeneity and inherent anisotropy Sarout et al. (2014); Favero et al. (2018); Minardi et al. (2016). The embedded layering orientation in these sedimentary rocks results in a high anisotropy factor along the material’s principal axes. Therefore, based on fluid percolation tests conducted on cubic clay samples, the dependency of frac- ture openings on the imposed stress state as well as sample orientation were investigated (Fig.5a). A similar approach has been followed for salt rock samples by Kamlot (2009). A true triaxial device is used to apply and control the principal stresses along the axes. Pressurised oil using a syringe pump is used to perform the hydraulic fracturing. The experiments are conducted under two main boundary conditions: the direction of the applied borehole pressure perpendicular or parallel to the embedded layering orientations. In both cases the layering orientation controls the fracking paths and flow surfaces. The inherent anisotropy in Opalinus Clay rocks also affects the material response during the wetting and drying processes Minardi et al. (2016). Thin-section samples were prepared and placed inside a desiccator under different total suction values (Fig. 5b). The suction applied by means of saturated salt solutions ranged from 3.2 up to 367.5 MPa, which results in shrinkage and swelling of the samples. One sample was used to determine the surface linear strains along two different orientations, one parallel to the embedded lay- ering orientation and the other one perpendicular to it. A second sample was used to determine the weight change during the entire process. The room temperature was kept constant around 20 ◦ C during the test period. Equilibrium inside the desiccator is defined once the weight change of the sample during three consecutive readings equalled zero, where tolerance is set to be 1%. Along a drying path, the results indicate higher strains perpendicular to the embedded layers. When the suction exceeded 150 MPa, the strains in the perpendicular direction were almost 4.5 times larger than those parallel to the bedding. Similar behaviour was also observed by Minardi et al. (2016), however, the conducted experiment here was extended to cover a wider range of suc- tion values. Previous results from Minardi et al. (2016) were limited to a suction range of 9.8 up to 139  MPa and did not provided results for air-entry pressure values. Interestingly, along the wetting path, the differences between the strain gauges were much less pronounced. According to the water content data, the air-entry pressure for a sandy Opalinus Clay was found to be around 25 MPa. In‑situ experiments (URL) The in-situ part of the experimental work is formed by the collaborating underground research labs (URLs) from which the sample material was obtained and distributed accord- ing to the defined workflow (Fig. 3) for use in the labora- tory experiments described above. As an example, Opalinus Clay samples have been provided from the AD experiment in Mt. Terri for fluid percolation experiments and simulations Bossart et al. (2017). In addition, several in-situ experiments have been analysed within the GeomInt project such as the CD experiment in Mont Terri Ziefle et al. (2019) or the large wellbore test in Springen Nest and Naumann (2020). These analyses were possible due to links to related pro- jects in which the partners of the present consortium were involved, such as the Mt. Terri and DECOVALEX projects as well as the Geo:N project StimTec Schmidt et al. (2020). URL highlight: In‑situ experiments in Mont Terri The Mont Terri rock laboratory provides unique direct access to the Middle Jurassic Opalinus Clay formation in the Swiss Jura Mountain Belt. It branches off from the security gallery of the motorway tunnel near the town of St. Ursanne (NW Switzerland). Depending on the topography, it lies at 230 m to 320 m depth below ground Heitzmann and Bossart (2001). The Opalinus Clay in the rock laboratory offers a thickness of about 130 m Hostettler et al. (2018), the sedi- mentary bedding plane is dipping with about 40◦ towards SE. Within the Mont Terri project, 21 institutions from nine Fig. 5 Influence of embedded layers on a percolation orientation as well as the flow surfaces, and b suction-induced shrinkage strains Environmental Earth Sciences (2021) 80:509 1 3 509 Page 8 of 20 countries are currently investigating the Opalinus Clay as a potential host rock for the storage of nuclear waste. Since 1996, a total of 1400 m of galleries and niches have been excavated in the Mont Terri rock laboratory (Fig. 17a) and more than 140 experiments have been carried out. Almost 50 of them are still running Bossart and Milnes (2018). Some of them are noteworthy within the context of the GeomInt-Project. Serving the need for core samples for the GeomInt labo- ratory investigation program, the AD-experiment has been initiated. Therein, two fully cored boreholes with a diameter of 131 mm were drilled parallel and perpendicular to the sedimentary bedding, respectively. The 15.35 m long hori- zontal borehole is located entirely in the lower sandy facies. The second borehole is oriented perpendicular to the bed- ding (with a dip of 43◦ ), thus crossing the upper shaly, lower sandy, carbonate-rich sandy-facies and lower shaly facies of the Opalinus Clay. The cores have been geologically char- acterised Galletti and Jaeggi (2019) and analysed by thin section-microscopy and mineralogical-geochemical tech- niques Kneuker and Furche (2021). Furthermore, geophysi- cal in-situ-borehole-measurements have been conducted Furche (2020). These analysis bring out that the Opalinus Clay is characterised by a high intrafacies variability perpen- dicular to bedding, especially for the sandy and carbonate- rich parts, but a low variability parallel to bedding, except for the carbonate-rich sandy facies Kneuker et al. (2020); Kneuker and Furche (2021). The Model-Experiment-Exercise MEX 1-4 is based on the CD/LP experiment, which investigates the long-term cyclic deformation (CD) in a niche of the rock laboratory, as depicted in Fig. 6. The related long-term coupled hydro- mechanically effects are characterised by the interaction of various effects like the redistribution of stresses due to the excavation, swelling and shrinkage of the host rock, related changes of material parameters, as well as plastic and vis- cous effects. Within the GeomInt project these effects have been modelled with a fully coupled hydro-mechanically approach showing consistent results. Nevertheless, the complex interaction of these effects is not yet fully under- stood and is subject of the follow-up experiment CD-A and the GeomInt II project. The Fault Slip (FS) experiment addresses the fault reactivation due to pressure-induced percolation in a low-permeability, large-scale discontinu- ity in the Mont Terri rock laboratory. Within the GeomInt project it was shown that the implemented hydro-mechan- ically coupled model approach, using the Lower-Interface- Element (LIE)-method, is principally capable to reproduce the measured evolution of the deformation at the fault and the injected water volume. URL highlight: In‑situ experiments in the salt mine Kalibergwerk Springen/Merkers The large-scale test site Merkers benefits from the unique mining situation in the bedded salt mass of the Werra salt formation (z1, Zechstein sequence) where two potash seams were mined in a room-and-pillar system at 275 m (1st floor, potash seam “Hessen”, z1KH) and 360 m (2nd floor, potash seam “Thüringen”, z1KTh) depth, respectively. The evapo- rite rocks of the Zechstein formation were laid down dur- ing the Permian period around 250 million years ago. The intact mineral deposit was locally disturbed between 14 and 25 million years ago by tertiary volcanism, leading to the mutation of some potash salts to sylvinite, and the creation of pockets of CO2 under high pressure. The in-situ tests for pressure-driven percolation are con- ducted between the two potash seams in the very homoge- neous Middle Werra rocksalt (z1Na). It consists mostly of very pure halite layers intersected by thin anhydrite lines or bands of rock salt with finely distributed anhydrite acces- sories indicating the sedimentary bedding. Because the test results depend mainly on the acting stress field, i.e. the minimum principal stress distribution in the rock mass around the test area, it has been measured and characterised by hydro-frac measurements, and is thus well-known. The minimal stresses in the contour increase with progressive distance from the underground openings until reaching a constant value at a depth of around 15 m. The measured value of an undisturbed stress state of around 8  MPa corresponds fairly well to calculated lithostatic stresses of 7.8 to 8.8 MPa. The main facility is a large borehole, see Fig. 7, of nearly vertical 60 m height and 1.3 m diameter. It was drilled upwards from the second floor, ending about 20 m beneath the first floor. For access to the later sealed volume an 85 mm pilot hole has been drilled from the upper floor. The borehole was closed by a 20 m thick MgO-concrete plug. Fig. 6 EZ-B Niche in the Mont Terri Rock Laboratory with 3D-Joint- meter Girardin and Nussbaum (2006) Environmental Earth Sciences (2021) 80:509 1 3 Page 9 of 20 509 For monitoring of micro-seismic events, e.g. due to crea- tion of an excavated damage zone around underground open- ings or fluid flow-driven damage, a highly sensitive acoustic emission (AE) network was installed in observation bore- holes, which were drilled parallel to the main borehole. This network has constantly been updated and extended over the past years. Signals of magnitude M < −4 can be detected in a frequency range from 1 kHz to 200 kHz corresponding to intergranular microcracking along grain boundaries on a millimeter- to centimeter-scale. As an example, we present an experiment performed in 2018, where the borehole was loaded with brine. Figure 8 shows the last three pressure spikes over a span of 50 days. They were created by pumping brine with rates of 2 l/h, 3.3 l/h and 5.5 l/h, respectively. Each time, a rapid drop in pressure was observed consistent with the high bulk mod- ulus of brine. This behaviour is very different from gases (think natural gas/nitrogen blanket), where the pressure decreases very slowly. The central part of the figure shows, that the decay of the pressure still continues 6 weeks after the pressure spike, showing the slow expansion of the brine into the surrounding area. The final spike was an attempt to determine the highest pressure obtainable with available equipment. Numerical platform In order to render the gathered experimental data useful for the analysis of current and future geosystems and geotech- nologies, it needs to be integrated into software used for the design and assessment of such applications. An essential scientific goal of the GeomInt project was, therefore, the analysis of potentials and limitations of different numeri- cal approaches for the modelling of discontinuities in the rocks under consideration of multi-physical processes. An important prerequisite is the improved understanding of the diverse set of methods available to the analyst and their syn- ergies with regard to theoretical and numerical fundamen- tals. The following numerical methods were included in the comparison and development efforts: – discontinuous “Lattice Element Method” (LEM), – discontinuous “Discrete Element Method” (DEM), – Lagrangian “Smoothed Particle Hydrodynamics” (SPH), – imaging-based “Forces on Fracture Surfaces” (FFS), – continuum-based “Phase-Field Method” (PFM), – continuum-based “Non-Local Damage” (NLD) formula- tion, – co-dimensional continua-based method using “Lower- Interface-Elements” (LIE) and the “Hybrid-Dimensional Finite-Element-Method” (HDF). The methods were systematically investigated and appro- priately extended based on experimental results. The spe- cific arrangement of the numerical methods along scales represents typical applications within the GeomInt project. In general, numerical methods can be applied at different scales depending on the specific purpose of investigation. For example, smeared (phase field, non-local dam- age) and discrete (co-dimensional LIE) representations of fractures were compared for the case of hydraulic fracture propagation in Yoshioka et al. (2019). While all methods were able to simulate fracture propagation in a toughness- dominated regime, they strongly differ in several aspects: – The type and number of material parameters used: The PFM model for brittle fracture requires a discretisation- related length-scale parameter and otherwise well defined material parameters such as Young’s modulus, Poisson’s ratio and fracture toughness. Other choices are implic- Fig. 7 Schematic overview of the large borehole underground labora- tory Fig. 8 Sequence of three pressure applications over a span of 50 days Environmental Earth Sciences (2021) 80:509 1 3 509 Page 10 of 20 itly made, e.g., by the choice of degradation function and phase-field driving forces. The LIE approach with a cohesive-zone model requires stiffness and damage parameters of the lower-dimensional fracture zone, which in turn can be linked to quantities such as fracture tough- ness (critical energy release rate). – The impact of discretisation and how it can be com- pensated for: For example, an effective critical energy release rate can be determined for the PFM approach that maintains energetic equivalence to the sharp interface approach for a given choice of regularisation length and discretisation. – The accessibility of important physical measures of the fractured system: While methods such as LIE or HDF readily allow a quantification of the fracture aperture, smeared approaches such as PFM and NLD have to determine these values based on additional criteria and post-processing effort. This becomes particularly relevant for coupled processes, where the fracture aperture deter- mines flow and transport through the fractured porous medium. NUM highlight: Phase‑field method for fracturing processes In order to facilitate comparisons such as this one and to advance the development of the individual methods in light of the experimental work packages, so-called model- experiment exercises were introduced, cf. Section “Model- Experiment-Exercises (MEX)”. One of the exercises is con- cerned with pressure-driven fluid percolation and hydraulic fracturing in salt rock samples under anisotropic confining stresses. Cubic samples with a side length of 100 mm were prepared with a borehole drilled in the centre to a depth of 50 mm. The wellbore is cased with a steel tubing with a diameter of 4 mm, leaving the bottom 20 mm open for fluid entrance (Fig. 9). Pressurised oil is injected through the wellbore in the sample until eventual failure. Two different configurations of the confining stresses were considered as shown in Fig. 10. As the vertical stress applied is the minimum principal stress (compression positive) in the first configuration, a crack on the horizontal plane was induced (Fig. 10a). Similarly for the second configuration, a vertical crack is developed according to the confining stress configuration (Fig. 10b). The variational phase-field model bases its origin on the generalised Griffith’s theory Francfort and Marigo (1998) and one of its main advantages is to account for arbitrary numbers of cracks without imposing geometrical restrictions on their topology. It represents a crack set Γ by a scalar phase-field variable, v ∶ � ↦ [0, 1] which transitions from an intact ( v = 1.0 ) to a fully broken state ( v = 0.0 ) continu- ously. The transition length is controlled by a regularisation parameter with the dimension of a length, �s > 0 Bourdin et al. (2000). In the computation, we seek a pair of displace- ment � and phase-field v that minimise the following energy functional given loading in a quasi-static manner Bourdin et al. (2012); Yoshioka et al. (2019) where W is the strain energy density, Gc is the fracture toughness, and cn is a normalisation parameter defined as cn ∶= ∫ 1 0 (1 − �)n∕2d�. The two confining pressure cases above were simulated by the variational phase-field model (Fig. 10). Though the crack propagation direction is not prescribed, the stress con- centration is correctly identified through the energy mini- misation process and the cracks propagated in the expected direction with the expected orientation. NUM highlight: Hybrid‑dimensional finite‑element‑method (HDF) for high resolution investigations on the influence of transient fracture volume changes Another modelling exercise is concerned with the identifica- tion of hydro-mechanical effects throughout characterisation of a single, fluid-filled fracture embedded in crystalline rock. Transient perturbations of the fracture’s equilibrium state by (1) F 𝓁 = ∫ � v2WdV − ∫ ��N � ⋅ � dS + ∫ � p� ⋅ ∇v dV + 1 4cn ∫� Gc ( (1 − v)n 𝓁 + 𝓁|∇v|2 ) dV , Fig. 9 Schematic of the pressure-driven fluid percolation experiment in rock salt Environmental Earth Sciences (2021) 80:509 1 3 Page 11 of 20 509 Fig. 10 Two different confining stress states considered in the fluid percolation experiments in rock salt Kamlot (2009) and simulation results (simulated stresses are displayed in MPa) Environmental Earth Sciences (2021) 80:509 1 3 509 Page 12 of 20 systematic alterations of the fluid pressure induce fluid flow and changes of the fracture volume. The strong interaction of flow, transient geometrical properties such as the fracture’s permeability and changes in fracture volume are studied on a cylindrical sample with a radius of r = 15 mm and height h = 75mm, separated by a single fracture (Fig. 11). The fluid pressure within the fracture is varied by pressure controlled experiments in a triaxial cell through a borehole located at the samples centre. The water-filled fracture is stimulated by harmonic oscillations of the fluid pressure to study hydro-mechanical effects on the pressure p - flow q transients. Variation of the pressure amplitude lead to varying mechanically induced volumetric changes and flow characteristics of the fracture resulting in a non-linear p − q relationship. These effects are sensitive regarding the fracture’s volume, since small abso- lute variations in the � m range might induce large relative changes and is continuously influencing the flow character- istics within a single stimulation period. These phenomena are governed by a hybrid-dimensional element formulation, which is motivated by pressure driven, Poiseuille-type creeping flow conditions and low Reynolds numbers. The lower dimensional flow domain is not explic- itly accounting for the fracture’s geometry, but considers initial fracture apertures in terms of a material parameter and transient changes by means of the deformation state of the surrounding rock matrix. Consistent evaluation of the balance equations results in where p is the fluid pressure and � the fracture aperture. Eq. 2 consists (from left to right) of a) transient term , b) (2) �p �t − �2 12 ��R ∇p ⋅ ∇p − � 12 ��R �� ∇� ⋅ ∇p − � 12 ��R �� ∇ ⋅ (�2 ∇p) + 1 � �� �� �t = qlk quadratic term, c) convective term, d) diffusive term, e) volumetric coupling term and f) leak-off term. Due to their minor contribution to the overall solution terms b) and c) are mostly neglected in numerical investigations. The proposed discrete fracture formulation possesses a high numerical resolution of fracture volume changes and its influence on characteristic diffusion processes. Results displayed in Fig. 12 indicate the correlation of increasing pressure amplitudes and hydro-mechanically induced non-linearities between pressure and flow transients. The proposed model is capable to consider small volumetric changes of the fracture domain and reproduces their influ- ence on the flow processes in transient settings. Model‑Experiment‑Exercises (MEX) The basic idea of Model-Experiment-Exercises (MEX) is to link modelling and experimental work from the very begin- ning – i.e. in the conceptual phase of a project. Due to the inherent complexity of each part in the systems analysis, this combination could sometimes not be maintained as planned. Moreover, both models and experiments require highly sophisticated tools and equipment as well as highly special- ised professionals. This means that adequate measures and incentives for collaboration across institutions and between experimental as well as modelling personnel are required. GeomInt is introducing the MEX concept exactly for this purpose. Therefore, the following MEX were a dominant part of the GeomInt project work and supply the majority of the publications with research material. In order to illustrate the MEX concept, Fig. 2 (upper right) depicts the dependencies between laboratory (LAB) and field experiments (URL) as well as modelling work (MOD). Lab experiments (LAB) will be analysed by mod- els (MOD) in order to calibrate parameters and validate the Fig. 11 Set-up of the stimulated single fracture embedded in a cylindrical sample and visualisation of the applied harmonic fluid pressure: experimental set-up (left), applied boundary conditions (mid), axi-symmetric model set-up (right) Environmental Earth Sciences (2021) 80:509 1 3 Page 13 of 20 509 models [1]. This step in the analysis loop should demonstrate that the models are able to represent the experiments (vali- dation). Validated models will be used for planning experi- mental work in order to improve the general process under- standing [2] (experimental design). Both lab experiments and models must be scaleable to field experiments [3] where it is demonstrated that they can be applied in an engineering- scale context. Finally, the MEX concept was also used as a framework for comparing the different numerical approaches introduced in section “Numerical platform”. Table 1 lists an overview of the main Model-Experiment- Exercises (MEX) performed in GeomInt combining both theoretical work and experiments at various scales. Model- Experiment-Exercises (MEX) are organised mainly along the projects work packages as well as experiments (EXP: LAB/URL) and various modelling approaches according to the methods of the numerical platform presented in Sec- tion “Numerical platform”. Table 1) summarises a matrix of conducted MEX, where processes (WP1 swelling/shrink- ing, WP2 fluid percolation, WP3 stress redistribution) are arranged in lines and numerical approaches (MOD) in col- umns. Gray shaded cells indicated research within GeomInt. A full description of MEX can be found in Kolditz et al. (2020). Here, we present three selected model-experiment-exer- cises (MEX) which emphasise the following aspects: (i) comparison of various numerical methods against experi- mental data (section “Three point bending tests (MEX 0-1a)”), (ii) an example with a new numerical approach for shear tests (section “Direct shear tests using a basalt sample (MEX 3-2)”), and (iii) the role of in-situ experiments for model testing and experimental design (section “Model- Experiment-Exercises in Underground Research Laborato- ries (URL)”). Three point bending tests (MEX 0‑1a) Three numerical methods for fracturing - Lattice Element Method (LEM) Sattari et  al. (2017), Discrete Element Method (DEM), and Variational Phase-Field model (VPF) discretised by the Finite Element Method Yoshioka et al. (2019) - are compared against three point bending tests per- formed on Rockville Granite samples. Three point bend- ing tests were conducted by applying a point load on the top of the sample whose rate is controlled by keeping the Crack Mouth Opening Displacement (CMOD) rate fixed at 0.05 µm/s (Fig. 13a). The notch width is 1.2 mm and the ratio of notch length to sample height is 0.2 Tarokh et al. (2017). The thickness (w), the span (L), and the height (H) of the sample, and Young’s modulus (E), Poisson’s ratio ( � ), the uniaxial compressive strength (UCS), the uniaxial tensile strength ( �T ), and the fracture toughness ( Gc ) are listed in Table 2. Experimental force-displacement results from two samples are shown together with the simulation results in Fig. 13b. The test results exhibit a predominantly linear elas- tic response until the onset of the failure where they show hardening behaviour followed by softening. All the simula- tion results reproduce the linear elastic responses. For the variational phase-field model, this is straightforward as it takes the elastic properties directly. However, for lattice ele- ment and discrete element methods, these bulk properties need to be derived from the micromechanical properties. The main difference between the lattice/discrete element methods and the variational phase-field model is observed in the post-failure behaviour. While the variational phase-field model used in this study is based on brittle fracture and not equipped with a ductile softening behaviour, lattice/discrete element methods are capable of simulating a gradual soften- ing behaviour of the bulk material. Fig. 12 Pressure and flow transients (left) showing an increasing non-linear p − q relationship obtained by numerical investigations of the pro- posed boundary value problem (right) Environmental Earth Sciences (2021) 80:509 1 3 509 Page 14 of 20 Direct shear tests using a basalt sample (MEX 3‑2) At the rock mechanical laboratory of TU Bergakademie Freiberg, direct shear tests using a constant normal stiffness (CNS) have been conducted on samples with pre-existing discontinuities. Large-size intact drill cores with diameter of approximately 150 mm where split parallel to their axis in preparation for shear testing. After removing strips at the boundaries of the discontinuities where irregular and non- typical fracture patterns were observed due to the proximate splitting impact, net shear planes of approximately 125 x 175 mm could be generated. A scheme of a CNS test and a photograph of the shear testing device can be seen in Fig. 14. A basalt sample was sheared multiple times using increasing stiffness values Kn . The shear forces and the ver- tical movement of the sample were recorded. The simula- tion approach explicitly calculates the forces on the surface geometry. For that purpose, 3D scan data were obtained. Basic rock parameters were investigated in standard tests. In order to simulate shear behaviour of rock joints in crys- talline rock a MATLAB-code using a Forces on Fracture Surfaces (FFS) approach was developed. Its main objec- tive is to simulate shear behaviour by only using physical- founded input parameters and realising minimal calculation times. Starting point of the calculation algorithm is a matrix of the actual surface geometry prior to the direct shear test. A quadratic grid is used which is created from the raw scan data by linear interpolation of measurement data to the neighbouring points. The surface is then duplicated and shearing is done by shifting the matrices against each other at consecutive steps of one pixel. Thus, the grid constant is the increment of the shear dis- placement. At each shear step grid contact forces and grid geometry are updated based on classical approaches like Hooke’s formulation for linear elasticity, Coulomb’s fric- tion law and Mohr’s shear strength model. In Fig. 15, a visual comparison of the damaged areas of the rock surface proves that the code is able to detect areas which are damaged and therefore strongly influence the shear behaviour. The simulation results for a selected stiffness value can be seen in Fig. 16. In this case, the destroyed edge of the sam- ple was removed before starting the simulation. The overall quality of the numerical result is considered fair. It clearly can be seen that the dilatation and shear stress curves are tightly correlated. This implies that the key for adequate results is a good approach to model the vertical movements of the sample. In laboratory tests, a negative normal dis- placement at the first stage of the shear tests is measured which can be explained by not perfectly matching sample halves. The code is so-far not able to simulate this behaviour as it considers a perfect match. Finally it can be said that the simple, physically motivated approach works properly in terms of detecting the zones of importance. A general trend for the shape of the shear curves can be displayed but the difference to the laboratory data can be significant. The mentioned perfectly matching parts and the inability to destroy larger parts of the sample at Table 1 Model-Experiment-Exercises (MEX) matrix. Conducted experiments and modelling are coloured in grey Legend: (EXP) experimental data from literature (LIT), LAB and field experiments (URL). Various modelling approaches (MOD): lattice (LEM), distinct (DEM), finite elements (FEM), hybrid-dimensional (HDF), and forces on fracture surfaces (FFS) are developed and compared against experimental results MEX TOP EXP MOD WP LEM DEM FEM HDF FFS 0-1a Bending fracture test LIT 0-1b Bending fracture test (aniso) LAB 0-2 Humidity controlled bending Concept 1-1a Swelling of clay LAB 1-1b Swelling of clay LAB 1-2 Shrinkage of clay LAB 1-3 Desiccation of clay Concept 1-4 CD/LP experiment URL 2-1a Pressure driven percolation LIT 2-1b Pressure driven percolation LAB 2-2 Healing / closure LAB 2-3 Compressible fluids LIT 2-4 URL Springen URL 3-1 CNL test LAB 3-2 CNS test LAB 3-3 Cyclic loading LAB Environmental Earth Sciences (2021) 80:509 1 3 Page 15 of 20 509 once are the main limitations at the moment. The code and the input data will be published and a further development focusing on this problems is possible. For further details see Fruehwirt et al. (2020). Model‑Experiment‑Exercises in Underground Research Laboratories (URL) The Underground Research Laboratory Mont Terri plays an extraordinary role for model testing and experimental design Bossart and Milnes (2018). Figure 17a illustrates some of the ongoing model-experiment activities which are also integrated into a virtual reality context (VR Task) Rink et al. (2020). As an example, the study concerning the CD/CD-A experiments are briefly presented in the fol- lowing (Fig. 17b). Cyclic‑Deformation experiment (MEX 1‑4) The Cyclic-Deformation (CD/LP) experiment focuses on the deformation of a niche, caused by seasonally influenced desaturation of the Opalinus Clay. This experiment was used twofold: (i) for code verification and (ii) for experimental design of the planned CD-A follow-up experiment (see below). In particular, it has been used for the verification of a newly implemented model within the open-source finite element code OpenGeoSys 6 (OGS-6) for strongly hydro-mechanically coupled problems with consideration of desaturation by comparison against the former weakly- coupled version in OGS-5 Ziefle et al. (2017). Hereby, the new implementation in OGS-6 has proven its practical appli- cability under real conditions and its capacity for high per- formance computing. It was thus used for further validation against observations. The CD-A experiment has been prepared in recent years, to distinguish between deformation processes due to stress redistribution and seasonal variations in air humidity that cause saturation (swelling) and desaturation (shrinkage) of the rock and stress redistribution alone. To this end, two identical niches were excavated, one sealed towards the gal- lery and with a high humidity inside to minimise desatu- ration and one open to the general air circulation of the rock laboratory. The measurement campaign was started in October 2019 Ziefle et al. (2019). Pre-calculations with OGS-5 have been a major factors for the dimensioning of Table 2 Parameters used for the simulation of three point bending tests w L H E � UCS �T Gc 30 mm 127 mm 50.8 mm 27.5 GPa 0.175 106 MPa 8.1 MPa 100 Pa-m Fig. 13 Three-point-bending test: numerical interpretation on experimental data Fig. 14 Schematic (left) and laboratory view (right) of a direct shear test with constant normal stiffness boundary conditions. One part of the sam- ple is fixed and the other part sheared off. (from Nguyen et al. (2014), Konietzky et al. (2012)) Environmental Earth Sciences (2021) 80:509 1 3 509 Page 16 of 20 the experimental design, e.g. the distance between the two niches (Fig. 17b). Currently, a numerical 3D model of the CD-A experiment for OGS-6 has been set up and is being evaluated, cf. Fig. 18. The high-performance computing capacity of OGS-6 will assure a reasonable run-time even for this problem taking into account the complex non-linear deformation behaviour of the hydro-mechanically coupled system. Synthesis and outlook As a result of the GeomInt research project a broad com- bined experimental and numerical platform for the inves- tigation of discontinuities due to swelling and shrinkage processes, pressure-driven percolation and stress redistri- bution for important reservoir and barrier rocks (clay, salt, crystalline) has been developed. Model comparisons for damage and fracture phenomena driven by different physi- cal processes provide information on the optimal areas of application of the numerical methods as well as their limita- tions (see Section “Numerical platform”). A comprehensive validation of the platforms (“Proof-of- Concept”) and in particular their integration across institu- tions and work groups was carried out by “Model-Experi- ment-Exercises” (MEX). Each MEX targeted damage and fracture processes driven by different mechanisms along the lines of each work package (see section “Model-Experiment- Exercises (MEX)”). The MEX concept is the central syn- thesis element of GeomInt as it is directly linking models (MOD) with lab experiments (LAB) and paving the way towards the analysis of in-situ experiments (URL) which has been started in the current project phase and will be continued in future activities (see section “Outlook–Future directions”). The project results published in dedicated articles allow an improved understanding of the processes, the methods used and their application-oriented relevance on appropri- ate time and length scales in order to make the planning and implementation of geotechnical uses of the underground safer, more reliable and efficient. An important part of future work is the demonstration of transferability of the experi- mental-numerical concepts and methods to other geotech- nological applications (e.g. deep geothermal energy, energy storage, repository problems, methods for hydraulic stimula- tion, conventional and unconventional resource extraction or tunnel construction). Ongoing activities also aim at intensi- fying the internationalisation started in the current project in cooperation with complementary research projects (e.g. EURAD, DECOVALEX 2023, Mont Terri project). Synthesis–Numerical methods capabilities Significant emphasis in GeomInt was directed on the development of numerical methods for modelling of dis- continuities in various rock types as well as the analysis of Fig. 15 View of a shear plane in a basalt specimen after the CNS-test – comparison of abrasion after the experiment with model calculation result. Damaged areas (as calculated with the model) are marked with red colour overlay and coincide well with visibly, worn parts of the shear plane Fig. 16 Results of the new FFS model approach for normal stiffness Kn = 2.0 MPa/mm. For comparison the results of the lab experiment (Lab) and the classical Barton-Bandis model (Barton et al. 1985 Barton et al. (1985)) (Bar-Ban) are presented alongside with the new model results Environmental Earth Sciences (2021) 80:509 1 3 Page 17 of 20 509 advantages / disadvantages of particular methods for specific fracture pattern at various scales. The numerical methods are described in more detail in the GeomInt book Kolditz et al. (2020). Here the capabilities of the numerical methods used will be compactly summarised in the following tables. Table 3 provides an overview of the capabilities of the numerical methods for simulation of fracturing processes, i.e. (i) crack representation by elements explicitly/separative or by damage variables, (ii) failure criteria, (iii) used mesh types and (iv) fracture aperture calculation. An evaluation of the capabilities of all used methods concerning the analysis of discontinuities is provided in Table 4. An important part of the MEX concept was the application of multiple methods to these exercises (e.g. for MEX 0-1A three methods, LEM, DEM and FEM-VPF have been applied). For the analysis of lab-scale experiments all present methods (discontinuous and continuous) could be applied successfully. The advantages of the discontinuous methods (LEM, DEM) lies in the simulation of small scale processes, e.g. fracture initiation and propagation, and because they rely on fundamental fracture mechanical phe- nomena. Discontinuous methods exhibit some difficulties when it comes to coupled processes and at larger scales where continuum methods (FEM#) are rather strong. We also added 3D capabilities and HPC implementation status to the table. Benchmarking is a typical method in geosciences for evaluating and assessing the validity of experimental and modelling approaches. Related projects and initiatives e.g. are DECOVALEX Birkholzer et al. (2018, 2019), SeS Bench Steefel et al. (2015), InterComp Maxwell et al. (2014), and recently BenVaSim to mention a few. Those benchmarking initiatives have been focused on providing a best fit between measured data and simulated results. GeomInt follows the proven path but providing a systematic study of fracture mechanics methods at different scales, additionally. Fig. 17 Model-Experiment- Exercises in Underground Research Laboratories (URL) Mont Terri. a The left figures shows a plane view on the URL and experiments for which mod- els have been established (blue spheres). The model type (e.g. THM - thermo-hydro-mechan- ical coupled) is indicated by symbols with a tetrahedron. b A detail of the model for experimental design of the CD-A experiment is shown in the right figure, i.e. the moisture distribution for un-/ventilated niches, respectively Fig. 18 Desaturation around the tunnel contour as evaluated by the design simulations for the cyclic deformation experiment. Simulations were performed as part of the GeomInt and Mont Terri projects: Inflow velocity into the tunnel (left), fluctuations of humidity (saturation) (right) Environmental Earth Sciences (2021) 80:509 1 3 509 Page 18 of 20 Outlook – Future directions The main concept of a unified modelling framework for various rock type will be pursued also in the future. After a strongly methodically oriented first phase, the continuation is intended to demonstrate practical applicability under in-situ conditions. In this context, the focus of experimental investi- gations will shift significantly to the realisation and evalua- tion of in-situ experiments in various underground laborato- ries (URLs). For this purpose, the international cooperation with the Mont Terri project and DECOVALEX-2023 will be further intensified. Focused laboratory experiments will be performed particularly where knowledge gaps limit the ability to study the URL experiments. In addition, aspects of digitalisation will be emphasised in two regards, computation and virtualisation. For the efficient numerical simulation of coupled mechanical, thermal and hydraulic processes in the formation and development of discontinuities on the scale under consideration, the adapta- tion of models and algorithms will be further developed in future using high-performance computing (HPC) capabili- ties. The descriptive presentation of structural, experimental and model results in the real context (e.g. URLs) is to be carried out within the framework of an integrated visual data analysis (virtualisation). Acknowledgements The funding of the GeomInt project by the Federal Ministry for Education and Research (BMBF) under grant 03G0866A is highly acknowledged. We thank the Projektträger Jülich (PTJ) for the continuous support throughout the project. We are grateful to the Mont Terri consortium for supporting this research in particular by provi- sion of new samples, access and insight into the underground research laboratory. We thank Leslie Jakobs for support in project management and maintenance of the project webpage. We appreciate very much the thorough evaluation of the manuscript and the revised version by the reviewers based on which the paper could been significantly improved. Funding Open Access funding enabled and organized by Projekt DEAL. Open Access This article is licensed under a Creative Commons Attri- bution 4.0 International License, which permits use, sharing, adapta- tion, 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 Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted 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 licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. References Armand G, Leveau F, Nussbaum C, de La Vaissiere R, Noiret A, Jae- ggi D, Landrein P, Righini C (2014) Geometry and properties of the excavation-induced fractures at the Meuse/Haute-Marne URL drifts. Rock Mech Rock Eng 47(1):21–41 Armand G, Bumbieler F, Conil N, de la Vaissière R, Bosgiraud J-M, Vu M-N (2017) Main outcomes from in  situ thermo-hydro- mechanical experiments programme to demonstrate feasibility of Table 3 Numerical methods comparison Methods Crack representation Failure criteria Mesh Fracture apert. LEM Element separation Element strength Lattice Direct DEM Element separation Element strength Discrete Direct SPH Particle based NA NA Particle based FEM-LIE Element explicit Cohesive law Conforming Direct FEM-VPF Damage variable Fracture mechanics Non-conf. Indirect FEM-NLD Damage variable Stress based Non-conf. Indirect FEM-HDF Element explicit Cohesive law Non- & conf. Direct Table 4 Numerical methods competencies Legend: PD - pre-defined, UD - undefined, Visc. - viscous, diss. - dissipation, 3D - three-dimensional simulation, HPC - High-Performance- Computing http://creativecommons.org/licenses/by/4.0/ Environmental Earth Sciences (2021) 80:509 1 3 Page 19 of 20 509 radioactive high-level waste disposal in the Callovo-Oxfordian claystone. J Rock Mech Geotech Eng 9(3):415–427 Armero F, Garikipati K (1996) An analysis of strong discontinuities in multiplicative finite strain plasticity and their relation with the numerical simulation of strain localization in solids. Int J Solids Struct 33(20–22):2863–2885 Barton N, Bandis S, Bakhtar K (1985) Strength, deformation and conductivity coupling of rock joints. Int J Rock Mech Min Sci Geomech Abstr 22(3):121–140 Bauer S, Class H, Ebert M, Feeser V, Götze H, Holzheid A, Kolditz O, Rosenbaum S, Rabbel W, Schäfer D, Dahmke A (2012) Modeling, parameterization and evaluation of monitoring methods for co2 storage in deep saline formations: The co2-mopa project. Environ Earth Sci 67(2):351–367 Birkholzer JT, Bond AE, Hudson JA, Jing L, Tsang C-F, Shao H, Kolditz O (2018) Decovalex-2015: an international collabora- tion for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (thmc) processes in geologi- cal systems. Environ Earth Sci 77(14) Birkholzer JT, Tsang C-F, Bond AE, Hudson JA, Jing L, Stephans- son O (2019) 25 years of DECOVALEX - Scientific advances and lessons learned from an international research collabora- tion in coupled subsurface processes. Int J Rock Mech Min Sci 122(2018):103995 Bossart P, Jaeggi D, Nussbaum Ch (2017) Experiments on thermo- hydro-mechanical behaviour of Opalinus Clay at Mont Terri rock laboratory, Switzerland. J Rock Mech Geotech Eng 9(3):502–510 Bossart P, Bernier F, Birkholzer J, Bruggeman C, Connolly P, Dewonck S, Fukaya M, Herfort M, Jensen M, Matray J-M, Mayor JC, Moeri A, Oyama T, Schuster K, Shigeta N, Vietor T, Wieczorek K (2017) Mont terri rock laboratory, 20 years of research: introduc- tion, site characteristics and overview of experiments. Swiss J Geosci 110(1):3–22 cited By 44 Bossart P, Milnes AG (eds) (2018) Mont Terri Rock Laboratory, 20 Years, volume 5 of Swiss Journal of Geosciences Supplement. Birkhäuser Basel, 1 edition Bourdin B, Francfort GA, Marigo J-J (2000) Numerical experiments in revisited brittle fracture. J Mech Phys Solids 48(4):797–826 Bourdin B, Chukwudozie C, Yoshioka K (2012) A variational approach to the numerical simulation of hydraulic fracturing. Proc SPE Annu Tech Conf Exhib 2:1442–1452 Buscarnera G (2012) A conceptual model for the chemo-mechanical degradation of granular geomaterials. Géotech Lett 2(3):149–154 Chang KW, Segall P (2016) Injection-induced seismicity on basement faults including poroelastic stressing. J Geophys Res Solid Earth 121(4):2708–2726 Döhmann MJEA, Brune S, Nardini L, Rybacki E, Dresen G (2019) Strain localization and weakening processes in viscously deform- ing rocks: numerical modeling based on laboratory torsion experi- ments. J Geophys Res Solid Earth 124(1):1120–1137 Duretz T, Souche A, Borst R, Le Pourhiet L (2018) The Benefits of Using a Consistent Tangent Operator for Viscoelastoplastic Computations in Geodynamics. Geochem Geophys Geosyst 19(12):2018GC007877 Favero V, Ferrari A, Laloui L (2018) Anisotropic behaviour of opalinus clay through consolidated and drained triaxial testing in saturated conditions. Rock Mech Rock Eng 51(5):1305–1319 Field B, Barton B, Funnell R, Higgs K, Nicol A, Seebeck H (2018) Managing potential interactions of subsurface resources. Proc Inst Mech Eng Part A J Power Energy 232(1):6–11 Francfort GA, Marigo J-J (1998) Revisiting brittle fracture as an energy minimization problem. J Mech Phys Solids 46(8):1319–1342 Frühwirt T, Pötschke D, Konietzky H (2021) Simulation of direct shear tests using a forces on fracture surfaces (FFS) approach. Environ Earth Sci 80:312. https:// doi. org/ 10. 1007/ s12665- 021- 09606-6 Furche M (2020) AD experiment: High resolution ERT measurements in borehole BAD-2. Technical Report TN 2020-64, Federal Office of Topography (swisstopo), Wabern, Switzerland, www.mont- terri.ch Galletti M, Jaeggi D (2019) AD experiment: Drillcore documentation of boreholes BAD-1 and BAD-2. Technical Report TN 2019-06, Federal Office of Topography (swisstopo), Wabern, Switzerland, www.mont-terri.ch Gens A (2010) Soil-environment interactions in geotechnical engineer- ing. Géotechnique 60(1):3–74 Gens A, Guimaräes LN, Olivella S, Sánchez M (2010) Modelling thermo-hydro-mechano-chemical interactions for nuclear waste disposal. J Rock Mech Geotech Eng 2(2):97–102 Ghanbarzadeh S, Hesse MA, Prodanovi M, Gardner JE (2015) Deformation-assisted fluid percolation in rock salt. Science 350(6264):1069–1072 Girardin C, Nussbaum C (2006) CD experiment (cyclic deformations): Installation and procedures of measurement. Technical Report TN 2006-64, Geotechnical Institute Ltd., St-Ursanne, Switzerland Guglielmi Y, Cappa F, Avouac J-P, Henry P, Elsworth D (2015) Seis- micity triggered by fluid injection-induced aseismic slip. Science 348(6240):1224–1226 Hansbo A, Hansbo P (2004) A finite element method for the simulation of strong and weak discontinuities in solid mechanics. Comput Methods Appl Mech Eng 193(33–35):3523–3540 Harrington JF, Cuss RJ, Talandier J (2017) Gas transport properties through intact and fractured Callovo-Oxfordian mudstones. Geol Soc Lond Spec Publ 454(1):131–154 Heitzmann P, Bossart P (2001) Das mont-terri-projekt untersuchungen über den opalinuston im internationalen felslabor. Bulletin für Angewandte Geologie 6:183–197 Hobbs BE, Mühlhaus H-B, Ord A (1990) Instability, softening and localization of deformation. Geol Soc Lond Spec Publ 54(1):143–165 Hostettler B, Reisdorf AG, Jaeggi D, Deplazes G, Bläsi H, Morard A, Feist-Burkhardt S, Waltschew A, Dietze V, Menkveld-Gfeller U (2018) Litho-and biostratigraphy of the opalinus clay and bound- ing formations in the mont terri rock laboratory (switzerland). In Mont Terri Rock Laboratory, 20 Years, pages 23–39. Springer Hudson JA, Stephansson O, Andersson J, Tsang C-F, Jing L (2001) Coupled T-H-M issues relating to radioactive waste repository design and performance. Int J Rock Mech Min Sci 38(1):143–161 Jirásek M (2007) Mathematical analysis of strain localization. Revue Européenne de Génie Civil 11(7–8):977–991 Kabuth A, Dahmke A, Beyer Ch, Bilke L, Dethlefsen F, Dietrich P, Duttmann R, Ebert M, Feeser V, Görke U-J, Köber R, Rabbel W, Schanz T, Schäfer D, Würdemann H, Bauer S (2017) Energy storage in the geological subsurface: dimensioning, risk analysis and spatial planning: the ANGUS+ project. Environ Earth Sci 76(1):23 Kamlot WP (2009) Habilitationsschrift: Gebirgsmechanische Bew- ertung der geologischen Barrierefunktion des Hauptanhydrits in einem Salzbergwerk. des Instituts für Geotechnik der Technischen Universität Bergakademie Freiberg Kneuker T, Furche M (2021) Capturing the structural and compo- sitional variability of Opalinus Clay: constraints from mul- tidisciplinary investigations of Mont Terri drill cores (Swit- zerland). Environ Earth Sci 80:421. https:// doi. org/ 10. 1007/ s12665- 021- 09708-1 Kneuker T, Vowinckel B, Ziefle G, Maßmann J (2020) Opalinus clay- stone from mont-terri, switzerland. In GeomInt - Geomechanical integrity of host and barrier rocks - experiment, modeling and analysis of discontinuities, Terrestrial Environmental Sciences. Springer, Heidelberg Kolditz O, Görke UJ, Konietzky H, Maßmann J, Nest M, Steeb H, Wuttke F, Nagel Th, (eds) (2021) GeomInt - Geomechanical https://doi.org/10.1007/s12665-021-09606-6 https://doi.org/10.1007/s12665-021-09708-1 https://doi.org/10.1007/s12665-021-09708-1 Environmental Earth Sciences (2021) 80:509 1 3 509 Page 20 of 20 integrity of host and barrier rocks - experiment, modeling and analysis of discontinuities. In: Terrestrial Environmental Sciences, page ca. 300. Springer, Heidelberg Konietzky H, Frühwirt T, Luge H (2012) A new large dynamic rockmechanical direct shear box device. Rock Mech Rock Eng 45(3):427–432 Lei Q, Latham JP, Tsang C-F (2017) The use of discrete fracture net- works for modelling coupled geomechanical and hydrological behaviour of fractured rocks. Comput Geotech 85:151–176 Liu P, Yan X, Yao J, Sun S (2019) Modeling and analysis of the acidiz- ing process in carbonate rocks using a two-phase thermal-hydro- logic-chemical coupled model. Chem Eng Sci 207(June):215–234 Loret B, Hueckel T, Gajo A (2002) Chemo-mechanical coupling in saturated porous media: elastic-plastic behaviour of homoionic expansive clays. Int J Solids Struct 39(10):2773–2806 Maxwell RM, Putti M, Meyerhoff S, Delfs J-O, Ferguson IM, Ivanov V, Kim J, Kolditz O, Kollet SJ, Kumar M, Lopez S, Niu J, Pani- coni C, Park Y-J, Phanikumar MS, Shen C, Sudicky EA, Sulis M (2014) Surface-subsurface model intercomparison: a first set of benchmark results to diagnose integrated hydrology and feed- backs. Water Resour Res 50(2):1531–1549 McCartney JS, Sánchez M, Tomac I (2016) Energy geotechnics: advances in subsurface energy recovery, storage, exchange, and waste management. Comput Geotech 75:244–256 McDermott CI, Fraser-Harris A, Sauter M, Couples GD, Edlmann K, Kolditz O, Lightbody A, Somerville J, Wang W (2018) New experimental equipment recreating geo-reservoir conditions in large, fractured, porous samples to investigate coupled thermal, hydraulic and polyaxial stress processes. Sci Rep 8(1):14549 Minardi A, Crisci E, Ferrari A, Laloui L (2016) Anisotropic volumetric behaviour of opalinus clay shale upon suction variation. Géotech Lett 6:1–5 Minkley W, Knauth M, Fabig T, Farag N (2015) Stability and integrity of salt caverns under consideration of hydro-mechanical loading. In Mechanical Behaviour of Salt VIII, pages 217–227. CRC Press Nest M, Naumann D (2020) Rock salt samples. In GeomInt - Geome- chanical integrity of host and barrier rocks - experiment, mod- eling and analysis of discontinuities, Terrestrial Environmental Sciences, pages 22–25. Springer, Heidelberg. in Vorbereitung Nguyen V-M, Konietzky H, Frühwirt T (2014) New methodology to characterize shear behavior of joints by combination of direct shear box testing and numerical simulations. Geotech Geol Eng 32(4):829–846 Parisio F, Tarokh A, Makhnenko R, Naumov D, Miao X-Y, Kolditz O, Nagel T (2019) Experimental characterization and numeri- cal modelling of fracture processes in granite. Int J Solids Struct 163:102–116 Parisio F, Vilarrasa V, Wang WQ, Kolditz O, Nagel T (2019) The risks of long-term re-injection in supercritical geothermal systems. Nat Commun 10(1):4391 Quintal B, Jänicke R, Rubino JG, Steeb H, Holliger K (2014) Sensitiv- ity of S-wave attenuation to the connectivity of fractures in fluid- saturated rocks. Geophysics 79(5):WB15–WB24 Regenauer-Lieb K, Veveakis M, Poulet T, Wellmann F, Karrech A, Liu J, Hauser J, Schrank Ch, Gaede O, Fusseis F, Trefry M (2013) Multiscale coupling and multiphysics approaches in earth sci- ences: applications. J Coupled Syst Multiscale Dyn 1(3):281–323 Ringrose PS, Meckel TA (2019) Maturing global CO2 storage resources on offshore continental margins to achieve 2DS emis- sions reductions. Sci Rep 9(1):17944 Rink K, Kalbacher T, Kolditz O (2012) Visual data exploration for hydrological analysis. Environ Earth Sci 65(5):1395–1403 Rink K, Bilke L, Raith F, Zänker S, Jaeggi D, Kunz H, Maßmann J, Naumov D, Schefer S, Scheuermann G, Shao H, Wang WQ, Ziefle G, Bossart P, Nussbaum Ch, Nagel Th, Kolditz (2020) A virtual exploration of the underground rock laboratory mont terri. In EuroGraphics/EuroVis 2020 WISDOME Contest, pages 1–3 Rizvi ZH, Sembdner K, Suman A, Giri Prasad MJ, Wuttke F (2019) Experimental and numerical investigation of thermo-mechanical properties for nano-geocomposite. Int J Thermophys 40(5): Sarout J, Esteban L, Delle Piane C, Maney B, Dewhurst DN (2014) Elastic anisotropy of opalinus clay under variable saturation and triaxial stress. Geophys J Int 198(3):1662–1682 Sattari AS, Rizvi ZH, Motra HB, Wuttke F (2017) Meso-scale mod- eling of heat transport in a heterogeneous cemented geomaterial by lattice element method. Granular Matter 19(66) Sattari A, et al (2020) Thermo-mechanical fracture behaviour of opa- linus clay - experimental and numercial analysis. preparation for the Geo:N Topical Collection in Environmental Earth Sciences Schmidt P, Dutler N, Renner J, Steeb H (2020) Hydro-mechanical analysis of fracture flow and deformation. J Geophys Res Solid Earth submitted. to be submitted Steeb H, Kurzeja PS, Schmalholz SM (2014) Wave propagation in unsaturated porous media. Acta Mech 225(8):2435–2448 Steefel CI, Yabusaki SB, Mayer KU (2015) Reactive transport bench- marks for subsurface environmental simulation. Comput Geosci 19(3):439–443 Tarokh A, Makhnenko RY, Fakhimi A, Labuz JF (2017) Scaling of the fracture process zone in rock. Int J Fract 204(2):191–204 Ucar E, Berre I, Keilegavlen E (2017) Postinjection normal closure of fractures as a mechanism for induced seismicity. Geophys Res Lett 44(19):9598–9606 Volchko Y, Norrman J, Ericsson LO, Nilsson KL, Markstedt A, Öberg M, Mossmark F, Bobylev N, Tengborg P (2020) Subsurface plan- ning: towards a common understanding of the subsurface as a multifunctional resource. Land Use Policy 90(2019):104316 Vowinckel B, Alfarra A, Maßmann J, Guevara Morel C, Naumov D, Zielfe G, Kolditz O, Nagel T (2020) Computational modeling of hydraulically induced fault slip in argillaceous rock using lower- dimensional interface elements with local enrichment of the dis- placement field. preparation for the Geo:N Topical Collection in Environmental Earth Sciences Yoshioka K, Parisio F, Naumov D, Lu R, Kolditz O, Nagel T (2019) Comparative verification of discrete and smeared numerical approaches for the simulation of hydraulic fracturing. GEM Int J Geomath 10(1):13 Ziefle G, Bossart P, Costabel S, Czaikowski O, Furche M, Graupner B, Hesser J, Jaeggi D, Kolditz O, Königer F, Maßmann J, Rink K, Schuhmann R, Schuster K, Vowinckel B, Wieczorek K (2019) Influence of Humidity on Cyclic and Longterm Deformations: CD-A Experiment, Mont Terri Rock Laboratory. In Berichts- band zur 10. CMM Tagung, “Material – Prozesse – Systeme” in Karlsruhe, Germany Ziefle G, Matray J-M, Maßmann J (2017) Coupled hydraulic-mechan- ical simulation of seasonally induced processes in the mont terri rock laboratory (switzerland). Swiss J Geosci 110(195) Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. GeomInt: geomechanical integrity of host and barrier rocks–experiments, models and analysis of discontinuities Abstract Introduction Methodology Experimental platform Laboratory experiments (LAB) LAB highlight: Anisotropy of clay rocks In-situ experiments (URL) URL highlight: In-situ experiments in Mont Terri URL highlight: In-situ experiments in the salt mine Kalibergwerk SpringenMerkers Numerical platform NUM highlight: Phase-field method for fracturing processes NUM highlight: Hybrid-dimensional finite-element-method (HDF) for high resolution investigations on the influence of transient fracture volume changes Model-Experiment-Exercises (MEX) Three point bending tests (MEX 0-1a) Direct shear tests using a basalt sample (MEX 3-2) Model-Experiment-Exercises in Underground Research Laboratories (URL) Cyclic-Deformation experiment (MEX 1-4) Synthesis and outlook Synthesis–Numerical methods capabilities Outlook – Future directions Acknowledgements References