[S08P-18] Dynamic rupture simulation of the 2011 Mj7.0 Fukushima Hamadori, Japan, earthquake
2011年Mj7.0福島県浜通りの地震 (浜通り地震) は, 単一の断層だけでなく, 複数の断層が立て続けに破壊する複雑な過程を経て発生したと考えられている. 現在までに, InSAR解析で得られた地表面変位 (例: Kobayashi et al., 2012; Fukushima et al., 2013) や強震波形 (例: Hikima, 2012) を用いた逆解析によって浜通り地震発生時の複数断層上での破壊すべりの時空間分布を調べる試みが行われてきた. しかし一方で, 複数断層の連続的な破壊の発生過程を動的破壊シミュレーションによって調べる試みは, 投稿者の知る限りでは現在まで行われていない. そこで本研究では, 浜通り地震に関する動的破壊シミュレーションを実施し, それを通して浜通り地震発生時に生じた複数断層の連続的な破壊の発生過程を明らかにすることを目指す.
動的破壊シミュレーションの実施時に設定する必要のある要素は主に「断層形状」「広域応力場」「断層面上に課す摩擦則」の3つである. 断層形状に関しては, 浜通り地震発生後の余震分布および浜通り地震発生時に地表面上に現れた塩ノ平断層, 井戸沢断層, 湯之岳断層の3断層それぞれに対応する変位不連続 (Kobayashi et al., 2012) を利用して制約する. 余震分布を利用した断層形状の制約の際に用いる震源カタログとしては, Kato et al. (2013) にてdouble-difference tomography method (Zhang and Thurber, 2003) を用いて構築された浜通り地震の震源域における震源カタログとその更新版 (解析期間: 2015年12月31日まで延長) を選択する. 広域応力場に関しては, Uchide et al. (2022) にて構築された日本全国内陸部の地殻内応力マップを参照して, 浜通り地震の震源域における広域応力場の設定を行う. また, 断層面上に課す摩擦則は線形すべり弱化摩擦則 (Ida, 1972) を仮定する. 以上のように各要素を設定した上で, FDP-BIEM (Ando, 2016) を用いて動的破壊シミュレーションを実施する.
本発表では, 上記のようにして実施する動的破壊シミュレーションの予備的結果を報告し, 浜通り地震発生時に生じたとされる複数断層の連続的な破壊の発生過程について議論する予定である.
動的破壊シミュレーションの実施時に設定する必要のある要素は主に「断層形状」「広域応力場」「断層面上に課す摩擦則」の3つである. 断層形状に関しては, 浜通り地震発生後の余震分布および浜通り地震発生時に地表面上に現れた塩ノ平断層, 井戸沢断層, 湯之岳断層の3断層それぞれに対応する変位不連続 (Kobayashi et al., 2012) を利用して制約する. 余震分布を利用した断層形状の制約の際に用いる震源カタログとしては, Kato et al. (2013) にてdouble-difference tomography method (Zhang and Thurber, 2003) を用いて構築された浜通り地震の震源域における震源カタログとその更新版 (解析期間: 2015年12月31日まで延長) を選択する. 広域応力場に関しては, Uchide et al. (2022) にて構築された日本全国内陸部の地殻内応力マップを参照して, 浜通り地震の震源域における広域応力場の設定を行う. また, 断層面上に課す摩擦則は線形すべり弱化摩擦則 (Ida, 1972) を仮定する. 以上のように各要素を設定した上で, FDP-BIEM (Ando, 2016) を用いて動的破壊シミュレーションを実施する.
本発表では, 上記のようにして実施する動的破壊シミュレーションの予備的結果を報告し, 浜通り地震発生時に生じたとされる複数断層の連続的な破壊の発生過程について議論する予定である.