EMRP1.6 | Fault deformation across scales: from laboratory observations to numerical simulations
Fault deformation across scales: from laboratory observations to numerical simulations
Co-organized by SM4/TS5
Convener: Nathalie Casas | Co-conveners: Chiara Cornelio, Pierre Romanet, Federica Paglialunga, Carolina Giorgetti
Orals
| Wed, 17 Apr, 16:15–18:00 (CEST)
 
Room -2.20
Posters on site
| Attendance Tue, 16 Apr, 16:15–18:00 (CEST) | Display Tue, 16 Apr, 14:00–18:00
 
Hall X2
Orals |
Wed, 16:15
Tue, 16:15
The upscaling of laboratory results to regional geophysical observations is a fundamental question and a current challenge in geosciences. Indeed, earthquakes are non-linear and multi-scale problems, whose dynamics depend strongly on the geometry and the physical properties of the fault and its surrounding medium. To reproduce realistic boundary conditions in the laboratory, fault mechanisms are often scaled down to examine the physical and mechanical characteristics of earthquakes. Small-scale experiments are a powerful tool to study friction and bring to light new insights into weakening or dynamic rupture processes. However, it is not evident how the observed mechanisms can be extrapolated to large-scale observations, and this is where numerical simulations can help to bridge the gap in scale. Laboratory experiments, numerical simulations, and geophysical observations are complementary and necessary to understand fault mechanisms across the different scales. In this session, we aim to convene contributions dealing with multiple aspects of earthquake mechanics, such as:
(i) the thermo-hydro-mechanical processes associated with all the different stages of the seismic cycle, e.g., healing, nucleation, co-seismic fault weakening;
(ii) multidisciplinary studies combining laboratory and numerical experimental results;
(iii) bridging the gap between the different scales of fault deformation mechanisms.

We particularly welcome novel observations and/or innovative approaches to study earthquake faulting. Contributions from early career scientists are highly solicited.

Orals: Wed, 17 Apr | Room -2.20

Chairpersons: Nathalie Casas, Chiara Cornelio, Pierre Romanet
16:15–16:20
16:20–16:30
|
EGU24-15273
|
ECS
|
On-site presentation
Mathias Lebihain, Thibault Roch, Marie Violay, and Jean-François Molinari

Earthquake nucleation is traditionally described using cascading or slow pre-slip models. In the latter, nucleation occurs as the sudden transition from quasi-static slip growth to dynamic rupture propagation. This typically occurs when a region of the fault of critical size Lc, often called nucleation length, is sliding. This transition is relatively well-understood in the context of homogeneous faults. Yet, faults exhibit multiple scales of heterogeneities that may emerge from local changes in lithologies or from its self-affine roughness. How these multiscale heterogeneities impact the overall fault stability is still an open question.

Combining the nucleation theory of [Uenishi and Rice, JGR, 2003] and concepts borrowed from statistical physics, we propose a theoretical framework to predict the influence of brittle/ductile asperities on the nucleation length Lc for simple linear slip-dependent friction laws. Model predictions are benchmarked on two-dimensional dynamic simulations of rupture nucleation along planar heterogeneous faults. Our results show that the interplay between frictional properties and the asperity size gives birth to three (in)stability regimes: (i) a local regime, where fault stability is controlled by the local frictional properties, (ii) an extremal regime, where it is governed by the most brittle asperities, and (iii) a homogenized regime, in which the fault behaves at the macroscale as if it was homogeneous and the influence of small-scale asperities can be averaged.  

Using this model, we explore the overall stability of rough faults, featuring multiscale distributions of frictional properties. We also investigate the stability of velocity-neutral faults that features brittle asperities. Overall, our model provides a theoretical basis to discriminate which heterogeneity scales should be explicitly described in a comprehensive modelling of earthquake nucleation, and which scales can be averaged.

How to cite: Lebihain, M., Roch, T., Violay, M., and Molinari, J.-F.: Impact of multiscale heterogeneities on the nucleation of earthquakes, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15273, https://doi.org/10.5194/egusphere-egu24-15273, 2024.

16:30–16:40
|
EGU24-19397
|
On-site presentation
Fabian Barras, Einat Aharonov, and François Renard

Observations suggest that large earthquakes often propagate as self-healing slip pulses but the mechanical reason of this ubiquity remains debated. Pulse-like ruptures differ from the classical crack-like dynamics by the fact that the slipping portion of the fault is limited to the immediate vicinity of the propagating tip. In this work, we first propose a minimal model describing the dynamics of large earthquakes. In its simplest form, the model contains only two free parameters: a dimensionless stress parameter characterizing the initial state of stress along the fault and a ratio of elastic moduli. The model illuminates how self-healing slip pulses can be produced by the paucity of elastic strain energy that arises once the rupture dynamics interplays with the finite geometry of fault zones—even in the absence of additional mechanisms such as rate-dependent friction.

Next, we discuss the example of faults surrounded by a damage zone whose reduction in elastic wave velocity restricts the flow of strain energy to the rupture tip and promotes pulse-like rupture. Using the proposed model, we demonstrate how the contrast in wave velocities and the initial stress level in the fault zone mediate the propagation mode of the earthquake.

How to cite: Barras, F., Aharonov, E., and Renard, F.: The pulse-like dynamics of large earthquakes illuminated by a minimal elastodynamic model, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-19397, https://doi.org/10.5194/egusphere-egu24-19397, 2024.

16:40–16:50
|
EGU24-7508
|
On-site presentation
Guilhem Mollon, Adriane Clerc, Amandine Ferrieux, Lionel Lafarge, and Aurelien Saulot

Seismic faults are often represented using two different and self-excluding conceptual models. In the first representation, seismic faults are seen as the interface between two surfaces of bare rock, with a roughness extending at all scales. These surfaces interact mechanically through a certain number of “asperities” which constitute the “real contact area”. When adopting this view, attention is paid on the statistics of the asperities population in the fault plane. Faults are thus considered as 2D objects, since their thickness is disregarded.

In the second representation, seismic faults are seen as mathematical planes separated by a certain thickness of granular gouge created by abrasive wear of the surfaces during previous slips. This view is analogous to the tribological “third body” theory, and is supported by field observations and experimental evidences of gouge creation in rotary shear and triaxial experiments. It is convenient to adopt this perspective when weakening phenomena within the gouge are to be spatially resolved in the direction orthogonal to the fault plane. Variations along this plane are then ignored, as well as fault roughness, and faults are mostly seen as 1D objects.

Unification of these two representations requires a better understanding of the interactions between geometrical asperities and a layer of gouge, and in particular of the phenomena that lead to the creation of the latter through the wear of the former. In this communication, we present a numerical model which aims at reproducing lab tests of millimetric single-asperity friction and wear. The model is essentially granular in order to represent the progressive degradation of the asperity along sliding, the separation of powdery matter, its successive ejection and reinjection by the contact (thanks to a periodicity in boundary conditions), and the build-up of a gouge layer. It also includes a coupling with continuum mechanics in order to maintain a meaningful stress field in the asperity beyond the region of degradable rock.

Numerical results show that: (i) the rate of wear of the asperity and the counterface are directly linked to the normal load applied to the contact; (ii) an established layer of gouge develops in the interface and controls the friction coefficient; (iii) a constant level of surface roughness is established after a sufficient sliding distance, both for the asperity and the counterface; (iv) an accurate control of the asperity boundary conditions is necessary in order to obtain repeatable friction and wear. These results are a first step towards a better understanding of the wear kinetics as a function of asperity geometry, load, and roughness, before the introduction of thermal aspects (including melting) in a future version of the model.

How to cite: Mollon, G., Clerc, A., Ferrieux, A., Lafarge, L., and Saulot, A.: A granular numerical model for the friction and wear of a lab-scale fault asperity, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7508, https://doi.org/10.5194/egusphere-egu24-7508, 2024.

16:50–17:00
|
EGU24-3912
|
ECS
|
On-site presentation
Xiaotian Ding, Shiqing Xu, Eiichi Fukuyama, and Futoshi Yamashita

In recent years, an intriguing feature of back-propagating rupture (BPR) has been reported during some earthquakes (Ide et al., 2011; Houston et al., 2011; Hicks et al., 2020; Okuwaki et al., 2021; Vallée et al., 2023). The occurrence of BPR challenges the classical interpretation of rupture propagation as a “forward” problem, while remaining less understood by the earthquake science community. Here, using fracture mechanics, we first argue that BPR is nothing but an intrinsic component of rupture propagation; however, its observability is usually masked by the superposition effect of interfering waves behind the primary, forward-propagating rupture front. We then suggest that perturbation to an otherwise smooth rupture process can break the superposition effect and hence can make BPR observable. To test our idea, we report results of mode-II rupture propagation from both numerical simulations and laboratory observations. By introducing a variety of perturbations (lateral variation in bulk or interfacial properties along the fault, singular stress drop, coalescence of two rupture fronts, etc.) during rupture propagation, we show that prominent phases of BPR indeed can be successfully excited. We further classify BPR into two modes: higher-order rupture or interface wave, depending on whether the already-ruptured fault is quickly healed and whether additional stress drop can be produced. Lastly, we propose several application potentials for BPR, such as constraining the velocity structure of fault zones, probing the mechanical state of faults, and studying the stability of perturbed slip along a homogeneous or bimaterial interface. Our study refines the understanding of the nature and complexity of rupture process, and can help improve the assessment of earthquake hazards.

How to cite: Ding, X., Xu, S., Fukuyama, E., and Yamashita, F.: Back-Propagating Rupture: Nature, Excitation, and Applications, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3912, https://doi.org/10.5194/egusphere-egu24-3912, 2024.

17:00–17:10
|
EGU24-10998
|
ECS
|
On-site presentation
Dong Liu and Nicolas Brantut

During natural and induced seismic activities, pore fluid pressure within fault zones and their surrounding rock may respond differently to stress variations, introducing additional complexities to seismic hazard assessment. While theoretical investigations have recognized the influence of such poroelastic heterogeneity on fault instability, incorporating phenomena like slip-induced dilation or compaction, the chosen poroelastic properties in these studies lack robust constraints from experimental measurements. Addressing this gap, our study focuses on quantifying the heterogeneity of poroelastic properties in the presence of a fresh fault, aiming to elucidate the coupling between poroelasticy and fault dilatancy during fault slip.

In our experimental investigation, we examined the evolving dynamics of pore pressure both on- and off-fault in initially intact Westerly granite samples. Applying confining stress of 100 MPa and a pore pressure of 60 MPa at two sample ends to replicate crustal settings, we induced a sliding fault plane through loading to failure under a constant strain rate. In the faulted samples, we measured the pore pressure response under sudden step loading in the direction of the maximum compression σ1. Each loading step of around 5 MPa was imposed incrementally increasing the differential stress from 5 MPa to approximately 80 MPa (frictional resistance) after achieving pore pressure equilibrium. Detailed measurements, including displacement, bulk deformation, differential stress, local pore pressure and acoustic emissions were recorded throughout these tests. A spring-slider model coupled with 1-D fluid diffusion was used to try to simulate experimental observations.

Our results indicate that both the shear zone and the bulk exhibit a diminishing Δp/Δσ1 with increasing differential stress. Measurements within the fault zone consistently yield positive values, surpassing those off the fault, with the discrepancy more pronounced at lower stress levels. In regions farther away from the shear zone, the off-fault response Δp/Δσ1 presents a smaller value compared to locations proximal to the fault zone and may even exhibit slight negativity. During fault slip, on-fault measurements exhibit an instantaneous increase upon step loading followed by a gradual decrease, as a result of the interplay between poroelasticity and fault dilatancy. These observations were effectively reproduced by the numerical model integrating the poroelastic measurements and rate-and-state fault friction with slip-dependent dilatancy. The implications of this investigation extend to an enriched understanding of the heterogeneity in poroelastic responses between fault zones and host rocks, serving as valuable benchmarks for informing future numerical simulations, particularly in the context of naturally formed fresh faults. 

How to cite: Liu, D. and Brantut, N.: Poroelastic heterogeneity in the presence of a fresh fault: experimental insights and numerical modelling, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10998, https://doi.org/10.5194/egusphere-egu24-10998, 2024.

17:10–17:20
|
EGU24-11200
|
On-site presentation
Francois Passelegue, Pierre Dublanchet, Nicolas Brantut, and Hervé Chauris

A growing amount of evidence indicate that aseismic transients driven by overpressure play an important role in the triggering of induced seismicity. Understanding the physical control on aseismic slip development is thus important for seismic hazard assessment. We conducted an investigation into the propagation dynamics of a fluid-driven slip front along a laboratory frictional interface composed of granite. The experiments were carried out under a confining pressure of 90 MPa, with an initial uniform fluid pressure of 10 MPa. Fault reactivation was initiated by injecting fluids through a borehole directly connected to the fault.

Our findings reveal that the peak fluid pressure at the borehole leading to reactivation exhibits an increase proportionate to the injection rate. Employing three fluid pressure sensors and eight strain gauges strategically positioned around the experimental faults, we performed an inversion analysis to image the spatial and temporal evolution of (i) hydraulic diffusivity and (ii) kinematic fault slip during each injection experiment. Our inversion methods integrated both deterministic and Bayesian procedures, facilitating the tracking of the fluid pressure front along the fault interface and the subsequent propagation of the slip front over time.

The migration pattern shares many similarities with natural slow slip events suspected to play a role in the development of natural and induced earthquake swarms or aftershock sequences.  We demonstrate that increasing the fluid injection rate induces a transition from a quasistatic propagation of the slip front correlated with the increase in fluid pressure to a dynamic scenario where the slip front outgrows the fluid pressure front, accelerating during its propagation. Furthermore, we establish that temporarily shutting off fluid pressure during injection induces the propagation of a pore-pressure back-front, which halts the propagation of the slip front, aligning with theoretical expectations.

How to cite: Passelegue, F., Dublanchet, P., Brantut, N., and Chauris, H.: Influence of Injection rate and slip-induced dilatancy on the propagation of fluid-driven slip front, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11200, https://doi.org/10.5194/egusphere-egu24-11200, 2024.

17:20–17:30
|
EGU24-13433
|
On-site presentation
Michele Fondriest, Fabio Arzilli, Benoit Cordonnier, Michael Carroll, and Mai-Linh Doan

The propagation of earthquake fault ruptures in the crust involve the generation of unloading stress pulses sufficiently large to induce dynamic failure of water-saturated rocks under tensional stresses and hydrofracturing. Similar processes are also activated during underground rock mass excavation activities in mines and tunnels. The current knowledge about rock fracturing via dynamic unloading is mainly limited to empirical records and numerical simulations, while there is a general paucity of experimental studies, due to difficulties in reproducing large instantaneous decompressions on rock samples using standard triaxial rigs. Until now rapid decompression and fracturing of large rock samples in dry conditions was reported only by using an unconventional gas-confined vessel.

Here, we report rock-fracture results for newly conceived rock decompression experiments, completed through the innovative use of a “cold-seal pressure vessel” (CSPV) apparatus which is routinely employed in experimental petrology. We applied instantaneous large decompressions on water-saturated rock samples equilibrated at high confinement (up to 200 MPa) and temperatures (up to 540°C). The tested rock samples were fine-grained Westerly granite, coarse-grained tonalite and micritic limestone. During the decompressions the rock samples hydrofractured due to the confinement dropping faster than the pore pressure within the rock. Porosity measurements, SEM imaging and X-ray µCT acquired before and after the tests suggest that the magnitude of dynamic fracturing not only positively correlates with the pressure drops but it mostly increases when the decompression is associated to a phase change of the pore water (e.g. supercritical fluid to subcritical gas) . Water vaporization or degassing imply an instantaneous volume expansion (up to 70 times) which critically enhances dynamic fracture propagation along rock grain boundaries. The induced fractures span from mm-long transgranular cracks to microcracks with submicrometric aperture. Therefore, synchrotron light high-resolution microtomography (final pixel resolution of 0.3 µm) was employed to fully resolve and quantify the 3D fracture networks of these deformed rock samples. Such unique dataset allowed us to determine at different scales the fracture intensity, aperture and connectivity of the dynamically induced fracture networks and to assess the key contribution of pore-water physical state changes on the initial stages of dynamic fracturing in rocks at crustal conditions. Such results will contribute to close a current knowledge gap in rock mechanics.

How to cite: Fondriest, M., Arzilli, F., Cordonnier, B., Carroll, M., and Doan, M.-L.: The effect of pressure drop and fluid expansion during rock fracturing by dynamic unloading, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13433, https://doi.org/10.5194/egusphere-egu24-13433, 2024.

17:30–17:40
|
EGU24-2384
|
ECS
|
On-site presentation
Songlin Shi, Meng Wang, Yonatan Poles, and Jay Fineberg

Earthquake-like ruptures disrupt the frictional interface between contacting bodies and initiate frictional motion (stick-slip). The interfacial slip (motion) immediately resulting from a rupture during each stick-slip event is usually much smaller than the total slip recorded during the duration of the event. Slip after the onset of friction is generally attributed to the continuous motion of global ‘dynamic friction’. Here, we demonstrate that numerous hitherto invisible secondary ruptures are initiated immediately after each initial rupture by directly measuring the contact area and slip at the frictional interface. Each secondary rupture generates incremental slip that, when not resolved, may appear as steady sliding of the interface. Each slip increment is linked, via fracture mechanics, to corresponding variations of contact area and local strain. Cumulative interfacial slip can only be described if the effects of these secondary ruptures are taken into account. These weaker slip sequences can also be observed in bimaterial interfaces and exhibit strong directional effects. In addition, the seismic moments we estimate based on slip sequences are consistent with the Gutenberg-Richter (G-R) law. These results have important implications for our fundamental understanding of frictional motion and the important role of aftershocks within natural faults in generating earthquake-mediated slip/afterslip.

How to cite: Shi, S., Wang, M., Poles, Y., and Fineberg, J.: Frictional slip sequences in homogeneous and bimaterial interfaces, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2384, https://doi.org/10.5194/egusphere-egu24-2384, 2024.

17:40–17:50
|
EGU24-5187
|
On-site presentation
Fabio Corbi, Giacomo Mastella, Elisa Tinti, Adriano Gualandi, Laura Sandri, Matthias Rosenau, Silvio Pardo, and Francesca Funiciello

Modeling the seismic cycle requires multiple assumptions and parameters. Providing a quantitative assessment of the model behavior is pivotal for determining the degree of similarity between different scales and modeling strategies and for exploring dependencies with respect to selected parameters. Here we compare stick-slip ruptures nucleating spontaneously in scaled seismotectonic models (i.e., laboratory experiments capturing the first-order physics of the seismic cycle of subduction megathrusts) with slow earthquakes in nature. We rely on two non-dimensional parameters, namely the Ruina number (Ru) and system dimension (D) to quantify model behavior. Ru is proportional to the ratio of the asperity size to the critical nucleation size. Within the rate- and state friction framework, for velocity weakening asperities Ru controls the behavior of the system, which can be either periodic or not, and it can exhibit both slow and fast ruptures. D measures how complicated the system evolution is. D reveals how many variables are required to describe the seismic cycle because it tells us the minimum dimension needed to embed the observed dynamics. 

By coupling the Simulated Annealing algorithm and quasi-dynamic numerical simulations, we retrieve rate and state friction parameters characterizing single asperity models with different lateral extent of the velocity weakening patch. Similarly to slow earthquakes, we found optimal rate and state parameters indicative of low (< 4) Ru. We also found a direct proportionality between Ru and the lateral extent of the asperity. 

Next, we implement tools from non-linear time-series analysis and Extreme Value Theory to compute D from models of different sizes, materials, deformation rates and frictional configurations (single or twin asperities along strike). Our analysis supports the existence of a low dimensional attractor (D<5) describing the dynamics of scaled seismotectonic models. In particular, our models display D=3.0-4.2, which is remarkably similar to D=3.2 of slow earthquakes identified along the Cascadia subduction zone. Under the explored conditions, D appears more affected by the material behavior of the analog upper plate (i.e., gelatin vs. foam rubber) than by the lateral frictional segmentation of the megathrust.

Despite the different spatio-temporal scales, our results support a scenario where scaled seismotectonic models and slow earthquakes share similar dynamics.



How to cite: Corbi, F., Mastella, G., Tinti, E., Gualandi, A., Sandri, L., Rosenau, M., Pardo, S., and Funiciello, F.: The similarity between ruptures in scaled laboratory seismotectonic models and slow earthquakes , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5187, https://doi.org/10.5194/egusphere-egu24-5187, 2024.

17:50–18:00
|
EGU24-11536
|
Highlight
|
On-site presentation
Men-Andrin Meier, Domenico Giardini, Stefan Wiemer, Massimo Cocco, Florian Amann, Elena Spagnuolo, Paul Selvadurai, Elisa Tinti, Luca Dal Zilio, Alba Zappone, Giacomo Pozzi, Mohammadreza Jalali, and Valentin Gischig and the FEAR science team

Our understanding of earthquake rupture processes is generally limited by the resolution of available observations. In all but exceptional cases, earthquake observations are made at comparatively large distances from the rupture itself, which puts a limit on what spatial scales can be resolved. At the same time, it is clear that small scale processes may play a crucial, if not dominant, role for various seismogenic processes, including rupture nucleation, co-seismic weakening and stress re-distribution.

The Fault Activation and Earthquake Rupture ('FEAR') project aims at collecting and interpreting a multitude of earthquake-relevant observations from directly on and around the process zone of an induced earthquake. To this end, we attempt to activate a natural granitic fault zone in the BedrettoLab, at a depth of ~1km, after instrumenting the fault zone with a multi-domain and multi-scale monitoring system. The goal is to observe and study earthquake rupture phenomena in a natural setting, from unusually close distance.

In this talk, we outline the project status, the science goals, and the plans for the main experiments, which are scheduled for the years 2024 - 2026. Notable milestones we report on include

  • the identification and detailed characterisation of the target fault zone
  • the beginning of niche and tunnel excavations
  • laboratory experiments that characterise the frictional and mechanical behaviour of both gauge material and host rock of the target fault zone
  • development of numerical models for 2D and 3D dynamic rupture propagation
  • development of tailored monitoring methods for seismicity, strain, temperature, pressure, bio-geo-chemistry and other relevant observables
  • development of remote experiment control methods
  • test stimulations in a nearby rock volume of similar geology, with an already existing monitoring system, where we tested the influence of pre-conditioning injection protocols
  • similar test stimulations in the same volume where we aim at triggering a larger event (target Mw~0)
  • active seismic experiments in an underground salt mine, to calibrate the very- to ultra-high frequency (1k Hz - 500k Hz) acoustic emission sensors

Together, these and other efforts constitute the necessary ingredients we need for interpreting the near-source observations that we will collect during the fault activation experiments.

How to cite: Meier, M.-A., Giardini, D., Wiemer, S., Cocco, M., Amann, F., Spagnuolo, E., Selvadurai, P., Tinti, E., Dal Zilio, L., Zappone, A., Pozzi, G., Jalali, M., and Gischig, V. and the FEAR science team: Fault activation from up close, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11536, https://doi.org/10.5194/egusphere-egu24-11536, 2024.

Posters on site: Tue, 16 Apr, 16:15–18:00 | Hall X2

Display time: Tue, 16 Apr, 14:00–Tue, 16 Apr, 18:00
Chairpersons: Pierre Romanet, Carolina Giorgetti, Nathalie Casas
Laboratory experiments: thermo-hydro-mechanical processes associated with all the different stages of the seismic cycle, e.g., healing, nucleation, co-seismic fault weakening;
X2.144
|
EGU24-18100
|
ECS
Nico Bigaroni, Julian Mecklenburgh, and Ernest Rutter

During interseismic periods a fault at depth can experience non-constant effective normal stress due to fluctuations in the pore-fluid pressure. Pore-pressure oscillations may influence the healing capability of the fault and ultimately affect its reactivation. Thus, studying the behaviour of faults during interseismic periods is a critical factor in understanding the seismicity. Triaxial tests were conducted using saw-cut (45o) samples of Pennant Sandstone to investigate the influence of pore-pressure oscillations during slide-hold-slide (SHS) tests (th = 900 – 7300s) on its frictional behaviour and fault reactivation. The cylindrical samples were hydrostatically compacted at 30 MPa and pore-pressurized with argon gas at 5, 10 and 18 MPa resulting in effective normal stress (σ’n) 25, 20 and 12 MPa, respectively. Then the saples were deformed at a constant shear displacement rate ≈ 4.5 μm/s. To overcome the displacement hardening tendency of the sample geometry, we servo-controlled the confining pressure so that the resolved normal stress on the sliding surface is kept constant. Experimental observations revealed a significant influence of pore-pressure oscillation on the frictional behaviour resulting in an increase in both frictional healing and creep relaxation. Moreover, this effect was enhanced as the effective normal stress was increased further. To understand better the underling mechanism(s) that influences these time-dependent processes we coupled the frictional results with permeability measured using the oscillating pore pressure method during the SHS tests. Finally, we tested how the pore-pressure oscillation affected the fault reactivation by conducting creep experiments at constant shear stress while the fault was brought to reactivation via progressive increase in fluid pressure. Our results demonstrated how non-constant effective normal stress history during interseismic periods deeply affects the fault behaviour, with important implications for natural and human-induced seismicity.

How to cite: Bigaroni, N., Mecklenburgh, J., and Rutter, E.: Frictional Response of Clay-rich Sandstone to Pore-Pressure Oscillation Throughout Interseismic Periods, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18100, https://doi.org/10.5194/egusphere-egu24-18100, 2024.

X2.145
|
EGU24-9541
Marco Scuderi, Nathalie casas, Giuseppe Volpe, and Cristiano Collettini

Mineralogy, fabric, and frictional properties are fundamental aspects of natural and experimental faults that concur in controlling the fault strength and the fault slip behavior. Mineralogy controls the fabric evolution influencing the micro-mechanisms at play during fault deformation and needs an in-depth investigation to better understand and foresee the frictional response of experimental faults. Classically, this investigation has been conducted by relating the fault frictional behavior to the post-experimental microstructures. However, this “classical” approach provides a direct but static view of the fault deformation where the evolution of fabric with deformation can be only speculated.

To investigate in “real-time” the deformation micro-mechanisms at play during the experiments, the recording and analyses of Acoustic Emissions (AEs) produced by the deforming fault gouge can provide new insights.

In this study, we present a systematic study of microstructural, mineralogical, frictional, and AEs analysis coming from a suite of frictional experiments in a double direct shear configuration (biaxial apparatus, BRAVA2). We conducted experiments on gouges made of bi-disperse and layered mixtures of quartz and phyllosilicate. These experiments were performed at a constant normal stress of 52MPa and under 100% humidity. The friction evolves with the phyllosilicate content from µ ~ 0.6 for 100% quartz to µ ~ 0.4 for 100% phyllosilicates. At the end of the experiments samples were carefully collected and prepared for microstructural analysis. The fabric of the experimental samples show an evolution from localized to distributed and foliated fabric with increasing amount of phyllosilicate content.

We then integrate specific features of AEs, such as amplitude and AE rate, to unveil the micro-mechanisms at play during the experimental fault deformation. Our results show that the overall AE behavior is controlled by mineralogy. Deformation of quartz gouge produces the largest number of AEs whereas phyllosilicates are almost not producing AEs. Furthermore, the AE behavior of bi-disperse mixtures of quartz and phyllosilicates is strongly controlled by the amount of phyllosilicates. In fact, increasing the amount of phyllosilicate, the number, the rate, and the amplitude of AEs decrease. This behavior could be explained by the lubricant role of phyllosilicates which hinder the interaction between quartz grains favoring foliation sliding as main deformation mechanism and thus reducing the frictional strength. These results suggest that for bi-disperse mixtures the AEs reflect the frictional behavior of the mixture. Layered quartz-phyllosilicates mixtures show instead a non-trivial acoustic emission behavior which cannot be directly related to the measured frictional strength of the layered mixture: friction is controlled by the frictionally weaker mineral phase, whereas the AEs are probably dependent by the interplay between the stronger and weaker phase of the layered mixture.

Our results show that fault fabric together with mineralogy strongly control the micro-mechanisms at play during deformation and therefore the frictional response. Our findings support the use of the AE analysis as a new tool for the investigation of the micro-mechanisms at play during deformation, improving our interpretation of the mechanical behavior of fault gouges.

How to cite: Scuderi, M., casas, N., Volpe, G., and Collettini, C.: Unraveling the micro-mechanics of shear deformation through acoustic attributes of quartz-muscovite mixtures, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9541, https://doi.org/10.5194/egusphere-egu24-9541, 2024.

X2.146
|
EGU24-15205
|
ECS
Federica Paglialunga, Francois Passelegue, and Marie Violay

Many aspects of earthquake physics are still not completely understood given its intrinsically complex nature. Among the others, the nucleation process; when and where an earthquake will occur, as well as its magnitude. Seismology is a commonly used method for studying earthquakes, but it faces challenges in accessing precise information about the physical processes taking place on the fault plane.

Here, we show how laboratory seismology can directly shed light on fault plane dynamics. Our approach involves reproducing in the laboratory on a large biaxial apparatus with a fault length of 2.5 m generated by two analog (PMMA) samples brought into contact. The experimental setup allows to impose both a heterogeneous loading distribution through the use of independent pistons loading the fault in the normal direction and specific boundary conditions (i.e. by modifying stopper and puncher dimensions). The stress state is measured through strain gauges at high frequency (40 KHz) along 15 locations along the fault. The experiments provide insights into two crucial aspects of laboratory earthquakes: (i) the nucleation location of ruptures and (ii) the complexity of the seismic cycle.

Our findings reveal that the initial on-fault stress distribution plays a significant role in both aspects. We observe that ruptures consistently nucleate in locations where the stress ratio τ/σn is the highest. Notably, such values change among experiments, challenging the widespread notion that a friction coefficient solely governs the onset of instability. Furthermore, we demonstrate how the heterogeneity of the initial prestress distribution along the fault controls the complexity of the seismic cycle. In certain cases, the seismic cycle manifests as system-size events with complete ruptures occurring regularly in time, devoid of precursors. Conversely, other initial stress distributions generate more complex cycles, characterized by multiple precursors before a main rupture, predominantly occurring in zones of elevated τ/σn (referred to as 'friction asperity'). The complexity of the seismic cycle can be described in terms of the number of precursory events, inter-event time, and the size of finite ruptures.

This study, carried out in a long laboratory fault, highlights the complexities that emerge when heterogeneous, hence more realistic, stress conditions are applied, providing valuable insights into the physics of natural earthquakes.

How to cite: Paglialunga, F., Passelegue, F., and Violay, M.: Initial stress distribution dictates nucleation location and complexity of the seismic cycle of long laboratory faults, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15205, https://doi.org/10.5194/egusphere-egu24-15205, 2024.

X2.147
|
EGU24-11407
|
ECS
Sofia Michail, Paul Antony Selvadurai, Markus Rast, Antonio Felipe Salazar Vásquez, Patrick Bianchi, Claudio Madonna, and Stefan Wiemer

Faults in nature exhibit complex surface characteristics with patches of the fault (asperities) that may slip dynamically while other sections are more prone to creep (Beeler et al., 2011). Asperities forming in nature may be due to the geometric interactions between surfaces within a fault that contribute to complex stress states that are not well understood. Fault roughness is believed to play an important role in the control of the contact conditions established by asperities, directly affecting its potential to slip unstably. How the asperities are formed and how their seismogenic properties evolve due to wear is an important question with implications to slip budget and earthquake potential.

In this study, we performed a triaxial experiment at sequentially increasing confining pressures (Pc = 60, 80, 100 MPa) on a saw-cut sample of Carrara marble. We analysed the quasi-static frictional response that benefited from novel arrays of distributed strain sensors (DSS) obtained using fiber optics. This sensor offered unique insight into the axial strain with a spatial resolution of 2 mm. The frictional behaviour during the first confining pressure step exhibited a dynamic instability in the form of a stick-slip event (SS) that produced a measurable stress drop. In the subsequent confining pressure stages, where an increase in confining pressure translated to increased normal stress, the fault behaved in a stable manner and no dynamic instabilities were produced. This observation is inconsistent with frictional stability theory (e.g. Rubin and Ampuero, 2005) and required pre- and post-mortem campaigns into the surface characteristics and their evolution to explain this abnormal behaviour. Therefore, we employed experimental techniques (pressure sensitive film (PSF), optical and stylus profilometry) along with finite element (FE) model in ABAQUS to characterize the pressure and roughness.

The DSS array showed extensional axial strain closer to the edges of the fault, while only compression was expected in this triaxial loading test. The pre-experimental profilometry revealed an asperity located at the centre of the fault with a curvature ratio of h/L=0.1% inherited from the hand-lapping preparation, which dominated the initial contact conditions prior to the SS and explained the DSS observations. The DSS results were confirmed using a FE model which justified the effect of the fault geometry (h/L) on the strain response. After the SS, wear and smoothening of the central asperity was seen in roughness measurements. The profilometric measurements showed that gouge was deposed adjacent to the high normal stress asperity center (PSF) and were characterized by increased RMS roughness. These small amounts of gouge on the fault surface were sufficient to suppress the seismic response of the asperity. These findings show that the seismic potential of a carbonate (softer) asperity, may be highly influenced by the debris produced during wear. Its impact on earthquake nucleation could provide insight into large-scale earthquake preparation processes on carbonate faults in nature.

 

References:

  • Beeler, M., Lockner, D. L. and Hickman, S. H. (2001), Bull. Seis. Soc. Am., 91 (6): 1797–1804
  • Ampuero, J.-P. and Rubin, A. M. (2008), J. Geophys. Res., 113, B01302

How to cite: Michail, S., Selvadurai, P. A., Rast, M., Salazar Vásquez, A. F., Bianchi, P., Madonna, C., and Wiemer, S.: Laboratory Insight into the Evolution of the Seismic Potential of an Asperity due to Wear, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11407, https://doi.org/10.5194/egusphere-egu24-11407, 2024.

X2.148
|
EGU24-14566
Weifeng Qin, Lu Yao, Tongbin Shao, Wei Feng, Jianye Chen, and Shengli Ma

The frictional properties of faults are primarily controlled by their mineral composition, as well as ambient and deformation conditions, such as temperature, pore fluid, normal stress, and slip displacement. While many studies have been conducted to decipher how temperature and pore fluid may affect the frictional behavior of faults, less attention has been paid to the slip displacement effects, especially under hydrothermal conditions. By employing a rotary shear apparatus equipped with an externally-heated hydrothermal pressure vessel, we conducted large-displacement (up to 521 mm) friction experiments on chlorite under temperature (T) of 25 to 400℃ and pore water pressure (Pp) of 30MPa. The imposed effective normal stresses were 200 MPa and the slip rates ranged from 0.4 to 10 μm/s. The experiments unveiled significant slip strengthening in chlorite within the temperature range of 25 to 400 °C. Moreover, with increasing temperatures, there was an overall increasing trend in both the rate of slip strengthening and the ultimate frictional strength. For example, under T = 25 °C, the friction coefficients at displacements of 5, 90, and 521 mm were 0.33, 0.49, and 0.59, respectively, in contrast to 0.46, 0.79, and 0.88, respectively, at the same three displacements under T =400 °C. Under all the temperature and displacement conditions, chlorite exhibited velocity strengthening behavior without discernible temperature dependence, although the velocity-dependence parameter (a-b) increased with slip displacement. Microstructural analysis revealed that, the entire layer of the chlorite gouge experienced pervasive and intense shear deformation after slip of 521 mm, with extremely remarkable grain-size reduction. The thermogravimetrical and FTIR data of the deformed chlorite samples, together with the microstructural data, suggest that the dehydroxylation and the distortion of crystal structure of chlorite might occur during the friction experiments conducted at T ≥ 200 °C. Such changes may explain the more pronounced slip strengthening of chlorite with increasing temperatures towards 400 °C. This explanation can be further demonstrated by a comparative experiment conducted under varying temperatures (400°C for the first 100 mm of slip, followed by 25°C for the rest of 100 mm slip), wherein the friction coefficient at T = 25°C during the latter stage of slip remains as high as that at T = 400°C. These findings highlight the importance of slip displacement in controlling the frictional strength and its variations of chlorite-bearing faults at depths, and have profound implications for understanding the fault slip behaviors and earthquake mechanisms in subduction zones.

How to cite: Qin, W., Yao, L., Shao, T., Feng, W., Chen, J., and Ma, S.: Frictional behavior of chlorite in large-displacement experiments under hydrothermal conditions, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14566, https://doi.org/10.5194/egusphere-egu24-14566, 2024.

X2.149
|
EGU24-9941
|
ECS
Michele De Solda, Michele Mauro, Federico Pignalberi, and Marco Scuderi

In the last decades, rock mechanics laboratory experiments have allowed framing earthquake physics as a frictional problem. When the accumulated stress on a fault exceeds the frictional forces holding it in place, a rapid acceleration occurs. This movement can be stable or unstable, involving phases of adhesion (stick) and rapid sliding (slip). In these terms, an earthquake results from the release of mechanical energy during one of these slip phases.

 

Modern friction theories propose that the frictional forces holding the fault in place are controlled by small asperities defining the real contact area (RCA). Therefore, understanding the mechanics of contacts on the fault and their evolution under stress and velocity changes can shed light on the microphysical processes underlying earthquakes.

 

In the laboratory, it is now possible to investigate the dynamics of experimental faults, predicting their instability behavior based on Rate and State Friction theory and its experimentally obtainable parameters (a-b, Dc). However, these parameters lack an explicit relationship with contact mechanics, necessitating additional measurements complementing the system's state information. One of the most widely used techniques for studying RCA during laboratory experiments involves investigating changes in acoustic transmissivity (velocity, amplitude) of generated and recorded ultrasonic waveforms (UW) passing through the sample during the deformation. At a given wavelength, analytical expressions for these quantities depend on the elastic properties and densities of the fault portion crossed by the wave. Simultaneous knowledge of stress conditions and elastic properties allows the formulation of constitutive laws for the evolution of contacts between fault asperities.

 

In double direct shear experiments (DDS) within biaxial apparatuses, the sample dimensions (gouge) impose stringent limits on the spatial and temporal resolution of the signal. These limits highlight the current sensor technology's deviation from the ideal behavior.

 

Here, we present a methodology and a waveform recording and synchronization protocol

implemented on the biaxial apparatus BRAVA2 in the Rock Mechanics and Earthquake Physics laboratory at Sapienza University of Rome. We focus on the types of sensors used and their specifications to provide accurate measurements of the deformation processes occurring within the gouge layers.

 

Several studies have conducted DDS experiments using UW, but they rarely take into account the characterization of the impulse signal, various reflections in the sample assembly, and conversion modes of the generated waveforms. These are all essential components to identify the interaction of the experimental system with the propagation of the ultrasonic waves, to exploit the received signal in its entirety.

We believe that a careful signal characterization is necessary to fully understand the physical processes during deformation within the sample and, consequently, to attempt upscaling to natural earthquakes.

How to cite: De Solda, M., Mauro, M., Pignalberi, F., and Scuderi, M.: Strengthen and limitations of ultrasonic wave testing: examples from Double Direct Shear experiments on gouge, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9941, https://doi.org/10.5194/egusphere-egu24-9941, 2024.

X2.150
|
EGU24-8381
Lei Zhang, Changrong He, and Sylvain Barbot

To investigate the frictional behavior of basalt under hydrothermal conditions, we apply sliding experiments using basalt gouge under the temperature of 100-600ºC, effective normal stress of 150MPa, and fluid pressure of 30MPa and 100MPa, respectively. Experiment results under 30MPa pore pressure show that basalt exhibits velocity-strengthening behavior at 100-200ºC and changes to velocity-weakening behavior at 400-600ºC; meanwhile, at 400ºC, velocity dependence of basalt evolves with slip from initial velocity weakening to velocity-strengthening. Results under 100MPa fluid pressure show a similar transition of velocity dependence at 300ºC; however, at higher temperatures of 400-600ºC, velocity strengthening behavior occurs, accompanied by strong slip weakening behavior at the slowest loading rate (0.04μm/s). During the velocity step, the experiment exhibits a rate-dependent creep without transient evolution with slip. Microstructure observation reveals significant differences between samples sheared under 30MPa and 100MPa fluid pressure. At higher fluid pressure and temperatures of 400-600ºC, the porosity of the basalt gouge layer is significantly reduced, and deformation is characterized by pervasive shear with no apparent localization. Such results suggest that the healing process/plastic deformation is activated at higher fluid pressure, leading to slip stability transition and slip-weakening of frictional strength.

How to cite: Zhang, L., He, C., and Barbot, S.: Transition from Unstable Slip to Rate-Dependent Creep Controlled by High Fluid Pressure, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8381, https://doi.org/10.5194/egusphere-egu24-8381, 2024.

X2.151
|
EGU24-8887
|
ECS
Jianhua Huang, Bo Zhang, Junjie Zou, Honglin He, Jiaxiang Dang, and Jinjiang Zhang

        Abundant cherts (nodules and bands) were discontinuously hosted by dolostones of the Mesoproterozoic group Strata (∼1.5 Ga) in the Shanxi graben system, North China, where earthquakes are common. Measurements of the shear strength and stability of granular quartz reveal that quartz is a typical tectosilicate which exhibits high frictional strength and velocity-weakening properties. Conversely, dolomite is usually frictionally weak but velocity strengthening. The two minerals also behave differently during coseismic slip. Due to the high temperature generated by frictional heating, the thermal decomposition of dolomite usually results in calcite, periclase nanoparticles and carbon dioxide. However, quartz melts by friction at high temperatures. In order to investigate the role of quartz in dolomite fault rock during the process of coseismic slip, high velocity shear experiments were conducted on the quartz-bearing dolomite fault gouge taken from Yuguang Basin South Margin Fault (YBSMF), northeast of the Shanxi graben system. Also, we carried out high velocity experiments with the synthetic quartz-dolomite gouge with different mass ratio. For a slip velocity ≥ 0.1 m/s and normal stresses from 1.0 to 1.5 MPa, the friction values of the gouge decrease exponentially from a peak value of more than 0.5 to a steady-state value from 0.1 to 0.4. This high-velocity weakening feature was observed in the synthetic quartz-dolomite gouge as well as in the YBSMF gouge. With the increase of quartz content, the slip weakening distance (Dw) increases from 4.27 to 13.24 m, and the steady-state friction coefficient increases from 0.2 to 0.4 at 1.0 MPa normal stress and 1.0 m/s slip velocity. The textures of the gouge are characterized by grain comminution, R shear planes and localized deformation zone in the friction weakening stage. The slip surfaces are characterized by mirror-like smooth surface and nanoparticle aggregates. The theoretical calculation results show that the temperature inside the gauge layer did not exceed 300 °C during the experiments. However, the microstructures present that the dolomite experienced thermal decomposition, indicating that the temperature at the asperity exceeds 550 ℃. We suggest that thermal decomposition together with flash heating may lead to the slip-weakening behavior of quartz-bearing dolomite gauge, and the addition of quartz will increase of the strength of the dolomite gouge.

How to cite: Huang, J., Zhang, B., Zou, J., He, H., Dang, J., and Zhang, J.: High-velocity frictional behavior and microstructure evolution of quartz-bearing dolomite fault gouge, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8887, https://doi.org/10.5194/egusphere-egu24-8887, 2024.

X2.152
|
EGU24-9165
|
ECS
Wei-Hsin Wu, Wei Feng, Rodrigo Gomila, Telemaco Tesei, Marie Violay, Anette K. Mortensen, and Giulio Di Toro

Fault’s frictional strength and particularly its ability to heal during the interseismic period (fault frictional healing Δμ) is critical to understand the seismic cycle, yet the understanding of temperature and phase-dependent healing characteristics of natural geothermal conditions remains limited. Here we examined the frictional healing of both simulated fresh and chlorite-altered basaltic gouge from Krafla geothermal field (Iceland) under realistic geothermal conditions of water temperature Tf = 100-400 ˚C and pressure Pf = 10-30 MPa (water in liquid, vapor and supercritical state) by performing Slide-Hold-Slide (SHS) experiments. All experiments were performed under a constant effective normal stress of 10 MPa and initiated with a 5-mm run-in slip at a loading point slip rate V of 10 mm/s before the SHS sequence. For each SHS sequence, shearing was held from 3 s to 10,000 s, separated by a slip interval of 1mm. Our mechanical results indicate that frictional healing, the difference between peak friction reached upon re-shear and the steady-state friction before the hold, increases with increasing logarithm of hold time in all experiments, as suggested by previous studies. Meanwhile, frictional healing rate (β = Δμ/log(1+ thold/tcutoff)), commonly regarded as the quantification of the rate of healing, increases with increasing temperature for both fresh and altered basalt. For fresh basalt, β increases from 0.007 at Tf = 100 ˚C to 0.060 at Tf = 300 ˚C (liquid) before dropping to 0.036 at Tf = 400 ˚C (vapor) and eventually increases to 0.096 at Tf = 400 ˚C (supercritical). For altered basalt, β  increases continuously from 0.003-0.007 at Tf = 100 ˚C to 0.013-0.022 at Tf = 300 ˚C and reaches its maximum value of β = 0.024-0.035 at Tf = 400 ˚C (vapor) and β = 0.030 at Tf = 400 ˚C (supercritical). Besides this temperature-dependent relationship, the dramatic decrease of β in fresh basalt to values similar to those of altered basalt when water changed from liquid to vapor state also suggests that the physical state of water can control the healing rate. Subsequent microanalytical analyses (XRPD, XRF, SEM-EDS) performed on the deformed gouges from altered basalts suggest an increase in hydrothermal alteration with increasing temperature, as shown by a depletion in K2O at Tf ≥ 300 ˚C. SEM-BSE images of fine platy matrices in shear bands formed at Tf = 400 ˚C point towards a dissolution of quartz, pyroxene and plagioclase. Therefore, we suggest that the healing rate of both fresh and altered basalt not only scales with the ambient temperature but is also affected by the physical state of water, particularly in the case of fresh basalt, potentially related to more intense fluid-rock interactions with increasing temperature.

Keywords: frictional healing, frictional healing rate, hydrothermal fluids, basaltic gouge, Krafla geothermal field

How to cite: Wu, W.-H., Feng, W., Gomila, R., Tesei, T., Violay, M., Mortensen, A. K., and Di Toro, G.: Temperature and Physical State of Water Controls Frictional Healing of Basaltic Gouges from Krafla (Iceland), EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9165, https://doi.org/10.5194/egusphere-egu24-9165, 2024.

X2.153
|
EGU24-12986
|
ECS
Elvira Latypova, Fabio Corbi, Giacomo Mastella, Jonathan Bedford, and Francesca Funiciello

The short seismic record with respect to the return time of large subduction earthquakes and the spatial fragmentation of available geophysical data represent unfavourable conditions for robust hazard assessment. Over the last decade, data from scaled seismotectonic models have become useful in filling the observational gaps of seismic and geodetic networks. Such models allow reproducing hundreds of analogue seismic cycles in a few minutes of experimental time and with the advantage of known and controllable boundary conditions. 

Here we present experimental results from Foamquake – an established 3D seismotectonic model, which simulates megathrust subduction. Recent technical advances in experimental monitoring have allowed us to include into our research a high-frequency camera to record model surface deformation at 50 Hz and a network of 5 accelerometers (located on the model surface) that measure the three components of acceleration at 1 kHz. To analyse the camera data, we used particle image velocimetry (PIV) to derive surface displacements, such as in a dense, homogeneously distributed geodetic network spanning updip to scaled depths that are often offshore and, therefore, typically under-monitored in natural subduction zones.

We performed 33 experiments exploring 10 different geometrical configurations of asperities along the analog megathrust. In particular, we varied the number of asperities, their size, location, and extra normal load. We observed that the rupture pattern of analogue earthquakes predictably changes as the extra normal load varies and the distribution of asperity configurations becomes more complex. Depending on the number and size of the asperities and the size of the barrier between them, we noticed different ratios between full and partial ruptures with different recurrence time (Rt) intervals. In some experiments we detected cascades of ruptures. We used the coefficient of variation (CoV) of recurrence time to quantify analog earthquakes periodicity. Most of our models display quasi-periodic analog earthquakes recurrence with CoV<0.5, but multi-asperity experiments with variable-size and extra normal load lean toward random behaviour as testified by CoV~0.8.

Future investigations include the following steps – exploring this great volume of data using machine learning, looking for spatial and temporal relationships between accelerometer and PIV displacements, and tracking in detail the aseismic processes that may precede and follow earthquake rupture.

How to cite: Latypova, E., Corbi, F., Mastella, G., Bedford, J., and Funiciello, F.: Exploring earthquake recurrence and nucleation processes with Foamquake and a variety of asperity configurations, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12986, https://doi.org/10.5194/egusphere-egu24-12986, 2024.

X2.154
|
EGU24-6284
|
ECS
Raphael Affinito, Derek Elsworth, and Chris Marone

Elevated pore fluid pressures are frequently implicated in governing fault zone seismicity. While substantial evidence from geodetic and geological studies supports this notion, there is a notable scarcity of experimental observations of how fluid pressure influences fault stability during shear. Understanding the precise interplay between porosity, fault slip rate, and frictional stability is pivotal for assessing the significance of processes like dilational strengthening or thermal pressurization in the context of seismic hazards. Here, we prepare fault gouges from the Utah FORGE enhanced geothermal field injection well 16A at depths corresponding to seismic events (between 2050 – 2070m). Experiments were conducted inside a pressure vessel and loaded under a true-triaxial stress state, replicating in-situ stress conditions observed at the Utah FORGE site. The applied fault normal stress and during the experiments were held constant at 44MPa. Pore fluid pressure was varied between successive experiments (13, 20, and 27 MPa) to span a range of effective stresses to examine impacts on fault dilation/compaction and the successive frictional stability. Different fluid pressure boundary conditions: constant volume or pressure were applied to explore how changes in shearing rate influence gouge stability thought the fault drainage state. Our data indicate that the Utah FORGE samples are velocity-neutral and transition to velocity-weakening behavior at elevated pore pressure and shear strains >7. We find dilatancy coefficients e = ∆f/∆ln(v), where f is porosity and v is fault slip velocity, consistent with quartz-feldspathic-rich rocks ranging from 5–12^10-4, indicating a conditionally unstable regime. Furthermore, our results demonstrate that the boundary conditions for pore fluids influence frictional stability viachanges in effective normal stress. For example, when pore volume has zero flux, an expansion in the void volume during slip results in a decrease in pore pressure, transitioning the system towards frictional stability. Our results indicate that the connectivity of pore conduits may be more important than the imposed pore pressure conditions when considering the impact on fault stability. We suggest that the interplay between fault slip and fluid mobility within a fault is a delicate balance for predicting and managing seismic hazards.

How to cite: Affinito, R., Elsworth, D., and Marone, C.: Fault drainage state and frictional stability in response to shearing rate steps in natural gouge, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6284, https://doi.org/10.5194/egusphere-egu24-6284, 2024.

X2.155
|
EGU24-10969
Matt Ikari and Marianne Conin

The slip behavior of crustal faults is known to be controlled by the mineralogic composition of the fault gouge. The exact properties determining the frictional behavior of geologic materials, including diverse remains an important question. Here, we use a geochemical approach considering the role of water-rock interactions. As a mechanism, we suspect that the mineral surface charge allows attractive and repulsive forces (Van Der Waals type), and that those forces may influence the static mechanical behavior of clays (cohesion, static friction).  On the other hand, we suspect that the water bound to the mineral surfaces may play a role during shearing.  To address these ideas, we measured the cation exchange capacity (CEC) of 10 different rock and mineral types, including non-clays and a range of phyllosilicate minerals, using CEC as a proxy for the mineral surface charge and the ability to bind water to the mineral surfaces.  For these materials, we conducted laboratory shearing experiments measuring the pre-shear cohesion, peak friction coefficient, residual friction coefficient, post-shear cohesion, and velocity-dependent friction parameters under 10 MPa effective normal stress.  
Our results show that low CEC materials (< 3 mEq/100g) tend to exhibit high friction, low cohesion, and show velocity-weakening frictional behavior. The phyllosilicate minerals exhibit larger CEC values up to 78 mEq/100g and correspondingly lower friction coefficients, higher cohesion, and velocity-strengthening frictional behavior. Zeolite exhibits a relatively high CEC value typical of phyllosilicates, but its strength and frictional characteristics are that of a non-clay with low CEC. This suggests that grain shape and contact asperity size may be more important for non-phyllosilicates. For phyllosilicates, we suggest that the systematic patterns in strength and frictional behavior as a function of CEC could be explained by water bound to the mineral surfaces, creating bridges of hydrogen or van der Waals bonds when the particles are in contact. Such bonding explains the large cohesion values for high-CEC materials under zero effective stress, whereas surface-bound water trapped between the particles under load explains low friction.  Beyond the results of this study, CEC appears to be a controlling factor for other properties such as permeability and even the amount of bound DNA in sediments.

 

How to cite: Ikari, M. and Conin, M.: Cation Exchange Capacity Quantifies the Link Between Mineral Surface Chemistry and Frictional-Mechanical Behavior of Simulated Fault Gouges, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10969, https://doi.org/10.5194/egusphere-egu24-10969, 2024.

Multidisciplinary studies combining laboratory and/or numerical experimental results
X2.156
|
EGU24-13117
|
ECS
Ludovico Manna, Giovanni Toscani, Matteo Maino, Leonardo Casini, and Marcin Dabrowski

The 2D, plane strain, Finite Element Method-based linear elastic model that I present aims to assess the differential stress response to variations in the geometric configuration of a system of multiple collinear elliptic cracks intercepting a body of rock undergoing elastic deformation. The assumption underlying this simulation is that a collection of thin voids in a continuum medium can replicate the features observed in a system consisting of rough fault profiles in partial contact subjected to shear. The linear elastic model is designed to reproduce the stress and displacement fields around a rough fault, with a specific focus on stress concentration around its contact asperities. The model also allows to record the principal stress field on the domain for a wide range of scales and geometric properties of the system of collinear cracks embedded in the deforming rock. Analyzing the dependence of differential stress on parameters describing the geometry of rough fractures allows for considerations on the primary factors influencing brittle failure. Additionally, the examination of principal stresses around the tips of the cracks helps evaluate the potential orientation of new fracture patterns that may emerge when the yield strength of the deforming material is locally exceeded. The magnitude and orientation of the principal stresses are also crucial for the understanding of fracture coalescence and frictional reactivation of shear cracks in an elastic rock, which in turn is one of the main factors that govern the seismic cycle of natural faults. Furthermore, a comparison of the results of the present model with recent wing crack models of brittle creep suggest that our code may also be useful to obtain estimates of the critical distance between cracks for their interaction to coalesce into larger fractures. The process is assumed to indefinitely continue at greater scales, which offers the chance to propose a model for fault formation and propagation.

How to cite: Manna, L., Toscani, G., Maino, M., Casini, L., and Dabrowski, M.: A model for the formation and propagation of faults from the coalescence of smaller-scale systems of cracks: Finite Element Method-based numerical approach, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13117, https://doi.org/10.5194/egusphere-egu24-13117, 2024.

X2.157
|
EGU24-14003
Sylvain Barbot

Establishing a constitutive law for fault friction is a crucial objective of earthquake science. However, the complex frictional behavior of natural and synthetic gouges in laboratory experiments eludes explanations. Here, we present a constitutive framework that elucidates the slip-rate, state, temperature, and normal stress dependence of fault friction under the relevant sliding velocities and temperatures of the brittle lithosphere during seismic cycles. The competition between healing mechanisms explains the low-temperature stability transition from steady-state velocity-strengthening to velocity-weakening as a function of slip-rate and temperature. In addition, capturing the transition from cataclastic flow to semi-brittle creep accounts for the stabilization of fault slip at elevated temperatures. The brittle behavior is controlled by the real area of contact, which is a nonlinear function of normal stress, leading to an instantaneous decrease of the effective friction coefficient upon positive normal stress steps. The rate of healing also depends on normal stress, associated with an evolutionary response. If these two effects do not compensate exactly, steady-state friction follows a nonlinear dependence on normal stress. We calibrate the model using extensive laboratory data covering various relevant tectonic settings. The constitutive model consistently explains the evolving frictional response of fault gouge from room temperature to 600º for sliding velocities ranging from nanometers to millimeters per second, and normal stress from atmospheric pressure to gigapascals. The frictional response of faults can be uniquely determined by the in situ lithology and the prevailing hydrothermal conditions.

How to cite: Barbot, S.: Constitutive behavior of rocks during the seismic cycle in non-isothermal, non-isobaric conditions, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14003, https://doi.org/10.5194/egusphere-egu24-14003, 2024.

X2.158
|
EGU24-10564
Yu-Han Wang and Elías Rafn Heimisson

The interactions between fluids and fault structures play a pivotal role in understanding fault slip behavior. Over the years, various numerical methods have been developed to simulate these interactions. Volume-based methods, like the finite difference method (FDM), excel in their capacity to handle the intricacies of real-world fault structures, including material heterogeneity. On the other hand, the spectral boundary integral method (SBIM) is renowned for its computational efficiency. Recently, a hybrid approach has garnered significant attention, offering the benefits of both volume-based and SBI methods. This hybrid method allows for the consideration of fault structures' heterogeneity while maintaining computational efficiency. In this study, we introduce a novel hybrid method that bridges the SBIM and the FDM to model fluid migration in fault structures. Through rigorous model verification, we establish that our hybrid method can achieve a remarkable speedup of up to one thousand times compared to the FDM. Furthermore, we conducted two parametric studies to address open questions in fluid migration modeling within fault structures. First, we investigate the mobility contrast ratio between the host rock and the damage zone to determine the limits under which we can assume a zero-leak-off interface. Second, we explore the role of fault zone width in maintaining the validity of this zero-leak-off assumption. Building upon these foundational investigations, we demonstrate the possibility of extending the numerical framework to describe fault-fluid interactions considering poroelastic coupling.

How to cite: Wang, Y.-H. and Rafn Heimisson, E.: An efficient hybrid SBI-FD method for modeling fluid migration and fault-fluid interactions, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10564, https://doi.org/10.5194/egusphere-egu24-10564, 2024.

X2.159
|
EGU24-14504
Luyuan Huang, Luca Dal Zilio, and Elías Rafn Heimisson

Understanding earthquake rupture propagation across fault stepovers is pivotal for assessing the seismic hazard, offering vital insights into dynamic rupture processes within intricate fault geometries. However, the role of poroelastic effects within strike-slip fault systems featuring stepovers remains unexplored in dynamic models simulating Sequences of Earthquakes and Aseismic Slip (SEAS). Many existing models neglect poroelastic effects, and among those that consider them, a typical standard value of 0.8 is adopted for Skempton's coefficient B. Furthermore, a single dynamic rupture simulation is unable to address the frequency at which ruptures propagate through the stepover. Instead, these simulations only provide a binary status, indicating whether the ruptures jump or arrest. Thus, the investigation into how poroelasticity influences the likelihood of an earthquake jumping through a stepover emerges as a significant area of study. In response, we introduce a quasi-dynamic boundary element model that simulates 2D plane-strain earthquake sequences. This model incorporates undrained pore pressure responses affecting the fault's clamping and unclamping mechanisms and is governed by rate-and-state friction, with state evolution defined by the aging law. We first illustrate that dynamic rupture occurring in either left-lateral or right-lateral fault stepovers leads to a dynamic decrease (unclamping) or increase (clamping) in the effective normal stress. Dynamic variations of the effective normal stress depend on Skempton's coefficient. Consequently, higher Skempton's coefficients can promote rupture jumping across fault segments even for larger stepover distances. We then conduct a thorough parameter space study, evaluating the effects of Skempton's coefficient variations and stepover width on fault interactions within a fluid-filled porous environment. The likelihood of rupture jumping involves a trade-off between Skempton's coefficient and stepover width. We validate the numerical model by comparing it to an analytical solution that involves a plane strain shear dislocation on a leaky plane within a linear poroelastic, fluid-saturated solid. This validation demonstrates that a simple analytical solution, primarily dependent on fault dislocation and Skempton's coefficient, has the potential to effectively predict the pore pressure change. The critical jumping width for 50% chance of rupture jumping predicted by our model explains the threshold dimension of the fault step, above which ruptures do not propagate. This study highlights the significance of incorporating poroelastic effects on- and off-fault in understanding the dynamic variations of the effective normal stress, which could significantly alter the overall length of fault rupture.

How to cite: Huang, L., Dal Zilio, L., and Rafn Heimisson, E.: The role of poroelasticity in rupture dynamics across fault stepovers, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-14504, https://doi.org/10.5194/egusphere-egu24-14504, 2024.

X2.160
|
EGU24-7803
|
ECS
Betti Hegyi, Taras Gerya, Luca Dal Zilio, and Whitney Behr

The role of fluid flow in triggering earthquakes in subduction zones is a critical yet complex aspect in seismology. Despite extensive study through geological, geophysical observations, and laboratory experiments, fully understanding and modelling these processes within a coupled solid-fluid interaction framework remain challenging. This study employs a coupled seismo-hydro-mechanical code (i2elvisp) to simulate fluid-driven earthquake sequences in a simplified subduction megathrust environment. We incorporate non-uniform grid resolution, enhancing the resolution of seismic events within the subduction channel. The code integrates solid rock deformation with fluid dynamics, solving mass and momentum conservation equations for both phases, alongside gravity and temperature-dependent viscosity effects. Brittle/plastic deformation is modelled through a rate-dependent strength formulation, with slip instabilities governed by compaction-induced pore fluid pressurisation. Our approach demonstrates the significant impact of fluid pressurisation on deformation localization, achieving slip rates up to metres per second in a fully compressible poro-visco-elasto-plastic medium. By refining the vertical model resolution in the subduction channel to less than or equal to 200 metres, we ensure convergence in terms of event recurrence interval and slip velocity. The models successfully replicate various slip modes observed in nature, ranging from regular earthquakes (including partial and full ruptures) to transient slow slip phenomena and aseismic creep. This research focuses on the parameters influencing the dominant slip mode, their distributions, and interactions along a modelled subduction interface. Our findings indicate that the dominant slip mode and the earthquake sequences are significantly influenced by porosity, permeability, and temperature-dependent viscosity. We explore two distinct viscosity gradients in the subduction channel to represent subduction zones with differing thermal profiles. In 'hot' subduction models, the brittle-ductile transition commences at shallower depths than in 'cold' subduction cases, influencing the nucleation depth of seismic events. These viscosity variations markedly impact model evolution; regular earthquakes exhibit higher velocities and slip rates in 'hot' scenarios, which are also more conducive to hosting aseismic creep or slow slip events. In conclusion, our study elucidates the pivotal role of fluid pressure evolution in seismicity within subduction zones and provides deeper insights into earthquake source processes. Through comprehensive modelling and analysis, we enhance understanding of the complex dynamics governing fluid-induced seismic activity and contribute to the broader field of earthquake source processes. 

How to cite: Hegyi, B., Gerya, T., Dal Zilio, L., and Behr, W.: Fluid driven seismic cycle modelling in subduction zones, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7803, https://doi.org/10.5194/egusphere-egu24-7803, 2024.

X2.161
|
EGU24-5767
|
ECS
Pierre Romanet, Marco Scuderi, Stéphanie Chaillat, Jean-Paul Ampuero, and Frédéric Cappa

Numerical modeling of Discrete Fracture Networks (DFNs) is commonly used to assess the behavior and properties of hydraulic diffusion and seismicity in the Earth’s crust within a network of fractures and faults, and to study the hydromechanical evolution of fractured reservoirs stimulated by hydraulic injection and production. The modelling of such fractures is typically carried out under a quasi-static approximation, and occasionally accounting for elasto-dynamics in single-rupture studies that assume a slip-weakening friction law. 

In this work, we develop a 2D DFN model to simulate fluid-induced seismicity that couples hydraulic diffusion and slip governed by rate-and-state friction on multiple interacting faults. The main goal of this numerical model is to establish a connection between laboratory derived friction parameters and field observations, enabling the inference of the long-term evolution of fractured reservoirs and crustal fault systems undergoing multiple earthquakes and (slow) slip events induced by fluid pressure perturbations.

In the model, the elastic interactions are computed with a boundary element method, accelerated by the hierarchical matrix method. We assessed the convergence of the method at fracture junctions and verified it does not create unphysical singularities. The use of rate-and-state friction makes it possible to model several seismic events over the injection duration.

The simulations will be later used to fit measurements of permeability and friction collected in laboratory experiments, in-situ observations of fault slip and opening from fluid injection experiments at decametric scale, and finally, induced seismicity at reservoir scale.

 

How to cite: Romanet, P., Scuderi, M., Chaillat, S., Ampuero, J.-P., and Cappa, F.: Towards a 2D model of Discrete Fracture Network with permeability and friction evolution for modeling fluid-induced seismicity , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5767, https://doi.org/10.5194/egusphere-egu24-5767, 2024.