[STT51-P03] Application of H-matrices method to the calculation of the post-seismic relaxation
When we consider the post-seismic stress change after an earthquake, the effect of viscoelastic deformation of the crust will be important. Recently, a new method based on BIEM is proposed by Barbot and Fialko (2010) in which stress change due to an inelastic strain is calculated as the solution of the inhomogeneous Navier’s equation with equivalent body forces of the inelastic strain. Then, using the stress-strain greenfunction in an elastic medium, we can take into account the inelastic effect.
In this study, we employ their method to evaluate the stress change due to the Nankai/Tonankai earthquakes. Their method requires the computational cost and memory storage of O(N2), where N is the number of discretized cells of the inelastic medium. We reduce the computational cost by applying the fast computation method of H-matrices method. With H-matrices method, a dense matrix is divided into hierarchical structure of submatrices, and each submatrix is approximated to be low rank. When we divide the viscoelastic medium into N = 8,640 or 69,120 uniform cuboid cells and apply the H-matrices method, the required storage memory for the matrices of stress-strain greenfunction are reduced to 0.17 times or 0.05 times of those for the uncompressed original matrices with enough accuracy. In this study, using this method, we show the time development of the stress change at the volcanoes around the Nankai/Tonankai earthquakes, assuming the simple viscos structure. We also discuss the discretized cells and the accuracy for the stress calculation.