In this work, we consider simulation of heat fl w in the shallow subsurface. As sunlight heats up the surface of soil, the thermal energy received dissipates downward into the ground. This process can be modeled using a partial differential equation known as the heat equation. The spatial distribution of soil thermal conductivities is a key factor in the modeling process. Prior to this study, t...