Simulation of a fully coupled 3D glacial isostatic adjustment – ice sheet model for the Antarctic ice sheet over a glacial cycle

• Published in: Geoscientific Model Development 2023-09, Vol.16 (18), p.5473-5492
• Main Authors: van Calcar, Caroline J, van de Wal, Roderik S. W, Blank, Bas, de Boer, Bas, van der Wal, Wouter
• Format: Article
• Language: English
• Subjects: Analysis; Antarctic ice sheet; Bedrock; Coupling; Deformation; Deglaciation; Earth; Earth mantle; Environmental Sciences; Evolution; Feedback; Finite element method; Glaciation; Glaciology; Ice; Ice cover; Ice loads; Ice melting; Ice sheet models; Ice sheets; Ice shelves; Ice thickness; Lithosphere; Mantle; Mathematical models; Meltwater; Modelling; Polar environments; Response time; Rheology; Sea level; Sheet modelling; Simulation; Stabilizing; Temporal resolution; Thickness; Three dimensional models; Topography; Viscoelasticity; Viscosity
• ISSN: 1991-9603; 1991-959X; 1991-962X; 1991-9603; 1991-962X; 1991-959X
• DOI: 10.5194/gmd-16-5473-2023

Glacial isostatic adjustment (GIA) has a stabilizing effect on the evolution of the Antarctic ice sheet by reducing the grounding line migration following ice melt. The timescale and strength of this feedback depends on the spatially varying viscosity of the Earth's mantle. Most studies assume a relatively long and laterally homogenous response time of the bedrock. However, the mantle viscosity is spatially variable, with a high mantle viscosity beneath East Antarctica and a low mantle viscosity beneath West Antarctica. For this study, we have developed a new method to couple a 3D GIA model and an ice sheet model to study the interaction between the solid Earth and the Antarctic ice sheet during the last glacial cycle. With this method, the ice sheet model and GIA model exchange ice thickness and bedrock elevation during a fully coupled transient experiment. The feedback effect is taken into account with a high temporal resolution, where the coupling time steps between the ice sheet and GIA model are 5000 years over the glaciation phase and vary between 500 and 1000 years over the deglaciation phase of the last glacial cycle. During each coupling time step, the bedrock elevation is adjusted at every ice sheet model time step, and the deformation is computed for a linearly changing ice load. We applied the method using the ice sheet model ANICE and a 3D GIA finite element model. We used results from a regional seismic model for Antarctica embedded in the global seismic model SMEAN2 to determine the patterns in the mantle viscosity. The results of simulations over the last glacial cycle show that differences in mantle viscosity of an order of magnitude can lead to differences in the grounding line position up to 700 km and to differences in ice thickness of the order of 2 km for the present day near the Ross Embayment. These results underline and quantify the importance of including local GIA feedback effects in ice sheet models when simulating the Antarctic ice sheet evolution over the last glacial cycle.