Higher fidelity seafloor topography and crustal structure models have become available with accumulation of observation data. Previous studies have shown that the consideration of such high-fidelity models produces significant effects, in some cases, on crustal deformation results that are used as inputs for tsunami analysis. However, it is difficult to apply high-fidelity model of crustal deformation computations to tsunami computations because of large computational costs. In this paper, we propose a new crustal deformation computation method for estimating inputs for tsunami computations, which is based on a finite element analysis method with remarkable reduction of computation costs by efficient use of the arithmetic space and the solution space. This finite element analysis method enables us to conduct [Formula: see text]-times crustal deformation computations using high-fidelity models with a degree of freedom on the order of [Formula: see text] for the 2011 Tohoku earthquake example. Tsunami computations with typical settings are conducted as an application example to present the advantages and characteristics of the proposed method. Comparisons between results of the proposed and the conventional method reveal that large shallow fault slip around the trench axis may lead to significant differences in tsunami waveforms and inundation height distributions in some cases.
Fast Finite Element Analysis Method Using Multiple GPUs for Crustal Deformation and its Application to Stochastic Inversion Analysis with Geometry Uncertainty
