1Department of Civil and Environmental Engineering, University of Trento, Trento, Italy
2Environment Canada – National Hydrology Research Centre, Saskatoon, Saskatchewan, Canada
3Department of Geography, University of Zurich, Winterthurerstrasse 190 Zurich, Switzerland
*now at: Mountain-eering srl, Via Siemens 19 Bolzano, Italy
Abstract. In this paper we provide a method for solving the energy equation in freezing soil. The solver is linked with the solution of Richards equation, and therefore able to approximate water movement near the liquid-solid phase transition. The equations show non-linear characteristics causing oscillatory behavior in the solution close to the phase transition, when normal methods of iterative integration, as Newton or Picard, are used. Thus, a globally convergent Newton method has been implemented to achieve convergence. The method is tested by comparison with an analytical solution to the Stefan problem and by comparison with experimental data derived from literature.