Soil information, including their spatial variability is vital for devising various soil and land management policies and strategies as well as assessing the environmental impact of different land use changes. Among the various soil forming factors, terrain plays a vital role in determining the spatial distribution of soil properties especially in hilly terrains due to its impact on soil processes as well as the redistribution of soil materials. The study was carried out for assessing the applicability and efficiency of various terrain parameters for prediction as well as spatial mapping of different soil properties using Random Forest (RF) modeling approach, in a mid-Himalayan watershed located in Tehri Garhwal district, Uttarakhand. Field data collection was done by collecting soil samples from 0-15cm depth (surface layer) at 68 different locations in the study area, following a transect survey method. Collected samples were pre-processed and analyzed in the laboratory for estimating various soil properties such as sand, silt and clay contents as well as soil nutrients such as Soil Organic Carbon (SOC) and nitrogen (N). Fifteen terrain parameters derived by digital terrain analysis of CartoDEM of 10m spatial resolution were used as predictor variables during RF modeling. Pearson correlation analysis revealed poor linear relationship of soil properties with various terrain parameters. Variable importance ranking revealed that compound topographic index, elevation, heat load index, integrated moisture index and mean curvature are the most important variables, explaining maximum variability in different properties. The performance of RF models for predicting various soil properties were evaluated based on the coefficient of determination (R²), Root Mean Square Error (RMSE) as well as Lin’s Concordance Correlation Coefficient (CCC). The models were found to perform well with R² and CCC values of 0.47 and 0.50 for sand, 0.52 and 0.55 for silt, 0.44 and 0.57 for clay, 0.41 and 0.36 for organic carbon and 0.55 and 0.34 for nitrogen, respectively. The validated models were further employed for generating the spatial distribution maps of various soil properties. Such attempts combining various environmental variables with machine learning based models may aid us in mapping the spatial distribution of various soil properties in the hilly mountainous terrains with considerable accuracies.