Seismic Compression-Quiescence Anomalies as Earthquake Precursors: Prospective Prediction for the Kashmir Himalayan Fault Belt and Validation Across Four Tectonic Regimes

1. Probleemstelling

De Kashmir–Himalaya-boog vormt een van de meest seismisch gevaarlijke regio's wereldwijd, gedreven door de noordwaartse convergentie van de Indiase en Euraziatische platen. Deze convergentie slaat aanzienlijke elastische rekenergie op langs parallelle breuksystemen (Main Frontal Thrust, Main Boundary Thrust, Panjal Thrust en Main Central Thrust). Ondanks historische catastrofes zoals de M 7.6 Kashmir-aardbeving in 2005, kampen huidige methoden voor aardbevingvoorspelling in deze regio met drie specifieke beperkingen:

• Het gebrek aan een enkele kwantitatieve drempelwaarde die geldig is in verschillende tektonische settings.
• De afwezigheid van prospectief verifieerbare voorspellingsgegevens specifiek voor Kashmir.
• Het onvermogen om belastingssignalen toe te schrijven aan specifieke breuklijnen in plaats van aan brede geografische regio's.

Bestaande methoden, zoals seismische quiescentie (Z-statistiek), b-waarde mapping en versnellende momentvrijgave, hebben deze doelen niet simultaan bereikt.

2. Methodologie

De studie introduceert een datagestuurd kader gecentreerd rond de Magnitude Compressie Index ($C$) en het C-drop signaal, gevalideerd over vier verschillende tektonische regimes (Kashmir, Japan, Turkije, Italië).

Kernmetrieken:

• Magnitude Compressie Index ($C$): Een parameter-vrije maatstaf gedefinieerd als $C(k) = \max(0, 1 - \sigma[W] / \sigma_0)$. Het kwantificeert de vernauwing van de magnitudeverdelingen van aardbevingen naarmate een breuk spanning accumuleert. Naarmate de stressbelasting de lokale seismiciteit homogeniseert, neemt de standaarddeviatie ($\sigma$) van de magnitudes af, waardoor $C$ stijgt richting 1.
• De C-drop: Het primaire precursore signaal. Dit wordt gedefinieerd als een systematische, aanhoudende post-piek daling in $C$ (minstens 0,050 vanaf een piek $\ge$ 0,25). Deze daling treedt op wanneer de breuk begint te falen, waardoor grotere aardbevingen in het observatievenster terechtkomen en de magnitude-variantie toeneemt.
• Seismische Quiescentie Index ($Q$): Een secundaire maatstaf voor de reductie in frequentie ($Q = \max(0, 1 - \lambda[W] / \lambda_0)$). Verhoogde $Q$ gecombineerd met verhoogde $C$ duidt op een kritisch vergrendelde breuk.

Gegevensverwerking & Zone-identificatie:

• Rolling Window: Berekeningen maken gebruik van een rolling window van de $W=20$ meest recente aardbevingen in een zone, uitgesmeerd door een 12-gebeurtenissen bewegend gemiddelde.
• Baseline: Referentiewaarden ($\sigma_0, \lambda_0$) worden afgeleid van de periode 2000–2009.
• Data-gestuurde Zonering: In tegenstelling tot eerdere studies die willekeurige coördinaten gebruikten, zijn de centrumpunten van de zones geplaatst op natuurlijke pieken van seismische dichtheid, geïdentificeerd via Gaussische smoothing van de ISC-catalogus (2000–2026). De radii van de zones vangen 90% van de lokale seismische clusters op.
• Zes-criteria Score-systeem: Om fout-positieven te filteren (bijv. naschok-sequenties), wordt een zone alleen als KRITIEK geclassificeerd als deze gelijktijdig aan ten minste 3 van de 6 criteria voldoet:
  1. Compressie: $C \ge 0,25$.
  2. Quiescentie: $Q \ge 0,20$.
  3. C-drop: Bevestigde daling vanaf de piek.
  4. b-waarde: $b < 0,75$ (wijst op stressconcentratie).
  5. Lage Activiteit: $<15$ gebeurtenissen in de laatste 90 dagen.
  6. Versheid: De piek vond plaats binnen de laatste 180 dagen.

3. Belangrijkste Bijdragen

• Het C-drop Signaal: Vestigt de post-piek daling in magnitude-compressie als een primair, fysiek onderbouwd precursor voor naderende significante aardbevingen, onderscheidend van eenvoudige frequentie-gebaseerde quiescentie.
• Data-gestuurde Zone-decompositie: Demonstreert dat samengestelde seismische zones kunnen worden gedecomposeerd om specifieke actieve breuklijnen te identificeren. Zo isoleert de studie het actieve belastingssignaal in de Anantnag-regio tot de Panjal Thrust (Ramban–Banihal) in plaats van de ontlaste Main Central Thrust kern.
• Prospectieve Verificatie: Het artikel presenteert een formele voorspelling gearchiveerd op Zenodo (18 juni 2026) voordat de activerende gebeurtenissen plaatsvonden, wat een verifieerbaar record biedt van de toepassing van de methode.
• Cross-tektonische Validatie: De methode wordt retrospectief getest over vier verschillende plaatgrenstypen (botsing, subductie, strike-slip en compressieve boog), wat de universaliteit van het C-drop mechanisme aantoont.

4. Resultaten

Retrospectieve Validatie:
Over 53 bevestigde C-drop cycli in vier tektonische regimes trad een aardbeving met een doel-magnitude op binnen 90 dagen in 83,0% van de gevallen (44/53).

• Japan (Subductie): 88,9% detectiegraad.
• Turkije (Transform): 91,7% detectiegraad.
• Italië (Compressie): 75,0% detectiegraad.
• India/Kashmir (Botsing): 73,3% detectiegraad.
• High-Confidence Subset: Voor cycli met $C_{peak} \ge 0,30$ en een drop $\ge 0,055$ bereikte de detectiegraad 100% over alle regio's.

Status van de Kashmir-zone (per 20 juni 2026):

• Zone Z21 (Anantnag): Score 5/6. $C_{drop} = 0,129$. Twee M 4.8 gebeurtenissen vonden plaats op 11 april en 5 juni 2026, tijdens de actieve C-drop fase. De cyclus blijft open met $C=0,260$ (67% van de piekbelasting behouden).
• Zone Z15b (Noord-Kashmir): Score 5/6. $C_{drop} = 0,092$. Hoogste belasting in 26 jaar ($C_{peak}=0,585$). Historische detectiegraad voor M $\ge$ 4,5 is 96,3%. Voorspelling M 5,0–6,0 voor juli–september 2026.
• Zone Z15 (Muzaffarabad): Score 4/6. $C=0,342$ met uitzonderlijk diepe quiescentie ($Q=0,795$). Bevestiging van C-drop (0,050 drempel) wordt verwacht binnen 2–4 weken.
• Zone FZ-NPB (Gilgit): Score 4/6. Hoogste $C$-waarde in de studie (0,530). De C-drop is net begonnen; een opkomende voorspelling is gemarkeerd voor augustus–oktober 2026.

Validatiegebeurtenissen:
Drie gebeurtenissen valideerden het prospectieve voorspellingsvenster:

1. M 3.3 (19 juni 2026): Vond 18 uur na het Zenodo-archief plaats, na een stilte van 13 dagen in Zone Z21.
2. M 4.8 (11 april 2026) & M 4.8 (5 juni 2026): Beiden vonden plaats tijdens de actieve C-drop fase van Zone Z21, wat staat voor gedeeltelijke spanningsontlading.

5. Betekenis en Claims

Het artikel claimt een prospectieve, prospectief verifieerbare aardbevingvoorspelling te bieden voor de Kashmir Himalaya breukbelt, waarmee een kritiek gat in de huidige seismologische praktijk wordt gedicht.

• Fysisch Mechanisme: De auteurs stellen dat de C-drop een universeel breukmechanisme reflecteert waarbij de magnitude-variantie afneemt onder belasting en toeneemt bij falen, onafhankelijk van de tektonische setting.
• Precisie: Door gebruik te maken van data-gestuurde zone-decompositie verfijnt de studie de ruimtelijke attributie van spanning, waarbij de Panjal Thrust wordt geïdentificeerd als de actieve bron voor het Anantnag-signaal, in plaats van de voorheen aangenomen Main Central Thrust.
• Risicobeoordeling: De studie benadrukt dat ongeveer 10 miljoen mensen in de getroffen zones leven in Jammu en Kashmir, door Pakistan beheerd Kashmir en Gilgit-Baltistan.
• Oproep tot Actie: De auteurs nodigen de internationale seismologische gemeenschap (inclus van NCS en CSEP) uit om de voorspelling onafhankelijk te monitoren en te evalueren, waarbij zij benadrukken dat de methode zelfconsistent is in het onderscheiden van geladen (voorspelbare) en ontlaste (niet-voorspelbare) breuktoestanden.

Het artikel concludeert dat de C-drop methode, gecombineerd met data-gestuurde zonering en een zes-criteria filter, een robuust kader biedt voor het identificeren van kritieke breukbelastingsstaten en het uitgeven van prospectieve voorspellingen met een hoge statistische betrouwbaarheid.