Preprints
https://doi.org/10.5194/egusphere-2026-788
https://doi.org/10.5194/egusphere-2026-788
17 Mar 2026
 | 17 Mar 2026
Status: this preprint is open for discussion and under review for The Cryosphere (TC).

Observational data provide valuable insights for glacier thickness reconstruction in High Mountain Asia

Gillian M. A. Smith, Daniel N. Goldberg, Guillaume Jouvet, James R. Maddison, and Hamish D. Pritchard

Abstract. Mountain glaciers provide an irreplaceable water resource in High Mountain Asia, with a significant proportion of water input to rivers coming from glacial meltwater. However, the volume of water held in these glaciers and their evolution over the coming decades is subject to great uncertainty. The reliability of existing glacier ice thickness estimates in High Mountain Asia is limited by the use of low-order models, known to be locally unreliable on mountain glaciers, to describe the relationship between ice velocity and thickness, and by the scarcity of measured thicknesses available for constraint and validation at the time those estimates were produced. We use the Instructed Glacier Model (v2.2.3), a deep-learning-based high-order ice flow model with the capability to invert observed glacier surface velocity for ice thickness, to construct an estimated thickness map of Bhote Kosi glacier. Our thickness inversion is constrained using data collected via a novel airborne radar method for measuring ice thickness. We perform an in-depth case study, carefully justifying inversion parameter choices and quantifying the accuracy of our results. We demonstrate that in the absence of thickness observations, results can be optimized via the use of an L-curve to select the regularization parameter, with significant bias in the unconstrained results, but comparable accuracy to leading thickness estimates. We find that while thickness-constrained inversions are able to correct the modelled thickness field where there is limited information from observed surface velocity, cross-validation experiments demonstrate that the "interpolative power" of thickness observations is weak.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Share
Gillian M. A. Smith, Daniel N. Goldberg, Guillaume Jouvet, James R. Maddison, and Hamish D. Pritchard

Status: open (until 28 Apr 2026)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Gillian M. A. Smith, Daniel N. Goldberg, Guillaume Jouvet, James R. Maddison, and Hamish D. Pritchard

Data sets

Bhote Kosi inversion experiment input files and parameter files Gillian Smith https://doi.org/10.5281/zenodo.18495260

Model code and software

IGM v2.2.3 (extended) IGM developers; Gillian Smith https://doi.org/10.5281/zenodo.18484149

Gillian M. A. Smith, Daniel N. Goldberg, Guillaume Jouvet, James R. Maddison, and Hamish D. Pritchard
Metrics will be available soon.
Latest update: 17 Mar 2026
Download
Short summary
We estimate the thickness of a large glacier in the Himalaya, using recently available thickness measurements and a computational model which uses higher-order physics to match estimated thickness with satellite-observed velocity. Our velocity inversion achieves similar accuracy to leading thickness estimates, while thickness-constrained inversions show increased accuracy, but limited interpolative power. We make recommendations for future measurement locations and for choosing model parameters.
Share