High-fidelity elastic Green’s functions for subduction zone models consistent with the global standard geodetic reference system

<jats:title>Abstract</jats:title><jats:p>Green’s functions (GFs) for elastic deformation due to unit slip on the fault plane comprise an essential tool for estimating earthquake rupture and underground preparation processes. These estimation results are often applied to generate important information for public such as seismic and tsunami hazard assessments. So, it is important to minimize the distortion of the estimation results on the numerical models used for calculating GFs to guarantee assessment reliability. For this purpose, we here calculated GFs based on a numerical model that is of high fidelity to obtain realistic topography and subsurface structural models of the Earth. We targeted two well-known subduction zones in Japan, the Nankai Trough and the Japan Trench. For these subduction zones, databases for realistic topography and subsurface structural models of the Earth are available in the “Japan integrated velocity structure model version 1”, which was proposed for earthquake hazard assessments conducted by the Japanese government. Furthermore, we eliminated the inconsistency in processing calculated GFs and space geodetic observation data for surface displacements, which is often overlooked, by using the same coordinate system. The ellipsoidal shape of the Earth, which is often approximated with a projected plane or a spherical shape, was also incorporated by faithfully following the definitions of the coordinate systems in Geodetic Reference System 1980, which is the global standard for space geodesy. To calculate elastic GFs based on such high-fidelity subduction zone databases with the ellipsoidal shape of the Earth, we introduced the finite element (FE) method. In the FE meshes, the resolution of the topography and subsurface structure is the same as that of the original databases. Recent development of the state-of-the-art computation techniques for the rapid calculation of crustal deformation using large-scale FE models allows for GF calculation based on such a high-fidelity model. However, it is generally not easy to perform such calculations. Thus, we released a library for the GFs calculated with 1-km grid spacing on the ground surface in this study to the geoscience community on a web server, aiming to contribute more reliable seismic and tsunami hazard assessment.</jats:p>