EMRP1.3 | Multiscale rock damage in geology, geophysics and geo-engineering systems
EDI
Multiscale rock damage in geology, geophysics and geo-engineering systems
Co-organized by GM3/NH3
Convener: Federico Agliardi | Co-conveners: Carolina GiorgettiECSECS, Anne Voigtländer, Christian Zangerl, Patrick Baud, Sergio Vinciguerra
Orals
| Mon, 24 Apr, 14:00–18:00 (CEST)
 
Room -2.21
Posters on site
| Attendance Tue, 25 Apr, 10:45–12:30 (CEST)
 
Hall X3
Orals |
Mon, 14:00
Tue, 10:45
Rock mass deformation and failure at different stress levels (from the brittle regime to the brittle-ductile transition) are controlled by damage processes occurring on different spatial scales, from grain (µm) to geological formation (km) scale. These lead to a progressive increase of micro- and meso-crack intensity in the rock matrix and to the growth of inherited macro-fractures at rock mass scale. Coalescence of these fractures forms large-scale structures such as brittle fault zones, rockslide shear zones, and excavation damage zones (EDZ) in open pit mining and underground construction. Diffuse or localized rock damage have a primary influence on rock properties (strength, elastic moduli, hydraulic and electric properties) and on their evolution across multiple temporal scales spanning from geological time to highly dynamic phenomena as earthquakes, volcanic eruptions, slopes and man-made rock structures. In subcritical stress conditions, damage accumulation results in brittle creep processes key to the long-term evolution of geophysical, geomorphological and geo-engineering systems
Damage and progressive failure processes must be considered to understand the time-dependent hydro-mechanical behaviour of fault damage zones and principal slip zones, and their interplay (e.g. earthquakes vs aseismic creep), volcanic systems and slopes (e.g. slow rock slope deformation vs catastrophic rock slides), as well as the response of rock masses to stress perturbations induced by artificial excavations (tunnels, mines) and loading. At the same time, damage processes control the brittle behaviour of the upper crust and are strongly influenced by intrinsic rock properties (strength, fabric, porosity, anisotropy), geological structures and their inherited damage, as well as by the evolving pressure-temperature with increasing depth and by fluid pressure, transport properties and chemistry.
In this session we will bring together researchers from different communities interested in a better understanding of rock deformation and failure processes and consequence, as well as other related rock mechanics topics. We welcome innovative and novel contributions on experimental studies (both in the laboratory and in situ), continuum / micromechanical analytical and numerical modelling, and applications to fault zones, reservoirs, slope instability and landscape evolution, and engineering applications.

Orals: Mon, 24 Apr | Room -2.21

Chairpersons: Federico Agliardi, Anne Voigtländer, Patrick Baud
14:00–14:05
Experimental studies on rock mechanics
14:05–14:15
|
EGU23-1518
|
EMRP1.3
|
On-site presentation
Ik Woo

The pore size and the distribution of individual or connected pores contribute to the porosity in a rock which is closely related to rock weathering degree and rock strength. The chemical reaction is normally higher for the larger specific surface area which is closely related to the pore size distribution in a rock. The variation of pore size distribution in sedimentary rocks from Gyeongsan basin in Korea was determined by the laboratory artificial acceleration weathering experiment using peristatic pumps. The pore size distribution of rock specimens was measured by the nitrogen gas adsorption method using BELSORP-max II of Microtrac MRB. The pore characteristics were measured on the outer surface and the innermost part of rock samples to determine the variation of pore size distribution since the outer surface was directly affected by weathering processes while the innermost part was not. The high-purity nitrogen gas is used to evaluate the pore size distribution with different methods such as BET, BJH, and HK. The overall pore volume and size have been increased by the weathering experiment for the tested sedimentary rocks-sandstone, conglomerate, and shale. The increase of macropore in sandstone by weathering experiment leads mainly to the increase in pore volume, while the rise of micropore and mesopore in conglomerate drives the increase of pore volume. 

Acknowledgments

This research was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education(NRF-2020R1F1A107576412).

How to cite: Woo, I.: Pore Size Redistribution by Laboratory Weathering Tests on Sedimentary rocks, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-1518, https://doi.org/10.5194/egusphere-egu23-1518, 2023.

14:15–14:25
|
EGU23-8050
|
EMRP1.3
|
ECS
|
On-site presentation
Manuel Asnar, Christoph Sens-Schönfelder, Audrey Bonnelye, and Georg Dresen

In rocks and concrete, dynamic excitation leads to a fast softening of the material, followed by a slower recovery process where the material recovers part of its initial stiffness as a logarithmic function of time. This requires us to exit the convenient framework of time independent elastic properties, linear or not, and investigate non-classical, non-linear elastic behavior.

These phenomena can be observed during seismic events in affected infrastructure as well as in the subsurface. Since the transient material changes are not restricted to elastic parameters but also affect hydraulic and electric parameters as well as material strength, as documented for instance by long lasting changes in landslide rates, it is of major interest to characterize the softening and recovery phases. It may help us gain more insight in hazard prediction from both a geological and engineering perspective.

The underlying physics behind those non-classical, non-linear effects, sometimes referred to as Nonlinear Mesoscopic Elasticity”, are not agreed upon. There is a lack of experiments that would allow us to discriminate between the existing models.: we aim to contribute to filling that knowledge gap.

Our experiments are made on a sample of Bentheim sandstone, initially dry and then fully saturated, in a triaxial cell. We subject the sample to loading and holding cycles in the microstrain range, while also varying confining pressure and pore pressure. Active acoustic measurements during those loading cycles with an array of 14 piezoelectric sensors allow us to monitor relative velocity changes during the experiment by using Coda Wave Interferometry (CWI).

We observe the dynamic softening as well as the recovery processes in the sample during repeated loading phases of different durations. We find that characteristics of the observed velocity changes vary depending on the observed sensor combination, indicating spatial variability of the response, as well as depending on the lapse time and frequency content of the acoustic measurements that we perform the CWI on.

These experiments serve to estimate the exact capabilities of our experimental setup in terms of signal quality, signal stability and lapse time dependent decorrelation of coda waves. We expect our results to inform a future series of similar but more refined experiments addressing the pore pressure dependence of the non-classical response of rocks.

How to cite: Asnar, M., Sens-Schönfelder, C., Bonnelye, A., and Dresen, G.: Non-classical, non-linear elasticity in rocks: experiments in a triaxial cell with pore pressure control, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8050, https://doi.org/10.5194/egusphere-egu23-8050, 2023.

14:25–14:35
|
EGU23-4578
|
EMRP1.3
|
On-site presentation
Marie Violay, Corentin Noel, and Patrick Baud

The presence of water causes a dramatic reduction of the strength of most rocks. Under compressive stress conditions, fracture mechanics models show the strength of a rock sample is in particular controlled by frictional parameters and the fracture toughness of the material. Previous studies suggested that these parameters could change significantly in the presence of water, but there is a paucity of data quantifying this. Here, we report fracture toughness, frictional and uniaxial compression tests performed on five sandstones and five limestones under dry and water-saturated conditions, that provide new insight into the mechanical influence of water on sedimentary rock strength. Our new data showed that on both sandstones and limestones, the presence of water causes a reduction of both the fracture toughness (from 0 to 50%) and the static friction coefficient (from 0 to 40%), suggesting that water weakening in these sedimentary rocks is mostly due to a reduction of these two parameters under the relatively high strain rate conditions investigated here. While for sandstone we found a reduction of the Uniaxial Compressive Stress between 0 to 35%, it was less variable in limestone, in most cases around 40%. The measured fracture toughness and frictional parameters were then introduced into two well-known micro-mechanical models (the pore-emanating cracks model and the wing crack model), which provide simple theoretical expressions for the Uniaxial Compressive Strength. We found that the predicted water-weakening based on our toughness and friction parameter measurements is in overall agreement with our strength measurements on dry and wet samples.

How to cite: Violay, M., Noel, C., and Baud, P.: Effect of water on sandstone and limestone, fracture toughness, frictional parameters and brittle strength., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-4578, https://doi.org/10.5194/egusphere-egu23-4578, 2023.

14:35–14:45
|
EGU23-14779
|
EMRP1.3
|
On-site presentation
Fabio Trippetta, Roberta Ruggieri, Hem B. Motra, and Cristiano Collettini

Pressure, temperature, and infilling fluids influence the petrophysical properties and the associated damaging processes of rocks at all scales. Moreover, each fluid-rock system possesses peculiar mechanical behaviours being these particularly complex in carbonate rocks hosting fluids. In this work, we analyze the laboratory results of deformed clean and hydrocarbon-filled limestones under varying pressure and temperature, providing links between recorded physical properties (seismic velocity), fluid behavior, and damaging. We focus on carbonate-bearing reservoir (Bolognano Formation) rocks, sampled in the Majella massif (Central Italy) that represents a very good analogue for buried carbonate reservoirs. This reservoir is composed by calcarenites with connected porosity of about 20% saturated by hydrocarbon in the solid state at the outcrop conditions. We performed hydrostatic, triaxial and true-triaxial deformation tests up to a temperature of 100º C and a confining pressure up to 100 MPa on both clean and naturally hydrocarbon-filled limestone samples. Results show increasing seismic velocity and Young’s modulus with increasing confining pressures for both clean and saturated samples as expected. However, different results are observed when the temperature is increased. At low temperatures saturated samples show larger seismic velocity and rigidity with respect to clean samples whilst at higher temperatures the opposite occurs. In particular, when temperature is rised up to 100º C the Young’s modulus of the saturated samples dramatically decreases, being this coupled by a clear volume reduction even during hydrostatic tests (no differential stress applied). Accordingly, microstructural observations highlight grain crushing related to a large amount of randomly distributed cracks within saturated samples. On the contrary, clean samples are characterized by few microfractures, pointing out the primary role played by liquid hydrocarbons. These observations are in good agreement with meso and microstructural features observed on outcropping hydrocarbon-filled carbonate-bearing faults. The presence of fluid hydrocarbons (high temperature) severely weakens the rock promoting fracturing whilst at lower temperature the presence of solid hydrocarbons increases the mechanical properties of hydrocarbon-bearing rocks. These observations have a large impact for the petrophysical characterization of reservoirs and for the understanding microscale to mesoscale mechanisms of deformation and fluids movement along deformed rock volumes.

How to cite: Trippetta, F., Ruggieri, R., Motra, H. B., and Collettini, C.: Mechanical behavior of porous carbonates as a function of pressure, temperature, and fluid content from laboratory experiments and correlation with larger scale structures, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14779, https://doi.org/10.5194/egusphere-egu23-14779, 2023.

14:45–14:55
|
EGU23-5667
|
EMRP1.3
|
ECS
|
On-site presentation
Adarsh Tripathi, Mohammad Mohasin Khan, Ashok K. Singh, and Anindya Pain

Rocks are very much susceptible to deformation in tension, especially under elevated temperatures. Therefore, the study of the dynamic tensile behaviour of rock exposed to high temperature is highly significant to understand the tensile deformation behaviour in dynamic loading conditions which will be proved useful in a variety of engineering problems such as quantifying the blast load impact in fire affected underground/opencast coal mine regions; assessment of ground subsidence due to coalmine fire coupled with blast loading etc. The Jharia coalfield region, known as the coal capital of India, is affected by pervasive underground coalmine fire for decades resulting in small to large-scale surface fracturing. So,the present study focuses on the effect of high temperature on dynamic tensile behaviour and its relation with micro-mineralogical properties of subsurface coal-bearing sandstone samples from a fire-affected mine.  To achieve the objective, the prepared samples were kept in the furnace for 24h with a heating rate of 5°C/min and then allowed to cool down naturally within the furnace. Samples were divided into nine groups based on the thermal treatment at 25 °C, 100 °C, 200 °C, 300 °C, 400 °C, 500 °C, 600 °C, 700 °C, and 800 °C. Using the Split-Hopkinson Pressure Bar (SHPB), the indirect dynamic tensile strength was measured for each group. Based on the obtained results, the indirect dynamic tensile strength of heat-treated specimens is characterized into three zones; viz.: Zone 1 (25-400°C), Zone 2 (400-600°C) and Zone 3 (600-800°C). In zone 1, an increase in average indirect dynamic tensile strength is observed with elevated temperature. However, in zone 2, a sharp decreasing trend in indirect dynamic tensile strength was observed with increasing temperature. This zone is characterised by a progressive increase in thermal cracks and porosity which is possibly the prime reason for a sharp transition in thermal properties. An overall reduction in indirect dynamic tensile strength is observed within zone 3, however, the rate of reduction is gentle. The plasticity that occurred due to high temperature was responsible for a slow rate of reduction in indirect dynamic tensile strength.

How to cite: Tripathi, A., Khan, M. M., K. Singh, A., and Pain, A.: Dynamic tensile behaviour of Barakar sandstone under high-temperature conditions, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5667, https://doi.org/10.5194/egusphere-egu23-5667, 2023.

14:55–15:05
|
EGU23-16739
|
EMRP1.3
|
ECS
|
On-site presentation
Filip Adamus, Ashley Stanton-Yonge, Tom Mitchell, Dave Healy, and Phil Meredith

3D printing is a rapidly evolving technology that has proven useful for a wide variety of disciplines and industries. However, knowledge of its applicability to the fields of rock and soil physics remains limited. 3D printing allows the design of samples with any desired microstructural composition, enabling independent control of properties such as pore space fabric, size, and density; a feat impossible to accomplish with naturally occurring geomaterials. The use of 3D printed samples is therefore highly attractive for relating the effective properties of heterogeneous materials to their microstructural arrangement, a key subject in the fields of rock and soil physics.

This study aims to characterize the physical properties of 3D printed materials (i.e., elasticity parameters, porosity, permeability) and evaluate whether they are suitable to be used as proxies for heterogeneous geomaterials. Two distinct 3D printing technologies were employed for this purpose: the Fused Deposition Modelling (FDM), and Stereolithography (SLA) methods. The FDM method constructs 3D objects by superposing layers of polymer-based filament through a heated nozzle, whereas the SLA method, also known as resin 3D printing, uses a laser light source to cure liquid resin into hardened plastic.

Samples with a variety of pore shapes (sphere, needle, penny shaped), sizes, and pore densities were designed and printed as cylindrical samples of 25 mm diameter and 62.5 mm height. Samples were then subjected to uniaxial compression to measure their effective elastic parameters (elasticity modulus and Poisson’s ratio), and these measurements were compared with theoretical predictions. Preliminary results indicate that the FDM printing method is inadequate for representing a heterogeneous solid composed of an isotropic matrix and void space, due to the intrinsically anisotropic fabric resulting from the layer-by-layer printing method. Additionally, samples with a porous microstructure appear to be effectively stiffer than the intact material, which is attributed to enhanced material sintering surrounding the edges of the void spaces. On the other hand, SLA printing appears to hold more promise and be able to represent a composite material composed of an isotropic matrix with a heterogeneous void space. Further measurements need to be made to confirm these preliminary findings, and this work is currently in progress. 

How to cite: Adamus, F., Stanton-Yonge, A., Mitchell, T., Healy, D., and Meredith, P.: Physical properties of 3D printed materials and their applicability as proxies for heterogeneous geomaterials, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16739, https://doi.org/10.5194/egusphere-egu23-16739, 2023.

15:05–15:15
|
EGU23-11447
|
EMRP1.3
|
ECS
|
On-site presentation
Patrick Bianchi, Paul Antony Selvadurai, Antonio Salazar Vásquez, Luca Dal Zilio, Claudio Madonna, Taras Gerya, and Stefan Wiemer

At present, a reliable method for forecasting earthquakes has not been developed yet, as the physical mechanisms that generate them are very complex and still not completely understood. To overcome the difficulties of retrieving direct observations and measurements in the field, here we employ laboratory and numerical models to investigate and better understand strain localization preceding mainshocks.

We perform a failure test on an intact and dry sample of Berea sandstone confined at 20 MPa with a triaxial machine (LabQuake). Employing in-house developed, conical-type and fully calibrated piezo-electric transducers (PZT), we are able to investigate the acoustic emission (AE) clouds by relocating the single events and by computing their focal mechanisms and scalar seismic moments. The PZT sensors are also used actively to allow for the construction of inhomogeneous and anisotropic velocity models. We further employ distributed strain sensing (DSS) with optical fibers to capture the heterogeneous spatial distribution of the surface strain by gluing the fibers on the sample surface. We observe AE clustering in two regions located at the top and bottom of the rock specimen throughout the majority of the experiment. As the test approaches the main failure, AE localize at one side of the sample in the lower half before obliquely propagating upwards by forming a macro-fracture. Surface strain heterogeneities are detected during the experiment, and regions of higher extensional strain correlate in time and space with rock volumes experiencing high AE activity. Numerical simulations, which are conducted using a two-dimensional continuum-based and fully coupled seismo-hydro-mechanical poro-visco-elasto-plastic modelling tool (H-MEC), are validated with both AE and DSS data. The combination of laboratory and numerical investigations allows us to individuate and study physical mechanisms (e.g., visco-plastic compaction of pores and shear banding) that explain the processes responsible for both surface strain concentration and the generated AE clouds. These findings suggest that the deformation in the interior of the sample is mainly occurring inelastically and is localized along an obliquely forming shear band. We estimate the partitioning between seismic and total deformation to be ~0.09 % and this effectively confirms the previous evidences related to the irreversible localization of strain within the rock specimen. 

How to cite: Bianchi, P., Selvadurai, P. A., Salazar Vásquez, A., Dal Zilio, L., Madonna, C., Gerya, T., and Wiemer, S.: Evidence of Strain Localization Preceding Rock Failure: Insights From Laboratory and Physics-Based Poroelastic Models, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11447, https://doi.org/10.5194/egusphere-egu23-11447, 2023.

15:15–15:25
|
EGU23-14869
|
EMRP1.3
|
ECS
|
On-site presentation
Federico Franzosi, Chiara Crippa, Roberto Garzonio, Stefano Casiraghi, Marc-Henry Derron, Michel Jaboyedoff, and Federico Agliardi

The reliable assessment of the fracturing state of rock masses is a fundamental step towards the evaluation of their geomechanical quality and the quantification of their hydraulic and mechanical properties. Traditional field discontinuity mapping techniques remain fundamental to collect statistical populations of discontinuity attributes and characterize rock mass structure and quality. However, point-like field surveys are strongly biased by scale and orientation. The development of 3D surveys allowed to partly overcome this problem by providing high-resolution point clouds. These allow a robust characterization of fracture geometry but require significant mapping efforts. Here we proposed a quantitative contactless approach to rock mass fracturing assessment by the use of Infrared Thermography (IRT).

IRT is increasingly used in rock-mechanics to characterize rock porosity/fracturing and to monitor rock mass stability, by measuring the thermal response of rock materials to heating or cooling. However, existing IRT applications to the geomechanical study of rock masses are mostly qualitative and lacking sound theoretical and experimental foundations. Starting from the laboratory scale, studying the thermal behaviour of rock samples with different fracture degrees, we propose a quantitative approach to quantify rock mass fracturing, that combines IRT rock temperature monitoring during cooling with the quantification of different descriptors of fracturing state suitable for different analysis scales (laboratory vs in situ).

As a field laboratory we used the Mount Gorsa porphyry quarry (Trento, North Italy), characterized by a homogeneous rock type but strongly variable fracturing states related by complex structurally-controlled and progressive slope damage processes.

We performed a field campaign on quarry front making a) Geometrical UAV surveys and b) field Geological Strength Index (GSI) evaluation on typical spots, c) we carried out IRT monitoring during night cooling using FLIRT1020 thermal camera.

During data processing d) thermal data acquired were corrected by environmental effects (blue sky radiation, slope inclination etc.) adopting original and ad hoc calibrated filters to skim the thermal response from geometrical and external biases. Finally we try to find a correlation between the thermal response of rock-mass outcrop to their quality index.

Our results support the possibility to upscale the analysis to field conditions in order to account for the radiative characteristics of natural environments, the limitations of the technique and upscaling issues typical of fractured rock-mass, taking into account that fracturing metrics (used in laboratory phase) at rock-mass scale, should influence block size distributions, which is fundamental in the evaluation of quality indices, e.g the Geological Strength Index (GSI) widely used in engineering applications.

Emphasising all these issues, the goal of our work is to investigate the relationship between the thermal response of rock mass quality index, through an experimental method developed at laboratory scale and upscaled to in situ conditions.

How to cite: Franzosi, F., Crippa, C., Garzonio, R., Casiraghi, S., Derron, M.-H., Jaboyedoff, M., and Agliardi, F.: In situ assessment of rock mass fracturing using infrared thermography, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14869, https://doi.org/10.5194/egusphere-egu23-14869, 2023.

15:25–15:35
|
EGU23-3096
|
EMRP1.3
|
On-site presentation
Céline Mallet, Gautier Laurent, and Mohamed Azaroual

Heterogeneities control rock properties, especially hydraulic and geophysical properties. Complex systems typically include multiple porosities at embedded scales, from the micro/meso cracks and pores to geological macro-fractures and karsts. This complex network play coupled roles and introduces difficulties in the characterization of the whole formation.

In order to constrain these coupled effects, we use seismic to acoustic data to characterize a multi-scale double porosity network and to understand the corresponding flow and mechanical properties of a shallow aquifer reservoir. The study focusses on the platform “Observatoire des transferts dans la Zone Non-Saturée” (O-ZNS, Orléans, France), an artificial excavation in the karstified and fractured limestone formation of Beauce aquifer. It is composed by an exceptional well (20 m-depth, 4 m-diameter) surrounded by 8 cored boreholes.

Two seismic refraction profiles crossing the O-ZNS site were carried out to determine P-wave velocities. The profiles delineated three main geological units: (i) a clayey soil (0-2 m), (ii) a weathered and karstified limestone layer (2-7 m), and (iii) massive limestone down to the underlying Molasse du Gâtinais layer at a depth of 25 m. In consistence with the lithological log, a thin layer of more massive limestone is highlighted around 5 m-depth. In addition, we also observed that the increase in P-wave velocity slows down after 15 m. This effect is consistent with the increasing fracture density and karst development observed on the direct log imagery and on the well 3D scan. In the massive thin limestone layer of 5 m-depth, the interpreted relative crack density is low, around 0.08. However, in the last layer from 15 to 20m-depth, the relative crack density is much more important, even so discrepant, with maximal values around 0.4.

In parallel to large scale field investigation, mechanical tests and elastic wave velocities have been measured on representative core samples. A strong discrepancy is observed, whatever the property. For example, at 16 m-depth, P-wave velocities are distributed from 3,650 to 5,700 m.s-1 and the corresponding mechanical parameter of crack density ranges from 0 to 0.5. In addition, extreme values of crack density, above 1 are observed around 19 m-depth. These large discrepancies and crack density values are consistent with mechanical behavior and microstructure observation made directly on core samples, even though some samples are more porous than cracked and the distinction need to be kept. Samples are then classified through image processing in three categories: the porous ones, the cracked ones, and the mixed ones allowing to discuss and organize the heterogeneity distribution of the O-ZNS.

To complete the study, an intermediate characterization is running on metric blocs sampled at different depth from the O-ZNS well. Their analyses include 3D external scans at high resolution (as for the surface of O-ZNS well) and P-wave velocity measurements at intermediate frequencies. These blocs are then iterativelly, cut into smaller blocs and re-characterized in order to obtain the distribution of heterogeneity size and characteristics with depth targeting the determination of a REV (Relative Elementary Volume) for future modeling developments.

How to cite: Mallet, C., Laurent, G., and Azaroual, M.: Seismic to acoustic characterization of geomechanical and microstructural properties of a vadose zone heterogeneous limestone formation, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3096, https://doi.org/10.5194/egusphere-egu23-3096, 2023.

15:35–15:45
|
EGU23-6698
|
EMRP1.3
|
ECS
|
On-site presentation
Reinhard Gerstner, Erik Kuschel, Christine Fey, Klaus Voit, Gerald Valentin, and Christian Zangerl

How to implement rock bridges and rock bridge failure in slope stability analysis is an ongoing discussion within the rock mechanic and landslide community. Although there has been intensive research over several decades, there is still a lack of knowledge on how to measure intact rock bridges on rock slopes, how to quantify their impact on rock mass strength, and how they affect the initial failure mechanism. Therefore, we present the analysis of a rock fall case study located in the alpine environment of southern Salzburg (Austria), where a rock slope composed of a polymetamorphic rock mass hosted three rock fall events in the year 2019. The primary aim of this study is the reconstruction of the multiphase failure event and the investigation of the influence of the discontinuity network with its intact rock bridges on the initial failure mechanism.

In our study, we performed a detailed reconstruction of the rock fall process by helicopter-borne event documentation. Moreover, we identified the rock fall failure mechanism by analysing a video capturing the first rock fall event.

Furthermore, we developed a high-resolution digital surface model of the complex post-failure topography by unmanned aerial vehicle photogrammetry (UAV-P) with real-time kinematics (RTK). Based on this model, we map the location, orientation and persistence of pre-existing discontinuities and identify failed intact rock bridges on the rupture surface of the unstable rock slope.

Additionally, we conducted point load and direct shear tests in the rock mechanic laboratory. We applied the former on block specimens to derive the uniaxial compressive strength of the intact rock. The latter allowed us to estimate the Mohr-Coulomb shear strength properties of intact rock and of failure planes, which formed sub-parallel to foliation planes in course of the test procedure.

After the third rock fall event of 2019, a ground-based interferometric synthetic aperture radar (GbInSAR) was installed for 166 days to monitor the actual deformation of the rock slope. We analysed the obtained deformation data at mm resolution to detect zones of ongoing slope movements.

Finally, we integrate the topographical and geological model, the structural inventory, and the geomechanical properties into a 2D numerical model based on the distinct element method (UDEC). We use Voronoi tessellation to allow the development of any failure path within intact rock bridges. By varying the persistence of pre-existing discontinuities and the shear-strength properties of rock bridges, we study the impact of rock bridge location, spatial distribution, and strength on the initial failure mechanisms of the rock slope. We validated the distinct element model by comparing its outcome with the essential characteristics of the rock fall observed in the event reconstruction and deformation monitoring.

By this integrated approach of methods applied to a polyphase rock fall process, we show that the initial rock fall failure mechanism is sensitive to the spatial distribution of rock bridges and their assigned shear strength properties.

How to cite: Gerstner, R., Kuschel, E., Fey, C., Voit, K., Valentin, G., and Zangerl, C.: Rock bridge control on the failure mechanism of a rock fall in a metamorphic rock mass, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6698, https://doi.org/10.5194/egusphere-egu23-6698, 2023.

Experimental studies on fault mechanics
Coffee break
Chairpersons: Sergio Vinciguerra, Carolina Giorgetti, Christian Zangerl
16:15–16:25
|
EGU23-8807
|
EMRP1.3
|
On-site presentation
Elena Spagnuolo, Chiara Cornelio, Stefano Aretusini, Giacomo Pozzi, Massimo Cocco, Paul Selvadurai, and Giuseppe Di Stefano

We present a novel apparatus designed to investigate the mechanical and chemical processes active during the nucleation and the subsequent propagation of a seismic rupture. The earthquake is experimentally represented by the sudden frictional sliding of two blocks caused by either: i) the passage of a rupture front from a nearby seismogenic source at prescribed slip velocity, or ii) by the sudden release of strain energy cumulated during the slow (tectonic) loading stage preceding the nucleation of seismic rupture. M.E.E.R.A. (Mechanics of Earthquake and Extended Rupture Apparatus) is a biaxial horizontal machine installed at the laboratories of the Istituto Nazionale di Geofisica e Vulcanologia of Rome (Italy) thanks to a grant funded by the Italian Dipartimento di Protezione Civile. MEERA works on two blocks sized 320x80x50 mm3 put in frictional contact under a normal load up to 30 MPa. Blocks can be either rocks or analogue materials. The normal load and the shear stress are supplied by 6 hydraulic piston cylinders. One piston applies the tangential force up to 150 kN and up to 40 mm/s of slip rate. The other 5 cylinders modulate the normal force on the 320 x 50 mm2 contact surface. The 6 pistons are mounted on a rigid stainless-steel vessel that can be closed by a top built in plexiglass, which enables the environmental chamber for fluid confinement. The plexiglass top resists up to 6 MPa of fluid pressure exerted and controlled by using two ISCO pumps.  MEERA is designed following the outline described in McLaskey and Yamashita (2017) and introduces three novelties: the control in displacement and displacement rate of the tangential piston up to 1kS/s; the environmental chamber; the rigid stainless-steel frame. MEERA is designed to study how the tectonic loading of a frictional interface composed of natural rocks determine the stress state and shear stress evolution governing seismogenic processes. To this end, the simulated fault in MEERA is equipped with acoustic sensors, strain gauges, optical fibers and high velocity cameras to measure and constrain rupture nucleation processes and earthquake source parameters, including directivity and rupture velocity, the dynamics of seismic ruptures and the earthquake energy budget at different scales. We aim at comparing the laboratory observations and the signals collected by MEERA with those collected by the newly developed on-fault observatory of the ERC FEAR project in the Bedretto Underground Laboratory for Geosciences and Geoenergies (BULGG, Swiss Alps) to provide novel insights in earthquake mechanics.

 

How to cite: Spagnuolo, E., Cornelio, C., Aretusini, S., Pozzi, G., Cocco, M., Selvadurai, P., and Di Stefano, G.: A novel apparatus to study the mechano-chemical processes active during the nucleation and propagation of earthquakes (MEERA)., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8807, https://doi.org/10.5194/egusphere-egu23-8807, 2023.

16:25–16:35
|
EGU23-11449
|
EMRP1.3
|
On-site presentation
Yuval Boneh, Markus Ohl, Oliver Plümper, Greg Hirth, and Matěj Peč

Talc is an important product of several hydration and dehydration reactions in deep faults and subduction zones. The unique weakness of talc along its basal planes makes it an essential component in understanding various fault slip behaviors (e.g., episodic vs continuous slip, seismic vs aseismic) or realistic geodynamic models. A recent experimental study by Boneh et al. (2023) on talc mechanical behavior at high P-T conditions highlighted: (i) talc’s low friction coefficient under all conditions (<0.14), with thermal weakening down to µ~0.01 at 700 °C. (ii) Grain-scale microstructures demonstrate a component of fracturing and microcracking under all conditions tested. And (iii) pressure-dependence of talc strength decreases at higher temperatures, where there is also a greater tendency for localization. A vital part of depicting mineral rheology is the understanding of their underlying mechanisms of deformation associated with the observed bulk mechanical and microstructural behavior. To reveal the underlying deformation mechanism/s we analyzed the deformed samples through high-resolution transmission electron microscopy (TEM) at Utrecht university of samples prepared using a focus ion beam (FIB). Five talc samples were examined – an undeformed sample, and samples deformed at 400, 600, and 700°C under 1 GPa, and at 400 °C under 1.5 GPa.

Seven FIB lamellae sampled areas adjacent to the main fracture (if exists) or high damage zones. The starting material shows talc flakes with a thickness of ~100-400 nm without a sample-scale preferred alignment. The sample deformed under 400°C and 1.5 GPa exhibits distributed deformation with opening cracks along talc basal planes and pervasive kinking normal to the basal planes. The sample deformed at 400°C and lower pressure (1.0 GPa) exhibits thin lamination (~50 nm) well oriented with the orientation of the main fracture plane. The sample deformed at 600°C exhibits crystal delamination along the basal cleavage (forming grain fragments <10 nm in width) along the main fracture. The sample deformed at 700°C exhibits more areas of high damage, possibly due to the similar basal-cleavage delamination. A key incentive is to relate the observed nano-scale crystal defects with the bulk mechanical behavior and with processes that might promote the localization of deformation. Pressure-dependent strength can be accounted for by kinking and kinking-induced porosity while thermal weakening can be related to temperature-dependent mobility of crystal defects leading to delamination along the basal cleavage. We will discuss possible physical mechanisms of talc deformation and the prospect of extrapolating the mechanical behavior of talc achieved at the lab to the range of conditions expected in natural settings.

How to cite: Boneh, Y., Ohl, M., Plümper, O., Hirth, G., and Peč, M.: The Weakest Link – Revealing the microphysical deformation mechanisms of talc under P-T conditions associated with fault creep and slow slip events, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11449, https://doi.org/10.5194/egusphere-egu23-11449, 2023.

16:35–16:45
|
EGU23-9823
|
EMRP1.3
|
On-site presentation
Marianne Conin, Emilio AbiAad, Olivier Deck, and Jana Jaber

Morphology of rock joints (faults, fractures) has been recognized as the key factor controlling their mechanical behaviour, including the pre-peak and post peak phases as well as dilatancy. It also appears playing a significant role in two mechanical behaviours that have been observed on natural fracture during shearing: (i) joint dilation, or (ii) joint closure in association with asperities crushing, and rock matrix plastic deformation. We examine how (i & ii) occur in the joint, discussing their relationship with normal stress, joint morphology and intact matrix mechanical properties. To do this, two innovative methodologies based on 3DP technologies using a sand and phenolic binder on one side and a polymer (PA12) and binder jetting technology on the other one are applied to built fractures in a weak matrix, and in a strong matrix respectively. Joint surface roughness are built as fractal property with a self–affine replication, in accordance with natural observations. Results of direct shear tests under constant normal stress reveal that the mechanical behavior of the joints is first controlled by the mechanical parameters of the material (UCS/σn ratio), then by the joint geometry. In the case where the UCS/σn ratio is high (>40) corresponding to a strong material compares to the mechanical solicitation, no significant damage is notice on the joint and the maximal dilation angle is controlled by the steepest angles of the shorter wavelength asperities, which may only represent a small percentage of the surface roughness. In a the case of a weak material the joint behaviour is more complex, and is controlled by a specific range of asperities sizes. Three behaviours were observed depending on the applied normal stress: (i) at low normal stress the larger wavelengths asperities cause dilation since they are not sheared off; (ii) at normal stress over 40% of the UCS value, tensile and/or slip cracks were observed around those asperities, leading to their crushing and beheading; (iii) at normal intermediate stress, the two mechanisms were conjointly observed. In the second case (ii) joint closure is observed and the permeability increases in the surrounding matrix. Those results implies that the UCS/σn ratio plays an key role in fault shear behaviour, off-fault damage propagation and fluid circulation.

How to cite: Conin, M., AbiAad, E., Deck, O., and Jaber, J.: Investigation of fault behaviour during shear process for weak and strong materials using 3D printing technologies, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9823, https://doi.org/10.5194/egusphere-egu23-9823, 2023.

16:45–16:55
|
EGU23-5402
|
EMRP1.3
|
ECS
|
On-site presentation
Luc Illien, Jens M. Turowski, Christoph Sens-Schönfelder, Clement Berenfeld, and Niels Hovius

Earthquakes introduce long-lasting transient mechanical damage in the subsurface that can take years to recover to a new elastic steady-state. The associated transient perturbation of the elastic moduli can cause postseismic hazards such as enhanced landsliding.  This dynamics is linked to relaxation, a phenomenon observed in a wide class of materials after straining perturbations. In this study, we analyze the successive effect of two large earthquakes (the 2017 Mw7.7 Tocopilla and the 2014  Mw8.2 Iquique earthquakes) on ground properties through the monitoring of seismic velocity from ambient noise interferometry in the Atacama desert in Chile. The absence of rainfall in this area allows study of the mechanical state of the subsurface by limiting the potential effect of variations in groundwater content. We show that relaxation timescales are a function of the current state of the subsurface when perturbed by earthquakes, rather than ground shaking intensity. Our study highlights the predictability of earthquake damage dynamics in the Earth's near-surface and potentially other materials. We propose to reconcile this paradigm with existing physical frameworks by considering the superposition of different populations of damaged contacts. 

How to cite: Illien, L., Turowski, J. M., Sens-Schönfelder, C., Berenfeld, C., and Hovius, N.: Predictable healing rates in near-surface materials after earthquake damage in Chile, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5402, https://doi.org/10.5194/egusphere-egu23-5402, 2023.

16:55–17:05
|
EGU23-7857
|
EMRP1.3
|
On-site presentation
Wei Liu and Michel Jaboyedoff

The propagation of micro-cracks will cause the change of rock infrared radiation (IR) information, which provides the possibility to study the rock damage behavior and failure precursor using IR. In this paper, a new quantitative characterization method of rock damage evolution using IR is proposed. Firstly, the maximum classes square error and median filter methods are used to separate the temperature increment caused by crack development in IR images. On this basis, a new index, Damage Infrared Energy Response (DIER), is proposed to describe the crack evolution state and recognize the failure precursor of rock. It is found that the change characteristics of DIER and Acoustic Emission (AE) count are consistent: the DIER remains at the level in the compaction and elastic stages, rises gradually in the stable crack propagation stage, and increases sharply in the unstable crack propagation stage and fluctuates with the appearance of AE count “peak”. The change characteristics of DIER in unstable crack propagation stages can be regarded as the failure precursors of rock, about 84.10% of peak stress. Then, according to the continuum damage mechanics theory, the DIER is used to establish a theoretical characterization of the damage variable for rock, which can accurately describe the damage evolution process of the rock under uniaxial compression. The research results can provide experimental and theoretical support for monitoring slope and rock engineering stability by IR.

How to cite: Liu, W. and Jaboyedoff, M.: Theoretical damage characterization and failure precursor recognition of the rock under uniaxial compression using infrared radiation, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7857, https://doi.org/10.5194/egusphere-egu23-7857, 2023.

Modeling studies
17:05–17:15
|
EGU23-6256
|
EMRP1.3
|
ECS
|
On-site presentation
Marthe Grønlie Guren, Henrik Anderson Sveinsson, Anders Malthe-Sørenssen, Razvan Caracas, and Francois Renard

At the nanoscale, fracturing creates surface area and flow pathways, which control the rates of fluid-rock interactions. However, how fractures form at the nanoscale remains enigmatic. Here, we implement molecular dynamics simulations to reproduce fracture propagation in quartz and basalt. These simulations require large systems and long simulation times and are therefore currently depending on interatomic potentials. In the recent years, machine learning approaches have been established as a way to fit interatomic potentials, where the potentials are trained with quantum-mechanical data obtained from ab initio molecular dynamics simulations. We have developed machine-learned interatomic potentials for silica and basalt that allow using molecular dynamics simulations to simulate fracture propagation at the nanoscale. The interatomic potentials reproduce the mechanical properties of bulk silica and basalt sand have also been trained to account for fracture propagation. First, we trained a potential on silica to verify the fitting procedure, and then we used the same procedure to train an interatomic potential for basalt. By training the potential with water and carbon dioxide as fluids, we aim to study how a dynamic fracture damage basaltic glass and how the water and carbon dioxide enter these fractures in the wake of rupture. Our results are relevant for carbon mineralization where a coupling between dissolution of the basalt and precipitation of carbonate minerals can lead to nanofracturing of the rock.

How to cite: Guren, M. G., Sveinsson, H. A., Malthe-Sørenssen, A., Caracas, R., and Renard, F.: Machine-learned interatomic potentials for modelling nanoscale fracturing in silica and basalt, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6256, https://doi.org/10.5194/egusphere-egu23-6256, 2023.

17:15–17:25
|
EGU23-13293
|
EMRP1.3
|
On-site presentation
Fahimeh Najafi, Henrik Andersen Sveinsson, Christer Dreierstad, and Anders Malthe-Sørenssen

Neural networks have proved to be able to capture the relevant and informative features of a wide variety of data types and predict the desired output for different regression or classification problems. Finding a mapping between materials’ structure and a given physical property of those systems is an example of a problem that could be approached with machine learning methods like neural networks. Especially when we are dealing with systems with a very large design space where using classical computational methods like molecular dynamics can be very time and resource-consuming for the study of a very large number of systems, a well-trained neural network can be greatly faster and more efficient for computing the relevant properties. In this work, we study α-quartz crystals with one porous layer with simplex noise as the shape of porosity. Simplex noise is a gradient based procedural algorithm that can produce irregular geometries with surface morphology resembling what is observed in nature. The property that we want the neural network to learn is the yield stress of these systems under both shear and tensile deformation. Molecular dynamics simulations are used for a randomly selected sample of possible structures in order to generate the ground truth to be used as the training data. We employ deep convolutional neural networks (CNN) which are commonly used when dealing with image or image-like data since the input data for the problem in hand is a binary 2-D structure of the porous layer of the systems. The trained model is compared with a basic polynomial fit of stress versus porosity. The trained CNN is successful in predicting the yield stress of a system based on the geometry of that given system, with lower variability and higher precision compared to the base polynomial regression method. The saliency maps created with the trained model show the model to be successful in capturing the physics of the problem when compared with the stress fields calculated using molecular dynamics simulations. This method of modeling materials can be further developed for the inverse design of structures with desired properties without the need for a huge number of simulations on a wide domain of possible systems.

How to cite: Najafi, F., Andersen Sveinsson, H., Dreierstad, C., and Malthe-Sørenssen, A.: Modeling the relationship between mechanical yield stress and material geometry using convolutional neural networks, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-13293, https://doi.org/10.5194/egusphere-egu23-13293, 2023.

17:25–17:35
|
EGU23-12979
|
EMRP1.3
|
ECS
|
On-site presentation
Fengchang Bu, Michel Jaboyedoff, Marc-Henri Derron, and Lei Xue

The Universal Distinct Element Code (UDEC) has been rising much prevalence and is a validated technique applied to simulations in varied branches of geotechnics. Nevertheless, the characterization of water-weakening effects on a rock model remains elusive, resulting in a poor constraint referring to water-induced simulations by UDEC. In this context, previous research has been attempting to conduct an intuitive link between modelling parameters and saturation degrees, Sr, to implement a water-weakening process in UDEC, leading to a detrimental potential devoid of the basic logic that modelling parameters determine macroscopic properties of a rock model in contrast to dominance by Sr owing to the discrepancy in a physical sense and spatial scale. To fill in this gap, a new methodology coupled with a parametric study is first proposed with procedures that macroscopic properties of actual rock with different Sr are input into parametric relations to acquire predicted modelling parameters, which will be sequentially calibrated and adjusted until simulations are in line with actual tests. Utilizing this methodology, water-weakening effects on macroscopic properties, mechanical behaviours, and failure configurations of numerical models in UDEC are commensurate with tested ones to the utmost with noticeable computational expediency harnessing the benefit of the parametric study, indicating the feasibility and simplicity of the methodology. Thus, in the implementation of a water-weakening process in UDEC, we suggest converting an embarking from an intuitive link between modelling parameters and Sr into a parametric analysis to determine modelling parameters according to macroscopic properties under different Sr.

How to cite: Bu, F., Jaboyedoff, M., Derron, M.-H., and Xue, L.: Parametric study to characterize water-weakening effects in UDEC, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-12979, https://doi.org/10.5194/egusphere-egu23-12979, 2023.

17:35–17:45
|
EGU23-254
|
EMRP1.3
|
ECS
|
Virtual presentation
Özge Dinç Göğüş, Elif Avşar, Kayhan Develi, and Ayten Çalık

The deformation and failure processes of rocks under stress are primarily induced by microcracking. Detecting this micro-interaction phenomenon before the ultimate failure has paramount importance for predicting the post-failure rock damage characteristics. In this study, we aim to quantify the evolution of microcracking through fractal analyses of scanning electron microscope (SEM) images, captured from three different rock types subjected to uniaxial loading at various stress levels. In terms of uniaxial compressive (UCS) and tensile strength (UTS) values, the rocks range from the strongest to the weakest as being diabase, ignimbrite, and marble, respectively.  All rock samples are uniaxially loaded up to critical stress thresholds as crack initiation (σci), crack damage (σcd), and peak stress (σp) levels, considering their pre-defined characteristic stress-strain curves. Using the box-counting technique, the fractal dimension values (DB) of cracking intensity, induced by loading are determined for all these three stages. Here, it should be noted that higher fractal dimensions represent more intense microcracking according to the fractal theory. The results show that the DB values are increasing with the increasing amount of microcracks and the greatest DB values are calculated for Diabase due to its highest strength ratio (UCS/UTS). Although the marble has the weakest strength values, it presents a higher DB value than that of ignimbrite (DBmarble = 1.215 and DBignimbrite = 1.133) once the σcd stress threshold is reached. Furthermore, the DBmarble value is also greater than the DBignimbrite value for the σp stress level. It is because marble has a higher UCS/UTS ratio than the ratio of ignimbrite. Our results highlight the important role of rock texture on brittleness which exerts a primary control on fractal dimensions (DB). A decrease in volumetric rigidity is more dramatic in marble than in ignimbrite with incremental loading. The insights provide a better understanding of the microcracking process that leads to macro-scale deformations in rock engineering.

How to cite: Dinç Göğüş, Ö., Avşar, E., Develi, K., and Çalık, A.: Progressive failure characteristics of different rock types through fractal analysis, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-254, https://doi.org/10.5194/egusphere-egu23-254, 2023.

17:45–17:55
|
EGU23-4620
|
EMRP1.3
|
Virtual presentation
|
John Rudnicki

Recent industrial processes that involve injection of fluids, such as geothermal stimulation, disposal of waste water from hydraulic fracturing and carbon sequestration, have induced seismicity that has caused concern and resulted in discontinuation of the activity. Although field observations are the ultimate test of the effects of pore fluid on failure, their interpretation is complicated by heterogeneity of hydrologic and mechanical structure, and pumping and loading history. In particular circumstances, well-designed field tests can overcome some of these limitations. Laboratory experiments, despite their limited size and time scales, provide a more controlled environment that can yield an understanding of fundamental processes. Simple models that simulate the experiments can assess whether the mechanisms included in the models are sufficient to describe well the response or more complex formulations are needed. In addition, simulations can extend results for parameter values and loading programs beyond those achievable in experiments and aid in extrapolation to field applications.

This work uses a spring-block model and rate and state friction to simulate experiments conducted in a double direct shear apparatus on simulated carbonate fault gouge (Scuderi et al., EPSL, 2017) and on a shale bearing rock (Scuderi and Collettini, JGR, 2018). Both sets of experiments used the same loading protocol and injected pore fluid under creep conditions. When velocity strengthening rate and state friction is used to simulate the experiments on the simulated carbonate fault gouge the results agree well with the observed onset of tertiary creep in the experiment. Thus, the simulation reinforces the observation that pore fluid injection can induce rapid slip even when the friction relation is velocity strengthening. The rate and state framework provides an interpretation alternative to the standard one of the Mohr's circle moving to the left as pressure increases. In the rate and state framework, the friction coefficient must increase with pore pressure increase. The shale has a very low nominal friction coefficient (0.28) and is much more velocity strengthening than the carbonate. The simulation agrees with the observations that increases in pore pressure induce an increase in slip velocity but the magnitudes reach only about 100 µm/s by the end of the experiment. The simulation predicts reasonably well the times at which representative values of the slip velocity and displacement occur but the overall agreement of simulation and observation is not as good as for the carbonate. Mechanisms other than rate and state friction, for example, direct dependence of the friction coefficient on slip and porosity changes, may be significant.

How to cite: Rudnicki, J.: Rate and State Simulation of Two Experiments with Pore Fluid Injection Under Creep Conditions, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-4620, https://doi.org/10.5194/egusphere-egu23-4620, 2023.

17:55–18:00

Posters on site: Tue, 25 Apr, 10:45–12:30 | Hall X3

X3.14
|
EGU23-13723
|
EMRP1.3
|
ECS
Even Marius Nordhagen

Static friction is the force required to impose sliding on a rested body. The force depends on material properties and external factors such as normal pressure and temperature, but also a time dependent component is important. The frictional aging effect is at origination of the difference between static and dynamic friction, and is also believed to be responsible for the velocity-weakening of sliding friction. Despite immense effort, how microscopic processes affect the macroscopic aging is still not fully understood. We have performed molecular dynamics simulations where we demonstrate that high surface diffusion may provoke rapid contact area growth of an asperity-substrate interface, inducing a strong frictional aging effect. This mechanism differs from elastic and plastic creep in the sense that it occurs even at no normal pressure. The growth of contact area was found to be nearly logarithmic due to an exponentially decaying diffusivity. Furthermore, when applying a normal stress the aging effect is enhanced due to plastic creep. Our work suggests a new explanation of the logarithmic nature of aging and helps bridging the gap between empirical macroscale friction laws and the microscale behavior. While aging due to plastic and elastic creep is well-known and incorporated into most friction laws, diffusion aging has yet to be considered. The ultimate goal is to design or redesign friction laws taking the microscopic behavior into account and conceivably improve the accuracy of the laws. In the long term, this may contribute to improved earthquake forecast.

How to cite: Nordhagen, E. M.: Atomic scale frictional aging in silicon carbide due to diffusion and creep, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-13723, https://doi.org/10.5194/egusphere-egu23-13723, 2023.

X3.15
|
EGU23-10164
|
EMRP1.3
|
ECS
Derek Leung, Florian Fusseis, and Ian Butler

Fatigue and damage accumulation in granitoids are classical, but poorly characterised, rock mechanics problems. In order to explore these phenomena, we consider colliding curling stones as a rock physics experiment. Curling stones are made using granitoids from either Ailsa Craig (Scotland) or Trevor (North Wales), which are chosen for their uniformity, strength, and durability. During a curling game, stones are slid over an ice sheet and made to collide along a circumferential striking band. From a rock physics perspective, the collision of curling stones can be modelled as unconfined uniaxial compression of two convex surfaces under well defined boundary conditions. A curling stone experiences about 2900 collisions per season and is played for 10-15 years before refurbishment, which provides a unique long-term opportunity to study fatigue and damage accumulation under cyclic loading.

Here, we first determine the stress magnitudes and strain rates of head-on curling stone impacts using a series of on-ice experiments involving a high speed camera and pressure-sensitive films. We then characterise the observed damage that these collisions produce on the centimetre and micrometre scale using photogrammetry, synchrotron microtomography, optical microscopy, and backscattered electron imaging. We show that during each impact, a curling stone is stressed to at least 300-680 MPa (for a maximum-velocity scenario of 2.9±0.1 ms-1), which exceeds the unconfined compressive strength of the rocks (232-395 MPa; Nichol, 2001, J. Gemm. 27/5). Over its lifetime, a curling stone thus experiences thousands of impacts that will cause damage. The strain rates of these impacts (24±4 s-1) most closely resemble seismic magnitudes, suggesting that the impacts are dynamic in nature. This is supported by the type of damage observed in aged curling stones: (1) Hertzian cone fractures, (2) ejection of rock powder during collisions, and (3) minor spalling microcracks. Most samples show damage being confined to macroscopic Hertzian cone fractures and their immediate collet zones in the relatively narrow striking band. Within the striking band, the circumferential density of cone fractures is limited to about 2-2.5 fractures/cm. Surprisingly, damage does not appear to extend beyond about 3-5 cm into the stones along a radial direction.

Our observations allow us to formulate a model for damage evolution in curling stones. We infer that high-velocity/high-stress impacts initiate cone fractures up to a specific spatial density. As they mature over repeated impacts in the same regions of the striking band, cone fractures progressively propagate and coarsen with subsequent collisions, concentrating and channelling the accumulating damage. This damage geometry is surprisingly effective in shielding the rest of the stones from the reaching critical stress levels for damage. Our findings are significant for applications where rocks are exposed to large numbers of high-stress impacts and suggest that a relatively narrow damage zone can dampen even high-impact stresses over a relatively moderate network of fractures.

How to cite: Leung, D., Fusseis, F., and Butler, I.: Where curling stones collide with rock physics: Cyclical damage accumulation and fatigue in granitoids, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10164, https://doi.org/10.5194/egusphere-egu23-10164, 2023.

X3.16
|
EGU23-8188
|
EMRP1.3
Sergio Vinciguerra, Thomas King, Guido M. Adinolfi, and Philip Benson

Parametric analysis of laboratory Acoustic Emission (AE) during rock deformation laboratory experiments has revealed periodic trends and precursory behaviour of the rupture source, as crack damage nucleates, it grows and coalesces into a fault zone. Due to the heterogeneity of rocks and the different effective pressures, finding a full prediction of rupture mechanisms is still an open goal.

4x10cm cylindrical samples of Alzo granite were triaxially deformed at confining pressures of 5-40 MPa, while AE are recorded by an array of twelve 1MHz Piezo-Electric Transducers. AE are then post-processed to derive attributes and parameters. We aim to identify what are our most important parameters, and more interestingly, when they are most relevant for predicting when the rock will fail.

Time Delay Neural Networks (TDNN) have shown promise in forecasting failure when using AE-derived parameters. We trained a TDNN with 5 key parameters: 1) AE event rate, i.e. the number of events obtained during the incremental deformation (strain); 2) AE amplitude, i.e. maximum amplitude of S-waves, 3) AE source mechanisms inferred by the source radiation patterns to categorize events and obtain source orientations of mixed-mode type mechanisms; 4) Seismic scattering, i.e. the ratio between the low frequency (LF, 50-500 kHz) and high frequency (HF, 500-1000 kHz) peak delay (PD) values for individual AE and 5) Bulk elastic S-wave velocity measured at intervals throughout the experiment along the ray-paths created by transmitters and receivers. As each parameter investigates a specific mechanical aspect, taken together they provide information on deformation, fracturing and the evolving state of the background medium as failure is approached. These timeseries are then classified by the TDNN as variations in stress and strain (target parameters).

We are currently assessing the importance of individual parameters by omitting one at a time from the training routine. The more important the omitted parameter, the larger the misfit will be when comparing the network output and the target timeseries. The omission analysis determines what are the most important parameters to use when training a neural network to predict dynamic failure. Results are strongly dependent on the methods used to define the training parameters, but several trends are emerging. Event rate and amplitude differently influence predictions of stress and strain. Event rate appears relevant only in the early deformation phases, while amplitude seems much more significant during the coalescence/propagation phase. Seismic scattering and source mechanisms also show an early relevance, interpreted as due 1) to the breakup of low frequency surface waves as microcracks begin to coalesce and 2) bursts of tensile events in the enucleation phase and an increase at ~80% UCS, likely related to the crack propagation. Similarly, there is a clear pivot in the importance of seismic velocity during the early stage, but it emerges a progressive increase ~40% UCS whose origin is unclear. We are currently determining if these variations are directly related to the mechanics of the fault zone or are simply an artifact of the processing.

How to cite: Vinciguerra, S., King, T., Adinolfi, G. M., and Benson, P.: Using AE based Machine Learning Approaches to Forecast Rupture during Rock Deformation Laboratory Experiments, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8188, https://doi.org/10.5194/egusphere-egu23-8188, 2023.

X3.17
|
EGU23-2397
|
EMRP1.3
Patrick Baud, Fanbao Meng, Lingcao Huang, and Teng-fong Wong

Understanding compaction localization in porous limestone in the laboratory is significantly more challenging than in sandstone because of the lack of consistent acoustic emission activity in carbonate samples. Previous studied have therefore relied on X-ray Computed Tomography imaging (CT). The first unambiguous evidence of compaction band development in limestone was provided by Huang et al. (2019), who performed synchrotron in situ CT imaging during shear-enhanced compaction in a sample of Leitha limestone. This sample was deformed in the HADES rig at the European Synchrotron Radiation Facility, in dry conditions and at a confining pressure of 20 MPa. In this study, we analysed this data set using Digital Volume Correlation (DVC). Not only could we use DVC to characterize quantitatively the spatiotemporal development of displacement and strain, we were also able to compare with direct observations to assess the stress-induced damage in multiple scales. Our new results confirm that inelastic compaction occurred in two stages in Leitha limestone: macropore collapse first and then sequential growth of compaction bands. In the pore collapse stage, DVC reveals complex and heterogeneous grain-scale strains, implying significant heterogeneity in the internal stress field. Such complexity is to be accounted for if one were to connect micromechanical and continuum models. At higher stresses, we have obtained further quantitative constraints on the spatial distribution of volumetric and shear strain during the growth of compaction bands. Our results demonstrate that compaction banding in Leitha limestone can be analysed as a bifurcation phenomenon, that would typically occur preferentially in zones of high porosity. The displacement field inferred from DVC revealed that the bands showed mostly normal displacement discontinuities, as expected for compaction bands. DVC analysis also gave more constraints on band geometric attributes. Analysis of the autocorrelation function for the strain suggested that the decay and rebound of the autocorrelation as a function of the axial separation may provide proxies for the mean width and spacing of compaction bands. The 2D autocorrelation function on the band planes also provides relevant clues on the complex sequential growths of the compaction bands.

How to cite: Baud, P., Meng, F., Huang, L., and Wong, T.: Compaction localization in 4D imaged by X-ray Computed Tomography and Digital Volume Correlation, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2397, https://doi.org/10.5194/egusphere-egu23-2397, 2023.

X3.18
|
EGU23-6624
|
EMRP1.3
|
ECS
Thomas Alcock, Sergio Vinciguerra, and Phillip Benson

Stromboli volcano, located in the north-easternmost island of the Aeolian archipelago (Southern Italy) and well known for its persistent volcanic activity, has experienced at least four sector collapses over the past 13 thousand years. The most recent activity resulted in the formation of the Sciara del Fuoco (SDF) horseshoe-shaped depression and a tectonic strain field believed to have promoted flank collapses and formed a NE / SW trending weakness zone across the SDF and the western sector of the island. The tectonic strain field interplayed with dyking and fracturing appears to control the episodes of instability and the onset of slip surfaces. This study presents new data identifying areas of damage that could promote fracturing via remote sensing and rock friction measurements taken on rocks around the SDF and the coupled “weak” zone. We have carried out a multiscale approach by integrating satellite and microscale observations with frictional tests carried out in triaxial configuration on cm scale slabs.

 

Key units have been sampled on the field (Paleostromboli, Vancori and Neostromboli) with reference to SDF and the weak zone. Direct-shear tests in triaxial configuration were carried out to explore the frictional and seismic properties using rectangular basalt slabs at 5 – 15 MPa confining pressure in dry and saturated conditions, while recording acoustic emissions (AE) via two Piezo-Electric Transducers. The sliding velocity was changed to acquire rate and state friction parameters (RSF). Preliminary results show a variation in the friction coefficient (m) between 0.55 and 0.9 with a general m decrease with increasing confining pressure and saturation. RSF parameters a-b (0.1 < a-b < 0.1) and steady state friction coefficient (mss) (0.6 < mss < 0.9) are controlled by changing sliding velocity, confinement and by the physical properties of each unit, in particular the porosity.  AE key attributes, such amplitude, frequency and duration and their evolution confirm the relation to sliding velocity, confinement and porosity. Ongoing post mortem SEM analysis are aiming to assess the impact that textural features, such porosity, crystal distribution and glass groundmass for the different units have on the evolution of crack damage and their control on the frictional properties. Quantitative crack density analysis will be carried out using the Matlab tool box FracPaQ on the microstructures to quantify fractures properties and highlight which mechanical features (for example crystals or pores) control the development of asperities/stress concentration. This finding can be related to the field scale fracture density analysis, providing quantitative support for the identification of structurally weak zones across the SDF and constraint the mechanical behaviour of the fractured zones prone to instability.

 

How to cite: Alcock, T., Vinciguerra, S., and Benson, P.: Multiscale analysis of physical rock properties at Stromboli Volcano: what controls the frictional properties?   , EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6624, https://doi.org/10.5194/egusphere-egu23-6624, 2023.

X3.19
|
EGU23-1727
|
EMRP1.3
Wen-Bin Doo and Hsueh-Fen Wang

A zone of significant high-amplitude magnetic anomalies is observed without a comparable gravity high along the Cascadia margin and is spatially correlated with the low-velocity fore-arc mantle wedge, which is understood to be serpentinized fore-arc mantle and is further considered to be the main source of the high-amplitude magnetic anomalies. To test this concept, the magnetization-density ratio (MDR) is estimated along the Cascadia margin to highlight the physical characteristics of serpentinization (reduce density and increase in magnetization). Interestingly, high MDR values are found only in central Oregon, where slab dehydration and fore-arc mantle serpentinization (50%-60% serpentinization) are inferred in conjunction with sparse seismicity. This result may indicate either a poorly serpentinized fore-arc mantle or that the fore-arc mantle is deeper than the Curie temperature isotherm for magnetite in northern and southern Cascadia. This finding means that magnetic anomaly highs and serpentinized fore-arc mantle may not be completely positively related in subduction zones.

How to cite: Doo, W.-B. and Wang, H.-F.: Relationship between the high-amplitude magnetic anomalies and serpentinized fore-arc mantle in the Cascadia subduction zone, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-1727, https://doi.org/10.5194/egusphere-egu23-1727, 2023.

X3.20
|
EGU23-7718
|
EMRP1.3
Chung-Liang Lo, Shu-Kun Hsu, Shiao-Shan Lin, Ching-Hui Tsai, Wen-Bin Doo, Song-Chuen Chen, and Pin-Ju Su

As an active mountain building, Taiwan Island is deduced from the oblique collision between the Eurasian and the Philippine Sea plate, and its northern part and offshore region are under post-collision collapse. The magmatism induced from the post-collision collapse therefore distributes in the northern Taiwan and its offshore area. A series of submarine volcanoes and igneous rock isles are rooted in the area. For investigation of the volcanic and igneous arrangement, we have collected the magnetic data over the past few decades to combine and compile a map of regional magnetic anomalies. A pronounced magnetic high largely dominates the area of most submarine volcanoes and extends eastward, while the adjacent areas to the north and west are lower. To better understand the magnetic features for the submarine volcanic area, the magnetization for an equivalent magnetic layer thickness was calculated. The result shows that a high magnetization concentrated on the SV7 and extends northwestward that could be a magnetic dipole combining with its northeastern low part. To the southward, the submarine volcanoes SV1, SV3-SV6 locate between this high and another low magnetization. We also applied the enhanced analytic signal technique from the same magnetic data to evaluate the magnetic source strength distribution. Except for the SV2, SV5 and SV6, the 0th degree of enhanced analytical signal shows that most signal high concentrated on the submarine volcanic areas. For higher degree of enhanced analytical signal, the highest magnitude focus on the Ks and SV1, and PV, SV3 and SV4 are slightly minor.

How to cite: Lo, C.-L., Hsu, S.-K., Lin, S.-S., Tsai, C.-H., Doo, W.-B., Chen, S.-C., and Su, P.-J.: Geomagnetic characteristics of submarine volcanic area off the coast of northern Taiwan, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7718, https://doi.org/10.5194/egusphere-egu23-7718, 2023.

X3.21
|
EGU23-10329
|
EMRP1.3
A new method for identifying and evaluating bitumen-bearing reservoirs based on conventional logging curves
(withdrawn)
Jue Hou, Lun Zhao, Xing Zeng, Xueke Wang, Yiqiong Zhang, and Congge He
X3.22
|
EGU23-10226
|
EMRP1.3
|
ECS
Che-Ming Yang, I-Lin Chuang, and En-Lun Zhang

Deep-seated landslides may become long-term creeping or transform to catastrophic landslides. Both serious threat to mountainous roads, villages, tourist area, and reservoir areas, which belongs to long-term and extensive effects. Many historical catastrophic landslides have caused the devastating disasters, such as Tsaoling landslide induced by the 1999 Chichi Earthquake, the Hsiaolin village landslide induced by extreme rain of the 2009 Typhoon Morakot, and the long-term large-scale landslides (creeping) of Lushan and Lishan have affected on the environment of the adjacent areas for decades or even longer. In Taiwan, there are many regions of widespread dip slope landform with potential planar failures, study area ranges the right bank of Chishan River from south to north (~35km in length). Most potential planar failure areas had delineated, and several platforms or gentle surfaces on the slope represent the deposits of paleo planar slides or old landslides in study area. But few cases of the geomorphologic evolution are investigated. However, it is difficult to estimate the slope stability of potential planar failures without geomorphologic evolution model. Therefore, the main purpose of this study is to use field investigation and topographic analysis to establish engineering geological models, then propose the geomorphologic evolution model for evaluation of slope stability. The methods of this study include: (1) identifies the microstructure of landslides by high-resolution LiDAR data, (2) performs the geological investigation to verify topographic interpretation and records occurrences of outcrops, (3) collects historical orthoimages to evaluate the activity of slope, (4) use high-precision aerial photography to establish digital surface model and analysis point cloud data to obtains the discontinuous plane state, and (5) the failure mechanism would be analyzed by the stereographic projection. In study area, 160 platforms are identified and area ranges from 220 m2 to 82386 m2. The largest two platforms, Tianziding (TZD) and Mujiliao (MZL) platforms, are investigated by field survey and drone. The attitudes of interbeds under TZD platform are surveyed along gullies, which is obvious gentler than the strata of dip slope. The front of MZL platform occurred slope failure during 2016 Typhoon Megi, therefore, the exposed rock mass of platform can be identified as deposits of planar slide. The preliminary results can infer that numerous paleo planar slides are exist in study area and the geological profiles of TZD and MZL are plotted. The structures of rock masses and attitudes of discontinuities by field survey and point cloud analysis need to interpreted carefully. Then, the geological model and geomorphological evolution can be proposed. 

How to cite: Yang, C.-M., Chuang, I.-L., and Zhang, E.-L.: Re-failure and geomorphological evolution of paleo planar slide, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10226, https://doi.org/10.5194/egusphere-egu23-10226, 2023.

X3.23
|
EGU23-14394
|
EMRP1.3
|
ECS
|
Highlight
|
Luigi Massaro, Federica Rauseo, Melania De Falco, Ermanno Marino, Giovanni Forte, and Antonio Santo

Ischia (Southern Italy) is a volcanic island of the Phlegrean Volcanic District that was recently affected by multiple geological hazards, including floodings, landslides, rockfall and earthquakes.

In this study, rockfall stability is analysed, assuming as a case study a 400m-wide cliff made of Green Tuff and located on the western area of Mt. Epomeo. The two outcrops studied are located at 280 and 420 m a.s.l., above the site of Frassitelli, Forio d’Ischia, which is an area of high residential, touristic and agricultural importance. The former is a high-angle outcrop affected by tens of meters-long faults, whereas the latter is characterised by high-dip pinnacles.

We analysed the fracture systems affecting the examined formation to compute the kinematic analysis of the potential rupture mechanisms and to perform numerical simulations of potential rockfall scenarios. The data acquisition was carried out by means of classical geological field surveys and structural analysis on Virtual Outcrop Models (VOM) obtained from images acquired by drones. The VOMs were analysed with ‘CloudCompare v2.10.2’ and ‘OpenPlot’ software. The former allowed the automatic digitalisation of the exposed discontinuities by applying the ‘Facets’ plugin, based on a least-square fitting algorithm (Fernández, 2005). ‘OpenPlot’ enabled the extraction of the geostructural information from the VOM, by computing the best-fit planes of the polylines manually drawn along the interference between the geological surface and the outcrop topography (Tavani et al., 2011).

The measured and the extracted features were classified following their attitude. Three main sets were defined, striking N-S, NW-SE and NE-SW. The fracture dataset was used to perform a kinematic analysis with ‘DIPS’ software on the surface discontinuities extracted from ‘Facets’ plugin. The 'wedge sliding' resulted the most critical potential rupture mechanism to occur on the analysed outcrops. Successively, numerical simulations of rockfall scenarios were computed based on the acquired structural information. The latter permitted us to identify the maximum run out of the potential blocks and draw some consideration on the rockfall hazard of the area.

How to cite: Massaro, L., Rauseo, F., De Falco, M., Marino, E., Forte, G., and Santo, A.: Characterization of rockfall mechanisms and run-out in active volcano-tectonic areas: a case study from Ischia Island, Southern Italy, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14394, https://doi.org/10.5194/egusphere-egu23-14394, 2023.