REVIEW 4 major objections 6 minor 21 references
Relativity of Time in Earthquake Physics
T0 review · 4 major / 6 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read Two foreshock phases of the Tohoku earthquake emerge when events are counted as the source's own clock.
desk verdict A small empirical case study that re-labels event order as 'proper time' and finds a two-phase foreshock pattern in Tohoku, but the phase split is fixed by the M7.3 and the statistics are too weak to prove a regime change. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central mechanism is the 'underground clock': each catalogued shock (foreshock or aftershock) is treated as one unit of the source's proper time x, so the clock ticks at every event. The deactivation coefficient sigma is recovered from the data through the relation sigma = d ln T / dx, where T is the mean interval between successive events in universal time. In the phenomenological theory, sigma is the coefficient in the hyperbolic decay equation dn/dt + sigma n^2 = 0, so that the standard aftershock decay law corresponds to constant sigma; the clock's purpose is to make this comparison possible without assuming any universal-time law in advance.
What would settle it
Take the same foreshock catalog and recompute the deactivation coefficient using equal windows of world time instead of event number, or randomize the event order and re-run the fit; if the straight-line/exponential split at event 104 disappears or the coefficients change widely across shuffles, the claimed two phases depend on the chosen clock rather than on source physics.
Extended reading notes
Core claim
On its own terms, the paper establishes a two-phase evolution of the Tohoku foreshocks. Numbering foreshocks by x = 1..166 and plotting their occurrence times t(x) gives a straight line t = 46.5x - 485 for x = 1..103 (R^2 = 0.98) and an exponential trend exp(0.065 x) for x = 104..166 (R^2 = 0.97). Equivalently, the deactivation coefficient sigma = d ln T / dx is zero in the first phase and 0.065 in the second. The authors interpret this as a sudden transition from a degenerate constant-rate regime to a classical hyperbolic decay regime, and note the similarity to the aftershock pattern: the initial aftershock phase also has sigma = 0 before it jumps to a positive constant. The proper time of
Load-bearing premise
The entire analysis assumes that simply numbering successive earthquakes creates a physical clock, so that each event marks one equal unit of the source's proper time; if the event number is not a meaningful time variable, the two phases are an artifact of how the data are reordered.
Editorial extensions
If this is right
- The Tohoku foreshock sequence is not stationary: it switches abruptly from a constant-rate regime (sigma = 0) to an exponentially slowing regime (sigma = 0.065) at the 104th foreshock.
- The decelerating phase likely reflects a cascade of aftershocks triggered by the magnitude-7.3 foreshock 51 hours before the main shock, meaning the main shock was preceded by a detectable secondary sequence.
- The same ordering by event number reveals an analog in aftershock relaxation, where sigma also starts at zero and then jumps to a positive constant, suggesting a common regime-change pattern in source evolution.
- The uneven flow of proper time relative to world time gives a direct experimental measure of source non-stationarity: sigma(t) is not constant when expressed in world time.
- Applying the procedure to other 'complete classical triad' earthquakes is the next planned test of whether this two-phase foreshock structure is generic.
Reading between the lines
- If the pattern reproduces on other events, the event-count clock would provide a simple, parameter-free time axis for comparing foreshock and aftershock evolution across different magnitudes and tectonic settings.
- The transition at x = 104 could be detected in real time if the event count and magnitude are known, offering a possible short-term precursor indicator: a switch from steady to hyperbolic slowing before a large event.
- The exponential segment is fit with a single coefficient; a cleaner falsification is to test on shuffled event order or to compare the exponential fit against a two-parameter power law on the same interval.
- The vortex, or Umov-energy, reading of the epicenter path is speculative, but it implies that foreshock epicenters should trace a systematically curved path in space-time, which is measurable.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a phenomenological concept of 'proper time' for an earthquake source, defined operationally as the ordinal number x of events in a catalog ('underground clock'). For the 2011 Tohoku earthquake, the authors plot world-time occurrence times t(x) for 166 foreshocks and identify two phases: a linear relation t(x)=46.5x−485 for x=1–103 (Eq. 5) and an exponential relation t(x)∝exp(0.065x) for x=104–166 (Eq. 6). They report deactivation coefficients σ=0 and σ=0.065 for the two phases (Table 1) and compare these with earlier aftershock phases. The paper also presents preliminary spatial analyses (foreshock convergence, aftershock divergence, and a 'vortex movement' interpretation). The central claim is that ordering events by this proper time reveals a previously unnoticed regime change in foreshock evolution.
Significance. If the two-phase foreshock structure is real, it would complement the authors' earlier three-phase aftershock relaxation picture and could indicate a sudden change in source deactivation about 51 hours before the main shock. The paper is clearly written, uses a public catalog, provides explicit functional fits, and honestly acknowledges that the result may be specific to this earthquake. The analogy to relativity is explicitly stated to be symbolic, avoiding overclaiming on that front. However, the statistical support for the central regime-change claim is currently weak: the phase split is imposed at the time of the known M7.3 foreshock rather than estimated, no competing model is tested, and the reported σ values are derived from the same fits. The paper is better viewed as a suggestive phenomenological description than as a statistically validated discovery.
major comments (4)
- [Section 2, Eqs. (5)–(6) and Table 1] The two-phase claim is not statistically validated. The breakpoint is fixed at the time of the M7.3 foreshock (x≈103/104) rather than estimated from the data, and the paper does not compare the two-phase model against a single-phase alternative (e.g., one exponential or a power-law acceleration). Because t(x) is a cumulative sequence of waiting times, consecutive points are strongly autocorrelated, so high R² values (0.98 and 0.97) are expected for many smooth curves and do not by themselves establish a regime change. The authors should estimate the breakpoint (or at least test sensitivity to its location), provide confidence intervals for the slopes/rates, and use a model-comparison criterion such as AIC or a likelihood-ratio test.
- [Section 2, Eq. (7) and Table 1] The deactivation coefficients in Table 1 are not independent measurements; they are direct consequences of the fits in Eqs. (5) and (6). For the first phase, t(x) is linear, so T=dt/dx is constant and Eq. (7) gives σ=0 by construction. For the second phase, t(x)∝exp(0.065x), so T∝exp(0.065x) and σ=d ln T/dx=0.065 exactly. Thus Table 1 restates the phase fits rather than confirming them. The paper should be explicit that σ is a fitted parameter of the piecewise model, not an independently estimated quantity.
- [Section 2, definition of proper time] The 'underground clock' is defined so that each event increments proper time by one unit. Under this definition, the statement that 'proper time flows unevenly relative to world time' is a restatement of the fact that inter-event times vary. The physical content of the paper therefore reduces to the empirical claim that the inter-event time pattern changes at x≈103/104. This may be a legitimate empirical finding, but the manuscript should not present it as a measurement of an independent source clock without justifying that ordinal numbering carries physical meaning. A concrete test would be to apply the same procedure to other earthquake sequences and show that the piecewise structure is reproducible and preferable to a single-phase null model.
- [Section 3, Figs. 6 and 7] The spatial results (foreshock convergence, aftershock divergence, 'vortex movement') are presented without quantitative uncertainty or significance testing. The moving average over 20 points is not accompanied by confidence bounds, and the curvature/torsion of the epicenter path in Fig. 7 is not measured. Since these claims are secondary to the main temporal phase analysis, they could be relegated to clearly labeled preliminary observations, or supported with appropriate statistical measures.
minor comments (6)
- [References, Ref. [6]] The journal title contains a typo: 'Solid Evarth' should be 'Solid Earth'.
- [Section 2, Eq. (1) and Fig. 3] The symbol τ is used in Eq. (1), but the main text and figures use x for proper time. Please define the relation between τ and x explicitly.
- [Figures 4 and 5 captions] Fig. 4 refers to a green dot marking 'a distinct change' but the caption does not specify its coordinates or how it was determined. Please clarify whether this dot is the M7.3 foreshock.
- [Section 2, Eq. (5)] The intercept -485 in Eq. (5) is dimensionally ambiguous. Please specify the reference time (e.g., days before the main shock) and the time units.
- [Section 2, data description] The catalog source is given as USGS/NEIC, but the retrieval date and version are not stated. This would aid reproducibility.
- [Section 3, Fig. 6] The moving-average window is described as '20 points' in the text, but it would be helpful to state whether this is a centered or trailing window and whether the smoothing changes results for different window lengths.
Circularity Check
No significant circularity; the two-phase result is an empirical fit, not a prediction derived from the theory.
full rationale
The paper's derivation chain is observational: it defines proper time x as the event ordinal, plots the catalog occurrence times t_k against x, fits a straight line and an exponential to different segments, and defines the deactivation coefficient sigma from T = dt/dx via Eq. (7). Nothing in this chain is equivalent to its inputs by construction: the clock definition (event ordinal as proper time) does not force t(x) to be piecewise linear/exponential, and the fitted slopes are not predicted from Eqs. (1)–(4); they are estimated from the USGS/NEIC catalog. The Table 1 sigma values are transformations of the fitted curves, but the paper does not present them as an independent prediction from theory. The split at x=103/104 is chosen from the known M7.3 foreshock, which is a statistical/validation weakness (no breakpoint estimation, no model comparison), not a circularity. Self-citations to [1,2,6,7] supply the phenomenological framework and prior aftershock results, but the key equations are stated and applied directly rather than imported as an unexamined black box. Therefore no circular step is present.
Assumptions & free parameters
free parameters (6)
- Slope of first foreshock phase (a1) =
46.5 (per event, units unspecified)
- Exponential growth rate of second foreshock phase (beta) =
0.065 per event
- Deactivation coefficient for first phase (sigma1) =
0
- Deactivation coefficient for second phase (sigma2) =
0.065
- Aftershock deactivation coefficients =
0 and 0.0014
- Moving-average window length =
20 points
assumptions (4)
- domain assumption The aftershock (and foreshock) frequency n is related to the deactivation coefficient sigma by the differential equation dn/dt + sigma*n^2 = 0 (eq. 2).
- domain assumption The source can be described as a dynamic system with a well-defined deactivation coefficient sigma(t).
- ad hoc to paper Each event in the catalog serves as a tick of the source's proper time, i.e., the ordinal number x is a uniform measure of proper time.
- domain assumption The deactivation coefficient can be estimated from the logarithm of the inter-event time gradient: sigma = d ln T / dx (eq. 7).
invented entities (4)
-
Proper time of the earthquake source (tau or x)
-
Hypothetical underground clock
-
Deactivation coefficient sigma
-
Umov energy flow and vortex movement
Cite this review
Pith. "Pith review of Relativity of Time in Earthquake Physics." pith.science (2026). https://pith.science/paper/5MOHKORT
@misc{pith2026250904858,
author = {Pith},
title = {Pith review of: Relativity of Time in Earthquake Physics},
year = {2026},
howpublished = {\url{https://pith.science/paper/5MOHKORT}},
note = {Machine review of arXiv:2509.04858}
}
read the original abstract
In this paper we develop an understanding of the proper time of the Tohoku earthquake source. The paper is dedicated to the 120th anniversary of Einstein's theory of relativity, but the dedication is symbolic, since we are investigating a purely non-relativistic geophysical object. We still found it possible to borrow the terms time relativity and proper time from the theory of relativity. The concept of the proper time of the source is developed by us within the framework of the phenomenological theory of earthquakes. The paper describes a procedure for measuring proper time using observation data of foreshocks and aftershocks. We synchronized the imaginary underground clock that counted the proper time and the clock that showed world time. The idea of a hypothetical underground clock turned out to be effective. We have shown that the proper time of the source flows unevenly relative to the flow of world time. Two phases of the evolution of the source before the main shock of the earthquake were discovered. This complements the picture of three-phase relaxation of the source after the main shock, which we established earlier. The uneven flow of proper time relative to world time indicates the non-stationarity of the source parameters. The overall conclusion is that the concept of the proper time of the source has enriched the possibilities of experimental study of earthquakes.
Reference graph
Works this paper leans on
-
[1]
On the elementary theory of earthquake aftershocks
Guglielmi A.V., Zotov O.D. On the elementary theory of earthquake aftershocks // arXiv:2505.10414 [physics.geo-ph]
-
[2]
Zotov O.D. Guglielmi A.V . Evolution of the Tohoku earthquake aftershocks in the framework of the phenomenological theory of aftershocks // arXiv:2508.04323 [physics.geo-ph]
-
[3]
On he Omori Law in the physics of earthquakes // Appl
Zavyalov A., Zotov O., Guglielmi A., Klain B. On he Omori Law in the physics of earthquakes // Appl. Sci. 2022, vol. 12, issue 19, 9965. https://doi.org/10.3390/app12199965. 11
-
[4]
The fundamentals of a phenomenological theory of earthquakes // J
Guglielmi, A.V ., Klain, B.I., Zavyalov, A.D., Zotov, O.D. The fundamentals of a phenomenological theory of earthquakes // J. V olcanol. Seismol., 2023, vol. 17, no. 5, pp. 428–438
work page 2023
-
[5]
Guglielmia A.V ., Zavyalov A.D., Zotov O.D., Klain B.I. A method for processing and analysis of aftershocks due to a tectonic earthquake: A new look at an old problem // J. V olcanology and Seismology. 2025. V . 19. No. 2. P . 196–202
work page 2025
-
[6]
Interpretation of the Omori Law
Guglielmi A.V. Interpretation of the Omori Law // Izv. Phys. Solid Evarth. 2016. V . 52. P . 785–786 // arXiv:1604.07017 [physics.geo-ph]
work page Pith review arXiv 2016
-
[7]
Omori's law: a note on the history of geophysics // Phys
Guglielmi A.V. Omori's law: a note on the history of geophysics // Phys. Usp. 2017. V. 60. P . 319–324. DOI: 10.3367/UFNe.2017.01.038039
-
[8]
On the aftershocks of earthquake // J
Omori F. On the aftershocks of earthquake // J. Coll. Sci. Imp. Univ. Tokyo. 1894. V . 7. P. 111–200
Show all 21 references
-
[9]
Investigation of aftershocks of the Great Kanto earthquake at Kumagaya // Kishoshushi
Hirano R. Investigation of aftershocks of the Great Kanto earthquake at Kumagaya // Kishoshushi. Ser. 2. 1924. V . 2. P. 77–83 (in Japanese)
1924
-
[10]
Aftershocks and periodicity in earthquakes // Gerlands Beitr
Jeffreys H. Aftershocks and periodicity in earthquakes // Gerlands Beitr. Geophys. 1938. V
1938
-
[11]
The centenary of the Omori formula for a decay law of aftershock activity // J
Utsu T., Ogata Y., Matsu’ura R.S. The centenary of the Omori formula for a decay law of aftershock activity // J. Phys. Earth. 1995. V . 43. № 1. P. 1–33. http://dx.doi.org/10.4294/jpe1952.43.1
1995 doi
-
[12]
Space–time ETAS models and an improved extension // Tectonophysics
Ogata Y., Zhuang J. Space–time ETAS models and an improved extension // Tectonophysics
-
[13]
Lagrangian formulation of Omori’s law and analogy with the cosmic Big Rip, Eur
Faraoni V . Lagrangian formulation of Omori’s law and analogy with the cosmic Big Rip, Eur. Phys. J. C., 2020, 80(5): 445. DOI:10.1140/epjc/s10052-020-8019-2
2020 doi
-
[14]
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
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. 12
-
[15]
On the possibility of reproducing Utsu’s Law for earthquakes with a spring-block SOC model // Entropy
Salinas-Martínez A., Perez-Oregon J., Aguilar-Molina A.M., Muñoz-Diosdado A., Angulo- Brown F . On the possibility of reproducing Utsu’s Law for earthquakes with a spring-block SOC model // Entropy. 2023. V . 25. P. 816. https://doi.org/10.3390/e25050816
2023 doi
-
[16]
Instanton Representation of Foreshoсk— Aftershock Sequences // Izv
Bogomolov L.M., Rodkin M.V ., Sychev V.N. Instanton Representation of Foreshoсk— Aftershock Sequences // Izv. Phys. Solid Evarth. 2025, № 2. P . 43–57. DOI: 10.31857/S0002333725020045
2025 doi
-
[17]
Dependence of the source deactivation factor on the earthquake magnitude // arXiv:2108.02438 [physics.geo-ph]
Guglielmi A.V., Zotov O.D. Dependence of the source deactivation factor on the earthquake magnitude // arXiv:2108.02438 [physics.geo-ph]
-
[18]
On the three laws of earthquake physics // Journal of V olcanology and Seismology
Guglielmi A.V., Zavyalov A.D., Zotov O.D., Klain B.I. On the three laws of earthquake physics // Journal of V olcanology and Seismology. 2025, vol. 19, # 5 // Cornell University Library: arXiv:2506.06504
2025 arXiv
-
[19]
On the spatial-temporal structure of aftershock sequences // In: Yanovskaya T
Zotov O.D, Zavyalov A.D., Klain B.I. On the spatial-temporal structure of aftershock sequences // In: Yanovskaya T. et al. (eds). Problems of Geocosmos–2018, Springer Proceedings in Earth and Environmental Sciences. Springer, Cham. 2020. P . 199. DOI: 10.1007/978-3-030-21788-4_16
2018 doi
-
[20]
Seismicity localization before large Kamchatkan earthquakes // Volcanology and Seismology
Zavyalov A.D., Nikitin Yu.V. Seismicity localization before large Kamchatkan earthquakes // Volcanology and Seismology. 2000, vol.21, pp.525-534
2000
-
[2006]
413, Iss
V . 413, Iss. 1–2. P. 13-23
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.