Normal mode simulation of prompt elastogravity signals induced by an earthquake rupture

As soon as an earthquake starts, the rupture and the propagation of seismic waves redistribute masses within the Earth. This mass redistribution generates in turn a long-range perturbation of the Earth gravitational field, which can be recorded before the arrival of the direct seismic waves. The recent first observations of such early signals motivate the use of the normal mode theory to model the elastogravity perturbations recorded by a ground-coupled seismometer or gravimeter. Complete modelling by normal mode summation is challenging due to the very large difference in amplitude between the prompt elastogravity signals and the direct P-wave signal. We overcome this problem by introducing a two-step simulation approach. The normal mode approach enables a fast computation of elastogravity signals in layered self-gravitating Earth models. The fast and accurate computation of gravity perturbations indicates instrument locations where signal detection may be achieved, and may prove useful in the implementation of a gravity-based earthquake early warning system.