Exploring the Influencing Factors in Identifying Soil Texture Classes Using Multitemporal Landsat-8 and Sentinel-2 Data

Soil texture is a key soil property driving physical, chemical, biological, and hydrological processes in soils. The rapid development of remote sensing techniques shows great potential for mapping soil properties. This study highlights the effectiveness of multitemporal remote sensing data in identifying soil textural class by using retrieved vegetation properties as proxies of soil properties. The impacts of sensors, modeling resolutions, and modeling techniques on the accuracy of soil texture classification were explored. Multitemporal Landsat-8 and Sentinel-2 images were individually acquired at the same time periods. Three satellite-based experiments with different inputs, i.e., Landsat-8 data, Sentinel-2 data (excluding red-edge parameters), and Sentinel-2 data (including red-edge parameters) were conducted. Modeling was carried out at three spatial resolutions (10, 30, 60 m) using five machine-learning (ML) methods: random forest, support vector machine, gradient-boosting decision tree, categorical boosting, and super learner that combined the four former classifiers based on the stacking concept. In addition, a novel SHapley Addictive Explanation (SHAP) technique was introduced to explain the outputs of the ML model. The results showed that the sensors, modeling resolutions, and modeling techniques significantly affected the prediction accuracy. The models using Sentinel-2 data with red-edge parameters performed consistently best. The models usually gave better results at fine (10 m) and medium (30 m) modeling resolutions than at a coarse (60 m) resolution. The super learner provided higher accuracies than other modeling techniques and gave the highest values of overall accuracy (0.8429), kappa (0.7611), precision (0.8378), recall rate (0.8393), and F1-score (0.8398) at 30 m with Sentinel-2 data involving red-edge parameters. The SHAP technique quantified the contribution of each variable for different soil textural classes, revealing the critical roles of red-edge parameters in separating loamy soils. This study provides comprehensive insights into the effective modeling of soil properties on various scales using multitemporal optical images.