Analytical approach for the toroidal relaxation of viscoelastic earth

Abstract

This paper is concerned with post-seismic toroidal deformation in a spherically symmetric, non-rotating, linear-viscoelastic, isotropic Maxwell earth model. Analytical expressions for characteristic relaxation times and relaxation strengths are found for viscoelastic toroidal deformation, associated with surface tangential stress, when there are two to five layers between the core-mantle boundary and Earth's surface. The multilayered models can include lithosphere, asthenosphere, upper and lower mantles and even low-viscosity ductile layer in the lithosphere. The analytical approach is self-consistent in that the Heaviside isostatic solution agrees with fluid limit. The analytical solution can be used for high-precision simulation of the toroidal relaxation in five-layer earths and the results can also be considered as a benchmark for numerical methods. Analytical solution gives only stable decaying modes-unstable mode, conjugate complex mode and modes of relevant poles with orders larger than 1, are all excluded, and the total number of modes is found to be just the number of viscoelastic layers between the core-mantle boundary and Earth's surface - however, any elastic layer between two viscoelastic layers is also counted. This confirms previous finding where numerical method (i.e. propagator matrix method) is used. We have studied the relaxation times of a lot of models and found the propagator matrix method to agree very well with those from analytical results. In addition, the asthenosphere and lithospheric ductile layer are found to have large effects on the amplitude of post-seismic deformation. This also confirms the findings of previous works. © 2006 The Authors Journal compilation © 2006 RAS.