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

1. समस्या विवरण (Problem Statement)

कश्मीर-हिमालयन आर्क (Kashmir–Himalayan arc) वैश्विक स्तर पर सबसे अधिक भूकंपीय खतरों वाले क्षेत्रों में से एक है, जो भारतीय और यूरेशियाई प्लेटों के उत्तर की ओर अभिसरण (convergence) द्वारा संचालित है। यह अभिसरण समानांतर भ्रंश प्रणालियों (मुख्य अग्रवर्ती थ्रस्ट, मुख्य सीमावर्ती थ्रस्ट, पंजल थ्रस्ट और मुख्य केंद्रीय थ्रस्ट) के साथ महत्वपूर्ण लोचदार विकृति ऊर्जा (elastic strain energy) को संचित करता है। 2005 के M 7.6 कश्मीर भूकंप जैसी ऐतिहासिक आपदाओं के बावजूद, इस क्षेत्र में भूकंप भविष्यवाणी के वर्तमान तरीके तीन विशिष्ट सीमाओं का सामना करते हैं:

• विभिन्न विवर्तनिक (tectonic) परिवेशों में मान्य एकल मात्रात्मक दहलीज (quantitative threshold) का अभाव।
• विशेष रूप से कश्मीर के लिए भविष्योन्मुखी रूप से सत्यापन योग्य (prospectively verifiable) भविष्यवाणी रिकॉर्ड की अनुपलब्धता।
• व्यापक भौगोलिक क्षेत्रों के बजाय विशिष्ट भ्रंश धागों (fault strands) को लोडिंग संकेतों के रूप में जोड़ने में असमर्थता।

मौजूदा विधियाँ, जैसे कि भूकंपीय शांति (Z-statistic), b-वैल्यू मैपिंग, और त्वरित मोमेंट रिलीज (accelerating moment release), इन लक्ष्यों को एक साथ प्राप्त करने में विफल रही हैं।

2. कार्यप्रणाली (Methodology)

यह अध्ययन मैग्नीट्यूड कंप्रेशन इंडेक्स ($C$) और C-ड्रॉप (C-drop) सिग्नल पर केंद्रित एक डेटा-संचालित ढांचे को पेश करता है, जिसे चार अलग-अलग विवर्तनिक क्षेत्रों (कश्मीर, जापान, तुर्की, इटली) में सत्यापित किया गया है।

मुख्य मेट्रिक्स:

• मैग्नीट्यूड कंप्रेशन इंडेक्स ($C$): एक पैरामीटर-मुक्त माप जिसे $C(k) = \max(0, 1 - \sigma[W] / \sigma_0)$ के रूप में परिभाषित किया गया है। यह भूकंपीय परिमाण वितरण के संकुचन को मापता है क्योंकि एक भ्रंश तनाव संचित करता है। जैसे-जैसे तनाव लोडिंग स्थानीय भूकंपीयता को समरूप (homogenize) बनाता है, परिमाण का मानक विचलन ($\sigma$) कम हो जाता है, जिससे $C$ बढ़कर 1 की ओर बढ़ता है।
• C-ड्रॉप (C-drop): प्राथमिक पूर्ववर्ती संकेत (precursory signal)। इसे एक व्यवस्थित, निरंतर पीक के बाद होने वाली गिरावट के रूप में परिभाषित किया गया है (एक पीक $\ge$ 0.25 से कम से कम 0.050 की गिरावट)। यह गिरावट तब होती है जब भ्रंश विफल होने लगता है, जिससे बड़े भूकंप अवलोकन विंडो में प्रवेश कर पाते हैं और परिमाण भिन्नता (magnitude variance) बढ़ जाती है।
• भूकंपीय शांति सूचकांक ($Q$): एक माध्यमिक माप जो दर में कमी को दर्शाता है ($Q = \max(0, 1 - \lambda[W] / \lambda_0$)। उच्च $Q$ के साथ उच्च $C$ एक गंभीर रूप से लॉक (critically locked) हुए भ्रंश का संकेत देता है।

डेटा प्रोसेसिंग और ज़ोन पहचान:

• रोलिंग विंडो: गणनाएँ एक ज़ोन में $W=20$ सबसे हालिया भूकंपों के रोलिंग विंडो का उपयोग करती हैं, जिसे 12-इवेंट मूविंग एवरेज द्वारा सुचारू (smooth) किया जाता है।
• बेसलाइन: संदर्भ मान ($\sigma_0, \lambda_0$) 2000-2009 की अवधि से प्राप्त किए गए हैं।
• डेटा-संचालित ज़ोनिंग: मनमाने निर्देशांकों का उपयोग करने वाले पिछले अध्ययनों के विपरीत, ज़ोन केंद्रों को ISC कैटलॉग (2000-2026) के गाऊसी स्मूथिंग (Gaussian smoothing) के माध्यम से पहचाने गए प्राकृतिक भूकंपीय घनत्व शिखरों पर रखा गया है। ज़ोन त्रिज्या स्थानीय भूकंपीय क्लस्टर्स के 90% हिस्से को कवर करती है।
• छह-मानदंड स्कोरिंग प्रणाली: गलत सकारात्मक परिणामों (जैसे कि आफ्टरशॉक अनुक्रम) को फ़िल्टर करने के लिए, एक ज़ोन को केवल तभी क्रिटिकल (CRITICAL) वर्गीकृत किया जाता है यदि वह एक साथ 6 में से कम से कम 3 मानदंडों को पूरा करता है:
  1. कंप्रेशन: $C \ge 0.25$।
  2. क्विएसेंस (शांति): $Q \ge 0.20$।
  3. C-ड्रॉप: पीक से पुष्टि की गई गिरावट।
  4. b-वैलue: $b < 0.75$ (तनाव संकेंद्रण का संकेत)।
  5. कम गतिविधि: पिछले 90 दिनों में $<15$ घटनाएं।
  6. ताजगी (Freshness): पीक पिछले 180 दिनों के भीतर हुआ हो।

3. प्रमुख योगदान (Key Contributions)

• C-ड्रॉप सिग्नल: आगामी महत्वपूर्ण भूकंपों के लिए एक प्राथमिक, भौतिक रूप से आधारित पूर्ववर्ती के रूप में मैग्नीट्यूड कंप्रेशन में पोस्ट-पीक गिरावट को स्थापित करता है, जो सरल दर-आधारित शांति (rate-based quiescence) से अलग है।
• डेटा-संचालित ज़ोन अपघटन (Decomposition): यह अध्ययन प्रदर्शित करता है कि समग्र भूकंपीय क्षेत्रों को विशिष्ट सक्रिय भ्रंश धागों में विभाजित किया जा सकता है। उदाहरण के लिए, यह अध्ययन अनंतनाग क्षेत्र में सक्रिय लोडिंग सिग्नल को डिस्चार्ज किए गए मुख्य केंद्रीय थ्रस्ट के बजाय पंजल थ्रस्ट (रम्बन-बनिहाल) के लिए अलग करता है।
• प्रोस्पेक्टिव वेरिफिकेशन (Prospective Verification): यह शोध पत्र एक औपचारिक भविष्यवाणी प्रस्तुत करता है जो जून 18, 2026 को जेनेडो (Zenodo) पर आर्काइव की गई है, जो सक्रिय घटनाओं के होने से पहले की है, जिससे यह एक सत्यापन योग्य रिकॉर्ड प्रदान करता है।
• क्रॉस-टेक्टोनिक सत्यापन: इस पद्धति का परीक्षण चार प्लेट बाउंड्री प्रकारों (टकराव, सबडक्शन, स्ट्राइक-स्लिप, और संपीड़न आर्क) में पूर्वव्यापी रूप से किया गया है, जो C-ड्रॉप तंत्र की सार्वभौमिकता को प्रदर्शित करता है।

4. परिणाम (Results)

पूर्वव्यापी सत्यापन (Retrospective Validation):
चार विवर्तनिक क्षेत्रों में 53 पुष्ट C-ड्रॉप चक्रों में से, एक लक्षित-परिमाण वाला भूकंप 90 दिनों के भीतर 83.0% मामलों (44/53) में घटित हुआ।

• जापान (सबडक्शन): 88.9% डिटेक्शन दर।
• तुर्की (ट्रांसफॉर्म): 91.7% डिटेक्शन दर।
• इटली (कंप्रेशन): 75.0% डिटेक्शन दर।
• भारत/कश्मीर (कोलिजन): 73.3% डिटेक्शन दर।
• उच्च-विश्वास उपसमूह: उन चक्रों के लिए जहाँ $C_{peak} \ge 0.30$ और ड्रॉप $\ge 0.055$ था, सभी क्षेत्रों में डिटेक्शन दर 100% तक पहुँच गई।

कश्मीर ज़ोन की स्थिति (20 जून, 2026 तक):

• ज़ोन Z21 (अनंतनाग): स्कोर 5/6। $C_{drop} = 0.129$। 11 अप्रैल और 5 जून, 2026 को दो M 4.8 की घटनाएं सक्रिय C-ड्रॉप चरण के दौरान हुईं। चक्र अभी भी खुला है जिसमें $C=0.260$ (पीक लोडिंग का 67%) बना हुआ है।
• ज़ोन Z15b (कश्मीर उत्तर): स्कोर 5/6। $C_{drop} = 0.092$। 26 वर्षों में उच्चतम लोडिंग ($C_{peak}=0.585$)। M $\ge$ 4.5 के लिए ऐतिहासिक डिटेक्शन दर 96.3% है। जुलाई-सितंबर 2026 के लिए M 5.0–6.0 की भविष्यवाणी की गई है।
• ज़ोन Z15 (मुज़ाफ़राबाद): स्कोर 4/6। $C=0.342$ के साथ असाधारण गहरी शांति ($Q=0.795$)। C-ड्रॉप पुष्टि (0.050 थ्रेशोल्ड) 2-4 सप्ताह के भीतर अपेक्षित है।
• ज़ोन FZ-NPB (गिलगित): स्कोर 4/6। अध्ययन में उच्चतम $C$ मान (0.530)। C-ड्रॉप अभी शुरू हुआ है; अगस्त-अक्टूबर 2026 के लिए एक उभरती हुई भविष्यवाणी को फ्लैग किया गया है।

सत्यापन घटनाएं (Validation Events):
तीन घटनाओं ने भविष्योन्मुखी भविष्यवाणी विंडो को सत्यापित किया:

1. M 3.3 (जून 19, 2026): ज़ोन Z21 में 13 दिनों की चुप्पी को तोड़ते हुए, जेनेडो आर्काइव के 18 घंटे बाद घटित हुआ।
2. M 4.8 (अप्रैल 11, 2026) और M 4.8 (जून 5, 2026): दोनों ज़ोन Z21 के सक्रिय C-ड्रॉप चरण के दौरान हुए, जो आंशिक तनाव रिलीज का प्रतिनिधित्व करते हैं।

5. महत्व और दावे (Significance and Claims)

यह शोध पत्र दावा करता है कि यह कश्मीर हिमालयन फॉल्ट बेल्ट के लिए एक प्रोस्पेक्टिव, भविष्योन्मुखी रूप से सत्यापन योग्य भूकंप भविष्यवाणी प्रदान करता है, जो वर्तमान भूकंप विज्ञान अभ्यास में एक महत्वपूर्ण अंतर को संबोधित करता है।

• भौतिक तंत्र: लेखक तर्क देते हैं कि C-ड्रॉप एक सार्वभौमिक भ्रंश यांत्रिकी प्रक्रिया को दर्शाता है जहाँ भार (load) के तहत परिमाण भिन्नता कम हो जाती है और विफलता (failure) पर बढ़ जाती है, जो विवर्तनिक सेटिंग से स्वतंत्र है।
• परिशुद्धता (Precision): डेटा-संचालित ज़ोन अपघटन का उपयोग करके, यह अध्ययन तनाव के स्थानिक एट्रिब्यूशन को परिष्कृत करता है, जिससे यह पहचान में आता है कि सक्रिय स्रोत पूर्व-माना गया मुख्य केंद्रीय थ्रस्ट के बजाय पंजल थ्रस्ट है।
• जोखिम मूल्यांकन: अध्ययन रेखांकित करता है कि लगभग 1 करोड़ लोग जम्मू और कश्मीर, पाकिस्तान-प्रशासित कश्मीर और गिलगित-बाल्टिस्तान के प्रभावित क्षेत्रों में रहते हैं।
• कार्यवाही का आह्वान: लेखक अंतर्राष्ट्रीय भूकंप विज्ञान समुदाय (NCS और CSEP सहित) को इस भविष्यवाणी की स्वतंत्र रूप से निगरानी और मूल्यांकन करने के लिए आमंत्रित करते है, इस बात पर जोर देते हुए कि यह पद्धति स्वयं-सुसंगत (self-consistent) है जो लोड किए गए (अनुमान लगाने योग्य) और डिस्चार्ज किए गए (गैर-अनुमान लगाने योग्य) भ्रंश अवस्थाओं के बीच अंतर करती है।

शोध पत्र निष्कर्ष निकालता है कि C-ड्रॉप विधि, डेटा-संचालित ज़ोनिंग और छह-मानदंड फिल्टर के साथ मिलकर, उच्च सांख्यिकीय विश्वास के साथ भविष्योन्मुखी पूर्वानुमान जारी करने के लिए एक मजबूत ढांचा प्रदान करती है।