The study of land use/land cover (LULC) has become an increasingly important stage in the development of forest ecosystems strategies. Hence, the main goal of this study was to describe the vegetation change of Azrou Forest in the Middle Atlas, Morocco, between 1987 and 2017. To achieve this, a set of Landsat images, including one Multispectral Scanner (MSS) scene from 1987; one Enhanced Thematic Mapper Plus (ETM+) scene from 2000; two Thematic Mapper (TM) scenes from 1995 and 2011; and one Landsat 8 Operational Land Imager (OLI) scene from 2017; were acquired and processed. Ground-based survey data and the normalized difference vegetation index (NDVI) were used to identify and to improve the discrimination between LULC categories. Then, the maximum likelihood (ML) classification method was applied was applied, in order to produce land cover maps for each year. Three classes were considered by the classification of NDVI value: low-density vegetation; moderate-density vegetation, and high-density vegetation. Our study achieved classification accuracies of 66.8% (1987), 99.9% (1995), 99.8% (2000), 99.9% (2011), and 99.9% (2017). The results from the Landsat-based image analysis show that the area of low-density vegetation was decreased from 27.4% to 2.1% over the past 30 years. While, in 2017, the class of high-density vegetation was increased to 64.6% of the total area of study area. The results of this study show that the total forest cover remained stable. The present study highlights the importance of the image classification algorithms combined with NDVI index for better understanding the changes that have occurred in this forest. Therefore, the findings of this study could assist planners and decision-makers to guide, in a good manner, the sustainable land development of areas with similar backgrounds.

Machine learning (ML) models are commonly used in solar modeling due to their high predictive accuracy. However, the predictions of these models are difficult to explain and trust. This paper aims to demonstrate the utility of two interpretation techniques to explain and improve the predictions of ML models. We compared first the predictive performance of Light Gradient Boosting (LightGBM) with three benchmark models, including multilayer perceptron (MLP), multiple linear regression (MLR), and support-vector regression (SVR), for estimating the global solar radiation (H) in the city of Fez, Morocco. Then, the predictions of the most accurate model were explained by two model-agnostic explanation techniques: permutation feature importance (PFI) and Shapley additive explanations (SHAP). The results indicated that LightGBM (R² = 0.9377, RMSE = 0.4827 kWh/m², MAE = 0.3614 kWh/m²) provides similar predictive accuracy as SVR, and outperformed MLP and MLR in the testing stage. Both PFI and SHAP methods showed that extraterrestrial solar radiation (H₀) and sunshine duration fraction (SF) are the two most important parameters that affect H estimation. Moreover, the SHAP method established how each feature influences the LightGBM estimations. The predictive accuracy of the LightGBM model was further improved slightly after re-examination of features, where the model combining H₀, SF, and RH was better than the model with all features.

Soil erosion is an increasingly issue worldwide, due to several factors including climate variations and humans’ activities, especially in Mediterranean ecosystems. Therefore, the aim of this paper is: (i) to quantify and to predict soil erosion rate for the baseline period (2000–2013) and a future period (2014–2027), using the Revised Universal Soil Loss Equation (RUSLE) and the Soil and Water Assessment Tool (SWAT) model in the R’Dom watershed in Morocco, based on the opportunities of Remote Sensing (RS) techniques and Geographical Information System (GIS) geospatial tools. (ii) we based on classical statistical downscaling model (SDSM) for rainfall prediction. Due to the lack of field data, the model results are validated by expert knowledge. As a result of this study, it is found that both agricultural lands and bare lands are most affected by soil erosion. Moreover, it is showed that soil erosion in the watershed was dominated by very low and low erosion. Although the area of very low erosion and low erosion continued to decrease. Hence, we hereby envisage that our contribution will provide a more complete understanding of the soil degradation in this study area and the results of this research could be a crucial reference in soil erosion studies and also may serve as a valuable guidance for watershed management strategies.

Calcul De L’indice De Qualité De L’eau (IQE) Pour L’évaluation De La Qualité Physico-Chimique Des Eaux Superficielles De L’Oued Moulouya (NE, Maroc)