📄 earth_science

A Statistical Investigation of Four-Point Relay Seismicity Patterns Across Four Regional Datasets

This study utilizes a statistical framework and iterative traceback algorithm on 20,693 earthquake records from 1960–2026 to demonstrate that four-point relay seismicity chains occur with a significantly non-random 41.30% completion probability within a 168-hour window, revealing distinct regional susceptibilities and energy attenuation patterns that offer new insights into stress transfer and intermediate-term seismic hazard assessment.

Original authors: Chenghao Zhu

Published 2026-08-04

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Chenghao Zhu

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the Earth's crust not as a solid, unyielding rock, but as a giant, slightly wobbly Jell-O mold sitting on a table. When you poke it in one spot, the whole thing wiggles. Sometimes, that wiggle is so strong it triggers a wobble in a completely different spot. In the world of earthquake science, this is called "triggering." Usually, scientists study how one big quake shakes up its immediate neighborhood, creating a swarm of smaller "aftershocks" nearby. But what if the Earth's Jell-O is connected in a way that a poke in Japan could eventually cause a wiggle in Thailand, which then nudges a spot in Myanmar, which finally sets off a tremor in the Himalayas? This idea of a "relay race" of earthquakes, where stress hops from one distant region to another like a baton, is the mystery this study tries to solve.

Why does this matter? Because if we can spot the pattern of these distant hand-offs, we might get a few days' warning before a big event happens in a place that seems quiet right now. It's like noticing that the first three runners in a relay race have already started sprinting; you might guess the fourth runner is about to take off, even if you can't see them yet. This isn't about predicting the exact second an earthquake will happen, but rather understanding if the Earth has a secret rhythm where trouble in one place makes trouble in another place more likely a few days later.

The Great Earthquake Relay Race

In this study, researcher Chenghao Zhu decided to play detective with a massive pile of earthquake data. They looked at over 20,000 earthquakes that happened between 1960 and 2026 across four very different corners of the world: the Northwest Pacific (near Japan), the Taiwan region, Southeast Asia (around Myanmar), and Central/South Asia (the Hindu Kush mountains).

The goal was to see if these four regions were playing a game of "four-point relay." The rules of this game were simple but strict:

An earthquake happens in Region A.
Within a certain number of hours, an earthquake happens in Region B.
Then, one in Region C.
Finally, one in Region D.

If this happened in the right order, it was a "chain." The big question was: Is this just a random coincidence, like flipping a coin and getting heads four times in a row? Or is there a real, physical connection where the first three quakes are actually "priming" the fourth one?

The Magic Number: 7 Days

To find the answer, the researcher had to figure out the perfect time limit for the race. How long does the baton stay in the air? If the time window is too short (like 24 hours), you might miss the connection. If it's too long (like a month), you might accidentally link two earthquakes that have nothing to do with each other.

They tested windows ranging from 1 day up to a week and a half. The results were surprisingly clear. The "sweet spot" turned out to be exactly 168 hours, or 7 days.

When they looked at the data from the year 2000 onward (when earthquake detectors got much better and more reliable), they found something fascinating. If three regions had already had their turn in the relay race, there was a 41.30% chance that the fourth region would join the party within that 7-day window.

To put that in perspective, if these earthquakes were just random noise, you'd only expect this to happen about 15% of the time. The fact that it happened nearly 41% of the time suggests that these regions are indeed talking to each other. It's like if you saw three friends high-five each other, and then the fourth friend high-fived them too, not just by chance, but because they were all part of the same group hug.

The Quality of the Evidence

The study also noticed that the quality of the data mattered a lot. When they looked at older records from before the year 2000, the "relay" seemed to happen less often. Why? Because older earthquake detectors weren't as good at hearing the smaller, quieter quakes. It's like trying to hear a whisper in a noisy room; if you miss the whispers, you can't hear the conversation. Once they switched to the high-quality data from 2000 to 2026, the signal got stronger, and the probability of the chain completing went up. This tells us that the pattern is real, but we needed better ears to hear it clearly.

Who Wins and Who Loses?

Not all regions were equally eager to join the relay. The study found that the region in Southeast Asia (Location 3) was the most likely to be the fourth runner, joining the chain about 25.75% of the time. On the other hand, the Central/South Asia region (Location 4) was the most reluctant, joining only about 17.39% of the time. It seems some parts of the Earth's crust are just more "connected" to the others than others.

The Energy Fizzle

One of the most interesting findings was about the size of the earthquakes. The researchers compared the size (magnitude) of the final earthquake in the chain to the average size of the first three. They found that the final earthquake was almost always slightly smaller—about 0.04 to 0.09 units smaller on the magnitude scale.

Think of it like a game of "telephone" where the message gets a little quieter with each person who repeats it. Or imagine a line of dominoes: the first few might be big and heavy, but by the time the energy reaches the last one, it's a little weaker. This suggests that as the stress travels from one region to the next, some of that energy is lost along the way, perhaps due to friction or the way the rock absorbs the shock.

What This Means for Us

So, what does a 41% chance and a 7-day window actually mean for us? The author is careful to say this isn't a crystal ball. It doesn't mean that if you see three quakes, the fourth one will happen. It just means that the odds are much higher than they would be by pure luck.

This discovery gives scientists a new tool for "intermediate-term" hazard assessment. It suggests that if three distant regions have been active in the last week, the fourth region might be in a state of higher stress and could be more likely to shake soon. It's a statistical clue, a hint that the Earth's crust is a connected system where a ripple in one ocean can eventually reach the shore of another. While we can't stop the earthquakes, understanding these relay patterns helps us understand the hidden conversations happening deep underground.

Technical Summary: A Statistical Investigation of Four-Point Relay Seismicity Patterns

Problem Statement
While traditional seismological research has extensively characterized mainshock-aftershock sequences and doublet earthquakes, the possibility of extended, cascading seismic sequences involving multiple distinct geographic regions remains under-explored. Conventional models, such as the Epidemic-Type Aftershock Sequence (ETAS), typically emphasize magnitude-dependent triggering within localized proximity. This study addresses the gap in understanding regional-scale stress transfer that may manifest as sequential activation across distant tectonic provinces. The core problem is to determine whether "four-point relay seismicity chains"—where events in four distinct regions occur in a strict temporal cascade—exist with statistical significance above random expectations, and to quantify the temporal and magnitude characteristics of such patterns.

Methodology
The study employs a statistical framework utilizing 20,693 earthquake records (magnitude ≥ 4.5) from 1960 to 2026 across four distinct tectonic regions: the Northwest Pacific (Japan-Kuril Trench), the Taiwan region, Southeast Asia (Myanmar Arc), and Central/South Asia (Hindu Kush-Pamir).

Algorithm: An iterative traceback algorithm was developed to identify chains. For a designated final event ( $E_4$ ) in a specific region, the algorithm searches backward in time for an event ( $E_3$ ) in any of the other three regions within a variable time window ( $\Delta T$ ). It then recursively searches for $E_2$ and $E_1$ in the remaining untested regions. To ensure conservative estimates, if multiple candidates exist within the window, the temporally earliest event is selected.
Metric: The study calculates the conditional probability ( $P_{4|3}$ ) of chain completion: the ratio of four-point chains to three-point chains (where the three-point chain excludes the final response region).
Optimization: A two-stage scanning process (coarse: 24–168 hours; fine: 150–186 hours) was used to identify the optimal $\Delta T$ .
Quality Control: Analyses were stratified by data quality, comparing the full dataset (1960–2026) against subsets starting from 1980 and 2000 to account for improvements in global seismic monitoring.
Magnitude Analysis: For completed chains, the relative magnitude ( $\Delta M$ ) was computed as the difference between the final event's magnitude and the average magnitude of the preceding three events.

Key Contributions

Framework Definition: The paper introduces a specific statistical definition for "four-point relay seismicity chains" that relies on temporal succession across predetermined regions without requiring spatial proximity assumptions, establishing a baseline for detecting remote correlations.
Data Quality Sensitivity: The study systematically demonstrates that catalog completeness significantly impacts detection rates, showing that earlier, incomplete data dilutes rather than enhances the signal of these patterns.
Quantification of Regional Susceptibility: The research provides location-specific susceptibility metrics, identifying distinct variations in how different regions act as the final link in a relay chain.
Magnitude Attenuation Evidence: The study quantifies a consistent trend of magnitude attenuation along relay paths, offering a statistical basis for energy decay in cascading seismic events.

Results

Optimal Time Window: The conditional probability of chain completion increases monotonically with the time window, peaking at 168 hours (7 days). At this optimum, the post-2000 dataset yields a completion probability of 41.30%. This rate drops at longer windows (e.g., 39.90% at 186 hours), confirming a well-defined optimum rather than a saturation effect.
Data Quality Impact: The probability of completion increases systematically with data quality: 20.65% (all data), 21.13% (≥1980), and 22.11% (≥2000). This indicates that the signal is robust but was previously obscured by catalog incompleteness.
Regional Variability:
- Location 3 (Southeast Asia) exhibits the highest susceptibility as the fourth point, with a completion probability of 25.75% in the post-2000 dataset.
- Location 4 (Central/South Asia) shows the lowest susceptibility (17.39%), suggesting relative isolation from the inter-regional triggering network.
- Location 2 (Taiwan region) displays a unique characteristic: its mean relative magnitude is near zero (0.00 ± 0.45), indicating that events terminating in this region do not experience the same magnitude attenuation observed in other locations.
Magnitude Attenuation: Across all configurations, the mean relative magnitude ( $\Delta M$ ) is negative, ranging from -0.04 to -0.09 magnitude units. This suggests that response events are systematically smaller than the average of the preceding chain events.

Significance and Claims
The paper claims that the observed 41.30% completion probability at the 7-day window substantially exceeds the estimated random baseline of ~15%, providing strong statistical evidence for non-random temporal clustering across these regions.

The author posits that this 7-day timescale is consistent with physical mechanisms such as post-seismic relaxation (stress diffusion via viscoelastic or fluid-mediated processes), slow slip events, and regional stress field adjustments. The study concludes that these findings offer a quantitative basis for intermediate-term seismic hazard assessment. Specifically, the author suggests that if three regions demonstrate active seismicity within a one-week period, there is an elevated probability of activation in the fourth region relative to baseline rates.

The paper maintains a modest tone regarding its predictive power, explicitly stating that the 41.30% figure is a statistical estimate of chain closure frequency, not a forecasting probability for individual sequences. The author emphasizes that these results serve as a benchmark for evaluating physical models of regional stress transfer and suggest that the method could be generalized to other datasets, potentially refined by incorporating spatial distance constraints.