Spatial distribution of soil data is a current issue with which most of areas are being faced. Imba-Numa watershed, Japan, is also one among those areas which lack of spatial data of soil parameters. This study aims to generate spatially distribution of soil properties which consist of soil particles, total carbon, total nitrogen and Bulk density in Imba-Numa watershed. Three geostatistical interpolation methods: Ordinary Kriging (OK), Universal Kriging (UK), and Inverse Distance Weighting (IDW), were applied to interpolate soil properties into spatially continuous data. To evaluate the performance of methods to obtain the best method which has the minimum error, Root Mean Square Error (RMSE), Mean Absolute Error (MAE), and Mean Error Percentage (MEPE) were used. Cross Validation comparison was also used to compare and validate method. The results show that Universal Kriging with Hole-effect Model is the best method to interpolate soil parameters (except Bulk density) in Imba-Numa watershed. Due to limitation of available data, the interpolated maps could not perfectly provide the satisfactory map and consists of error, but it is still acceptable.