Soil heating and evaporation under extreme conditions: Forest fires and slash pile burns

Heating any soil during a sufficiently intense wild fire or prescribed burn can alter soil irreversibly, resulting in many significant and well known, long term biological, chemical, and hydrological effects. To better understand how fire impacts soil, especially considering the increasing probability of wildfires that is being driven by climate change and the increasing use of prescribe 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, improving understanding of heat and mass transport during such extreme conditions should also provide insights into the associated transport mechanisms under more normal conditions as well. Here I describe the development of a new model designed to simulate soil heat and moisture transport during fires where the surface heating often ranges between 10,000 and 100,000 Wm^-2 for several minutes to several hours. Model performance is tested against laboratory measurements of soil temperature and moisture changes at several depths during controlled heating events created with an extremely intense radiant heater. The laboratory tests employed well described soils with well known physical properties. The model, on the other hand, is somewhat unusual in that it employs formulations for temperature dependencies of the soil specific heat, thermal conductivity, and the water retention curve (relation between soil moisture and soil moisture potential). It also employs a new formulation for the surface evaporation rate as a component of the upper boundary condition, as well as the Newton-Raphson method and the generalized Thomas algorithm for inverting block tri-diagonal matrices to solve for soil temperature and soil moisture potential. Model results show rapid evaporation rates with significant vapor transfer not only to the free atmosphere above the soil, but to lower depths of the soil, where the vapor re-condenses ahead of the heating front. Consequently the trajectory of the solution (soil volumetric water content versus soil temperature) is very unusual and highly nonlinear, which may explain why more traditional methods (i.e., those based on finite difference or finite element approaches) tend to show more numerical instabilities than the Newton-Raphson method when used to model these extreme conditions. But, despite the intuitive and qualitative appeal of the model's numerical solution, it underestimates the rate of soil moisture loss observed during the laboratory trials, although the soil temperatures are reasonably well simulated.