Robust reconstruction of glacier beds using transient 2D assimilation with Stokes

Initialising model glaciers such that they match well with their real counterparts and are thus able to make more accurate predictions is an ongoing challenge in glacier modelling. We set out a data-assimilation approach using an ensemble Kalman filter in a 2D flowline example that provides one possible solution to this problem. We show that our approach is valid across a range of parameters and scenarios, including deliberately data-deficient or inaccurate ones, and leads to robust retrieval of the glacier bed. We also provide some suggestions for how best to use data assimilation within a mountain-glacier context.