[S08P-04] Estimating effective normal stress during SSE in Bungo channel from slip velocities and shear stress variations
１．はじめに
スロースリップ(SSE)は、地震波を出さずにプレート境界等を10日～数年かけて数十cmゆっくりすべる現象である。SSEはGNSS観測以降発見され、プレート境界でのエネルギー蓄積・解放の新たな1つの過程として位置づけられる。これを調べることはプレート境界のテクトニクス等を知る上で重要である。 速度と状態に依存した摩擦則を用いた理論解析により、速度弱化に伴う準静的なすべり挙動が示されている。Kobayashi and Sato(GRL, 2021, https://doi.org/10.1029/2021GL095690)はGNSSのデータから房総半島沖で発生した６回のSSEの時空間すべり分布を推定し、剪断応力変化とすべり速度の関係をグラフに示した。すると、その軌跡のほとんどは理論的に安定な軌跡と類似していた。これは、房総沖のSSEは準静的な現象であり、すべり加速時の軌跡の傾きはSSEを引き起こしたプレートの摩擦係数（a, b）と有効法線応力（σn）の積（ (a-b)σn ）で表されることが示唆された。本研究では、Kobayashi and Sato(GRL, 2021)で用いた手法を豊後水道で発生したSSEに当てはめ、有効法線応力の絶対量を推定する。
２．解析方法
用いたデータは国土地理院による日々の座標値(F5解)である。SSEによる変動を抽出するために、SSEとは無関係な線形トレンドや季節変動等を推定し、除去する。線形トレンドや季節変動等を除去したデータをフィッティングによってなめらかな曲線で繋ぎ、１ヶ月毎の変動に分割する。次に、１ヶ月毎の変動データからABICを用いたインバージョン解析によって１ヶ月毎のすべり分布を推定する。そしてUSGS(United States Geological Survey)のソフトウェアCoulomb3.3を用いて、プレート境界面を長方形の微小パッチに分割してすべり分布に伴う剪断応力変化をそれぞれのパッチで計算する。求めた剪断応力変化とすべり速度の関係をグラフに表し、加速時の傾きから有効法線応力の絶対量を推定する。
３．結果
GNSSデータを処理し、豊後水道SSEのすべりの時空間分布を推定したところ、すべりは徐々に加速し、そして減速することが示された。現在、剪断応力変化を推定するために、長方形の微小パッチを用いてプレート境界面を分割し、すべり分布に伴う剪断応力変化をそれぞれのパッチで計算している。発表では、この結果を用いて、有効法線応力の絶対量について議論する。
謝辞
解析にあたり、国土地理院の日々の座標値（F5解）を使いました。応力の計算では、米国地質調査所のCoulomb3.3を使いました。ここに記して感謝します。本研究は科研費（23K03541）の補助を受けました。
スロースリップ(SSE)は、地震波を出さずにプレート境界等を10日～数年かけて数十cmゆっくりすべる現象である。SSEはGNSS観測以降発見され、プレート境界でのエネルギー蓄積・解放の新たな1つの過程として位置づけられる。これを調べることはプレート境界のテクトニクス等を知る上で重要である。 速度と状態に依存した摩擦則を用いた理論解析により、速度弱化に伴う準静的なすべり挙動が示されている。Kobayashi and Sato(GRL, 2021, https://doi.org/10.1029/2021GL095690)はGNSSのデータから房総半島沖で発生した６回のSSEの時空間すべり分布を推定し、剪断応力変化とすべり速度の関係をグラフに示した。すると、その軌跡のほとんどは理論的に安定な軌跡と類似していた。これは、房総沖のSSEは準静的な現象であり、すべり加速時の軌跡の傾きはSSEを引き起こしたプレートの摩擦係数（a, b）と有効法線応力（σn）の積（ (a-b)σn ）で表されることが示唆された。本研究では、Kobayashi and Sato(GRL, 2021)で用いた手法を豊後水道で発生したSSEに当てはめ、有効法線応力の絶対量を推定する。
２．解析方法
用いたデータは国土地理院による日々の座標値(F5解)である。SSEによる変動を抽出するために、SSEとは無関係な線形トレンドや季節変動等を推定し、除去する。線形トレンドや季節変動等を除去したデータをフィッティングによってなめらかな曲線で繋ぎ、１ヶ月毎の変動に分割する。次に、１ヶ月毎の変動データからABICを用いたインバージョン解析によって１ヶ月毎のすべり分布を推定する。そしてUSGS(United States Geological Survey)のソフトウェアCoulomb3.3を用いて、プレート境界面を長方形の微小パッチに分割してすべり分布に伴う剪断応力変化をそれぞれのパッチで計算する。求めた剪断応力変化とすべり速度の関係をグラフに表し、加速時の傾きから有効法線応力の絶対量を推定する。
３．結果
GNSSデータを処理し、豊後水道SSEのすべりの時空間分布を推定したところ、すべりは徐々に加速し、そして減速することが示された。現在、剪断応力変化を推定するために、長方形の微小パッチを用いてプレート境界面を分割し、すべり分布に伴う剪断応力変化をそれぞれのパッチで計算している。発表では、この結果を用いて、有効法線応力の絶対量について議論する。
謝辞
解析にあたり、国土地理院の日々の座標値（F5解）を使いました。応力の計算では、米国地質調査所のCoulomb3.3を使いました。ここに記して感謝します。本研究は科研費（23K03541）の補助を受けました。