Material versus local incompressibility and its influence on glacial-isostatic adjustment

We present an analytical form of the layer propagator matrix for the response of a locally incompressible, layered, linear-viscoelastic sphere to an external load assuming that the initial density stratification rho (o)(r) within each layer is parametrized by Darwin's law. From this, we show that the relaxation of a sphere consisting of locally incompressible layers is governed by a discrete set of viscous modes. The explicit dependence of the layer propagator matrix on the Laplace transform variable allows us to determine the amplitudes of the viscous modes analytically. Employing Darwin's parametrization, we construct three simplified earth models with different initial density gradients that are used to compare the effects of the local incompressibility constraint, div (rho (o)u)=0, and the material incompressibility constraint, div u=0, on viscoelastic relaxation. We show that a locally incompressible earth model relaxes faster than a materially incompressible model. This is a consequence of the fact that the perturbations of the initial density are zero during viscoelastic relaxation of a locally incompressible medium, so that there are no internal buoyancy forces associated with the continuous radial density gradients, only the buoyancy forces generated by internal density discontinuities. On the other hand, slowly decaying internal buoyancy forces in a materially incompressible earth model cause it to reach the hydrostatic equilibrium after a considerably longer time than a locally incompressible model. It is important to note that the approximation of local incompressibility provides a solution for a compressible earth model that is superior to the conventional solutions for a compressible earth with homogeneous layers because it is based on an initial state that is consistent with the assumption of compressibility.