Soil mass and heat transfer processes were represented in three levels of model complexities to understand soil freeze–thaw mechanisms. Results indicate that coupled mass and heat transfer models considerably improved simulations of the soil hydrothermal regime. Vapor flow and thermal effects on water flow are the main mechanisms for the improvements. Given the explicit consideration of airflow, vapor flow and its effects on heat transfer were enhanced during the freeze–thaw transition period.