In this study, the heat distribution throughout the profile of unsaturated multilayered soil is determined using finite difference method while its thermal diffusivity varies with time and depth. First, the input parameters such as water content, dry density and sand content of the soil profile...