Quickly and correctly mapping soil nutrients significantly impact accurate fertilization, food security, soil productivity, and sustainable agricultural development. We evaluated the potential of the new PRISMA hyperspectral sensor for mapping soil organic matter (SOM), available soil phosphorus (P₂O₅), and potassium (K₂O) content over a cultivated area in Khouribga, northern Morocco. These soil nutrients were estimated using (i) the random forest (RF) algorithm based on feature selection methods, including feature subset evaluation and feature ranking methods belonging to three categories (i.e., filter, wrapper, and embedded techniques), and (ii) 107 soil samples taken from the study area. The results show that the RF-embedded method produced better predictive accuracy compared with the filter and wrapper methods. The model for SOM showed moderate accuracy ($R_{val}^2$ = 0.5, RMSEP = 0.43%, and RPIQ = 2.02), whereas that for soil P₂O₅ and K₂O exhibited low efficiency ($R_{val}^2$ = 0.26 and 0.36, RMSEP = 51.07 and 182.31 ppm, RPIQ = 0.65 and 1.16, respectively). The interpolation of RF-residuals by ordinary kriging (OK) methods reached the highest predictive results for SOM ($R_{val}^2$ = 0.69, RMSEP = 0.34%, and RPIQ = 2.56), soil P₂O₅ ($R_{val}^2$ = 0.44, RMSEP = 44.10 ppm, and RPIQ = 0.75), and soil K₂O ($R_{val}^2$ = 0.51, RMSEP = 159.29 ppm, and RPIQ = 1.34), representing the best fitting ability between the hyperspectral data and soil nutrients. The result maps provide a spatially continuous surface mapping of the soil landscape, conforming to the pedological substratum. Finally, the hyperspectral remote sensing imagery can provide a new way for modeling and mapping soil fertility, as well as the ability to diagnose nutrient deficiencies.

Integrating satellite data at different resolutions (i.e., spatial, spectral, and temporal) can be a helpful technique for acquiring soil information from a synoptic point of view. This study aimed to evaluate the advantage of using satellite mono- and multi-sensor image fusion based on either spectral indices or entire spectra to predict the topsoil clay content. To this end, multispectral satellite images acquired by various sensors (i.e., Landsat-5 Thematic Mapper (TM), Landsat-8 Operational Land Imager (OLI), Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), and Sentinel2-MultiSpectral Instrument (S2-MSI)) have been used to assess their potential in identifying bare soil pixels over an area in northeastern Tunisia, the Lebna and Chiba catchments. A spectral index image and a spectral bands image are generated for each satellite sensor (i.e., TM, OLI, ASTER, and S2-MSI). Then, two multi-sensor satellite image fusions are generated, one from the spectral index images and the other from spectral bands. The resulting spectral index and spectral band images based on mono-and multi-sensor satellites are compared through their spectral patterns and ability to predict the topsoil clay content using the Multilayer Perceptron with backpropagation learning algorithm (MLP-BP) method. The results suggest that for clay content prediction: (i) the spectral bands’ images outperformed the spectral index images regardless of the used satellite sensor; (ii) the fused images derived from the spectral index or bands provided the best performances, with a 10% increase in the prediction accuracy; and (iii) the bare soil images obtained by the fusion of many multispectral sensor satellite images can be more beneficial than using mono-sensor images. Soil maps elaborated via satellite multi-sensor data fusion might become a valuable tool for soil survey, land planning, management, and precision agriculture.

The potential of Visible Near-Infrared and Short-Wave Infrared (VNIR–SWIR, 400 nm–2500 nm) hyperspectral imagery for use in multivariate approaches and geostatistical techniques for mapping topsoil properties has been previously demonstrated. However, the use of VNIR–SWIR hyperspectral imagery remains costly, which limits the spatial scales over which it can be applied. This paper aims to evaluate the potential for substituting the more accessible Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) VNIR–SWIR multispectral data for hyperspectral imagery in mapping surface soil clay contents. This study used ASTER multispectral data (nine bands in the VNIR–SWIR spectral domain) acquired over the Cap-Bon region in Tunisia (2000 km²) and 262 surface soil samples collected within the ASTER scene that were subjected to laboratory analysis of the clay fraction (soil particles less than 2 µm). The approach followed two steps: i) estimation of surface soil clay contents for bare soil areas using a Multiple Linear Regression (MLR) model built from the 9 ASTER VNIR–SWIR bands and ii) spatial interpolation (co-kriging) of the soil sampling of measured points and the ASTER-estimates over the whole study area. The MLR model for estimating clay contents using ASTER multispectral data performed correctly ( $R_{\mathrm{v}\mathrm{a}\mathrm{l}}^2$ = 0.60). In addition, this performance is only slightly lower than that obtained using hyperspectral imagery (specifically, an Airborne Imaging Spectrometer for Applications (AISA-DUAL) dual hyperspectral sensor) in a previous study. Moreover, the co-kriging process appeared to yield encouraging results for capturing the large range of variability of clay content values, although it was not able to represent the short scale variability ( $R_{\mathrm{C}\mathrm{V}}^2$ = 0.43). Finally, the ASTER multispectral data, despite being underused in the mapping of soil properties, may open up new ways to perform more extensive mapping of surface soil properties in semi-arid contexts characterized by extensive bare and dry soil surfaces.

L'article présente les résultats de l'étude de l'évolution spatiale et temporelle de l'occupation du sol dans le bassin de l'Oued Chiba (200 km2) situé au Nord-Est de la Tunisie. Cette analyse est fondée sur une classification de données de télédétection (images Landsat TM des années 1992, 2011 et une image ASTER de l'année 2001) et de photographies aériennes à l'échelle de 1/20 000 (des années 1962, 1974 et 1984). La méthode de classification des images satellites utilisées est un classifieur supervisé de type séparateur à vastes marges (SVM) avec un noyau gaussien.L'analyse de l'occupation des sols sur une période de 49 ans se traduit par une augmentation de l'espace agricole (classe des cultures / jachères et la classe d'arboricultures) de 24 km2 en 1962 à 75 km2 en 2011. Ces changements d'occupation et d'utilisation du sol sont liés directement à la création du barrage de l'Oued Chiba en 1963.