[Predicting soil property in hilly regions by using landscape and multiscale micro-landform features]

To assess the high-resolution digital soil mapping method for small watersheds in hilly areas, we explored the role of landscape classification and multiscale micro-landform features in predicting soil pH, soil clay content (SCC), and cation exchange capacity (CEC). Geomorphons (GM) terrain classification method was used to create landform units. The traditional digital elevation model (DEM) derivatives and remote sensing variables were employed for different combinations with landscape and micro-landform classification variables, with further compa-rison and analysis being conducted. In addition, three machine learning techniques, including support vector machine (SVM), partial least squares regression (PLSR), and random forest (RF), were used to build prediction models. The best method was then selected, and then combined with regression kriging by modeling spatial structure of the model residuals. The results showed that the application of landscape and multiscale micro-landform classification variables effectively improved the prediction accuracy of pH, SCC, and CEC by 18.8%, 8.2% and 8.7%, respectively. The map of landscape classification that contained vegetation coverage information had greater model contribution than land use data. The GM classification map with 5 m resolution was more suitable for high-precision DSM than those with lower resolution. The composite model of RF performed the best in predicting SCC, while the pH and CEC were not suitable for adding the residual regression kriging on the basis of RF model. Finally, the combination of landscape and multiscale micro-landform classification variables, DEM derivatives and remote sensing variables had the highest prediction accuracy for all the three soil properties. This result indicated that multivariable contained more effective soil information than single data source for rolling areas. The landscape variables composed of GM and surface classified data explained about 40% of the spatial variation of tested soil attributes in hilly area. Therefore, multi-resolution GM and landscape classified variables could be included into the construction of prediction model in research of soil mapping.