Quantifying the changes in elastic properties of rocks during deformation is an important task. Effective Medium Theory (EMT), as formulated by Sayers & Kachanov (1995) relates the crack fabric (or damage) to the elastic properties. EMT has been successfully applied in the forward sense to predict the evolution of elasticity and related acoustic velocities in response to prescribed changes in crack density; and in the inverse sense to recover crack densities from laboratory measurements of acoustic velocities.  However, EMT fails to predict an important observation from laboratory studies of rock deformation: cyclic loading under uniaxial and conventional triaxial loads of rock samples can produce significant increases in Poisson’s ratio. These increases correlate with increasing number of cycles and with increasing crack density. This phenomenon has been known since the work of Walsh (1965), Brace et al. (1966) and Zoback & Byerlee (1975). More recent work by Heap & Faulkner (2008) and Heap et al. (2009; 2010) has extended the findings across a range of different lithologies.

Published EMT equations predict Poisson’s ratios that stay constant or decrease with increasing crack density. Resolving this discrepancy is important because Poisson’s ratio may play a key role in producing stress rotations in the damage zones of faults, thereby making them ‘weak’ and prone to slip even when the normal stress is high e.g. the San Andreas Fault (Faulkner et al., 2006; Healy, 2008). Building on the work of David et al. (2012 & 2020) incorporating the effects of crack closure, sliding on cracks (Kachanov, 1992) and grain boundaries (Sayers, 2018) during loading, and delayed back-sliding during unloading, closed form micromechanical equations have been derived to describe increases of Poisson’s ratio with increasing number of cycles. Critically, increases in Poisson’s ratio are predicted even without including the effects of new cracks. Examples are shown comparing the predicted changes in Poisson’s ratio using the newly derived equations to data from uniaxial and triaxial laboratory tests on cracked rocks.

How to cite: Healy, D.: Increases in Poisson’s ratio due to cyclic deformation in cracked rock, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3893, https://doi.org/10.5194/egusphere-egu24-3893, 2024.

Understanding the damage processes in clay-bearing rocks is a decisive factor in geological engineering, and for instance considering nuclear waste deep geological repositories. But, more generally they may also contribute to localized deformation, and thus the rupture of fault gauges in seismic zones. However, owing to their complex mineralogy, multiscale microstructures and anisotropy, the mechanisms of clay-rich rock damage and their chronology are not yet well understood..

Here we focus on the impact of micro-damage on ultrasonic wave propagation velocity, which is confronted with the corresponding full deformation fields calculated by digital image correlation (DIC). 

The aim is to associate the acoustic signature with the active deformation mechanisms identified by DIC. To this end, an integrated experimental approach is proposed to  characterize localization and to identify the related deformation micro-mechanisms  during uniaxial compression of natural clayey rock samples (Tournemire shales) with two simultaneous measurements: 1) the evolution of P-wave velocity within the sample by active acoustics, 2) the development of the 2D mechanical full field by digital image correlation.

Both experimental techniques are well known, but the innovation of our approach is to combine simultaneously both measurements. Deformation localization is a multiscale problem, which obviously occurs at the sample scale, but also at the fines scales of the microstructure. Therefore, we developed two different experimental setups. On the one hand, during uniaxial compression with a standard MTS loading frame the macro-scale localization patterns are characterized by optical observations, which image resolution is well suited to the cm sample scale (sample diameter: 3.6 cm and double in length). On the other hand, in order to characterize the initiation of micro-damage at the microstructure scale of the composite type of rock, the same loading protocol is reproduced (while keeping the acoustic diagnosis) on smaller scale mm-sized specimens (sample diameter : 8 mm, double in length), using a home-designed miniature loading frame fit for an environmental scanning electron microscope (ESEM). The latter analysis is carried out under  controlled relative humidity  of RH = 80%, hence preventing the samples to dry out due to the high vacuum

A similar acoustic signature is identified at both scales of observation, in spite of the variations of experimental conditions imposed by the environmental SEM. We are therefore confident to be able to understand the fracturing process from micro-cracking initiation (microscale) to sample failure (macroscale), and to assess its impact on ultrasonic wave propagation.

How to cite: Lusseyran, M., Dimanov, A., Bonnelye, A., Fortin, J., and Tanguy, A.: Multi-scale experimental deformation and damage initiation of clay-rich rocks  : Coupling ultrasonic wave propagation and  full field deformation measurements by digital image correlation (DIC), EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9811, https://doi.org/10.5194/egusphere-egu24-9811, 2024.

In the context of the global environmental crisis and the urgent need for energy transition and efficient energy storage solutions, salt caverns have gained attention as promising reservoirs for hydrogen. However, current literature predominantly focuses on deriving macroscopic constitutive relations, lacking crucial insights into the underlying physical mechanisms of deformation and damage active at various microscopic scales. This study addresses this gap by undertaking qualitative and quantitative investigations into the micro-mechanisms of rock salt, employing advanced micro-scale observation techniques. Natural rock salt from diverse mines and re-synthetic salts, produced through the cold compaction of grinded natural halite powder, are used to encompass a wide range of microstructural morphologies. Initial microstructure characterization involves SEM, EBSD, and CT, followed by classic uniaxial compressive tests coupled to optical microscopy monitoring. High-resolution images of the sample surface are continuously captured during testing, allowing for 2D full field measurements by subsequent application of digital image correlation techniques : the analysis of relative displacements of markers randomly distributed on the sample surface enables the retrieval of surface displacement fields and the calculation of the corresponding local strain fields over statistically representative domains. Segmentation of digital images and quantitative identification, specifically focusing on crystal slip plasticity and grain boundary sliding using an in-house computation program, reveal the complex local interactions of different micro-mechanisms. The estimation of the relative contributions of these mechanisms to global deformation all along the loading path, along with an analysis of the impact of salt grain size, provides insights into physically grounded micromechanical constitutive relations. These findings are essential for the safety assessment of industrial applications involving rock salt caverns with respect to short-term mechanical loading conditions relevant to daily hydrogen filling and withdrawal.

How to cite: Li, X., Dimanov, A., Bornert, M., Hallais, S., and Gharbi, H.: Multiscale experimental investigation of crystal plasticity and grain boundary sliding in rock salt using digital image correlation, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10099, https://doi.org/10.5194/egusphere-egu24-10099, 2024.

Recent observations of large earthquakes document the progressive localization of rock damage around future rupture zones that is also coupled with the spatial migration of foreshock sequences (Kato & Ben-Zion, 2020). This implies that the precursory deformation may act as a potential tracer for preparatory process that result in large earthquakes. It has also been observed that self-organization of the localized damage regions can govern the eventual macroscopic brittle failure in geomaterials (Renard et al., 2019). How the presence of fluid controls the self-organized precursory deformation along localized damage zone remains an open question. In this study, we have performed two triaxial compression experiments on dry and water saturated Berea sandstone, using distributed strain sensing (DSS) technology to visualize the strain field on the sample surface (Salazar Vásquez et al., 2022) with high spatial resolution. By tracking components of the strain field, specifically the region on the sample that sustained the largest incremental change in strain, we tested the effect of fluid on the predictability of phase transition between intact and failed state, under the context of critical hypothesis. Strain was progressively localized around the eventual faulting region for both samples, while a slow faulting was observed in the wet sample accompanied by a diffuse deformation pattern and unstable crack nucleation at failure. The results showed that, the failure in the dry sample was preceded by a critical power law acceleration of the largest increment, thus the dynamic faulting occurred in a well-defined singularity. The strain distribution also provided evidence for a predictable evolution of precursors. In contrast, the wet test showed evidence for a first-order transition with an exponential increase in largest increment, leading to an abrupt failure with a transient increase of strain. We interpreted this abrupt transition to be due to the increasing dominance of fluid-driven subcritical crack growth in the faulting. In this process, the local stress at crack tips decreases with crack lengthening, hence impeding the crack interaction and leading to an abrupt development of fault network. Our observation unravels the mechanisms of precursory deformation with fluid-assisted subcritical cracking, which has important implication in forecasting large earthquakes in nature.

References:

Kato, A., & Ben-Zion, Y. (2020). The generation of large earthquakes. Nature Reviews Earth & Environment, 2(1), 26–39. https://doi.org/10.1038/s43017-020-00108-w

Renard, F., McBeck, J., Kandula, N., Cordonnier, B., Meakin, P., & Ben-Zion, Y. (2019). Volumetric and shear processes in crystalline rock approaching faulting. Proceedings of the National Academy of Sciences, 116(33), 16234–16239. https://doi.org/10.1073/pnas.1902994116

Salazar Vásquez, A., Rabaiotti, C., Germanovich, L. N., & Puzrin, A. M. (2022). Distributed Fiber Optics Measurements of Rock Deformation and Failure in Triaxial Tests. Journal of Geophysical Research: Solid Earth, 127(8). https://doi.org/10.1029/2022JB023997

How to cite: Chen, H., Selvadurai, P., Salazar, A., Bianchi, P., Michail, S., Rast, M., Madonna, C., and Wiemer, S.: The influence of fluid pressure on the phase transition of brittle faulting, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13010, https://doi.org/10.5194/egusphere-egu24-13010, 2024.

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. We trained a Time Delay Neural Networks (TDNN) on key seismic attributes derived from AE, such: Event rate, i.e. the negative log time difference between successive events; Amplitude, i.e. the average max amplitude of all waveforms for each single event AE; Source mechanism estimated from first-motion polarity spheres (King et al., JGR, 2021); Seismic scattering, i.e the ratio between high and low frequency peak delays (King et al., GJI, 2022); Vp/Vs ratios from vertical P-wave velocities and horizontal S-wave velocities for individual AE (King et al., GJI, 2023).

These timeseries are then classified by the TDNN as variations in stress and strain (target parameters). TDNN require continuous, regularly sampled data but AE are discrete and irregular. To transform the training data for the TDNN, parameters are smoothed in a weighted moving window of 100 AE events, where weighting is given towards high amplitude events that occur close in space together. Data processing is applied to waveform data from all experimental condition. Despite the inherent complexity in the raw data, clear increasing or decreasing trends are repeated at different experimental conditions.

Hyperparameters for the neural network are optimised using a Genetic Algorithm (GA) by evaluating the misfit between training target (mechanical data) and model output. Each model is trained on the 10 MPa dataset and validated on the 40 MPa dataset. Roles are reversed and the results summed. This approach ensures consistent trends in the training data (waveform parameters) whilst reducing bias towards a particular dataset. We then investigate 120 configurations for the training data following a ‘leave-one-out’ strategy. E.g., a model is trained on 5, 10 and 20 MPa datasets whilst omitting the Event rate parameter. The model is then validated on the 40 MPa dataset.

Model output on validation datasets demonstrate that the TDNN can classify AE-derived parameters as increasing variations in stress and strain. 10 and 40 MPa demonstrate the best fit and are likely linked to the GA optimisation, highlighting biases driven by the training data. Forecasting results for strain and stress reveal notable over- and under-estimations of values. However, both 10 and 40 MPa are generally accurate to within 20% further highlighting the feasibility of using a TDNN for forecasting the development of new fracture under conventional triaxial conditions.

How to cite: Vinciguerra, S., King, T., and Benson, P.: Using AE based Machine Learning Approaches to Forecast Rupture during Rock Deformation Laboratory Experiments, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5440, https://doi.org/10.5194/egusphere-egu24-5440, 2024.

Lineaments are elongated elements in spatial data such as valleys and ridges on topographic maps, or linear lows and highs in aeromagnetic data. Topographic linear depressions (topolineaments) are generally considered as the morphological expressions of easily erodable, elongated rock bodies situated within a mechanically stronger rock mass. In most circumstances topolineaments are even directly interpreted as faults and fractures, which forms the basis for lineament analysis study to understand brittle deformation patterns. However, topolineaments may also be formed by other tectonic and non-tectonic causes, such as alternating layers, foliation traces, dikes etc., or river- and glacial erosion not controlled by any bedrock features. This mix of potential causes begs the question: “How robust is actually lineament analysis for characterizing and quantifying faults and fractures?”. Testing the topolineament vs. fracture-relationship is not straight forward since topolineaments are usually occupied by rivers or creeks and covered with colluvium, which prevents direct observation of rock types and bedrock structures.

Underground excavations allow for continuous logging of bedrock types, rock mass quality and fracture density and orientations, which is done routinely at tunnel face during tunnel construction. Here we make use of such underground data from a large dataset of Norwegian road tunnels to compare the position of topolineaments spatially and statistically with rock fracture density and orientations in the subsurface. The tunnel dataset comprises data from across Norway in areas with widely varying bedrock geology, tectonic evolution, and geomorphology (e.g. etched surface, alpine, lowlands), which allow for an evaluation of the robustness of lineament analysis in various settings. Topolineaments are acquired using a newly developed algorithm (OttoDetect) run on both 10x10m and 50x50m resolution digital elevation models. The algorithm ensures that tunnel data is compared to a homogeneous and reproducible lineament dataset without operator or hillshade illumination biases.

Preliminary results from tunnels in areas with etched geomorphology show that c. 75% of all topolineaments correspond to weakness zones in the bedrock (i.e. very high fracture densities/very low rock mass quality compared to the surroundings). Hit rate increases for longer lineaments, which generally correspond to thicker fault zones. At the same time, only up to c. 60% of all weakness zones mapped at tunnel face can be spatially associated to a topolineaments, which demonstrate that significant brittle deformation is not expressed as topolineaments. Further analysis will be carried out to build a statistically robust dataset on the validity of lineament analysis.

How to cite: Torgersen, E., Arctander, K., Redfield, T. F., Svendby, A. K., Dichiarante, A. M., and Arntsen, M. L.: Are all lineaments the surface expression of faults and fractures? – A novel analysis using tunnel face mapping data from Norwegian road tunnels, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16840, https://doi.org/10.5194/egusphere-egu24-16840, 2024.

Rock weathering is ubiquitously observed at or near Earth’s surface as a fundamental component in many landscape evolution process. In arid landscapes, where limited moisture availability restricts the rate and effectiveness of chemical and biological weathering – salt weathering (regarded herein as the physical disintegration of rocks in the presence of salts) is commonly acknowledged as an especially effective mechanism for progressive weathering of rocks. While volumetric expansion and contraction of salts in response to changes in ambient moisture conditions are broadly recognized as the primary drivers of salt weathering, our understanding of the environmental conditions that produce such moisture dynamics in otherwise extremely dry settings, such as hyper-arid deserts, remains largely unknown.

Here, we present preliminary results from field-based acoustic emission (AE) measurements for boulders with salt-laden cracks perched on abandoned shorelines of the hypersaline Dead Sea. Continuous measurements since April 2023 revealed daily fracturing activity displaying a bi-modal distribution with AE activity peaks during the early predawn and afternoon hours when T changes are minimal and RH fluctuations reach maximum or minimum values, respectively. Time-lapse photography revealed a recurring pattern of salts that crystalize along the rock cracks during the afternoon AE peak hours and subsequently disappear towards the predawn AE peak hours. The appearance of salt crystals during lowest RH conditions (warmest afternoon hours) and their disappearance during highest RH conditions (coldest predawn hours) suggests that stresses induced by repeated cycles of salt deliquescence/efflorescence in response to daily fluctuations in atmospheric RH are most likely responsible for the bi-modal distribution of daily fracturing activity. This suite of new field-based measurements of salt-weathering activity in natural hyper-arid settings suggests that atmospheric RH fluctuations and the volumetric changes they induce in hygroscopic salts can be key drivers of progressive rock fracturing in extremely dry and salt-rich environments on Earth and possibly Mars where other moisture sources are limited to effectively non-existent.

How to cite: Mushkin, A., Boroda, R., Malik, U., Lensky, N., Haggai, E., Muravin, B., and Amit, R.: Salt-driven Progressive Fracturing of Alluvial Boulders Along the Hyper-arid Shoreline of the Dead Sea , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14362, https://doi.org/10.5194/egusphere-egu24-14362, 2024.

Time- and stress-dependency of elastic properties are features particularly observed in a variety of complex solids, ranging from steel, polymers, and cracked structures, to rocks and concrete. Recently, considerable effort has been made to understand the underlying physics of these phenomena commonly regarded as Nonlinear Mesoscopic Elasticity (NME) in laboratory setups.

As a result, various models have been suggested to explain a range of NME phenomena like hysteresis, dynamic softening, and slow dynamics, among others. Due to the high sensitivity of NME to the presence of imperfections or internal damage on solids, there is a growing interest in taking the current models and applying them to construction materials for damage assessment.

Intending to observe and incorporate these models into real-condition structures, we carried out a 1-day multifrequency vibration experiment in a 24-meter-long reinforced concrete test bridge equipped with a pretension system, to investigate the possible presence of internal damage with vibration-based methodologies.

We used the pretension system to subject the specimen to eight compression states in its longitudinal direction (forces of 400kN at the highest, and 280kN at the lowest). At every compression state, we struck the structure in the vertical direction three times on the north and south sides of the bridge with an impulse drop weight. Throughout the whole experiment, we recorded ambient seismic noise at different frequency bands with a 14-six-component sensor array to measure the acceleration in the conventional translational components and the angular velocity (rotation rate), a 14-geophone array of 4.5 Hz of natural frequency, and four pairs of embedded ultrasound transducers were used to estimate relative velocity changes (dv/v) by applying the Coda Wave Interferometry (CWI) stretching technique. internal temperature of the concrete was also recorded to correct our measurements by first-order thermal effects.

At the material scale (ultrasound regime) we observe stress-dependent dv/v at four different locations in the specimen and describe them by using the acoustoelastic effect concept regarded as a classical nonlinear phenomenon. We also analyze the relative velocity drop and the subsequent healing process in the concrete triggered by the action of the drop weight. We used the model of Snieder and Sens-Schönfelder (2017) to numerically describe the relaxation process happening at different time scales in the specimen through a deterministic inversion procedure. The north side of the structure showed to have a higher acoustoelastic effect and higher velocity drops, as well as longer relaxation times, it is important to mention that there is evidence of external cracking in this span of the bridge.

We present preliminary results in the seismic frequency band (structural scale), where we expect to observe the influence of the vertical beams that support the bridge on the spatial distribution of changes in dv/v. Changes in the fundamental frequency of the structure as a function of the stress level are also expected.

With this work, we point towards the development of new nondestructive testing methodologies highly sensitive to small cracks and imperfections using conventional and non-conventional seismic instruments, and linear and nonlinear wave propagation models.

How to cite: Dominguez-Bureos, M., Hadziioannou, C., Niederleithinger, E., and Sens-Schönfelder, C.: Time- and stress-dependent elastic properties in a concrete structure; spotting internal damage footprints, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16652, https://doi.org/10.5194/egusphere-egu24-16652, 2024.

Deep-seated, high-alpine rock slides frequently occur in highly schistose, fractured, anisotropic rock masses. Many studies have shown that pre-existing geological structures are decisive for a rock slide’s initiation and kinematics, as they provide weakness zones that may be reactivated in the rock slide process. Besides this structural pre-disposition, internal deformation processes by brittle rock mass fracturing play an important role in the evolution of a rock slide. Nonetheless, the effect of multiscale rock mass fracturing due to the rock slide process is yet to be fully understood. Especially, as it is challenging to measure, characterize, and to numerically model these processes. In our contribution, we present three deep-seated rock slides located in the European Alps in heavily foliated, fractured rock masses with failure volumes above 500.000 m³ each. Focusing on these case studies, we investigate the internal deformation processes with a combined approach, comprising field mapping, laboratory testing, remote sensing, and numerical modelling.

During extensive geological field surveys, we mapped the geomorphological rock slide features and characterized the structural framework of each study site, yielding geometrical models of the rock slides. This provided the basis for our 2D distinct element modelling studies using UDEC, backed by lithological and rock mechanical laboratory investigations.

Whilst UDEC allows for modelling large displacement of blocks bounded by pre-existing discontinuities, it lacks the capability to simulate fracture initiation and propagation of new failure paths within intact blocks, thus neglecting brittle rock mass fracturing. We circumvent this constraint by tessellating the intact rock mass into random polygons – referred to as Voronoi elements. Here, we adapted the original Voronoi technique by assigning an asymmetry to the Voronoi elements, characterized by an elongated axis to consider rock mass anisotropy related to schistosity. By applying this approach, we modelled the fractured, anisotropic, metamorphic rock masses as a combination of pre-existing, field-related structures within a matrix of small, asymmetric Voronoi elements.

In order to confirm the model outputs, we used terrain and deformation data derived from various remote sensing techniques – e.g. satellite based synthetic aperture radar, terrestrial laser-scanning (Riegl VZ 4000) and several campaigns of unmanned aerial vehicle photogrammetry.

In our study, we were able to reproduce the failure mechanism and kinematics of all three rock slides in accordance with our remote sensing deformation data. Thereby, the asymmetric Voronoi tessellation proved to be feasible in reproducing the brittle rock mass fracturing processes in remarkable agreement with our observations in the field. Thus, our results show, how the formation and kinematics of deep-seated rock slides are controlled by the reactivation of pre-existing geological structures and brittle rock mass fracturing. In doing so, our integrated field, laboratory, and numerical modelling approach further contributes to a better understanding of rock slide initiation and kinematics in complex geological media.

How to cite: Gerstner, R., Frießenbichler, M., Avian, M., Maschler, A., Fey, C., Valentin, G., Keuschnig, M., Mair, V., Goldschmidt, F., and Zangerl, C.: High-alpine rock slides controlled by pre-existing geological structures and brittle rock mass fracturing, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-1670, https://doi.org/10.5194/egusphere-egu24-1670, 2024.

When excavating a tunnel, the stresses are distributed asymmetrically along the tunnel cross-section. Other factors, particularly slope friction force and excavation speed, can also contribute to the deformation and displacement of a tunnel. Despite this, several authors have used the complex potential method to predict the ground deformation surrounding the tunnel. However, their applicability to the ground response caused by the asymmetric stress distribution around the mine wall is analyzed in this context. This project, therefore, proposes an approximate solution on the slope to predict the mine cross-section deformation. The solution is based on the complex potential method to predict analytically and numerically the ground deformation around the tunnel. However, two variables called the “complex potential functions” for the Laurent series expansion are used for the stress redistribution to the tunnel boundary conditions. Data from the Datong mine case are used to justify the proposed analytical solutions. The solution is an essential guide for analyzing deformations in complex geological conditions and structures, such as steeper slopes.

How to cite: Mabe Fogang, P. and Huo, B.: Contribution of programming language to novel mine risk assessment project, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3021, https://doi.org/10.5194/egusphere-egu24-3021, 2024.