Comparison of theoretical and observed temperature profiles in Devon Island ice cap, Canada

Summary. A non-steady-state theoretical model is used to predict the present variation of temperature with depth in two boreholes in the Devon Island ice cap, Arctic Canada. The boreholes are 300 m apart and one of them reaches bedrock. The heat transfer equation is solved numerically with the record of past temperatures obtained from measurements of the variations of oxygen—isotope ratio with depth in the cores as surface boundary condition. The effects of ice advection, refreezing of meltwater percolating from the surface (the amount of which is recorded in the cores), heating due to firn compaction and ice deformation, and heat flow in the bedrock below the ice sheet are all included in the model. The free parameters are geothermal heat flux, present surface temperature and heat loss at the surface which depends on the depth of meltwater penetration and other factors. Agreement between observed and predicted temperature—depth profiles is very close. Latent heat released by percolating meltwater is the predominating factor in determining the temperature distribution in the upper half of each borehole. The temperature distribution is insensitive to the value of the factor used to convert oxygen—isotope ratio to temperature.