Numerical study of the error sources in the experimental estimation of thermal diffusivity: an application to debris-covered glaciers

In tectonically active mountain regions, the thinning of alpine glaciers due to climate change favors the development of debris covered glaciers. This debris layer significantly modifies a glacier’s melt depending on the debris thickness and therefore modifies its evolution. Debris thermal conductivity is a critical parameter for calculating ice melt beneath a debris layer. The most commonly used method to calculate apparent thermal conductivity of supraglacial debris layers is based on an estimate of volumetric heat capacity of the debris and simple heat diffusion principles presented by Conway and Rasmussen (2000). The analysis of heat diffusion requires a vertical array of temperature measurements through the supraglacial debris cover. This study explores the effect of the temporal and spatial sampling interval, and method on the thermal diffusivity values derived using this method. Results show that increasing temporal and spatial sampling intervals increase truncation errors and therefore systematically underestimate values of thermal diffusivity. Also, the thermistor precision, the shape of the diurnal temperature cycle, and vertical thermistor displacement result in systematic errors. Overall these systematic errors would result in an underestimation of glacier ice melt under a debris layer. We have developed a best practice guideline to help other researchers to investigate the effect of the sampling interval on their calculated sub-debris ice melt and better plan future measurement campaigns.