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

1. 問題の所在

カシミール・ヒマラヤ弧は、インドプレートとユーラシアプレートの北向きの収束によって駆動される、世界で最も地震災害のリスクが高い地域の一つである。この収束は、並行する断層系（主前縁断層、主境界断層、パンジャル断層、および主中央断層）に沿って、膨大な弾性歪エネルギーを蓄積させている。2005年のマグニチュード（M）7.6のカシミール地震のような歴史的な大災害にもかかわらず、この地域における現在の地震予測手法には、以下の3つの具体的な限界が存在する。

• 異なるテクトニクス設定において有効な、単一の定量的閾値の欠如。
• カシミール地域に特化した、将来予測として検証可能な記録の不在。
• 広域的な地理的領域ではなく、特定の断層系に負荷信号を帰属させる能力の欠如。

地震静穏期（Z統計量）、b値マッピング、加速モーメント放出などの既存の手法は、これらの目標を同時に達成できていない。

2. 手法

本研究は、**マグニチュード圧縮指数（$C$）とCドロップ（C-drop）**信号を中心としたデータ駆動型のフレームワークを導入し、4つの異なるテクトニクス体制（カシミール、日本、トルコ、イタリア）を用いて検証を行っている。

中核となる指標:

• マグニチュード圧縮指数（$C$）: パラメータフリーの尺度であり、$C(k) = \max(0, 1 - \sigma[W] / \sigma_0)$ と定義される。これは、断層が応力を蓄積するにつれて地震のマグニチュード分布がどのように狭まるかを定量化するものである。応力負荷が局所的な地震活動を均質化するにつれ、マグニチュードの標準偏差（$\sigma$）が減少し、これにより$C$は1に向かって上昇する。
• Cドロップ（C-drop）: 主要な前駆現象信号である。これは、ピークからの系統的かつ持続的な$C$の減少（ピークから少なくとも0.050の減少）として定義される。この減少は、断層が破壊を開始し、より大きな地震が観測ウィンドウに入り込むことで、マグニチュードの分散が増大することによって起こる。
• 地震静穏指数（$Q$）: 二次的な指標であり、発生率の減少を示す（$Q = \max(0, 1 - \lambda[W] / \lambda_0)$）。高まった$Q$と高まった$C$の組み合わせは、臨界状態にあるロックされた断層を示唆する。

データ処理およびゾーン特定:

• ローリングウィンドウ: 計算には、ゾーン内の直近$W=20$個の地震を用いたローリングウィンドウを使用し、12イベントの移動平均によって平滑化を行う。
• ベースライン: 基準値（$\sigma_0, \lambda_0$）は、2000年から2009年までの期間から導出される。
• データ駆動型ゾーニング: 任意の座標を用いる従来の研究とは異なり、ゾーンの中心はISCカタログ（2000–2026）のガウス平滑化によって特定された、自然な地震密度ピークに配置される。ゾーンの半径は、局所的な地震クラスターの90%を捉える。
• 6つの基準によるスコアリングシステム: 偽陽性（例：余震列）を排除するため、以下の6つの基準のうち少なくとも3つを同時に満たした場合のみ、そのゾーンを**クリティカル（CRITICAL）**と分類する。
  1. 圧縮: $C \ge 0.25$。
  2. 静穏: $Q \ge 0.20$。
  3. Cドロップ: ピークからの確定的な減少。
  4. b値: $b < 0.75$（応力集中を示す）。
  5. 低活動: 直近90日間で15イベント未満。
  6. 新鮮さ: ピークが過去180日以内に発生。

3. 主な貢献

• Cドロップ信号: 単純な発生率ベースの静穏期とは異なる、差し迫った重大な地震の前駆現象として、ポストピークのマグニチュード圧縮の減少を、物理的根拠に基づいた主要な信号として確立した。
• データ駆動型ゾーン分解: 複合的な地震ゾーンを分解することで、特定の活動的な断層系を特定できることを示した。例えば、本研究ではアナントナグ地域の活動的な負荷信号を、放電した主中央断層コアではなく、パンジャル断層（ランバン–バニハル）に特定している。
• 将来予測の検証可能性: 本研究は、事象の発生前に（2026年6月18日に）Zenodoにアーカイブされた正式な予測を提示しており、手法の適用に関する検証可能な記録を提供している。
• クロス・テクトニクス検証: 本手法を4つの異なるプレート境界タイプ（衝突、沈み込み、横ずれ、および圧縮弧）に対して遡及的にテストし、Cドロップメカニズムの普遍性を実証した。

4. 結果

遡及的検証:
4つのテクトニクス体制における53件の確定したCドロップサイクルにおいて、ターゲットとなるマグニチュードの地震が90日以内に発生した割合は83.0%（44/53）であった。

• 日本（沈み込み帯）: 検知率 88.9%。
• トルコ（トランスフォーム断層）: 検知率 91.7%。
• イタリア（圧縮帯）: 検知率 75.0%。
• インド／カシミール（衝突帯）: 検知率 73.3%。
• 高信頼度サブセット: $C_{peak} \ge 0.30$ かつ ドロップ $\ge 0.055$ のサイクルでは、全地域を通じて検知率は**100%**に達した。

カシミール・ゾーンのステータス（2026年6月20日時点）:

• ゾーン Z21（アナントナグ）: スコア 5/6。$C_{drop} = 0.129$。2026年4月11日と6月5日にM 4.8のイベントが2回発生した。これらは能動的なCドロップ局面中に発生した。サイクルは継続中であり、$C=0.260$（ピーク負荷の67%を保持）である。
• ゾーン Z15b（北カシミール）: スコア 5/6。$C_{drop} = 0.092$。過去26年間で最高の負荷（$C_{peak}=0.585$）。M $\ge$ 4.5の歴史的検知率は96.3%である。2026年7月〜9月の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年8月〜10月の予測がフラグ立てされている。

検証された事象:
3つの事象が将来予測ウィンドウを検証した：

1. M 3.3（2026年6月19日）: Zenodoアーカイブの18時間後に発生し、ゾーンZ21における13日間の沈黙を破った。
2. M 4.8（2026年4月11日）および M 4.8（2026年6月5日）: いずれもゾーンZ21の能動的なCドロップ局面中に発生しており、部分的な応力解放を表している。

5. 意義および主張

本論文は、現在の地震学の実践における決定的なギャップを埋める、将来予測可能で、将来的に検証可能な地震予測をカシミール・ヒマラヤ断層帯に対して提供すると主張している。

• 物理的メカニズム: 著者らは、Cドロップが、テクトニクス設定に依存せず、負荷の下ではマグニチュードの分散が減少し、破壊に伴って増加するという、普遍的な断層力学プロセスを反映していると論じている。
• 精度: データ駆動型のゾーン分解を用いることで、空間的な帰属を精緻化し、アナントナグの信号の活動的なソースを、以前想定されていた主中央断層ではなくパンジャル断層へと特定している。
• リスク評価: 本研究は、ジャム・カシミール、パキスタン管理下のカシミール、およびギルギット・バルティスタンにわたる影響圏に、約1,000万人の人々が居住していることを強調している。
• 行動喚起: 著者らは、国際的な地震学コミュニティ（NCSおよびCSEPを含む）に対し、この予測を独立して監視し、評価することを求めている。これは、本手法が、負荷された（予測可能な）状態と放電した（予測不可能な）状態を区別する上で自己整合的であることを強調しているものである。

結論として、Cドロップ法は、データ駆動型のゾーニングと6つの基準によるフィルターと組み合わせることで、決定的な断層の負荷状態を特定し、高い統計的信頼度を持って将来の予測を発行するための堅牢な枠組みを提供する。