2025/04/03 by David Hötten, Hötten, David, Jenny Niebsch +5
Earth and Planetary Sciences · Environmental Science · #65M22 #65M60 #Climate change and permafrost #FOS: Mathematics #Numerical Analysis (math.NA) #Plant Water Relations and Carbon Dynamics #Tree-ring climate responses
paper · pdf · doi:10.48550/arxiv.2504.02475
openalex publication_date 2025/04/03 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We consider the problem of heat conduction with phase change, that is essential for permafrost modeling in Land Surface Models and Dynamic Global Vegetation Models. These models require minimal computational effort and an extremely robust solver for large-scale, long-term simulations. The weak enthalpy formulation of the Stefan problem is used as the mathematical model and a finite element method is employed for the discretization. Leveraging the piecewise affine structure of the nonlinear time-stepping equation system, we demonstrate that this system has a unique solution and provide a solver that is guaranteed to find this solution in a finite number of steps from any initial guess. Comparisons with the Neumann analytical solution and tests in the Lund-Potsdam-Jena managed Land vegetation model reveal that the new method does not introduce significantly higher computational costs than the widely used DECP method while providing greater accuracy. In particular, it avoids a known nonphysical artifact in the solution.