Cross-Regional Digital Soil Mapping of OpenLandMap-Derived USDA Soil Texture Classes Using Remote Sensing and Machine Learning

Soil is a fundamental environmental resource that supports food production and sustainable agricultural systems. Soil’s texture is a key soil property that affects its capacity to hold water, infiltration rate, fertility, erosion risk, and crop productivity. Traditionally, it is determined using hydrometer or pipette analysis. These laboratory procedures are costly and time-consuming, typically requiring approximately two working days, which limits their practicality for large-scale assessments. This paper proposes a cross-regional digital soil texture classification system based on remote sensing data and machine learning in Malaysia and Pakistan. The two countries were chosen because they represent contrasting environmental conditions: Malaysia has humid tropical conditions, and Pakistan has arid and semi-arid conditions. The OpenLandMap United States Department of Agriculture (USDA) soil texture classes served as the reference labels, and the spectral bands from Landsat, together with spectral indices and Shuttle Radar Topography Mission (SRTM) elevation, were used as predictor variables. The three ensemble models, including Random Forest (RF), XGBoost, and CatBoost, were evaluated for multi-class classification of USDA soil texture. The main Landsat-based Malaysia–Pakistan dataset produced the best results, with RF achieving the highest accuracy of 71.45%, XGBoost 71.31%, and CatBoost 67.33%. Additional experiments using 2025 Landsat 8/9 and Sentinel-2 imagery achieved lower accuracies of 63.00% and 56.10%, respectively, while five-fold spatial block cross-validation produced a mean accuracy of 59.27 ± 1.84%. To strengthen dataset reliability, representative OpenLandMap-derived soil texture classes were externally validated using PCRWR sand–silt–clay fractions plotted on the USDA soil texture triangle and independent Malaysian soil resource references. The results of the feature importance analysis showed that elevation, Bare Soil Index (BSI), Normalized Difference Water Index (NDWI), and Normalized Difference Vegetation Index (NDVI) were important predictors for classifying OpenLandMap-derived USDA soil texture classes. These show that the Landsat-based spectral and terrain features, when used in conjunction with ensemble learning, could assist cross-regional soil texture classification for the USDA. However, class imbalance, spectral overlap between soil texture classes, temporal variability, and resolution mismatch between satellite predictors and soil labels remain important challenges. The findings of this study provide a useful foundation for large-scale soil texture mapping and future spatial validation across different environmental regions.