Most models used in land surface hydrology, vadose zone hydrology, and hydro-climatology require an accurate representation of soil thermal properties (soil thermal conductivity and volumetric heat capacity). Various empirical relations have been suggested to estimate soil thermal properties. However, they require many input parameters such as soil texture, mineralogical composition, porosity and water content, which are not always available from laboratory experiments and field measurements. In this paper, to overcome the above challenge, a hybrid numerical method, Genetic Algorithm- Finite Difference (GA-FD), is proposed to estimate soil thermal properties using land surface temperature (LST) as the only input. The genetic algorithm (GA) optimization method coupled with the finite difference (FD) modeling technique is a viable hybrid approach for estimating soil thermal properties. The finite difference method is employed to solve the heat diffusion equation and simulate LST, while a robust optimization technique (GA) is used to retrieve soil thermal properties by minimizing the difference between observed and simulated LST. Furthermore, a generalization of the hybrid model is developed for inhomogeneous soil, in which soil thermal properties are not constant throughout the soil slab. The proposed model is applied to the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE). The results show that the proposed hybrid numerical method is able to estimate soil thermal properties accurately, and therefore effectively eliminate the need for the unavailable soil parameters which are required by empirical methods for determining the soil thermal conductivity and volumetric heat capacity. Remarkably, the temporal variation of the retrieved soil thermal conductivity is consistent with the volumetric water content, even though no water content information is used in the model.