Exponential Distribution for Assessing Earthquake Probabilities (Mw >6.7) within a 100 km Radius in Northeastern Crete

技術要約：北東クレタ島における地震確率評価のための指数分布

問題提起
本研究は、ギリシャの北東クレタ島から半径100 km圏内における主要地震（モーメントマグニチュード $M_w > 6.7$）の局所的な地震ハザードおよび再来間隔を評価するという課題に取り組んでいる。主な目的は、これらの事象の時系列分布が「記憶のない（memoryless）」確率過程に従うかどうかを判断することであり、これにより、歴史的データに基づいた将来の地震活動の条件付き確率を算出することを可能にする。

データ選択およびカタログの構築
本研究では、西暦1500年から1935年までの間に発生した、$N = 10$ 個の主要または壊滅的な地震事象からなる均質なデータセットを利用している。歴史的アーカイブに見られる不完全性を軽減し、データの信頼性を確保するため、著者らは時間窓を限定し、厳格なマグニチュード閾値（$M_w > 6.7$）を適用した。

• 情報源: 歴史的事象の詳細は、欧州歴史地震アーカイブ（AHEAD）から抽出され、計測データはPapazachos & Papazachos (2000) によるカタログと照合された。
• 空間領域: 北東クレタ島を中心とした半径約100 kmの円形領域。
• 完全性: カタログは、記録が明確に残っている1935年2月25日の地震をもって終了している。

手法
著者らは、本震の時系列分布を、主震は記憶のない挙動を示すという標準的な地震学の仮定（余震は適切にフィルタリングされているものとする）の下で、確率的更新過程を用いてモデル化した。

• 分布モデル: 事象間隔（$\Delta t$）について、確率密度関数 $f(\Delta t) = \lambda e^{-\lambda \Delta t}$ で定義される指数分布（ここで $\lambda$ は平均年間発生率）への適合性を検証した。
• パラメータ推定: 発生率パラメータ $\lambda$ は、最大尤度推定（MLE）法を用いて推定された。
• 適合度検定: 実測データと理論的な指数モデルとの間の適合性を検証するため、2つの非パラメトリック検定を実施した：
  1. コルモゴロフ–スミルノフ（K-S）検定。
  2. アンダーソン–ダリング（A-D）検定（小標本に対するブラウン近似を利用）。

主な結果

• 再来率: 推定された平均年間発生率（$\lambda$）は 0.02108 であり、標準誤差は 0.00701 である。$\lambda$ の95%信頼区間は 0.007338 から 0.03482 の範囲にある。
• 統計的妥当性: 適合度検定の結果、高いp値が得られ、データが指数分布に従うという帰無仮説は棄却できないことが示された：
  • K-S検定: $p = 0.959$
  • A-D検定: $p = 0.802$
    これらの結果は、当該地域の再来間隔が指数分布に密接に従っていることを裏付けており、この地域における主震発生の記憶のないポアソン過程という仮定を支持している。
• 確率推定: 累積分布関数（CDF） $P(T \le t) = 1 - e^{-\lambda t}$ を用いて、後続の主要地震の条件付き確率を算出した。最後の事象が1935年に発生したことを踏まえると、本モデルは、90年以上の時間軸（すなわち2025年/2026年まで）において、主要地震（$M_w > 6.7$）が発生する確率は 85% ($0.85$) であると推定している。

意義および主張
本論文は、「簡潔かつ堅牢な統計的枠組み」を地域的な地震ハザード評価に提供すると主張している。北東クレタ島における主要地震の事象間隔が指数分布に従うことを検証することで、本研究は、記憶のないポアソン過程がこの特定の空間的・時間的クラスターに対して適切なモデルであることを示している。主な貢献は、1935年の地震からの経過時間と算出された平均再来率から直接導き出された、将来の主要事象に関する高い条件付き確率（85%）の定量的推定である。著者らは、これらの知見を、新たな物理的メカニズムや将来の実験プロトコルを提案するものではなく、歴史的な地震データに対する統計的更新過程の直接的な応用として提示している。