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 are provided. These data are coupled with the theoretical approaches to estimate thermal properties of soil such as thermal conductivity and thermal diffusivity of multilayered soil. Second, finite difference method is used to model heat distributions in soil profile taking into account the initial and boundary conditions. A continuity of heat flux between each layer is performed as a condition in the numerical model. A comparison of estimated temperature within time throughout the profile with the thermal probe measurements shows a satisfactory capacity of the numerical model. Finally, different cases of nonhomogeneous and homogeneous soil show that thermal response of homogeneous and nonhomogeneous soils are almost similar at average value of thermal diffusivity where hydrothermal characteristics of each soil layer (such as water content, dry density, and soil texture) are required to calculate this average value. (C) 2014 Elsevier Ltd. All rights reserved.