VIS-NIR spectroscopy and environmental factors coupled with PLSR models to predict soil organic carbon and nitrogen

Soil profile organic carbon (OC) and total nitrogen (TN) are influenced by topographic attributes, and land use. The visible and near-infrared (Vis-NIR) spectroscopy method can be used for the prediction of OC and TN because it is reliable, nondestructive, fast, and cost-effective. VIS-NIR soil spectral and environmental data were combined with the Partial least squares regression (PLSR) model to examine the effect of topography attributes and land use on topsoil and subsoil OC and TN stocks. After this, based on the soil depth, 114 soil samples were collected from 0 to 20 cm (topsoil) and 20–50 cm (subsoil) under three land uses, as well as OC and TN, along with several soil properties including soil particles (sand, silt, clay), pH, and bulk density in both topsoil and subsoil samples were measured. A DEM with a resolution of 30 m was used to derive the topography factors and remote sensing data was used to calculate the vegetation index. Soils (0–50 cm) under orchard land use had the highest stock of SOC (7.4 kg m−2) as well as TN (2.4 kg m−2). There was a significant increase in the organic matter stock of soils located on the south aspect (8.3 kg m−2) compared to soils located on other aspects, particularly on the north aspect (3.9% increase). Soils on the south aspect contain higher soil-water contents and lower temperatures, resulting in a decrease in the decomposition of soil organic matter. A strong positive correlation was demonstrated between topography wetness index (0.57–0.63) and topography TN stocks (0.54–0.66) as well as the highest loading score among terrain attributes, suggesting that topography is the primary factor controlling SOC stocks, particularly subsoil stocks. Additionally, we found that soils on the south-facing aspects (N aspects) had the highest spectra. Additionally, the PLSR, which showed an R2 of 0.82, a RMSE of 0.15 %, and a RPD of 0.39 indicated excellent prediction capabilities for the OC content. We concluded that the PLSR model coupled with Vis-NIR spectroscopy is able to predict topsoil and subsoil OC and N content under different aspect slopes. •Topography attributes and land use on topsoil and subsoil OC and TN.•Orchard land use had the largest SOC (7.4 kg m−2), and TN (2.4 kg m−2) stocks.•TWI was the main topographic feature controlling OC stocks.•North aspects had the highest spectra compared to other aspects.•PLSR with R2 = 0.82 showed excellent prediction for OC content.