An adjoint-based FEM optimization of coseismic displacements following the 2011 Tohoku earthquake: new insights for the limits of the upper plate rebound

Abstract

The 11 March 2011 Tohoku earthquake was the strongest event recorded in recent historic seismicity in Japan. Several researchers reported the deformation and possible mechanism as triggered by a mega thrust fault located offshore at the interface between the Pacific and the Okhotsk Plate. The studies to estimate the deformation in detail and the dynamics involved are still in progress. In this paper, coseismic GPS displacements associated with Tohoku earthquake are used to infer the amount of slip on the fault plane. Starting from the fault displacements configuration proposed by Caltech-JPL ARIA group and Geoazur CNRS, an optimization of these displacements is performed by developing a 3D finite element method (FEM) model, including the data of GPS-acoustic stations located offshore. The optimization is performed for different scenarios which include the presence of topography and bathymetry (DEM) as well as medium heterogeneities. By mean of the optimized displacement distribution for the most complete case (heterogeneous with DEM), a broad slip distribution, not narrowly centered east of hypocenter, is inferred. The resulting displacement map suggests that the beginning of the area of subsidence is not at east of MYGW GPS-acoustic station, as some researchers have suggested, and that the area of polar reversal of the vertical displacement is rather located at west of MYGW. The new fault slip distribution fits well for all the stations at ground and offshore and provides new information on the earthquake generation process and on the kinematics of Northern Japan area.