Modelling sea-level fingerprints of glaciated regions with low mantle viscosity

• Published in: Earth System Dynamics
• Main Authors: A. Bartholet, G. A. Milne, K. Latychev
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2021
• Subjects: Science; Q; Geology; QE1-996.5; Dynamic and structural geology; QE500-639.5; Antarctic; Canada; Antarc*; glacier; glacier*; glaciers; Alaska
• Online Access: https://doi.org/10.5194/esd-12-783-2021; https://doaj.org/article/ca9d514b1aa642f7a249e6376e30085c

Global patterns of sea-level change – often termed “sea-level fingerprints” – associated with future changes in ice/water mass re-distribution are a key component in generating regional sea-level projections. Calculation of these fingerprints is commonly based on the assumption that the isostatic response of the Earth is dominantly elastic on century timescales. While this assumption is accurate for regions underlain by mantle material with viscosity close to that of global average estimates, recent work focusing on the West Antarctic region has shown that this assumption can lead to significant error where the viscosity is significantly lower than typical global average values. Here, we test this assumption for fingerprints associated with glaciers and ice caps. We compare output from a (1D) elastic Earth model to that of a 3D viscoelastic model that includes low-viscosity mantle in three glaciated regions: Alaska, southwestern Canada, and the southern Andes (Randolph Glacier Inventory (RGI) regions 1, 2, and 17, respectively). This comparison indicates that the error incurred by ignoring the non-elastic response is of the order of 1 mm in most areas (or about 1 % of the barystatic signal) over the 21st century with values reaching the centimetre level in glaciated regions. However, in glaciated regions underlain by low-viscosity mantle, the non-elastic deformation can result in relative sea-level changes with magnitudes of up to several tens of centimetres (or several times the barystatic value). The magnitude and spatial pattern of this non-elastic signal is sensitive to variations in both the projected ice history and regional viscosity structure, indicating the need for loading models with high spatial resolution and improved constraints on regional Earth viscosity structure to accurately simulate sea-level fingerprints in these regions. The anomalously low mantle viscosity in these regions also amplifies the glacial isostatic adjustment signal associated with glacier changes during the 20th century, ...