A rapidly converging spin-up method for the present-day Greenland ice sheet using the GRISLI ice-sheet model.

Providing reliable projections of the ice-sheet contribution to future sea-level rise has become one of the main challenges of the ice-sheet modelling community. To increase confidence in future projections, a good knowledge of the present-day state of the ice flow dynamics, which is critically dependent on basal conditions, is strongly needed. The main difficulty is tied to the scarcity of observations at the ice-bed interface at the scale of the whole ice sheet, resulting in poorly constrained parameterisations in ice-sheet models. To circumvent this drawback, inverse modelling approaches can be developed and validated against available data to infer reliable initial conditions of the ice sheet. Here, we present a spin-up method for the Greenland ice sheet using the thermo-mechanical hybrid GRISLI ice-sheet model. Our approach is based on the adjustment of the basal drag coefficient that relates the sliding velocities at the ice-bed interface to basal shear stress in unfrozen bed areas. This method relies on an iterative process in which the basal drag is periodically adjusted in such as way that the simulated ice thickness matches the observed one. The process depends on three parameters controlling the duration and the number of iterations. The best spin-up parameters are chosen according to two criteria to minimize errors in sea-level projections: the final difference between the simulated and the observed Greenland ice volume as well as the final ice volume trend which must both be as low as possible. To increase confidence in the inferred parameters, we also make sure that the final ice thickness root mean square error from the observations is not greater than a few tens of meters. Our best results are obtained after only 420 years of simulation, highlighting a rapid convergence and demonstrating that our method can be used for computationally expensive ice sheet models. [ABSTRACT FROM AUTHOR]