Past and Future Land Use/Land Cover Changes in the Ethiopian Fincha Sub-Basin

: The increasing human pressure on African regions is recognizable when looking at Land Use Land Cover (LULC) change maps, generally derived from satellite imagery. Using the Ethiopian Fincha watershed as a case study, the present work focuses on (i) identifying historical LULC change in the period 1989–2019; (ii) estimating LULC in the next thirty years, combining Geographical Information Systems (GIS) with Land Change Modelling (LCM). Landsat 5/8 images were combined with ﬁeld evidence to map LULC in three reference years (1989, 2004, 2019), while the Multi-Layer Markov Chain (MPL-MC) model of LCM was applied to forecast LULC in 2030, 2040, and 2050. The watershed was classiﬁed into six classes: waterbody, grass/swamp, built-up, agriculture; forest; and shrub. The results have shown that, in the past 30 years, the Fincha watershed experienced a reduction in forest and shrubs of about − 40% and − 13%, respectively, mainly due to ever-increasing agricultural activities, and such a trend is also expected in the future. In fact, for the period 2019–2050, LCM simulated a signiﬁcant decrease in both forest and shrubs (around − 70% and − 20%, respectively), in favor of more areas covered by grass (19%) and built-up (20%). It is worth noting that a decrease in natural forests can drive an increase in soil erosion, fostering siltation in the water reservoirs located in the sub-basin. The study pointed out the urgency of taking actions in the sub-basin to counteract such changes, which can eventually lead to a less sustainable environment.


Introduction
Land use is defined as how land is utilized by human beings and their habitats, usually with an accent on the practical role of land for economic activities, whereas land cover is a physical characteristic of the Earth's surface or attributes of a part of the Earth's land surface and immediate subsurface, including biota, soil, topography, surface and groundwater, and human structures [1][2][3][4]. As it is strictly connected with representing the hydrological cycle [5,6], land use and land cover (LULC) change has been one of the most widely used methods to understand past land uses, types of changes estimated, the forces behind such changes, and the perceptible transformations of the Earth's surface [1]. LULC changes could involve critical issues such as biodiversity degradation and negative impact on human life [7,8]. The study of LULC change has attracted growing interest in recent years and is a complex issue that involves physical, environmental, and socioeconomic facts. According to Lambin et al. [9], the modeling of land cover processes can answer questions such as (i) which are the main environmental and cultural variables that contribute most to the observed changes, and why? (ii) within a geographical region, which locations can be affected by land cover changes, where and at what rate does land cover change, and when?
Prediction using time series data is important for the future management plan of LULC, and it is frequently employed as a diverse appropriateness measure as a proxy of human influence on land change processes [4,10]. Analysis of the historical trends of LULC is paramount in modeling future LULC, as past information generally represents a good proxy of human influence on land processes [10,11]. To adequately predict future scenarios, Four main seasons characterize the region: Summer, from June to August, with heavy rainfalls; Autumn from September to November (called harvest season); Winter, from December to February (the dry season with frost in the morning, especially in January); Spring, from March to May (occasional showers and the hottest season). The annual rainfall in the study area ranges between 1367 and 1842 mm with the minimum rainfall occurring in the northern lowlands and maximum rainfall greater than 1500 mm in the southern and western highlands. June to September is the main rainy season of the catchment, with an average of 1604 mm and maximum rainfall between July and August.
Natural resources, such as the Fincha, Amarti, and Nashee lakes (see Figure 1c), contribute to the national economy by generating hydroelectric power but are also used for irrigating large fields devoted to sugar cane. The area is of interest for national and international hydro-politics due to its downstream connection to the Nile basin and its intense agriculture.

Dataset
The study was performed using freely available satellite imagery and a Digital Elevation Model (DEM). The latter, having a resolution of 30 m and referring to 2019, was acquired from the GIS and Remote Sensing Department, Ministry of Water, Irrigation and Energy of Ethiopia [42].
Landsat-5 TM (L5, for the years 1989 and 2004) and Landsat-8 OLI-TIRS (L8, for the year 2019) data were downloaded from the United States Geological Survey (USGS) website (earthexplorer.usgs.gov). As the Landsat-5 mission started in March 1984, it was not possible to acquire the images every 10 years, so the first reference year was set to 1989. The images referred to January, when there is a clear sky corresponding to the dry season, and were atmospherically corrected via QGIS (qgis.org). To cover the whole watershed area, a composite of Landsat images from different paths/rows was created, ensuring that the images refer to the same season (Table 1). Field surveys have been conducted to assist the LULC classification of the satellite images. In addition, key informant interviews (KII) and focal group discussions were performed to obtain socio-economic support data, as this is paramount to understanding how locals interact with the environment [39][40][41][42][43][44][45]. KII were conducted with elders, as they have known the area for at least 30 years and have good knowledge of past LULC changes. Focal Group Discussions (FGD) were conducted with experts from zonal and district offices of agriculture; natural resources management; environment and climatic change; land use administration, and with local people's representatives. Ground truth data were collected using GPS and digital cameras to evaluate the current LULC.
Open-ended questions about LULC's significant changes in the study area, the connection between the biophysical environment, institution, socioeconomic activity, and demography were utilized during both KII and FGD. To learn more from a management point of view, assess the efforts made towards resource management, and identify obstacles, discussions on the practices and regulations that affect land management and policies in the area were held. The topic of land degradation and the most urgent problems that need solutions were also covered. The major goals of the discussion and interviews were to gather enough information on the historical and present trends of LULC changes, identify their fundamental causes, and assess their effects on regional socio-economic life and the environment. In detail, farmers were asked to describe the areas of the landscape that have altered and to recognize the reasons behind those changes. Moreover, they were questioned on the effects of the modifications to their way of life, their surroundings, and their working environment. In addition, farmers were asked to describe how their socioeconomic activity affects the change in land usage.
Based on a checklist created in advance to assess the situation in the watershed, field observations were conducted, and images of significant sites were obtained.

Land Use Land Cover Classification of Historical Data
To map LULC, satellite images should be classified, assigning predefined LULC classes to some pixels. As pointed out by Jemberie et al. [44], this phase could be affected by various factors such as classification methods, algorithms, collecting of training sites, and the quality (correctness) of the classification should be assessed via field evidence [46,47].
The study was performed by classifying three reference years (1989,2004, and 2019) and considering six LULC classes. The selection of these classes was performed based on past studies [7], field evidence, and information coming from local farmers and experts, as well as the personal considerations of the authors. Based on this, the Fincha sub-basin was classified into six classes, namely waterbody, built-up (urban and rural settlement), agriculture, forest (dense forest and sparse/desert forest), grass, and swamps ( Table 2). In the process of classification, it is difficult to differentiate some LULC classes' spectral properties from other classes. For instance, sugar cane, currently planted for the Fincha sugar factory, and grass; urban and rural settlements, roads, industries, and infrastructures have the same properties. This happens because the spectral properties of some LULC classes appear identical to others. To simplify the complexity and reduce the number of LULC classes, some related LULC classes were merged to form one class. For example, urban and rural settlements, roads, industries, and infrastructures were aggregated as built-up, while grass, swamps, and land that is covered by sugar cane, were commonly named grass/swamp.
Following ample literature on this topic e.g., [46][47][48][49][50], the maximum likelihood supervised classification method was applied via ArcGIS by creating training sites. For the L8 image of 2019, such training sites were defined using 100 ground truth points, while, for the two older L5 images, training signature sites were defined via unsupervised classification, ancillary data (Google Earth and Copernicus data), and KII information and literature data [7]. To improve image quality, quality assessments were used by taking a total of 50 ground truth points (20 agriculture, 5 waterbody, 5 built-up, 10 forest, 5 shrub, and 5 grass/swamps). The points were uniformly distributed across the watershed to guarantee a proper classification.
To quantitatively assess the accuracy, statistical methods such as overall accuracy and kappa value were applied. Based on this, random sampling data were prepared to check the overall accuracy OA and to determine the kappa coefficient K. Comparing the total corrected samples TCS and the total samples TS, OA provides an idea of how many sites are correctly classified (Equation (1)), and ranges from 0 (corrected samples) to 1 (very accurate classification).
The kappa coefficient K (Equation (2)) is generated from a statistical test and describes the accuracy of a classification compared to a random classification [51,52]. Its value varies between 0 and 1, where 0 indicates a total accidental classification, while 1 indicates a very accurate classification. According to Gidey et al. [53], good classifications have K > 0.8, while bad classifications have K < 0.4.
where the matrix columns indicate the correspondence between ground truth data and the pixel location, while the matrix rows indicate to which class the is pixel assigned.

Prediction of Future LULC and Associated Driving Forces
To manage natural resources (biodiversity) influences, and to analyze and forecast spatial LULC changes, the Land Change Modeler (LCM) in TerrSet (formerly known as IDRISI) software was developed [54,55]. LCM is an ArcGIS-integrated suite of tools for the assessment of future LULC changes, detecting gains and losses, net change, persistence, and identification of transitions between LULC classes [56]. To map future LULC scenarios, LCM utilizes historical LULC maps and a series of driving forces ( Table 3). The Markov chain projection is performed by creating matrixes to estimate the transition probability and the area of each LULC class for future dates [57,58]. In this study, LCM was applied to forecast the future LULC in three scenarios (2030, 2040, 2050), via a few main steps: (i) analysis of historical LULC maps (1989,2004, and 2019) and associated changes, (ii) creation of transition probability matrixes, (iii) model validation, iv) prediction of future LULC maps, accounting for possible driving forces. In this work, we define the probability transitional matrix as a matrix showing the transfer direction of LULC types from one category to other categories in a given year [10].
In LCM, there are two options for modeling algorithms that are used to model the transition variables: logistic regression and Multi-Layer Perceptron (MLP) neural network [59,60]. MLP uses minimal parameters, is more easily approachable, and has been extensively enhanced to offer an automatic mode that requires no user intervention. Therefore, in the present study, the MLP neural network has been employed.
To evaluate the capability of LCM in predicting future LULC, a predicted map of 2019 was created based on 1989 and 2004 LULC, and then compared with the actual 2019 map. To evaluate the quality of the 2019 predicted map against the 2019 reference map, the TerrSet validation module was used [61], mimicking the approach proposed in similar studies [62]. In TerrSet, two tools are available to assess the fit of the model to the sample data. First, the cross-validation tool iteratively removes a sample data point and interpolates a new value for the location. A table is produced to show the difference between the predicted attributes and the known attributes at those locations. Second, a variance image is produced that shows the spatial variation of uncertainty as a result of the fitted model. The variance image provides information to assist in identifying the problem areas where the relationship between the fitted model and the sample data points is poor [61].
Kappa indices, such as kappa for no information (Kno), kappa for location (Klocation), and kappa standards (Kstandards)l are used to identify potential errors [2,63,64]. Kappa values vary between 0 and 1, with values >0.8 meaning an almost perfect agreement. In detail, Kstandards is an index of agreement that attempts to account for the expected agreement due to random spatial reallocation of the categories in the comparison map; Kno is identical to Kstandards if both the quantity and allocation of categories in the comparison map are selected randomly; Klocation represents the extent to which the maps agree in terms of location of each LULC category.
To corroborate the study outcomes, a series of statistics were considered [65]: agreement due to chance (agreement chance), agreement due to quantity (agreement quantity), agreement due to the location at the grid cell level (agreement grid cell), disagreement due to the location at the grid cell level (disagree grid cell), and disagreement due to quantity (disagree quantity) were calculated to indicate how well the comparison map agrees with the reference map [40].
Driving forces are the factors that affect LULC changes at the local scale, and therefore they should be locally investigated and addressed [8,[66][67][68]. The driving variables (Table 3) were selected based on the actual literature and past studies, selecting the most important. In fact, there are still some other factors that are difficult to quantify, such as the population in the area. In simulating future LULC, LCM differentiates between static and dynamic variables, where the first are stable in time while the latter change temporally, and are therefore recalculated at each time step.
The type of land cover is strongly correlated with anthropogenic disturbance, for example, the local population can access resources more conveniently while changing the land use because of the distance from the stream. The ease with which land can transition to urban usage depends on the distance from urban centers, which can be a highly powerful force for change. One of the key factors in drawing more urban uses and expansion is the distance from roads, which determines accessibility. The primary topographic component known to affect LULC change is elevation. In addition, it seems reasonable to use the evidence likelihood, a quantitative variable that reflects the likelihood of discovering change between agricultural land and all other land classes at the relevant pixel, given that the annual pace of agricultural expansion was considerable. The watershed slope influences the spatial trends of land cover change, leading one to assume that changes in land use are highly influenced by the terrain slope: gains in urban land are primarily concentrated on relatively flat slopes and deforestation declines as the slope's gradient increases.
It is worth remembering that the selection of variables and indicators, to a certain extent, may cause some differences in the simulation results or model parameters, which will have effects on the prediction of LULC change. For example, for distance from the road, if the forest is very close to the road, the rate of forest clearance (deforestation) is very high, and vice versa, i.e., if there is road availability, the people living nearby can clear the forest for agriculture or other purposes. This is also applicable to other driving factors.
The Cramer's V Coefficient (CVC), sometimes called Cramer's V strategy [66], was used to assess the correlations between the various driving variables. In statistics, CVC is a measure of association between two categorical variables, giving a value between 0 and +1, and it is based on Pearson's chi-squared test [67]. According to Eastman [68], variables that have a Cramer's V > 0.40 are good and these drivers will have the greatest impact on the modification process and its spatial distribution [69][70][71]. One has to remember that CVC does not recognize interaction effects between the explanatory variables and land cover categories, while it only helps to determine whether to include a specific variable as a driving factor of LULC changes.

LULC Detection
LULC changes were detected via a few parameters: magnitude of change C, rate of change R, and change percentage P, using the following equations [2,72,73]: where i represents the LULC class, B i and L i are the areas [ha] with the earliest and latter LULC, respectively. The period between B i and L i is T [year] and determines the rate of change R i . Positive values of P i mean an increase in a specific LULC in the study period T

Historical LULC Maps
Three reference years (1989,2004, and 2019) were considered to evaluate historical LULC via a maximum likelihood supervised classification ( Figure 2). As reported in Table 4, in 1989 most of the study area was covered by agriculture (32%), grass/swamps (24%), and shrub (22%), with only a very minor part occupied by settlements (0.4%). Similar LULC was observed in 2004, with agriculture (34%), grass/swamp (24%), and shrub (18%) being the most dominant LULC classes and just a small increase in the area covered by built-up (1%). In 2019, the class distribution remained similar, with an increase in the built-up area (1.7%). In summary, in the past, agriculture was always the most dominant LULC class in the Fincha watershed, followed by grass/swamp and shrub.  The results reported in Figure 2 are in agreement with Dibaba et al. [39], who pointed out that the Fincha watershed is characterized by an expansion of agriculture and built-up LULC, resulting in a decline of natural vegetation.

Accuracy Assessment for Historical LULC
The overall accuracies OA and kappa values K were 82.80%, 85.57%, and 89.82% and 80.51%, 82.54%, and 87.84%, respectively, for the three reference years (Tables 5 and A1,  Tables A2 and A3). These results indicate that the accuracy of the classifications improved from 1989 to 2019, also thanks to the higher quality of the satellite data used. The accuracy of a map could be different for users and map developers. The user's accuracy indicates how often a specified class on the map is present on the ground, while the producer's (mapmaker) accuracy shows the probability that a certain land cover is classified according to field evidence.
Hailu et al. [40] defined the kappa statistics <40%, 40-75%, and >75% as poor, good, and excellent, respectively. Using this approach, from Table 5 one can notice that the statistics of the Fincha watershed were excellent, meaning a very good agreement between the classification maps and the reference information.

Historical LULC Changes and Transition Probability Matrix
Comparing the three reference years, it is possible to observe a considerable reduction in the area covered by forest and shrubs during the observation period (Table 6). In detail, yearly, around 639 ha of forest and 280 ha of shrubs were cleared in favor of other LULC classes. As anticipated, human pressure contributed to changing the environment, as recognizable by the increase in areas covered by agricultural fields, built-up, grass/swamp, and waterbodies, which yearly gained around 616 ha, 125 ha, 7 ha, and 171 ha, respectively. Waterbodies increased significantly during the last 30 years, mainly because of human intervention. In fact, in 1989, the Amerti reservoir, one of the reservoirs located in the Fincha watershed, was not filled, while it was filled in 2004. In 2019, another dam was constructed over the Nashe River [2]. The study pointed out small changes in terms of grass/swamps, at least in terms of net variation. In fact, as is visible from Figure 3 and Table 6, the majority of the Fincha watershed was affected by variations in LULC that include this class.
The results presented here are in line with the existing literature on LULC in the Fincha watershed [17,39]. All the authors agreed that the shift from natural LULC towards more anthropized environments could threaten biodiversity and decrease the total values of ecosystem services [74].  It is worth recalling that the probability transitional matrix is the transfer direction of LULC types from one category to other categories in the given year [10]. The nature change can be distinguished from the Markov transitional matrices for historical LULC over the period (1989-2004 and 2004-2019). The nature of change can be distinguished from the trend as depicted in the Markov transition matrices over the period between 1989 and 2019. The diagonal values represent the probability that each land cover class remains persistent (constant) from earlier to later years. The other values represent a given land cover land class undertaking transition to another land cover land class.
Between 1989 and 2004, the highest and the lowest persistent LULC classes were waterbody and grass/swamps, characterized by a percentage of stability of 91% and 42%, respectively. During the period 2004-2019, the most and the least stable LULC class categories were waterbody and grass/swamp, which accounted for around 80% and 31%, respectively. Over the entire temporal horizon observed (1989 to 2019), a large part of the forest was converted to agriculture and grass/swamps (see Tables A4-A6 for the detailed LULC transition matrixes).
According to information obtained during field investigations (KII and field evidence), waterbodies increased after the construction of the Nashe Dam. During the construction of the reservoir, many farmers along and downstream of the Nashe stream were displaced to other agricultural places or towns. The abandonment of fields and the need for resettling in other areas caused a decrease in forests and an increase in built-up areas. In addition, poorly planned and long-term urban development and agricultural management strategies contribute to negatively affecting natural resources, causing a significant decline in the last decades ( Figure 3).

Model Validation
The LULC map of 2019, predicted from the 1989 and 2004 data, has been validated with the classified LULC map of the very same year ( Table 7), showing that the LCM model can effectively forecast LULC changes. The capability of LCM in predicting the 2019 LULC was assessed via K-indexes and other statistics in TerrSet (Table 8). All the values of the computed K-indexes (>80%) indicate good agreement between the projected and the actual LULC map [2]. The Disagree Quantity (0.0742) is greater than the Disagree Gridcell (0.0268), indicating that the model has a higher ability to predict the LULC in location (spatial) than in quantity for the Fincha watershed.

Future LULC
To forecast future LULC changes, it is necessary to account for the most important driving variables ( Table 9). As shown, all variables but the slope should be included in LCM, as the Cramer's V value of the slope is very low. The LULC maps for 2030, 2040, and 2050 were created via LCM, using the historical maps as a basis ( Figure 4). As observed in the past, an increase in areas covered by agriculture, built-up, grass/swamp, and waterbody is forecast, while a drastic decrease in forest and shrub should be expected, with a slower rate of deforestation in the decade 2040-2050 ( Figure 5 and Table 10). The decrease in natural forests, in combination with climate change, is likely to negatively affect the hydrological cycle of the whole Fincha sub-basin, as already shown by preliminary studies [7]. Moreover, the lower availability of wood (mostly timber) for construction will also impact the local economy.   In terms of transition probability (Tables A7-A9), areas covered by forest and shrubs are more prone to be converted into agricultural land, while built-up areas should be expected on the grass/swamp zones. This indicates that, in the future, agriculture and built-up zones will expand at a high rate since the other LULC classes will be converted to them. Conversely, forests and shrubs will decline at a significant rate.

Discussion and Policy Implications
As suggested in the review carried out by Regasa et al. [3], most of the studies of LULC in Ethiopia are on a local scale, analyzing how past variations affected water resources and socio-economic conditions in the region. This was mostly connected with the difficulty in obtaining proper field evidence (e.g., photos, description of LULC changes observed by local inhabitants). Therefore, very few works tried to forecast future LULC changes at the watershed scale, eventually providing new insights that can be useful for developing future basin-wide management strategies [74]. The present study was developed to fill this gap, estimating the Fincha sub-basin LULC for the next three decades (2030, 2040, and 2050) based on past information (LULC in 1989(LULC in , 2004, and 2019), to infer trends to be used in multiple ways. The LCM results point out that, in the coming decades, significant changes in LULC should be expected, mostly because of the ever-increasing pressure of humans in need of more land for settlements and cropland. Indeed, the local population is growing, and more natural resources are needed to satisfy their needs for food, energy, and construction materials [2,7,[75][76][77]. Apart from the direct consequences on the environment, the ongoing deforestation in the Fincha sub-basin is also causing an intensification of soil erosion, triggering sediment siltation in the various reservoirs located in the region, eventually leading to a decrease in reservoir efficiency in terms of water availability and hydropower production. Therefore, future studies should integrate the analysis of LULC changes with the simulation of soil erosion and sediment transport, to help local authorities better plan adequate management strategies to reduce siltation and guarantee the sustainability of local water resources, biodiversity, and socio-economy.
The reduction in forest cover pointed out in the present study compares well with similar research performed at the Ethiopian level [2,7,17,28,39,74,[78][79][80]. This trend was historically caused by the political situation, as, during the 1970s, the military regime proclaimed the nationalization of all rural land, abolishing private property [28,75,81,82]. However, the state was unable to adequately manage the land, and the majority of the forest was converted to settlements, agricultural land, and highly degraded areas because of the low level of land management practices. A similar approach has been taken since, but nowadays the land is property of the nations, nationalities, and people of Ethiopia according to Article 40(3) of the constitution endorsed in 1994 [83]. Because of the present policies, an increasing number of farmers are pushing for the expansion of their fields at the expense of the natural environment, in search of better socio-economic conditions driven by market-oriented choices [84][85][86].
As pointed out by Tolessa et al. [86], the reduction in land coverage and ecosystem services connected to natural resources such as forests are strictly connected with the present policy. As also shown in past investigation [87] and in the present work, the trend of decreasing forests to provide room for agriculture is very likely to continue in the future, also following the Ethiopian legislation. In fact, in Ethiopia, farmers have more legal rights over their land if they convert forest land into farmland, as the law stipulates natural forests as the property of the government. This forces farmers (legally or illegally) to convert forests into agricultural areas, as this guarantees them the use of the land for an indefinite period. As in other Ethiopian watersheds, the Fincha sub-basin was also affected by the 1975 land reform as, after it, grazing lands and forests were freely accessed for various uses. Tefera and Sterk [88] showed that, of the land potentially available for community use, cropland covered 77% of the whole region in 2001, indicating that further expansion cannot accommodate new farm families. This corroborates the hypotheses made in the present work, where future LULC was simulated assuming a trend similar to recent decades.
As with other Ethiopian basins [89], in the Fincha sub-basin, population pressure is a top driver of LULC change. To address this, in addition to policy changes, the local community should also start reforestation of the degraded forest area to assure the sustainability of the environment in the future, as also stressed by Senbeta [90].
To provide both decision-makers and local stakeholders such as farmers with more evidence on the importance of safeguarding water resources and ecosystem services at the basin scale, additional studies are needed, also taking advantage of information remotely acquired, such as satellite imagery, or considering different simulation algorithms and expanding the dataset of field evidence. Indeed, acquiring ground truth data and information from the local population, paramount for assuring a proper calibration/validation of the model, in locations that are hardly reachable or with unstable political situations such as the Fincha sub-basin, could be very difficult and expensive.

Conclusions
The present study investigated the historical LULC (years 1989, 2004, and 2019) in the Ethiopian Fincha watershed via a combination of satellite imagery and field support data. Based on such analysis, the Land Change Modeller was applied to forecast LULC over the next three decades (years 2030, 2040, and 2050). The 2019 LULC map was used for validating the LCM approach, comparing the forecast situation with the actual one derived from satellite images, indicating that the used Multi-Layer Perceptron (MLP) neural network of Markov chain (MC) has enough capability to predict future LULC.
The study results have shown that, over the last thirty years, the forest covering the Fincha watershed was mostly converted to agricultural and grass/swamp areas. An increase in zones covered by waterbody and built-up was also observable, mainly connected to the increasing human pressure and the construction of new hydropower reservoirs. This trend is recognizable not only in the study sub-basin but also in many Ethiopian basins [3], showing that LULC changes represent a major problem in the country.
As pointed out by the modeling results, in the future, a similar trend is more than probable. Indeed, if management strategies are not changed towards a more sustainable approach, also via proper reforms at the national level, an even more significant decrease in forest coverage should be expected in favor of new settlement areas and cropland. This change could help locals in sustaining their livelihood in the short term, but, in the medium/long term, the reduction in areas covered by forest will contribute to decreasing biodiversity and ecosystem services, as well as fostering soil erosion, with detrimental consequences such as reservoir siltation.
It is worth remembering that, to corroborate the results presented here and to reduce the uncertainties, additional data should be included in the study, mostly deriving from laborious and expensive field investigations. However, due to the current conditions of the study area, obtaining such information in the coming months could be very challenging. On the other hand, the increasing availability of commercial high-resolution satellite images could partially help in enlarging the dataset of field evidence, pointing out LULC changes happening at a more detailed scale. Therefore, in the future, a deeper analysis of satellite information is planned. Data Availability Statement: Additional data are available from the corresponding author, and partially available at https://dataportal.igf.edu.pl/dataset/land-use-land-cover-changes-percentage-inethiopia (accessed on 1 June 2022).
Acknowledgments: A great thank you goes to all the authorities and stakeholders involved in the field survey, which was paramount for validating our results. Moreover, we would like to thank the three anonymous reviewers for their useful comments, which helped us to better detail our research. We also thank the editors for their work in handling our manuscript and the revision process.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Accuracy Assessment of Historical LULC
The accuracy of the historical maps of LULC was assessed separately, showing the following results:

Appendix B. Transition Matrixes Historical LULC
Tables A4-A6 summarize the transition matrix between the three reference years. Detailed analysis and a discussion of these results are reported in Section 3.3.

Appendix C. Transition Matrices for Future LULC
Tables A7-A9 summarize the transition matrix of the forecasted LULC (see Section 3.5).