Estimating the ice thickness of the Müller Ice Cap using an inversion of the shallow ice approximation

The ice thickness of the Müller Ice Cap, Arctic Canada, is estimated using regression parameters obtained from an inversion of the shallow ice approximation by the use of a single Operation IceBridge flight line in combination with the glacier outline, surface slope, and elevation. The model is compared with an iterative inverse method of estimating the bedrock topography using PISM as a forward model. In both models the surface elevation is given by the Arctic Digital Elevation Model. The root mean squared errors of the ice thickness on the ice cap is 131 m and 139 m for the shallow ice inversion and the PISM model, respectively. Including the outlet glaciers increases the root mean squared errors to 136 m and 396 m, respectively.

The simplicity of the shallow ice inversion model, combined with the good results and the fact that only remote sensing data is needed, means that there is a possibility of applying this model in a global glacier thickness estimate by using the Randolph Glacier Inventory. Most global glacier estimates only provide the volume and not the ice thickness of the glaciers. Hence, global ice thickness models is of great importance in quantifying the potential contribution of sea level rise from the glaciers and ice caps around the globe.