Realistic Post-seismic deformation in a spherically symmetric Maxwell Earth

The realistic post-seismic deformation in a spherically symmetric earth is obtained for the first time. The calculation is done by the new algorithm which makes it possible to consider the self-gravitation, the compressibility, and a continuously varying structure of the Earth without any approximations. This makes a remarkable contrast to the previous methods which have neglected such effects for the intrinsic numerical difficulties. The essential point of the new algorithm is to evaluate numerically the Laplace integration without taking the sum of the innumerable poles. Using this method, a complete set of the Green's function is shown; time variations of displacement, gravity, geoid height at the surface for a strike-slip, dip-slip, horizontal and vertical tensile point dislocation. As an earth model, we employ the 1066A model and the standard viscosity profiles. The result shows a diverse spatial pattern due to a viscous structure or a source depth. The considerations on the result suggests that what determines those patterns is mainly the thickness of the lithosphere and a relative location of the source. Of particularly interested is that the near filed deformation is mostly prevented by the lithosphere. This makes the post-seismic deformation be distinctive in an epicentral distance of a few hundreds km. This indicates that a post-seismic gravity change is possibly detected by the satellite mission because the wavelength exceeds 100 km, which strongly encourages us to observe it to infer the viscosity under the lithosphere.