pith. sign in

arxiv: 2509.04909 · v3 · submitted 2025-09-05 · ⚛️ physics.geo-ph

Fault volume digital twin to reproduce the full slip spectrum, scaling and statistical laws

M. Almakari , N. Kheirdast , C. Villafuerte , M. Y. Thomas , P. Dubernet , J. Cheng , A. Gupta , P. Romanet
show 2 more authors
S. Chaillat H. S. Bhat
This is my paper

Pith reviewed 2026-05-18 19:07 UTC · model grok-4.3

classification ⚛️ physics.geo-ph
keywords fault damage zoneslip spectrumOmori lawGutenberg-Richter scalingrate-and-state frictiondigital twinquasi-dynamic simulationsoff-fault fractures
0
0 comments X

The pith

A main fault surrounded by power-law distributed off-fault fractures reproduces the full range of slip speeds and seismic laws with uniform friction.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper establishes that modeling a fault zone with an off-fault damage zone of distributed fractures can explain diverse slip behaviors and statistical patterns seen in real earthquakes and geodesy. The fractures follow power-law size and density distributions, are oriented optimally or parallel to the main fault, and all obey the same rate-and-state friction law without any added spatial variations in properties. A sympathetic reader would care because this single setup generates a natural continuum from slow to fast ruptures along with Omori laws, Gutenberg-Richter scaling, and moment-duration relations. It also produces localization of seismicity before main events and migration patterns during slow slip that resemble diffusion fronts, all emerging from the fault volume itself.

Core claim

Incorporating an off-fault damage zone characterized by distributed fractures surrounding a main fault can reproduce many key features observed in seismic and geodetic data, including a natural continuum from slow to fast ruptures, Omori law, inverse Omori law, Gutenberg-Richter scaling, and moment-duration scaling, without introducing spatial heterogeneities in frictional properties. All fractures follow rate-and-state friction with parameters enabling slip instabilities, and quasi-dynamic boundary integral simulations show tremors, very low frequency earthquakes, low frequency earthquakes, slow slip events, and earthquakes emerging naturally within this framework.

What carries the argument

off-fault damage zone modeled as power-law distributed fractures surrounding the main fault, all obeying uniform rate-and-state friction

If this is right

  • A natural continuum from slow to fast ruptures emerges despite spatially uniform frictional properties.
  • Omori law, inverse Omori law, Gutenberg-Richter scaling, and moment-duration scaling are all reproduced by the same setup.
  • Seismicity localizes toward the main fault before nucleation of main-fault events.
  • During slow slip events, off-fault seismicity migrates in patterns resembling fluid diffusion fronts without any fluids present.
  • Tremors, very low frequency earthquakes, low frequency earthquakes, slow slip events, and regular earthquakes all arise within the single fault-volume model.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The framework could serve as a controlled testbed for inverting mechanical properties that are hard to measure directly in the field.
  • Similar damage-zone effects might unify slow and fast slip behaviors in other tectonic settings beyond the 2D shear fault examined here.
  • Extending the model to three dimensions could reveal whether the same statistical laws persist when fracture orientations and interactions become more complex.

Load-bearing premise

Off-fault cracks must follow power-law size and density distributions, be oriented optimally or parallel to the main fault, and obey rate-and-state friction parameters that enable slip instabilities.

What would settle it

A quasi-dynamic simulation of an isolated main fault with the same uniform friction but no surrounding power-law damage zone would fail to produce the continuum of slip rates or the full set of statistical laws.

Figures

Figures reproduced from arXiv: 2509.04909 by A. Gupta, C. Villafuerte, H. S. Bhat, J. Cheng, M. Almakari, M. Y. Thomas, N. Kheirdast, P. Dubernet, P. Romanet, S. Chaillat.

Figure 1
Figure 1. Figure 1: Hierarchical structure of fault systems over a wide range of length scales. (a) Fault map of Southern California (Fletcher et al. 2014). Black lines indicate fault traces. Stars and col￾ored lines indicate the epicenters and rupture traces of historical earthquakes. (b) Fault map and rupture traces (in red) associated with the 1992 Landers earthquake (modified from Sowers et al. (1994)). (c) Smaller-scale … view at source ↗
Figure 2
Figure 2. Figure 2: Fault volume geometry of a case study: a) sketch of fault volume geometry (not to scale): the main fault in black and the off-fault fractures in red b) off-fault fracture density per￾pendicular to the main fault (linear scale) c) probability density function of length distribution of the faults, including the main fault d) rose diagram of off-fault fracture orientations and principal stress direction e) cr… view at source ↗
Figure 3
Figure 3. Figure 3: Time series of the moment rate M˙ 0. The black curve represents the contribution com￾ing from the main rough fault, and the red curve represents the summation of the moment rate released by all the fractures. Panel a) The entire seismic cycle. Top x-axis is absolute time in years, and bottom x-axis is time normalized by recurrence time of earthquakes when considering only the main rough fault, without a da… view at source ↗
Figure 4
Figure 4. Figure 4: Sequences of slip rate profiles on the main fault. Top x-axis is average slip accumulated during the specific time sequence. Left y-axis represents the along-strike distance on the main fault. Horizontal lines show the specific durations for different phases of the rupture sequence, with blue lines indicating slow phases and red lines indicating fast phases. Nucleation and after￾slip phases are indicated b… view at source ↗
Figure 5
Figure 5. Figure 5: (a) Omori Law: Seismicity rate decreases over time following the mainshock. Each gray curve represents one earthquake cycle with one mainshock on the main fault, the blue curve represents the stacked sequences of all the seismic cycles in this simulation case (b) Magnitude￾frequency distribution of the cataloged fast ruptures follows the Gutenberg-Richter distribution. Red and blue colors represent fast an… view at source ↗
Figure 6
Figure 6. Figure 6: Scaling laws of inferred seismicity: Moment-duration scaling with a threshold of detec￾tion of slow ruptures at (a) 10−6 m/s and (b) 10−8 m/s. Events are color-coded based on rupture velocity vr . In Figure 6a, we delve into the moment-duration scaling, a crucial aspect of seismic rup￾ture behavior. Observations in natural seismicity have highlighted a cubic relation be￾tween moment and duration for fast r… view at source ↗
Figure 7
Figure 7. Figure 7: Effect of slip rate threshold (Vth), or signal to noise ratio (SNR), on the inference of the magnitude and duration of a slow slip event. The slip rate is plotted with four different minimum detection levels: (a) 10−9 m/s, (b) 10−8 m/s, (c) 10−7 m/s, and (d) 10−6 m/s. (e) Same as (a) with off-fault seismicity. Circles represent ruptures on the off-fault fractures, and are color-coded with respect to their … view at source ↗
Figure 8
Figure 8. Figure 8: Localization and delocalization of off-fault seismic activity, measured by the standard deviation of hypocentral distances from the main fault in sequential 500-event batches, λy. The timing of earthquakes is indicated by yellow stars. The standard deviation values are normalized by their maximum value over all batches. The gray region highlights the time window shown in [PITH_FULL_IMAGE:figures/full_fig_… view at source ↗
Figure 9
Figure 9. Figure 9: Zoomed-in view of the gray-hatched window shown in [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Off-fault Seismicity during two examples of slow-slip ruptures. a-c) Slip rate profiles of two slow slip events. b-d) The temporal migration of off-fault seismicity follows a distance from epicenter (r)–time (t) evolution (black dashed curve) similar to that of a diffusion front (r ∝ D √ t). The relative distance is calculated as the distance between the off-fault event and the detected nucleation point o… view at source ↗
Figure 11
Figure 11. Figure 11: Moment-rate functions (MRFs) (left panels), and their spectra (right panels), across different Mw ranges for both off-fault and main fault events. (a) Selected MRFs from slow and fast events in the off-fault region within a narrow range of −0.9 ≤ Mw ≤ −0.8. (b) Same as (a), but for −0.3 ≤ Mw ≤ −0.2. (c) MRFs from main fault events, and their spectra, with a broader range of 2.2 ≤ Mw ≤ 3.8 and 2.4 ≤ Mw ≤ 4… view at source ↗
Figure 12
Figure 12. Figure 12: MRF shapes for slow and fast events after normalization by centroid time. Panels (a) and (b) correspond to slow events in the off-fault region for the same magnitude ranges as in Figures 11a and 11b. (c) shows slow events on the main fault. Panels (d) and (e) correspond to fast off-fault events with the same Mw ranges as in (a) and (b). (f) Fast events on the main fault. Color-coded by Mw. events across a… view at source ↗
Figure 13
Figure 13. Figure 13: (a) Effect of off-fault fracture orientation on the Omori law: The average Omori decay curve over a period of 30 days is calculated from the stacked curves of different fracture orien￾tation, and the average p coefficient is reported (b) Effect of off-fault fracture orientation on the b-value found from stacked catalog of fast ruptures with same off-fault fracture orientation (c) Effect of off-fault fract… view at source ↗
Figure 14
Figure 14. Figure 14: Moment-duration scaling for the comprehensive catalog compiled from all simula￾tions : a) with a threshold of detection of slow ruptures at 10−6 m/s, showing a M0 ∝ T 3 relation for fast and slow ruptures; b) with a threshold of detection of slow ruptures at 10−9 m/s. Events are color-coded based on rupture velocity vr . events with confidence reinforces the presence of this gap, raising the possibility o… view at source ↗
Figure 15
Figure 15. Figure 15: Summary of the digital twin model, SAFAR, that reproduces broadband slip dynamics, statistics and scaling laws of a fault volume. frequency distribution, Omori and inverse-Omori behavior of seismicity rates before and after mainshocks, moment–duration and moment–area scaling relationships, and the ob￾served localization-delocalization of seismicity around an event. We also observe that off￾fault seismicit… view at source ↗
read the original abstract

Seismological and geodetic observations of fault zones reveal diverse slip dynamics, scaling, and statistical laws. Existing mechanisms explain some but not all of these behaviors. We show that incorporating an off-fault damage zone-characterized by distributed fractures surrounding a main fault-can reproduce many key features observed in seismic and geodetic data. We model a 2D shear fault zone in which off-fault cracks follow power-law size and density distributions, and are oriented either optimally or parallel to the main fault. All fractures follow rate-and-state friction with parameters enabling slip instabilities. We do not introduce spatial heterogeneities in frictional properties. Using quasi-dynamic boundary integral simulations accelerated by hierarchical matrices, we simulate slip dynamics and analyze events produced both on and off the main fault. Despite spatially uniform frictional properties, we observe a natural continuum from slow to fast ruptures, as seen in nature. Our simulations reproduce the Omori law, inverse Omori law, Gutenberg-Richter scaling, and moment-duration scaling. We observe seismicity localizing toward the main fault before nucleation of main-fault events. During slow slip events, off-fault seismicity migrates in patterns resembling fluid diffusion fronts, despite the absence of fluids. We show that tremors, Very Low Frequency Earthquakes, Low Frequency Earthquakes, Slow Slip Events, and earthquakes can all emerge naturally within this fault volume framework, making it an ideal digital twin for testing hypotheses, performing ground-truth inversions, and probing mechanical properties inaccessible with natural observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript proposes a 2D fault-zone model in which a main fault is embedded in an off-fault damage zone populated by distributed fractures obeying power-law size and density distributions. Fractures are restricted to optimal or main-fault-parallel orientations and all surfaces (main fault and off-fault cracks) obey rate-and-state friction with spatially uniform parameters chosen to permit slip instabilities. Quasi-dynamic boundary-integral simulations accelerated by hierarchical matrices are used to generate slip events both on and off the main fault. The central claim is that this geometry alone, without frictional heterogeneity, produces a natural continuum from slow to fast ruptures together with Omori and inverse-Omori laws, Gutenberg-Richter scaling, moment-duration scaling, pre-nucleation seismicity localization, and fluid-like migration patterns.

Significance. If the quantitative matches and robustness claims are substantiated, the work would provide a unified mechanical explanation for the full slip spectrum and associated statistical laws emerging from off-fault damage rather than from frictional heterogeneity. The framework could serve as a useful digital twin for hypothesis testing and ground-truth inversions. The technical implementation (hierarchical-matrix acceleration of boundary-integral methods) is a positive methodological contribution.

major comments (3)
  1. [Results] Results section (and associated figures): the manuscript asserts reproduction of the Omori law, inverse Omori law, Gutenberg-Richter scaling, and moment-duration scaling, yet supplies no quantitative goodness-of-fit metrics, fitted exponents with uncertainties, or direct comparison statistics between simulated and target distributions. Without these measures it is impossible to judge how well the simulated catalogs actually match the claimed laws.
  2. [Methods, Results] Methods and Results sections: the central claim that the observed behaviors emerge generically from any off-fault damage zone depends on the specific power-law exponents for crack size and density, the restriction to optimal or parallel orientations, and the particular rate-and-state parameters. No sensitivity tests varying these exponents, densities, or randomizing orientations are reported, leaving open the possibility that the statistical laws are tied to these geometric and parametric choices rather than arising generically.
  3. [§3] §3 (numerical method): while hierarchical-matrix acceleration is used to handle the large number of off-fault cracks, the manuscript does not quantify the approximation error introduced by the hierarchical compression for the quasi-dynamic evolution or demonstrate convergence with respect to the number of cracks. This is load-bearing for the reliability of the reported slow-to-fast continuum and migration patterns.
minor comments (2)
  1. [Figures] Figure captions and text: several panels show event catalogs or migration patterns; the distinction between main-fault and off-fault events should be made visually clearer (e.g., by symbol or color) and the time windows used for each panel should be stated explicitly.
  2. [Methods] Notation: the symbols used for crack density, size-distribution exponent, and the rate-and-state parameters (a, b, Dc) are introduced without a consolidated table; a single parameter table would improve readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed review. The comments identify key areas where additional quantification and robustness checks will strengthen the manuscript. We respond to each major comment below and commit to revisions that address the concerns without altering the core claims.

read point-by-point responses

  1. Referee: [Results] Results section (and associated figures): the manuscript asserts reproduction of the Omori law, inverse Omori law, Gutenberg-Richter scaling, and moment-duration scaling, yet supplies no quantitative goodness-of-fit metrics, fitted exponents with uncertainties, or direct comparison statistics between simulated and target distributions. Without these measures it is impossible to judge how well the simulated catalogs actually match the claimed laws.

    Authors: We agree that quantitative metrics are necessary for rigorous evaluation. In the revised manuscript we will add fitted exponents with uncertainties (e.g., b-value for Gutenberg-Richter and p-value for Omori’s law) together with goodness-of-fit statistics such as R² for log-log regressions and Kolmogorov-Smirnov p-values comparing simulated versus theoretical distributions. These will appear in the Results section, in new table entries, and in updated figure captions. revision: yes

  2. Referee: [Methods, Results] Methods and Results sections: the central claim that the observed behaviors emerge generically from any off-fault damage zone depends on the specific power-law exponents for crack size and density, the restriction to optimal or parallel orientations, and the particular rate-and-state parameters. No sensitivity tests varying these exponents, densities, or randomizing orientations are reported, leaving open the possibility that the statistical laws are tied to these geometric and parametric choices rather than arising generically.

    Authors: The chosen power-law exponents and orientations are directly motivated by field observations of damage zones; the rate-and-state parameters are uniform and lie within the unstable regime. Nevertheless, we recognize the value of explicit sensitivity checks. We will add a dedicated subsection (or appendix) presenting results for modest variations in the size and density exponents (±0.2) and will explicitly discuss the physical rationale for restricting orientations to optimal and main-fault-parallel directions. Full randomization of all orientations remains computationally prohibitive for the largest catalogs but will be addressed qualitatively. revision: partial

  3. Referee: [§3] §3 (numerical method): while hierarchical-matrix acceleration is used to handle the large number of off-fault cracks, the manuscript does not quantify the approximation error introduced by the hierarchical compression for the quasi-dynamic evolution or demonstrate convergence with respect to the number of cracks. This is load-bearing for the reliability of the reported slow-to-fast continuum and migration patterns.

    Authors: We performed internal convergence and error analyses during code development but omitted the quantitative details for brevity. In the revision we will expand §3 with explicit error metrics (relative L2-norm differences versus direct summation on benchmark problems) and will demonstrate convergence of key observables (event rates, migration speeds, slip-rate spectra) with increasing crack number and hierarchical rank. These additions will directly support the reliability of the slow-to-fast continuum and fluid-like migration patterns. revision: yes

Circularity Check

0 steps flagged

No significant circularity; forward simulation with independent inputs

full rationale

The paper describes a forward-modeling study: off-fault cracks are assigned explicit power-law size/density distributions and fixed orientations (optimal or parallel), all fractures receive uniform rate-and-state parameters chosen to allow instabilities, and quasi-dynamic boundary-integral simulations are run. The observed continuum of slip speeds, Omori laws, GR scaling, and moment-duration scaling are reported as emergent outputs of these simulations. No parameter is fitted to the target statistics and then relabeled as a prediction; no self-citation supplies a load-bearing uniqueness theorem; no ansatz is smuggled in; and the central claim does not reduce by construction to its own inputs. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 1 invented entities

The model rests on standard rate-and-state friction, power-law distributions for crack size and density, and the assumption that these choices suffice to produce the observed behaviors without additional heterogeneity.

free parameters (2)
  • power-law exponents for crack size and density
    Chosen to characterize the off-fault damage zone; specific values not stated in abstract but required to set up the fracture population.
  • rate-and-state friction parameters enabling slip instabilities
    Selected so that fractures can produce both slow and fast slip; these are free parameters tuned to allow instabilities.
axioms (2)
  • domain assumption Rate-and-state friction governs slip on all fractures
    Invoked to allow slip instabilities on both main fault and off-fault cracks.
  • domain assumption Quasi-dynamic approximation is sufficient for the slip dynamics
    Used in the boundary integral simulations.
invented entities (1)
  • off-fault damage zone as digital twin no independent evidence
    purpose: To reproduce the full slip spectrum and statistical laws
    The damage zone is the central modeling construct introduced to generate the observed behaviors.

pith-pipeline@v0.9.0 · 5842 in / 1531 out tokens · 31421 ms · 2026-05-18T19:07:28.619347+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

10 extracted references · 10 canonical work pages

  1. [1]

    Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard

    Abdelmeguid, M. & A. Elbanna (2022). “Sequences of seismic and aseismic slip on bimaterial faults show dominant rupture asymmetry and potential for elevated seismic hazard”.Earth and Planetary Science Letters593, p. 117648.DOI: 10.1016/j.epsl.2022.117648. Abercrombie, R. E. & J. R. Rice (2005). “Can observations of earthquake scaling constrain slip weaken...

    work page doi:10.1016/j.epsl.2022.117648 2022

  2. [2]

    Characterization of barriers on an earthquake fault

    AGU.DOI:10.1029/ME004p0566. Aki, K. (1979). “Characterization of barriers on an earthquake fault”.J. Geophys. Res.84.B11, pp. 6140–6148.DOI:10.1029 /JB084iB11p06140. Amezawa, Y., T. Maeda & M. Kosuga (2021). “Migration diffusivity as a controlling factor in the duration of earthquake swarms”.Earth, Planets and Space73.1.DOI:10.1186/s40623-021-01480-7. Amp...

    work page doi:10.1029/me004p0566 1979

  3. [3]

    Dynamic plane-strain shear rupture with a slip-weakening friction law calculated by a boundary integral method

    Time-dependent case”.J. Geophys. Res.86.B11, pp. 10821–10834.DOI: 10.1029/JB086iB11p10821. Andrews, D. J. (1985). “Dynamic plane-strain shear rupture with a slip-weakening friction law calculated by a boundary integral method”.Bulletin of the Seismological Society of America75.1, pp. 1–21.DOI:10.1785/bssa0750010001. Andrews, D. J. (2005). “Rupture dynamic...

    work page doi:10.1029/jb086ib11p10821 1985

  4. [4]

    Earthquake nucleation on faults with rate-and state-dependent strength

    Experimental results and constitutive equations”.Journal of Geophysical Research: Solid Earth84.B5, pp. 2161–2168.DOI:10.1029/jb084ib05p02161. Dieterich, J. H. (1992). “Earthquake nucleation on faults with rate-and state-dependent strength”.Tectonophysics211.1-4, pp. 115–134.DOI:10.1016/0040-1951(92)90055-b. Dor, O., Y. Ben-Zion, T. K. Rockwell & J. N. Br...

    work page doi:10.1029/jb084ib05p02161 1992

  5. [5]

    Slow earthquake scaling reconsidered as a boundary between distinct modes of rupture propagation

    Ide, S. & G. C. Beroza (2023). “Slow earthquake scaling reconsidered as a boundary between distinct modes of rupture propagation”.Proceedings of the National Academy of Sciences120.32.DOI:10.1073/pnas.2222102120. Ide, S., G. C. Beroza, D. R. Shelly & T. Uchide (2007). “A scaling law for slow earthquakes”.Nature447.7140, pp. 76–79. DOI:10.1038/nature05780....

    work page doi:10.1073/pnas.2222102120 2023

  6. [6]

    Community-Driven Code Comparisons for Three-Dimensional Dynamic Modeling of Sequences of Earthquakes and Aseismic Slip

    Jiang, J., B. A. Erickson, V . R. Lambert, J. Ampuero, R. Ando, S. D. Barbot, C. Cattania, L. D. Zilio, B. Duan, E. M. Dunham, A. Gabriel, N. Lapusta, D. Li, M. Li, D. Liu, Y. Liu, S. Ozawa, C. Pranger & Y. van Dinther (2022). “Community-Driven Code Comparisons for Three-Dimensional Dynamic Modeling of Sequences of Earthquakes and Aseismic Slip”.Journal o...

    work page doi:10.1029/2021jb023519 2022

  7. [7]

    The physics of earthquakes

    Kanamori, H. & E. E. Brodsky (2004). “The physics of earthquakes”.Reports on Progress in Physics67.8, pp. 1429–1496.DOI: 10.1088/0034-4885/67/8/r03. Kato, A. & Y. Ben-Zion (2020). “The generation of large earthquakes”.Nature Reviews Earth & Environment2.1, pp. 26–39. DOI:10.1038/s43017-020-00108-w. Kato, N., M. Ohtake & T. Hirasawa (1997). “Possible Mecha...

    work page doi:10.1088/0034-4885/67/8/r03 2004

  8. [8]

    Dynamics, Radiation, and Overall Energy Budget of Earthquake Rupture With Coseismic Off-Fault Damage

    Okubo, K., H. S. Bhat, E. Rougier, S. Marty, A. Schubnel, Z. Lei, E. E. Knight & Y. Klinger (2019). “Dynamics, Radiation, and Overall Energy Budget of Earthquake Rupture With Coseismic Off-Fault Damage”.Journal of Geophysical Research: Solid Earth124.11, pp. 11771–11801.DOI:10.1029/2019jb017304. Okubo, K., E. Rougier, Z. Lei & H. S. Bhat (2020). “Modeling...

    work page doi:10.1029/2019jb017304 2019

  9. [9]

    Characteristics of Slow Slip Events Explained by Rate-Strengthening Faults Subject to Periodic Pore Fluid Pressure Changes

    Perez-Silva, A., Y. Kaneko, M. Savage, L. Wallace & E. Warren-Smith (2023). “Characteristics of Slow Slip Events Explained by Rate-Strengthening Faults Subject to Periodic Pore Fluid Pressure Changes”.Journal of Geophysical Research: Solid Earth128.6.DOI:10.22541/essoar.167327976.60438163/v1. Perfettini, H., J. Schmittbuhl, J. R. Rice & M. Cocco (2001). “...

    work page doi:10.22541/essoar.167327976.60438163/v1 2023

  10. [10]

    Fault rocks and fault mechanisms

    Sibson, R. H. (1977). “Fault rocks and fault mechanisms”.Journal of the Geological Society (London, United Kingdom)133.3, pp. 191–213.DOI:10.1144/gsjgs.133.3.0191. Sibson, R. H. (1986). “Rupture interaction with fault jogs”.Earthquake Source Mechanics. Ed. by S. Das, J. Boatwright & C. H. Scholz. AGU Geophys. Monogr., pp. 157–167.DOI:10.1029/gm037p0157. S...

    work page doi:10.1144/gsjgs.133.3.0191 1977