A Structural-State Effective-Friction Closure for Oceanic-Plate Subduction Interfaces

技術要約：海洋プレート沈み込み界面における構造状態有効摩擦閉鎖モデル

問題提起
海洋プレートの沈み込み界面は、埋没し、水和し、熱的に進化し、かつ幾何学的に複雑であるため、その摩擦抵抗をフルスケールで直接測定することは不可能である。この界面は、固着したアスペリティ、クリープ領域、流体、堆積物、および変質鉱物の不均質な混合物で構成されている。したがって、見かけの抵抗は単純な実験室レベルの摩擦対ではなく、「有効」なフィールドスケールの現象である。当面の観測上の課題は、アクセス不可能な構造状態を、熱流量、有効摩擦係数、地殻変動結合、および地震モーメント放出といった観測量へと投影する、有効なフィールドスケールの閉鎖式（クロージャ）として機能する有用な公式を開発することである。本研究では、他のアクセス不可能な大規模摩擦問題に適用されてきた「構造状態摩擦アーキテクチャ」が、沈み込みの観測データに対して意味のある投影が可能かどうかを検証する。

手法
本研究では、構造状態閉鎖モデルを現実世界のデータに対して検証するために、二段構えの観測監査を実施した。

1. 熱流量および熱モデル監査：

   • 様々な沈み込み帯における熱流量および熱モデルから導出された有効摩擦係数（$\mu'$）をまとめた、文献コード化された制約条件のデータセットを編纂した。
   • これらの制約条件は、通常、低い有効摩擦値（0.03から0.07の範囲）を示しており、これは乾燥した実験室レベルのバイアリー摩擦（$\mu \approx 0.6$）よりも著しく低い。
   • 縮小された構造状態有効摩擦閉鎖モデルをこれらの区間制約に適合させ、一定の有効摩擦ベースラインと比較した。

2. USGS地震モーメント監査：

   • 1900年から現在までのマグニチュード $\ge$ 7の地震カタログデータを、USGS FDSN APIを介して12の代表的な沈み込み領域からダウンロードし、合計715件のイベントを取得した。
   • 各領域の総地震モーメント（$M_0$）を、長期的なスティック・スリップによるモーメント放出の粗い投影として算出した。
   • 以下のモデルを含む様々な投影モデルを、このモーメントデータに対してテストした：
     • 収束速度のみ。
     • 熱状態のみ。
     • 単純な「構造的負荷指数」。
     • 速度、熱状態、およびスラブの傾斜（ディップ）を組み合わせた補正モデル。
   • パフォーマンスは、リージョン間の転移性を評価するために、leave-one-out (LOO) 交差検証を用いて評価された。

構造状態閉鎖モデル
提案された閉鎖モデルは、有効摩擦（$\mu_{eff}$）をプレートの収束速度（$V$）と、縮小された熱・流体弱化状態（$\Theta$）の関数として定義する。モデルは以下の4段階のアーキテクチャで構成されている：

1. 物理的実行可能集合： 沈み込み力学によって課される（例：$0 \le \mu_{eff} \le \mu_{Byerlee}$）。
2. 構造座標： 固着、粗さ、流体含有量、変質鉱物、クリープ、および幾何学を統合した複合変数（$\Phi$）。
3. 構造応答： 構造座標に対するシステムの応答を表す、生産損失形式。
4. 投影： 状態を測定可能な観測量（熱流量、結合、スロー・スリップ、地震モーメント）へとマッピングする。

主な結果

• 熱流量監査： 構造状態閉鎖モデルは、編纂されたすべての有効摩擦区間（0.03–0.07）内に適合することに成功した。これは、最良の一定有効摩擦ベースライン（0.00004225）と比較して、わずかに小さい区間損失（0.00003544）をもたらした。両モデルとも、乾燥した実験室レベルの摩擦値0.6を正しく拒絶しており、沈み込みの観測量が、流体および温度によって修飾された弱い有効係数を表していることを裏付けている。
• 地震モーメント監査：
  • 速度のみの失敗： 収束速度のみでは、総地震モーメントの予測因子としては不十分である（LOO RMSE = 0.8293）。例えば、トンガ・ケルマデック帯は収束速度が高いが、ペルー・チリや日本のような低速の領域よりも総モーメントが低い。
  • 単純な構造指数の失敗： 単純な構造的負荷指数は性能が悪く（LOO RMSE = 0.7210）、この公式が普遍的な速度のみの摩擦則として扱うことはできないことを示している。
  • 熱・幾何学的成功： 熱状態とスラブの幾何学を組み込んだモデルは、著しく優れた性能を示した。「速度・熱・傾斜」投影は、非カタログモデルの中で最も低いLOO RMSE（0.6005）を達成し、速度のみおよび熱のみのベースラインを上回った。「structural_interaction」モデルはさらに高い誤差（LOO RMSE = 0.8876）を示した。

意義および主張
本論文は、熱力学的沈み込みモデル、速度・状態摩擦理論、または地震サイクル・シミュレーションに取って代わることを目的としていないと明記しており、その範囲を控えめに保っている。代わりに、特定のアーキテクチャ的アプローチのパイロットレベルの検証を提供している。

最も強く支持される結論は以下の通りである：

1. 適合性： 沈み込みの有効摩擦は、構造状態の閉鎖モデルと適合する。このモデルは、熱流量研究で見られる低い摩擦区間を正常にカバーしている。
2. 投影の必要性： この閉鎖の観測可能な投影には、熱・流体状態およびスラブの幾何学を含めることが不可欠である。収束速度や単純な構造指数のみに基づく公式では、リージョンごとの地震モーメント放出を説明するには不十分である。
3. 普遍的な係数ではなくアーキテクチャ： 本研究は、沈み込み摩擦の「普遍的」な側面とは、固定された数値係数ではなく、一つのアーキテクチャであることを示唆している。すなわち、硬い実行可能な沈み込み状態 $\rightarrow$ 構造座標 $\rightarrow$ 生産損失応答 $\rightarrow$ 熱・流体・幾何学的投影 $\rightarrow$ 観測量 というプロセスである。

要約すると、本研究は、熱・流体および幾何学的投影項によって補正される場合に限り、構造状態有効摩擦閉鎖モデルが沈み込み界面に対して実行可能な枠組みであることを示す証拠を提供している。外部への普遍性は未証明であり、現在の証拠は決定的な普遍法則というよりも、パイロットレベルの閉鎖モデルを支持している。