Modelled glacier dynamics over the last quarter of a century at Jakobshavn Isbræ

Abstract. Observations over the past 2 decades show substantial ice loss associated with the speed-up of marine-terminating glaciers in Greenland. Here we use a regional three-dimensional outlet glacier model to simulate the behaviour of Jakobshavn Isbræ (JI) located in western Greenland. Our approach is to model and understand the recent behaviour of JI with a physical process-based model. Using atmospheric forcing and an ocean parametrization we tune our model to reproduce observed frontal changes of JI during 1990–2014. In our simulations, most of the JI retreat during 1990–2014 is driven by the ocean parametrization used and the glacier's subsequent response, which is largely governed by bed geometry. In general, the study shows significant progress in modelling the temporal variability of the flow at JI. Our results suggest that the overall variability in modelled horizontal velocities is a response to variations in terminus position. The model simulates two major accelerations that are consistent with observations of changes in glacier terminus. The first event occurred in 1998 and was triggered by a retreat of the front and moderate thinning of JI prior to 1998. The second event, which started in 2003 and peaked in the summer 2004, was triggered by the final break-up of the floating tongue. This break-up reduced the buttressing at the JI terminus that resulted in further thinning. As the terminus retreated over a reverse bed slope into deeper water, sustained high velocities over the last decade have been observed at JI. Our model provides evidence that the 1998 and 2003 flow accelerations are most likely initiated by the ocean parametrization used but JI's subsequent dynamic response was governed by its own bed geometry. We are unable to reproduce the observed 2010–2012 terminus retreat in our simulations. We attribute this limitation to either inaccuracies in basal topography or to misrepresentations of the climatic forcings that were applied. Nevertheless, the model is able to simulate the previously observed increase in mass loss through 2014.