Open AccessMathematical modeling and numerical simulation of polythermal glaciers
Author(s) -
Aschwanden A.,
Blatter H.
Publication year - 2009
Publication title -
journal of geophysical research: earth surface
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.67
H-Index - 298
eISSN - 2156-2202
pISSN - 0148-0227
DOI - 10.1029/2008jf001028
Subject(s) - latent heat , enthalpy , stagnation enthalpy , mechanics , heat flux , sea ice growth processes , geology , thermodynamics , physics , sea ice , heat transfer , climatology , cryosphere , sea ice thickness
A mathematical model for polythermal glaciers and ice sheets is presented. The enthalpy balance equation is solved in cold and temperate ice together using an enthalpy gradient method. To obtain a relationship between enthalpy, temperature, and water content, we apply a brine pocket parameterization scheme known from sea ice modeling. The proposed enthalpy formulation offers two advantages: (1) the discontinuity at the cold‐temperate transition surface is avoided, and (2) no treatment of the transition as an internal free boundary is required. Fourier's law and Fick‐type diffusion are assumed for sensible heat flux in cold ice and latent heat flux in temperate ice, respectively. The method is tested on Storglaciären, northern Sweden. Numerical simulations are carried out with a commercial finite element code. A sensitivity study reveals a wide range of applicability and defines the limits of the method. Realistic temperature and moisture fields are obtained over a large range of parameters.