A Dual-Permeability Approach for Modeling Soil Water Flow and Heat Transport during Freezing and Thawing

Abstract

Preferential flow may become significant in partially frozen soils because infiltration can occur through large, initially air-filled pores surrounded by a soil matrix with limited infiltration capacity. The objectives of this study were to develop and evaluate a dual-permeability approach for simulating water flow and heat transport in macroporous soils undergoing freezing and thawing. This was achieved by introducing physically based equations for soil freezing and thawing into the dual-permeability model MACRO. Richards’ equation and the heat flow equation were loosely coupled using the generalized Clapeyron equation for the soil micropore domain. Freezing and thawing of macropore water is governed by a first-order equation for energy transfer between the micropore and macropore domains. We assumed that macropore water was unaffected by capillary forces, so that water in macropores freezes at 0°C. The performance of the model was evaluated for four test cases: (i) redistribution of water in the micropore domain during freezing, (ii) a comparison between the first-order energy transfer approach and the heat conduction equation, (iii) infiltration and water flow in frozen soil with an initially air-filled macropore domain, and (iv) thawing from the soil surface during constant-rate rainfall. Results show that the model behaves in accordance with the current understanding of water flow and heat transport in frozen macroporous soil. To improve modeling of water and heat flow in frozen soils, attention should now be focused on providing experimental data suitable for evaluating models that account for macropore flow.