Rough Subducting Seafloor Reduces Interseismic Coupling and Mega-Earthquake Occurrence: Insights From Analogue Models
Pith reviewed 2026-05-25 19:22 UTC · model grok-4.3
The pith
Rough subduction interfaces reduce megathrust coupling and limit mega-earthquake size in analogue models.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Models characterized by a very rough interface have lower interface frictional strength and lower interseismic coupling than models with a smooth interface. Overall, ruptures in the rough models have smaller rupture area, duration and mean displacement. Individual slip distributions indicate a segmentation of the subduction interface by the rough geometry. Flexure of the overriding plate is proposed as one mechanism contributing to the heterogeneous strength distribution.
What carries the argument
Seismotectonic analogue models comparing rough and smooth subduction interfaces, which demonstrate heterogeneous strength distribution and reduced coupling due to roughness.
If this is right
- Rough interfaces produce lower interseismic coupling than smooth ones.
- Ruptures have smaller area, duration, and mean displacement.
- The rough geometry segments the interface, leading to heterogeneous slip distributions.
Where Pith is reading between the lines
- Subduction zones with smoother seafloor may be more prone to mega-earthquakes if roughness controls coupling.
- Plate flexure effects on strength could be isolated in follow-up models that vary overriding plate thickness.
- Field measurements of coupling differences between rough and smooth natural zones would test the scaling from the models.
Load-bearing premise
The analogue models accurately reproduce the stress state and scaling of natural subduction megathrusts, including the effects of plate flexure on interface strength.
What would settle it
Direct comparison of interseismic coupling measurements from GPS data in natural subduction zones with rough versus smooth subducting seafloor.
Figures
read the original abstract
The roughness of the subduction interface is thought to influence seismogenic behavior in subduction zones, but a detailed understanding of how such roughness affects the state of stress along the subduction megathrust is still debated. Here, we use seismotectonic analogue models to investigate the effect of subduction interface roughness on seismicity in subduction zones. We compared analogue earthquake source parameters and slip distributions for two roughness endmembers. Models characterized by a very rough interface have lower interface frictional strength and lower interseismic coupling than models with a smooth interface. Overall, ruptures in the rough models have smaller rupture area, duration and mean displacement. Individual slip distributions indicate a segmentation of the subduction interface by the rough geometry. We propose that flexure of the overriding plate is one of the mechanisms that contribute to the heterogeneous strength distribution, responsible for the observed seismic behavior.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses seismotectonic analogue models to compare subduction megathrust behavior for smooth versus very rough interface end-members. It reports that rough-interface models exhibit lower interface frictional strength and interseismic coupling, with ruptures showing smaller area, duration, and mean displacement; slip distributions indicate interface segmentation by roughness, and overriding-plate flexure is proposed as one mechanism producing heterogeneous strength.
Significance. If the analogue scaling holds, the work supplies direct experimental evidence that seafloor roughness can segment the megathrust, reduce coupling, and limit mega-earthquake size—observations that are difficult to obtain from natural data alone. The ability to visualize full slip distributions across repeated events is a clear strength of the analogue approach and could inform interpretations of coupling maps and paleoseismic records in subduction zones with varying basement roughness.
major comments (2)
- [Abstract and Discussion] The central claim that rough models have lower frictional strength and coupling rests on the assumption that the analogue setup reproduces the stress state and scaling of natural megathrusts, including plate-flexure effects invoked in the abstract. No dimensionless analysis, stress-state comparison, or scaling checks are described that would confirm this reproduction; without them the reported differences cannot be confidently extrapolated beyond the laboratory.
- [Results] Results: the reported differences in coupling coefficient, rupture area, duration, and mean displacement between the two roughness end-members are presented without error bars, standard deviations across realizations, or statistical tests. This makes it impossible to judge whether the differences are robust to model variability or specific to the chosen material properties and loading rates.
minor comments (2)
- [Figures] Figure captions should explicitly state the roughness amplitude and wavelength values used for each end-member and the number of experimental runs performed for each case.
- [Methods] The term 'interseismic coupling' is used for the analogue models; a brief clarification of how this quantity is defined and measured from the laboratory time series would aid readers unfamiliar with the method.
Simulated Author's Rebuttal
We thank the referee for their constructive comments, which highlight important aspects of our analogue modeling study. We respond to each major comment below.
read point-by-point responses
-
Referee: [Abstract and Discussion] The central claim that rough models have lower frictional strength and coupling rests on the assumption that the analogue setup reproduces the stress state and scaling of natural megathrusts, including plate-flexure effects invoked in the abstract. No dimensionless analysis, stress-state comparison, or scaling checks are described that would confirm this reproduction; without them the reported differences cannot be confidently extrapolated beyond the laboratory.
Authors: We agree that an explicit dimensionless analysis and stress-state comparison were not included in the original manuscript. In revision we will add a dedicated scaling section that presents the relevant non-dimensional numbers (e.g., force ratios governing gravitational, frictional and viscous stresses) and compares the modeled interface stresses and overriding-plate flexure magnitudes to typical natural values reported in the literature. The flexure mechanism itself is observed directly in the experiments and is offered as one contributing process; the added scaling discussion will clarify the conditions under which the laboratory results can be extrapolated. revision: yes
-
Referee: [Results] Results: the reported differences in coupling coefficient, rupture area, duration, and mean displacement between the two roughness end-members are presented without error bars, standard deviations across realizations, or statistical tests. This makes it impossible to judge whether the differences are robust to model variability or specific to the chosen material properties and loading rates.
Authors: The original manuscript reports results from repeated experiments for each roughness end-member, yet quantitative variability measures were omitted. In the revised version we will include standard deviations (or error bars) on the reported parameters and add a short statistical comparison (e.g., two-sample t-tests) between the smooth and rough populations to demonstrate that the observed differences are robust across realizations. revision: yes
Circularity Check
No significant circularity; purely experimental comparison
full rationale
The paper reports direct observations from analogue models comparing two roughness end-members. Central claims (lower frictional strength, lower coupling, smaller ruptures in rough models) are empirical outcomes of the physical experiments, with no derivations, fitted parameters renamed as predictions, or load-bearing self-citations. The study is self-contained against external benchmarks as an experimental result set; no equations or ansatzes reduce the findings to inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Analogue models can be scaled to reproduce the stress state and rupture behavior of natural subduction megathrusts
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Models characterized by a very rough interface have lower interface frictional strength and lower interseismic coupling... ruptures... smaller rupture area, duration and mean displacement... segmentation of the subduction interface by the rough geometry... flexure of the overriding plate
-
IndisputableMonolith/Foundation/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We use seismotectonic analogue models... 3D-printed geometry... stick-slip behavior
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
-
[1]
https://doi.org/10.1130/G37258.1 Husen, S., Kissling, E., & Quintero, R. (2002). Tomographic evidence for a subducted seamount beneath the Gulf of Nicoya, Costa Rica: The cause of the 1990 Mw = 7.0 Gulf of Nicoya earthquake. Geophysical Research Letters, 29(8), 79-1-79–4. https://doi.org/10.1029/2001GL014045 Kelleher, J., & McCann, W. (1976). Buoyant Zone...
-
[2]
https://doi.org/10.1130/0091-7613(1997)025<0487:TEOSSO>2.3.CO;2 Singh, S. C., Hananto, N., Mukti, M., Robinson, D. P., Das, S., Chauhan, A., … Harjono, H. (2011). Aseismic zone and earthquake segmentation associated with a deep subducted seamount inSumatra. Nature Geoscience, 4(5), 308–311. https://doi.org/10.1038/ngeo1119 Thielicke, W., & Stamhuis, E. J....
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.