Several Hundred Finite Element Analyses of an Inversion of Earthquake Fault Slip Distribution using a High-fidelity Model of the Crustal Structure

To improve the accuracy of inverse analysis of earthquake fault slip distribution, we performed several hundred analyses using a 108-degree-of-freedom finite element (FE) model of the crustal structure in an efficient way based on capacity computing. we developed a meshing method and an efficient computational method for these large FE models, which methods are specialized for the problem setting. We applied the model to the inverse analysis of coseismic fault slip distribution for the 2011 Tohoku-oki Earthquake. The high resolution of our model provided a significant improvement of the fidelity of the simulation results compared to existing computational approaches.

[1]  Takane Hori,et al.  Afterslip distribution following the 2003 Tokachioki earthquake: An estimation based on the Green’s functions for an inhomogeneous elastic space with subsurface structure , 2010 .

[2]  Tsuyoshi Ichimura,et al.  Fast numerical simulation of crustal deformation using a three-dimensional high-fidelity model , 2013 .

[3]  Yuji Yagi,et al.  Introduction of uncertainty of Green's function into waveform inversion for seismic source processes , 2011 .

[4]  Charles L. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[5]  Motoyuki Kido,et al.  Trench‐normal variation in observed seafloor displacements associated with the 2011 Tohoku‐Oki earthquake , 2011 .

[6]  Akira Asada,et al.  Displacement Above the Hypocenter of the 2011 Tohoku-Oki Earthquake , 2011, Science.

[7]  Walter D. Mooney,et al.  Poroelastic stress-triggering of the 2005 M8.7 Nias earthquake by the 2004 M9.2 Sumatra–Andaman earthquake , 2010 .

[8]  Arthur Raefsky,et al.  A simple and efficient method for introducing faults into finite element computations , 1981 .

[9]  Y. Kaneda,et al.  Effect of elastic inhomogeneity on the surface displacements in the northeastern Japan: Based on three-dimensional numerical modeling , 2007 .

[10]  Motoyuki Kido,et al.  Frontal wedge deformation near the source region of the 2011 Tohoku‐Oki earthquake , 2011 .

[11]  O. C. Zienkiewicz,et al.  A novel boundary infinite element , 1983 .

[12]  Gene H. Golub,et al.  Inexact Preconditioned Conjugate Gradient Method with Inner-Outer Iteration , 1999, SIAM J. Sci. Comput..

[13]  Brendan J. Meade,et al.  Geodetic imaging of coseismic slip and postseismic afterslip: Sparsity promoting methods applied to the great Tohoku earthquake , 2012 .

[14]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[15]  Masaya Matsuura,et al.  Geodetic data inversion using a Bayesian information criterion for spatial distribution of fault slip , 1992 .

[16]  Thomas J. R. Hughes,et al.  Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies , 1985 .