Deformation Due to a Pair of Infinite and a Pair of infinite-finite planar faults in Fractional Standard Linear Solid

技術要約：分数階標準線形固体モデルにおける一対の無限および一対の無限・有限平面断層による変形

問題提起
地震活動が活発な地域には、ある断層に沿った滑りが隣接する断層の応力蓄積や分布に大きな影響を与える、相互作用する複数の断層が頻繁に含まれている。従来の研究では、粘弾性媒体（例：マックスウェルモデルや標準線形固体モデル）における整数階微分を用いて断層間の相互作用をモデル化してきたが、これらの伝統的な手法は、地質学的材料に固有の「記憶」効果を見落としがちである。本論文では、以下の2つの特定のシナリオにおける、クリープする横ずれ断層の相互作用に伴う地殻の変形を扱うことで、この問題に対処する。

1. ケースI： 二つの無限埋没断層間の相互作用。
2. ケースII： 一つの無限埋没断層と、一つの有限な地表破断断層との相互作用。

本研究は、地球の媒体が**分数階標準線形固体（FSLS）**モデルに基づく粘弾性半空間として振る舞うと仮定しており、これは整数階モデルと比較して、地殻岩石の時間依存的なレオロジー挙動と記憶効果をより適切に捉えるものである。

手法
著者らは、変位、応力、および歪み成分の解析解を導出するために、理論的な数学モデルを開発した。その手法は、以下の主要なステップで構成される。

• 数学的定式化： 媒体は、カプト（Caputo）分数階微分演算子（${}_{0}D^{\alpha}_{t}$、ここで $\alpha \in (0, 1]$）を含む分数階標準線形固体構成方程式を用いて定義される。応力－歪み関係には、有効剛性（$\mu$）および有効粘性（$\eta$）が組み込まれる。
• 支配方程式： 解析は、慣性力を無視した準静的な変形を仮定している。応力の運動方程式および境界条件は、無限および有限の断層幾何学の両方に対して確立されている。
• 解法： 時間依存の分数階微分方程式を扱うために、ラプラス変換を採用している。解は以下を用いて導出される：
  • 対応原理（Rybicki, 1971）。
  • 断層面における転位条件を考慮するための修正グリーン関数法（Maruyama, 1966）。
  • 分数階演算子の逆ラプラス変換から自然に現れるミッタッハ・レフラー関数（$E_{\alpha}$）。
• 分析されたシナリオ：
  • 断層クリープの「前」（テクトニクス応力の蓄積）の状態と、「後」（滑りによる転位）の状態に関する解を導出している。
  • クリープ速度（$V_1, V_2$）および深さ依存のクリープ関数を適用し、断層 $F_1$ および $F_2$ の動きをシミュレートしている。
  • 数値計算は、地質学的データ（例：剛性 $\mu = 3.5 \times 10^{10}$ N/m²、粘性 $\eta = 5 \times 10^{19}$ Pa·s）から得られたパラメータを用いて行われ、様々な断層傾斜角（$\theta$）、分数階（$\alpha$）、およびクリープ速度に対してグラフ結果を生成している。

主な貢献

1. 断層相互作用の分数階モデリング： 本論文は、分数階微積分を組み込むことで既存の断層相互作用モデルを拡張し、リソスフェアの記憶依存的な粘弾性挙動をより正確に表現することを可能にした。
2. 複雑な幾何学形状に対する解析解： 本研究は、一対の無限埋没断層および、無限と有限の地表破断断層の混合システムという二つの異なる構成について、変位、応力、および歪み場の閉じた形式の解析解を提供している。
3. 記憶効果の定量化： 本研究は、分数階の次数 $\alpha$ が応力の蓄積と解放にどのように影響するかを明示的に示し、数学的パラメータと地質学的媒体の物理的な「記憶」を結びつけている。

結果および数値的知見
グラフ表現に裏付けられた数値解析により、以下の具体的な観察結果が得られている。

• 変位の大きさ： 断層付近の表面変位の変化は、年間 $10^{-3}$ m のオーダーである。これは、観測された地球物理学的データ（例：USGS報告による 2～6 mm/年）とよく一致しており、従来の整数階モデルが予測していた $10^{-5}$ の範囲とは乖離している。
• 断層幾何学の影響：
  • ケースI（二つの無限断層）では、第二の断層付近の変位量は第一の断層付近よりも大きい。
  • ケースII（無限および有限の断層）では、有限の断層付近の変位量は無限の断層付近よりも小さい。
  • 変位は断層トレース付近でピークに達し、距離が増すにつれて減少し、最終的にゼロに収束する。
• 分数階の次数（$\alpha$）の影響： 研究の結果、応力の蓄積および解放の速度は $\alpha$ に強く依存することが判明した。$\alpha$ が低い値（強い記憶効果を示す）の場合、せん断応力の蓄積量が高くなり、緩和が遅くなる。$\alpha$ が 1 に近づく（整数挙動に近づく）につれて、応力解放の速度が増加する。
• 駆動力と応力分布： 断層間の総駆動力は、クリープの開始後に急速に解放される。この解放の大きさは、断層の傾斜角およびクリープ速度によって変化する。等値線図は、応力の蓄積が自由表面および断層転位帯の近くで最も顕著であり、応力の減少は断層の位置に対応する深さで顕著であることを示している。
• 歪み解析： せん断歪み（$e_{12}$）の変動は $10^{-4}$ の範囲で見られ、これは 2011 年の東北沖地震のような実際の事象のデータと一致している。

意義と主張
著者らは、自らの分数階モデルが、従来の整数階モデルと比較して、より一般的かつ正確な断層相互作用ダイナミクスの表現を提供すると主張している。具体的には：

• 分数階のアプローチは、標準的なモデルでは見落とされがちな、粘弾性媒体の記憶効果と非線形性を捉えることに成功している。
• モデルによる変位速度（$10^{-3}$）の予測は、観測された地球物理学的データと強い一致を示しており、分数階微積分が地震挙動や地下変形を予測するための優れたツールであることを示唆している。
• 知見は、断層間の相互作用が単一の断層よりも広範囲にストレス場に影響を与えること、および分数階の次数 $\alpha$ が応力解放の速度を決定する重要なパラメータであることを示している。

本論文は、これらの理論的モデルが、地下変形およびそれが断層ダイナミクスや地震発生に与える影響の理解を深め、地震活動領域におけるアセイスミック（無地震）期間を分析するための堅牢な枠組みを提供すると結論付けている。