Shake and Bake: Structural or Subduction Control on Sumatran Volcanism? A Statistical Approach

技術要約：シェイク・アンド・ベイク（Shake and Bake）：スマトラ火山活動における構造的制御か、それとも沈み込みによる制御か？統計学的アプローチ

問題提起
スマトラにおける火山活動の起源については、沈み込みダイナミクスと右横ずれのスマトラ断層（SF）の相対的な影響力に関して、依然として議論の対象となっている。ホロセーン（完新世）の火山活動の約70%は収束境界で発生しているが、スマトラ弧は、斜め沈み込みがSFを駆動している複雑なテクトニクス体制を呈している。先行研究では相反する見解が示されている。すなわち、火山とSFの近接性は偶然であると主張するものもあれば、特定の火山は構造的に制御されていると示唆するものもある。課題は、主にマグマプロセス（沈み込み駆動型）によって形成された火山と、構造的制御（断層駆動型）の影響を受けた火山を、個々のスマトラ火山の包括的な地球化学データセットの可用性が限られている状況下で、いかに区別するかにある。

手法
著者らは、19のアクティブまたは潜在的に活動的な火山（ホロセーンおよび選定された更新世）を分析するために、物理学に基づいた多変量統計フレームワークを開発した。形態計測の信頼性を確保するため、カルデラおよび高度に浸食された火山の構造体は除外された。

1. データ収集： 本研究では、形態計測パラメータ（MORVOLCアルゴリズムを用いた10m解像度のDEMから派生）と地質学的/構造的パラメータを統合した。主要な変数は以下の通りである：

   • 形態計測： 火山体総体積、平均傾斜、火口の方向、頂部対基部幅比、高さ対基部幅比、および火口数。
   • テクトニクス/地質力学： 海溝からの距離、スラブの深さ、スマトラ断層からの距離、火口の配列、および$\Delta\beta$値。
   • 地震学： $\Delta\beta$値は、火山下のグーテンベルク・リヒターの$\beta$値の偏差を表し、地殻の応力状態や損傷のプロキシ（代理指標）として機能する。

2. 統計解析：

   • 主成分分析 (PCA)： パラメータを標準化し、PCAを用いてデータセットの共分散構造を分解した。
     総分散の分解には、PCAによって得られた主成分スコアを、それぞれの分散で重み付けして合計した分散加重影響行列を作成した。
   • 重み付けスキーム： 断層制御の診断的であるとみなされた特定のパラメータ（断層からの距離、$\Delta\beta$、火口数、頂部/基部幅比、高さ/基部幅比）には1.5の重みを割り当て、その他のパラメータには1.0を割り当てた。
   • 断層親和性指数 (FAI)： 各火山のFAIは、SFに制御されていることが独立して認められている3つの火山（Lubukraya、Sibualbuali、Talang）によって定義された「ベンチマーク・シグネチャ」に対して、各火山の加重影響プロファイルを比較することによって導出された。正のFAI値は、SF制御のシグネチャとの類似性を示し、負の値は非類似性を示す。

3. 地球化学的および地質力学的コンテキスト： マグマの強制力のみが形態の変化を説明できるかどうかを検証するため、著者らは、直交する収束速度に基づくモンテカルロ法（Goltz et al., 2026に基づく）を用いて、マントル融解計算を行い、弧に沿ったマグマフラックスの変動を推定した。また、既存のコンピレーションから沿弧的な微量元素系（La/Sm, Sm/Yb, Th/Yb, Ce/Pb）を分析した。

主な結果

• 統計的クラスタリング： PCAの結果、最初の3つの主成分が全分散の67.7%を説明した。火山はPC1の負の側に一貫したグループとしてクラスター化しており、これらはSFへの距離が短く、中〜高の$\Delta\beta$、整列した火口、大きな体積、および特定の幅の比率を特徴としている。
• FAIによる分類：
  • 高い正のFAI（強い構造的親和性）： ベンチマーク火山（Lubukraya, Sibualbuali, Talang）のSF親和性を確認し、新たに4つの候補を特定した：Dempo, Sorikmarapi, Seulawah Agam, Malintang。これらの構造体は、断層制御と一致する形態計測学的および地震学的シグネチャを示している。
  • 負のFAI（弱い/なしの構造的親和性）： Rajabasa, Peuet Sague, Telong Bur ni, Sinabung, Singgalang, Tandikat, Marapi, Talakmau などの火山は負のFAIを示した。それらの形態（例：急傾斜、大きな断層距離、SFに対して垂直な火口の配列）は、SFの影響とは一致しない。
  • 中間/中立： Kerinci と Sumbing は、ゼロに近い、あるいは緩やかな正のFAIを示し、限定的または混合的な影響を示唆している。
• マグマの強制力： マントル融解モデルは、北スマトラから南スマトラにかけてのマグマフラックスの変動が極めて小さいこと（~55から62 km³/Myr）を示しており、これは比較的均一なスラブの深さと一致している。微量元素系は、抑制された一次的なマグマ強制力を示した。著者らは、マグマの変動性のみでは観察された火山形態の多様性を説明できないと結論付けている。
• 火山の年代： 噴火年代（ホロセーン vs 更新世）とFAIとの間に系統的な関係は見られなかった。これは、構造的シグネチャが、利用可能な解像度内における最新の噴火のタイミングとは独立していることを示唆している。

意義および主張
本論文は、ハイブリッドなテクトニクス環境において、構造的に制御された火山と、主にマグマプロセスによって形成された火山を区別するための、再現可能で定量的な根拠を提供すると主張している。マグマフラックスの変動が形態的な差異を説明するには不十分であることを示すことで、本研究は、スマトラ断層が少なくとも一部のスマトラ火山に対して、重要な構造的影響を及ぼしているという仮説を補強している。

「構造的シグネチャ」の特定は、火山ハザードの再評価を可能にする。著者らは、SFへの構造的親和性が高い火山は、短周期の地震トリガリングに対して特有の感受性を持つ可能性があると示唆している。具体的には、巨大地震がSFに荷重をかけ、それに続く断層の滑りイベントがSF関連の火山におけるマグマ貯留層を乱し、噴火を誘発する可能性がある。この二重のハザード経路（直接的なメガスラスト応力 ＋ 断層制御による応力再分配）は、スマトラにおける火山リスクにおける、潜在的に過小評価されている要素である。

限界
著者らは、本研究が探索的なものであり、比較的規模の小さいデータセットに基づいていることを認めている。火山固有の地球化学データの体系的な欠如により、マグマ源の特性を統計モデルに直接統合することが困難となっている。さらに、FAIはベンチマークに対する類似性の相対的な尺度であり、因果関係の制御に関する決定的な確率ではない。形態計測パラメータは、浸食や構造体の保存状態に関連する不確実性の影響を受ける。本研究は、これらの関係を精緻化するために、地球物理学的イメージング、岩石学的モデリング、およびフィールドベースの構造マッピングを統合した今後の研究を求めている。