Spatial information of soil carbon storage at national and global level is essential for soil quality and environmental management. Improved knowledge on the amount and spatial distribution of the carbon stock in soils is crucial in estimating changes in the terrestrial carbon dynamics and management options for carbon-storing. A study was conducted to map the soil organic carbon stock (SOC) over 56,763 km2 area of Western Ghats of south India using a digital soil mapping approach. Landsat data, terrain attributes, and bioclimatic variables were used as covariates. Equal-area quadratic splines were fitted to soil profile datasets to estimate soil organic carbon stock at six standard soil depths (0–5, 5–15, 15–30, 30–60, 60–100 and 100–200 cm) and Quantile Regression Forest (QRF) algorithm was used to predict the SOC stocks. Prediction of SOC stock was better for surface layer (R2 = 31–43%) and the performance was decreasing with depth (R2 = 7–21%). The modal performance was also compared with SoilGrids products. Although the spatial patterns were similar, the present predicted SOC maps outperformed SoilGrids products in terms of both R2 and RMSE. The predicted total soil organic stock in the Western Ghats ranged from 7.1 kg m−2 to 30.9 kg m−2 and the total estimated SOC was 917 Tg. The present high resolution SOC maps help to assess and monitor the soil health and preparation of proper land use planning.