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

The central idea behind diffusional earthquakes is that fault slip can be driven not by the abrupt rupture of frictional patches but by the gradual, pressure-driven migration of fluids and stress through the fault zone. In this framework, a fault zone is treated as a porous, deformable medium in which pore pressure evolves by diffusion. When a perturbation in pore pressure or stress travels along the fault, it can progressively weaken the fault surface and generate slip as it moves. The result is an event that looks, in many respects, like an earthquake, but whose physics is governed by diffusion rather than by the elastic waves that dominate ordinary seismic rupture.

One of the most consequential outcomes of the study is the reported scaling relationship between slip and distance. In classical earthquake physics, the average slip on a fault surface is related to the rupture length, and this relationship has been used for decades to estimate the magnitude of ancient earthquakes from the size of exposed fault scarps. Diffusional earthquakes, according to the new work, follow their own scaling in which slip accumulates as the disturbance propagates. That means that unlike ordinary ruptures, which tend to show slip that depends strongly on the total area of the ruptured patch, diffusional events appear to involve slip that grows systematically with the distance the disturbance has traveled along the fault.

This slip-distance relationship matters because it provides a diagnostic signature. If scientists can measure how slip accumulates with distance in observed fault-slip events, they can, in principle, distinguish between processes driven by elastic rupture and those governed by diffusion. That distinction is not merely academic. Diffusional slip may be associated with swarms of small earthquakes, with slow-slip episodes detected by geodesy, and with creep events on faults that never generate destructive shaking. Understanding which scaling law a given event obeys could therefore tell researchers something fundamental about the mechanism beneath the observed motion.

Fluids occupy a central role in this emerging picture. Pore fluids within a fault zone carry pressure that both reduces the effective normal stress clamping the fault shut and transports stress perturbations through the medium. Because the transport is diffusive, the characteristic timescales and length scales of the resulting slip depend on the hydraulic diffusivity of the fault-zone rock. Highly permeable damage zones can transmit pressure changes over considerable distances in relatively short times, while low-permeability regions trap fluids and localize deformation. In the diffusional-earthquake model, these hydraulic properties effectively set the pace of fault motion, replacing the control that frictional instability dynamics exert in conventional seismicity.

The mathematical treatment presented in the study blends continuum poromechanics with fault-friction concepts. Rather than treating the fault as a simple sliding surface, the approach accounts for the coupled evolution of deformation, fluid pressure and slip along the fault zone. The analysis shows that a self-sustaining front of slip can emerge, and that the scaling of slip with distance follows directly from the balance between the elastic response of the surrounding rock and the diffusive transport of pore pressure through the fault. In effect, the geometry of slip is imprinted by diffusion, and the authors demonstrate that the resulting scaling differs in a testable way from the linear crack-like scaling familiar from classical seismology.

The implications reach into several active areas of earthquake science. Seismic swarms in geothermal fields, volcanic regions and fluid-injection sites have long been suspected to involve migrating fluid pressure, and models of pressure diffusion are routinely used to explain how such swarms expand over days to months. A well-defined slip-distance scaling for diffusional slip gives researchers a quantitative tool for these interpretations. If the cumulative slip of migrating events can be inferred from geodetic measurements or from seismological observations, the scaling law can be checked directly against field data, turning what has been a qualitative association between fluids and swarm seismicity into a testable quantitative prediction.

The framework may also speak to slow-slip phenomena observed along subduction zones and other major fault systems. Slow-slip events release strain over hours to weeks, far slower than ordinary earthquakes, and they often migrate along strike at rates that some researchers have compared to diffusion. Whether the migration of slow slip is controlled by dilatancy, by viscous rheologies or by fluid diffusion remains debated, but a scaling law derived from diffusive physics offers a way to discriminate among these possibilities. Events that follow the diffusional scaling reported in the study would point toward pore-pressure transport as the controlling process, while departures from that scaling would implicate other mechanisms.

There are also practical consequences for hazard assessment. Fault motion that is driven by diffusion tends to be slower and gentler than ordinary rupture, but it can still load adjacent locked patches of fault that may eventually fail seismically. Estimating how much slip accumulates in diffusional events, and over what distances, therefore contributes to a fuller accounting of how strain is redistributed through fault networks. The new scaling provides a compact way to make such estimates, potentially improving models of how slow, fluid-driven slip feeds into the earthquake cycle on larger faults.

As with any new theoretical framework, the test now is confrontation with observations. Field campaigns that combine dense seismic arrays, geodetic networks and measurements of fault-zone hydraulic properties will be needed to verify the predicted slip-distance relationship in natural settings. Laboratory experiments on fluid-saturated fault gouge could offer another route to testing the scaling under controlled conditions. If the diffusional scaling holds up, it will give seismologists a second, distinct signature to look for in fault-slip data, alongside the classical scaling of ordinary earthquakes, and a sharper lens through which to view the hidden, fluid-mediated processes that shape the earthquake cycle.

Subject of Research: Diffusional fault-slip events driven by fluid diffusion and their scaling of slip with propagation distance.

Article Title: Diffusional earthquakes and their slip-distance scaling

Article References: Sato, D. S., & Yoshida, K. (2026). Diffusional earthquakes and their slip-distance scaling. Communications Earth & Environment. https://doi.org/10.1038/s43247-026-04043-4

Image Credits: AI Generated

DOI: 10.1038/s43247-026-04043-4

Keywords: diffusional earthquakes, fault slip, slip-distance scaling, pore pressure diffusion, fluid-driven seismicity, slow slip, earthquake swarms, poromechanics, seismology, fault mechanics, Diffusional, earthquakes

The appeal of the roughness hypothesis has always been its visual plausibility. On bathymetric maps, the regions that hosted the largest recorded events, such as the 2004 Sumatra-Andaman earthquake, the 2011 Tohoku-oki earthquake and the 1960 Chile earthquake, seemed to correspond with segments of the trench where the incoming plate appeared comparatively smooth, while rougher margins like those off parts of Central America or Japan’s older, sediment-starved trenches appeared to host only smaller ruptures. The hypothesis gained traction because it offered a forecasting shortcut: map the bumps on the incoming seafloor before it subducts, and you have a proxy for how big the next earthquake might be. Given how difficult it is to observe the actual fault surface at depth, a measurable feature of the incoming plate seemed like a gift to hazard assessors.

Yet the new analysis finds that the correlation does not survive careful, systematic testing. Rather than selecting a few celebrated case studies, the researchers assembled a globally consistent dataset of subduction interface roughness measurements and paired them with the maximum magnitudes that each margin has produced, constrained by both the instrumental record and, where available, longer historical and paleoseismic evidence. When roughness is quantified uniformly, using the same topographic measures of relief and spectral character of the incoming plate across all margins, the expected relationship between rough interfaces and smaller maximum earthquakes largely dissolves. Margins with prominently rough seafloor have still hosted very large ruptures, and some comparatively smooth margins have not delivered the giant events the hypothesis predicts.

The authors point to several reasons why the roughness argument fails as a control on maximum magnitude. First, roughness measured on the seafloor before subduction is a poor guide to roughness on the fault at seismogenic depths. As the plate descends, sediments smear into the topographic lows, hydrothermal alteration and mineralization smooth the interface, and the accretionary wedge redistributes material. A seamount that looks like a formidable asperity at the trench may be largely underplated or subducted within a weaker sedimentary layer by the time it reaches the depth range where most seismic moment is released. The fault that ruptures in a great earthquake is therefore not the same rough surface that geodesists and marine geologists can map from shipboard sonar.

Second, and more fundamentally, the maximum magnitude of an earthquake is a geometric property of how large a patch can rupture in a single event, and that is governed by the lateral and downdip extent of coherent locking, the segmentation of the plate boundary, and the accumulated slip deficit, not necessarily by centimeter- to kilometer-scale frictional heterogeneity. A rupture can jump or creep past small roughness elements if the surrounding fault is sufficiently stressed and strongly coupled. Conversely, a smooth interface segmented by major structural boundaries such as fracture zones or tear faults may still be unable to host a rupture longer than a few hundred kilometers. In other words, the features that actually arrest or release ruptures may be structural discontinuities with tens of kilometers of offset, not the relief on the incoming plate.

The study also re-examines the physical reasoning behind the roughness hypothesis itself. Laboratory and theoretical work shows that roughness influences frictional behavior most strongly at small scales, affecting the onset of slip and the distribution of aftershocks, but its effect on the total possible rupture area at the scale of a magnitude 9 earthquake is weak. Scale considerations matter: an earthquake of magnitude 9 ruptures a fault patch on the order of 1,000 kilometers long. Bumps a few kilometers across are simply too small to act as persistent barriers to a rupture front carrying enormous elastic strain energy. The comparison is often made to tearing a sheet of paper: small wrinkles in the paper do not determine where the tear stops if the sheet is being pulled hard enough.

The implications for seismic hazard assessment are significant and uncomfortable. If roughness cannot be used as a proxy for maximum magnitude, then several regional hazard models that factor seafloor roughness into estimates of maximum credible earthquakes may need revision. The study suggests that hazard assessors should weight other evidence more heavily, including geodetic measurements of plate coupling, the distribution of past ruptures inferred from historical accounts, tsunami deposits and coral microatolls, and the structural segmentation of the margin. It also argues against using roughness to declare any margin incapable of producing a giant earthquake, a conclusion with direct consequences for coastal communities and infrastructure planning along subduction margins worldwide.

The researchers emphasize that their findings do not make roughness irrelevant. Interface roughness still shapes where within a rupture the largest slip occurs, how strong ground shaking is distributed, and possibly the frequency of smaller-to-moderate events. What the study removes is the assumption that roughness imposes a hard ceiling on rupture size. Distinguishing between influences on the distribution of slip and controls on maximum magnitude is, the authors argue, an essential refinement that the field has too often glossed over. Roughness may sculpt the earthquake, but it does not appear to cap it.

The work also speaks to a broader lesson in earthquake science: the danger of building predictive frameworks on a handful of spectacular case studies. The great earthquakes of the past century are few, and any global comparison involving fewer than a dozen giant events is statistically fragile. By compiling a comprehensive, uniformly processed dataset, the new study reduces the risk of overfitting a narrative to memorable examples. The result is a more sober picture of subduction zones: nearly any of them, given enough time to accumulate strain, may be capable of producing ruptures larger than their instrumental records suggest, and the size of the next great earthquake may be far less predictable from the shape of the seafloor than scientists once hoped.

Subject of Research: The relationship between subduction interface roughness and maximum earthquake magnitude

Article Title: Maximum earthquake magnitude unlikely to be controlled by subduction interface roughness

Article References: Yang, X., Bell, R. E., Whittaker, A. C., Xu, H., Han, X., Knowlson, A. R., & Locher, V. A. (2026). Maximum earthquake magnitude unlikely to be controlled by subduction interface roughness. Nature Geoscience. https://doi.org/10.1038/s41561-026-02093-z

Image Credits: AI Generated

DOI: 10.1038/s41561-026-02093-z

Keywords: subduction zones, earthquake magnitude, fault roughness, seismic hazard, plate boundaries, rupture segmentation, great earthquakes, seismology, geodesy, marine geophysics, plate coupling, Nature Geoscience

Nejla Essaddi, Imen Ziadi, and Mongi Besbes, affiliated with the University of Tunis El Manar, SUP’COM, and the Higher Institute of Information and Communication Technologies at the University of Carthage, set out to answer a question that has hovered over the quantum machine learning community for a decade: does quantum computing offer any genuine advantage for practical classification problems, or is it all theoretical promise? Rather than testing on toy datasets, the team turned to the Mediterranean, one of the most tectonically complicated regions on Earth, where the African and Eurasian plates collide in a patchwork of subduction zones, strike-slip faults, and continental collision. They drew their data from the United States Geological Survey catalog, a publicly accessible record of global seismicity, and framed the task as a binary classification problem: does an event exceed local magnitude 4.0, the threshold above which earthquakes begin to pose genuine hazards?

The scale of the benchmark is what distinguishes this study from much of the quantum machine learning literature. The researchers evaluated 261 distinct experimental configurations, systematically varying the temporal windows of seismic features fed into each model, the dimensionality of those feature vectors, and the choice of classifier. On the classical side, they tested logistic regression, the workhorse of interpretable statistics; random forests, ensembles of decision trees renowned for robustness; and support vector machines, which separate data classes by finding optimal boundaries in transformed feature spaces. On the quantum side, they deployed two fundamentally different architectures: the variational quantum classifier, a hybrid quantum-classical circuit whose parameters are tuned by a classical optimizer, and the quantum kernel support vector machine, which computes distances between data points in an exponentially large quantum feature space and feeds those similarities to a classical SVM.

The technical distinction between these two quantum approaches matters enormously for interpreting the results. Variational quantum classifiers function like quantum neural networks: data is encoded into qubit states through parameterized rotation gates, a measurement produces a prediction, and the parameters are iteratively adjusted to minimize a loss function. Quantum kernel methods, by contrast, bypass training of the quantum circuit altogether. Instead, each pair of data points is loaded into a quantum circuit whose evolution depends on the data values, and the overlap between the resulting quantum states—a quantity requiring exponentially many classical operations to compute exactly—serves as a kernel in a classical support vector machine. When the encoding circuit is chosen well, this kernel can capture patterns that are classically hard to represent, which is precisely the kind of advantage quantum machine learning theorists have been hunting for.

The headline finding is nuanced but striking. Across the full Mediterranean dataset, the classical random forest remained the most dependable global performer, reaching accuracies as high as 0.833. But in specific tectonic settings—particular zones combined with particular temporal histories of seismic activity—the entangled quantum kernel SVM achieved dramatically better results. Its peak accuracy reached 0.947, a 28.5 percent relative improvement over the best classical result in that same zone-history combination, achieved with only six qubits. In an era when quantum hardware is noisy, limited, and expensive, the fact that a six-qubit model could dominate a tuned random forest on any slice of real geophysical data is a genuinely remarkable result.

The authors, however, are careful not to oversell the finding, and their statistical honesty is part of what makes the study valuable. That peak performance did not reach statistical significance in their pre-specified fixed-model comparison; after applying the Holm–Bonferroni correction for multiple comparisons, the adjusted p-value was 1.0. The reasons are practical rather than mysterious: test sets in narrow tectonic zones are small, and the standout quantum result emerged from post-hoc model selection—examining many configurations and highlighting the best one—which inherently inflates apparent performance. The team explicitly frames these results as exploratory, a signal worth pursuing rather than a settled proof of quantum advantage. This kind of methodological transparency is rare in a field often criticized for hype, and it sets a benchmark for how quantum machine learning claims should be reported.

The study also delivered a cautionary tale about variational quantum classifiers. While quantum kernel methods thrived in favorable geometries, the VQC models suffered severe training instability as qubit counts increased—a phenomenon well known in the quantum computing literature as the barren plateau problem, where gradients of the loss function vanish exponentially with circuit size, leaving optimizers wandering a nearly flat landscape with no useful direction. McClean and colleagues first characterized this pathology in 2018, and the Mediterranean earthquake data confirms it in practice: pushing the variational approach to higher qubit counts degraded rather than improved results. Intriguingly, however, the VQC showed unexpected niche competitiveness in data-scarce regimes, hinting that different quantum architectures may suit different data availability conditions—a finding with real implications for seismically active but poorly instrumented regions.

What does it mean for a quantum feature space to align with geophysical data geometry? The Mediterranean’s seismotectonic zones produce feature distributions shaped by interacting fault systems, depth-dependent attenuation, and regional magnitude scales. When the embedding circuit entangles features in a way that mirrors these physical correlations, the quantum kernel can draw decision boundaries that classical kernels approximate only crudely. The 0.947 accuracy in specific configurations suggests that, at least locally, the structure of earthquake catalog features resonates with the expressivity of a modest quantum circuit. The authors argue this points toward targeted, hybrid quantum-classical early warning systems, in which classical models handle most of the workload and quantum kernels are deployed selectively where their strengths apply.

The practical roadmap emerging from this research is one of selective hybridization rather than wholesale replacement. Earthquake early warning is a domain where seconds matter and false alarms carry real economic and social costs, so any improvement in classification reliability, even in narrow regimes, is consequential. The Tunisian team’s work demonstrates that quantum kernels deserve a place in the toolkit—not as a universal solution, but as a specialized instrument whose deployment should be guided by the geometry of the data and the tectonic character of the region. As quantum hardware matures and qubit counts grow beyond the six used here, the boundary between classical and quantum competitiveness will inevitably shift. For now, the Mediterranean has provided the testing ground where quantum machine learning took a measurable, statistically honest step from theory toward the seismic frontier.

Subject of Research: Benchmarking quantum kernel machine learning against classical classifiers for earthquake magnitude classification in the Mediterranean region.

Article Title: Harnessing quantum kernels for robust earthquake classification: a Mediterranean case study

Article References: Essaddi, N., Ziadi, I., & Besbes, M. (2026). Harnessing quantum kernels for robust earthquake classification: a Mediterranean case study. Quantum Machine Intelligence, 8(2), Article 102. https://doi.org/10.1007/s42484-026-00443-z

Image Credits: AI Generated

DOI: 10.1007/s42484-026-00443-z

Keywords: quantum machine learning, earthquake classification, quantum kernel SVM, variational quantum classifier, Mediterranean seismicity, USGS catalog, random forest, support vector machine, early warning systems, seismology, six-qubit circuits, Holm-Bonferroni correction