There is a great challenge in combining soil proximal spectra and remote sensing spectra to improve the accuracy of soil organic carbon (SOC) models. This is primarily because mixing of spectral data from different sources and technologies to improve soil models is still in its infancy. The first objective of this study was to integrate information of SOC derived from visible near-infrared reflectance (Vis-NIR) spectra in the laboratory with remote sensing (RS) images to improve predictions of topsoil SOC in the Skjern river catchment, Denmark. The second objective was to improve SOC prediction results by separately modeling uplands and wetlands. A total of 328 topsoil samples were collected and analyzed for SOC. Satellite Pour l’Observation de la Terre (SPOT5), Landsat Data Continuity Mission (Landsat 8) images, laboratory Vis-NIR and other ancillary environmental data including terrain parameters and soil maps were compiled to predict topsoil SOC using Cubist regression and Bayesian kriging. The results showed that the model developed from RS data, ancillary environmental data and laboratory spectral data yielded a lower root mean square error (RMSE) (2.8%) and higher R² (0.59) than the model developed from only RS data and ancillary environmental data (RMSE: 3.6%, R²: 0.46). Plant-available water (PAW) was the most important predictor for all the models because of its close relationship with soil organic matter content. Moreover, vegetation indices, such as the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI), were very important predictors in SOC spatial models. Furthermore, the ‘upland model’ was able to more accurately predict SOC compared with the ‘upland & wetland model’. However, the separately calibrated ‘upland and wetland model’ did not improve the prediction accuracy for wetland sites, since it was not possible to adequately discriminate the vegetation in the RS summer images. We conclude that laboratory Vis-NIR spectroscopy adds critical information that significantly improves the prediction accuracy of SOC compared to using RS data alone. We recommend the incorporation of laboratory spectra with RS data and other environmental data to improve soil spatial modeling and digital soil mapping (DSM).

After decades of mining and industrialization in Qatar, it is important to estimate their impact on soil pollution with toxic metals. The study utilized 300 topsoil (0–30 cm) samples, multi-spectral images (Landsat 8), spectral indices and environmental variables to model and map the spatial distribution of arsenic (As), chromium (Cr), nickel (Ni), copper (Cu), lead (Pb) and zinc (Zn) in Qatari soils. The prediction model used condition-based rules generated in the Cubist tool. In terms of R² and the ratio of performance to interquartile distance (RPIQ), the models showed good predictive capabilities for all elements. Of all of the prediction results, Cu had the highest R² = 0.74, followed by As > Pb > Cr > Zn > Ni. This study found that all of the models only chose images from January and February as predictors, which indicates that images from these two months are important for soil toxic metals’ monitoring in arid soils, due to the climate and the vegetation cover during this season. Topsoil maps of the six toxic metals were generated. The maps can be used to prioritize the choice of remediation measures and can be applied to other arid areas of similar environmental/socio-economic conditions and pollution causes.