Material versus isobaric internal boundaries in the Earth and their influence on postglacial rebound

Most previous earth models used to calculate viscoelastic relaxation after the removal of the Late Pleistocene ice loads implicitly assume that there is no exchange of mass across the mantle density discontinuities on periods of tens of thousands of years (the material boundary formulation). In the present study, simple incompressible models are used to determine the Earth's behaviour in the case where the density discontinuity remains at a constant pressure rather than deforming with the material (the isobaric boundary formulation). The calculation of the movement of the boundary is more rigorous than in earlier studies and uses the local incremental pressure calculated at the depth of the boundary and allows for the vertical deformation caused by the change in volume as material changes phase. It is shown that the buoyancy modes associated with the density discontinuities decrease in strength and increase in relaxation time analogous to what results when the density contrast is reduced. Also, two viscoelastic modes arise from an isobaric boundary, which is also predicted when there is a contrast in rigidity or viscosity across a material boundary. The difference in predicted radial deformation between the isobaric boundary model and the material boundary model is largest for long-wavelength loads for which the material incremental pressure at depth is largest. If the isobaric boundary model is appropriate for the treatment of the mineral phase changes in the mantle on glacial rebound timescales, then previous inferences of the deep-mantle to shallow-mantle viscosity ratio based on large-scale deformation (spherical harmonic degree < 10) of the Earth and including data from the early part of the glacio-isostatic uplift are too small.