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

講演情報

[J] ポスター発表

セッション記号 S (固体地球科学) » S-SS 地震学

[S-SS13] 強震動・地震災害

2019年5月27日(月) 17:15 〜 18:30 ポスター会場 (幕張メッセ国際展示場 8ホール)

コンビーナ:栗山 雅之(一般財団法人 電力中央研究所 地球工学研究所 地震工学領域)、染井 一寛(一般財団法人地域地盤環境研究所)

[SSS13-P07] 強震波形記録を用いて推定された平成30年北海道胆振東部地震の震源過程

*久保 久彦1鈴木 亘1岩城 麻子1青井 真1関口 春子2 (1.国立研究開発法人防災科学技術研究所、2.京都大学防災研究所)

キーワード:平成30年北海道胆振東部地震、震源過程、強震動

The 2018 Hokkaido Eastern Iburi earthquake (Mjma 6.7) struck southern Hokkaido, Japan on September 6, 2018 (JST). This earthquake caused strong ground motions over southern Hokkaido with a maximum seismic intensity of 7 on the JMA scale and a maximum PGA of over 1500 cm/s/s. In this study, we estimate the source process of this earthquake using strong motion waveforms.
For source-process analysis of this event, we develop a curved fault model based on the spatial distribution of aftershocks. This curved fault model consists of multiple planes, each with a width of 20 km, a common dip angle (65 degree), and a common top depth (approximately 22 km). The common dip angle refers to the moment tensor analysis of F-net. The strike angle of the planes smoothly changes along the strike direction: the strike angle is 15° at the north part, 340° at the central part, and 40° at the south part. The rupture starting point is set to the hypocenter determined by JMA. The curved fault model is divided into 15 subfaults along the strike direction and 10 subfaults along the dip direction.
We use velocity waveforms at 20 stations of K-NET and KiK-net, F-net of NIED. The waveforms are band-pass filtered between 0.05 and 0.5 Hz, resampled to 10 Hz, and windowed from 1 s before S-wave arrival for 25 s.
Green's functions are calculated with the discrete wavenumber method (Bouchon 1981) and the reflection/transmission matrix method (Kennett and Kerry 1979) assuming 1-D velocity structure models. The structure models are obtained for each station from the 3-D structure model (Fujiwara et al. 2009). Logging information is also used for the KiK-net station. To consider the rupture propagation effect, 25 point-sources are uniformly distributed over each subfault in the calculation of Green’s functions.
The multi-time-window linear waveform inversion method (Olson and Apsel 1982; Hartzell and Heaton 1983) is used in this study. The slip time history of each subfault is discretized using eight smoothed ramp functions (time windows) progressively delayed by 0.4 s and having a duration of 0.8 s each. The triggering velocity of first time window is set to 1.4 km/s to minimize the data-fit residual. Two orthogonal slips of each time window at each subfault are derived by minimizing the difference between the observed and synthetic waveforms using the non-negative least-squares scheme (Lawson and Hanson 1974). The slip angle is allowed to vary centered at 107°, which is the rake angle of the F-net analysis. In addition, we impose a spatiotemporal smoothing constraint on slips (Sekiguchi et al. 2000).
In the estimated source model, the large slips are found at approximately 30 km depth in the up-dip direction from the hypocenter. In the first 6 s, the rupture grows slowly around the hypocenter with small slips. Then, the rupture developed with the large moment release between 6 s and 12 s in the large slip area. The large slip area does not overlap with the active aftershock area