Between 2006 and 2017, North Korea carried out six underground nuclear tests at the Punggye-ri facility, which lies beneath Mt. Mantap in the country’s northeast. The final and largest of these explosions, conducted in September 2017, had an estimated yield of 100 to 250 kilotons. The United States Geological Survey registered the event as a magnitude 6.3 earthquake, and satellite radar measurements revealed substantial deformation of Mantap’s summit, evidence that the blast had physically warped the rock above the test cavity. For years afterward, most scientific attention focused on the standard questions of nuclear test monitoring: where the explosion occurred, when it happened, and how large it was.

The new research, led by Xingli Fan and colleagues, takes a broader view. Rather than characterizing individual explosions, the team set out to understand how the tests changed the seismic behavior of the region as a whole. To do so, the researchers analyzed seismic data recorded in China and South Korea since 2008, drawing on stations located between 80 and 200 kilometers from the test site. That distance matters: instruments far from the epicenter record the region’s background seismicity rather than just the violent, short-lived signals of the detonations themselves, allowing the team to track subtle changes in earthquake activity over nearly two decades.

The scale of what they found surprised even the authors. By combing through the continuous data streams, Fan and colleagues identified 1,399 local earthquakes around Mt. Mantap between 2008 and 2025, far more than previous earthquake catalogs had recorded. Many of these events were small enough to escape earlier detections, but together they paint a picture of a crust that has been fundamentally unsettled by the testing program. The sheer number of events suggests that the mountain and its surroundings experienced a level of seismic disturbance that conventional monitoring approaches had substantially underestimated.

What makes the sequence remarkable is its timing. After most underground nuclear explosions, seismologists observe a familiar pattern: a burst of aftershocks in the immediate vicinity of the blast cavity that decays rapidly over days or weeks, much like the aftershock sequences that follow natural earthquakes. Mt. Mantap refused to follow the script. Instead of decaying, seismic activity after the September 2017 test began roughly three weeks after the detonation and then continued to grow in both frequency and magnitude through 2025, intensifying years after the last explosion. This is not the signature of a crust settling back into equilibrium; it is the signature of an ongoing process.

High-precision relocation of the earthquakes revealed an equally striking organizational pattern. The events were not scattered randomly through the rock but concentrated along two roughly north-northwest–trending fault structures, some of which had existed before the testing began while others had gone unrecognized until now. That alignment indicates persistent and organized fault reactivation rather than the chaotic shattering one might expect from simple blast damage. The faults, in other words, appear to have been switched on as coherent structures, slipping repeatedly along their length in the years following the detonations.

The physical explanation proposed by the researchers is a gradual one. According to the study, the repeated nuclear explosions progressively damaged the shallow crust around Mt. Mantap and altered its internal stress field. In a region where many faults already sat close to failure, even modest perturbations to the stress balance could tip them into sliding. Rather than failing all at once, the faults became active gradually over several years as stresses redistributed through the damaged rock, producing the slow escalation in seismicity that the team documented. The 2017 explosion, with its enormous yield and the pronounced ground deformation it caused, likely delivered the decisive push.

The implications extend well beyond a single mountain on the Korean Peninsula. The study demonstrates that underground nuclear explosions have the potential to reactivate faults that were previously seismically quiet, awakening geological structures that no monitoring program would have flagged as hazardous. Earthquakes triggered in this way may be difficult to distinguish from naturally occurring tectonic activity, since they occur on real faults and can continue for years with no obvious temporal link to the explosion that caused them. For organizations tasked with verifying compliance with nuclear test bans, that ambiguity is a serious concern: a cluster of earthquakes near a former test site could represent natural tectonics, lingering explosion effects, or something else entirely.

Monitoring former nuclear test sites, the authors suggest, will require longer horizons and finer tools than have typically been applied. Dense regional seismic networks capable of detecting and precisely locating small events, combined with sustained observation over many years, appear essential to capturing delayed sequences like the one at Mt. Mantap. The study’s catalog of 1,399 earthquakes, assembled from stations hundreds of kilometers away, shows what such sustained analysis can reveal even at considerable distance from the source. As more nations conduct and then abandon underground nuclear testing, the geological legacies they leave behind may keep generating earthquakes, and scientific attention, long after the political headlines have moved on. For now, Mt. Mantap stands as the clearest demonstration yet that the consequences of nuclear testing are written not only in treaties and diplomacy but in the slow, patient mechanics of the Earth’s crust itself.

Subject of Research: Delayed earthquake activity caused by nuclear explosion–induced fault reactivation at North Korea's Mt. Mantap test site

Article Title: Mt. Mantap nuclear testing triggered delayed seismicity via fault reactivation

Article References: Mt. Mantap nuclear testing triggered delayed seismicity via fault reactivation. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: Mt. Mantap, Punggye-ri, nuclear test monitoring, fault reactivation, delayed seismicity, induced earthquakes, seismology, North Korea, intraplate faults, crustal stress, Science journal, nuclear explosions

Injection-induced earthquakes, or IIEs, have become a growing concern in regions where oil and gas operations pump large volumes of fluid into the subsurface. Most induced events are small, but some have reached magnitudes capable of causing damage, prompting jurisdictions to adopt operational safeguards. The most common of these is the traffic-light system, which calls for injection rates to be reduced or halted altogether when earthquakes of increasing size occur near an operation. The logic behind such systems rests on two assumptions: that smaller foreshocks can serve as warning signs of a larger event to come, and that stopping injection quickly enough can prevent that larger event from occurring. Yet until now, scientists have lacked a clear explanation for why some induced earthquakes are preceded by swarms of foreshocks while others arrive with little or no advance notice, and how foreshock behavior connects to the physical processes that culminate in a mainshock.

The stakes of this gap in understanding are not abstract. In Western Canada, a surge of induced seismicity over the past decade has drawn sustained attention from regulators, operators, and researchers alike. Events felt by residents near communities such as Fox Creek in Alberta drew international scientific scrutiny after being linked to hydraulic fracturing operations, and similar concerns have arisen in other jurisdictions where wastewater disposal or stimulation activities have coincided with earthquake swarms. The economic and social consequences are tangible: operational shutdowns carry significant costs for industry, while residents living near active wells face repeated shaking, concerns about structural damage, and questions about long-term risk. A more reliable way to anticipate when small events might escalate into damaging ones would therefore carry value well beyond the scientific community.

To address those questions, Bei Wang and colleagues turned to a decade of high-quality seismic recordings from western Canada, a region that has experienced a well-documented rise in induced seismicity tied to fluid injection. From that dataset, the team examined foreshock activity preceding 77 IIEs of magnitude 3 or greater. The scale and quality of the records allowed the researchers to detect and characterize foreshocks that would have gone unnoticed in less sensitive monitoring, providing an unusually detailed statistical foundation for analyzing how foreshock sequences develop in space and time before larger ruptures.

The results were striking in their consistency. Roughly 92% of the earthquakes in the study — 71 of the 77 events — were preceded by identifiable foreshocks. In other words, detectable precursory seismicity was the norm rather than the exception for injection-induced earthquakes of significant magnitude in this region. That finding carries immediate practical weight, because it supports the basic premise underlying traffic-light protocols: that smaller earthquakes accompanying injection operations are, far more often than not, part of a build-up toward a larger event rather than isolated, self-limiting rumblings.

The near-ubiquity of foreshocks also reframes how small events during injection should be interpreted. Under many current protocols, a magnitude threshold triggers a response — caution, reduced injection, or shutdown — but the underlying assumption has often been that most small earthquakes will simply die out on their own. The new results suggest that in this region, at least, small events during injection are statistically likely to be early expressions of a developing rupture sequence. That does not mean every foreshock presages a damaging mainshock, but it does mean the base rates that inform protocol design should be revisited with this high foreshock-precedence figure in mind.

But the study went beyond simply counting foreshocks. Wang and colleagues found that foreshock productivity — how many foreshocks occur, and how vigorous the sequences are — along with their spatial and temporal patterns, reflected the interplay of three key factors: the fluid injection itself, the seismogenic index of the affected volume of rock, and the stress state of the faults involved. In essence, the character of foreshock activity encodes information about how close a fault is to failure, how susceptible the surrounding rock is to seismic rupture, and how the injected fluid is perturbing the system. That interplay helps explain why foreshock behavior varies from one induced earthquake to another: different combinations of injection history, rock properties, and fault stress conditions produce different precursory signatures.

The seismogenic index, a concept drawn from statistical seismology, is worth unpacking for readers unfamiliar with it. In broad terms, it captures how prone a given volume of rock is to producing seismicity in response to a given amount of fluid-induced perturbation. Two sites receiving similar injection volumes can behave very differently: one may respond with abundant small earthquakes, while the other remains nearly silent because its faults are oriented unfavorably or its stress conditions are far from failure. By tying foreshock productivity to this index, the study provides a physical rationale for why precursory activity is so variable — and why a warning threshold tuned in one basin may misfire in another.

Perhaps the most consequential contribution of the work is the identification of three distinct rupture nucleation models describing how injection-induced earthquakes get underway. The first is a fluid-driven cascade, in which injected fluid drives a chain of small seismic events that progressively leads to a larger rupture. The second is fluid-driven pre-slip with an intact source asperity, in which fluids promote aseismic — silent, non-radiating — slip on parts of a fault while a locked patch, or asperity, remains intact until it finally breaks in the mainshock. The third is injection-weakened pre-slip, in which fluid injection weakens the fault sufficiently that pre-slip develops and ultimately cascades into rupture. Together, these models illustrate the diverse ways in which fluids can promote aseismic slip and transfer stress through the fault system before a mainshock strikes.

The distinction among these models matters because each implies a different warning signature. A fluid-driven cascade should announce itself through an accelerating sequence of detectable small earthquakes clustered near the injection point. Pre-slip models, by contrast, can unfold largely beneath the detection threshold of seismic networks, with foreshocks appearing only sporadically as the silent slip loads stress onto locked patches. An injection-weakened fault may show a more gradual drift toward instability. Recognizing which model best describes a developing sequence — and the study shows that more than one pathway can operate even within a single region — is a step toward interpreting real-time monitoring data with physical rather than purely statistical rules.

The recognition that fluids can drive aseismic processes in the run-up to an induced earthquake is particularly important. Foreshocks recorded at the surface are only the audible portion of a broader deformation process; silent slip, accelerated by fluid pressure changes, can load stress onto locked fault patches without generating any seismic signal at all. This means that seismic monitoring alone may provide an incomplete picture of the state of a fault during injection operations. The authors accordingly conclude that risk mitigation for induced earthquakes should be informed by combined seismic and geodetic monitoring — pairing earthquake detection with measurements of ground deformation that can reveal aseismic slip in progress.

Geodetic techniques such as satellite-based radar interferometry and continuous GPS stations can detect millimeter-to-centimeter scale deformation of the ground surface, deformation that often accompanies slow slip on faults at depth. Integrating such measurements with dense seismic arrays would give operators and regulators a two-channel view of fault behavior: the seismic channel capturing radiating events, and the geodetic channel capturing the silent loading that may precede them. The challenge, as the authors acknowledge, is that geodetic data are not routinely collected or ingested in real time by most current monitoring programs, making this recommendation as much an operational roadmap as a scientific conclusion.

The implications for traffic-light protocols are significant. The finding that more than 90% of significant induced earthquakes in western Canada had foreshocks validates the use of small-event thresholds as warning triggers. At the same time, the study shows that foreshock productivity and patterns vary according to the seismogenic index and fault stress state, meaning that a uniform, one-size-fits-all threshold may perform differently from one region, formation, or fault to another. A traffic-light scheme calibrated for one geological setting could under-warn in another, or impose costly operational shutdowns where precursory activity does not actually signal a large rupture. The authors argue that IIE risk mitigation should instead be tailored to regional and local conditions.

The study also carries implications for fundamental earthquake science. Foreshocks preceding natural tectonic earthquakes remain an active and contested research area, with debate over whether they reflect a genuine nucleation process or are simply part of a self-similar earthquake cascade. Injection-induced sequences offer a valuable natural laboratory for these questions because the timing and location of the forcing — the injection — are comparatively well constrained. By classifying foreshock sequences into three physically distinct nucleation models, the study demonstrates that multiple pathways to rupture can operate even within a single region and dataset, underscoring the complexity of how earthquakes begin.

Several caveats accompany the findings. The analysis draws on a decade of recordings from western Canada, and the extent to which the 92% foreshock-precedence figure and the three nucleation models apply to other regions with different geology, injection practices, or monitoring networks remains to be established. Detecting foreshocks also depends heavily on the sensitivity and density of the seismic network; regions with sparser instrumentation would likely identify a lower fraction of events with precursors, so cross-regional comparisons must be made carefully. In addition, while the study links foreshock patterns to the seismogenic index and fault stress state, translating those insights into operational decision rules for individual injection sites will require further work, including real-time integration of geodetic data that many current monitoring programs do not routinely collect.

Nevertheless, the research marks a meaningful advance in the science of induced seismicity. By demonstrating that foreshocks are nearly ubiquitous before significant injection-induced earthquakes in western Canada, and by articulating concrete physical models for how fluids shepherd faults toward rupture, Wang and colleagues have provided both a diagnostic framework and a practical directive. The combination of seismic and geodetic monitoring, tailored to the seismogenic character of each region and fault, offers a path toward traffic-light systems that are more predictive and less blunt — a development with clear stakes for communities and industries operating where the subsurface is increasingly called upon to absorb vast quantities of injected fluid.