On four representations of the law of aftershock evolution

A.V. Guglielmi; O.D. Zotov

REVIEW 4 major objections 5 minor 2 references

On four representations of the law of aftershock evolution

T0 review · 4 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash

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

Pith's one-line read The paper argues that aftershock activity decays exponentially with the earthquake source's 'proper time' — the clock defined by the shocks themselves — making Omori's and Utsu's laws only partial approximations.

desk verdict The exponential decay law is a definitional identity; the empirical support is too thin to carry the claim. read the letter →

arxiv 2607.29381 v1 pith:HOKGPVFL submitted 2026-07-31 physics.geo-ph

A.V. Guglielmi , O.D. Zotov This is my paper

classification physics.geo-ph

keywords aftershockdecayOmorilawUtsusourcepropertimedeactivationcoefficientundergroundclockTohokuearthquakerelaxation

verification ladder T0 review T1 audit T2 compute T3 formal

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 tries to establish that an earthquake source relaxes after a mainshock in a way that is invisible when aftershocks are timed by ordinary clocks: if the shocks themselves are used as ticks of an 'underground clock,' the interval between successive aftershocks grows exponentially, and aftershock frequency decays exponentially with the source's proper time. The authors compare four representations — Omori's hyperbolic law, the Hirano-Utsu power law, and their own linear and discrete forms — and argue that only the linear and discrete forms are complete, while the traditional laws hold only during an initial 'Omori epoch' or fail altogether. If correct, this reframes aftershock analysis from fitting a universal-time power law to measuring a single state parameter, the deactivation coefficient, which characterizes the source and appears to decrease with mainshock magnitude.

What carries the argument

The deactivation coefficient σ, a state variable of the source, and the source's proper time, defined as the accumulated integral of σ over universal time. In the discrete representation, the aftershock sequence itself serves as an 'underground clock' whose ticks mark proper time. The load-bearing identity is the synchronization function t(x) = t1 exp(σ(x−1)); its logarithmic derivative gives σ, and it converts a nonlinear Omori-type decay into the linear equation dn/dχ + n = 0.

What would settle it

Examine the early hours of the Tohoku sequence with a completeness correction or a local borehole catalog that detects every early aftershock: if the inter-event intervals do not form a geometric progression once completeness is accounted for, or if the σ from the discrete representation no longer matches the σ from the linear representation, the central claim is refuted. Alternatively, find an aftershock sequence with p far from 1 that still shows the claimed exponential proper-time decay; if the exponential fails there, the law is not general.

Watch

Extended reading notes

Core claim

The central claim is that aftershock activity decays exponentially with the source's proper time, not with universal time. In the discrete representation, the excitation times of successive aftershocks, plotted against their natural numbering, form an exponential synchronization curve t(x) = t1 exp(σ(x−1)); for the Tohoku sequence the paper reports t1 = 7.2 h, σ = 0.0014, and R² = 0.97, with the same σ recovered from the linear representation. This means the frequency obeys the first-order linear equation dn/dχ + n = 0, where χ = ∫σ dx is the source's proper time. The authors present this as a holistic law: Omori's hyperbolic law holds only while σ is constant, and the Hirano-Utsu power law

Load-bearing premise

The load-bearing premise is that the catalog supplies a complete, unbroken enumeration of aftershock excitation times, so that numbering them x_i = i is faithful and the interval between successive aftershocks stands for 1/n; missing early aftershocks or contaminating non-aftershocks would make the exponential synchronization an artifact, and the paper does not specify the magnitude cutoff, declustering, completeness period, or uncertainty for the 4,537 Tohoku events.

Editorial extensions

If this is right

If the claim is correct, the Omori and Hirano-Utsu laws are not fundamental; they are first-stage approximations valid only while the deactivation coefficient stays constant, and the ubiquitous p ≈ 1.1 scatter is an artifact of using universal time instead of source proper time.
The deactivation coefficient becomes a directly measurable, source-specific physical quantity, with an observed monotonic decrease as mainshock magnitude increases, making it a candidate state variable for earthquake source physics.
Aftershock evolution is two-stage: a predictable Omori epoch followed by a bifurcation into non-monotonic, apparently chaotic variation of σ, so forecasts based on a single power-law fit should be expected to break down at the bifurcation.
Discrete proper-time synchronization offers a simple extrapolation tool: once the exponential t(x) is established, future aftershock times can be estimated until the epoch ends.
Proper-time chronometry reveals spatial order invisible in universal time: foreshocks converge toward the mainshock epicenter and aftershocks diverge from it when plotted against source proper time.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the exponential synchronization is robust, the most important test is catalog completeness: missing early aftershocks or contamination by non-aftershocks would make the numbering x_i = i unfaithful, so the exponential law should be rechecked with completeness-corrected or borehole records.
The same deactivation/activation framework could be extended to foreshocks, where the left branch of g(t) gives a source activation coefficient; a unified triad law for activation and relaxation might follow, though the paper only sketches this.
A testable extension is to apply the underground-clock method to induced seismicity or laboratory acoustic emission, where event detection is nearly complete, to see whether exponential proper-time decay is a universal source property or an artifact of catalog selection.
If proper time effectively stops during background activity, the source clock freezes between earthquake cycles; this could have consequences for recurrence-time models, but the paper leaves that implication unstated.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Colleague,

The paper's central 'law' is a tautology. With their definitions, Eq. (4) follows for any smooth, decreasing aftershock sequence—σ = -n'/n² and τ = ∫σ dt make dn/dτ + n² = 0 an identity, not a discovery. The discrete representation is no stronger: t(x) is fitted as an exponential, then the exponential decay law is read off the same curve. So the headline claim—exponential decay in source proper time—is not independently established.

Credit where it is due. The paper is clearly written and honest enough to label its axioms as postulates, and it explicitly concedes it cannot derive them. The source-centered time idea is a genuinely useful way to think about aftershocks, and the reported two-stage relaxation for Tohoku, with a constant σ in the first month and erratic behavior after, is a pattern worth checking. Figure 4, on foreshock/aftershock distances versus proper time, is a new and thought-provoking visual.

The soft spots are load-bearing, though. There are no uncertainty estimates for σ; R²=0.97 on a two-parameter exponential fit to 4,537 ordered intervals is weak evidence without residual analysis. The catalog completeness problem is serious: if early aftershocks are missing, numbering events by index can make a power-law sequence look geometric. The paper gives no magnitude cutoff, declustering protocol, or completeness period, and the 'more than a hundred events' claim is unsubstantiated. The dismissal of Hirano-Utsu rests on a hand-wavy monotonicity argument, not a statistical test. Much of the mathematical content comes from the authors' previous papers; the new material here is the comparative framing and the distance figure.

Who should read this? Seismologists interested in reparameterizations and anyone who teaches what counts as a physical law. As a claim about earthquake physics, I would not treat it as established. The authors are serious and there may be a real empirical regularity underneath, but this preprint does not make the case.

My recommendation: if the journal can find a referee who will demand catalog completeness analysis, error bars, and an out-of-sample test, it deserves peer review. Otherwise, it is a desk reject with an invitation to resubmit with real data handling.

Referee Report

4 major / 5 minor

Summary. The paper presents four representations of aftershock decay: the classical Omori hyperbolic law, the Hirano–Utsu power law, and two new representations called 'linear' and 'discrete.' In the linear representation the deactivation coefficient is defined by σ = dg/dt with g = 1/n − 1/n0, and the source proper time is τ = ∫σ dt; from these definitions the paper derives dn/dτ + n² = 0. In the discrete representation, aftershocks are numbered x_i = i and fitted by t(x) = t1 exp(σ(x−1)), from which the paper derives n(χ) = n0 e^{−χ} and dn/dχ + n = 0. The paper concludes that Omori's law holds only during an 'Omori epoch,' that the Hirano–Utsu law contradicts observations, and that the holistic law is exponential decay of aftershock activity with respect to the source's proper time. The main empirical support is a single Tohoku aftershock sequence (4,537 events) shown in Fig. 3, plus superposed-epoch analyses in Figs. 2 and 4.

Significance. If established, the claimed exponential relaxation in source proper time would be a substantial reorientation of aftershock physics, replacing power-law fits with a state-variable description of the source and introducing an 'underground clock.' The paper has the virtue of being explicit and internally transparent: the definitions of σ, τ, and χ are given in closed form, and the algebraic steps are easy to check. However, the central equation of the linear representation is an identity following from the definitions, and the discrete exponential law is a restatement of an exponential fit to one unvalidated catalog. The empirical evidence lacks uncertainty estimates, completeness analysis, or an independent test. Thus the paper's contribution is better described as a proposed re-parameterization of aftershock sequences than as an empirically established law.

major comments (4)

[Linear representation, Eqs. (3)–(5)] Equation (4) is not a physical law; it is an identity under the paper's own definitions. With g = 1/n − 1/n0 and σ = dg/dt, the proper time τ = ∫σ dt equals g up to a constant. Therefore dn/dτ = (dn/dt)/(dg/dt) = (−σ n²)/σ = −n² for any smooth function n(t). The text itself notes that Eq. (3) imposes no restrictions on σ(t). Consequently, the derivation of (4) and the 'Omori epoch' behavior in Fig. 1 carry no independent empirical content: σ = const is a smoothing/fitting statement, not a tested law. This is load-bearing because the paper's central claim of a holistic law is presented as following from (3)–(5).
[Discrete representation, Eqs. (8)–(11)] The discrete exponential law is a restatement of the fitted exponential (9). Since σ = d ln t/dx by (8), an exponential t(x) immediately gives inter-event intervals T = dt/dx ∝ t and hence n = 1/T ∝ e^{−σx}; Eq. (11) is then the derivative of the fitted curve. More importantly, x_i = i assumes the catalog lists the complete, unbroken sequence of source excitation events. The paper gives no magnitude cutoff, declustering procedure, completeness period, or uncertainty for the 4,537 Tohoku events in Fig. 3. With t1 = 7.2 h, early missing aftershocks would censor the smallest inter-event intervals and can produce an apparent geometric progression even if the physical sequence follows a power law. The agreement σ = 0.0014 between the discrete and linear fits is therefore not independent confirmation: both values are derived from the same unvalidated event-time series.
[Power representation / Hirano–Utsu rejection] The rejection of the Hirano–Utsu law rests on a qualitative monotonicity argument: from Eq. (2), σ(t) is monotonic for constant p, whereas the paper says the measured σ is non-monotonic in the second stage. This is not a quantitative test. It relies on the same smoothed σ from Fig. 1 and on the claim of 'more than a hundred events,' which is not documented in the manuscript. No uncertainty, formal goodness-of-fit comparison, or handling of the well-known p ≠ 1 deviation is provided. Thus the statement that the Hirano–Utsu representation 'contradicts observations' is not supported by the evidence presented.
[Conclusion / three postulates] The third postulate, 'aftershock activity decays exponentially with the source's proper time,' is the theorem that the paper derives from definitions (3)–(8) and the empirical fit (9), rather than an independent axiom. The paper also concedes that the axiom system is incomplete and that the postulates cannot be derived from general physics. As formulated, the theory is unfalsifiable: σ(t) is defined from the data, and τ (or χ) is constructed so that Eqs. (4) and (11) hold. A testable version would require a priori restrictions on σ(t), a prediction from one catalog tested on another, or a completeness-corrected analysis showing that the exponential form survives when missing early events are accounted for.

minor comments (5)

[Introduction] The reference to 'Utsu, 1861, 1962' contains a typo: the first date should be 1961.
[Appendix] In the final paragraph, 'The experiment showed that 1σ = 1' appears to be a typo; the intended statement is σ ≪ 1.
[Fig. 3] The red dot marking the end of the Omori epoch is not defined by a quantitative criterion; please state the procedure used to set this boundary.
[General] The terms 'holistic' and 'relevant' are used without definitions. They should be replaced or precisely specified if they are meant to carry technical meaning.
[Eqs. (10)–(11)] The notation χ is introduced in (10) but the lower limit of the integral ∫σ dx is not specified; clarify the relationship between χ and the discrete index x.

Circularity Check

3 steps flagged · score 8.0 of 10

The exponential law is built into the definitions: σ and proper time are constructed from n(t), so Eqs. (3)-(4) and (8)-(11) are identities that restate the fitted curves rather than predict them.

self definitional [Linear representation, Eqs. (3)-(4)]
"We postulate the law of aftershock evolution in the form of the simplest differential equation dg/dt=σ, ... It is easy to verify that for σ=const the solution to equation (3) is equivalent to Omori's law (1). However, the dynamic equation (3) is remarkable precisely because, unlike the classical Omori law, it imposes no restrictions on the time dependence of the deactivation coefficient."

The coefficient σ is not independently specified; setting σ(t)=d(1/n−1/n0)/dt makes Eq. (3) true for any observed n(t). Then τ=∫σ dt makes Eq. (4), dn/dτ+n²=0, an algebraic identity for every smooth aftershock sequence. The paper itself says (3) 'imposes no restrictions.' Thus the 'holistic law' contains no dynamical content unless σ=const is imposed empirically; that constancy is a fit, not a consequence of the equation.
fitted input called prediction [Discrete representation, Eqs. (8)-(11) and Fig. 3]
"The set of points is well approximated by the exponential function t(x)=t1 exp(σ(x−1)) ... with t1=7.2 h and σ=0.0014. The coefficient of determination is close to unity: R²=0.97. ... To conclude this section, let us state in words the law of aftershock evolution in a discrete representation: The interval between successive aftershocks increases exponentially over proper time. Accordingly, the frequency of aftershocks decreases exponentially with the passage of proper time: n(χ)=n0 exp(−χ), χ=∫σ dx. (10) ... (11)."

Equation (9) is a fit to the empirical t_i-versus-i points. With σ=d ln t/dx and χ=∫σ dx, Eq. (9) is χ=ln(t/t1); defining the discrete frequency as n=1/T≈1/(dt/dx)=1/(σt) makes n(χ)=n0 e^{−χ} an algebraic rewrite of the same fitted exponential. The exponential relaxation law is therefore not independently predicted or tested: it is the fitted input re-expressed in new variables. The only non-tautological residue is the goodness of fit (R²=0.97).

1 more flagged steps

other [Discrete representation, paragraph after Fig. 3]
"The proper functioning of our special-purpose clock is evidenced by the fact that the deactivation coefficient σ=0.0014 matches the value we obtained from the linear representation of the Tohoku aftershock evolution."

This 'match' is internal consistency, not independent confirmation. Both σ values are derived from the same Tohoku catalog: the discrete σ is the slope of log t versus x in Eq. (9), and the linear σ is d(1/n)/dt. Under the discrete identification n=1/T≈1/(σt), these are the same parameter by construction. Hence the agreement is built into the processing and cannot validate the clock or the law.

full rationale

The paper's central derivation is largely definitional. In the linear representation, σ is introduced via g=1/n−1/n0, so Eq. (3) is an identity and Eq. (4) follows by the definition of proper time; any smooth aftershock sequence satisfies these equations. The empirical content must come from the claim that σ is constant in the Omori epoch and small, which is a fitted observation, not a derivation. In the discrete representation, an exponential curve is fitted to the catalog event index versus time, and then the exponential decay law is read off the same curve after defining χ=∫σ dx. The law is a reformulation of the fit, not a prediction. The paper's additional check that the discrete and linear σ values agree is also internal to the same processing. I do not score the two-stage relaxation observation itself as circular: it is an empirical pattern, though its robustness is weakened by the absence of any documented completeness, declustering, or magnitude-cutoff analysis. The self-citations are numerous, but the central circularity is mathematical/definitional rather than a self-citation chain. Because the headline 'exponential decay with proper time' is forced by the way σ and χ are defined, the score is 8 rather than lower.

Assumptions & free parameters 3 free parameters · 6 assumptions · 3 invented entities

The mathematical content reduces to a change of variable: with σ = -n'/n², τ = ∫σ dt = 1/n - 1/n0, and Eq. (4) is an identity for any smooth n. The empirical laws rest on postulates 1-3 and on a single fitted exponential.

free parameters (3)

σ (deactivation coefficient, Tohoku Omori epoch) = 0.0014
Chosen so that the exponential t(x)=t1 exp(σ(x-1)) fits the Tohoku aftershock times in Fig. 3; the same value is obtained from the smoothed linear representation. It is a per-event empirical parameter, not derived from a theory.
t₁ (first aftershock time offset) = 7.2 h
Free parameter in the exponential fit Eq. (9); sets the origin of the discrete proper-time scale.
Smoothing window / binning choices = hourly bins; 25-point smoothing (Fig. 2); aftershock sequence selection for Fig. 3
Regularization choices in Eq. (7) and Fig. 2 that determine whether σ(t) appears constant during the 'Omori epoch' and chaotic afterward; no sensitivity analysis is given.

assumptions (6)

ad hoc to paper Eq. (3): dg/dt = σ(t) is postulated as the law of aftershock evolution.
§Linear representation: 'We postulate the law of aftershock evolution in the form of the simplest differential equation...' The paper explicitly calls it a postulate.
domain assumption There exists a state function of the source (formalized by σ(t)).
Conclusion: 'There exists a function representing the state of the earthquake source as a dynamic system.'
ad hoc to paper There exists a proper time of the source, a functional of the state function.
Postulate 2 in Conclusion; the paper states it as an axiom.
ad hoc to paper Aftershock activity decays exponentially with source proper time.
Postulate 3 in Conclusion; this is the central claim, stated as an unproved postulate.
domain assumption The aftershock frequency n(t) is a smooth, slowly decaying function of time.
Appendix uses this to conclude σ << 1; it is the 'first postulate of the classical theory.'
domain assumption Catalog completeness and representativeness of the USGS/NEIC events.
The discrete representation requires t_i to be the actual aftershock excitation times; no completeness or declustering discussion is given.

invented entities (3)

Source proper time τ (linear) / χ (discrete)
purpose: A reparametrized time coordinate in which aftershock decay becomes Omori/exponential.
Defined only from the aftershock frequency itself (τ = ∫σ dt = 1/n - 1/n0), so no observable outside the data used to define it.
Underground clock / discrete proper time x
purpose: Measures source time by counting aftershocks (tremors as ticks).
The clock's rate is synchronized using the same aftershock times whose evolution it is supposed to explain; no independent standard.
Deactivation coefficient σ(t)
purpose: State variable describing source relaxation; central parameter in both representations.
Computed from smoothed aftershock frequency; no independent measurement (stress, strain, etc.) is provided.

how reviews work

0 comments

Cite this review

Pith. "Pith review of On four representations of the law of aftershock evolution." pith.science (2026). https://pith.science/paper/HOKGPVFL

@misc{pith2026260729381,
  author       = {Pith},
  title        = {Pith review of: On four representations of the law of aftershock evolution},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HOKGPVFL}},
  note         = {Machine review of arXiv:2607.29381}
}

read the original abstract

The relaxation of the earthquake source following the main shock has been investigated both theoretically and experimentally, based on data regarding the evolution of aftershocks. It is shown that Omori picture of aftershock evolution is not holistic, whereas the Hirano-Utsu representation contradicts observations. Linear and discrete representations are proposed. They are based on original concepts regarding deactivation and the proper time of the source. The difference between the linear and discrete representations lies in the procedure for processing experimental data on aftershocks. In the linear representation, proper time is calculated based on the frequency of aftershocks. In a discrete representation, the source's proper time is measured by an imaginary underground clock, the rate of which differs radically from that of universal time. In the linear representation, the frequency of aftershocks is described by the simplest nonlinear differential equation. In the discrete representation, the aftershock frequency is described by a first-order linear differential equation. It has been found that aftershock activity decays exponentially with the source's proper time. Keywords: earthquake source, foreshocks, mainshock, aftershock theory, Omori law, Utsu's law, proper time, underground clock, deactivation coefficient, source relaxation.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages

[7]

P. 111–200. Rodrigo M.R. A spatio -temporal analogue of the Omori -Utsu law of aftershock sequences // Cornell University Library: arXiv:2111.02955v1 [physics.geo -ph], submitted on 21 Oct 2021, pp. 1-12. https://doi.org/10.48550/arXiv.2111.02955 Utsu T. A statistical study on the occurrence of aftershocks // Geophys. Mag. 1961. V. 30. P. 521–605. Salinas...

work page Pith review arXiv doi:10.48550/arxiv.2111.02955 2021
[1938]

V . 56. P. 111–139. Ogata Y., Zhuang J. Space –time ETAS models and an improved extension // Tectonophysics. 2006. V. 413, Iss. 1–2. P. 13-23. Omori F . On the aftershocks of ear thquake // J. Coll. Sci. Imp. Univ. Tokyo. 1894. V

2006

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.

[1] [7]

P. 111–200. Rodrigo M.R. A spatio -temporal analogue of the Omori -Utsu law of aftershock sequences // Cornell University Library: arXiv:2111.02955v1 [physics.geo -ph], submitted on 21 Oct 2021, pp. 1-12. https://doi.org/10.48550/arXiv.2111.02955 Utsu T. A statistical study on the occurrence of aftershocks // Geophys. Mag. 1961. V. 30. P. 521–605. Salinas...

work page Pith review arXiv doi:10.48550/arxiv.2111.02955 2021

[2] [1938]

V . 56. P. 111–139. Ogata Y., Zhuang J. Space –time ETAS models and an improved extension // Tectonophysics. 2006. V. 413, Iss. 1–2. P. 13-23. Omori F . On the aftershocks of ear thquake // J. Coll. Sci. Imp. Univ. Tokyo. 1894. V

2006