Crustal heterogeneity effects on coseismic deformation: numerical simulation of the 2008 MW 7.9 Wenchuan earthquake

Abstract

Coseismic deformation of large earthquakes causes significant property damages and fatalities, which requires quantitative research of multiple disciplines such as geodesy, geological investigation, seismic tomography, and seismic dislocation theory. The finite element method accounts for material heterogeneity and geometric complexity, making it suitable for studying the coseismic deformation of large earthquakes. This paper develops a parallel elastic finite element program that utilizes split nodes and high-performance parallel computing technology on the FELAC software platform to study the coseismic deformation of large earthquakes. We verify the accuracy of the parallel elastic finite element program by comparing its results with the analytical solutions from seismic dislocation theory for four ideal earthquake cases. Finally, we established parallel elastic finite element models to study the coseismic deformation of the 2008 Wenchuan earthquake. The simulation results are consistent with the GPS and InSAR data. Coseismic surface deformation results are significantly influenced by medium regional heterogeneity with different layered structures besides the Longmenshan fault. The finite element program lay the foundation for the inversion of the coseismic fault rupture process based on the heterogeneous medium model and complex geometric model.