Topology-Phase Monitoring of Kilauea Eruption Preparation: A Reproducible Seismic Audit and Forecasting Potential Assessment

技術要約：キラウエア噴火準備におけるトポロジー・フェーズ・モニタリング

問題提起
高精度な火山噴火予測は、依然として大きな課題である。なぜなら、噴火の準備過程は、単一のスカラー変数（例：地震回数の急増）によって示されることは稀だからである。むしろ、それは地震活動、地殻変動、ガス輸送、およびマグマと破砕系の連結性の、分散的な再編成を伴うものである。従来のモニタリングは、しばしば孤立したチャネルに依存しており、それらは曖昧になりやすい。例えば、地震群（スウォーム）が必ずしも噴火につながるわけではなく、イベント数の増加は、切迫した活動ではなく応力調整を反映している場合もある。本論文は、噴火の準備とは、地震場の「トポロジー状態」における遷移、すなわち、拡散した背景活動から、連結し、移動し、かつ異方性を持つ状態への移行として理解されるべきであると提唱している。

手法
本研究は、USGS FDSN地震カタログデータを用い、キラウエア火山系に対する再現可能な、地震のみによる監査を実施した。

• データ範囲: キラウエア（19.4069°N, 155.2834°W）の半径80 km圏内の、マグニチュード $M \ge 1.0$ の地震を対象とした。
• 時間窓: データは、2014年から2019年までの同一暦期間（1月1日から5月5日まで）のウィンドウに集計された。2014–2017年のウィンドウは訓練参照分布として用いられ、2018年のウィンドウは、下部東リフトゾーン（LERZ）噴火（開始日：2018年5月3日）のターゲット・テストケースとして用いられた。
• トポロジー状態ベクトル: 7日間のローリングウィンドウを、以下の7つの特徴量からなるベクトルに変換した。
  1. イベント数 ($N$)
  2. マグニチュード・エネルギー・プロキシ ($E$)
  3. 最大連結成分の割合 ($C_{max}$)
  4. グラフのエッジ密度 ($\rho_G$)
  5. 正規化された空間エントロピー ($S_{sp}$)
  6. イベント・クラウドの異方性 ($A$)
  7. セントロイドの移動 ($\Delta r_c$)
  • グラフ構築: イベント間の投影空間上の距離が $\le 8$ km かつ、時間間隔が $\le 2$ 日以内である場合に、それらを連結したものとする。
• 指標の構築: 複合トポロジー・フェーズ指標（TPS）は、標準化された構成要素の平均として算出された。空間エントロピーは、拡散した状態から集中した組織化への移行を反映するため、負の値とした。閾値は、2014–2017年の訓練分布の99パーセンタイルに設定された。

主な結果
2014–2017年のデータで学習されたロックされたパラメータを用いた監査により、2018年のLERZ噴火に関して以下の知見が得られた。

• 早期警戒: 複合トポロジー・フェーズ指標は、2018年4月12日に訓練データの99パーセンタイル閾値を最初に超え、5月3日の噴火開始に対して21日間のリードタイムを提供した。
• スカラー指標との比較:
  • イベント数: 2018年5月1日（噴火開始の2日前）にのみ、99パーセンタイルを超えた。
  • エネルギー・プロキシ: 噴火開始前に閾値を超えなかった。
• 構成要素分析: 早期警戒は、協調的な構造的再編成によってもたらされた。グラフ密度は4月12日に閾値を突破し、最大連結成分の割合と空間エントロピーは4月17日にそれぞれの閾値を突破した。移動（Migration）は3月に早期の突破を示したが、訓練最大値を超えるのは5月1日であった。異方性は1月に早期に突破したが、単独の予測因子としては不十分であると判断された。
• ピーク値: トポロジー指標は2018年5月1日に4.0411のピーク値に達し、訓練の99パーセンタイル（2.0462）および訓練最大値（2.8166）の両方を上回った。

主な貢献

1. トポロジー状態の検出: 本研究は、十分に監視された玄武岩質システムにおいて、噴火の準備が、単純なイベント数やエネルギーの増加とは異なる、地震場のトポロジー的な遷移として検出可能であることを示している。
2. 再現可能な監査: 著者らは、特定のスクリプトとデータファイルを含む、完全に再現可能なワークフローを提供しており、これにより「ロックされたパラメータ」による監査の独立した検証が可能となっている。
3. リードタイムの改善: 連結性、密度、および空間的組織化の指標を組み合わせることで、従来のスカラー・ベースラインと比較して、噴火準備状態を検出するためのリードタイムを大幅に延長できることが示唆されている。

意義と主張
本論文は、控えめかつ具体的な主張を行っている。すなわち、トポロジー状態の物理学が、すべての火山噴火の正確な日付、噴火口の位置、または規模を決定論的に予測できると断言するものではない。

• 確率的な予兆: 主要な貢献は、一つの「噴火準備状態」の確率的な、マルチチャネルによる検出を検証したことにある。この手法は、システムが歴史的なベースラインと比較して、異常に組織化されたトポロジー状態に入ったことを特定するものである。
• 運用の目標: 近接した目標は、決定論的な予測ではなく、トポロジー状態のモニタリングによるリードタイムの向上と、スカラーによる誤報の削減である。
• 限界と今後のステップ: 著者らは、現在の監査が地震のみに基づいていること、および単一の強力な事例（キラウエア2018）に依存していることを認めている。手法が高精度な運用レベルへと進展するためには、以下が必要であるとしている。
  • 誤報率を推定するための、独立した非噴火の不安定期における検証。
  • マルチチャネル・データ（地殻変動、ガス、熱、傾斜）の統合。
  • 標準的なモデル（b値、ETASなど）との比較。
  • 新しい事象に対するテスト前の厳格なパラメータ・ロッキング。

結論として、本論文は、キラウエアの2018年の噴火準備が、地震学的なトポロジーの再編成として検出可能な形で現れたという仮説を支持しており、これは、確率的な状態検出のための潜在的な経路を提示している。