Simulating the evolution of the Antarctic Ice Sheet including 3D GIA feedback during the Last Glacial Cycle

Changes in ice load over time deform the Earth’s crust and mantle. This effect, Glacial Isostatic Adjustment (GIA), induces vertical deformation of the bedrock of the Antarctic continent and impacts the grounding line position which is critical for the dynamical state of the Antarctic Ice Sheet (AIS). GIA introduces a negative feedback and stabilizes the ice sheet evolution, hence GIA modelling is important for transient studies. Most ice dynamic models use a two-layer flat Earth approach with a laterally homogenous relaxation time or a layered Earth approach with a laterally homogenous viscosity (1D) to compute the bedrock deformation. However, viscosity of the Earth’s interior varies laterally (3D) and radially with several orders of magnitude across the Antarctic continent. Here we present a new coupled 3D GIA – ice dynamic model which can run over hundred thousands of years with a resolution of 500 years. The method is applied using various 1D and 3D rheologies. Results show that the present-day ice volume is 3 % lower when using a 1D viscosity of 10²¹ Pa·s than using a 3D viscosity. However, local differences in grounding line position maybe up to a hundred kilometres around the Ronne and the Ross Ice Shelfs, and ice thickness differences are up to a kilometre for present day conditions when comparing 1D rheologies and 3D rheologies. The difference between the use of various 3D rheologies is significantly smaller. These results underline and quantify the importance of including local GIA feedback effects in ice dynamic models when simulating the Antarctic Ice Sheet evolution over the Last Glacial Cycle.