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

講演情報

[E] 口頭発表

セッション記号 S (固体地球科学) » S-CG 固体地球科学複合領域・一般

[S-CG40] Science of slow-to-fast earthquakes

2024年5月29日(水) 09:00 〜 10:15 コンベンションホール (CH-B) (幕張メッセ国際会議場)

コンビーナ:加藤 愛太郎(東京大学地震研究所)、山口 飛鳥(東京大学大気海洋研究所)、濱田 洋平(国立研究開発法人海洋研究開発機構)、野田 朱美(気象庁気象研究所)、座長:松澤 孝紀(国立研究開発法人 防災科学技術研究所)、佐藤 大祐(海洋研究開発機構)

09:15 〜 09:30

[SCG40-22] Granular simulation for earthquakes in geological scale structure evolution

*古市 幹人1、Chen Jian1、西浦 泰介1新井 隆太1山本 由弦2井出 哲3 (1.海洋研究開発機構 、2.神戸大学、3.東京大学)

キーワード:地震モデル、スラスト形成、シームレスモデル、DEM、HPC

We present virtual earthquake simulation using a numerical granular rock box experiment based on the Discrete Element Method (DEM). This simulation extends the numerical sandbox test by incorporating cohesive contact forces. We simulate a thin granular rock layer sized 100km in convergence direction, 250m laterally, and about 2.0km in gravity direction, using about 5.6 million DEM elements. The horizontal shortening of this layer generates sequential thrust formation like that of an accretionary wedge [1]. Intermittent fast motions of elements generating seismic waves, are observed through convergence tests following the development of geological-scale structures [2]. These results suggest that the granular rock box simulation seamlessly reproduced seismological and geological-scale phenomena. The hypocenter, which is not predefined by the model, emerges within the active fault damage zone and shifts as the geological structures evolved. The seismic event shows the 3D elementwise rotation and the double-couple behavior, and Gutenberg–Richter law regarding kinetic energy. Quantitative analysis of a seismic event shows that the fracture propagation speed is ~2.6 km/s, and the P-wave and S-wave velocities are ~3.7 and ~2.7 km/s, respectively. The slip distance and the mean stress drop at the hypocenter are ~1.8 cm and ~0.46 MPa, respectively. These characteristics demonstrate the feasibility of granular rock box simulation for reproducing multiscale mechanisms that bridge seismological and geological phenomena. On the other hand, we acknowledge the limitations of the model in terms of fast convergence speed, large element size, and element shape.

[1] Furuichi, M., Chen. J., Nishiura, D., Arai, and R., Yamamoto, Y., Tectonophysics, 862, 5, 229963 (2023). 654 DOI:10.1016/j.tecto.2023.229963.
[2] Furuichi, M., Chen. J., Nishiura, D., Arai, and R., Yamamoto, Y., Ide, S., Tectonophysics, 874, 11, 230230 (2024). DOI:10.1016/j.tecto.2024.230230