The choice of the physical model of postglacial rebound plays a decisiverole to derive information about the mantle rheology and viscosity from observed data. In models for the mantle rheology, an incompressible Maxwell material is often assumed in spite of seismic observations showing that the Earth's mantle is composed of compressible material.
In this study, in order to assess the influence of compressibility on glacial isostatic adjustment (GIA), the spectral-finite element approach proposed by Martinec is extended to incorporate the effect of compressibility.