測地データと物理モデルを融合した断層摩擦特性の推定手法の構築

点击翻译按钮获取中文摘要

測地学的に観測される余効変動は主に断層における余効すべりとマントルの粘弾性応力緩和の寄与の和であると考えられている。従って、余効変動の観測データに基づき、断層の摩擦特性やマントルのレオロジーを明らかにすることができると考えられる。本研究では、応力によって駆動される余効すべりと粘弾性緩和の物理モデルを用いて、断層の摩擦構成則のパラメータやマントルのレオロジーのパラメータの空間分布を測地データに基づいて推定する手法の構築を行っている。パラメータの最適値のみならず、その不確実性を明らかにするために、本研究ではパラメータ推定の逆問題をベイズ的に定式化し、パラメータの事後確率分布を推定することを目指している。しかし、パラメータの空間分布を推定する問題では、未知パラメータ数が非常に多くなるため、事後確率分布推定は計算コストが非常に高く困難である。事後確率分布を現実的な計算コストで推定するために、iterative ensemble smootherに基づく近似推定手法を開発した。この手法では、事前確率分布からサンプルを観測データとモデル計算結果に基づいて反復的に更新し、最終的に事後確率分布を近似するガウス分布からのサンプルを得る。この手法の性能を評価するために、モデルを用いて人工的なGNSS時系列データを作成し、このデータからモデルのパラメータ（断層の摩擦構成則のパラメータ、マントルのレオロジーのパラメータ、地震時の応力変化等の空間分布）の事後確率分布を推定した。その結果、推定された事後確率分布の平均が真値を良く再現できることを確認した。また、地震時の応力変化が大きい場所でパラメータの不確実性（事後確率分布の標準偏差）が小さいという妥当な結果が得られた。