Next Article in Journal
Empirical Study of PEFT Techniques for Winter-Wheat Segmentation
Previous Article in Journal
Trainable Noise Model as an Explainable Artificial Intelligence Evaluation Method: Application on Sobol for Remote Sensing Image Segmentation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Assessing the Impact of Climate Change on Seasonal Variation in Agricultural Land Use Using Sentinel-2 and Machine Learning †

by
Musa Mustapha
*,‡ and
Mhamed Zineddine
School of Digital Engineering and Artificial Intelligence, Euro-Mediterranean University of Fes, Fes 30030, Morocco
*
Author to whom correspondence should be addressed.
Presented at the 5th International Electronic Conference on Remote Sensing, 7–21 November 2023; Available online: https://ecrs2023.sciforum.net/.
These authors contributed equally to this work.
Environ. Sci. Proc. 2024, 29(1), 51; https://doi.org/10.3390/ECRS2023-16365
Published: 6 November 2023
(This article belongs to the Proceedings of ECRS 2023)

Abstract

:
The Fez region in Morocco has experienced changes in agricultural land use as a result of climate change. These changes include erratic rainfall, rising temperatures, and evapotranspiration. The objective of this research is to investigate the impact of these changes on agricultural land use between 2018 and 2022 using remote sensing data (sentinel-2 and MODIS), climate data, drought index (Vegetation Condition Index (VCI)) and two machine learning algorithms (Random Forest (RF) and Gradient Tree Boost (GTB). The RF and GTB algorithms were trained and tested, and their performance was analyzed, revealing that the GTB algorithm is more efficient than the RF, with a Kaffa coefficient of 91% and overall accuracy of 93%. The analysis of climate change on land use and land cover (LULC) variations revealed a significant (54%) reduction in rainfall. Furthermore, agricultural land use and water were reduced by 41% and 17%, respectively. Conversely, barren land and built-up areas increased by 58% and 4%, respectively, and the annual mean VCI decreased from 39.72 in 2018 to 19.9 in 2022. The study concluded that climate change had a significant impact on the region’s agricultural land cover, and decreases in rainfall directly affect agricultural land use.

1. Introduction

Agriculture plays a crucial role in the Moroccan economy, contributing between 14 and 20% to the nation’s Gross Domestic Product (GDP) [1]. It serves as a major source of employment, accounting for 43% of job opportunities, particularly in rural areas, while maintaining food security [1,2]. In addition, agriculture acts as a source of raw materials for both domestic and international industries. In recent times, the agricultural sectors in the Middle East and North Africa (MENA) region have encountered a range of challenges, notably the impact of climate change on rainfall patterns and the depletion of groundwater resources [3,4]. Furthermore, salinization in soils and groundwater causes a detrimental impact on agricultural product growth and quality [5]. Incorporating decision-making tools in Moroccan agriculture can significantly assist with evaluating and tackling climate-related shifts. Employing remote sensing technology offers a quantitative approach to comprehending environmental changes, resulting in reduced expenses and time spent analyzing land use and cover. This can be a useful tool in assessing and addressing the impact of climate change on agriculture in Morocco [6,7]. The popularity of machine learning (ML) methods for remote sensing has significantly increased, especially in the classification of land use and land cover (LULC). ML techniques have shown remarkable effectiveness in reducing processing time and improving accuracy, enabling the efficient detection of changes in LULC.
When collecting remote sensing data for LULC classification, two main satellite platforms are utilized: Sentinel-2 and Landsat. Although both have their advantages, Landsat provides a more extensive dataset that goes back to 1972. However, Sentinel-2’s spatial resolution is significantly higher, at 10 meters, compared to Landsat’s 60 meters. Therefore, it is preferred to use Sentinel-2 for analyzing LULC changes from 2016 to the present. Additionally, Sentinel-2 offers 13 spectral bands, while Landsat 9 has only 9 bands that cover the visible, near-infrared (NIR) and shortwave infrared (SWIR) wavelengths. To extract relevant information on surface features such as vegetation, urban areas, soil or water, spectral indices (SI) are utilized, significantly enhancing the accuracy of LULC classification. The accuracy of LULC classification is contingent on the selection of appropriate spectral bands and indices, adequate training samples, and the application of a suitable machine-learning algorithm. Numerous scholars have thoroughly investigated the utilization of spectral bands and indices to evaluate changes in LULC, thereby facilitating comprehension of the impact of climate change on changes in urban, agricultural and water bodies. A study conducted by Beroho et al. (2021) examined the changes in land use and land cover (LULC) in a Mediterranean watershed in Morocco from 1998 to 2018. The study used a Markov Chain and Cellular Automaton (CA-Markov) model [8]. The outcomes indicate that urbanization has significantly reduced agricultural land in the Drâa Valley. Similarly, Karmaoui et al. (2021) investigated LULC changes in Errachidia province, Morocco, from 2005 to 2020 using NDVI, NDWI and EVI [9]. The findings suggest that urban development has had an adverse impact on groundwater, soil quality and natural ecosystems. The study also highlights the importance of remote sensing in detecting hydrological droughts and aiding in the conservation of land and water resources.
The objective of this research was to analyze the effects of temperature, drought and rainfall patterns on agricultural land located in the Fez prefecture of Morocco during the rainy seasons of 2018 and 2022. To classify the land-use/land-cover (LULC) map into four categories, namely builtup, agriculture, water and bareland, machine learning models were employed based on Sentinel-2 satellite data. The LULC changes were compared with climate and drought changes using MODIS, CHIRPS and AgERA5 databases. The results indicate that machine learning methods can efficiently be employed to monitor changes in agricultural land.

2. Materials and Methods

2.1. Study Overview

Figure 1 depicts the study flowchart, which delineates the sequential procedures involved in LULC classification and change analysis. Remote sensing data were initially collected from the Sentinel-2 satellite. Following that, spatial and temporal filtering techniques were used to obtain cloud-free imagery while also limiting the selection of dates within the scope of the investigation. After that, the spectral indices listed in Table 1 were computed using ground truth data. Next, Random Forest and Gradient Tree Boosting were trained with the training set. The validation set was used to evaluate the models’ performance using a variety of performance metrics. The LULC results and climate data from the CHIRPS dataset were used to analyze the changes.

2.2. Site Description

Fez is one of the oldest cities, with origins dating back to the 14th century. It is situated in the northern part of Morocco as depicted in Figure 2. Fez is renowned for its prominence in the fields of tourism and agriculture. The climate of the region is classified as semi-arid, featuring an annual precipitation of approximately 500 mm and an average temperature of 18 degrees Celsius.

2.3. Data Collection

2.3.1. Sentinel-2 and MODIS Data

Using the Google Earth Engine platform, Remote sensing imagery data from the Sentinel-2 satellite and MODIS sensor were processed for the years 2018 and 2022. The data underwent cloud filtration to ensure clarity and accuracy, and were further limited to a specific study area. This process enabled the collection of high-quality, relevant data for analysis and research purposes.

2.3.2. Climate Data

The study area’s annual historical aggregated time series data for rainfall, maximum and minimum temperature, and evapotranspiration for the years 2018 and 2022 were obtained from the CHIRPS [10] and AgERA5 [11] datasets available at (https://aquastat.fao.org/climate-information-tool/climate-data?lat=33.8843&lon=-5.7375&year=2022, accessed on 2 September 2023).

2.3.3. Ground Truth Data

We collected primary data from the field by conducting on-site visits. During these visits, we carefully recorded the geographical coordinates and assigned appropriate labels corresponding to the LULC classes (water, land, agriculture and bare land) for each specific location. Using the Google Earth Engine platform, ground data were subsequently utilized to accurately pinpoint places within the study area and generated 1126 samples (908 training and 218 validation sets).

2.3.4. Spectral Indices (SI)

Spectral indices refer to linear combinations of multiple spectral bands from remote sensing satellites, which are utilized to extract additional information and enhance the identification and differentiation of specific features on the Earth’s surface. In this study, we employed four spectral indices presented in Table 1 to augment the feature set.
Table 1. List of spectral indices used for LULC classification.
Table 1. List of spectral indices used for LULC classification.
IndicesNameFormulaReference
NDVINormalized difference vegetation index(B8 − B4)/(B8 + B4)[12]
NDWINormalized difference water index(B3 − B8)/(B3 + B8)[13]
NDBINormalized difference built-up index(B11 − B8)/(B11 + B8)[14]
BSIBare soil index(B11 − B12 − B8)/(B11 + B12 + B08)[15]

2.3.5. Drought Index (VCI)

The vegetation condition index (VCI) is a metric derived from NDVI using Equation (1) to evaluate vegetation health. It is used to classify drought severity into three classes: no drought (VCI > 50), drought (35 <= VCI <= 50) and severe drought (VCI < 35) [16].
V C I = N D V I N D V I m i n N D V I m a x N D V I m i n × 100

2.4. Classifiers

2.4.1. Random Forest (RF)

The RF algorithm is a supervised machine-learning technique that is utilized for both classification and regression tasks [17]. It falls under the category of controlled nonparametric methods. The RF algorithm utilizes an ensemble of decision trees and combines their predictions through a majority voting technique. RF has been utilized in numerous studies to address various LULC classification problems [18], consistently demonstrating superior performance and yielding optimal outcomes. The fundamental components of random forests encompass the number of trees, variables per split, bag fraction (BF), maximum nodes and minimum tree leaves.

2.4.2. Gradient Tree Boost (GTB)

The GTB model is a supervised algorithm that is commonly used to solve classification and regression problems [19]. GTB aggregates and produces a more precise final result by using an ensemble of weak individual decision trees. GTB can avoid overfitting by fitting the residuals of the regression tree at each iteration with negative gradient loss values. Several authors have reported on the use of GTB for LULC classification, revealing that the algorithm also produces good results [20,21].

2.5. Performance Evaluation Metrics

The evaluation of the machine learning algorithm’s accuracy for the classification of LULC was conducted using a confusion matrix. This assessment included metrics such as overall accuracy (OA), producer’s accuracy (PA), consumer’s accuracy (CA), and kappa coefficient.

3. Results and Discussion

3.1. Classification Accuracy

The classification performance of the RF and GTB algorithms for LULC is presented in Table 2. Both algorithms demonstrated exceptional performance in generating accurate LULC maps for the study, with an Overall Accuracy and kappa coefficient exceeding 80%.
The assessment of the classification models’ performance was based on the metrics outlined in Section 2.5. The GTB model outperformed RF in terms of overall accuracy, achieving a score of 93% and a kappa coefficient of 91%. In contrast, RF achieved an accuracy rate of 92% and a kappa coefficient of 89%, which is comparatively lower.

3.2. Analysis of LULC, Climate and VCI Changes

The best technique for examining land-cover changes and comprehend how changes occur within classes is to compare classification results in detail. Figure 3A and Table 3 provide the LULC map comparison results for 2018 and 2022. Figure 3B displays the VCI for the month of March in both 2018 and 2022. The graph clearly illustrates the noticeable effect of drought in 2022. Figure 4 illustrates the VCI time-series plot for the years 2018 and 2022, with a noticeable decrease in the trend observed in 2022.
As indicated in Section 2.3.2, meteorological data for the specified period were acquired and are presented in Table 4. It is observed that there was a significant increase in rainfall of 392 mm/year (46%) between the years 2018 and 2022. Furthermore, there was a subsequent increase in evapotranspiration of 173 mm (11.5%). The minimum and maximum temperatures both experienced increases of 0.5 and 1.4 degrees Celsius, respectively.
Based on these findings, it is evident that a reduction in rainfall has a direct impact on the decline in agricultural land coverage and VCI. Additionally, the drop in rainfall contributes to an increase in the rate of evapotranspiration and VCI. Consequently, there exists an inverse correlation between rainfall and evapotranspiration. Secondly, there exists a classification overlap between the built-up region and bareland class as a result of the resemblance in their features. The bareland class encompasses terrains such as mountains and rocks, which are utilized to produce construction materials like marble and tiles.

4. Conclusions

In this research, two classifiers, namely RF and GTB, were utilized for the purpose of LULC classification. The findings of the study indicated that the GTB model had superior performance compared to the RF model in terms of both overall accuracy and the Kaffa coefficient. Based on our findings, it can be inferred that climate change has a direct influence on agricultural land cover. Therefore, it is of utmost importance to monitor changes in agricultural land cover to identify factors contributing to degradation, particularly with regard to climate and food security. The scope of this research is limited to the analysis of LULC changes specifically during the rainfed season. Further investigation of the irrigation season is essential to extrapolate the results and assess the broader changes in agricultural land and the state of food security within the region.

Author Contributions

Conceptualisation, M.M.; methodology, investigation, M.M. and M.Z.; writing—original draft preparation, M.M.; writing—review and editing, supervision, M.M. and M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data utilized in this manuscript are open-access data and is available and publicly accessible online.

Acknowledgments

Thanks to all anonymous reviewers and editors for taking the time to review this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Epule, T.E.; Chehbouni, A.; Chfadi, T.; Ongoma, V.; Er-Raki, S.; Khabba, S.; Etongo, D.; Martínez-Cruz, A.L.; Molua, E.L.; Achli, S.; et al. A Systematic National Stocktake of Crop Models in Morocco. Ecol. Model. 2022, 470, 110036. [Google Scholar] [CrossRef]
  2. Abdelmajid, S.; Mukhtar, A.; Baig, M.B.; Reed, M.R. Climate Change, Agricultural Policy and Food Security in Morocco. In Emerging Challenges to Food Production and Security in Asia, Middle East, and Africa; Springer International Publishing: Cham, Switzerland, 2021; Book Section Chapter 7; pp. 171–196. [Google Scholar] [CrossRef]
  3. Ortega-Pozo, J.L.; Alcalá, F.J.; Poyatos, J.M.; Martín-Pascual, J. Wastewater Reuse for Irrigation Agriculture in Morocco: Influence of Regulation on Feasible Implementation. Land 2022, 11, 2312. [Google Scholar]
  4. Eddoughri, F.; Lkammarte, F.Z.; El Jarroudi, M.; Lahlali, R.; Karmaoui, A.; Yacoubi Khebiza, M.; Messouli, M. Analysis of the Vulnerability of Agriculture to Climate and Anthropogenic Impacts in the Beni Mellal-Kh&eacute;nifra Region, Morocco. Sustainability 2022, 14, 13166. [Google Scholar]
  5. Oumara, N.G.A.; El Youssfi, L. Salinization of Soils and Aquifers in Morocco and the Alternatives of Response. Environ. Sci. Proc. 2022, 16, 65. [Google Scholar]
  6. Seyam, M.M.H.; Haque, M.R.; Rahman, M.M. Identifying the land use land cover (LULC) changes using remote sensing and GIS approach: A case study at Bhaluka in Mymensingh, Bangladesh. Case Stud. Chem. Environ. Eng. 2023, 7, 100293. [Google Scholar] [CrossRef]
  7. Sohl, T.; Sleeter, B. Role of Remote Sensing for Land-Use and Land-Cover Change Modeling. In Remote Sensing of Land Use and Land Cover; CRC Press: Boca Raton, FL, USA, 2012; pp. 225–239. [Google Scholar] [CrossRef]
  8. Beroho, M.; Briak, H.; Cherif, E.K.; Boulahfa, I.; Ouallali, A.; Mrabet, R.; Kebede, F.; Bernardino, A.; Aboumaria, K. Future Scenarios of Land Use/Land Cover (LULC) Based on a CA-Markov Simulation Model: Case of a Mediterranean Watershed in Morocco. Remote Sens. 2023, 15, 1162. [Google Scholar] [CrossRef]
  9. Karmaoui, A.; Ben Salem, A.; El Jaafari, S.; Chaachouay, H.; Moumane, A.; Hajji, L. Exploring the land use and land cover change in the period 2005–2020 in the province of Errachidia, the pre-sahara of Morocco. Front. Earth Sci. 2022, 10, 962097. [Google Scholar] [CrossRef]
  10. Funk, C.; Peterson, P.; Landsfeld, M.; Pedreros, D.; Verdin, J.; Shukla, S.; Husak, G.; Rowland, J.; Harrison, L.; Hoell, A.; et al. The climate hazards infrared precipitation with stations—A new environmental record for monitoring extremes. Sci. Data 2015, 2, 150066. [Google Scholar] [CrossRef] [PubMed]
  11. Copernicus, C.C.S. Agrometeorological Indicators from 1979 up to 2019 Derived from Reanalysis. 2019. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/10.24381/cds.6c68c9bb?tab=doc (accessed on 2 September 2023).
  12. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring vegetation systems in the Great Plains with ERTS. NASA Spec. Publ. 1974, 351, 309. [Google Scholar]
  13. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  14. Zha, Y.; Gao, J.; Ni, S. Use of normalized difference built-up index in automatically mapping urban areas from TM imagery. Int. J. Remote Sens. 2003, 24, 583–594. [Google Scholar] [CrossRef]
  15. Nguyen, C.T.; Chidthaisong, A.; Kieu Diem, P.; Huo, L.Z. A Modified Bare Soil Index to Identify Bare Land Features during Agricultural Fallow-Period in Southeast Asia Using Landsat 8. Land 2021, 10, 231. [Google Scholar] [CrossRef]
  16. Liu, W.T.; Kogan, F.N. Monitoring regional drought using the Vegetation Condition Index. Int. J. Remote Sens. 1996, 17, 2761–2782. [Google Scholar] [CrossRef]
  17. Breiman, L.; Friedman, J.; Stone, C.; Olshen, R. Classification and Regression Trees; Taylor & Francis: Abingdon, UK, 1984. [Google Scholar]
  18. Sheykhmousa, M.; Mahdianpari, M.; Ghanbari, H.; Mohammadimanesh, F.; Ghamisi, P.; Homayouni, S. Support vector machine versus random forest for remote sensing image classification: A meta-analysis and systematic review. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 6308–6325. [Google Scholar] [CrossRef]
  19. Friedman, J.H. Greedy Function Approximation: A Gradient Boosting Machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  20. Hamedianfar, A.; Gibril, M.B.A.; Hosseinpoor, M.; Pellikka, P.K. Synergistic use of particle swarm optimization, artificial neural network, and extreme gradient boosting algorithms for urban LULC mapping from WorldView-3 images. Geocarto Int. 2022, 37, 773–791. [Google Scholar] [CrossRef]
  21. Sahin, E.K. Assessing the predictive capability of ensemble tree methods for landslide susceptibility mapping using XGBoost, gradient boosting machine, and random forest. SN Appl. Sci. 2020, 2, 1308. [Google Scholar] [CrossRef]
Figure 1. Research flow chart.
Figure 1. Research flow chart.
Environsciproc 29 00051 g001
Figure 2. Study area map.
Figure 2. Study area map.
Environsciproc 29 00051 g002
Figure 3. LULC and VCI change maps (2018 and 2022).
Figure 3. LULC and VCI change maps (2018 and 2022).
Environsciproc 29 00051 g003
Figure 4. VCI time series chart (2018 and 2022).
Figure 4. VCI time series chart (2018 and 2022).
Environsciproc 29 00051 g004
Table 2. Validation accuracies of GTB and RF algorithm.
Table 2. Validation accuracies of GTB and RF algorithm.
AlgorithmClassOverall
Accurary
Consumer’s
Accurary
Producer’s
Accurary
Kaffa
Coeffiecient
RFBuilt-Up Area0.920.950.890.89
Water0.90.82
Agriculture0.930.98
Bareland0.890.94
GTBBuilt-Up Area0.930.960.890.91
Water0.960.93
Agriculture0.910.99
Bareland0.860.91
Table 3. The percentage of LULC-changes for each thematic class in 2018 and 2022.
Table 3. The percentage of LULC-changes for each thematic class in 2018 and 2022.
Class2018 (km2)2022 (km2)Changes (km2)
Built-Up Area74.277.5+3.3
Water0.570.47−0.1
Agriculture141.582.8−58.7
Bareland95.1150.7+55.6
Table 4. Climate and VCI changes between 2018 and 2022.
Table 4. Climate and VCI changes between 2018 and 2022.
YearMin. Temp (°C)Max. Temp (°C)Rainfall (mm/yr)Evapotranspiration (mm/yr)Annual Mean VCI
20185.136.9729134339.72
20225.637.5337151619.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mustapha, M.; Zineddine, M. Assessing the Impact of Climate Change on Seasonal Variation in Agricultural Land Use Using Sentinel-2 and Machine Learning. Environ. Sci. Proc. 2024, 29, 51. https://doi.org/10.3390/ECRS2023-16365

AMA Style

Mustapha M, Zineddine M. Assessing the Impact of Climate Change on Seasonal Variation in Agricultural Land Use Using Sentinel-2 and Machine Learning. Environmental Sciences Proceedings. 2024; 29(1):51. https://doi.org/10.3390/ECRS2023-16365

Chicago/Turabian Style

Mustapha, Musa, and Mhamed Zineddine. 2024. "Assessing the Impact of Climate Change on Seasonal Variation in Agricultural Land Use Using Sentinel-2 and Machine Learning" Environmental Sciences Proceedings 29, no. 1: 51. https://doi.org/10.3390/ECRS2023-16365

APA Style

Mustapha, M., & Zineddine, M. (2024). Assessing the Impact of Climate Change on Seasonal Variation in Agricultural Land Use Using Sentinel-2 and Machine Learning. Environmental Sciences Proceedings, 29(1), 51. https://doi.org/10.3390/ECRS2023-16365

Article Metrics

Back to TopTop