How to cite: Bletery, Q. and Nocquet, J.-M.: Do earthquakes start with precursory slow aseismic slip?, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9961, https://doi.org/10.5194/egusphere-egu24-9961, 2024.

How to cite: Ju, F., Meng, H., Chen, X., and Yu, C.: Detection of Immediate Foreshocks Using Dense Seismic Array: A Case Study of the 2021 Ms 6.4 Yangbi aftershock sequence, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3922, https://doi.org/10.5194/egusphere-egu24-3922, 2024.

The acoustic emissions (AEs) produced during the shearing of granular materials reflect the accumulation and release of stress, offering valuable insights into the failure mechanisms of seismic faults and stick-slip-controlled landslides. While various characteristics such as amplitude, energy, counts, and frequency of AE signals generated by stick-slip have been studied, the stress changes corresponding to different frequency AEs at various stages of the stick-slip process remain unclear. This knowledge gap hinders our understanding of the precursory signals leading to stick-slip failure. In order to enhance our comprehension of the physical mechanisms underlying granular stick-slip, we conducted monitoring of both mechanical and AE signals using high-frequency (2 MHz) synchronous acquisition. This was done during the constant-speed shearing of packs containing uniform glass beads of different sizes under varying normal stresses. Our findings revealed an accelerated release rate of AE energy in tandem with sample volume dilatation. Additionally, the stress drop during stick-slip increased with higher normal stress and particle size. This study identified three distinct events during a single cycle of stick-slip: main slip, minor slip, and microslip. We analyzed the AE frequency spectra associated with each of these events. Main slip and minor slip correlated with stress drop, generating high-frequency AEs (approximately several hundred kHz). In contrast, microslip produced lower AE frequencies (around tens of kHz) and exhibited stress strengthening. These characteristics, overlooked in prior studies due to low-frequency acquisition, suggest that microslip is primarily a result of sliding on grain contacts, while main slip and minor slip arise from the breakage and reformation of force chains. The low-frequency AEs from microslip may serve as a crucial precursor to seismic faults and landslides, providing a deeper understanding of the granular stick-slip phenomenon.

How to cite: Gou, H. and Hu, W.: Detection of Stick-Slip Nucleation and Failure in Homogeneous Glass Beads Using Acoustic Emissions in Ring-Shear Experiments: Implications for Recognizing Acoustic Signals of Earthquake Foreshocks, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5252, https://doi.org/10.5194/egusphere-egu24-5252, 2024.

A growing number of observations made using geodetic approaches have been able to detect large preparatory regions that experience accelerated deformation prior to and in close proximity to an earthquake’s hypocenter. An uptick in localized seismicity has also been observed in these regions and represents an opposite end-member of the spectra of deformation, in both space and time, from the opposite broad and slow process. If and how these preparatory observations are linked are not well understood. To study this, we conducted a triaxial experiment on a granitic rock sample instrumented with calibrated acoustic emission (AE) sensors and a distributed strain sensing (DSS) method using fibre optics. These two technologies were sensitive to seismic (100 kHz to 1 MHz) and aseismic (DC to 0.4 Hz) deformation at our sample scale and these were monitored as it was loaded and experienced brittle shear failure. DSS measurement allowed us to visualize the emergence of slow, heterogeneous strain fields that localized well before the failure of the sample. In the early stages of localized deformation, the regions exhibiting preferential damage were growing and doing so without producing seismicity. However, when approaching failure, these regions accommodating slow deformation began to accelerate and now produced clusters of seismicity. The cumulative seismic moment of the precursory seismicity was a fraction of the total anelastic deformation (< 0.1%) precluding the runaway dynamic failure. We also examined the clustering and frequency-magnitude distribution of the seismicity with respect to the localized strain field. In the later stages, moments prior to nucleation, the b-value begins to drop and becomes anti-correlated to the rapidly accelerating average volumetric strain rate measured using the DSS array. This observation better constrains the hypothesis that dilation of the relatively large preparation zone can host larger precursory earthquakes therein. These findings can help constrain models that better replicate the physics associated with the large spectrum of brittle deformation and will in turn help with our understanding of preparatory earthquake processes.

How to cite: Selvadurai, P. A., Salazar Vasquez, A. F., Bianchi, P., Madonna, C., Germanovich, L., Puzrin, A. M., Rabaiotti, C., and Wiemer, S.: Localizing slow deformation holds crucial information related to seismicity patterns precluding brittle failure in crystalline rocks, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10887, https://doi.org/10.5194/egusphere-egu24-10887, 2024.

Faults in the brittle crust lie at any orientation to the far-field stress. However, laboratory experiments designed to investigate earthquake physics commonly simulate favorably oriented faults, potentially overlooking the complexity of natural fault behavior. Here, we assess the role of stress field orientation in fault reactivation and earthquake precursors by conducting triaxial saw-cut experiments with laboratory faults oriented at different angles to the maximum principal stress, ranging from 30° to 70°. The samples were instrumented with strain gauges and piezo-electric sensors. Laboratory well-oriented faults describe a rather simple system in which the elastic energy is stored via the deformation of the surrounding host rock during the inter-seismic period and released via on-fault slip during the co-seismic phase with associated precursor acoustic activity. Consistent with previous laboratory data, an abrupt increase in the on-fault acoustic emission rate occurs shortly before the laboratory earthquake. A more complex picture emerges when deforming laboratory misoriented faults. Particularly, acoustic emissions and strain gauge data indicate that when the fault is misoriented, off-fault permanent deformation occurs well before fault reactivation. The stress state in the host rock surrounding the fault is indeed far beyond the one required for the onset of inelastic deformation. In this case, acoustic activity distributed in the rock volume during the pre-seismic phase is associated with permanent deformation in the critically stressed host rock and is not a direct precursor to the following laboratory earthquake. Unlike well-oriented faults, laboratory mis-oriented faults lack detectable seismic precursors. The laboratory-observed increase in acoustic activity prior to, but not precursor of, mis-oriented fault reactivation impacts our understanding of earthquake precursors in natural faults.

How to cite: Giorgetti, C. and Brantut, N.: Fault orientation in earthquake seismic precursors: insights from the laboratory, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15878, https://doi.org/10.5194/egusphere-egu24-15878, 2024.

The Mw 7.8 earthquake of 2023 ruptured the southern (main) branch of the East Anatolian Fault (EAF), followed by the Mw 7.6 earthquake on the northern (Çardak Fault) branch of the EAF nine hours later. In March, we installed ten carbon-rod extensometers across segments of the main and northern branches of the ruptured faults, where potential slip deficits were considered possible, to investigate if afterslip continues. It is important to measure afterslip to understand the behaviour of a fault, if any, resulting from stresses associated with local coseismic slip deficits.  Seven of these extensometers recorded less than a few millimetres of slip since March. In Göksun near the western end of the Çardak fault, we recorded more than a 25 mm of accelerating afterslip preceding local aftershocks of magnitude ≤Mw5.1.

The Mw 6.8 Elazığ-Sivrice earthquake of January 24, 2020, and the Mw 7.8 Kahramanmaraş earthquake of February 6, 2023, stopped in the Pütürge region, where it is named as Pütürge Gap. To understand why these two earthquakes terminated there, an array of five extensometers were ultimately deployed. One of the extensometers, which is 52 m long, shows that slip > 3.8 mm/yr continues at depth. Extensometers spaced 45 km apart recorded an eastward propagating subsurface creep event in September 2023. Four cGPS stations recording at 5 Hz were installed in an array to better investigate the subsurface evolution of aseismic slip on the Pütürge Fault in the village of Taşmış.

How to cite: Ayruk, E. T., Turğut, M., Farımaz, İ., Köküm, M., Bilham, R., and Doğan, U.: AFTERSLIP of 6 FEBRUARY 2023 KAHRAMANMARAS EARTHQUAKE SEQUENCE : PRELIMINARY RESULTS, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-428, https://doi.org/10.5194/egusphere-egu24-428, 2024.

The interactions between aseismic slip and seismicity is mostly studied during postseismic afterslip associated with the generation of aftershocks following large earthquakes. However because of the large stress perturbations produced by the coseismic rupture, it remains difficult to distinguish the contribution of the coseismic stress perturbation and the effects of stresses induced by the afterslip on the triggering of aftershocks. Studying seismicity triggered by slow slip events enables us to understand the direct effect of aseismic slip on the generation of the seismicity. While most subduction slow slip events are deep-seated, at the down-dip edge of seismogenic zones accompanied by tectonic tremors, some are also observed at shallower depths associated with seismic swarms. Among them are the well-studied Boso slow slip events located on the Sagami trough, between 10 and 20 km depth. Recorded every ~4 years since 1996, they are always accompanied by swarms of Mw 1 to 5 earthquakes on their northern and western flanks. Being located right beneath the Boso Peninsula coastline, the kinematics of these slow slip events is particularly well imaged by dense GNSS and tiltmeter networks. In this study, we model the time dependent aseismic slip of the 2018 Boso slow slip event, with the largest moment released of all Boso slow slip events, by inverting time step by time step the slip on the fault with joint GNSS and tiltmeter data. We do not impose arbitrary temporal smoothing in the inversion and find that the well constrained fault slip is first growing and then migrating southwestward with a migration speed of ~ 2 km/day. In order to model the interaction with the seismicity, we compute the Coulomb stress change due to the transient slip on receiver faults located in 9 cubes centered in the seismicity swarm and parallel to the subduction interface. From June 2nd to June 18th, the seismicity is migrating up-dip at a rate of 1 km/day. This migration period coincides with the peak slip rate and with Coulomb stress produced by the slow slip migrating updip together with the seismicity, indicating a causal relationship. Adopting a rate and state friction formalism to explain the nucleation of the seismicity, we finally investigate the ensemble of parameters, in particular the constitutive parameter that relates changes in stress to logarithmic changes of slip velocity, the effective normal stress and the tectonic stressing rates, that can explain the seismicity rate. 

How to cite: Rousset, B., Inbal, A., Bürgmann, R., Uchida​​, N., Socquet, A., Marill, L., Matsuzawa, T., and Kimura, T.: Boso seismicity swarm propagation driven by slow slip stress change, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4561, https://doi.org/10.5194/egusphere-egu24-4561, 2024.

In earthquake research, the discovery and ongoing investigation of interseismic transient processes has revealed that faults are non-steady between large earthquakes. These transients are typically identified in continuously operating tectonic GNSS stations, whereby an acceleration away from the average interseismic rates of displacement can be identified with a variety of time series analysis methods. However, the features of these transients can vary depending on the processing strategy employed to derive displacement time series from the raw GNSS observables.

In the processing strategy, the definition of the geodetic datum is necessary to determine global terrestrial reference frames (TRFs), providing an accurate and stable absolute reference of Earth's locations. It is essential for comprehending the dynamic changes in Earth's geometry driven by factors like tidal and non-tidal loading, plate tectonic seismic activity, and ongoing climate change. Therefore, just as geodesy aims for accuracy and stability in the TRF, the datum definition—i.e., the realization of the TRF-defining parameters origin, orientation, and scale—may emerge as a critical factor in processing GNSS networks for geodynamic purposes.

The purpose of this study is to assess up to what extent the transient velocities obtained from GNSS-derived displacement time series change under different regional and global datum definitions for the Cascadia subduction zone and Hikurangi margin; regions with very well-known catalog of interseismic transient tectonic events. In our study, we process data from Cascadia to produce network solutions both NNR (No-Net-Rotation) and NNR+NNT (No-Net-Translation) constraint for regional and global datum definition, respectively. We employed dual-frequency ionosphere-free linear combination observations from 125 GNSS stations for the time between 2015 and 2020. The same GNSS processing strategy is then followed for the Hikurangi subduction zone using a set of 72 stations from the GeoNet project as well as the same control stations used for Cascadia spanning from 2002 to 2010.

For the Cascadia displacement time series, we find variations in transient velocities under different datum definitions emphasizing the need for a comprehensive understanding of its impact on dynamic geophysical processes. Processing and analyses of the New Zealand data is ongoing and results will be presented, along with recommendations for both regions on how to reduce the occurrence of likely non-tectonic transients in the displacement time series. Ultimately, our results may have implications for improving the estimate of the slip budget at plate boundaries that is released aseismically.

How to cite: Garcia, C., Bedford, J., Männel, B., and Glaser, S.: Impact of datum definition on transient velocities from GNSS displacement time series in Networks mode: A Case Study of Cascadia Subduction Zone and Hikurangi margin., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12979, https://doi.org/10.5194/egusphere-egu24-12979, 2024.

Topographic features such as seamounts can influence the buoyancy of the slab and the short- and long-timescale mechanical properties of the subduction interface. How seamounts in the trench interact with the upper plate accretionary wedge during subduction— their stress field, their potential for ‘decapitation’, and their ability to host large megathrust earthquakes— is not fully understood. We utilize exhumed rocks to investigate seamount–upper plate interactions at shallow subduction interface conditions. 

We focus on a 250-m-thick cross-section of deformed, weakly metamorphosed basalt, limestone, chert, argillite and greywacke exposed in the inboard part of the Chugach accretionary complex near Grewingk glacier, southern Alaska. Temperatures from Raman spectroscopy on graphite yield ~260°C, suggesting deformation and metamorphism down to ~15-20 km depth. Detrital zircon data from greywacke lenses within and outside the shear zone overlap within error suggesting emplacement over less than ~1 m.y. at ~167 Ma. 

Basalts in the shear zone are dismembered into ~3 slices up to 35 m thick, all of which contain limestone patches suggesting the basalt is derived from the seamount’s very top (limited decapitation). The basalt slices are bounded by high-strain melange-like shear zones up to 25 m thick, interpreted to represent décollements along which the seamount slices were underplated. These mélange belts exhibit a block-and-matrix texture with a macroscopically ductile argillite and chert matrix, and pervasively disaggregated and brittlely deformed greywacke and basalt lenses. Both the matrix and the blocks show several generations of dilational and shear veins, suggesting high fluid pressures and low differential stresses. Features suggesting deformation at fast (potentially seismic) strain rates include fluidized cataclasites, but these do not extend along strike for more than 0.25 m and do not occur within the larger (m-to-dm-scale) basalt lenses, suggesting that large-magnitude earthquakes were limited during seamount underplating. Instead, the observed mix of brittle and macroscopically ductile deformation at high fluid pressures is more consistent with a potential record of shallow tremor and slow slip.

Our findings support geophysical observations and numerical models that suggest relatively weak mechanical and seismic coupling between seamounts and the overriding plate, and are consistent with recent suggestions (e.g. for the Hikurangi margin) that sediment envelopes around subducting seamounts are conducive to slow slip and tremor.

How to cite: Behr, W., Akker, I. V., and Rast, M.: Deformation processes during seamount dismemberment and underplating along the shallow subduction interface: a case study from the Chugach Complex, Alaska, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8244, https://doi.org/10.5194/egusphere-egu24-8244, 2024.

A general correlation between maximum co-seismic fault slip and fault length has been well constrained by observations from numerous earthquakes occurring on many crustal faults. In spite of this global consistency, there are notable examples of co-seismic fault slip magnitudes that far exceed the expected maxima for the fault dimensions.  Two instances of extreme fault slip that occurred on short-to-moderate length upper-plate faults in subduction systems, where the megathrust lies at shallow depths are the 2016 Kaikoura (New Zealand) earthquake, and the 1855 Wairarapa (New Zealand) earthquake. In both cases, co-seismic fault slip was 5 to 10 times greater than expected from the rupture length - fault slip scaling relationships. The typical crustal fault model  is a fault with a brittle-to-ductile transition at depth; in this scenario, co-seismic slip is inhibited by viscous resistance from the deeper, ductile component of the fault.  In contrast, the tectonic characteristics of faults that experience extreme co-seismic slip, involve upper plate faults that truncate against the megathrust at seismogenic depths (i.e. are fully frictionally coupled over their entire depth extent). During a megathrust earthquake, the plate interface unlocks and upper-plate faults that extend to the ruptured plate interface transiently experience free-slip boundary conditions on both their upper (surface) and lower (megathrust) ends. As a result, such upper-plate faults can potentially experience full strain release (and therefore maximum slip), independent of their length. For appropriately oriented faults, this effect may be enhanced by co-seismic stress changes associated with the megathrust earthquake. 

Geologic evidence of large displacement (and/or displacement rate) upper-plate faults in other subduction systems indicates this process may commonly occur. One example is a set of upper-plate faults along the Cascadia margin (near Newport, Oregon), that have strike-slip geologic slip rates, averaged over 10s of thousands of years, exceeding tens of mm/yr and approaching local plate convergence rates. These upper-plate Cascadia faults are also located where the plate interface is sufficiently shallow and seismogenic, indicating that these high-slip, upper-plate faults are likely frictionally locked over their entire depth range. In spite of the high overall slip rates of these upper-plate faults, because they are locked along their entire depth extent between earthquakes,  they may be unrecognized by inter-seismic geodetic observations.

How to cite: Furlong, K. P., Herman, M. W., and McKenzie, K. A.: Big Slip on Small Faults - How Does Extreme Fault Slip Occur on Short-to-Moderate Length Faults, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4427, https://doi.org/10.5194/egusphere-egu24-4427, 2024.

The shallowest region of subduction megathrust accommodates deformation by a spectrum of seismic modes including continuous aseismic creep and peculiar seismic phenomena as slow slip events. However, the mechanisms behind these phenomena remain enigmatic because they are not explained by conventional frictional models. This because the shallowest regions of subduction zones are characterized by unconsolidated, clay-rich lithologies that, nominally, cannot nucleate seismic events due to their frictionally weak and rate-strengthening attributes. Here we present laboratory friction experiments showing that clay-rich experimental faults with bulk rate strengthening behavior and low healing rate can contemporaneously creep and nucleate slow slip events. These instabilities are self-healing, slow ruptures propagating within a thin shear zone and driven by structural and stress heterogeneities. We propose that the bulk rate-strengthening frictional behavior promotes the observed long-term aseismic creep whereas local frictional mechanism causes slow rupture nucleation and propagation. Our results illustrate the complex behavior of clay-rich lithologies, providing a new paradigm for the interpretation of the genesis of slow slip as well as significant implications for seismic hazard.

How to cite: Volpe, G., Collettini, C., Taddeucci, J., Marone, C., and Pozzi, G.: Complex Frictional Behavior of Clay and Implications for Slow Slip , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-5061, https://doi.org/10.5194/egusphere-egu24-5061, 2024.

Aseismic slip has been recognized over the last 50 years as one of the modes of elastic stress accommodation by large tectonic faults. Long strike slip faults have been imaged with InSAR and scrutinized with GNSS networks and creepmeters, revealing the strikingly ubiquitous occurrence of aseismic slip globally. From the 600-km-long creeping section of the Chaman to the 70 km-long Ismetpasa creeping section along the North Anatolian Fault, geodetic imaging illustrates the rich behavior of such aseismic slip, from mm-scale transient episodes of slip to seemingly continuously sliding fault segments.

Most models explaining the occurrence of aseismic slip along continental faults rely entirely on an ad hoc parameterization of the frictional rheology of the fault. While the friction law governing slip along faults has been determined from laboratory experiments, the inference of the constitutive parameters of such friction law entirely derives from reproducing geodetic data in most cases. In particular, most creeping sections are interpreted as the signature of rate-strengthening material, diffusing stress through stable sliding. However, most rocks at seismogenic depths exhibit a rate-weakening behavior and some even show transient episodes of slip incompatible with purely strengthening properties. In addition, other mechanisms, including complex geometric configuration of faults or fluid circulation may offer the conditions for slow slip. Therefore, the direct inference of constitutive properties of a fault zone from kinematic observations may not be simple.

We propose here a model in which upwelling of fluids sourced in the upper mantle through a vertical fault zone leads to the conditions for slow slip, irrespective of the fault constitutive properties. We map aseismic slip along three different fault zones, including the North Anatolian Fault (Turkiye), the San Andreas Fault (USA) and the Leyte fault (Philippines) and find a systematic relationship between the effective locking depth and the occurrence of aseismic slip. Our model explains this modulation of locking depth along strike and the subsequent modulation of surface shear stressing rate with the along strike variation in the mantle fluid source. This model applied to fault segments with relatively high mantle fluid source leads to low effective normal stress, large nucleation size of a frictional instability, and predicts occurrences of shallow aseismic slip. We perform numerical modeling to show that the critical parameter is the flux of upwelling fluid through the fault zone, which increase leads to the widening of the near-surface region of aseismic slip and transition to full-fault aseismic slip at large enough flux. We finally discuss the potential sources of fluids explaining such behavior.

How to cite: Jolivet, R., Garagash, D., Pierre, D., and Jara, J.: Creeping sections on continental strike slip faults as the signature of deep fluid upwelling, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9437, https://doi.org/10.5194/egusphere-egu24-9437, 2024.

How to cite: Thomas, M. Y. and S. Bhat, H.: Combined Effect of Brittle Off-Fault Damage and Fault Roughness on Earthquake Rupture Dynamics, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-12095, https://doi.org/10.5194/egusphere-egu24-12095, 2024.

Deformation within the upper crust is mainly accommodated through slip on fault systems. These systems can accommodate slip via different modes, going from aseismic creep (i.e., stable motion) to dynamic earthquake (i.e., unstable motion). Notably, a single fault is not confined to a specific slip mode, as recent geodetic observations have indicated that a single fault can exhibit both stable and unstable motions. The distinct slip behaviours have been attributed to fault spatial heterogeneity of the frictional properties, rheological transitions, or geometrical fault complexities.

To comprehensively characterize the impact of frictional heterogeneities, we deformed heterogeneous fault samples in a triaxial apparatus, at confining pressure ranging from 30 to 90 MPa. The fault planes, sawcut at a 30° angle from the sample axis, consisted of two materials: granite and marble. Experiments were conducted for both marble asperities embedded in granite and vice versa, alongside homogeneous fault samples of single lithology. The selection of granite and marble was based on their different mechanical and frictional characteristics, with granite exhibiting seismic behaviour, while marble demonstrated aseismic behaviour across the pressure range tested.

Our findings reveal that the stress drops of seismic events are dependent on fault composition, with faults containing higher granite content exhibiting larger seismic events. In addition, by coupling the inversion of the kinematic slip from strain-gauge measurements and the records of acoustic activity during experiments, we demonstrate that the nucleation and propagation of seismic events are significantly influenced by lithological heterogeneity on the fault plane. In the case of homogeneous faults, the seismic event nucleation is relatively straightforward, initiating in the highest stressed region and propagating uniformly. Conversely, heterogeneous faults display more intricate nucleation patterns, often featuring multiple nucleation regions converging into a major dynamic event. The dynamic event propagation is expedited when traversing granite areas and more restrained within the marble. Remarkably, our experiments demonstrate that heterogeneities are required in order to induce earthquake afterslip. These results emphasize the crucial role of fault heterogeneity in earthquake nucleation and propagation, highlighting that even minor lithological heterogeneities are sufficient to complicate laboratory earthquake dynamics.

How to cite: Noël, C., Dublanchet, P., and Passelègue, F.: Exploring the impact of frictional heterogeneities on the seismic cycle: Insights from laboratory experiments, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9595, https://doi.org/10.5194/egusphere-egu24-9595, 2024.

Many previous studies have explored the role of granular media in controlling friction of faults. A gap exists though in understanding the failure process and sliding of a fluid-saturated fault gouge. Here we use a coupled 2D DEM-fluid code to simulate fault-gouge as a layer of grains, sheared by a constant stress boundary. We explore and compare two scenarios: 1) a dry granular layer, in which shear stress on the top wall is incrementally increased, or 2) a fluid-saturated granular layer, into which fluid is injected, so that fluid pressure is incrementally increased. Once the applied stress/pressure is high enough, the layer fails and starts accelerating, until it reaches a steady-state sliding rate (determined by the layers’ velocity-strengthening friction). We next incrementally step-down the shear stress or fluid pressure. Consequently, the slip-rate is observed to slow down linearly with decreasing stress/pore-pressure, until the layer finally stops, at a stress/pressure lower than that required to initiate the failure. Both the dry and fluid-saturated granular systems exhibit two main behaviors: 1) velocity-strengthening friction, following the mu(I) rheology, 2) a hysteresis effect between friction and velocity, porosity and grain coordination numbers. The hysteresis and strain-rate dependence agree with previous experimental, numerical and theoretical results in dry granular media, yet our work suggests these behaviors extend to fluid-filled granular media. We theoretically predict the transient and steady-state observations for dry and fluid-saturated layers, using the mu(I) friction rheology with an added component of hysteresis. Importantly, we show that fluid-filled faults exhibit a process which is absent in dry systems: fluid-injected layers may exhibit failure delay, with some time passing between pressure rise and failure. We link this delay to pre-failure creeping dilative strain, interspersed by small dilative slip events. Our numerical and analytical results may explain: (i) field measurements of fault creep triggered by fluid pressure rise (e.g. via injection), (ii) fault motion which is triggered by fluid-injection but continues even after fluid pressure returns to its pre-injection level. (iii) observed delay prior to failure in fluid-injection experiments.

How to cite: Aharonov, E., Sarma, P., Toussaint, R., and Parez, S.: Fluid-induced failure and sliding of a gouge-filled fault zone: Hysteresis, creep, delay and shear-strengthening., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-10783, https://doi.org/10.5194/egusphere-egu24-10783, 2024.

In recent decades, there has been a proliferation of observations related to spatiotemporally intricate slip events occurring in fault systems. These events encompass a spectrum of transient energy releases, ranging from slow slip events to low-frequency earthquakes (LFEs) and tremors, in addition to the more familiar creep and fast ruptures. The prevailing focus in recent research has been to interpret these events by considering variations in frictional behavior along the fault plane.

However, it is crucial to acknowledge the inherent geometric complexity of fault systems across multiple scales. Recent studies have illuminated the significance of incorporating a fault volume or damage zone surrounding the fault in the analysis of slip dynamics. In the context of this study, we endeavor to investigate the influence of "realistic" fault geometry on the dynamics of slip events. To achieve this, we approach the problem from three interrelated perspectives:

• Forward Source Modeling: We employ forward source modeling techniques to simulate and understand the behavior of slip events.
• Bridging Source Modeling and Observations: We establish a connection between our source modeling and observed data by generating synthetic surface records that can be compared to actual observations.
• Energy Budget Analysis: We meticulously analyze the variations in the energy budget that occur throughout the earthquake cycles to gain insights into the mechanics of slip events.

Our primary objectives include deciphering how deformation within the volume is accommodated by both the off-fault damage zone and the primary fault. Specifically, we aim to determine the proportion of the supplied moment rate that is absorbed by off-fault fractures during an earthquake cycle. Additionally, we seek to unravel how the diverse sequences of complex behavior observed on the fault plane manifest in the signals recorded by seismic stations. This entails assessing the distinct contributions of the main fault and off-fault fractures to the radiated signals detected at the monitoring stations. Lastly, we delve into the evolution of the medium's energy budget throughout the earthquake cycles and evaluate the dissipative contribution of off-fault fractures to ascertain their energetic role in the context of earthquake cycles.

How to cite: Bhat, H., Almakari, M., Kheirdast, N., Villafuerte, C., and Thomas, M.: A mechanical insight into the continuous chatter of a fault volume, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-9007, https://doi.org/10.5194/egusphere-egu24-9007, 2024.