Numerical Modeling of Coupled Water Flow and Heat Transport in Soil and Snow

Abstract:

A numerical model is developed to calculate coupled water flow and heat transport in seasonally frozen soil and snow. Both liquid water flow and water vapor flow are included. The effect of dissolved ions on soil water freezing point depression is included by combining an expression for osmotic head with the Clapeyron equation and the van Genuchten soil water retention function. The coupled water flow and heat transport equations are solved using the Thomas algorithm and Picard iteration. Ice pressure is always assumed zero and frost heave is neglected. The new model is tested using data from a high-elevation rangeland soil that is subject to significant soil freezing and a mountainous forest soil that is snow-covered for about 8 months of the year. Soil hydraulic parameters are mostly based on measurements and only vegetation parameters are fine-tuned to match measured and calculated soil water content, soil & snow temperature, and snow height. Modeling statistics for both systems show good performance for temperature, intermediate performance for snow height, and relatively low performance for soil water content, in accordance with earlier results with an older version of the model.