Modeling soil heating and moisture transport under extreme conditions: Forest fires and slash pile burns

Abstract

Heating any soil during a sufficiently intense wildfire or prescribed burn can alter it irreversibly, causing many significant, long‐term biological, chemical, and hydrological effects. Given the climate‐change‐driven increasing probability of wildfires and the increasing use of prescribed burns by land managers, it is important to better understand the dynamics of the coupled heat and moisture transport in soil during these extreme heating events. Furthermore, improved understanding and modeling of heat and mass transport during extreme conditions should provide insights into the associated transport mechanisms under more normal conditions. The present study describes a numerical model developed to simulate soil heat and moisture transport during fires where the surface heating often ranges between 10,000 and 100,000 W m−2 for several minutes to several hours. Basically, the model extends methods commonly used to model coupled heat flow and moisture evaporation at ambient conditions into regions of extreme dryness and heat. But it also incorporates some infrequently used formulations for temperature dependencies of the soil specific heat, thermal conductivity, and the water retention curve, as well as advective effects due to the large changes in volume that occur when liquid water is rapidly volatilized. Model performance is tested against laboratory measurements of soil temperature and moisture changes at several depths during controlled heating events. Qualitatively, the model agrees with the laboratory observations, namely, it simulates an increase in soil moisture ahead of the drying front (due to the condensation of evaporated soil water at the front) and a hiatus in the soil temperature rise during the strongly evaporative stage of the soil drying. Nevertheless, it is shown that the model is incapable of producing a physically realistic solution because it does not (and, in fact, cannot) represent the relationship between soil water potential and soil moisture at extremely low soil moisture contents (i.e., residual or bound water: θ < 0.01 m3 m−3, for example). Diagnosing the model's performance yields important insights into how to make progress on modeling soil evaporation and heating under conditions of high temperatures and very low soil moisture content.