Virtual earthquakes in a numerical granular rock box experiment

Mikito Furuichi,Jian Chen,Daisuke Nishiura,Ryuta Arai,Yuzuru Yamamoto,Satoshi Ide
DOI: https://doi.org/10.1016/j.tecto.2024.230230
IF: 3.66
2024-02-03
Tectonophysics
Abstract:This study presents virtual earthquake simulation using a numerical granular rock box experiment based on the Discrete Element Method (DEM). The horizontal shortening of a granular rock layer with a thin 3D geometry and periodic boundary conditions generated sequential thrust formation like that of an accretionary wedge. Intermittent fast motions of elements along with faults generating seismic waves were observed through a convergence test following the development of geological-scale structures. This indicates that the granular rock box simulation seamlessly reproduced seismological and geological-scale phenomena. The hypocenter, which was not predefined by the model, emerged within the active fault damage zone and shifted as the geological structures evolved. The seismic event occurred with the shear motion of the fault, characterized by the 3D elementwise rotation and the double-couple behavior, and followed the Gutenberg–Richter law. The quantitative analysis of a virtual earthquake in the main frontal thrust estimated that the fracture propagation speed was ~2.6 km/s, and the P-wave and S-wave velocities were ~ 3.7 and ~ 2.7 km/s, respectively. In addition, the slip distance and the mean stress drop at the hypocenter were ~ 1.8 cm and ~ 0.46 MPa, respectively, and the major wave frequency was 3–4 Hz. These results show the feasibility of granular rock box simulation for reproducing realistic multiscale mechanisms bridging seismological and geological phenomena. The limitations of the model in terms of fast convergence speed ( 5×10−4 m/s), large element size ( ~12.5 m), and element shape (sphere) were discussed to achieve a more realistic model.
geochemistry & geophysics
What problem does this paper attempt to address?