Soil moisture represents a vital component of the ecosystem, sustaining life-supporting activities at micro and mega scales. It is a highly required parameter that may vary significantly both spatially and temporally. Due to this fact, its estimation is challenging and often hard to obtain especially over large, heterogeneous surfaces. This study aimed at comparing the performance of four widely used interpolation methods in estimating soil moisture using GPS-aided information and remote sensing. The DistanceWeighting (IDW), Spline, Ordinary Kriging models and Kriging with External Drift (KED) interpolation techniques were employed to estimate soil moisture using 82 soil moisture field-measured values. Of those measurements, data from 54 soil moisture locations were used for calibration and the remaining data for validation purposes. The study area selected was Varanasi City, India covering an area of 1535 km2. The soil moisture distribution results demonstrate the lowest RMSE (root mean square error, 8.69\%) for KED, in comparison to the other approaches. For KED, the soil organic carbon information was incorporated as a secondary variable. The study results contribute towards efforts to overcome the issue of scarcity of soil moisture information at local and regional scales. It also provides an understandable method to generate and produce reliable spatial continuous datasets of this parameter, demonstrating the added value of geospatial analysis techniques for this purpose. © 2019 by the authors.