The effects of compressibility on the GIA in southeast Alaska

• The regional glacial isostatic rebound in southeast Alaska is modeled.
• The effects of compressibility exceed observational precisions.
• Incompressible models can underestimate the mantle viscosity by approximately 30%.
• Flexural rigidity should be adjusted when using incompressible models.

Recent theoretical simulations on the glacial isostatic adjustment (GIA) have revealed that the model differences arising from considering mantle compressibility are not necessarily negligible if compared with the observation accuracy of present-day deformation rates. In this study, a compressible model is constructed for the GIA in southeast Alaska, and the uplift rate is compared with GPS data and the incompressible case for the first time. It is shown that, for Maxwell rheology, the incompressible model potentially underestimates the mean uplift rate by approximately 27% (4 mm/yr) with respect to the compressible case and the difference is detectable given observational precision. This difference between the compressible and incompressible models is reduced to 10% by matching the flexural rigidity of both earth models. When carrying out an inversion using incompressible models, this adjustment is important to infer a physically more correct viscoelastic structure.