日本地震学会2021年度秋季大会

講演情報

A会場

一般セッション » S15. 強震動・地震災害

PM-1

2021年10月14日(木) 13:30 〜 15:00 A会場 (A会場)

座長:野津 厚(港湾空港技術研究所)、吉田 昌平(大崎総合研究所)

14:00 〜 14:15

[S15-08] Estimate of the rupture process of the 2021 off-Fukushima earthquake based on the finite-fault source inversion Inferred from strong motion records

〇吉田 昌平1、津田 健一2、Chen Ji3、Ralph Archuleta3,4、佐藤 俊明1 (1.大崎総合研究所、2.清水建設、3.カリフォルニア大学サンタバーバラ校、4.RJA Ground Motion Analysis, Inc. )

A Mw 7.1 intra-slab earthquake occurred off the coast of the Fukushima prefecture, northeastern Japan, at 23:07 on 13th February 2021 (JST). This earthquake generated strong shaking (the seismic intensity 6+ and 6- in JMA scale) along the coast of the Fukushima and Miyagi prefecture. In this study, we have estimated the rupture process of this event by the waveform finite-fault inversion method developed by Ji et al. (2002).We use three-components strong-motion records at 16 K-NET and KiK-net surface stations in the Fukushima and Miyagi prefecture. To improve the quality of the inversion results, we have introduced an individual 1D velocity structure model for each station in calculating Green’s functions. First, we have constructed the 1D velocity model for each station. We interpolate the Japan Integrated Velocity Structure Model [JIVSM] (Koketsu et al.,2012) for its subsurface, the seismic tomography model of Matsubara et al.(2019) for crust, and the PREM (e.g., Shearer, 2018) for mantle. Then, we calibrate these models by forward simulation for the record of a Mw 5 aftershock occurring in March 17, 2021. Consequently, we got good agreements in both arrival time and amplitude for the frequency range between 0.05~0.5Hz. For the inversion of the mainshock, we assumed one fault plane dipping 33 degrees to the east based on the aftershock distribution and F-net focal mechanism (Fukuyama et al., 1998). The fault plane was divided into 2km by 2km sub-faults. We optimized five fault parameters (e.g. slip amplitude, rupture time, rake angle, starting and ending time of an assumed analytical slip-rate functions) for each individual sub-fault by minimizing an objective function defined in wavelet domain using simulated annealing (SA) inversion procedure. Our inverted results shown in the figure suggest that the fault rupture mainly propagated towards to the south-west. The area with large slip was detected on this side. The centroid time of the rupture of this asperity was about 10 sec from rupture initiation. Also, the area with large slip velocity is coincide with the area with large slip.