Coupled poroelastic and gravitational deformation of a layered spherical Earth – I: load Love numbers

In this paper, we propose a fully coupled two-phase poroelastic deformation theory for a spherically layered and self-gravitating Earth. The earth model consists of a solid inner core, a fluid outer core and poroelastic mantle and crust. The boundary-value problems are posed in the Laplace-transformed domain using the spherical system of vector functions, and analytical solutions are obtained in each layer using the dual-variable and position matrix method. As an application, the surface loading problem is considered. The undrained and drained limits are discussed with the corresponding governing equations, which are different from the conventional ones where body force is ignored. Numerical examples show that the poroelastic effect can be significant since the difference between the results for undrained and drained limits are large with obvious time-dependent deformation. This newly derived theory for the coupled boundary-value problem will have broad applications, such as displacement and gravity change due to groundwater depletion, poroelastic rebound after an earthquake and risk evaluation of earthquakes induced by water injection.