Slip along pre-existing faults in the Earth’s crust occurs whenever the shear stress resolved on the fault plane overcomes fault frictional strength, potentially generating catastrophic earthquakes. The coupling between shear stress and normal stress during fault loading depends on 1) the orientation of the fault within the stress field and 2) the tectonic setting. In compressional settings, a load-strengthening path occurs because along thrust faults the increase in shear stress is coupled with an increase in effective normal stress. On the contrary, in extensional settings, the increase in shear stress is coupled with a decrease in effective normal stress, resulting in load-weakening paths for normal faults.

Analytical approaches to evaluate the potential for fault reactivation are generally based on the assumption that faults are ideal planes, characterized by zero thickness and constant friction, embedded in homogeneous isotropic elastic media. However, natural faults typically host thick fault cores and highly fractured damage zones, which can compact or dilate under different loading paths (i.e., different coupling between normal and shear stress). In addition, in most laboratory friction experiments, the fault is loaded under constant or increasing normal stress and at optimal orientation for reactivation. Here, we present laboratory experiments simulating reactivation of thick gouge-bearing faults that experienced different loading paths.

Our results show that the differential stress required for reactivation strongly differs from theoretical predictions, and unfavourably oriented faults appear systematically weaker, especially when a thick gouge layer is present. Before reactivation fault zone compacts in load-strengthening paths whereas dilation is observed in load-weakening path. Upon fault reactivation at comparable normal stress, load-strengthening promotes stable creep  whereas load-weakening results in accelerated slip. Our study highlights the importance of fault thickness and loading path in fault hydromechanical coupling and stability with significant implications for fluid circulation within fault zones and earthquake mechanics.

How to cite: Giorgetti, C., Violay, M., and Collettini, C.: The role of loading path on fault reactivation: a laboratory perspective, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10016, https://doi.org/10.5194/egusphere-egu23-10016, 2023.

Catastrophic failure in brittle, porous materials initiates when structural damage, in the form of smaller-scale fractures, localises along an emergent failure plane or 'fault' in a transition from stable crack growth to dynamic rupture. Due to the extremely rapid nature of this critical transition, the precise micro-mechanisms involved are poorly understood and difficult to capture. However, these mechanisms are crucial drivers for devastating phenomena such as earthquakes, including induced seismicity, landslides and volcanic eruptions, as well as large-scale infrastructure collapse. Here we observe these micro-mechanisms directly by controlling the rate of micro-fracturing events to slow down the transition in a unique triaxial deformation experiment that combines acoustic monitoring with contemporaneous in-situ x-ray imaging of the microstructure. The results [1] provide the first integrated picture of how damage and associated micro-seismic events emerge and evolve together during localisation and failure and allow us to ground truth some previous inferences from mechanical and seismic monitoring alone. They also highlight where such inferences miss important kinematically-governed grain-scale mechanisms prior to and during shear failure.

The evolving damage imaged in the 3D x-ray volumes and local strain fields undergoes a breakdown sequence involving several stages: (i) self-organised exploration of candidate shear zones close to peak stress, (ii) spontaneous tensile failure of individual grains due to point loading and pore-emanating fractures within an emergent and localised shear zone, validating many inferences from acoustic emissions monitoring, (iii) formation of a proto-cataclasite due to grain rotation and fragmentation, highlighting both the control of grain size on failure and the relative importance of aseismic mechanisms such as crack rotation in accommodating bulk shear deformation. Dilation and shear strain remain strongly correlated both spatially and temporally throughout sample weakening, confirming the existence of a cohesive zone, but with crack damage distributed throughout the shear zone rather than concentrated solely in a breakdown zone at the propagating front of a pre-existing discontinuity.

Contrary to common assumption, we find seismic amplitude is not correlated with local imaged strain; large local strain often occurs with small acoustic emissions, and vice versa. The seismic strain partition coefficient is very low overall and locally highly variable. Local strain is therefore predominantly aseismic, explained in part by grain/crack rotation along the emergent shear zone. The shear fracture energy calculated from local dilation and shear strain on the fault is half of that inferred from the bulk deformation, with a smaller critical slip distance, indicating that less energy is required for local breakdown in the shear zone compared with models of uniform slip.

This improvement in process-based understanding holds out the prospect of reducing systematic errors in forecasting system-sized catastrophic failure in a variety of applications.

[1] Cartwright-Taylor et al. 2022, Nature Communications 13, 6169, https://doi.org/10.1038/s41467-022-33855-z

How to cite: Cartwright-Taylor, A., Mangriotis, M.-D., Main, I. G., Butler, I. B., Fusseis, F., Ling, M., Andò, E., Curtis, A., Bell, A. F., Crippen, A., Rizzo, R. E., Marti, S., Leung, D. D. V., and Magdysyuk, O. V.: Micromechanics of damage localisation and shear failure of a porous rock: sound and vision, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7933, https://doi.org/10.5194/egusphere-egu23-7933, 2023.

During earthquake propagation, a shock wave damages rocks at the rupture tip, creating numerous microfractures and altering the mechanical properties of fault zone rocks. This damage, which occurs dynamically at the millisecond time scale, controls rock strength during earthquake slip that occurs in the wake of rupture propagation. How the presence of water and the initial porosity of the rock control damage during high strain rate deformation remains an open question. We have performed a series of shock experiments using a split Hopkinson pressure bar apparatus installed at the European Synchrotron Radiation Facility. Using two ultra-fast cameras synchronized with the X-ray bunches of the synchrotron; we imaged deformation with microsecond time resolution on centimetre-scale core samples during shock wave damage. We deformed dry and water saturated low porosity Westerly granite and porous Berea sandstone samples. Several samples were surrounded by a thin aluminium jacket allowing recovering them after deformation and image them using X-ray microtomography with micrometre spatial resolution. Results confirm previous studies that have shown that rock pulverization occurs above a threshold strain rate produced by the shock wave. Water saturated samples are consistently weaker than dry samples as they pulverize under lower peak stress. Analyses of rock microstructure acquired using the ultrafast cameras and X-ray microtomography data shed light on the micro-mechanisms of damage production. Either the entire sample pulverized (Westerly granite) or a compaction of the sample occurred before shear zones were dynamically produced (Berea sandstone). These results demonstrate fundamental differences in dynamic damage production in crystalline and porous dry and wet rocks. Our data unravel mechanisms of gouge production before any significant slip has occurred on a fault, which control the shear strength during earthquake slip.

How to cite: Renard, F., Cordonnier, B., Doan, M.-L., Fondriest, M., Lukic, B., and Prastyani, E.: Dynamic damage in dry and wet rocks monitored by ultra-fast synchrotron imaging, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8778, https://doi.org/10.5194/egusphere-egu23-8778, 2023.

Preparatory earthquake processes such as slow preparatory slip (preslip) are connected to variations in frictional strength linked to frictional instabilities and appear in various scales across the Earth’s crust. For dry and bare surfaces, the fault surface characteristics affect the contact conditions. These conditions are established through asperities, which are topographical heights where the normal stress concentrates, imposing variations in fault strength. The effect of surface conditions on preslip can be studied in the laboratory where fault surface characteristics can be identified. Developing a more refined understanding of features controlling preslip (e.g., roughness) will lead to more realistic models describing frictional stability. In this study, we performed a triaxial test at sequentially increasing confining pressure steps (P_c = 60, 80, 100 MPa) on a saw-cut sample of Carrara Marble in dry and unlubricated conditions. Two types of technologies were used to study this frictional response in space and time: (1) an array of acoustic emission sensors monitored localized precursory seismicity and (2) quasi-static deformation in the fault-parallel strain was monitored using distributed strain sensing (DSS) with fiber optics. The differential stress was also measured throughout and allowed us to study the onset of frictional weakening/strengthening. In the first confining pressure step (P_c = 60 MPa), a single stick-slip event was observed with an associated 43 MPa static stress drop. In the subsequent confining pressure steps of P_c = 80 and 100 MPa, even though the normal stress on the fault was increased, no stick-slip events were observed, and the fault smoothly transitioned to sliding with smaller magnitude stress drops of 3 and 4 MPa, respectively. That suggests that a change in the frictional nature of the interface was incurred during the first rupture at P_c = 60 MPa. The high-density DSS array displayed a significant heterogeneous distribution of fault-parallel strain in time and space and experienced sudden reorganization at various phases of the experiment. Due to the high spatial resolution, DSS allowed us to investigate local deviations from an elastic response attributed to inelastic processes. A larger amount of local strain accumulation was needed to produce a stick-slip instability. At higher normal stress on the pre-ruptured fault, this level of locking was not possible in the subsequent confining pressure steps. Dissipative inelastic deformation was attributed to local frictional weakening that resulted in non-uniform preslip. Furthermore, priori measurements of contact pressure heterogeneities were obtained using a pressure sensitive film. These results showed regions of lower normal stress along the fault that correlated with regions that incurred more anelastic response on the DSS array. Post-mortem contact pressure measurements showed clear changes in the normal stress distribution that correlated to visual damage and wear. We believe that this contributed to the fault's inability to lock as before and mitigate dynamic rupture. Our results provide more insight into potential mechanisms controlling preslip distribution leading to dynamic and quasi-static frictional weakening.

How to cite: Michail, S., Selvadurai, P. A., Cebry, S. B. L., Salazar Vásquez, A. F., Bianchi, P., Rast, M., Madonna, C., and Wiemer, S.: Laboratory Observations linking Fault Surface Characteristics to Preparatory Earthquake Processes and Fault Stability, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8634, https://doi.org/10.5194/egusphere-egu23-8634, 2023.

We propose an approach that is aimed to enrich the catalogs of acoustic emission events recorded in laboratory experiments with such parameters as seismic moment and corner frequency. Because of the difficulty of separation of direct waves in experiments performed on small samples, we use the coda waves that are composed of the reverberation of the acoustic field in the tested sample. After multiple reverberations, the resulting wavefield can be approximated as nearly homogeneously distributed over the sample and with signal amplitudes decaying exponentially in time (linearly in a logarithmic scale).

Within the framework of this model, the frequency-dependent coda amplitude at any moment of time is described as combination of a source spectra, of a decay rate combining internal attenuation with reverberation losses, and of a sensor response. One of the main difficulties with the laboratory experiments is that acoustic sensors are very difficult to calibrate and their absolute response function in most of cases remains unknown. With the simple reverberation model, the logarithms of coda amplitudes at different times and sensors and for multiple events are described by a system of linear equations that we solve in a least-square sense to find frequency-dependent coda-decay rates, relative signal spectra and sensor responses. In a next step, we compute spectral ratios between spectra of different events to eliminate the sensor responses and to estimate main source parameters such as corner frequencies and relative seismic moments.

We provide details of our data analyses technique and present time-dependent corner frequency vs relative moment diagrams for two experiments on granite of the Voronezh massif and Berea sandstone under pseudo-triaxial loading. The dependence close to the cubic that is frequently estimated for tectonic earthquakes observed on the first stages of both experiments when confining pressure steps applied to the intact rock and therefore to the pre-existing inhomogeneties. After applying axial load changes in stress-drop is being observed: with higher stress-drops prevailing in granite and lower stress-drops in sandstone. Also there is a significant difference in Gutenberg-Richter relation in these two experimental conditions observed.

How to cite: Kartseva, T. and Shapiro, N.: Coda-Based Estimation of Source Parameters of Laboratory Acoustic-Emission Events, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-1330, https://doi.org/10.5194/egusphere-egu23-1330, 2023.

The elastic medium that hosts several, multi-scale, faults could be regarded as an energy reservoir that is charged by the far field stress rate and discharged by friction dissipation during earthquake slips on the faults.  In this study, we carefully analyze the energy budget variations that occur throughout a synthetic, 2D plane-strain, earthquake cycle on a fault system comprising of a main fault surrounded by a hierarchy of off-fault slip planes/fractures. We evaluate the rate of kinetic energy variation, stress power across the continuum, far field power supply, and the dissipation due to the rate-and-state friction on the faults given a spectrum of slips ranging from slow-slip to rapid ruptures. We study how the medium's energy budget evolves after these components have been determined.  Additionally, we compute the dissipation rate for a variety of slip rates to determine the contributions of so-called slow-slip events, low-frequency earthquakes (LFEs), and tremors to this budget. We also evaluate the share of off-fault fractures to determine their energetic role during earthquake cycles.

How to cite: Kheirdast, N., Almakari, M., Villafuerte, C., Thomas, M. Y., and Bhat, H. S.: Energy budget of quasi-dynamic earthquake cycle, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3982, https://doi.org/10.5194/egusphere-egu23-3982, 2023.

We propose a numerical model of laboratory earthquake cycle inspired by a set of experiments performed on a triaxial apparatus on sawcut Carrara marble samples. The model couples two representations of rock matter: rock is essentially represented as an elastic continuum, except in the vicinity of the sliding interface, where a discrete representation is employed. This allows to simulate in a single framework the storage and release of strain energy in the bulk of the sample and in the loading system, the damage of rock due to sliding, and the progressive production of a granular gouge layer in the interface. After independent calibration, we find that the tribosystem spontaneously evolves towards a stick-slip sliding regime, mimicking in a satisfactory way the behaviour observed in the lab. The model offers insights on complex phenomena which are out of reach in experiments. This includes the variability in space and time of the fields of stress and effective friction along the fault, the progressive thickening of the damaged region of rock around the interface, and the build-up of a granular layer of gouge accommodating shear. We present in detail several typical sliding events, we illustrate the fault heterogeneity, and we analyse quantitatively the damage rate in the numerical samples.

How to cite: Mollon, G., Aubry, J., and Schubnel, A.: On-fault damage evolution in laboratory earthquakes: a numerical perspective on fault complexity, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6517, https://doi.org/10.5194/egusphere-egu23-6517, 2023.

In addition to regular earthquakes, observations of spatiotemporally complex slip events have multiplied over the last decades. These slip events range along different time scales: from creep , slow slip events to LFEs and tremors. At present, these events are generally interpreted by imposed frictional heterogeneities along the fault plane. However, fault systems are geometrically complex in nature over different scales. We aim in this work to investigate the role of “realistic” fault geometry on the dynamics of slip events. We consider a fault system in a 2D quasi-dynamic setting. The fault system consists of a main self-similar rough fault, surrounded by a dense network of off-fault fractures. All fractures are frictionally homogeneous (rate weakening) and can potentially undergo dynamic slip. We aim to understand how the deformation in the volume is accomodated by the off-fault damage zone and the main fault. What fraction of the “supplied” moment rate is hosted by the off-fault fractures during an earthquake cycle?

How to cite: Bhat, H., Almakari, M., Kheirdast, N., Villafuerte, C., and Thomas, M.: Fault zone complexity naturally produces the full slip spectrum: Insights from numerical models, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16290, https://doi.org/10.5194/egusphere-egu23-16290, 2023.

Faults are complex systems embedded in an evolving medium fractured by seismic ruptures. This off-fault damage zone is shown to be thermo-hydro-mechano-chemically coupled to the main fault plane by a growing number of studies. Yet, off-fault medium is still, for the most part modelled as a purely elastic -- hence passive -- medium. Using a micromechanical  model we investigate the depth variation of dynamically triggered off-fault damage and its counter-impact on earthquake slip dynamics. We show that if the damage zone becomes narrower with depth, it is also denser and thus, unlike what is commonly believed, remains an energy sink even at depth. The results are in agreement with the complementary model by Okubo et al., 2019. In contrast to study cited above, our model accounts for the dynamics changes of elastic moduli related to crack growth. This lead to the dynamic creation of low-velocity zone that can trapped seismic waves and further impact the earthquake dynamics, even at greater depth. We therefore claim that the intertwined dynamics between the main fault plane and its surrounding medium must be including along the all seismogenic.

How to cite: Ferry, R., Thomas, M., and Jeandet, L.: Depth dependence of coseismic off-fault damage, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6974, https://doi.org/10.5194/egusphere-egu23-6974, 2023.

Major earthquake ruptures occur predominantly in thrust faults producing devastating events and tsunamis such as the 2011 Mw 9.0 Tohoku earthquake, the 2004 Mw 9.2 Sumatra earthquake and the 1999 Mw 7.7 Chi-Chi earthquake. Understanding the mechanics of earthquakes in thrust faults and the effect of the free surface is thus crucial to explain their large shallow slip, their asymmetric ground motion and their damage patterns surrounding the fault and the free surface. In this work, we carry out 2D dynamic rupture simulations on thrust faults to accurately characterize a possible unclamping effect, its responsible physical mechanism, and to produce dynamically activated off-fault fracture networks. To conduct the simulations, we use the software tool based on the Combined Finite-Discrete Element Method (FDEM), HOSSedu, developed by Los Alamos National Laboratory. Our dynamic rupture models in an elastic medium confirm that unclamping occurs in thrust faults and increases significantly as the rupture reaches the free surface and for the fault models with lower dip-angles. We show that this is a consequence of the torque mechanism induced in the hanging wall, and the release of this torque when the rupture reaches the free surface produces a “flapping” in the toe of the wedge where the most significant unclamping (possibly leading to fault opening) is taking place. Our results indicate that the free surface produces a considerable reduction of the compressive normal stress when the rupture is propagating up-dip that facilitates the extension and the amount of slip close to the trench as observed for large thrust earthquakes.This significant normal stress change is reflected in the orientation of the principal stresses before and after the rupture, where under certain conditions, the greatest principal stress changes from subhorizontal to almost vertical leading to a post-rupture tensional stress state in the hanging wall that has been confirmed by observations of recent in-situ, seismological and geodetic studies. Finally, we investigate whether this dramatic normal stress reduction stands when we allow for the activation of coseismic off-fault damage and explore its role in the rupture dynamics, the near-field deformation and radiation patterns.

How to cite: Villafuerte, C., Okubo, K., Rougier, E., Madariaga, R., and Bhat, H. S.: Dynamics and radiation of thrust earthquakes with coseismic off-fault damage, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7849, https://doi.org/10.5194/egusphere-egu23-7849, 2023.

In order to make seismic hazard estimates, it is necessary to assume some distribution of the number of earthquakes of a certain magnitude, i.e. a Gutenberg-Richter distribution. This is true for both natural seismicity as well as induced seismicity, but in both cases the number of historical earthquakes at the tail end of the distribution (i.e. the largest ones) is limited and often the maximum possible magnitude is unknown. In contrast, in a laboratory setting the maximum size of an unstable slip event (“stick-slip” or laboratory earthquake) is controlled by the size of the sample and the imposed stress. In our rotary shear apparatus, we can theoretically achieve unlimited fault displacement which allows us to produce earthquake-like distributions with thousands to tens of thousands event.

In this presentation, I will present results from experiments on simulated fault gouges which show unstable frictional behaviour at room temperature conditions. The results show that the event size distribution can change spontaneously, without any changes in the boundary conditions. Observations of fault gouge material after the experiment, suggest that wear of the granular material generates alternative surfaces for slip, which changes the macroscopic behaviour. Interestingly, the change in event size distribution is reversable, presumably because the fine-grained layers become disturbed with ongoing shear.

In an attempt to simulate the macroscopic behaviour, we have, for the first time, measured the rate-and-state frictional (RSF) properties on single grain contacts. Using these values in a numerical model for seismic slip (so-called “seismic cycle simulator”), we obtain maximum stress drops that are comparable to those obtained in the experiments, but with some differences. The differences are most likely due to the fact that the grains in our simulated fault are affected wear in previous slip events, which should change their RSF parameters. In addition, the normal stress at each individual grain contact is unknown in the experiment and could vary significantly from event to event.  This latter difference between model and experiment can be overcome by using a discrete element method with contact-scale RSF properties to simulate slip.

How to cite: Niemeijer, A. R., Korkolis, E., Mishra, T., Elbertsen, R., Jop, B., and Pires de Vasconselos, I.: Microstructural controls on seismicity distribution in simulated fault gouges, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2911, https://doi.org/10.5194/egusphere-egu23-2911, 2023.

Faults can act as conduits for the migration of hydrothermal fluids in the crust, affecting its mechanical behaviour and possibly leading to earthquake swarm activity. To date, there are still few constraints from the geological record on how fault-vein networks develop through time in high fluid-flux tectonic settings. Here, we describe small displacement (<1.5 m) epidote-rich fault-vein networks cutting granitoids in the exhumed Bolfin Fault Zone (Atacama Fault System, Chile). The epidote-rich sheared veins show lineated slickensides with scattered orientations and occur at the intersections with subsidiary structures in the fault damage zone. FEG-SEM cathodoluminescence (CL) reveals that magmatic quartz close to the sheared epidote-rich veins is affected by (i) thin (< 10 µm) interlaced deformation lamellae and (ii) a network of CL-dark quartz epitaxial veinlets sharply crosscutting the lamellae. EBSD maps of the deformed quartz indicate minor lattice distortion associated with the lamellae and an orientation nearly orthogonal to the c-axis. These deformation features disappear moving away into the host rock. The epidote-rich sheared veins (i) include clasts of magmatic quartz with both the deformation lamellae and the healed veinlets and (ii) show cyclic events of extensional-to-hybrid veining and localized shearing. We propose that the microstructures preserved in the quartz next to the sheared veins (i.e. deformation lamellae and epitaxial veinlets) record the high-strain rate loading associated with dynamic crack propagation and rapid micro-fracture sealing. On the other hand, the cyclic dilation and shearing within the epidote-rich veins is interpreted as the expression of a highly connected fault-vein network dominated by pore pressure oscillations leading to seismic swarm activity.

How to cite: Masoch, S., Fondriest, M., Gomila, R., Pennacchioni, G., Cembrano, J., and Di Toro, G.: Interplay between fluid flow and rock deformation in an exhumed hydrothermal fault-vein network, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-531, https://doi.org/10.5194/egusphere-egu23-531, 2023.

Fault frictional healing Δμ controls the storage of the elastic strain energy in the fault wall rocks and the re-occurrence of earthquakes in pre-existing faults. In the last 40 years, fault healing has been investigated with laboratory slide-hold-slide (SHS) experiments aimed at reproducing the seismic cycle. Experiments performed with different rocks types (e.g., granite, limestone, basalt) revealed that (1) Δμ increases with hold time t_h and, (2) the frictional healing rate β =Δμ/log t_h >0. This increase in fault frictional strength with t_h is interpreted as due to the increase (1) in the real area of contact or (2) of chemical bond strength. However, most of these experiments were conducted under room conditions, whereas natural earthquakes generally nucleate at ambient temperatures T  >150℃ and in the presence of pressurized fluids. Under these ambient conditions, fluid-assisted and thermally-activated processes (pressure-solution transfer, stress corrosion, etc.) may impact on the magnitude of Δμ and on β.

In this study, SHS experiments were performed on gabbro-built gouges (grain size <88 mm) in a rotary shear machine equipped with a pressurized vessel to explore frictional healing under hydrothermal conditions. All experiments were conducted at a constant effective normal stress (σ_eff =50MPa), and temperature (T) ranging from 25 to 400 ℃  under dry or pore fluid (deionized H₂O) pressure (P_f=30 MPa) conditions. In the SHS sequence, the imposed slip velocity was V=10 μm/s, and hold time t_h varied from 3 to 10000 s. For each experiment, two SHS sequences separated by a slip displacement interval of 40 mm were conducted.

Under dry conditions at all tested temperatures and under hydrothermal conditions but at T  <100℃, Δμ increases with t_h, consistent with previous experiments. Moreover, the Δμ and β values in the 2^nd SHS sequence are slightly higher than those in the 1^st sequence, possibly due to the smaller grain size at the larger displacement that promotes fault healing. By contrast, in the experiments performed under hydrothermal conditions but T >200℃, Δμ decreases and β switches to negative values (<0) when the hold is longer than a threshold hold time. In detail, at T=300℃: β= 0.0161±0.0017 for holds <300s and -0.0074±0.0043 for holds >300s, and at T=400℃: β= 0.0057±0.0020 for holds <100s and -0.0227±0.0042 for holds >100s.

The underlying mechanism responsible for the decrease in Δμ and the transition from β > 0 to β < 0 with the hold time, which could result in the transition from seismic to aseismic fault behavior in nature, is still poorly understood. However, high-resolution microstructural analyses conducted by scanning electron microscopy on experimental fault products rule out the formation of weak minerals (e.g., clays) in the gouge layer.  Consequently, the weakening of the fault is probably related to the decrease in bond strength at the asperity contacts.

The experimental data presented here suggest that fault healing of natural faults is controlled by the feedback of multiple physico-chemical processes associated with the slip history and type of fluid-rock interaction under hydrothermal conditions.

How to cite: Feng, W., Yao, L., Gomila, R., Ma, S., and Di Toro, G.: Healing of gabbro-built faults under hydrothermal conditions, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9616, https://doi.org/10.5194/egusphere-egu23-9616, 2023.

Pore fluids are ubiquitous throughout the lithosphere and are commonly cited as a major factor producing slow slip and complex modes of tectonic faulting. Here, we investigate the role of pore pressure on slow slip and the frictional stability transition and find that the mode of fault slip is largely unaffected by pore pressure once we account for effective stress. Ambient temperature experiments are done on synthetic fault gouge composed of quartz powder with a median grain size of 10μm with an average permeability of  8E-17m² – 6E-18m² from shear strains 0 - 26. We conduct constant velocity experiments at 20MPa σ_n’, with P_p/σ_n’ratios of λ from 0.05 to 0.28. Under these conditions, dilatancy strengthening is minimal and we find that slip rate dependent changes in the critical rate of frictional weakening are sufficient to explain slow slip.

How to cite: Affinito, R., Elsworth, D., and Marone, C.: The Stability Transition from Stable to Unstable Frictional Slip with Finite Pore Pressure, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16777, https://doi.org/10.5194/egusphere-egu23-16777, 2023.

In this paper, based on the model of thermal pressurization, we present a new way for the emergence of rate and state phenomenology (RSF, friction law) during the earthquake cycle. In the framework of fault mechanics, the common physical mechanism for the RSF phenomenology is slip and plastic deformation at the asperity contacts. We show that the fundamental physical mechanism of thermal pressurization together with viscosity inside the fault can also reproduce rate and state phenomenology.

More specifically, in our numerical analyses we model frictional weakening during large seismic slip due to thermal pressurization inside the fault. We introduce thermo-hydro-mechanical couplings to model thermal pressurization and a first order micromorphic Cosserat continuum, in order to avoid mesh dependence of the numerical results. Moreover, we introduce viscosity in the form of strain rate hardening. When we perform velocity stepping analyses, our numerical findings show that friction presents, initial peak over-strength and frictional oscillations around a residual value (see Figure). Our results, deriving from fundamental modeling assumptions, exhibit rate and state phenomenology, without the need to introduce the physical mechanism of slip at the asperity contacts.

Keywords: THM couplings; Viscosity; Cosserat continuum; Tribology; Earthquakes

How to cite: Stathas, A. and Stefanou, I.: Viscosity and thermal pressurization during large seismic slip lead to rateand state phenomenology, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14607, https://doi.org/10.5194/egusphere-egu23-14607, 2023.

The question "what arrests an earthquake rupture?" sits at the heart of any potential prediction of earthquake magnitude. Here, we present a one-dimensional, thin-elastic-strip, minimal model, to illuminate the basic physical parameters that control the arrest of large ruptures. The generic formulation of the model allows for wrapping various earthquake arrest scenarios into the variations of two dimensionless variables, valid for both in-plane and antiplane shear loading. Our continuum model is equivalent to the standard Burridge-Knopoff model, with an added characteristic length scale, that corresponds to either the thickness of the damage zone for strike-slip faults or to the thickness of the downward moving plate for subduction settings. We simulate the propagation and arrest of frictional ruptures and present closed-form expressions to predict rupture arrest under different conditions. Our generic model illuminates the different energy budget that mediates crack- and pulse-like rupture propagation and arrest. Despite its simplicity, this minimal model is able to reproduce several salient features of natural earthquakes that are still debated (e.g. various arrest scenarios, stable pulse-like rupture, back-propagating front, asymmetric slip profiles).

How to cite: Barras, F., Thøgersen, K., Aharonov, E., and Renard, F.: How do earthquakes stop? Insights from a minimal model of frictional rupture, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16320, https://doi.org/10.5194/egusphere-egu23-16320, 2023.

Although the physical mechanism of earthquake nucleation processes and the link with foreshocks are under debate, foreshocks are still considered as the most reliable earthquake precursors. Investigating the temporal and spatial evolution of foreshock sequences with high resolution and monitoring b-values in real time may shed light on these key issues. Many foreshock and aftershock sequences accompanying moderate mainshocks have been reported in the west of Yunnan Province, China, such as the 2016 Yunlong M 5.1 and 2021 Yangbi Ms 6.4 earthquake sequences. The recently improved coverage of seismic network in western Yunnan provides the opportunity to investigate how the foreshock sequence evolved and establish the temporal transient in b values. To find missing earthquakes and built more comprehensive earthquake catalogs, we carried out earthquake detection using the matched-filter detector. We used events in the standard catalog of China Earthquake Networks Center as templates to scan through continuous waveforms 3-6 months before and after the main shock. We then estimated the b-value and its temporal changes based on the newly developed catalogs. An obvious reduction in b-values before the major earthquake is observed in both the 2016 Yunlong and 2021 Yangbi sequences. We also found that the scattered spatial pattern of foreshocks exhibits a cascading manner and does not support the hypothesis of slow slip driving nucleation of mainshocks.

How to cite: Zhu, G. and Yang, H.: Foreshocks preceding moderate earthquakes in Western Yunnan, China, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10314, https://doi.org/10.5194/egusphere-egu23-10314, 2023.

Investigating clusters of events and geophysical screening often provides limited constraints on fault geometries. This imaging issue prevents the integration of realistic fault zone geometry into earthquake studies, which can affect our capacity to evaluate the role of pre-existing faults on seismicity. This study describes a modelling strategy accounting for realistic fault zone geometries. Our approach uses stochastic methods underpinned by quantitative fault zone parameterization, followed by an assessment of seismicity from simulations of rupture dynamics controlled by fault geometry. This method is used to investigate the role of fault maturity on seismicity for two case studies, including seismicity associated with the reactivation of a pre-existing fault network due to hydraulic fracturing in Harrison County, Ohio, from 2013 to 2015, as well as the natural seismicity associated with the Yushu-Ganzi left-lateral strike-slip fault system in central-eastern Tibet. In the Harrison County case, we analyze the effect of vertical variability in fault maturity and show how more mature faults in the deeper crystalline basement generate higher magnitude seismicity than shallow, immature faults in younger sedimentary rocks. In the Yushu-Ganzi case study, we show how lateral variability in structural maturity, arising from long-term fault propagation and strain rates, leads to different seismicity on individual fault segments. Our findings indicate that fault geometry determines seismic patterns, with rupture length controlled by fault zone geometry rather than fault lengths, and favour the adoption of a structural geological perspective for the integration of realistic fault geometry into seismicity prediction strategies.

How to cite: Roche, V., van der Baan, M., and Walsh, J.: Modelling seismicity based on fault geometry: maximum magnitudes and magnitude-frequency distributions., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5033, https://doi.org/10.5194/egusphere-egu23-5033, 2023.

Surface rupture produced by earthquakes can pose great threat on near-surface infrastructures and elevate damages. Accessing the potential of surface rupture along faults is critical to mitigating such hazards. It is commonly suggested that earthquakes with Mw>6.5 will break the surface. However, there are events with much smaller magnitudes rupturing the ground. To understand the potential controlling mechanisms, we first collect source parameters for earthquakes with  and  surface-breaching events in seismically active regions including west China, North America, Europe, Taiwan, Japan, and Iran. For strike-slip and normal events, almost all earthquakes with magnitudes over 6.7 broke the surface. In contrast, buried and surface-breaching events co-exist with moderate magnitude (6.0-6.7). For reverse events, there is no clear magnitude boundary, as thrust buried events can be quite large due to the downdip size of the seismogenic zone. The relocated hypocenter depths for moderate-to-large events are concentrated at depth of 5 -20 km with no systematic difference between buried and surface-breaching ruptures. Differently, all  surface-breaching events occurred at very shallow depths (<5 km). We also conduct dynamic rupture simulations and propose two conceptual models to explain whether or not ruptures may break the surface. The first model represents a fault with a continuous but heterogeneous seismogenic zone (velocity-weakening) that can hold moderate-to-large earthquakes. In this case, ruptures need to overcome the shallow velocity-strengthening zone (VS) with certain energy sink to reach the surface. Therefore, a thinner shallow RS zone and a higher stress drop of the earthquake can promote surface rupture, consistent with our observations. However, ruptures nucleating from different locations on heterogeneous faults may lead to different surface rupture patterns and final magnitudes, shedding lights on the diverse behaviors among moderate earthquakes. The second model is for small surface-breaching earthquakes. Those events are supposed to occur on shallow isolated velocity-weakening patches, consistent with the fact that usually no large earthquakes have been reported on the same fault zones. Such asperities may be formed on bodies with high-strength materials, leading to energetic ruptures with intense stress release. Our study contributes to the understanding of the surface rupture behaviors references for assessing near-surface damage in future earthquakes.

How to cite: Yao, S., Yang, H., and Tang, Z.: Investigating relationships between surface rupture and multiple source parameters of earthquakes, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7142, https://doi.org/10.5194/egusphere-egu23-7142, 2023.