Incorporating Horizontal Density Variations into Large‐scale Modelling of Ice Masses

• Published in: Journal of Geophysical Research: Earth Surface
• Main Authors: Schelpe, Camilla, Gudmundsson, Hilmar
• Format: Article in Journal/Newspaper
• Language: English
• Published: Wiley-Blackwell 2023
• Subjects: F700 Ocean Sciences; Antarctic; The Antarctic; West Antarctic Ice Sheet; Antarc*; Ice Sheet; Ice Shelves
• Online Access: https://nrl.northumbria.ac.uk/id/eprint/51277/; https://doi.org/10.1029/2022jf006744; https://nrl.northumbria.ac.uk/id/eprint/51277/15/JGR%20Earth%20Surface%20-%202023%20-%20Schelpe%20-%20Incorporating%20Horizontal%20Density%20Variations%20Into%20Large%E2%80%90Scale%20Modeling%20of%20Ice%20Masses.pdf; https://nrl.northumbria.ac.uk/id/eprint/51277/1/JGR%20Earth%20Surface%20-%202023%20-%20Schelpe%20-%20Incorporating%20Horizontal%20Density%20Variations%20into%20Large%E2%80%90scale%20Modelling%20of%20Ice%20Masses.pdf; https://nrl.northumbria.ac.uk/id/eprint/51277/2/ManuscriptFinal.pdf

Gravity-driven flow of large ice masses such as the Antarctic Ice Sheet (AIS) depends on both the geometry and the mass density of the ice sheet. The vertical density profile can be approximated as pure ice overlain by a firn layer of varying thickness, and for the AIS the firn thickness is not uncommonly 10 to 20% of the total thickness, leading to not insignificant variation in density. Nevertheless, in most vertically-integrated ice-flow models today the density is assumed constant, sometimes with an adjustment in thickness to compensate. In this study, we explore the treatment of horizontal density variations (HDVs) within vertically-integrated ice-sheet models. We assess the relative merits and shortcomings of previously proposed approaches, and provide new formulations for including HDVs. We use perturbation analysis to derive analytical solutions that describe the impact of density variations on ice flow for both grounded ice and floating ice shelves, which reveal significant qualitative differences between each of the proposed density formulations. Furthermore, by modelling the transient evolution of a large sector of the West Antarctic Ice Sheet (WAIS), we quantify the potential impact of HDVs on estimated sea level change. For the domain we considered, we find that explicitly including the horizontal density gradients in the momentum and mass conservation equations leads to about a 10% correction in the estimated change in volume above flotation over 40 years. We conclude that including horizontal density variations in flow modelling of the Antarctic Ice Sheet is important for accurate predictions of mass loss.