Pith. sign in

REVIEW 4 major objections 7 minor 3 references

Two parallel, overlapping faults with uniform rate-weakening friction can generate both slow slip events and earthquakes, whereas a single planar fault under the same conditions produces only earthquakes.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Two parallel faults with uniform rate-weakening friction generate slow slip events, ordinary earthquakes, and complex sequences, with slow slip only at intermediate geometric interaction strength.

T0 review reviewed 2026-08-02 challenge →

load-bearing objection Solid 3D step-over simulations show geometry alone can produce slow slip; the new Λ metric is a useful organizing device but its 2D basis needs direct validation against the 3D kernel. the 4 major comments →

arxiv 2602.16403 v2 pith:AY7B7ZBB submitted 2026-02-18 physics.geo-ph

Quantifying the Role of 3D Fault Geometry Complexities on Slow and Fast Earthquakes

J. Cheng , H. S. Bhat , M. Almakari , B. Lecampion , P. Dubernet This is my paper
classification physics.geo-ph
keywords slow slip eventsearthquake sequence simulationfault geometryfault interaction metricrate-and-state frictionstep-over faultsmoment-duration scaling3D quasi-dynamic model
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

The paper sets out to show that fault geometry alone—specifically, the interaction between two parallel, overlapping faults in a compressional step-over—can explain the coexistence of slow slip events and ordinary earthquakes, without invoking heterogeneous friction, fluids, or other extra physics. In 3D quasi-dynamic earthquake-cycle simulations with uniform velocity-weakening rate-and-state friction, the authors find four slip regimes: only slow slip, slow-slip-dominant, earthquake-dominant, and only earthquakes. A single planar fault under identical friction produces only earthquakes. They introduce a geometry-only interaction metric Λ, the maximum Coulomb stress induced on one fault by a unit uniform stress drop on the other, and find that slow slip occurs only for intermediate interaction strengths (0.03 < Λ < 1.2), with low strengths giving periodic earthquakes and high strengths giving complex earthquake sequences. The paper also reproduces observed moment-duration scaling and shows that its exponent depends on the slip-rate threshold used to detect events, implying a detection-sensitivity effect in real observations.

Core claim

On the paper's own terms, the discovery is that geometric interaction alone is enough to split a single slip behavior into the full observed spectrum: with two parallel, overlapping faults in a 3D quasi-dynamic rate-and-state model with spatially uniform velocity-weakening friction, slow slip events and earthquakes coexist, while an isolated planar fault with the same friction produces only periodic earthquakes. The mechanism is the time-dependent evolution of traction heterogeneities—'traction asperities'—created by one fault's stress field on the other. These asperities control where nucleation can complete: narrow high-ratio patches produce slow, partial ruptures, and low-ratio patches ar

What carries the argument

The load-bearing object is the dimensionless fault-interaction metric Λ, defined as the maximum Coulomb stress change induced on a receiving fault by a unit, spatially uniform stress drop on the neighboring fault. Λ is computed from the elastic stress field of a 2D shear crack via complex potentials, normalized so it depends only on geometry—fault length, width, overlap, and separation—and is independent of friction and nucleation length. Its role is to collapse the many geometric degrees of freedom of a step-over onto a single axis; the simulations show that the SSE moment-release ratio is controlled by Λ, with the boundaries at 0.03 and 1.2. The second mechanism is the evolving 'traction a

Load-bearing premise

The load-bearing premise is that the 2D crack-based metric Λ—computed from a unit uniform stress drop with a fixed friction coefficient and a crack half-length equal to the smaller of fault length and width—faithfully represents the 3D stress interaction between finite rectangular faults; if 3D finite-size or normal-stress coupling effects alter the interaction, the claimed intermediate window for slow slip could be an artifact of the tested geometries.

What would settle it

A single 3D simulation set that computes fault interaction with full 3D stress Green's functions rather than the 2D crack formula, run across the same geometries, would settle it: if slow slip events appear for Λ < 0.03 or Λ > 1.2, or fail to appear inside that window, the universal window is not general. In the field, a fault segment with Λ far outside 0.03–1.2 that nonetheless shows recurring slow slip would also contradict the claim.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Real segmented subduction faults can produce slow slip events without requiring heterogeneous friction or fluids; geometric interaction alone is sufficient.
  • The Λ metric, being friction-independent and purely geometric, can be applied to arbitrary fault systems to estimate where slow slip is likely to occur.
  • Observed variations in SSE moment-duration scaling across studies may reflect event-detection thresholds as much as physical differences in the faults.
  • Weakly interacting or strongly interacting fault segments will be earthquake-dominated, so only intermediate geometric coupling hosts slow slip.
  • Large fault width combined with small separation pushes a system into irregular, complex earthquake sequences rather than slow slip.

Where Pith is reading between the lines

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

  • Beyond the paper, Λ could be computed from imaged fault maps in subduction zones; segments falling in 0.03 < Λ < 1.2 should be preferentially slow-slip-prone, a testable field prediction for places like Cascadia and Hikurangi.
  • Beyond the paper, the sharp numerical boundaries 0.03 and 1.2 are likely sensitive to the choice of friction coefficient, loading rate, and fault aspect ratio; in nature the window may be broader or shifted, so the universal statement should be read as 'intermediate coupling' rather than exact numbers.
  • Beyond the paper, because Λ derives from a 2D crack solution, applying it to non-parallel or curved faults is an unproven extrapolation; a full 3D stress-interaction metric would either confirm or revise the window.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 7 minor

Summary. The paper presents 3D quasi-dynamic earthquake sequence simulations of two parallel, strike-slip-like faults in a compressional step-over, with spatially uniform rate-weakening rate-and-state friction. It reports four slip regimes—periodic earthquakes, SSE-dominant, earthquake-dominant, and complex sequences—whereas a single planar fault under identical friction produces only periodic earthquakes (Fig. S1). A scalar interaction metric Λ is introduced (Appendix A, Eq. A.6), defined via the maximum Coulomb stress induced by a unit uniform stress drop on a 2D crack, and the authors find that SSEs occur only in an intermediate interaction range, approximately 0.03 < Λ < 1.2, with SSE dominance for 0.05 < Λ < 0.75 (Fig. 3b). The paper also reports moment-duration scaling for SSEs that is threshold-dependent (M∼T^1.2 at 1e-8 m/s vs. M∼T^1.5 at 1e-6 m/s) and cubic scaling for earthquakes. The authors interpret the results as demonstrating that geometrical complexity alone can generate both slow and fast earthquakes via evolving traction heterogeneities.

Significance. If the central claim holds, the paper offers a mechanistically simple and observationally attractive route to slow slip: geometric fault interaction, without heterogeneous friction, fluids, or dilatant strengthening, can produce SSEs and their coexistence with earthquakes. The study's strengths include a relatively large parameter sweep (65 simulations), a single-fault control case, the use of a modern 3D quasi-dynamic BEM accelerated by hierarchical matrices, and public availability of code and data. The threshold sensitivity analysis in Section 3.2 is a useful caveat for interpreting scaling laws from synthetic and observed catalogs. However, the main conclusion rests on the validity of the 2D crack-based Λ metric as an ordering parameter for 3D fault interaction, and on the robustness of post hoc regime boundaries; these need stronger support before the universal intermediate-range claim can be accepted.

major comments (4)
  1. [Appendix A, Eq. (A.6); Section 4.1, Fig. 3] The load-bearing premise is that the 2D Muskhelisvili-Kolosov crack solution, evaluated with half-length a = min(Lf,W) and fs=0.6, faithfully represents the 3D stress interaction between finite rectangular faults. The simulations solve full 3D quasi-dynamic elastostatics with normal-stress coupling, yet no comparison is made between Λ computed from the 2D formula and the actual 3D BEM kernel. The decay of stress with distance differs between 2D cracks (~1/sqrt(r)) and 3D finite faults (closer to 1/r^2 in shadow zones), so Λ may misorder configurations, especially for large Lf/W. Supporting Figures S5-S6 validate only the 2D crack stress field, not its equivalence to the 3D transfer used in the simulations. Please add a direct scatter plot of Λ_2D vs. the maximum Coulomb stress computed from the 3D kernel for the same geometries, and test whether the 0.03 and 1.2 boundaries survive altern
  2. [Section 3.2, Fig. 2] The claim to reproduce observed moment-duration scaling is weakened by the simultaneous dependence on detection threshold and by the absence of uncertainties. The text reports 1925 SSEs with M∼T^1.2 at 1e-8 m/s and 640 SSEs with M∼T^1.5 at 1e-6 m/s, but no error bars, goodness-of-fit, or event-to-event scatter are quantified. Because the analysis uses only 10 single-fault recurrence intervals per simulation, the catalog is small, especially for the slower and rarer events that dominate the reported scaling. Please provide bootstrap confidence intervals for the fitted exponents and show the scaling for individual simulations or a clear aggregate with uncertainties.
  3. [Section 4.1, Fig. 3(b)] The regime boundaries Λ≈0.03 and Λ≈1.2 are assigned post hoc from the simulation data, and the supporting discussion of proportionality to 1/sqrt(D) and sqrt(W) is qualitative. As currently presented, the claim that SSEs occur only at intermediate interaction strength is a data-fitting statement, not a validated predictive law. The paper would be strengthened by an out-of-sample test—e.g., withholding a subset of the 65 geometries, fitting the boundaries on the rest, and checking classification—or by a derivation of the boundaries from nucleation criteria. Additionally, only three values of a/b are used; the statement that 'geometry remains the primary controlling factor' should be quantified with a regression or uncertainty measure.
  4. [Section 2; Section 3.2] The simulation duration is normalized to the time for ten earthquakes on a single planar fault. For SSE-dominant runs, which by definition contain few earthquakes, ten-equivalent cycles may be too short to sample the long-term statistics, particularly for recurrence variability and for the tail of the moment-duration relation. No convergence test or sensitivity to simulation length is reported. This is a load-bearing issue for the quantitative boundaries of the regimes, not merely a technical detail.
minor comments (7)
  1. [Plain Language Summary] “the transition between slip behaviors depends on how strongly the faults interact each other” should read “interact with each other.”
  2. [Figure 1(d) caption] “nomalized slip” is a typo for “normalized slip.”
  3. [Appendix A and Fig. S6] The friction coefficient is written fs in Eq. (A.6) and the main text, but µ in Fig. S6. Use a single symbol to avoid confusion with the shear modulus µ.
  4. [Throughout] “Muskhelisvili” is a misspelling of “Muskhelishvili.” Also, Λ is sometimes written as λ (e.g., Fig. S2, S3 captions); unify the notation.
  5. [Fig. 3(b)] The axis label “SSEs ratio” should read “SSE ratio” or “χ_sse” for consistency with the definition in Section 4.1.
  6. [Section 3.2] The text first says SSEs follow linear scaling and later says “near-linear” or “predominantly linear” with exponents 1.2–1.5. Please clarify whether linear scaling is claimed or whether the exponents are simply less than the earthquake cubic value.
  7. [Figure 4] The HRA/LRA (high/low traction-ratio asperity) concepts are central to the proposed mechanism, but the thresholds for “high” and “low” are not defined quantitatively. Please state the criterion used to classify a patch as HRA or LRA.

Circularity Check

0 steps flagged

No significant circularity: Λ is an independently defined geometric metric, SSE-regime boundaries are empirical, and self-citations are methodological.

full rationale

The paper's central claim—that two parallel, uniformly rate-weakening faults generate slow slip events only at intermediate interaction strength—is not circular. Λ (Eq. 3 / Appendix A Eq. A.6) is constructed from 2D crack complex potentials under a normalized unit stress drop, using only geometry plus an assumed f_s=0.6; it is defined before any simulated slip output is used. The SSE ratio χ_sse is measured from independent 3D quasi-dynamic simulations using a slip-rate threshold. Plotting χ_sse against Λ and reading empirical boundaries (0.03<Λ<1.2) is a correlation/order-parameter analysis, not a reduction of the outcome to the metric by construction; no RSF parameter is fit to the regime boundaries. The paper's self-citations are not load-bearing: FASTDASH (Cheng et al. 2025) is a published, code-reproducible numerical method; Romanet et al. (2018) is prior 2D motivation, not the proof of the 3D result; Cheng et al. (2026) is data availability. No uniqueness theorem is imported, and the complex-potential formulas are attributed to Scheel et al. (2021). The identified weakness—Λ uses a 2D crack with half-length min(L_f,W) instead of the full 3D BEM kernel—is a validity/robustness risk, not circularity, because Λ is not defined in terms of the simulated SSE behavior it is used to classify.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 1 invented entities

The central claim rests on standard rate-and-state friction plus a geometric interaction metric; no physical constants are fitted to target observations, but the regime boundaries and event thresholds are empirical.

free parameters (4)
  • Λ regime boundaries = Λ < 0.03, Λ > 1.2, SSE-dominant 0.05 < Λ < 0.75
    Determined after the fact from the 65 simulated catalogs (Fig. 3b) to separate only-EQ, SSE-dominated, and mixed regimes; not predicted before running the simulations.
  • Slip-rate event thresholds = EQ: ≥10^-3 m/s; SSE: 10^-8 to 10^-3 m/s (also 10^-6 threshold)
    Chosen by hand; the fitted moment-duration exponent changes from 1.2 to 1.5 when the SSE threshold changes, so the scaling result is threshold-dependent.
  • Friction coefficient fs in Λ (Coulomb stress) = 0.6
    Assumed constant in Appendix A Eq. A.6; standard Byerlee value, not fitted to simulations.
  • Fault aspect ratio Lf/W = 3 (fixed)
    Width and length are varied jointly while keeping aspect ratio fixed; the claim that Λ, not individual dimensions, controls behavior may depend on this choice.
axioms (4)
  • domain assumption Spatially uniform rate-weakening rate-and-state friction (a/b in {0.4, 0.6, 0.8}) with aging law
    Central to the claim that any slow slip is caused by geometry, not friction heterogeneity; real faults may have heterogeneous frictional properties.
  • domain assumption Quasi-dynamic approximation (radiation damping; no full wave propagation)
    Standard in cycle modeling; dynamic stress waves could affect complex sequences and fast earthquake arrest, but are omitted here.
  • domain assumption Two parallel finite planar faults in a homogeneous elastic medium subjected to uniform far-field loading
    Idealized geometry; real subduction interfaces are non-planar, segmented, and influenced by fluids.
  • ad hoc to paper The 2D Muskhelisvili-Kolosov crack solution (with half-length a = min(Lf, W)) adequately approximates 3D fault-fault stress transfer for the metric Λ
    Appendix A/Eq. A.6; no direct comparison to the 3D BEM stress kernel is shown, so the metric's sufficiency is assumed.
invented entities (1)
  • Traction asperity (HRA/LRA patches) no independent evidence
    purpose: To explain where SSEs and earthquakes nucleate and arrest on faults with uniform frictional properties: high shear-to-normal traction ratio patches nucleate events; low-ratio patches act as barriers.
    Introduced in §4.2 based on simulated traction fields; no external data or testable prediction is given, so it is a descriptive concept rather than a demonstrated physical entity.

reviewed 2026-08-02 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Quantifying the Role of 3D Fault Geometry Complexities on Slow and Fast Earthquakes." pith.science (2026). https://pith.science/paper/AY7B7ZBB

@misc{pith2026260216403,
  author       = {Pith},
  title        = {Pith review of: Quantifying the Role of 3D Fault Geometry Complexities on Slow and Fast Earthquakes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AY7B7ZBB}},
  note         = {Machine review of arXiv:2602.16403}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Traditional models of slow slip events (SSEs) oversimplify fault geometry, although imaging shows subduction faults are segmented and complex. We examine how fault interactions control slip behavior using 3-D quasi-dynamic earthquake sequence simulations of two parallel faults with uniform rate-weakening friction accelerated by hierarchical matrices. Four regimes emerge-periodic earthquakes, coexisting SSEs and earthquakes, only SSEs, and complex sequences-whereas with the same friction condition a single planar fault produces only earthquakes. We quantify interaction using the maximum Coulomb stress induced on a target fault by a spatially uniform unit stress drop on a neighboring fault. Because the stress drop is normalized, the metric depends only on geometry and is independent of friction, allowing extension to arbitrary fault systems. SSEs occur only at intermediate fault interaction strengths. At low interaction strengths, the system produces regular, periodic earthquakes. At high interaction strengths, fault interactions generate complex earthquake sequences with irregular recurrence and variable magnitudes. Simulations reproduce observed moment-duration scaling and show sensitivity to detection thresholds. These results demonstrate geometric complexity alone generates both slow and fast earthquakes through evolving traction heterogeneity.

Figures

Figures reproduced from arXiv: 2602.16403 by B. Lecampion, H. S. Bhat, J. Cheng, M. Almakari, P. Dubernet.

Figure 1
Figure 1. Figure 1: (a) Step-over fault configuration. The mesh is shown exaggerated for visualization purposes; all simulations use a grid spacing of Lb/3 (b-c) Example of SSE-dominant regime and earthquake-dominant regime. Time evolution of the maximum slip rate. The bottom x-axis shows real time, while the top x-axis shows time normalized by the recurrence interval of a single fault with the same frictional properties. Gre… view at source ↗
Figure 2
Figure 2. Figure 2: Moment-Duration scaling across all simulations. Red denotes earthquakes exhibiting cubic scaling. Blue and black represent slow slip events identified with slip rate thresholds of 10−6m/s and 10−8m/s, respectively, showing varying but predominantly linear scaling from M ∼ T 1.5 to M ∼ T 1.2 8 Feb. 19, 2026 at 01:43 [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Definition of Λ in relation to the geometric parameters L, D, and Lf (b) The relationship between the fault interaction metric Λ and the SSEs ratio χsse. Different colored dots represent various friction parameters a/b. The dashed lines indicate that SSEs emerge within a specific range of Λ. Grey area shows the results with SSE-dominant slip catalog. (c) Fault interaction metric Λ as a function of over… view at source ↗
Figure 4
Figure 4. Figure 4: Traction field evolution during four events (E, F, G, and H) on fault 2 in Model 2. The colormap indicates the ratio of shear to normal traction. The first column displays the traction field before the nucleation of each event, with red stars marking earthquake hypocenters and blue stars indicating SSE hypocenters. The second column shows the traction field after each event, with black contour lines repres… view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

3 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

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

    Almakari, M., N. Kheirdast, C. D. Villafuerte, M. Y . Thomas, P. Dubernet, J. Cheng, A. Gupta, P. Romanet, S. Chaillat & H. S. Bhat (2026). “Fault volume digital twin to reproduce the full slip spectrum, scaling and statistical laws”. In: under review J. Geophys. Res.DOI:10.48550/arXiv.2509.04909. arXiv:2509.04909. Ando, R., K. Ujiie, N. Nishiyama & Y . M...

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2509.04909 2026

  2. [3]

    Repeating slip events at a circular asperity: numerical simulation with a rate- and state-dependent friction law

    Kato, N (2003). “Repeating slip events at a circular asperity: numerical simulation with a rate- and state-dependent friction law”. In:Bull. Earthq. Res. Inst.78, pp. 151–166. Kirkpatrick, J. D., ˚A. Fagereng & D. R. Shelly (2021). “Geological constraints on the mechanisms of slow earthquakes”. In:Nature Reviews Earth Environment2.4, pp. 285–301.ISSN: 266...

    arXiv 2003

  3. [6]

    Fluid-induced aseismic fault slip outpaces pore-fluid migration

    Bhattacharya, P. & R. C. Viesca (2019). “Fluid-induced aseismic fault slip outpaces pore-fluid migration”. In:Science 364.6439, pp. 464–468.ISSN: 1095-9203.DOI:10.1126/science.aaw7354. Bletery, Q. & J.-M. Nocquet (2020). “Slip bursts during coalescence of slow slip events in Cascadia”. In:Nature Commu- nications11.1.ISSN: 2041-1723.DOI:10.1038/s41467-020-...

    arXiv 2019

This paper was first reviewed by deepseek-v4-flash on August 2, 2026.