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.

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.

Groundwater is a most important resource in arid and semi-arid regions and is required for drinking, irrigation and industrialization. Assessing the potential zone of groundwater recharge is extremely crucial for the protection of water quality and the management of groundwater systems. To identify the groundwater potential zone in the study area, thematic layers of lithology, slope, karst degrees, land cover, lineament and drainage density were generated using topographic maps, thematic maps, field data and satellite image, and were prepared, classified, weighted and integrated in a geographic information system (GIS) environment by the means of fuzzy logic. The fuzzy membership values have been assigned to different thematic layers according to their classification on respect for their contribution and their occurrence in groundwater. Based on the generated groundwater potential map, it was found that about 8% of the investigation area was categorized as very high potential for groundwater recharge, 31% as high, 28% as moderate, 17% as low and 16% as very low potential for groundwater recharge. Finally, the results were verified using well-yield data. The highest recharge potential area is located towards the downstream regions related to more fractured and karstified limestone.

Groundwater pollution poses a severe threat and issue to the environment and humanity overall. That is why mitigative strategies are urgently needed. Today, studies mapping groundwater risk pollution assessment are being developed. In this study, five new hybrid/ensemble machine learning (ML) models are developed, named DRASTIC-Random Forest (RF), DRASTIC-Support Vector Machine (SVM), DRASTIC-Multilayer Perceptron (MLP), DRASTIC-RF-SVM, and DRASTIC-RF-MLP, for groundwater pollution assessment in the Saiss basin, in Morocco. The performances of these models are evaluated using the Receiver Operating Characteristic curve (ROC curve), precision, and accuracy. Based on the results of the ROC curve method, it is indicated that the use of hybrid/ensemble machine learning (ML) models improves the performance of the individual machine learning (ML) algorithms. In effect, the AUC value of the original DRASTIC is 0.51. Furthermore, both hybrid/ensemble models, DRASTIC-RF-MLP (AUC = 0.953) and DRASTIC-RF-SVM, (AUC = 0.901) achieve the best accuracy among the other models, followed by DRASTIC-RF (AUC = 0.852), DRASTIC-SVM (AUC = 0.802), and DRASTIC-MLP (AUC = 0.763). The results delineate areas vulnerable to pollution, which require urgent actions and strategies to improve the environmental and social qualities for the local population.

Gully erosion is a serious threat to the state of ecosystems all around the world. As a result, safeguarding the soil for our own benefit and from our own actions is a must for guaranteeing the long-term viability of a variety of ecosystem services. As a result, developing gully erosion susceptibility maps (GESM) is both suggested and necessary. In this study, we compared the effectiveness of three hybrid machine learning (ML) algorithms with the bivariate statistical index frequency ratio (FR), named random forest-frequency ratio (RF-FR), support vector machine-frequency ratio (SVM-FR), and naïve Bayes-frequency ratio (NB-FR), in mapping gully erosion in the GHISS watershed in the northern part of Morocco. The models were implemented based on the inventory mapping of a total number of 178 gully erosion points randomly divided into 2 groups (70% of points were used for training the models and 30% of points were used for the validation process), and 12 conditioning variables (i.e., elevation, slope, aspect, plane curvature, topographic moisture index (TWI), stream power index (SPI), precipitation, distance to road, distance to stream, drainage density, land use, and lithology). Using the equal interval reclassification method, the spatial distribution of gully erosion was categorized into five different classes, including very high, high, moderate, low, and very low. Our results showed that the very high susceptibility classes derived using RF-FR, SVM-FR, and NB-FR models covered 25.98%, 22.62%, and 27.10% of the total area, respectively. The area under the receiver (AUC) operating characteristic curve, precision, and accuracy were employed to evaluate the performance of these models. Based on the receiver operating characteristic (ROC), the results showed that the RF-FR achieved the best performance (AUC = 0.91), followed by SVM-FR (AUC = 0.87), and then NB-FR (AUC = 0.82), respectively. Our contribution, in line with the Sustainable Development Goals (SDGs), plays a crucial role for understanding and identifying the issue of “where and why” gully erosion occurs, and hence it can serve as a first pathway to reducing gully erosion in this particular area.