Coseismic rotation changes from the 2004 Sumatra earthquake: the effects of Earth's compressibility versus earthquake induced topography

Polar motion is modelled for the large 2004 Sumatra earthquake via dislocation theory for an incompressible elastic earth model, where inertia perturbations are due to earthquake-triggered topography of density–contrast interfaces, and for a compressible model, where inertia perturbation due to compression-dilatation of Earth's material is included; density and elastic parameters are based on a multilayered reference Earth. Both models are based on analytical Green's functions, propagated from the centre to the Earth's surface. Preliminary and updated seismological solutions are considered in elucidating the effects of improving earthquake parameters on polar motion. The large Sumatra thrust earthquake was particularly efficient in driving polar motion since it was responsible for large material displacements occurring orthogonally to the strike of the earthquake and to the Earth's surface, as imaged by GRACE gravity anomalies over the earthquake area. The effects of earthquake-induced topography are four times larger than the effects of Earth's compressibility, for l = 2 geopotential components. For varying compressional Earth properties and seismic solution, modelled polar motion ranges from 8.6 to 9.4 cm in amplitude and between 117° and 130° east longitude in direction. The close relationship between polar motion direction, earthquake longitude and thrust nature of the event, are established in terms of basic physical concepts.