日本地球惑星科学連合2015年大会

講演情報

ポスター発表

セッション記号 M (領域外・複数領域) » M-TT 計測技術・研究手法

[M-TT40] 地球惑星科学データ解析の新展開:データ駆動型アプローチ

2015年5月27日(水) 18:15 〜 19:30 コンベンションホール (2F)

コンビーナ:*桑谷 立(東北大学大学院環境科学研究科)、駒井 武(東北大学大学院 環境科学研究所)、宮本 英昭(東京大学総合研究博物館)、小池 克明(京都大学大学院工学研究科 都市社会工学専攻地殻環境工学講座)、堀 高峰(独立行政法人海洋研究開発機構・地震津波海域観測研究開発センター)、長尾 大道(東京大学地震研究所)

18:15 〜 19:30

[MTT40-P01] データ同化手法による断層すべり数値モデルの状態・パラメータ推定

*福田 淳一1 (1.東京大学地震研究所)

キーワード:断層すべり, 摩擦パラメータ, 状態推定, パラメータ推定, データ同化, 状態空間モデル

Numerical models of the evolution of fault slip based on rate- and state-dependent frictions laws have been widely used to simulate variety of seismic and aseismic fault slip behavior during the earthquake cycle, including earthquakes, afterslip, slow slip events, and steady aseismic creep. These numerical models consist of the equations of motion and friction law for many subdivided fault patches and include parameters that describe frictional properties of the fault. Previous numerical studies have shown that these parameters, which are often called friction parameters, as well as initial conditions of the model, are one of the governing factors in determining the evolution of fault slip. However, it has been difficult to constrain the friction parameters and initial conditions for natural faults. In this study, we develop methods to simultaneously estimate the friction parameters and initial conditions of the model based on data assimilation methods.

We focus on afterslip, which is transient, decelerating, aseismic slip triggered by stress changes due to a large earthquake. We assume cumulative slip is observable and generate simulated data using the numerical model. A nonlinear state space model is used to relate the numerical model and observations to unknown parameters and dependent variables of the model.

We first adopt the particle filter and smoother to estimate probability distributions of the parameters and dependent variables of the model. In this case, the particle filter/smoother tends to degenerate even if a simplified model and relatively a large number of particles are employed, suggesting that it is difficult to correctly estimate the probability distributions of the unknowns with the particle filter/smoother. We then adopt the ensemble Kalman filter and smoother. We find that the ensemble Kalman filter/smoother does not degenerate even if the number of particles is small, indicating that this method is useful for our problem.

Our results indicate that our methods reproduce the true parameter values reasonably well if the initial conditions are fixed to the true values. In contrast, if the initial conditions are treated as unknowns, the state variable in the rate-state friction law cannot be constrained. In addition, the estimated friction parameters exhibit large uncertainties and a systematic bias. These results suggest that there are significant trade-offs between the friction parameters and initial conditions. We will discuss the details of the results, considerations on how to overcome the trade-offs, results from other data sets, and detailed comparisons of the two methods in the presentation.