Changes in soil temperature (ST) play an important role in the main mechanisms within the soil, including biological and chemical activities. For instance, they affect the microbial community composition, the speed at which soil organic matter breaks down and becomes minerals. Moreover, the growth and physiological activity of plants are directly influenced by the ST. Additionally, ST indirectly affects plant growth by influencing the accessibility of nutrients in the soil. Therefore, designing an efficient tool for ST estimating at different depths is useful for soil studies by considering meteorological parameters as input parameters, maximal air temperature, minimal air temperature, maximal air relative humidity, minimal air relative humidity, precipitation, and wind speed. This investigation employed various statistical metrics to evaluate the efficacy of the implemented models. These metrics encompassed the correlation coefficient (r), root mean square error (RMSE), Nash-Sutcliffe (NS) efficiency, and mean absolute error (MAE). Hence, this study presented several artificial intelligence-based models, MLPANN, SVR, RFR, and GPR for building robust predictive tools for daily scale ST estimation at 05, 10, 20, 30, 50, and 100cm soil depths. The suggested models are evaluated at two meteorological stations (i.e., Sulaimani and Dukan) located in Kurdistan region, Iraq. Based on assessment of outcomes of this study, the suggested models exhibited exceptional predictive capabilities and comparison of the results showed that among the proposed frameworks, GPR yielded the best results for 05, 10, 20, and 100cm soil depths, with RMSE values of 1.814°C, 1.652°C, 1.773°C, and 2.891°C, respectively. Also, for 50cm soil depth, MLPANN performed the best with an RMSE of 2.289°C at Sulaimani station using the RMSE during the validation phase. Furthermore, GPR produced the most superior outcomes for 10cm, 30cm, and 50cm soil depths, with RMSE values of 1.753°C, 2.270°C, and 2.631°C, respectively. In addition, for 05cm soil depth, SVR achieved the highest level of performance with an RMSE of 1.950°C at Dukan station. The results obtained in this research confirmed that the suggested models have the potential to be effectively used as daily predictive tools at different stations and various depths.