Discussion of ‘A Simple Thaw‐Freeze Algorithm for a Multi‐Layered Soil using the Stefan Equation’ by Xie and Gough (2013)

• Published in: Permafrost and Periglacial Processes
• Main Author: Kurylyk, Barret L.
• Other Authors: Natural Sciences and Engineering Research Council of Canada, The Killam Trusts
• Format: Article in Journal/Newspaper
• Language: English
• Published: Wiley 2015
• Subjects: Gough; Active layer thickness; permafrost; Permafrost and Periglacial Processes
• Online Access: http://dx.doi.org/10.1002/ppp.1834; https://api.wiley.com/onlinelibrary/tdm/v1/articles/10.1002%2Fppp.1834; https://onlinelibrary.wiley.com/doi/pdf/10.1002/ppp.1834

Abstract Xie and Gough (2013) developed an algorithm to expand the application of the well‐known Stefan equation in order to predict the freeze or thaw depth in multi‐layered soils. This discussion paper demonstrates that the development of this algorithm is premised on mathematical expressions that are not physically tenable. A previously proposed equation for applying the Stefan approach to multi‐layered soils is herein derived from first principles using the appropriate governing equations, boundary conditions and initial conditions. Numerical methods are employed to demonstrate the accuracy of this equation derived from first principles. This approach can be implemented in existing cold regions hydrology models to estimate the active layer thickness in permafrost regions underlain by layered soils. Copyright © 2015 John Wiley & Sons, Ltd.