Soil Suborder Discrimination Using Machine Learning Is Improved by SWIR Imaging Compared with Full VIS–NIR–SWIR Spectra

Rapid, standardised discrimination of soil taxonomic units remains challenging when relying solely on conventional field descriptions and laboratory analyses, particularly at high sampling densities. This study evaluated whether proximal spectroscopy and hyperspectral imaging can support the classification of Brazilian Soil Classification System (SiBCS) suborders and pedogenetic horizons when surface and subsurface spectra are treated separately. Six intact soil monoliths (0.12 × 1.60 m) were collected in Paraná State, southern Brazil, representing one Organossolo (Ooy), three Latossolos (LVd, LVd1, and LVd2) and two Argissolos (PVAd and PVd). For each monolith, 800 spectra were acquired per sensor with a non-imaging VIS–NIR–SWIR spectroradiometer (350–2500 nm), and 800 spectra per sensor per monolith were extracted from the SWIR hyperspectral images (1200–2450 nm). Principal component analysis (PCA) was used to summarise spectral variability, and supervised classification was performed via k-nearest neighbours, random forest, decision tree and gradient boosting for suborders (10-fold cross-validation), and a neural network was used for within-profile horizon classification. PCA indicated that most of the spectral variance was captured by a dominant axis, with clearer separation among suborders in the SWIR space than in the full VIS–NIR–SWIR range. With respect to suborder classification, subsurface spectra outperformed surface spectra, and SWIR outperformed VIS–NIR–SWIR: the best accuracies were 0.96 for subsurface SWIR (gradient boosting; AUC = 0.99; MCC = 0.95) and 0.89 for surface SWIR (k-nearest neighbours; AUC = 0.98; MCC = 0.87). Within-profile horizon classification via VIS–NIR–SWIR achieved accuracies of 0.84–0.97 with the Neural Network, with most misclassifications occurring between adjacent horizons. Overall, subsurface SWIR information provided the most reliable basis for taxonomic discrimination, whereas horizon classification was feasible but reflected gradual spectral transitions along the profile.