<figure><img src="image_1.png" /></figure>

HAL
open science

# Reclassification Scheme for Image Analysis in GRASS GIS Using Gradient Boosting Algorithm: A Case of Djibouti, East Africa

Polina Lemenkova

► To cite this version:

Polina Lemenkova. Reclassification Scheme for Image Analysis in GRASS GIS Using Gradient Boosting Algorithm: A Case of Djibouti, East Africa. Journal of Imaging, 2025, 11 (8), pp.249. 10.3390/jimaging11080249. [https://hal.science/hal-05178165v1](https://hal.science/hal-05178165v1)

HAL Id: hal-05178165

Submitted on 30 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Reclassification Scheme for Image Analysis in GRASS GIS Using Gradient Boosting Algorithm: A Case of Djibouti, East Africa

Polina Lemenkova ID

Department of Biological, Geological and Environmental Sciences, Alma Mater Studiorum-Università di Bologna,
Via Irnerio 42, Emilia-Romagna, 40126 Bologna, Italy; polina.lemenkova2@unibo.it; Tel.: +39-3446928732

## Abstract

Image analysis is a valuable approach in a wide array of environmental applications. Mapping land cover categories depicted from satellite images enables the monitoring of landscape dynamics. Such a technique plays a key role for land management and predictive ecosystem modelling. Satellite-based mapping of environmental dynamics enables us to define factors that trigger these processes and are crucial for our understanding of Earth system processes. In this study, a reclassification scheme of image analysis was developed for mapping the adjusted categorisation of land cover types using multispectral remote sensing datasets and Geographic Resources Analysis Support System (GRASS) Geographic Information System (GIS) software. The data included four Landsat 8-9 satellite images on 2015, 2019, 2021 and 2023. The sequence of time series was used to determine land cover dynamics. The classification scheme consisting of 17 initial land cover classes was employed by logical workflow to extract 10 key land cover types of the coastal areas of Bab-el-Mandeb Strait, southern Red Sea. Special attention is placed to identify changes in the land categories regarding the thermal saline lake, Lake Assal, with fluctuating salinity and water levels. The methodology included the use of machine learning (ML) image analysis GRASS GIS modules 'r.reclass' for the reclassification of a raster map based on category values. Other modules included 'r.random', 'r.learn train' and 'r.learn predict' for gradient boosting ML classifier and 'i.cluster' and 'i.maxlik' for clustering and maximum-likelihood discriminant analysis. To reveal changes in the land cover categories around the Lake of Assal, this study uses ML and reclassification methods for image analysis. Auxiliary modules included 'i.group', 'r.import' and other GRASS GIS scripting techniques applied to Landsat image processing and for the identification of land cover variables. The results of image processing demonstrated annual fluctuations in the landscapes around the saline lake and changes in semi-arid and desert land cover types over Djibouti. The increase in the extent of semi-desert areas and the decrease in natural vegetation proved the processes of desertification of the arid environment in Djibouti caused by climate effects. The developed land cover maps provided information for assessing spatial-temporal changes in Djibouti. The proposed ML-based methodology using GRASS GIS can be employed for integrating techniques of image analysis for land management in other arid regions of Africa.

**Keywords:** image processing; image analysis; remote sensing; satellite image; geoinformatics; machine learning

<figure><img src="image_3.png" /></figure>

### check for updates

Academic Editors: Chuang Niu, Qian Wang, Xin Cao and Shenghan Ren

Received: 14 June 2025

Revised: 16 July 2025

Accepted: 21 July 2025

Published: 23 July 2025

**Citation:** Lemenkova, P. Reclassification Scheme for Image Analysis in GRASS GIS Using Gradient Boosting Algorithm: A Case of Djibouti, East Africa. *J. Imaging* **2025**, *11*, 249. https://doi.org/10.3390/jimaging11080249

**Copyright:** © 2025 by the authors. License MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
# 1. Introduction

## 1.1. Background

Satellite images play a pivotal role in environmental monitoring for reconstructing the history of landscapes on Earth. Understanding how land cover types change over time gives insight into the development of human-environment interactions [1–3]. The main threat of environmental changes are linked to climate issues such as desertification, erosion, landscape fragmentation, coastal degradation and deforestation. Most of this information can be discovered on spaceborne images due to the fundamental principle of remote sensing (RS)—differences in spectral reflectance of pixels corresponding to diverse objects on the Earth. Using such benefits of RS, many environmental processes and phenomena can be recognised on satellite images and Earth Observation (EO) data, including, for instance, deforestation, soil erosion and desertification.

The important application of EO data for environmental monitoring consists of time series analysis. Such a technique is based on the comparison of changes visible in several multi-temporal images to evaluate trends in vegetation coverage during a given temporal period [4,5]. This enables the detection of land cover changes and the identification of landscape dynamics [6,7]. EO data processing requires advanced techniques that allow us to find optimal and time-efficient solutions for solving complex geospatial tasks in image analysis. While traditional methods of Geographical Information Systems (GISs) have limitations, programming algorithms have significant benefits in the time series analysis of satellite images.

Technologies of image analysis play an essential role in the development of environmental monitoring. Traditional tools for processing RS data rely on GISs [8,9]. However, due to the large volume of RS data, the automation of image analysis is needed to effectively perform time series analysis and extract geoinformation. Recent technological developments have enabled novel methods for extracting geoinformation from spaceborne data, including image enhancement [10], image differentiation, rationing or principal component analysis [11]. Machine learning (ML) methods have been found to improve the efficiency of image analysis through automated algorithms. Existing cases of ML algorithms in cartography and RS tend to be related to problems of image processing that are similar to those faced by decision algorithms. In each case, ML algorithms aim to identify the best solution strategy for image classification [12].

Examples of ML algorithms include Random Forest (RF) [13,14], Support Vector Machines (SVMs) [15] and Convolutional Neural Networks (CNNs) [16]. Such approaches are derived from programming and are adapted to image analysis [17,18]. Examples of novel techniques in image analysis include heuristic clustering [19–22], flexible control of variables [23], optimised pattern recognition [24,25], ensemble learning [26] and objective analysis [27,28].

Following such applications of ML in Earth and climate sciences, this study reports the case of a novel independent reclassification scheme for image analysis at the pixel level using ML coupled with Geographic Resources Analysis Support System (GRASS) GIS and gradient boosting algorithm. Its advantages include the optimisation approach, which considers previous iterative cycles of image classification for the generation of new models and the narrowing of pixel assignment options to various land cover classes.

## 1.2. Objectives and Goals

The main goal is to identify changes in land cover types of Djibouti using the time series of satellite images by the reclassification scheme. GRASS GIS version 8.4.1 is an advanced scripting cartographic and RS software with a powerful computational engine for geospatial processing [29]. Developed by the international GRASS developer community.
nity [30], it is used for processing cartographic data, maps, models and satellite image analysis across various applications in geoinformatics. Some of its important benefits are a built-in framework and a Python Application Programming Interface (API), which supports the automation of data processing through scripting. To achieve this goal, the objective is to map land cover changes in the coastal region of Djibouti using RS data. The scripts are included and commented for reference. The specific objectives of this research are as follows:

1. To produce land cover maps based on the classification of Landsat satellite images from 2015, 2019, 2021 and 2023;

2. To quantify and visualise dynamic changes in land cover types for the identification of landscape dynamics;

3. To assess the accuracy of classification for vegetation mapping in central and coastal zones around Lake Assal.

The rest of this manuscript is organised as follows. Section 1 presents the introduction with a related literature review and identifies the objectives and goals. Section 2 describes the study area, focusing on environmental and climate–hydrological settings that influence landscape types. Section 3 introduces the materials and numerical methods used in this study, including the dataset and workflow, which are described as the methodological approach. Major technical procedures are described in the subsections that follow. In Section 4, we present and discuss the results with cases of land cover changes, comment on uncertainties, and provide and interpret a series of maps. Finally, Section 5 concludes the manuscript, summarising the findings and providing recommendations for future studies.

## 2. Study Area

This study area is located in Djibouti and surrounding coastal areas in the Gulf of Tadjoura, Bab-el-Mandeb Strait (Figure 1). Situated in the region of the Horn of Africa and surrounded by waters of the southern Red Sea, Djibouti is a strategic point with geo-maritime importance [31,32]. Djibouti is one of the hottest and driest places on Earth with many issues related to climate and weather processes [33]. The thermal crater saline lake, Lake Assal, is located in this area, which is the lowest point on land in East Africa with extreme levels of evaporation. The topography affects temperature and precipitation, which vary according to the elevation and differ in coastal and mountainous regions [34].

This region experiences fluctuations in salinity and water levels due to seasonal and annual changes in temperature and precipitation. Dominating landscapes are represented by dry mountain vegetation typical of arid climates. The distribution and genera of flora and fauna are controlled by low precipitation and high temperatures in semi-desert and desert environments [35–37]. For example, the 2011 drought in East Africa affected rain-fed agriculture in Djibouti, which led to a regional famine [38–40]. Drought-induced crises endanger long-term investments and development, with catastrophic human costs and effects on vulnerable regions of the country [41–43].

The drought-prone regions of Djibouti have experienced strong effects from climate change, including fluctuations in the hydrological regime of thermal lakes. The arid climate has triggered environmental problems; for example, it has desiccated and abandoned agricultural lands, and salinated and blocked water sources due to sea level rise and flash floods [44,45]. Soil erosion and land degradation vary across topographically diverse landscapes, reflecting the prevailing climatic conditions. [46,47]. Moreover, environmental pollution has been reported along the coasts of the Gulf of Tadjoura [48,49].
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Topographic map of Djibouti with indicated study area (yellow rotated square) which shows the location of the Landsat images. Mapping software: Generic Mapping Tools (GMT) version 6.4.0. Data source: General Bathymetric Chart of the Oceans (GEBCO). Map source: Author.</figcaption></figure>

The surroundings of the Horn of Africa are notable for climate-hydrological fluctuations which cause land cover changes [50,51]. The coastal region located in the Bab-el-Mandeb region includes diverse types of xeric grasslands and shrubland typical for arid ecoregions (Figure 2). This region consists of various land cover classes typical for semi-desert areas: dry shrub steppe and grassland, including species of *Vachellia tortilis*, *Vachellia nubica*, and *Balanites aegyptiaca*. Many scarce vegetation spots include scattered trees or shrubs on the sandy plains of bare lands and semi-deserts.

The eastern strip along the Red Sea coast is part of the Eritrean coastal desert [52]. It includes vegetation along the muddy wadi mouths in diverse habitats [53], grasslands [54,55], and secondary growth areas in the urban areas. The floodplains along the wadi rivers include scarce vegetation typical for semi-desert ecosystems. Xeric vegetation, farmland, mangroves along the coasts of the Red Sea and valleys of wadi [56,57], coral reefs along the coasts [58–60] and urban areas collectively make up the region's landscape. Djibouti is one of the most seismically unstable zones of Africa with high tectonic-magmatic activities [61–64] due to the closeness of the Afar Triple Junction [65–67]. The thermal regime around the Red Sea is explained by its structural connection to the active zone of the East African Rift [68–72].

Such complex geological setting and arid environment have resulted in the deficit of water and favoured the distribution of thermal springs with high salinity, located in three regions—Lake Assal, Lake Hanle and Lake Abhe [73,74]. The craton origin of Lake Assal
is related to the Afar Triple Junction [75-77]. The geochemical composition of the Assal
Rift includes evaporitic halite deposits [78,79]. Its location in the arid zone with strong
evaporation and no riverine input contributed to the lake's increase in salinity, with a high
brine concentration above 100 g/L [80]. Fluctuations of this concentration are related to
evaporation, the hot climate and geothermal processes [81].

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Map of land cover types over Djibouti, showing the distribution of major vegetation and land cover types. Software: Quantum Geographic Information System (QGIS). Data source: Food and Agriculture Organisation (FAO). Map source: Author.</figcaption></figure>

**3. Materials and Methods**

Our methodology presents a case of image reclassification and analysis using ML
techniques for environmental mapping in Djibouti, East Africa. A classification scheme con-
sisting of 17 initial land cover classes was employed through a logical sequence of GRASS
GIS modules to extract land cover information of the study area in Djibouti. Reclassification
includes the sequence of changing the values in a raster dataset to new values. The process
is based on a set of rules which defined the criteria of updates in classes. This operation in
GRASS GIS was applied for the simplification of land cover polygons, standardisation of
assignment of pixels to classes and preparation of data for land cover analysis.

**3.1. Data**

The data used in this study include four images from the Landsat 8–9 Operational
Land Imager (OLI) and Thermal Infrared Sensor (TIRS): the training image from 2015, and
images with biannual intervals from 2019, 2021 and 2023 (Figure 3).

The Landsat images were selected as the data source due to their quality and reliability[82–84]. Earlier studies reported the effective use of the RS data in environmental studies which prove their utility and value [85–91]. Nowadays, satellite images are excellent sources of geoinformation for mapping arid regions [12,92,93]. Spaceborne data are successfully used in the analysis of land cover change and deforestation in tropical forests, and for reporting landscape dynamics [94–97]. The importance of Landsat data for environmental research has been described earlier [98–100].

In this study, four Landsat images were taken on the following dates: (1) 2 July 2015;
(2) 23 July 2017; (3) 19 August 2021; and (4) 16 July 2023. The metadata are summarised
in Table 1. The cloudiness percentage is less than 10% for all the images. The technical
parameters of the scenes are provided in Table 2. Common data for all four images are
as follows: Datum and Ellipsoid: World Geodetic System 1984 (WGS84); Product Map Projection L1: Universal Transverse Mercator (UTM) zone 38 (for Djibouti); Worldwide Reference System (WRS): Path—166, Row—52; Landsat Collection Number: 2; Category: T1; Station Identifier: LGN; and Sensor Identifier: OLI/TIRS. The images were taken during day time in nadir with a roll angle of 1° for the 2015 image and 0° for the images from 2017, 2021 and 2023.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Original Landsat 8–9 OLI/TIRS scenes of Djibouti from 2015, 2019, 2021 and 2023 used for image processing. Image source: United States Geological Survey (USGS). Compilation: Author.</figcaption></figure>

<table><caption>**Table 1.** Identifiers (IDs) and acquisition dates of the multispectral Landsat 8–9 OLI/TIRS images.</caption><thead><tr><td>Date</td><td>Landsat Product Identifier</td><td>Scene ID</td></tr></thead><tbody><tr><td>2 July 2015</td><td>LC08_L2SP_166052_20150702_20200909_02_T1</td><td>LC81660522015183LGN01</td></tr><tr><td>23 July 2017</td><td>LC08_L2SP_166052_20170723_20200903_02_T1</td><td>LC81660522017204LGN00</td></tr><tr><td>19 August 2021</td><td>LC08_L2SP_166052_20210819_20210827_02_T1</td><td>LC81660522021231LGN00</td></tr><tr><td>16 July 2023</td><td>LC09_L2SP_166052_20230716_20230719_02_T1</td><td>LC91660522023197LGN00</td></tr></tbody></table>

The image quality score was 9 for all scenes, according to the Landsat specifications. The remaining various metadata for each scene are summarised in Table 2. The geographic characteristics of the images are referenced in Table 3. The Landsat Scene Identifiers (IDs) are provided below.

The images have been filtered to reduce noise due to cloudiness and atmospheric effects. Images were systematically selected from the USGS EarthExplorer repository and
taken at the end of the dry season. During this period, cloud cover is generally low but vegetation appears bright green and contrasts with all other land cover types.

<table><caption>**Table 2.** Metadata of the Landsat 8–9 OLI/TIRS satellite images. Source: USGS (EarthExplorer).</caption><thead><tr><th>Dataset Attribute</th><th>Attribute Value</th><th>Attribute Value</th><th>Attribute Value</th><th>Attribute Value</th></tr></thead><tbody><tr><td>Date Acquired</td><td>16 July 2023</td><td>19 August 2021</td><td>23 July 2017</td><td>2 July 2015</td></tr><tr><td>Date Product Generated L2</td><td>19 July 2023</td><td>27 August 2021</td><td>3 September 2020</td><td>9 September 2020</td></tr><tr><td>Date Product Generated L1</td><td>16 July 2023</td><td>27 August 2021</td><td>3 September 2020</td><td>9 September 2020</td></tr><tr><td>Land Cloud Cover</td><td>0.02</td><td>0.19</td><td>0.45</td><td>0.08</td></tr><tr><td>Scene Cloud Cover L1</td><td>0.06</td><td>0.71</td><td>12.03</td><td>0.08</td></tr><tr><td>Ground Control Points Model</td><td>752</td><td>722</td><td>742</td><td>768</td></tr><tr><td>Ground Control Points Version</td><td>5</td><td>5</td><td>5</td><td>5</td></tr><tr><td>Geometric RMSE Model</td><td>3563</td><td>2963</td><td>2977</td><td>2098</td></tr><tr><td>Geometric RMSE Model X</td><td>2240</td><td>2009</td><td>1967</td><td>1089</td></tr><tr><td>Geometric RMSE Model Y</td><td>2771</td><td>2179</td><td>2235</td><td>1794</td></tr><tr><td>Processing Software Version</td><td>LPGS_16.3.0</td><td>LPGS_15.5.0</td><td>LPGS_15.3.1c</td><td>LPGS_15.3.1c</td></tr><tr><td>Sun Elevation L0RA</td><td>62.48488638</td><td>64.48715625</td><td>62.85521856</td><td>62.40770766</td></tr><tr><td>Sun Azimuth L0RA</td><td>65.62042782</td><td>84.79256539</td><td>68.51197578</td><td>61.89003236</td></tr><tr><td>TIRS SSM Model</td><td>N/A</td><td>FINAL</td><td>FINAL</td><td>ACTUAL</td></tr><tr><td>Scene Center Latitude</td><td>11.56786</td><td>11.56779</td><td>11.56783</td><td>11.56778</td></tr><tr><td>Scene Center Longitude</td><td>42.94139</td><td>42.97394</td><td>42.97879</td><td>42.96942</td></tr><tr><td>Corner Upper Left Latitude</td><td>12.60889</td><td>12.60924</td><td>12.60931</td><td>12.60921</td></tr><tr><td>Corner Upper Left Longitude</td><td>41.88103</td><td>41.91137</td><td>41.91688</td><td>41.90861</td></tr><tr><td>Corner Upper Right Latitude</td><td>12.62518</td><td>12.62530</td><td>12.62532</td><td>12.62529</td></tr><tr><td>Corner Upper Right Longitude</td><td>43.98441</td><td>44.01755</td><td>44.02308</td><td>44.01479</td></tr><tr><td>Corner Lower Left Latitude</td><td>10.50122</td><td>10.50151</td><td>10.50157</td><td>10.50149</td></tr><tr><td>Corner Lower Left Longitude</td><td>41.90427</td><td>41.93439</td><td>41.93986</td><td>41.93165</td></tr><tr><td>Corner Lower Right Latitude</td><td>10.51472</td><td>10.51483</td><td>10.51484</td><td>10.51482</td></tr><tr><td>Corner Lower Right Longitude</td><td>43.99199</td><td>44.02489</td><td>44.03037</td><td>44.02215</td></tr></tbody></table>

<table><caption>**Table 3.** Technical characteristics and geospatial location of the multispectral Landsat images.</caption><thead><tr><th>N-S</th><th>Data Values</th><th>E-W</th><th>Data Values</th></tr></thead><tbody><tr><td>North:</td><td>1,395,915.00</td><td>East:</td><td>393,015.00</td></tr><tr><td>South:</td><td>1,162,485.00</td><td>West:</td><td>164,085.00</td></tr><tr><td>Resolution:</td><td>30.00 m</td><td>Cells:</td><td>59,376,811</td></tr><tr><td>Rows:</td><td>7781</td><td>Columns:</td><td>7631</td></tr></tbody></table>

## 3.2. Workflow

The schematic workflow of the research is presented in Figure 4.

The image processing used to create a land cover map from Landsat OLI/TIRS imagery was carried out in six steps using the following techniques of GRASS GIS:

* Image import;

* Data exploration and analysis;

* Classification of images using MaxLike method;

* Reclassification scheme using 'echo' of Unix and 'r.reclass' module;

* Classification of images using gradient boosting algorithm of ML;

* Image post-processing and mapping;

* Accuracy assessment.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Workflow scheme. Source: Author.</figcaption></figure>

The workflow is based on scripting techniques of GRASS GIS for image analysis (software developed originally by the U.S. Army Construction Engineering Research Laboratories in 1982, now maintained by the developer team in Champaign, IL, USA). The topographic map is plotted using existing techniques in GMT version 6.4.0. [101–105] (software developed originally in 1988 by Paul Wessel and Walter H. F. Smith, in Lamont-Doherty Earth Observatory of Columbia University, Palisades, NY, USA).

### 3.3. Image Preprocessing

The Landsat 8–9 OLI/TIRS images were geometrically rectified according to EPSG coordinates: 2095, Universal Transverse Mercator (UTM) zone for Djibouti according to the USGS technical characteristics. In this study, a scripting method of GRASS GIS was applied, which uses several modules with code snippets explained and commented below.

The images were imported into the working folder using the 'r.import' module, then the same procedure was repeated for all the bands and the four images used in this study (2015, 2019, 2021 and 2023). Afterwards, the relief map was imported from GEBCO. The region was set up to the geospatial extent with the isolines obtained using the 'r.contour' module with a 200 m interval. To facilitate data processing, the best combinations of Landsat bands were stacked into coloured composites. The bands that produced the best separation for pixels were processed, excluding the panchromatic channels. The Tag Image File Format (TIFF) files were joined using the 'i.group' module and enhanced using preprocessing steps which included contrast stretching and atmospheric correction.
### 3.4. Clustering

The clustering was performed by the 'i.cluster' module of GRASS GIS, which generates a signature file and reports results using the k-means clustering algorithm, as explained in previous studies [106–109]. Signature file determines the categories (classes) of land cover types using the analysis of spectral reflectance of the pixels and contains clusters and covariance matrices for each image. The initial means for each class before data processing reported in Appendix A indicate the starting seed values of the multispectral band before iteration is assigned according to the signature files. The means are then calculated during the 'i.cluster' procedure and recalculated following the principles of maximal class separation (the contrast between the classes) and minimal class size (the details of classification with respect to a spatial resolution of 30 m). The table in Appendix A reports such values for each of the 10 classes and 7 multispectral bands of Landsat, starting from 1—'Coastal aerosol' to 7—'SWIR-2'.

The initial mean values for each land cover class for the years 2015, 2017, 2021 and 2023 are reported in Tables A1–A4; see Appendix A. Afterwards, the assignment of pixels was performed iteratively for each land cover class. The same was applied for all images for the years 2015, 2017, 2021 and 2023 with a biannual gap. The classification was performed using the maximum-likelihood discriminant analysis classifier by the 'i.maxlik' module [110,111]. The algorithm used the signature file calculated in the previous step. Using the signature file, the land cover types were classified using maximum-likelihood discriminant analysis (Figure 5).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Maps of land cover types based on the classified Landsat 8–9 OLI/TIRS images from 2015, 2019, 2021 and 2023 covering Djibouti. The approach involved the maximum-likelihood discriminant analysis classifier. Background relief: GEBCO. Software: GRASS GIS. Mapping source: Author.</figcaption></figure>
### 3.5. Accuracy Analysis

The analysis of accuracy (Figure 6) revealed the areas of misclassified pixels within water bodies, including the Bab-el-Mandeb Strait. The accuracy assessment, showing the probability of correctly assigning each pixel to its respective class across all images, is presented in Figure 6 and calculated with the averages of each calculated band. Afterwards, based on the classified ten categories in the images, two of them ('water in the sea' and 'water in shelf' areas) were merged as a single class 'water', and the names of the classes were assigned to each numerical cluster.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Accuracy analysis of image classification using reject threshold probability algorithm and the chi-squared test, which was applied to the reclassified images. Source: Author.</figcaption></figure>

After this step, ML was applied to all the images. The ML step started with the generation of training pixels from land cover classification performed in 2015, which were used as seed data. The extraction of training data was conducted using the 'r.random' module. Afterwards, the training pixels were used for classification on the Landsat images (example for 2023) by training a model using the 'r.learn train' module. The gradient boosting classifier was used as an algorithm to perform ML classification. Next, the prediction was performed using the 'r.learn predict' ML module, which applied a fitted estimator from the Python's Scikit-Learn library [112] to raster images in the imagery group. In the next step, the computed raster categories were automatically applied to the classification output and checked using the 'r.category' module.

### 3.6. Gradient Boosting

The images were processed using the gradient boosting ML method of image classification. The essential concept of this algorithm is setting the target outcomes of classified images for the next model based on previous iterations, which in the end minimises the
error of classification. Gradient boosting presents an advanced ensemble technique of ML, which combines the predictions of multiple classification steps by the principle of decision trees. Such an approach presents a novel method of satellite image processing based on the supervisor approach. The algorithm classifies images using the iterative approach which considers the previous and next classification steps as loops. It combines several decision trees as simplified learners to create a strong prediction mode. By minimising a loss function embedded in the algorithm, each new learner fixes the mistakes committed by the group of prior learners. In this way, the accuracy and predictive capacity of the classifier are increased iteratively, which is directed by the gradient of the loss function.

The core idea is that the algorithm finds the best possible next model when combined with previous ones. In this way, it minimises the overall prediction error. The advantage of such an approach is that it improves the accuracy of classification through the convergence of decision cycles and assigns cells that constitute the mosaic of pixels on a raster image to the defined categories. When several images are processed as a time series, this enables us to effectively monitor landscape dynamics and detect changes.

The class separability matrices were computed for each pair of years (2015–2017 and 2021–2023) and show the distinguishability between classes (Equations (A1) and (A2)). The iterative cycles of image classification are summarised in Appendix B in four statistical tables for each of the study periods, respectively. They report the pixels’ distribution by categories of land cover classes for each year: 2015, 2017, 2021 and 2023. The values indicate percent convergence, which shows at which values the cluster means of land cover categories become stable during iterative classification. Hence, iterative processes aim to achieve the maximal percentage of pixels that no longer move from cluster to cluster during iteration and reach the best possible values.

### 3.7. Reclassification

To improve accuracy, the images were reclassified using the principle of reclassification scheme. When clusters are being generated by the algorithm, the means of classified categories constantly change, because cells are assigned to these classes and the mean is then recalculated to include new pixels. The goal of the iteration is to maximise the distances between the classes through recalculation. To this end, after the creation of all clusters and the assignment of pixels into these groups, the algorithm ‘i.cluster’ changes cluster means through iteration processes. In this way, the algorithm attempts to increase the contrast between the classes, i.e., the numerical gap in values between the classes.

### 3.8. Generalisation of FAO Land Cover Map

The FAO land cover map of Djibouti was generalised from the classes of the Land Cover Classification System (LCCS) according to a technical report on the country. The scheme covering land cover categories comprising 10 classes was used and adopted for the study areas within the country (Table 4). This part of the work was performed using QGIS software, version 3.42.2 ‘Münster’. The QGIS was initially developed by G. Sherman and now a project of the Open Source Geospatial Foundation.

The reclassified image was generated for each year using the ‘r.reclass’ module of GRASS GIS and ‘echo’ Unix command. The raster map was reclassified based on the category values, and new raster maps for the years 2015, 2019, 2021 and 2023 were generated. The category values of the reclassified maps are based on an iteration of the categories using the information generated in ‘landusereclass.txt’ file. These categories were controlled using the ‘r.category’ module. The maps were visualised using modules ‘d.mon’ and ‘g.region’. The color palette was selected from the color tables. Afterwards, the maps were
displayed using a combination of 'd.rast' and 'd.vect' modules. Finally, a cartographic grid was added along with the legends by the 'd.legend' module.

<table><caption>**Table 4.** Training pixels for validation under LCCS for Djibouti, derived from FAO classification scheme of land cover types: spatial distribution of the training ground truth samples.</caption><thead><tr><th>ID</th><th>Landsat 8-9 OLI/TIRS Land Cover Class</th><th>FAO LCCS Codes</th></tr></thead><tbody><tr><td>1</td><td>Salt Plains and Hardpans</td><td>6020</td></tr><tr><td>2</td><td>Barren Land</td><td>6001, 6004</td></tr><tr><td>3</td><td>Water Bodies</td><td>7001 // 8001</td></tr><tr><td>4</td><td>Mangroves and Aquatic Vegetation on Regularly Flooded (Temporarily or Permanently) Fresh or Brackish Water</td><td>41,653-R1,R2 // 41,739-R1<br>41,521-R1,R2 // 41,597-R1,R2</td></tr><tr><td>5</td><td>Bushes and Shrubs</td><td>11,490 // 11,494</td></tr><tr><td>6</td><td>Farmland, Shrubs and Irrigated Trees</td><td>11,499 // 11,500 // 30,001,<br>11,491, 11,495</td></tr><tr><td>7</td><td>Built-up Areas, Artificial Surfaces and Associated Areas</td><td>0010</td></tr><tr><td>8</td><td>Broadleaved Semi-deciduous Forest // Woodland</td><td>21,496 // 21,497-15,048 //<br>21,496-12,134 // 21,497-12,940</td></tr><tr><td>9</td><td>Mosaic Cropland (Cultivated Areas)</td><td>0003-0004</td></tr><tr><td>10</td><td>Sparse Vegetation</td><td>20,049, 20,056 20,058</td></tr></tbody></table>

Information source: FAO (land use statistics).

The reclassification was performed since water surfaces in the Bab-el-Mandeb Strait and southern Red Sea had different colours. Nevertheless, they belong to the same class of water. The differences in the coastal regions were smooth at the time of image acquisition and appear as dark colours in all images. With each iteration, the values of the means shift to a higher percentage of pixels within each cluster. This is illustrated in the class separability matrix, which shows the distinction between the categories (Equations (A1) and (A2); Appendix C). As means never become completely static, a % convergence and a maximum number of iterations are defined to finalise the process of reclassification before the maximum number of cycles is reached. The final values of the computed class means for each land category are reported in Appendix D. Once the maximum number of reclassification cycles is reached, the optimal % convergence is achieved through the increased number of iterations during the 'i.cluster' procedure. Here, the number of cycles depends on the complexity of the terrain and land cover patterns.

# 4. Results and Discussion

## 4.1. Land Cover Change Analysis

The results of the classified satellite images are presented in Figure 5. The study area was categorised into 10 distinct land cover classes following the FAO classification scheme: (1) salt plains and hardpans; (2) barren land; (3) water bodies; (4) mangroves and aquatic vegetation on regularly flooded (temporarily or permanently) fresh or brackish water; (5) bushes and shrubs; (6) farmland, shrubs and irrigated trees; (7) built-up areas, artificial surfaces and associated areas; (8) broadleaved semi-deciduous forest and woodland; (9) mosaic cropland in cultivated areas; and (10) sparse vegetation in desert areas. Among these classes, the following categories included sub-clusters: water bodies comprise ponds, areas covered by marshes, rivers, estuaries, and coastal waterways; bushland and shrubland include non-classified types of vegetation such as grasslands, shrubs, rangelands, and savannah; and coastal areas surrounding Lake Assal are categorised by salt-tolerant vegetation and mangroves.
Over an estimated period of a decade, from 2013 to 2023, image analysis revealed certain changes in land cover types across the study area of Djibouti. Thus, it is noteworthy that the amount of mangroves along the coastal areas decreased slightly, shrinking by 0.35 km² from 6.28 km² to 5.93 km² during this period. Furthermore, a large area of 1392 km² that had previously been covered by bushes and shrubland was converted into other land use categories, resulting in a decrease in coverage from 1194 km². As a result of this reduction, the percentage of bush cover decreased to 15%. Over the course of the studied period, both salt-covered areas around Lake Assal and croplands showed downward tendencies, with farmland areas falling from 3.21 km² to 2.75 km² and salt coverage from 17 km² to 16.04 km².

Additionally, saline waters mixed with the fresh waters of the estuaries of the Gulf of Tadjoura have a lighter coloration than those of the open sea. Vegetation types such as mosaic croplands, sparse vegetation, xeric shrubland and grassland or bare soil areas are detected and indicated on the maps. Broadleaved deciduous vegetation is relatively scarce in the study area. The regions of shrub, and flooded brackish water areas are mostly occupied by mangrove forests, which have stronger spectral signals and appear as various shades of colors. The computational results of the pixels by land cover classes are summarised in Table 5.

<table><caption>**Table 5.** Distribution of pixels by land cover classes: biannual variations in study area of Djibouti over 10 land cover classes in summer periods from 2015 to 2023.</caption><thead><tr><th>Class</th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>2015</td><td>1531</td><td>251</td><td>314</td><td>755</td><td>964</td><td>761</td><td>699</td><td>616</td><td>448</td><td>591</td></tr><tr><td>2017</td><td>1452</td><td>361</td><td>120</td><td>746</td><td>1031</td><td>813</td><td>672</td><td>714</td><td>502</td><td>609</td></tr><tr><td>2021</td><td>588</td><td>1225</td><td>223</td><td>699</td><td>1007</td><td>798</td><td>608</td><td>644</td><td>495</td><td>644</td></tr><tr><td>2023</td><td>842</td><td>915</td><td>367</td><td>733</td><td>965</td><td>840</td><td>672</td><td>627</td><td>533</td><td>547</td></tr></tbody></table>

Xeric shrubland and sparse deciduous forests typical in Djibouti are located along the coasts of the Red Sea or Gulf of Tadjoura. Spectral signatures of this class differ from those of mosaic vegetation (semi-deciduous forests or sparse lands) found in inland regions of the country, on the border with Ethiopia and Eritrea. Thin coastal regions covered with mangrove plants along the coasts of the Bab-el-Mandeb Strait are identified as wetlands partially submerged in water. This type of vegetation shows little difference compared to the salt areas of Lake Assal, which experience significant changes.

There were also notable expansions in other land cover categories of landscapes in Djibouti. For example, the area of bare lands and in the interior deserts (western region of the country) increased significantly between the estimated period of 2015 and 2023. Specifically, the areas expanded from 5.56 km² to 6.24 km². This change highlights the processes of desertification of regional landscapes, which is caused by recent climate change and global warming. Hence, desertification processes have become notable over time, with arid regions occupying more areas in inner regions of the country.

## 4.2. Uncertainties and Sources of Error

Understanding possible uncertainties underlying the classification schemes in RS data processing is imperative. Accurate vegetation mapping using satellite image classification is not easy to achieve due to the spectral confusion between different vegetation types and similarity in spectral reflectance of various plant species [113–116]. Therefore, accuracy assessment is a crucial step in the evaluation process. In this case, the generated map was compared to the ground truth values, represented in this case by sample points taken from the FAO-based map. The overall accuracy (OA) and Kappa coefficient (KC) were computed
to evaluate the classified map's accuracy by measuring the similarity between the sample points and the correctly and incorrectly categorised classes (Table 6).

<table><caption>**Table 6.** Error matrix for validation of gradient boosting ML algorithm across land cover categories.</caption><thead><tr><th>Class</th><th>Salt Plains</th><th>Barren Land</th><th>Water</th><th>Mangrove</th><th>Shrub</th><th>Farmland</th><th>Built-Up Area</th><th>Forest</th><th>Cropland</th><th>Sparse Vegetation</th><th>Total</th></tr></thead><tbody><tr><td>1</td><td>3503</td><td>12</td><td>4</td><td>7</td><td>31</td><td>4</td><td>29</td><td>2</td><td>9</td><td>21</td><td>3622</td></tr><tr><td>2</td><td>1</td><td>1651</td><td>43</td><td>4</td><td>5</td><td>0</td><td>8</td><td>12</td><td>0</td><td>2</td><td>1726</td></tr><tr><td>3</td><td>0</td><td>1</td><td>1018</td><td>21</td><td>0</td><td>8</td><td>15</td><td>4</td><td>2</td><td>12</td><td>1081</td></tr><tr><td>4</td><td>2</td><td>3</td><td>14</td><td>945</td><td>1</td><td>0</td><td>38</td><td>51</td><td>27</td><td>12</td><td>1093</td></tr><tr><td>5</td><td>3</td><td>18</td><td>256</td><td>12</td><td>2751</td><td>12</td><td>48</td><td>35</td><td>0</td><td>0</td><td>3135</td></tr><tr><td>6</td><td>22</td><td>16</td><td>39</td><td>0</td><td>2</td><td>546</td><td>5</td><td>8</td><td>31</td><td>83</td><td>752</td></tr><tr><td>7</td><td>2</td><td>6</td><td>10</td><td>1</td><td>12</td><td>14</td><td>1132</td><td>44</td><td>6</td><td>13</td><td>1240</td></tr><tr><td>8</td><td>1</td><td>1</td><td>0</td><td>204</td><td>37</td><td>52</td><td>0</td><td>1439</td><td>2</td><td>0</td><td>1736</td></tr><tr><td>9</td><td>3</td><td>12</td><td>2</td><td>62</td><td>14</td><td>6</td><td>0</td><td>45</td><td>1968</td><td>0</td><td>2033</td></tr><tr><td>10</td><td>14</td><td>2</td><td>0</td><td>16</td><td>83</td><td>12</td><td>0</td><td>26</td><td>0</td><td>2382</td><td>2535</td></tr></tbody></table>

$$OA = 93.4\%; KC = 87.5\%$$

The sources of errors in automatic classification are also caused by the individual patterns of land cover categories on Earth; the data are not normally distributed in most cases. Moreover, the light signature of the deciduous broadleaved plants enabled us to detect occasionally growing plants on the mountainous slopes in the central regions of the country where precipitation is higher than in semi-desert areas. Thus, deciduous shrubland in all the images is attributed to the high backscatter coefficients of this vegetation type and spectral reflection effects. Therefore, validation using the computed KC was performed and is reported in Table 6.

In order to avoid uncertainties in classification caused by such cases, coastal types of vegetation were grouped in the same class with other plants having similar spectral properties (grassland). The classification accuracy was performed in order to evaluate possible sources of errors. To this end, images were processed using a rejection probability test, which examines the correctness of the pixel's classification using the chi-squared test. This made it possible to evaluate the classification of different types of land cover types over the study area using a traditional approach of GRASS GIS based on automated clustering (Figure 6).

### 4.3. Interpretation and Data Analysis

Major land cover types in Djibouti include xeric shrubland, bare land areas in mountainous deserts, estuary of the Gulf of Tadjoura, cropland areas and sparse vegetation. Multispectral bands of the images were analysed and compared using a spectral diagram and spatial analysis. The main land cover maps of Djibouti generated using reclassified scenes and ML-based analysis for the images is displayed in Figures 7 and 8.

The reclassified images are shown in Figure 7. Here, the areas of water are merged into one class, which improved the classification. Nevertheless, the neighbouring regions containing vegetation types with similar pixel reflectance were misclassified, which required the ML approach for image processing. A fluctuation in Lake Assal visible on the images is related to climate effects and geothermal activities, as is also mentioned in previous studies. This study revealed that similar to the fluctuations of saline sabkhas in the northern Sahara [117,118], the saline lake, Lake Assal, also witnessed slight changes in the water surface and extent of the lacustrine coasts. Nevertheless, the nature of such physiographic variations is more related to the geologic origin [119].
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Reclassified satellite images after reclassification used by the 'r.reclass' technique. Images are from 2019, 2021 and 2023 covering Djibouti. Background: GEBCO. Software: GRASS GIS. Mapping source: Author.</figcaption></figure>

Landscape change detection was performed through a comparison among the land cover types detected on the ML-based classified satellite images obtained on different dates. The classification of the images is based on an automated evaluation of the spectral reflectance values of the individual pixels on the images, which are assigned to different land cover classes. The comparative analysis of the classified images enabled us to detect landscape changes. The observed transitions of categories were the loss of intact vegetation areas; fluctuations in the saline lake, Lake Assal; gains in vegetation and coastal areas; and secondary vegetation growth. These transitions correspond to the predominant patterns of vegetation response to climate warming and growth annual temperatures, which are visible when comparing central and coastal regions of Djibouti, where the latter have been influenced by the marine climate, especially in the coastal areas of the Gulf of Tadjoura.

The traditional (maximal likelihood) and ML classification maps were compared for 2015, 2019, 2021 and 2023. Comparing the changes on classified images from different dates shows the dynamics of the landscapes over Djibouti from 2015 until 2023. The differences between the land cover or land use types are presented in Figures 5, 7 and 8. Figure 5 shows the results of the classification made using the maximal likelihood approach; Figure 6 shows the accuracy assessment; Figure 7 shows the reclassified maps; and Figure 8 shows the results of the classification performed using the ML approach.

The comparison of the classified maps revealed that natural vegetation, such as herbaceous and mosaic shrubland, cropland and grassland, has decreased while the sparse vegetation and bare soil areas typical of desert and semi-desert areas in Djibouti have increased. Changes in land cover types in the coastal region of Bab-el-Mandeb and central Djibouti were detected and visualised in the target landscapes over the period from 2015 to 2023. The images were compared based on different dates to assess the landscape dynamics, which is defined by the evaluation of land cover types through the identification of various land cover patches, in different regions of Djibouti.
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Image classification based on the ML gradient boosting classifier algorithm applied for Landsat scenes covering Djibouti in 2015, 2019, 2021 and 2023. Background relief: GEBCO. Software: GRASS GIS. Mapping source: Author.</figcaption></figure>

**5. Conclusions**

This paper demonstrated the efficiency of RS data processing using ML methods for visualising land cover changes. Using this approach, the dynamics of land cover types has been detected by identifying differences in the spectral reflectance of pixels on images assigned to different categories by ML algorithms of gradient boosting. This method was tested, and has been explained and presented in the form of a series of new maps covering Djibouti. Interdisciplinary problems in geographical sciences often require decisions to be made by diversified approaches. ML algorithms, such as the gradient boosting classifier, improved mapping workflow through programming to optimise RS data classification. When supplementary modules of GIS are used, such tools point out the shortcomings of traditional methods in geoinformatics. The development of cartographic tools for RS data processing using ML is a challenging task. As technologies mature along with the development of programming algorithms, cartographic modelling and the analysis of landscape changes improve by integrating applications of scripting methods. Here, we have demonstrated the use of such tools using Python's Scikit-Learn library, adapted to GRASS GIS.

Scripting cartographic tools enable us to highlight salient aspects of environmental dynamics through the automation of mapping and data visualisation. They support geospatial data processing, modelling, visualisation and interpretation. For such geologically complex areas, integrated programming methods assist in handling the issue of data processing in a spatially expressive manner. The effectiveness of ML is explained through the advanced algorithms of image classification and automated workflow. Hence, the use of ML tools is essential for implementing the workflow of processing geographic data to reveal and visualise environmental problems. Their applications demonstrate a prominent role of scripting for geospatial data handling and environmental analysis. As ML algorithms are developed and implemented in RS and mapping, their performance supports the computa-
tional complexity of the tasks involved in cartographic data processing. Many case studies on geographical analysis and mapping, which use conventional GISs, for example, require highly computational and costly mapping workload due to diversified cartographic tasks. In such cases, the use of ML, most notably gradient boosting, is an essential solution to optimising cartographic workflow.

Future research can use the thematic maps of land cover types and continue with the analysis of neighbouring regions of the Bab-el-Mandeb Straight for environmental monitoring. Moreover, it can enhance the thematic and topical research of East Africa, particularly in relation to the geological exploration of Djibouti, environmental analysis, and geophysical monitoring of tectonically active regions within the Afar Triple Junction. Furthermore, to continue this study, future works can adapt ML techniques of GRASS GIS for land cover monitoring of other regions over an extended period to analyse landscape dynamics using time series analysis in regions around the Red Sea, East Africa.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data are available in the GitHub repository: https://github.com/paulinelemenkova/GRASS_ML_Djibouti, accessed on 12 June 2025.

**Acknowledgments:** The author thanks the reviewers for reading and reviewing of this manuscript.

**Conflicts of Interest:** The author declares no conflicts of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

<table><tr><td>API</td><td>Application Programming Interface</td></tr><tr><td>CNN</td><td>Convolutional Neural Network</td></tr><tr><td>EO</td><td>Earth Observation</td></tr><tr><td>FAO</td><td>Food and Agriculture Organisation</td></tr><tr><td>GEBCO</td><td>General Bathymetric Chart of the Oceans</td></tr><tr><td>GIS</td><td>Geographic Information System</td></tr><tr><td>GMT</td><td>Generic Mapping Tool</td></tr><tr><td>GRASS GIS</td><td>Geographic Resources Analysis Support System GIS</td></tr><tr><td>KC</td><td>Kappa Coefficient</td></tr><tr><td>LCCS</td><td>Land Cover Classification System</td></tr><tr><td>OLI</td><td>Operational Land Imager</td></tr><tr><td>ML</td><td>Machine Learning</td></tr><tr><td>OA</td><td>Overall Accuracy</td></tr><tr><td>QGIS</td><td>Quantum Geographic Information System</td></tr><tr><td>RF</td><td>Random Forest</td></tr><tr><td>RS</td><td>Remote Sensing</td></tr><tr><td>SVM</td><td>Support Vector Machine</td></tr><tr><td>TIRS</td><td>Thermal Infrared Sensor</td></tr><tr><td>TIFF</td><td>Tag Image File Format</td></tr><tr><td>USGS</td><td>United States Geological Survey</td></tr><tr><td>UTM</td><td>Universal Transverse Mercator</td></tr><tr><td>WGS84</td><td>World Geodetic System 1984</td></tr><tr><td>WRS</td><td>Worldwide Reference System</td></tr></table>
# Appendix A. Initial Means for Each Multispectral Band Before Data Processing

<table><caption>**Table A1.** Initial mean values for pixels in spectral bands assigned to target land cover types in Djibouti for the image obtained on 2 July 2015, before iterative classification process.</caption><thead><tr><th rowspan="2">Class</th><th rowspan="2">1-Coastal Aerosol</th><th colspan="6">Initial Mean Values for Multispectral Bands of Landsat: 2015</th></tr><tr><th>2-Blue</th><th>3-Green</th><th>4-Red</th><th>5-NIR</th><th>6-SWIR-1</th><th>7-SWIR-2</th></tr></thead><tbody><tr><td>1</td><td>8103.97</td><td>8671.63</td><td>9562.73</td><td>9665.43</td><td>10,102</td><td>10,216.9</td><td>9598.71</td></tr><tr><td>2</td><td>8421.69</td><td>9028.36</td><td>10,061.5</td><td>10,389.8</td><td>10,972.8</td><td>11,250.5</td><td>10,522.6</td></tr><tr><td>3</td><td>8739.4</td><td>9385.08</td><td>10,560.2</td><td>11,114.2</td><td>11,843.6</td><td>12,284.2</td><td>11,446.6</td></tr><tr><td>4</td><td>9057.11</td><td>9741.81</td><td>11,059</td><td>11,838.6</td><td>12,714.5</td><td>13,317.8</td><td>12,370.5</td></tr><tr><td>5</td><td>9374.82</td><td>10,098.5</td><td>11,557.7</td><td>12,563</td><td>13,585.3</td><td>14,351.5</td><td>13,294.4</td></tr><tr><td>6</td><td>9692.54</td><td>10,455.3</td><td>12,056.4</td><td>13,287.4</td><td>14,456.2</td><td>15,385.1</td><td>14,218.3</td></tr><tr><td>7</td><td>10,010.2</td><td>10,812</td><td>12,555.2</td><td>14,011.8</td><td>15,327</td><td>16,418.8</td><td>15,142.3</td></tr><tr><td>8</td><td>10,328</td><td>11,168.7</td><td>13,053.9</td><td>14,736.2</td><td>16,197.9</td><td>17,452.4</td><td>16,066.2</td></tr><tr><td>9</td><td>10,645.7</td><td>11,525.4</td><td>13,552.7</td><td>15,460.6</td><td>17,068.7</td><td>18,486.1</td><td>16,990.1</td></tr><tr><td>10</td><td>10,963.4</td><td>11,882.1</td><td>14,051.4</td><td>16,185</td><td>17,939.6</td><td>19,519.7</td><td>17,914</td></tr></tbody></table>

<table><caption>**Table A2.** Initial mean values for pixels in spectral bands assigned to target land cover types in Djibouti for the image obtained on 23 July 2017, before iterative classification process.</caption><thead><tr><th rowspan="2">Class</th><th rowspan="2">1-Coastal Aerosol</th><th colspan="6">Initial Mean Values for Multispectral Bands of Landsat: 2017</th></tr><tr><th>2-Blue</th><th>3-Green</th><th>4-Red</th><th>5-NIR</th><th>6-SWIR-1</th><th>7-SWIR-2</th></tr></thead><tbody><tr><td>1</td><td>8074.96</td><td>8660.92</td><td>9599.33</td><td>9725.17</td><td>10,161.2</td><td>10,291.9</td><td>9670.42</td></tr><tr><td>2</td><td>8390.09</td><td>9009.56</td><td>10,079.2</td><td>10,419.5</td><td>11,011.6</td><td>11,300.7</td><td>10,555.8</td></tr><tr><td>3</td><td>8705.22</td><td>9358.2</td><td>10,559.1</td><td>11,113.9</td><td>11,861.9</td><td>12,309.4</td><td>11,441.2</td></tr><tr><td>4</td><td>9020.35</td><td>9706.85</td><td>11,039</td><td>11,808.3</td><td>12,712.3</td><td>13,318.1</td><td>12,326.7</td></tr><tr><td>5</td><td>9335.49</td><td>10,055.5</td><td>11,518.9</td><td>12,502.6</td><td>13,562.7</td><td>14,326.9</td><td>13,212.1</td></tr><tr><td>6</td><td>9650.62</td><td>10,404.1</td><td>11,998.8</td><td>13,197</td><td>14,413.1</td><td>15,335.6</td><td>14,097.5</td></tr><tr><td>7</td><td>9965.75</td><td>10,752.8</td><td>12,478.7</td><td>13,891.4</td><td>15,263.5</td><td>16,344.3</td><td>14,982.9</td></tr><tr><td>8</td><td>10,280.9</td><td>11,101.4</td><td>12,958.6</td><td>14,585.7</td><td>16,113.8</td><td>17,353.1</td><td>15,868.3</td></tr><tr><td>9</td><td>10,596</td><td>11,450.1</td><td>13,438.5</td><td>15,280.1</td><td>16,964.2</td><td>18,361.8</td><td>16,753.7</td></tr><tr><td>10</td><td>10,911.1</td><td>11,798.7</td><td>13,918.4</td><td>15,974.4</td><td>17,814.6</td><td>19,370.5</td><td>17,639.1</td></tr></tbody></table>

<table><caption>**Table A3.** Initial mean values for pixels in spectral bands assigned to target land cover types in Djibouti for the image obtained on 19 August 2021, before iterative classification process.</caption><thead><tr><th rowspan="2">Class</th><th rowspan="2">1-Coastal Aerosol</th><th colspan="6">Initial Mean Values for Multispectral Bands of Landsat: 2021</th></tr><tr><th>2-Blue</th><th>3-Green</th><th>4-Red</th><th>5-NIR</th><th>6-SWIR-1</th><th>7-SWIR-2</th></tr></thead><tbody><tr><td>1</td><td>8317.6</td><td>8806.66</td><td>9565.01</td><td>9718.97</td><td>10,254</td><td>10,232.6</td><td>9622.06</td></tr><tr><td>2</td><td>8578.74</td><td>9106.8</td><td>10,005.6</td><td>10,359.3</td><td>11,061</td><td>11,207.5</td><td>10,480.5</td></tr><tr><td>3</td><td>8839.87</td><td>9406.95</td><td>10,446.1</td><td>10,999.7</td><td>11,867.9</td><td>12,182.4</td><td>11,338.9</td></tr><tr><td>4</td><td>9101</td><td>9707.09</td><td>10,886.7</td><td>11,640.1</td><td>12,674.8</td><td>13,157.4</td><td>12,197.4</td></tr><tr><td>5</td><td>9362.13</td><td>10,007.2</td><td>11,327.3</td><td>12,280.4</td><td>13,481.8</td><td>14,132.3</td><td>13,055.8</td></tr><tr><td>6</td><td>9623.26</td><td>10,307.4</td><td>11,767.8</td><td>12,920.8</td><td>14,288.7</td><td>15,107.2</td><td>13,914.3</td></tr><tr><td>7</td><td>9884.39</td><td>10,607.5</td><td>12,208.4</td><td>13,561.1</td><td>15,095.7</td><td>16,082.1</td><td>14,772.7</td></tr><tr><td>8</td><td>10,145.5</td><td>10,907.7</td><td>12,649</td><td>14,201.5</td><td>15,902.6</td><td>17,057.1</td><td>15,631.2</td></tr><tr><td>9</td><td>10,406.7</td><td>11,207.8</td><td>13,089.5</td><td>14,841.9</td><td>16,709.6</td><td>18,032</td><td>16,489.6</td></tr><tr><td>10</td><td>10,667.8</td><td>11,507.9</td><td>13,530.1</td><td>15,482.2</td><td>17,516.5</td><td>19,006.9</td><td>17,348</td></tr></tbody></table>

<table><caption>**Table A4.** Initial mean values for pixels in spectral bands assigned to target land cover types in Djibouti for the image obtained on 16 July 2023, before iterative classification process.</caption><thead><tr><th rowspan="2">Class</th><th rowspan="2">1-Coastal Aerosol</th><th colspan="6">Initial Mean Values for Multispectral Bands of Landsat: 2023</th></tr><tr><th>2-Blue</th><th>3-Green</th><th>4-Red</th><th>5-NIR</th><th>6-SWIR-1</th><th>7-SWIR-2</th></tr></thead><tbody><tr><td>1</td><td>8504.75</td><td>9131.19</td><td>9993.94</td><td>10,161.4</td><td>10,600.6</td><td>10,554.5</td><td>9880.58</td></tr><tr><td>2</td><td>8747.11</td><td>9404.01</td><td>10,394.1</td><td>10,760.9</td><td>11,378.1</td><td>11,500.1</td><td>10,685.6</td></tr><tr><td>3</td><td>8989.47</td><td>9676.83</td><td>10,794.3</td><td>11,360.4</td><td>12,155.5</td><td>12,445.8</td><td>11,490.6</td></tr><tr><td>4</td><td>9231.83</td><td>9949.65</td><td>11,194.5</td><td>11,960</td><td>12,933</td><td>13,391.4</td><td>12,295.6</td></tr><tr><td>5</td><td>9474.19</td><td>10,222.5</td><td>11,594.7</td><td>12,559.5</td><td>13,710.4</td><td>14,337</td><td>13,100.6</td></tr><tr><td>6</td><td>9716.55</td><td>10,495.3</td><td>11,994.9</td><td>13,159.1</td><td>14,487.8</td><td>15,282.7</td><td>13,905.5</td></tr><tr><td>7</td><td>9958.91</td><td>10,768.1</td><td>12,395.1</td><td>13,758.6</td><td>15,265.3</td><td>16,228.3</td><td>14,710.5</td></tr><tr><td>8</td><td>10,201.3</td><td>11,040.9</td><td>12,795.3</td><td>14,358.1</td><td>16,042.7</td><td>17,173.9</td><td>15,515.5</td></tr><tr><td>9</td><td>10,443.6</td><td>11,313.7</td><td>13,195.5</td><td>14,957.7</td><td>16,820.1</td><td>18,119.6</td><td>16,320.5</td></tr><tr><td>10</td><td>10,686</td><td>11,586.6</td><td>13,595.7</td><td>15,557.2</td><td>17,597.6</td><td>19,065.2</td><td>17,125.5</td></tr></tbody></table>
## Appendix B. Iterative Cycles of Image Classification

<table><caption>Table A5. Iterative process of pixel assignment to target classes during image classification: 2015.</caption><thead><tr><th rowspan="2">2015</th><th rowspan="2">Stable Points</th><th colspan="10">Pixel Distribution by Categories of Land Cover Classes: 2015</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>Iteration 1</td><td>92.24%</td><td>1672</td><td>112</td><td>50</td><td>559</td><td>1110</td><td>860</td><td>613</td><td>390</td><td>635</td><td>929</td></tr><tr><td>Iteration 2</td><td>93.78%</td><td>1568</td><td>213</td><td>75</td><td>584</td><td>1057</td><td>850</td><td>622</td><td>494</td><td>647</td><td>820</td></tr><tr><td>Iteration 3</td><td>96.28%</td><td>1538</td><td>242</td><td>97</td><td>596</td><td>1026</td><td>843</td><td>658</td><td>541</td><td>626</td><td>763</td></tr><tr><td>Iteration 4</td><td>96.78%</td><td>1531</td><td>249</td><td>117</td><td>614</td><td>994</td><td>840</td><td>691</td><td>575</td><td>599</td><td>720</td></tr><tr><td>Iteration 5</td><td>96.97%</td><td>1531</td><td>249</td><td>137</td><td>622</td><td>980</td><td>836</td><td>714</td><td>610</td><td>568</td><td>683</td></tr><tr><td>Iteration 6</td><td>97.10%</td><td>1531</td><td>249</td><td>166</td><td>626</td><td>963</td><td>833</td><td>726</td><td>643</td><td>537</td><td>656</td></tr><tr><td>Iteration 7</td><td>97.13%</td><td>1531</td><td>250</td><td>194</td><td>632</td><td>945</td><td>834</td><td>732</td><td>669</td><td>512</td><td>631</td></tr><tr><td>Iteration 8</td><td>97.11%</td><td>1531</td><td>250</td><td>219</td><td>650</td><td>930</td><td>827</td><td>744</td><td>668</td><td>498</td><td>613</td></tr><tr><td>Iteration 9</td><td>97.52%</td><td>1531</td><td>250</td><td>243</td><td>658</td><td>936</td><td>819</td><td>732</td><td>674</td><td>480</td><td>607</td></tr><tr><td>Iteration 10</td><td>97.59%</td><td>1531</td><td>250</td><td>261</td><td>672</td><td>941</td><td>810</td><td>734</td><td>659</td><td>472</td><td>600</td></tr><tr><td>Iteration 11</td><td>97.75%</td><td>1531</td><td>251</td><td>279</td><td>681</td><td>954</td><td>798</td><td>727</td><td>650</td><td>462</td><td>597</td></tr><tr><td>Iteration 12</td><td>97.79%</td><td>1531</td><td>251</td><td>291</td><td>704</td><td>968</td><td>775</td><td>720</td><td>640</td><td>456</td><td>594</td></tr><tr><td>Iteration 13</td><td>97.95%</td><td>1531</td><td>251</td><td>305</td><td>729</td><td>964</td><td>768</td><td>710</td><td>629</td><td>450</td><td>593</td></tr><tr><td>Iteration 14</td><td>98.12%</td><td>1531</td><td>251</td><td>314</td><td>755</td><td>964</td><td>761</td><td>699</td><td>616</td><td>448</td><td>591</td></tr></tbody></table>

<table><caption>Table A6. Iterative process of pixel assignment to target classes during image classification: 2017.</caption><thead><tr><th rowspan="2">2017</th><th rowspan="2">Stable Points</th><th colspan="10">Pixel Distribution by Categories of Land Cover Classes: 2017</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>Iteration 1</td><td>91.34%</td><td>1698</td><td>165</td><td>76</td><td>574</td><td>1121</td><td>829</td><td>516</td><td>427</td><td>655</td><td>959</td></tr><tr><td>Iteration 2</td><td>93.52%</td><td>1584</td><td>263</td><td>117</td><td>620</td><td>1047</td><td>811</td><td>522</td><td>531</td><td>667</td><td>858</td></tr><tr><td>Iteration 3</td><td>95.37%</td><td>1518</td><td>316</td><td>137</td><td>669</td><td>996</td><td>788</td><td>553</td><td>605</td><td>640</td><td>798</td></tr><tr><td>Iteration 4</td><td>96.70%</td><td>1486</td><td>341</td><td>140</td><td>699</td><td>977</td><td>776</td><td>578</td><td>668</td><td>605</td><td>750</td></tr><tr><td>Iteration 5</td><td>97.15%</td><td>1468</td><td>357</td><td>134</td><td>718</td><td>972</td><td>771</td><td>611</td><td>707</td><td>575</td><td>707</td></tr><tr><td>Iteration 6</td><td>97.28%</td><td>1461</td><td>360</td><td>131</td><td>729</td><td>975</td><td>783</td><td>641</td><td>712</td><td>561</td><td>667</td></tr><tr><td>Iteration 7</td><td>97.51%</td><td>1456</td><td>362</td><td>123</td><td>741</td><td>986</td><td>798</td><td>651</td><td>730</td><td>528</td><td>645</td></tr><tr><td>Iteration 8</td><td>97.62%</td><td>1453</td><td>363</td><td>118</td><td>745</td><td>1016</td><td>795</td><td>670</td><td>725</td><td>510</td><td>625</td></tr><tr><td>Iteration 9</td><td>98.15%</td><td>1452</td><td>361</td><td>120</td><td>746</td><td>1031</td><td>813</td><td>672</td><td>714</td><td>502</td><td>609</td></tr></tbody></table>

<table><caption>Table A7. Iterative process of pixels' assignment to target classes during image classification: 2021.</caption><thead><tr><th rowspan="2">2021</th><th rowspan="2">Stable Points</th><th colspan="10">Pixels' Distribution by Categories of Land Cover Classes: 2021</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>Iteration 1</td><td>91.92%</td><td>1682</td><td>142</td><td>128</td><td>663</td><td>1062</td><td>784</td><td>485</td><td>371</td><td>626</td><td>988</td></tr><tr><td>Iteration 2</td><td>90.90%</td><td>1385</td><td>431</td><td>174</td><td>664</td><td>1017</td><td>766</td><td>497</td><td>467</td><td>659</td><td>871</td></tr><tr><td>Iteration 3</td><td>91.33%</td><td>1023</td><td>791</td><td>193</td><td>673</td><td>989</td><td>758</td><td>516</td><td>549</td><td>634</td><td>805</td></tr><tr><td>Iteration 4</td><td>93.64%</td><td>785</td><td>1029</td><td>198</td><td>684</td><td>985</td><td>752</td><td>540</td><td>587</td><td>621</td><td>750</td></tr><tr><td>Iteration 5</td><td>95.74%</td><td>666</td><td>1147</td><td>203</td><td>690</td><td>984</td><td>761</td><td>557</td><td>614</td><td>596</td><td>713</td></tr><tr><td>Iteration 6</td><td>96.90%</td><td>612</td><td>1201</td><td>208</td><td>694</td><td>984</td><td>771</td><td>573</td><td>641</td><td>556</td><td>691</td></tr><tr><td>Iteration 7</td><td>97.76%</td><td>596</td><td>1217</td><td>211</td><td>698</td><td>985</td><td>780</td><td>593</td><td>647</td><td>537</td><td>667</td></tr><tr><td>Iteration 8</td><td>97.88%</td><td>589</td><td>1224</td><td>216</td><td>701</td><td>990</td><td>795</td><td>603</td><td>645</td><td>513</td><td>655</td></tr><tr><td>Iteration 9</td><td>98.12%</td><td>588</td><td>1225</td><td>223</td><td>699</td><td>1007</td><td>798</td><td>608</td><td>644</td><td>495</td><td>644</td></tr></tbody></table>

<table><caption>Table A8. Iterative process of pixel assignment to target classes during image classification: 2023.</caption><thead><tr><th rowspan="2">2023</th><th rowspan="2">Stable Points</th><th colspan="10">Pixel Distribution by Categories of Land Cover Classes: 2023</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><td>Iteration 1</td><td>91.95%</td><td>1665</td><td>98</td><td>82</td><td>574</td><td>1139</td><td>882</td><td>567</td><td>370</td><td>632</td><td>1032</td></tr><tr><td>Iteration 2</td><td>91.27%</td><td>1393</td><td>367</td><td>114</td><td>611</td><td>1068</td><td>871</td><td>570</td><td>458</td><td>689</td><td>900</td></tr><tr><td>Iteration 3</td><td>92.44%</td><td>1134</td><td>623</td><td>138</td><td>641</td><td>1027</td><td>851</td><td>595</td><td>538</td><td>662</td><td>832</td></tr><tr><td>Iteration 4</td><td>94.08%</td><td>939</td><td>818</td><td>156</td><td>653</td><td>1012</td><td>837</td><td>618</td><td>588</td><td>644</td><td>776</td></tr><tr><td>Iteration 5</td><td>96.51%</td><td>872</td><td>885</td><td>173</td><td>661</td><td>1004</td><td>828</td><td>638</td><td>609</td><td>630</td><td>741</td></tr><tr><td>Iteration 6</td><td>97.26%</td><td>847</td><td>910</td><td>187</td><td>668</td><td>1003</td><td>820</td><td>662</td><td>606</td><td>625</td><td>713</td></tr><tr><td>Iteration 7</td><td>97.74%</td><td>842</td><td>915</td><td>200</td><td>676</td><td>996</td><td>829</td><td>663</td><td>621</td><td>606</td><td>693</td></tr><tr><td>Iteration 8</td><td>97.77%</td><td>842</td><td>915</td><td>216</td><td>681</td><td>996</td><td>825</td><td>671</td><td>635</td><td>583</td><td>677</td></tr><tr><td>Iteration 9</td><td>97.64%</td><td>842</td><td>915</td><td>237</td><td>683</td><td>994</td><td>829</td><td>666</td><td>648</td><td>569</td><td>658</td></tr><tr><td>Iteration 10</td><td>97.88%</td><td>842</td><td>915</td><td>250</td><td>700</td><td>983</td><td>830</td><td>666</td><td>656</td><td>558</td><td>641</td></tr><tr><td>Iteration 11</td><td>97.95%</td><td>842</td><td>915</td><td>264</td><td>711</td><td>982</td><td>826</td><td>667</td><td>655</td><td>558</td><td>621</td></tr><tr><td>Iteration 12</td><td>97.73%</td><td>842</td><td>915</td><td>282</td><td>719</td><td>982</td><td>822</td><td>668</td><td>655</td><td>557</td><td>599</td></tr><tr><td>Iteration 13</td><td>97.74%</td><td>842</td><td>915</td><td>303</td><td>721</td><td>980</td><td>832</td><td>664</td><td>650</td><td>549</td><td>585</td></tr><tr><td>Iteration 14</td><td>97.80%</td><td>842</td><td>915</td><td>332</td><td>710</td><td>980</td><td>837</td><td>672</td><td>639</td><td>545</td><td>569</td></tr><tr><td>Iteration 15</td><td>97.83%</td><td>842</td><td>915</td><td>352</td><td>723</td><td>966</td><td>844</td><td>672</td><td>634</td><td>533</td><td>560</td></tr><tr><td>Iteration 16</td><td>98.13%</td><td>842</td><td>915</td><td>367</td><td>733</td><td>965</td><td>840</td><td>672</td><td>627</td><td>533</td><td>547</td></tr></tbody></table>
## Appendix C. Computed Class Means for Land Cover Classes in the Region of Lake Assal, Djibouti, East Africa

<table><thead><tr><th></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><th>1</th><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>2</th><td>1.5</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>3</th><td>5.2</td><td>3.0</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>4</th><td>6.9</td><td>4.3</td><td>0.8</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>5</th><td>7.8</td><td>5.0</td><td>1.5</td><td>0.8</td><td>0</td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>6</th><td>7.9</td><td>5.3</td><td>2.2</td><td>1.7</td><td>0.9</td><td>0</td><td></td><td></td><td></td><td></td></tr><tr><th>7</th><td>6.2</td><td>4.5</td><td>2.4</td><td>2.0</td><td>1.5</td><td>0.8</td><td>0</td><td></td><td></td><td></td></tr><tr><th>8</th><td>7.4</td><td>5.5</td><td>3.4</td><td>3.1</td><td>2.6</td><td>1.8</td><td>0.8</td><td>0</td><td></td><td></td></tr><tr><th>9</th><td>8.6</td><td>6.7</td><td>4.7</td><td>4.5</td><td>4.0</td><td>3.2</td><td>1.9</td><td>1.1</td><td>0</td><td></td></tr><tr><th>10</th><td>11.0</td><td>8.6</td><td>6.5</td><td>6.5</td><td>6.0</td><td>5.0</td><td>3.3</td><td>2.4</td><td>1.2</td><td>0</td></tr></tbody></table>

(A1)

Equation (A1): Class separability matrices for classified images of Djibouti: 2015 and 2017.

<table><thead><tr><th></th><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th><th>10</th></tr></thead><tbody><tr><th>1</th><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>2</th><td>1.6</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>3</th><td>4.1</td><td>2.5</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>4</th><td>5.0</td><td>3.4</td><td>0.7</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>5</th><td>6.3</td><td>4.6</td><td>1.4</td><td>0.7</td><td>0</td><td></td><td></td><td></td><td></td><td></td></tr><tr><th>6</th><td>6.5</td><td>4.9</td><td>2.0</td><td>1.4</td><td>0.8</td><td>0</td><td></td><td></td><td></td><td></td></tr><tr><th>7</th><td>5.2</td><td>4.1</td><td>2.1</td><td>1.7</td><td>1.3</td><td>0.7</td><td>0</td><td></td><td></td><td></td></tr><tr><th>8</th><td>6.9</td><td>5.6</td><td>3.3</td><td>2.9</td><td>2.5</td><td>1.8</td><td>0.8</td><td>0</td><td></td><td></td></tr><tr><th>9</th><td>5.7</td><td>4.9</td><td>3.2</td><td>2.9</td><td>2.7</td><td>2.2</td><td>1.4</td><td>0.9</td><td>0</td><td></td></tr><tr><th>10</th><td>8.4</td><td>7.2</td><td>5.2</td><td>4.9</td><td>4.8</td><td>4.0</td><td>2.7</td><td>2.1</td><td>0.9</td><td></td></tr></tbody></table>

(A2)

Equation (A2): Class separability matrices for classified images of Djibouti: 2021 and 2023.

## Appendix D. Computed Class Means for Land Cover Classes in the Region of Lake Assal, Djibouti, East Africa

<table><caption>**Table A9.** Computed class means of 10 land cover types in Lake Assal, Djibouti. The computations are based on Digital Numbers (DNs) of pixels in the evaluated landscapes. Data: multispectral Landsat 8–9 OLI/TIRS images, band (1 to 7) for years 2015 to 2023.</caption><thead><tr><th>Class</th><th>Band 1</th><th>Band 2</th><th>Band 3</th><th>Band 4</th><th>Band 5</th><th>Band 6</th><th>Band 7</th></tr></thead><tbody><tr><td colspan="8">2015</td></tr><tr><td>1</td><td>7394.23</td><td>7878.92</td><td>8615.8</td><td>8244.47</td><td>8385.55</td><td>8786</td><td>8586.51</td></tr><tr><td>2</td><td>8998.3</td><td>9445.8</td><td>9981.23</td><td>9575.59</td><td>9422.92</td><td>9284.88</td><td>8887.37</td></tr><tr><td>3</td><td>9230.46</td><td>9913.19</td><td>11,138.2</td><td>11,958.2</td><td>12,864.5</td><td>12,862.9</td><td>11,816.5</td></tr><tr><td>4</td><td>9430.8</td><td>10,146.6</td><td>11,465.7</td><td>12,446</td><td>13,497.2</td><td>13,806.2</td><td>12,619</td></tr><tr><td>5</td><td>9547.38</td><td>10,303</td><td>11,752.4</td><td>12,912.6</td><td>14,092.7</td><td>14,641.3</td><td>13,360.3</td></tr><tr><td>6</td><td>9821.76</td><td>10,624.5</td><td>12,235.6</td><td>13,635.1</td><td>14,896.9</td><td>15,478.6</td><td>14,133.9</td></tr><tr><td>7</td><td>10,246.8</td><td>11,081.6</td><td>12,862.4</td><td>14,552</td><td>16,012.4</td><td>16,653.6</td><td>15,224.3</td></tr><tr><td>8</td><td>10,661.7</td><td>11,544</td><td>13,560.5</td><td>15,543.5</td><td>17,259.8</td><td>18,626.3</td><td>17,001.6</td></tr><tr><td>9</td><td>11,251.3</td><td>12,246.2</td><td>14,758.7</td><td>17,309</td><td>19,326.1</td><td>21,064.2</td><td>19,184.3</td></tr><tr><td>10</td><td>11,881.4</td><td>12,946.1</td><td>15,865.1</td><td>18,729.9</td><td>20,856</td><td>24,277.4</td><td>22,627.1</td></tr><tr><td colspan="8">2017</td></tr><tr><td>1</td><td>7371.5</td><td>7875.98</td><td>8664.69</td><td>8356.38</td><td>8565.28</td><td>8995.22</td><td>8778.1</td></tr><tr><td>2</td><td>8858.16</td><td>9357.22</td><td>9951.7</td><td>9564.73</td><td>9444.54</td><td>9371.12</td><td>8992.18</td></tr><tr><td>3</td><td>8945.89</td><td>9777.33</td><td>11,054.5</td><td>11,333.4</td><td>11,764</td><td>11,571.8</td><td>10,743</td></tr><tr><td>4</td><td>9321.25</td><td>10,018.9</td><td>11,324</td><td>12,242.3</td><td>13,244.4</td><td>13,468.5</td><td>12,268.8</td></tr><tr><td>5</td><td>9496.17</td><td>10,212.7</td><td>11,611.2</td><td>12,696.3</td><td>13,855.2</td><td>14,375.8</td><td>13,058.1</td></tr><tr><td>6</td><td>9715.96</td><td>10,479.6</td><td>12,049.1</td><td>13,373.4</td><td>14,635.2</td><td>15,185.6</td><td>13,795</td></tr><tr><td>7</td><td>10,087</td><td>10,877.3</td><td>12,605.2</td><td>14,223.3</td><td>15,670.3</td><td>16,380</td><td>14,930.2</td></tr><tr><td>8</td><td>10,526.8</td><td>11,378.8</td><td>13,341.7</td><td>15,248</td><td>16,971.7</td><td>18,216.1</td><td>16,582.3</td></tr><tr><td>9</td><td>11,127.8</td><td>12,082.7</td><td>14,519.5</td><td>16,951.3</td><td>19,046.9</td><td>20,773.9</td><td>18,771.8</td></tr><tr><td>10</td><td>11,728.2</td><td>12,813.2</td><td>15,676</td><td>18,406.1</td><td>20,794.3</td><td>24,015.7</td><td>22,082.7</td></tr></tbody></table>
<table><caption>Table A9. *Cont.*</caption><thead><tr><th>Class</th><th>Band 1</th><th>Band 2</th><th>Band 3</th><th>Band 4</th><th>Band 5</th><th>Band 6</th><th>Band 7</th></tr></thead><tbody><tr><td colspan="8">2021</td></tr><tr><td>1</td><td>7281.64</td><td>7688.52</td><td>8323.58</td><td>8019.43</td><td>8124.06</td><td>8480.67</td><td>8350.31</td></tr><tr><td>2</td><td>8578.45</td><td>8990.99</td><td>9497.29</td><td>9305.59</td><td>9360.4</td><td>9422.46</td><td>9107.76</td></tr><tr><td>3</td><td>9032.4</td><td>9610.46</td><td>10,681.6</td><td>11,350.9</td><td>12,512.9</td><td>12,392.4</td><td>11,404.3</td></tr><tr><td>4</td><td>9232.84</td><td>9855.5</td><td>11,063</td><td>11,913.4</td><td>13,238.9</td><td>13,353.2</td><td>12,168.5</td></tr><tr><td>5</td><td>9378.87</td><td>10,036.2</td><td>11,347.7</td><td>12,369.6</td><td>13,788.5</td><td>14,219</td><td>12,929</td></tr><tr><td>6</td><td>9563.91</td><td>10,256.7</td><td>11,715.6</td><td>12,988.2</td><td>14,561.6</td><td>15,110.5</td><td>13,758.7</td></tr><tr><td>7</td><td>10,037.5</td><td>10,785.2</td><td>12,429.3</td><td>13,969.7</td><td>15,629.4</td><td>16,072.7</td><td>14,651.6</td></tr><tr><td>8</td><td>10,332.3</td><td>11,120.3</td><td>12,974.7</td><td>14,809.1</td><td>16,765.2</td><td>17,972.9</td><td>16,439.1</td></tr><tr><td>9</td><td>10,934.6</td><td>11,846.4</td><td>14,132.9</td><td>16,407.4</td><td>18,587</td><td>20,167.3</td><td>18,278.5</td></tr><tr><td>10</td><td>11,319.5</td><td>12,358.9</td><td>15,074.5</td><td>17,683.4</td><td>20,102</td><td>23,286.8</td><td>21,439.7</td></tr><tr><td colspan="8">2023</td></tr><tr><td>1</td><td>7473.59</td><td>8079.42</td><td>8857.77</td><td>8532.88</td><td>8642.03</td><td>8940.82</td><td>8737.73</td></tr><tr><td>2</td><td>9179.37</td><td>9560.93</td><td>9889.56</td><td>9619.59</td><td>9502.14</td><td>9260.89</td><td>8893.31</td></tr><tr><td>3</td><td>9285.35</td><td>10,023</td><td>11,304.4</td><td>12,097.2</td><td>13,066.2</td><td>13,055</td><td>11,870.5</td></tr><tr><td>4</td><td>9260.65</td><td>9993.24</td><td>11,324.5</td><td>12,294.7</td><td>13,624</td><td>14,111.4</td><td>12,725.8</td></tr><tr><td>5</td><td>9612.99</td><td>10,407.7</td><td>11,850</td><td>12,914.8</td><td>14,185.6</td><td>14,445.7</td><td>13,030.6</td></tr><tr><td>6</td><td>9612.99</td><td>10,407.7</td><td>11,850</td><td>12,914.8</td><td>14,185.6</td><td>14,445.7</td><td>13,030.6</td></tr><tr><td>7</td><td>10,105.5</td><td>10,948.7</td><td>12,646.3</td><td>14,221.3</td><td>15,852.4</td><td>16,551.5</td><td>14,997.2</td></tr><tr><td>8</td><td>10,551</td><td>11,459.3</td><td>13,401.9</td><td>15,287.1</td><td>17,177.8</td><td>18,217.3</td><td>16,436.4</td></tr><tr><td>9</td><td>10,830.2</td><td>11,812.4</td><td>14,092.6</td><td>16,381.5</td><td>18,713.8</td><td>20,687.4</td><td>18,410.8</td></tr><tr><td>10</td><td>11,191.6</td><td>12,288.7</td><td>14,890.2</td><td>17,469.4</td><td>20,034.6</td><td>23,154.3</td><td>20,893.2</td></tr></tbody></table>

## References

1. Ndehedehe, C. Remote Sensing of Surface Vegetation. In *Satellite Remote Sensing of Terrestrial Hydrology*; Springer International Publishing: Cham, Switzerland, 2022; pp. 131–176. [CrossRef]

2. Foli, B.A.K.; Williams, I.K.; Boakye, A.A.; Azumah, D.M.Y.; Agyekum, K.A.; Wiafe, G. Earth Observation Services in Support of West Africa’s Blue Economy: Coastal Resilience and Climate Impacts. *Remote Sens. Earth Syst. Sci.* **2022**, *5*, 59–70. [CrossRef]

3. Ifejika Speranza, C.; Akinyemi, F.O.; Baratoux, D.; Benveniste, J.; Ceperley, N.; Driouech, F.; Helmschrot, J. Enhancing the Uptake of Earth Observation Products and Services in Africa Through a Multi-level Transdisciplinary Approach. *Surv. Geophys.* **2023**, *44*, 7–41. [CrossRef] [PubMed]

4. Springer, A.; Lopez, T.; Owor, M.; Frappart, F.; Stieglitz, T. The Role of Space-Based Observations for Groundwater Resource Monitoring over Africa. *Surv. Geophys.* **2023**, *44*, 123–172. [CrossRef]

5. Jia, N.; Li, Y.; Viña, A.; Cheng, J.; Dou, Y.; Song, Q.; He, L.; Liu, J. Long Image Time Series for Crop Extraction Based on the Automatically Generated Samples Algorithm. In Proceedings of the IGARSS 2023–2023 IEEE International Geoscience and Remote Sensing Symposium, Pasadena, CA, USA, 16–21 July 2023; pp. 3502–3505. [CrossRef]

6. Ngondo, J.; Mango, J.; Liu, R.; Nobert, J.; Dubi, A.; Cheng, H. Land-Use and Land-Cover (LULC) Change Detection and the Implications for Coastal Water Resource Management in the Wami–Ruvu Basin, Tanzania. *Sustainability* **2021**, *13*, 4092. [CrossRef]

7. Natarajan, K.; Latva-Käyrä, P.; Zyadin, A.; Pelkonen, P. New methodological approach for biomass resource assessment in India using GIS application and land use/land cover (LULC) maps. *Renew. Sustain. Energy Rev.* **2016**, *63*, 256–268. [CrossRef]

8. Borana, S.; Yadav, S. Chapter 10-Urban land-use susceptibility and sustainability—Case study. In *Water, Land, and Forest Susceptibility and Sustainability*; Chatterjee, U., Pradhan, B., Kumar, S., Saha, S., Zakwan, M., Fath, B.D., Fiscus, D., Eds.; Academic Press: Cambridge, MA, USA, 2023; Volume 2, pp. 261–286. [CrossRef]

9. Zhu, A.X.; Zhao, F.H.; Liang, P.; Qin, C.Z. Next generation of GIS: Must be easy. *Ann. GIS* **2021**, *27*, 71–86. [CrossRef]

10. Dyson, K.; Nicolau, A.P.; Saah, D.; Clinton, N. Image Manipulation: Bands, Arithmetic, Thresholds, and Masks. In *Cloud-Based Remote Sensing with Google Earth Engine: Fundamentals and Applications*; Springer International Publishing: Cham, Switzerland, 2024; pp. 97–114. [CrossRef]

11. Kochher, R.; Sharma, A. Improved principle component analysis based gray stretch algorithm for landsat image segmentation. In Proceedings of the 2016 2nd International Conference on Next Generation Computing Technologies (NGCT), Dehradun, India, 14–16 October 2016; pp. 765–771. [CrossRef]

12. Ferchichi, A.; Chihaoui, M.; Ferchichi, A. Spatio-temporal modeling of climate change impacts on drought forecast using Generative Adversarial Network: A case study in Africa. *Expert Syst. Appl.* **2024**, *238*, 122211. [CrossRef]
13. Balha, A.; Singh, C.K. Comparison of Maximum Likelihood, Neural Networks, and Random Forests Algorithms in Classifying Urban Landscape. In *Application of Remote Sensing and GIS in Natural Resources and Built Infrastructure Management*; Springer International Publishing: Cham, Switzerland, 2022; pp. 29–38. [CrossRef]

14. Lemenkova, P. Random Forest Classifier Algorithm of Geographic Resources Analysis Support System Geographic Information System for Satellite Image Processing: Case Study of Bight of Sofala, Mozambique. *Coasts* **2024**, *4*, 127–149. [CrossRef]

15. Zhang, D.; Zhao, T.; Xue, S.; Meng, Y.; Wang, W. Diagnosis of Dark Field Partial Slight Discharge of GIS Coaxial Busbar Using KPCA-SVM-Based Heuristic Sequence Comparison Method. In *Proceedings of the 2023 China Automation Congress (CAC)*, Chongqing, China, 17–19 November 2023; pp. 4487–4492. [CrossRef]

16. Suárez-Fernández, G.E.; Martínez-Sánchez, J.; Arias, P. Assessment of vegetation indices for mapping burned areas using a deep learning method and a comprehensive forest fire dataset from Landsat collection. *Adv. Space Res.* **2025**, *75*, 1665–1685. [CrossRef]

17. Huang, F.; Cao, Z.; Guo, J.; Jiang, S.H.; Li, S.; Guo, Z. Comparisons of heuristic, general statistical and machine learning models for landslide susceptibility prediction and mapping. *Catena* **2020**, *191*, 104580. [CrossRef]

18. Arabameri, A.; Arora, A.; Pal, S.C.; Mitra, S.; Saha, A.; Nalivan, O.A.; Panahi, S.; Moayedi, H. K-Fold and State-of-the-Art Metaheuristic Machine Learning Approaches for Groundwater Potential Modelling. *Water Resour. Manag.* **2021**, *35*, 1837–1869. [CrossRef]

19. Hara Sudhan, G.H.; Aravind, R.G.; Thanaraj, K.P. Multispectral Analysis of Satellite Images Using Heuristic Algorithm. In *Proceedings of the 2018 International Conference on Communication and Signal Processing (ICCSP)*, Chennai, India, 3–5 April 2018; pp. 1010–1013. [CrossRef]

20. Vitale, A.; Salvo, C.; Lamonaca, F. A Novel Geospatial Methodology for Measuring and Mapping Spatiotemporal Built-Up Dynamics Based on Google Earth Engine and Unsupervised K-Means Clustering of Multispectral Satellite Imagery. In *Proceedings of the 2024 IEEE International Workshop on Metrology for Living Environment (MetroLivEnv)*, Chania, Greece, 12–14 June 2024; pp. 57–62. [CrossRef]

21. Kirkby, S.; Badcock, J.; Eklund, P. Encoding heuristic knowledge for GIS classification. In *Proceedings of the 1997 IEEE International Conference on Intelligent Processing Systems (Cat. No.97TH8335)*, Beijing, China, 28–31 October 1997; Volume 2, pp. 1143–1145. [CrossRef]

22. Hong, Z.; Bian, F. A Heuristic Approach for Fast Mining Association Rules in Transportation System. In *Proceedings of the 2008 Fifth International Conference on Fuzzy Systems and Knowledge Discovery, Jinan, China, 18–20 October 2008; Volume 2*, pp. 541–545. [CrossRef]

23. Pham, V.D.; Nguyen, Q.H.; Nguyen, H.D.; Pham, V.M.; Vu, V.M.; Bui, Q.T. Convolutional Neural Network—Optimized Moth Flame Algorithm for Shallow Landslide Susceptible Analysis. *IEEE Access* **2020**, *8*, 32727–32736. [CrossRef]

24. Han, S.; Liao, S.; Gao, F.; Wang, B.; Yang, N. Pattern Recognition of UHF Partial Discharge Signals in GIS Based on Deep Learning. In *Proceedings of the 2021 International Conference on Advanced Electrical Equipment and Reliable Operation (AEERO)*, Beijing, China, 15–17 October 2021; pp. 1–5. [CrossRef]

25. Badr, G.; Oommen, B. On optimizing syntactic pattern recognition using tries and AI-based heuristic-search strategies. *IEEE Trans. Syst. Man, Cybern. Part B Cybernetics* **2006**, *36*, 611–622. [CrossRef] [PubMed]

26. Yang, J.; Zheng, X. Partial Discharge Pattern Recognition in GIS Based on Multiscale Dispersion Entropy and Stacking Ensemble Learning. In *Proceedings of the 2023 International Conference on Pattern Recognition, Machine Vision and Intelligent Algorithms (PRMVIA)*, Beihai, China, 24–26 March 2023; pp. 81–85. [CrossRef]

27. Sasmi Hidayatul, Y.T.; Djunaidy, A.; Muklason, A. Solving Multi-objective Vehicle Routing Problem Using Hyper-heuristic Method By Considering Balance of Route Distances. In *Proceedings of the 2019 International Conference on Information and Communications Technology (ICOIACT)*, Yogyakarta, Indonesia, 24–25 July 2019; pp. 937–942. [CrossRef]

28. Tang, Z.; Cao, Z. Application of Convolutional Neural Network Transfer Learning in Partial Discharge Pattern Recognition. In *Proceedings of the 2020 IEEE International Conference on High Voltage Engineering and Application (ICHVE)*, Beijing, China, 6–10 September 2020; pp. 1–4. [CrossRef]

29. Neteler, M.; Bowman, M.H.; Landa, M.; Metz, M. GRASS GIS: A multi-purpose open source GIS. *Environ. Model. Softw.* **2012**, *31*, 124–130. [CrossRef]

30. Team, G.D.; Landa, M.; Neteler, M.; Metz, M.; Petrášová, A.; Petráš, V.; Clements, G.; Zigo, T.; Larsson, N.; Kladivová, L.; et al. GRASS GIS, 2025. Available online: https://zenodo.org/records/14918754 (accessed on 12 June 2025).

31. McCabe, R. Policing the Seas: Building Constabulary Maritime Governance in the Horn of Africa–The Case of Djibouti and Kenya. *Afr. Secur.* **2019**, *12*, 330–355. [CrossRef]

32. Lawale, S.; Ahmad, T. The Role of Small States in Power Contestations in the Horn of Africa Case Study of Djibouti. *Asian J. Middle East. Islam. Stud.* **2023**, *17*, 325–340. [CrossRef]

33. Verdugo, D. Djibouti. In *Encyclopedia of Tourism*; Springer International Publishing: Cham, Switzerland, 2016; pp. 266–267. [CrossRef]
34. Haile, G.G.; Tang, Q.; Hosseini-Moghari, S.M.; Liu, X.; Gebremicael, T.G.; Leng, G.; Kebede, A.; Xu, X.; Yun, X. Projected Impacts of Climate Change on Drought Patterns Over East Africa. *Earth's Future* **2020**, *8*, e2020EF001502. [CrossRef]

35. Gebremeskel Haile, G.; Tang, Q.; Leng, G.; Jia, G.; Wang, J.; Cai, D.; Sun, S.; Baniya, B.; Zhang, Q. Long-term spatiotemporal variation of drought patterns over the Greater Horn of Africa. *Sci. Total. Environ.* **2020**, *704*, 135299. [CrossRef] [PubMed]

36. Lemenkova, P. Evapotranspiration, vapour pressure and climatic water deficit in Ethiopia mapped using GMT and TerraClimate dataset. *J. Water Land Dev.* **2022**, *54*, 201–209. [CrossRef]

37. AghaKouchak, A. A multivariate approach for persistence-based drought prediction: Application to the 2010–2011 East Africa drought. *J. Hydrol.* **2015**, *526*, 127–135. [CrossRef]

38. Agutu, N.; Awange, J.; Zerihun, A.; Ndehedehe, C.; Kuhn, M.; Fukuda, Y. Assessing multi-satellite remote sensing, reanalysis, and land surface models' products in characterizing agricultural drought in East Africa. *Remote Sens. Environ.* **2017**, *194*, 287–302. [CrossRef]

39. D'Ercole, R.; Casella, D.; Panegrossi, G.; Sanò, P. A high temporal resolution NDVI time series to monitor drought events in the Horn of Africa. *Int. J. Appl. Earth Obs. Geoinf.* **2024**, *135*, 104264. [CrossRef]

40. Tufa Dinku, P.C.; Connor, S.J. Challenges of satellite rainfall estimation over mountainous and arid parts of east Africa. *Int. J. Remote Sens.* **2011**, *32*, 5965–5979. [CrossRef]

41. Mishra, A.K.; Singh, V.P. A Review of Drought Concepts. *J. Hydrol.* **2010**, *391*, 202–216. [CrossRef]

42. Fava, F.; Vrieling, A. Earth observation for drought risk financing in pastoral systems of sub-Saharan Africa. *Curr. Opin. Environ. Sustain.* **2021**, *48*, 44–52. [CrossRef]

43. Bravar, L.; Kavvas, M. On the Physics of Droughts. I. A Conceptual Framework. *J. Hydrol.* **1991**, *129*, 281–297. [CrossRef]

44. Schulman, S. Climate Change Challenges and Djibouti. *RUSI J.* **2019**, *164*, 62–75. [CrossRef]

45. Carrão, H.; Naumann, G.; Barbosa, P. Mapping Global Patterns of Drought Risk: An Empirical Framework Based on Sub-National Estimates of Hazard, Exposure and Vulnerability. *Glob. Environ. Change* **2016**, *39*, 108–124. [CrossRef]

46. Fenta, A.A.; Tsunekawa, A.; Haregeweyn, N.; Poesen, J.; Tsubo, M.; Borrelli, P.; Panagos, P.; Vanmaercke, M.; Broeckx, J.; Yasuda, H.; et al. Land susceptibility to water and wind erosion risks in the East Africa region. *Sci. Total. Environ.* **2020**, *703*, 135016. [CrossRef] [PubMed]

47. Knight, J.; Abd Elbasit, M.A.M.; Adam, E. Land Degradation in Eritrea and Djibouti. In *Landscapes and Landforms of the Horn of Africa: Eritrea, Djibouti, Somalia*; Springer International Publishing: Cham, Switzerland, 2022; pp. 223–232. [CrossRef]

48. Boldrocchi, G.; Moussa Omar, Y.; Rowat, D.; Bettinetti, R. First results on zooplankton community composition and contamination by some persistent organic pollutants in the Gulf of Tadjoura (Djibouti). *Sci. Total. Environ.* **2018**, *627*, 812–821. [CrossRef] [PubMed]

49. Gajdzik, L.; Green, A.L.; Cochran, J.E.; Hardenstine, R.S.; Tanabe, L.K.; Berumen, M.L. Using species connectivity to achieve coordinated large-scale marine conservation efforts in the Red Sea. *Mar. Pollut. Bull.* **2021**, *166*, 112244. [CrossRef] [PubMed]

50. Billi, P. Climate Variability in the Horn of Africa Eastern Countries: Eritrea, Djibouti, Somalia. In *Landscapes and Landforms of the Horn of Africa: Eritrea, Djibouti, Somalia*; Springer International Publishing: Cham, Switzerland, 2022; pp. 1–39. [CrossRef]

51. Seka, A.M.; Guo, H.; Zhang, J.; Han, J.; Bayable, E.; Ayele, G.T.; Workneh, H.T.; Bayouli, O.T.; Muhirwa, F.; Reda, K.W. Evaluating the future total water storage change and hydrological drought under climate change over lake basins, East Africa. *J. Clean. Prod.* **2024**, *447*, 141552. [CrossRef]

52. Tesfai, R.T.; Parrini, F.; Owen-Smith, N.; Moehlman, P.D. African wild ass drinking behaviour on the Messir Plateau, Danakil Desert, Eritrea. *J. Arid. Environ.* **2021**, *185*, 104327. [CrossRef]

53. Medin, T.; Martínez-Navarro, B.; Rivals, F.; Libsekal, Y.; Rook, L. The late Early Pleistocene suid remains from the paleoanthropological site of Buia (Eritrea): Systematics, biochronology and eco-geographical context. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **2015**, *431*, 26–42. [CrossRef]

54. Ali, M.M.; Abd Ellah, R.G. Chapter 1-History and formation of African Lakes. In *Lakes of Africa*; El-Sheekh, M., Elsaied, H.E., Eds.; Elsevier: Amsterdam, The Netherlands, 2023; pp. 1–31. [CrossRef]

55. Idowu, T.E.; Lasisi, K.H. Seawater intrusion in the coastal aquifers of East and Horn of Africa: A review from a regional perspective. *Sci. Afr.* **2020**, *8*, e00402. [CrossRef]

56. Blanco-Sacristán, J.; Johansen, K.; Duarte, C.M.; Daffonchio, D.; Hoteit, I.; McCabe, M.F. Mangrove distribution and afforestation potential in the Red Sea. *Sci. Total. Environ.* **2022**, *843*, 157098. [CrossRef] [PubMed]

57. Beyin, A. Living by the land, gazing at the sea: Hominin occupation of near-coastal landscapes on the western periphery of the Red Sea. *Quat. Sci. Adv.* **2024**, *16*, 100247. [CrossRef]

58. Klaus, R. Coral Reefs and Communities of the Central and Southern Red Sea (Sudan, Eritrea, Djibouti, and Yemen). In *The Red Sea: The Formation, Morphology, Oceanography and Environment of a Young Ocean Basin*; Springer: Berlin/Heidelberg, Germany, 2015; pp. 409–451. [CrossRef]

59. Cowburn, B.; Samoilys, M.; Osuka, K.; Klaus, R.; Newman, C.; Gudka, M.; Obura, D. Healthy and Diverse Coral Reefs in Djibouti–A Resilient Reef System or Few Anthropogenic Threats? *Mar. Pollut. Bull.* **2019**, *148*, 182–193. [CrossRef] [PubMed]
60. Gapper, J.J.; Maharjan, S.; Li, W.; Linstead, E.; Tiwari, S.P.; Qurban, M.A.; El-Askary, H. A generalized machine learning model for long-term coral reef monitoring in the Red Sea. *Heliyon* **2024**, *10*, e38249. [CrossRef] [PubMed]

61. Chandrasekharam, D.; Lashin, A.; Al Arifi, N.; Al Bassam, A.A.; Varun, C. Evolution of geothermal systems around the Red Sea. *Environ. Earth Sci.* **2015**, *73*, 4215–4236. [CrossRef]

62. Moussa, N.; Bayon, G.; Dekov, V.; Yamanaka, T.; Shinjo, R.; Toki, T.; Le Gall, B.; Grassineau, N.; Langlade, J.; Awaleh, M.; et al. Mixed carbonate-siliceous hydrothermal chimneys ahead of the Asal propagating rift (SE Afar Rift, Republic of Djibouti). *J. Afr. Earth Sci.* **2023**, *197*, 104765. [CrossRef]

63. Ayele, A.; Luckett, R.; Baptie, B.; Whaler, K. The 2015 earthquake swarm in the Fentale volcanic complex (FVC): A geohazard risk for Ethiopia's commercial route to the Djibouti port. *J. Afr. Earth Sci.* **2024**, *213*, 105236. [CrossRef]

64. Galludin, M.; Razack, M. Analysis of pumping tests, with regard to tectonics, hydrothermal effects and weathering, for fractured Dalha and stratiform basalts, Republic of Djibouti. *J. Hydrol.* **1994**, *155*, 237–250. [CrossRef]

65. Rime, V.; Foubert, A.; Ruch, J.; Kidane, T. Tectonostratigraphic evolution and significance of the Afar Depression. *Earth-Sci. Rev.* **2023**, *244*, 104519. [CrossRef]

66. Gidafie, D.; Nedaw, D.; Azagegn, T.; Abebe, B.; Baba, A. Evaluation of the source and mechanisms of groundwater recharge for the southern sections of the western Afar rift margin and associated rift floor. *J. Hydrol. Reg. Stud.* **2024**, *56*, 102037. [CrossRef]

67. Darrah, T.H.; Tedesco, D.; Tassi, F.; Vaselli, O.; Cuoco, E.; Poreda, R.J. Gas chemistry of the Dallol region of the Danakil Depression in the Afar region of the northern-most East African Rift. *Chem. Geol.* **2013**, *339*, 16–29. [CrossRef]

68. Varet, J. Geomorphology of Afar. In *Landscapes and Landforms of the Horn of Africa: Eritrea, Djibouti, Somalia*; Springer International Publishing: Cham, Switzerland, 2022; pp. 81–124. [CrossRef]

69. Lemenkova, P. Tanzania Craton, Serengeti Plain and Eastern Rift Valley: Mapping of geospatial data by scripting techniques. *Est. J. Earth Sci.* **2022**, *71*, 61–79. [CrossRef]

70. Pinzuti, P.; Mignan, A.; King, G.C. Surface morphology of active normal faults in hard rock: Implications for the mechanics of the Asal Rift, Djibouti. *Earth Planet. Sci. Lett.* **2010**, *299*, 169–179. [CrossRef]

71. Lemenkova, P. Seismicity in the Afar Depression and Great Rift Valley, Ethiopia. *Environ. Res. Eng. Manag.* **2022**, *78*, 83–96. [CrossRef]

72. Meshesha, D.; Chekol, T.; Negussia, S. Major and trace element compositions of basaltic lavas from western margin of central main Ethiopian rift: Enriched asthenosphere vs. mantle plume contribution. *Heliyon* **2021**, *7*, e08634. [CrossRef] [PubMed]

73. Mologni, C.; Bruxelles, L.; Schuster, M.; Davtian, G.; Ménard, C.; Orange, F.; Doubre, C.; Cauliez, J.; Taezaz, H.B.; Revel, M.; et al. Holocene East African monsoonal variations recorded in wave-dominated clastic paleoshorelines of Lake Abhe, Central Afar region (Ethiopia & Djibouti). *Geomorphology* **2021**, *391*, 107896. [CrossRef]

74. D'Amore, F.; Giusti, D.; Abdallah, A. Geochemistry of the High-Salinity Geothermal field of Asal, Republic of Djibouti, Africa. *Geothermics* **1998**, *27*, 197–210. [CrossRef]

75. Houssein, D.E.; Axelsson, G. Geothermal resources in the Asal Region, Republic of Djibouti: An update with emphasis on reservoir engineering studies. *Geothermics* **2010**, *39*, 220–227. [CrossRef]

76. Awaleh, M.B.; Nishijima, J. Delineation of Geological Structures of Arta Geothermal Prospect in Djibouti Based on the Gravity Data Analysis and Interpretation. *Geothermics* **2024**, *117*, 102894. [CrossRef]

77. Bouh Houssein, D.; Chandrasekharam, V.C.; Galludin, M. Geochemistry of Thermal Springs Around Lake Abhe, Western Djibouti. *Int. J. Sustain. Energy* **2014**, *33*, 1090–1102. [CrossRef]

78. Castanier, S.; Perthuisot, J.P.; Rouchy, J.M.; Maurin, A.; Guelorget, O. Halite Ooids in Lake Asal, Djibouti: Biocrystalline Build-UPS. *Geobios* **1992**, *25*, 811–821. [CrossRef]

79. Gasse, F.; Fontes, J.C. Palaeoenvironments and Palaeohydrology of a Tropical Closed Lake (Lake Asal, Djibouti) Since 10,000 yr B.P. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **1989**, *69*, 67–102. [CrossRef]

80. Dekov, V.M.; Guéguen, B.; Yamanaka, T.; Moussa, N.; Okumura, T.; Bayon, G.; Liebetrau, V.; Yoshimura, T.; Kamenov, G.; Araoka, D.; et al. When a Mid-Ocean Ridge Encroaches a Continent: Seafloor-Type Hydrothermal Activity in Lake Asal (Afar Rift). *Chem. Geol.* **2021**, *568*, 120126. [CrossRef]

81. Zan, L.; Gianelli, G.; Passerini, P.; Troisi, C.; Haga, A.O. Geothermal Exploration in the Republic of Djibouti: Thermal and Geological Data of the Hanlé and Asal Areas. *Geothermics* **1990**, *19*, 561–582. [CrossRef]

82. Gong, Z.; Ge, W.; Guo, J.; Liu, J. Satellite remote sensing of vegetation phenology: Progress, challenges, and opportunities. *ISPRS J. Photogramm. Remote Sens.* **2024**, *217*, 149–164. [CrossRef]

83. Allu, A.R.; Mesapam, S. Impact of remote sensing data fusion on agriculture applications: A review. *Eur. J. Agron.* **2025**, *164*, 127478. [CrossRef]

84. Chanuwan Wijesinghe, D.; Chaminda Withanage, N.; Kumar Mishra, P.; Ranagalage, M.; Abdelrahman, K.; Fnais, M.S. An application of the remote sensing derived indices for drought monitoring in a dry zone district, in tropical island. *Ecol. Indic.* **2024**, *167*, 112681. [CrossRef]
85. Wang, Y.; Sun, Y.; Cao, X.; Wang, Y.; Zhang, W.; Cheng, X. A review of regional and Global scale Land Use/Land Cover (LULC) mapping products generated from satellite remote sensing. *ISPRS J. Photogramm. Remote Sens.* **2023**, *206*, 311–334. [CrossRef]

86. Darem, A.A.; Alhashmi, A.A.; Almadani, A.M.; Alanazi, A.K.; Sutantra, G.A. Development of a map for land use and land cover classification of the Northern Border Region using remote sensing and GIS. *Egypt. J. Remote Sens. Space Sci.* **2023**, *26*, 341–350. [CrossRef]

87. Farooq, B.; Manocha, A. Satellite-based change detection in multi-objective scenarios: A comprehensive review. *Remote Sens. Appl. Soc. Environ.* **2024**, *34*, 101168. [CrossRef]

88. Afuye, G.A.; Nduku, L.; Kalumba, A.M.; Santos, C.A.G.; Orimoloye, I.R.; Ojeh, V.N.; Thamaga, K.H.; Sibandze, P. Global trend assessment of land use and land cover changes: A systematic approach to future research development and planning. *J. King Saud Univ.-Sci.* **2024**, *36*, 103262. [CrossRef]

89. Chughtai, A.H.; Abbasi, H.; Karas, I.R. A review on change detection method and accuracy assessment for land use land cover. *Remote Sens. Appl. Soc. Environ.* **2021**, *22*, 100482. [CrossRef]

90. Velastegui-Montoya, A.; Escandón-Panchana, P.; Peña-Villacreses, G.; Herrera-Franco, G. Land use/land cover of petroleum activities in the framework of sustainable development. *Clean. Eng. Technol.* **2023**, *15*, 100659. [CrossRef]

91. Farhan, M.; Wu, T.; Amin, M.; Tariq, A.; Guluzade, R.; Alzahrani, H. Monitoring and prediction of the LULC change dynamics using time series remote sensing data with Google Earth Engine. *Phys. Chem. Earth Parts A/B/C* **2024**, *136*, 103689. [CrossRef]

92. Amazirh, A.; Chehbouni, A.; Bouras, E.; Benkirane, M.; Hssaine, B.; Entekhabi, D. Drought cascade lag time estimation across Africa based on remote sensing of hydrological cycle components. *Adv. Water Resour.* **2023**, *182*, 104586. [CrossRef]

93. Mekonnen, Y.G.; Alamirew, T.; Chukalla, A.D.; Malede, D.A.; Yalew, S.G.; Mengistu, A.F. Remote sensing in hydrology: A systematic review of its applications in the Upper Blue Nile Basin, Ethiopia. *HydroResearch* **2025**, *8*, 1–12. [CrossRef]

94. Parra-Paitan, C.; Meyfroidt, P.; Verburg, P.H.; zu Ermgassen, E.K. Deforestation and climate risk hotspots in the global cocoa value chain. *Environ. Sci. Policy* **2024**, *158*, 103796. [CrossRef]

95. Zhao, Z.; Ciais, P.; Wigneron, J.P.; Santoro, M.; Brandt, M.; Kleinschroth, F.; Lewis, S.L.; Chave, J.; Fensholt, R.; Laporte, N.; et al. Central African biomass carbon losses and gains during 2010–2019. *One Earth* **2024**, *7*, 506–519. [CrossRef]

96. Nzabarinda, V.; Bao, A.; Tie, L.; Uwamahoro, S.; Kayiranga, A.; Ochege, F.U.; Muhirwa, F.; Bao, J. Expanding forest carbon sinks to mitigate climate change in Africa. *Renew. Sustain. Energy Rev.* **2025**, *207*, 114849. [CrossRef]

97. Montero-Botey, M.; Kivuyo, E.; Sitati, N.; Perea, R. Deforestation and water availability as main drivers of human-elephant conflict. *Glob. Ecol. Conserv.* **2024**, *54*, e03068. [CrossRef]

98. Aksoy, S.; Yildirim, A.; Gorji, T.; Hamzehpour, N.; Tanik, A.; Sertel, E. Assessing the performance of machine learning algorithms for soil salinity mapping in Google Earth Engine platform using Sentinel-2A and Landsat-8 OLI data. *Adv. Space Res.* **2022**, *69*, 1072–1086. [CrossRef]

99. Nesari, H.; Shah-Hosseini, R.; Goodarzi, A.; Sobhanardakani, S.; Farzaneh, S. Integration of Landsat 8 (OLI) and MODIS images to monitor suspended particles and evaluate the spatial pattern of air pollution. *Atmos. Pollut. Res.* **2024**, *15*, 101989. [CrossRef]

100. Kumar, B.; Babu, K.; Anusha, B.; Rajasekhar, M. Geo-environmental monitoring and assessment of land degradation and desertification in the semi-arid regions using Landsat 8 OLI / TIRS, LST, and NDVI approach. *Environ. Chall.* **2022**, *8*, 100578. [CrossRef]

101. Wright, D.; Wood, R.; Sylvander, B. ArcGMT: A suite of tools for conversion between Arc/INFO® and Generic Mapping Tools (GMT). *Comput. Geosci.* **1998**, *24*, 737–744. [CrossRef]

102. Lemenkova, P. Mapping Climate Parameters over the Territory of Botswana Using GMT and Gridded Surface Data from TerraClimate. *ISPRS Int. J. Geo-Inf.* **2022**, *11*, 473. [CrossRef]

103. Köcher, A.; Markaj, A.; Fay, A. Toward a Generic Mapping Language for Transformations between RDF and Data Interchange Formats. In Proceedings of the 2022 IEEE 27th International Conference on Emerging Technologies and Factory Automation (ETFA), Stuttgart, Germany, 6–9 September 2022; pp. 1–4. [CrossRef]

104. Lemenkova, P. Console-Based Mapping of Mongolia Using GMT Cartographic Scripting Toolset for Processing TerraClimate Data. *Geosciences* **2022**, *12*, 140. [CrossRef]

105. Bawa, S.; Onotu, A.A.; Akomolafe, E.A.; Sa'i, U.I. Understanding the Topography of the Gulf of Guinea Seabed Using GMT Scripting and GEBCO Gridded Data. In Proceedings of the 2024 IEEE Mediterranean and Middle-East Geoscience and Remote Sensing Symposium (M2GARSS), Oran, Algeria, 15–17 April 2024; pp. 516–520. [CrossRef]

106. Tangkitjaroenmongkol, R.; Kaittisin, S.; Ongwattanakul, S. Inbound logistics cassava starch planning: With application of GIS and K-means clustering methods in Thailand. In Proceedings of the 2011 Eighth International Joint Conference on Computer Science and Software Engineering (JCSSE), Nakhonpathom, Thailand, 11–13 May 2011; pp. 204–209. [CrossRef]

107. Xing, F.; Zhao, Y.; Song, Z.; Liu, Z.; Chai, W.; Zhang, Y.; Li, M.; Jia, T.; Zhao, J.; Xu, M. GIS Gas Analysis Training Program Based on K-Means Clustering Method Training Auxiliary Analysis. In Proceedings of the 2023 International Conference on Educational Knowledge and Informatization (EKI), Guangzhou, China, 22–24 September 2023; pp. 78–82. [CrossRef]
108. Lemenkova, P. Mapping Wetlands of Kenya Using Geographic Resources Analysis Support System (GRASS GIS) with Remote Sensing Data. *Transylv. Rev. Syst. Ecol. Res.* **2023**, *25*, 1–18. [CrossRef]

109. Huo, X.J.; Moon, K.S.; Lee, S.H.; Seung, T.Y.; Kwon, K.R. Protecting GIS vector map using the k-means clustering algorithm and odd-even coding. In *Proceedings of the 2011 17th Korea-Japan Joint Workshop on Frontiers of Computer Vision (FCV), Ulsan, Republic of Korea*, 9–11 February 2011; pp. 1–5. [CrossRef]

110. Chowdhury, M.S. Comparison of accuracy and reliability of random forest, support vector machine, artificial neural network and maximum likelihood method in land use/cover classification of urban setting. *Environ. Chall.* **2024**, *14*, 100800. [CrossRef]

111. Lemenkova, P. Using Open-Source Software GRASS GIS for Analysis of the Environmental Patterns in Lake Chad, Central Africa. *Die Bodenkultur J. Land Manag. Food Environ.* **2023**, *74*, 49–64. [CrossRef]

112. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Müller, A.; Nothman, J.; Louppe, G.; et al. Scikit-learn: Machine Learning in Python. *arXiv* **2018**, arXiv:1201.0490. [CrossRef]

113. Wenting, X.; Bingfang, W.; Yichen, T.; Yuan, Z. Mapping plant diversity of broad-leaved forest ecosystem using Landsat TM data. In *Proceedings of the IGARSS 2004–2004 IEEE International Geoscience and Remote Sensing Symposium, Anchorage, AK, USA*, 20–24 September 2004; Volume 7, pp. 4598–4600. [CrossRef]

114. Dobrowski, S.Z.; Safford, H.D.; Cheng, Y.B.; Ustin, S.L. Mapping Mountain Vegetation using Species Distribution Modeling, Image-Based Texture Analysis, and Object-Based Classification. *Appl. Veg. Sci.* **2008**, *11*, 499–508. [CrossRef]

115. Moore, C.A.; Hoffmann, G.A.; Glenn, N.F. Quantifying Basalt Rock Outcrops in NRCS Soil Map Units Using Landsat-5 Data. *Soil Surv. Horizons* **2007**, *48*, 59–62. [CrossRef]

116. Masemola, C.; Cho, M.A.; Ramoelo, A. Assessing the Effect of Seasonality on Leaf and Canopy Spectra for the Discrimination of an Alien Tree Species, Acacia Mearnsii, From Co-Occurring Native Species Using Parametric and Nonparametric Classifiers. *IEEE Trans. Geosci. Remote Sens.* **2019**, *57*, 5853–5867. [CrossRef]

117. Lemenkova, P. Monitoring Seasonal Fluctuations in Saline Lakes of Tunisia Using Earth Observation Data Processed by GRASS GIS. *Land* **2023**, *12*, 1995. [CrossRef]

118. Lemenkova, P. A GRASS GIS Scripting Framework for Monitoring Changes in the Ephemeral Salt Lakes of Chotts Melrhir and Merouane, Algeria. *Appl. Syst. Innov.* **2023**, *6*, 61. [CrossRef]

119. Eldosouky, A.M.; Pham, L.T.; El-Qassas, R.A.Y.; Hamimi, Z.; Oksum, E. Lithospheric Structure of the Arabian–Nubian Shield Using Satellite Potential Field Data. In *The Geology of the Arabian-Nubian Shield*; Springer International Publishing: Cham, Switzerland, 2021; pp. 139–151. [CrossRef]

**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.