Quantitatively Assessing the Future Land-Use/Land-Cover Changes and Their Driving Factors in the Upper Stream of the Awash River Based on the CA–Markov Model and Their Implications for Water Resources Management

: Despite the rapid economic and population growth, the risks related to the current dynamics of land use and land cover (LULC) have attracted and the area of grassland decreased by 57.7%. In view of the greater increase in agricultural land and urbanization, as well as the decrease in shrubland, it means that the LULC of the region has changed. This research provides valuable information for water resources managers and land-use planners to make changes in the improvement of future LULC policies and development of sub-basin management strategies in the context of sustainable water resources and land-use planning and management.


Introduction
Land cover refers to natural structures such as vegetation and water surfaces, as well as manufactured structures that cover a certain region, whereas human activities that are related to the soil are called land use [1]. Change in land use and land cover (LULC) is an important factor in global environmental change, and it is crucial for regional development and land-use management to achieve sustainable development [2]. Due to rapid global economic and population growth and globalization, LULC is accelerating in many developing countries around the world [3]. As a result, land is being affected by competing demands and limited resources due to anthropogenic land degradation that severely disrupts various services and functions of the land [4].
Anthropogenic (human activity) affects the natural environment to a large extent through the dynamics of change (LULC). Consequently, the entire global ecosystem has been degraded due to the impact of human activities on natural resources [5]. Globally, in addition to population growth and global warming, the current dynamics of LULC are also of concern. For example, in East African countries and other regions strongly affected by population growth and climate change, human activities have greatly altered the natural landscape [6]. Future projections suggest that the world population will reach 9.1 billion by 2050, 34% higher than today, and agricultural production systems are expected to produce food for a growing global population [7,8]. Agriculture is a very important economic sector in many developing countries [9]. African forests play an important role in the development and transformation of a green economy, mitigating climate change and achieving the Sustainable Development Goals (SDGs) and the industrialization of the continent [10]. Ethiopia's economy is dominated by agriculture and is therefore highly sensitive to climate change [11]. Agriculture, forestry, and landscapes can significantly mitigate various aspects of sustainable development, including addressing the challenges of climate change, by increasing resilience and reducing vulnerability [12].
In Ethiopia, LULC has changed dramatically over the past few decades due to population growth, agricultural expansion, urbanization, and overgrazing. Similar to most other countries in sub-Saharan Africa, Ethiopia's population is rapidly expanding, with an estimated growth rate of around 2.8%, to a current population size of 108 million [11]. The service industry and agriculture are the main contributors to Ethiopia's GDP, both accounting for about 40%, and the industrial sector contributes the remaining 20% [11]. Smallholder farmers in Ethiopia face widespread problems of improper planting, overgrazing, and deforestation, leading to soil erosion and loss of soil fertility, water scarcity, lack of pasture and livestock feed, and a firewood crisis [13]. The main problem of deforestation is the overexploitation of natural resources; the clearing of agricultural land to meet the food needs of a growing population and the demand for fuel and building materials leads to a continuous reduction in forest area [14]. The rapid change in LULC is due to rapid population growth requiring more arable land, more trees for household firewood consumption, and more settlement area [15,16].
LULC changes have a significant impact on hydrology and water resources by affecting the hydrological cycle [17][18][19][20]. LULC and climate change are two major factors that are inherently anthropogenic and thus affect hydrological processes. Consequently, changes in hydrology and water components caused by climate change and anthropogenic activities, such as LULC, lead to many environmental problems [21,22]. LULC and climate change are the main drivers of land degradation by promoting soil erosion, especially heavy rainfall that accelerates soil erosion and bare vegetation cover [23]. Therefore, assessing LULC dynamics is important for effective environmental management, including effective water management practices.
In Ethiopia, due to changes in the environment and land use, soil erosion and consequent land degradation are considered major constraints to agricultural productivity and food security [24][25][26]. Likewise, water scarcity is increasingly common due to increased water demand, urbanization, economic development, land degradation, and climate change [27]. Climate and LULC changes are some of the main drivers of land degradation, which in turn will greatly affect water resources [19,23]. Therefore, LULC and climate-related natural disasters will affect people's lives, livestock, and economic development, which, in turn, will lead to water shortages in communities with different needs [28]. In Ethiopia, more than 350 million trees were planted in one day to minimize the above problems. This is part of the country's green legacy initiative to combat climate change and environmental degradation.
The main problem of the Awash River basin is land degradation, high population density, natural water degradation salinity, and wetland degradation [29]. This basin is the most important, intensively used, and environmentally vulnerable basin in Ethiopia [27,[30][31][32]. In addition to incorporating anthropogenic and natural changes in existing land degradation, new natural resource management strategies must be designed to address uncertainties associated with land use, land degradation, and climate change [33,34]. Any effort to address and reduce uncertainty in LULC and climate change assessments must support more effective adaptation and mitigation strategies for sustainable natural resource management [35]. A detailed and useful method for developing land-use classification maps is to use GIS and RS geospatial techniques [36]. GIS and RS are used to study LULC variation and mapping [37][38][39][40] and predict soil loss in agricultural watersheds [41,42], as well as in urban management [43], fire risk assessment [44], and hydrology and water resources management [45][46][47][48]. The cellular automata (CA) model can integrate remote-sensing and geospatial data, quantitative data (driving factors), and socioeconomic data [49]. Furthermore, a Land Change Modeler (LCM)-Markov chain model was applied to parameterize future land change predictions across multiple disciplines in different geographic regions. It has been evaluated and applied to simulate land-use change, and it has also been used in an Eastern Europe post-socialist land-use change model [20,21,[50][51][52][53], deforestation and forest monitoring [54,55], habitat and biodiversity [56,57], and extensively used in REDD+ project planning [58,59].
Ethiopia has fertile land and great agricultural potential; however, agriculture is still underdeveloped, and poverty persists, especially in rural areas [9]. Let us consider the current research portion of the USAB located within the Awash basin. Several studies have been conducted in the Awash River basin, but less attention has been paid to the assessment of changes in LULC in the study area. Thus, there is insufficient research in the USAB to identify and detect the role of specific types of LULC change and its drivers. In this context, there is an urgent need to estimate changes in LULC timing and predict future scenarios in the study area [60]. Currently, the sub-basin is particularly confronted with rapid population growth, land degradation, agricultural expansion, urbanization, and degradation of natural water bodies. Therefore, the main goal of this study was to analyze LULC dynamics and predict land in 2038 using spatial modeling (Markov chains and cellular automata) in the USAB's Land Change Modeler [61] by applying GIS and RS. The findings of this study provide useful information for water resource managers and land-use planners to develop sub-basin management strategies and improve future land-use policies within the framework of sustainable water resource management and land-use planning [62]. In addition, the results of this study can be used for future hydrological assessments of sub-catchments, with particular emphasis on surface runoff and river response.
In this study, it was hypothesized that the large changes in the LULC in the upper Awash River are associated with rapid population growth and human activities that meet the food security needs of the rapidly growing population through unplanned and inappropriate natural resource management practices.

Description of the Study Area
The Awash River is the second-largest river in Ethiopia and the main source of water for the rift valley in central Ethiopia. The study area is located in the upper part of Awash River basin, between 8 • 16 and 9 • 18 north latitude and between 37 • 57 and 39 • 17 east longitude with altitude varying between 1580 to 3396 m above mean sea level [63]. It lies upstream of the Koka dam and covers a total area of 9510 km 2 , which is considered as the major area where rapid population growth, urbanization, and industrialization expansion take place that force change in the LULC from time to time. The study area shares its boundaries with the rift-valley and upper-valley basins in the south, upper Blue Nile basin in the northwest, Omo Gibe in the west, and middle Awash River sub-basin in the east ( Figure 1). The basin received an annual mean rainfall of 1055 mm through the period of 1980-2017. Mean minimum, maximum, and average values of temperatures recorded at the basin during the season were 25.44 • C, 10.16 • C, and 17.80 • C, respectively [64]. Water in the sub-catchment area is mainly domestic water, mainly due to population growth and industrial expansion, and the demand for water has been increasing. Common types of land use in the upper reaches of the Awash River basin are cultivated land, shrubland, grassland, and woodland. Agriculture is the main land use in the study area. The main crops planted in this area are wheat, thrush, corn, sorghum, beans, and soybeans. Onions, tomatoes, and garlic are also some horticultural crops grown along river valleys and lake banks. In addition to growing crops, the livelihoods of communities in this area also depend on raising livestock, including cattle, sheep, goats, and poultry. The sub-basin also consists of different soil types. The predominant soil texture of the investigated area is sandy and silt clay loam. The most common soil types are Cambisols and Vertisols [65].

LULC Data Sources
Modeling of future LULC in the upper Awash River basin was based on baseline LULC data for the sub-basin (1984-2019). To model and evaluate the dynamics of LULC, multi-temporal and spatial data were collected. Digital elevation models (DEMs) and multi-temporal satellite imagery were used to develop future LULC thematic maps of the watershed based on historical LULC change patterns. Remote-sensing and geospatial data are reliable sources for identifying and understanding the dynamic drivers of LULC in any landscape [66]. Multiple cloud-free scenarios were obtained from the Global Land Cover Facility (GLCF) website (https://geog.umd.edu/ (accessed on 13 July 2020)) and the United States Geological Survey (USGS), website (https://earthexplorer.usgs.gov/ (accessed on 13 July 2020)). A basic summary of the multi-temporal satellite data used is shown in Table 1. The spatial resolution (30 × 30 m) of 1984, 2000, and 2019 Landsat imagery was used to create LULC maps and determine changes in the study area. In addition, various groundtruth data, document reviews, aerial photographs, and field observations were collected from the Ethiopian Survey and Mapping Agency (EMA), Ethiopian Ministry of Water and Energy (MoWE), and socioeconomic survey images from the validation satellite Landsat.

Landsat Data Acquisition
The Landsat TM, ETM+, and OLI images of 1984, 2000, and 2019 were selected, and GIS and IDRISI were used to analyze and process the changes in LULC, as shown in Figure 2. Remote-sensing image data of 1984, 2000, and 2019 were used, and atmospheric correction and radial calibration were performed. Relative geometric corrections were performed on the three images in order to remove geometric distortions caused by the Earth or sensor rotation [67]. Multi-temporal Landsat satellite imagery was used to analyze temporal and spatial changes to understand the dynamics of LULCs in the study area. Satellite imagery is geo-referenced using ground control points (GCPs) and projected using the Universal Transverse Mercator (UTM) system with WGS84 N 37 zones. In this study, false-color composite bands from Landsat images provided better visualization of surface features for further image processing ( Figure 2). In ArcGIS 10.4 software (Environmental Systems Research Institute (ESRI), Redlands, CA, the image subset setting process was performed using the image masking tool extraction based on the study area. The classification of LULC classes is based on using ground-truth validation data, Google Earth imagery, and visually interpreted detailed topographic maps to delineate and identify sample training locations. By defining a signature file and assigning the number of LULC classes, ArcGIS is used to apply the LULC classification through a supervised classification method (maximum likelihood algorithm), and then the geo-referenced terrain is used to verify that the LULC class map thus defined passes the ground-truth method. Supervised classification technique is more accurate than other existing image classification methods. Furthermore, it is the most commonly used pixel-wise method, which takes into account the spectral information of land-cover classes [68].

LULC Accuracy Assessments
In the process of image processing, accuracy evaluation is considered an important method for accurate image classification [69,70]. The verification/accuracy assessment of RS data is the final and one of the most essential steps to find the information value of the result data for the end user [40,71]. Using various statistical procedures to evaluate the accuracy of LULC classification can help understand the confidence of the results and determine whether the research goals have been achieved [40]. In this study, in order to verify the classification accuracy of LULC data in Google Earth, the most usually accepted stratified random point method was used. The accuracy of LULC is evaluated by comparing the classification map based on maximum likelihood classification with different ground-truth verification data. Map accuracy evaluation can be divided into user accuracy, producer accuracy, and overall accuracy and Kappa coefficient. Accuracy evaluation can be derived from various statistical processes of the error matrix, including user accuracy, producer accuracy, and the percentage (%) of the overall accuracy for resolving accidental errors [72]. Producer accuracy refers to the measurable map accuracy, which is used to indicate the probability that the reference pixel as the ground-truth feature is correctly displayed on the classified map or the probability that a certain land cover in a certain area of the ground is classified. The accuracy of the user is a measure of the commission error, which is calculated by dividing the number of correctly classified pixels in each category by the total number of pixels classified in that category. This represents the probability that the pixels classified on the map represent the ground category [73][74][75]. The accuracy of producers and users is based on the confusion matrix generated for accuracy evaluation using classified maps and ground-truth reference point data sources. The Kappa coefficient is a statistical measure of reliability or consistency between evaluators. It is a discrete multivariate technique used to present map accuracy evaluation research [74,76] and the Kappa statistics computed by the following Equation (1): where K is Kappa coefficient, N is the total number of observations, r is the number of rows and columns in the error matrix, x ii is the number of observations in row i and column i, and x i+ and x +i are the marginal totals of row i and column i, respectively [76][77][78].
The LULC transition matrix: The land-use transition matrix is widely used to evaluate the quantitative altering pattern of the transition of LULC changes particularly for RS and GIS studies [79]. The land-use transition error matrix refers to the LULC change into three stages, for instance, from 1984, 2000, and 2019 in this research. The transition matrix is computed using the following Equation (2) [80].
where p ij specifies the area in transition from landscape i to j, and each element in the transition matrix is categorized supposing (1) p ij is non-negative and (2) ∑ n j=1 p ij = 1.

LULC Change Detection
The area of the classified LULC map is calculated in square kilometers (km 2 ) and quantifies the change between each prescribed period. In addition, the absolute and relative changes of the classified LULC class are calculated by considering the difference between two consecutive LULC maps. The LULC change matrices between each prescribed period were developed to understand the temporal and spatial changes of LULC classes. In addition, the rate of change of LULC over time (km 2 /year) is calculated by dividing the area difference between the two time periods by the length of the year. Furthermore, the following Equation (3) is used to estimate the relative change in the LULC categories classified between each successive period: where X (%) is a relative change of LULC classes between earlier period X 1 and later period X 2 in km 2 , and X 3 is a time interval between X 2 and X 1 . The positive values imply an increase in area coverage of LULC class whereas the negative values indicate a decrease in area extent. This research adopted the post-classification change detection technology implemented in ArcGIS 10.4. The method framework for LULC classification and change detection is shown in Figure 3.

Prediction of Future Land Use/Land Cover
Long-term (2038)-scale changes in the LULC categories of the watershed were predicted based on the historical changes of the watershed from 1984 to 2019. To simulate future watershed LULC maps, the Land Change Modeler (LCM) module in IDRISI was used. TerrSet (formerly IDRISI) is an integrated remote-sensing software and geographic information system (GIS) developed by the Clark Laboratory at Clark University for the display and analysis of digital geospatial information [81]. We used LCM networks for predicting future LULC maps using Markov chains (MCs) and multilayer perceptron neural networks (MLPNNs) to determine historical transitions of LULCs developing in the basin from multi-temporal Landsat imagery. Based on past historical conversion potential, MC quantifies the LULC categories that are likely to transition from t 2 to the expected forecast period t 3 . The simulation process for future LULC classes is based on a transition probability matrix that shows the probability of each LULC class transitioning to other LULC classes. Let Y 1 and Y 2 be the states of the LULC classes at period t 1 and t 2 be the change matrix with elements T ij (ith row and jth column) which shows the transition Y 1 to Y 2 during t 1 and t 2 [82]. Then, the transition probability matrix can be obtained as: If t 3 is an even multiple of the time duration between t 1 and t 2 , then the transition probability matrix for t 3 can be found by powering the past transition probability matrix. If t 3 is not an even multiple of the periods between t 1 and t 2 , the transition probability matrix can be obtained by the interpolation technique [82].
The simulated LULC map can be used to quantify and visualize future LULC classes in the basin. The ability of LCM to predict future LULC maps in the study area was validated by both 2019 forecasts and historical LULC maps. The disagreement between the predicted and historical LULC maps was measured in terms of two parameters, viz., allocation disagreement and quantity disagreement [83]. Quantity difference calculates the area difference of LULC categories in historical and forecast LULC maps. The spatiotemporal allocation differences of the LULC classes were measured by using the parameter allocation divergence. The quantity and spatial disagreement space vary within [0, 1]. The sum of allocation disagreement and the spatial disagreement signifies the overall disagreement of the predicted LULC maps with the historical maps [82].

CA-Markov Model Approach
The Land Change Modeler (LCM) was built as an empirical parametric land-change prediction tool to support a wide range of planning activities [84]. In this study, the CA-Markov model in IDRISI was used to simulate future LULC in the upper Awash River basin. CA-Markov is one of the best tools for predicting future changes in land-use parameters [85]. The CA-Markov model combines Markov chains and cellular automata to predict trends and characteristics of LULC over time [86]. Furthermore, CA-Markov is one of the design support tools to analyze the spatial distribution and temporal variation of LULC [87]. In addition, the model is widely used to describe LULC dynamics, urban sprawl, plant growth, watershed management, and forest-cover modeling. Based on the analysis of historical LULC changes, this model develops an empirical model of the association between LULC transitions and a set of explanatory variables [62,86].
Markov chains have good statistical power to predict change probabilities, and cellular automata are considered to be a powerful method for reading spatial patterns of change [88][89][90][91]. The Markov chain model refers to the change of LULC from one time to another to predict future changes [91,92]. The prediction of land-use changes can be obtained as: where S(t) is the system status at the time of t, S(t + 1) is the system status at the time of t + 1, and P ij is the transition probability matrix [80,92]. The methodological framework applied for the CA-Markov model is shown in Figure 4. Transition suitability maps for each LULC category are derived using a series of factors and constraints. This is done through an expert's understanding of the interactions and effects of factors and a LULC-based multi-criteria analysis (MCE) [91]. Constraints, represented in the form of Boolean images, are criteria (0 to 1) that limit the variation of LULC, with a value of "1" representing the area suitable for suitability analysis, and "0" representing the area that controls the suitability analysis (e.g., water bodies) [62]. In contrast, the factor is usually a distance criterion that provides a degree of applicability to regional variation [93], for instance, distance to rivers, urban areas, or roads [62,[94][95][96]. Other factors include elevation (which is a good predictor of agricultural areas) and slope, which determine the usefulness of land to humans [49]. The constraints and factors used in this study are shown in Figure 5.  Factor and constraint maps-elevation (a); slope (b); roads (c); distance to rivers (d); distance to urban areas (e); distance to roads (f).

Model Calibration and Simulation
The scenario-driven CA-Markov model method was adopted and then calibrated and verified to simulate future LULC changes. Based on the Markov cellular automata model, the LULC changes were simulated under two different historical business-as-usual (BAU) scenarios. The images of 1984 and 2000 were used for the calibration and optimization of the Markov chain algorithm, and the images of 2019 were used to verify CA-Markov's predictions [97]. The LULC change between two time periods (t 1 and t 2 ) was modeled using real-land-cover maps to predict the land-cover map at t 3 and verify the model. We checked the simulated land-cover map (t 3 ) against the real map [98]. In addition, the land-cover map in 1984 is the earliest image (t 1 ), and the latest land-cover map in 2000 (t 2 ) was used to simulate the projection map of 2019 (t 3 ) and was checked against the real map in 2019. Compare the simulated LULC map with the actual map to verify the accuracy of the model. In order to predict the LULC maps in 2019 and 2038, the transition probability matrices for 1984-2000, 2000-2019, and 1984-2019 were calculated. The Markov and CA-Markov models were used to predict the LULC map in 2038.

CA-Markov Model Validation
Model validation is the main part of the CA-Markov model. It consists of test models on data not used to build the model [2]. The predicted LULC is considered consistent only when it is validated using existing ground data sets. Using the CA-Markov model and the verification model in the cellular automata, two LULC raster's were verified; the first was the actual LULC in 2019, and the second raster was predicted. The same process was applied when predicting the LULC map in 2038, but in this case, the iterations were given in order for the time difference between 1984 and 2000. In this study, the accuracy of the CA-Markov model was evaluated using Kappa statistics, including Kappa location (K location ), Kappa standard (K standard ), and Kappa for no information (K no ) [95,99], and Equations (7)-(9) [100] were used to calculate them. The level of the Kappa index is usually between 0 and 1. Therefore, the consistency level r of the Kappa index is as follows: Kappa < 0.5 indicates a rare agreement; 0.5 < Kappa < 0.75 indicates a moderate level of agreement; 0.75 < Kappa < 1 indicates a high degree of agreement; and Kappa = 1 indicates complete agreement [99,101]. We used the "validate" module to compare the 2019 classified LULC with the 2019 simulated LULC. Under acceptable Kappa statistics, the CA-Markov model is suitable for simulating the future LULC map in 2038, using the transition probability from 2000 to 2019 and the classification map in 2019 as the base map.
where N n is no information; M m , N m , and P m are medium grid-cell level information; and P p is perfect grid-cell level information across the landscape. However, when compared with the actual classification raster and the LULC prediction raster, the Kappa statistics give the degree of agreement between the raster and its probability.

Investigating the Drivers and Consequences of LULC Changes
In addition to remote-sensing data, sub-basin socio-economic data were collected and analyzed through key informant interviews, household surveys, and focus group discussions (FGDs) to verify the accuracy of the classification of LULC and further understand the possible main driving forces and consequences of watershed land use change. In this study, the multi-stage sampling method was used to select kebeles purposefully first because they represent three agro-ecological regions (Dega, Weyina Dega, and Kolla). FGD was made up of groups of elders, community leaders, women, and youth, who deliberately chose to participate in discussions to provide information related to current and past changes in LULC and to support the reliability of the survey. The socio-economic characteristics of the interviewees are shown in Table 2. We used a questionnaire with semi-structured questions to assess the local community's perception of LULC changes, their driving forces, and consequences. This study identified a total of 1800 households and determined a sample size of 10%, with 95% confidence and 10% accuracy. Therefore, a total of 180 sample households were selected from the list using the systematic random sampling technique [102]. Of the 180 households included in the sample, 80% were men, and the remaining 20% were women. The age of the head of the household interviewed was between 25 and 75 years old. In this study, the community's perception of the driving force of LULC changes focused on socio-economic characteristics, the reasons for LULC changes, and the impact of LULC changes requiring households to use "agree", "disagree", and "uncertainty", as well as "increase", "decrease", and "change", to assess the level of family perception. The interviewees were asked to explain how they viewed the LULC changes in the sub-catchment during the different time periods evaluated. In addition, qualitative data collection methods (personal observations, focus group discussions, and key informant interviews) were used to collect data on the local residents' views on the proximate and root causes of changes in LULC in their area. Finally, the main driving force, consequences, and direction of the LULC change were determined.

Land Use/Land Cover (LULC) Classification
The maximum likelihood classification algorithm in IDRISI17.02 was used for supervised classification, and the spatiotemporal pattern of LULCs in the sub-watershed was assessed. According to the LULC classification system of the Food and Agriculture Organization [103], six main types of LULC (agricultural land, forest, shrubland, grassland, water body, and urban) were identified. From 1984 to 2019, the supervised classification method (maximum likelihood algorithm) LULC classification analysis showed that the study area covered various land features (agricultural land, shrubs, forests, grassland, cities, and water bodies). Agricultural land refers to the area of cultivated land used to grow various crops such as corn, sorghum, teff grass, beans, barley, millet, and wheat. Shrubland types include dwarf tree species; sparsely located forest and multi-stem woody vegetation are considered as shrubland. Forest is a large area of land, covered with trees and plants, which have a closed canopy cover and appear greener throughout the year. Grassland is mainly composed of grassland, wetland, and non-herbous plants. Urban refers to built-up areas that include residential cities, towns, suburbs, commercial, industrial, transportation, and roads. Water bodies are areas covered by rivers, reservoirs, lakes, dams/ponds, streams, and seasonal waterlog areas. The obtained LULC maps for the three time periods of 1984, 2000, and 2019 are shown in Figure 6. The LULC classification map shows that the total land area of the upper reaches of the Awash River is 9510 km 2

Accuracy Assessment
The accuracy evaluation of the classified LULC map, especially the producer accuracy, user accuracy, overall accuracy, and Kappa coefficient, were derived from the statistical process of the error matrix. Table 4 shows the LULC map accuracy evaluation results in the three periods of 1984, 2000, and 2019. The overall map accuracy of using the error (confusion matrix) is 85%, 88%, and 90%, respectively. We calculated the Kappa coefficient of the LULC map for each category to measure the accuracy and confidence of the results. The Kappa index of the LULC classification map results in the three periods was 0.77, 0.79, and 0.89, respectively. This showed quite good overall accuracy and was accepted for subsequent change and detection analysis [104].

CA-Markov Model Validation
The LULC changes along the upper reaches of the Awash River are not uniform. In order to be able to use the CA-Markov model to simulate the future LULC in the upper reaches of the Awash River and validate the ability of the model, we first used the transition probabilities and the transition area between the LULCs of 1984-2019 to simulate the LULC in 2019 (Figure 7). Kappa variables were used to compare simulated and classified LULC in 2019. The simulation shows that the Kappa variation evaluation is as follows: 0.76 for K no , 0.88 for K location , and 0.78 for K standard , showing a high level of agreement. This shows the strength and reliability of the model in simulating future LULC changes in the upper reaches of the Awash River. In addition, visual analysis shows that the LULC category in the 2019 simulated LULC is more consistent with the 2019 classified LULC category. The calculated LULC change differences between the classified LULC and the simulated LULC in 2019 are shown in Table 5. The results show that the shrubland and urban areas were seriously overestimated, while the agricultural and grassland classification showed slightly lower results than the simulation. The calculated difference between the classified and simulated LULC shows that the agricultural, grassland, and water body categories are underestimated by approximately 28.06%, 26.65%, and 7.06%, respectively. On the other hand, the CA-Markov model proved to overestimate LULC categories, such as shrubland (222.23%), forest (23.50%), and urban (4.89%).

LULC Changes and Markov Probability Transition Matrices
This study identified six LULC categories in 1984, 2000, and 2019, namely, agricultural land, forest, grassland, shrubland, water body, and urban. The area and relative change of the LULC category on the temporal scale are shown in Figure 8a,b. Temporal analysis of the LULC map of the basin shows that the accelerated development activities implemented have led to a complex and violent LULC transition. The post-classification change matrix comparison technique clearly shows the area that has undergone a transition from one LULC class to another LULC class between t 1 and t 2 . Table 6 shows the Markov transition probabilities for the period 1984-2038, showing each category in the LULC map within that period (for example, 1984-2000, 2019-2019, and 2019-2038). That is, a given LULC category may change from one LULC category to any other category at any given time [62]. The LULC change matrix represents the transition between each LULC category during the time intervals t 1 and t 2 [49]. The change matrix in the vertical column is the area lost from one particular LULC class to another during t 1 and t 2 [105]. The area under each LULC class is represented by a diagonal line, and the area under each LULC class is the sum of the elements in the row at t 1 . The values on the diagonal represent the transition probability matrix of each land-use category that remains unchanged from time 0 to time 1 [106]. For instance, urban and agricultural land categories are highly resistant to changes in other categories, and the probability of remaining unchanged during 1984-2000, 2000-2019, and 2019-2038 is more than 70%. Between 2019 and 2038, the highest percentage (>85%) of the urban category that remains unchanged is observed. Classifications such as shrubland and grassland also show resistance, and more than 51% of the classifications may remain unchanged in the next period, while other classifications show less resistance, and thus they are more likely to change to other LULC classes in the following period. These include forests, grasslands, and water bodies, of which less than 38% and less than 60% remain unchanged for a period of time. On the contrary, the water body category showed a low probability and a high probability of remaining unchanged during the 36-year period, and more than 60% of the category remained unchanged during the period 1984-2000. Compared to showing persistence, the Markov transition probability matrix shows how each LULC class becomes another class [62]. For example, during the period 2019-2038, 45% of the water body category may be changed to the agricultural category. In addition, between 2000-2019 and 1984-2000, approximately 39% of grassland and 40% of shrubland are likely to be converted to agricultural land, which indicates that they are less likely to continue to exist during these periods.

Simulation of Future LULC Dynamics
Using the 2019 LULC as the base map, the transition suitability images and Markov transition probabilities generated from 2000 to 2019 simulated the LULC changes in the upper reaches of the Awash River basin in 2038, and the results are shown in Figure 9. Table 7 shows the dynamics of LULC changes in the spatiotemporal and area from 2019 to 2038. The results of the LULC change show that in 2019, 84% of the total land area of the sub-watershed was agricultural land, and it is expected to reach 78% by 2038. It is estimated that the level of shrubland will increase from 1025.64 km 2 (10.78%) in 2019 to 1524.54 km 2 (16.03%) in 2038. Similarly, the percentage of urban land is expected to double, intensifying from 163.08 km 2 (1.71%) in 2019 to 341.11 km 2 (3.59%) in 2038. The analysis results further indicate that the classification of agricultural land, forests, grasslands, and water bodies is expected to decrease in the next three decades. In particular, forest areal coverage will decrease by 59.79 km 2 (76.52%), and water bodies will decrease by 10.97 km 2 (51%). Similarly, grassland and agriculture have decreased by approximately 571.33 km 2 (7.14%) and 34.76 km 2 (15.42%), respectively.

Driving Forces of LULC Changes in the Upper Stream of the Awash River Basin
LULC dynamic analysis uses the views of local communities to accompany the interpretation of remote-sensing data. The results of key informant interviews, household surveys, focus group discussions, and field observations show that both natural processes and man-made LULC have changed. From a series of different driving factors, interviewees believe that nine activities related to natural and man-made processes are the main driving factors of LULC changes in the study area (Table 8). Therefore, more than 87.8% of the interviewees stated that the rapid growth of population pressure is one of the main driving forces for changes in the LULC in the study area. Similarly, approximately 92.8%, 67.8%, 70.0%, 80.6%, and 85.6% of the respondents argued that increasing agricultural activities, reduction of soil fertility, use of trees for firewood, charcoal, and construction, wood extraction, and settlement and urban expansion, respectively, were some of the perceived environmental consequences that caused LULC changes. Therefore, the driving force of LULC change in the sub-watershed is essentially related to human influence. Thus, the driving forces of LULC changes in the sub-basin are essentially linked to human-induced influences. The results show that between 1984 and 2019, the area of agricultural land and urban/built-up areas gradually increased, while forests, shrubs, and grasslands decreased. It is supposed that the impact of human activities such as forest clearing for agriculture and demand for house construction and firewood had a greater contribution to LULC changes in the socio-economic context of Ethiopia [102]. As shown in Table 8, 73.9% of the respondents indicated that overgrazing (livestock pressure) and the advantages of free-grazing systems are the main driving forces of LULC change and land degradation. In addition, 83.3% of the interviewees believe that the increases in drought and climate change are other factors leading to changes in LULC.
Changes in LULC were determined by local communities, and their responses are shown in Table 9. The perception of changes in LULC highlights that a large number of interviewees are aware of the long-term dynamic process and driving forces of LULC and its influence over the past three decades. Perception results confirm that 86.7% and 90.6% of the respondents believe that the area of farmland and urban land has increased since the 1980s. On the other hand, the respondents believe that the areas of forest (84.4%), shrubs (81.1%), grassland (69.4%), water bodies (78.9%), swamps (88.3%), and grassland (69.4%) have declined in the upper Awash River sub-basin. Thus, the indigenous knowledge perceived from the local interviewees confirmed the LULC change information extracted from the interpretation of the remote-sensing data set.
Therefore, understanding the causes and consequences of land-use changes is crucial for scholars, decision makers, and water and land managers, as it helps to take appropriate actions for future water management.

Validation of CA-Markov Model
The CA-Markov model was used to simulate the future LULC in the upper reaches of the Awash River basin. Model verification was conducted by comparing the classified and simulated LULC in 2019 to test the performance of the model in simulating future LULC changes. It is observed that the Kappa index estimate is between 0.76 and 0.88, which represents a high level of consistency [83,99] and the applicability of CA-Markov, which can accurately simulate the future LULC of the upstream Awash River watershed. However, the CA-Markov model proved to overestimate shrubland and urban areas while underestimating the LULC categories of agricultural land, forests, grasslands, and water bodies. Several researchers reported an acceptable CA-Markov model Kappa index [65,. Therefore, the CA-Markov model simulation is suitable for accurately predicting future changes in LULC in the upper reaches of the Awash River basin.

LULC Changes Observed from 1984 to 2019
The temporal LULC maps of the USAB were obtained by supervising maximum likelihood classification. The generated map helps to identify the dynamic pattern of visual LULC changes and quantify current and future changes. The study area sub-basin was categorized into six LULC classes, namely, agriculture land, forest, grassland, shrubland, water bodies, and urban. From the results of the above-mentioned LULC changes, it can be found that agriculture and urban areas have undergone great changes from time to time. The total area of urban and agricultural land from 1984 to 2019 was 3065.72 and 138.76 km 2 , respectively. Therefore, about 80-85% of Ethiopians are engaged in agriculture, mainly in self-sufficient and rain-fed agriculture and livestock production, and the country's agricultural land has increased [11,107]. The results show that during 1984-2019, agricultural area and urbanization increased significantly at the expense of other LULC types. From the results, it was found that agricultural land is the predominant land type of the sub-basin. In the past 36 years, the total area of forest cleared between 1984 and 2019 reached 616.85 km 2 . This means that about 88.7% of the existing forest coverage in 1984 was cleared. The largest expansion of agriculture and urbanization occurred between 1984 and 2019 and is estimated to be 62% and 570%, respectively. The results of the study show that the area of shrubland has dropped from −6% in 1984 to −68.6% in 2019. Similarly, grassland dropped from −47.5% in 1984 to −57.7% in 2019, and water bodies dropped from 41.5% in 1984 to 65.5% in 2019. It was found that the main driving forces of this transition were rapid population growth and urbanization, which led to deforestation, rainfall shortages, and farmland, which in turn exacerbated food security issues. The ever-increasing population has been driving the increase of agricultural land nationwide [5].
From 1984 to 2019, the drastic changes in the watershed were the expansion of agriculture and urbanization and the shrinking of forests, shrubs, grasslands, and water bodies. This shows that agriculture is expanding on a large scale. In Ethiopia, population growth has led to agricultural expansion, urbanization, and overgrazing; LULC has undergone tremendous changes in the past few decades [50,51,[108][109][110]. Similarly, in the past 36 years, the LULC of this sub-basin has undergone major changes. The results show that between 1984 and 2019, the urban area increased significantly. This shows that the LULC of the urban area category has undergone tremendous changes, exerting incredible pressure on non-urban areas, especially agricultural land. Through the construction of residential units, industrial and commercial units, road networks, sidewalks, ports and leisure facilities, and other impervious surfaces, the rapid expansion of urban areas has led to the continuous expansion of building surfaces in different corners of the city [36].
The change detection results show that the extent of LULC changes that occurred in shrubs, forests, and grasslands was concentrated in agricultural land. These changes indicate that the intensification of agricultural land is the result of population growth. Rapid population growth, human migration to urban areas, and detailed agricultural farming systems are the main drivers of changes in agricultural land use in the entire sub-catchment [16,38,[115][116][117][118][119].
Changes in land use from non-urban to urban or from forest to agricultural land have resulted in the loss of land cover for many different types of vegetation. The most serious consequence of changing land use through urbanization and agricultural expansion is the reduction of natural and expansion of agricultural land and the increase of hard surfaces in built-up areas [111]. This leads to soil erosion and loss of soil biodiversity, which leads to a decline in soil fertility, which in turn reduces agricultural productivity [112]. Changes in LULC are driven by human actions, which in turn drive changes that change the availability of human and livestock products and services [113].

Future LULC Changes for 2038
It is expected that the LULC in the upper reaches of the Awash River basin will continue to change, and the reduction and expansion of different LULC classifications may have an impact on the environment, mainly water resources. It is expected that by 2038, the types of shrubland and urban land will expand, while the types of agricultural land, forest land, grassland, and water bodies are expected to decrease. Agricultural land will be reduced and transformed into industries and urban areas if the current rate of industrialization and urbanization continues, which will have an impact on the water resources and environmental integrity of the sub-basin [67,[116][117][118][119][120][121][122][123][124]. In addition, from 2019 to 2038, urban and shrubland are expected to increase significantly to 109% and 49%, respectively. The growth of urban areas in the sub-watershed is a sign of population growth, mainly from rural to urban migration to Addis Ababa, because settlements in nearby towns and cities tend to promote urbanization, which is also observed in some countries [16,[114][115][116][117].
LULC losses are mainly due to the expansion of agriculture and urban areas. It is observed that the types of water bodies in the upper reaches of the Awash River basin are decreasing [60]. Changes in the water level of the sub-basin have been observed, and it is reported that the Koka dam has decreased water levels, which has exacerbated the shortage of hydropower and irrigation water supply in the surrounding communities [118,119]. Thus, LULC changes in the upper stream of the Awash River basin may disturb water resources such as runoff and streamflow [120] and could reduce the dam water levels [32]. As the population increases and LULC is expected to change in the next few decades, this sub-basin will face a high demand for water resources. Thus, there is a need for wellorganized and environmentally friendly industrialization and urbanization planning in both rural and urban areas.

Driving Forces of LULC Changes
As shown by the remote-sensing LULC change analysis, in the past 36 years (1984-2019), agricultural land and urban land in the upper reaches of the Awash River basin have increased significantly. Similarly, a large number of respondents believe that the main reasons for the changes in LULC are human disturbances mainly due to agricultural expansion (92.8%) and population growth (87.8%). Thus, the direction of LULC changes perceived by the respondents was consistent with the result obtained from remote-sensing image interpretation [121]. The driving force of LULC changes in the basin is essentially related to man-made influences. The results show that between 1984 and 2019, the area of agricultural land and urban/built-up areas gradually increased, while forests, shrubs, and grasslands decreased. Population growth has led to a substantial increase in the demand for food and the expansion of agricultural land by affecting the area of uncultivated forest land. It is speculated that in the socio-economic context of Ethiopia, the impact of human activities, such as agricultural deforestation and the demand for firewood and housing construction, will make a greater contribution to the changes in LULC.
Understanding the patterns and driving factors of LULC change is essential for the rational and specific planning of sustainable land management [122]. LULC is a complex and dynamic process that can be caused by many cooperative processes ranging from various natural factors to socio-economic dynamics [123]. Due to different driving factors, major LULC changes occurred in the study area (reduction of grassland, expansion of farmland and settlements, expansion of degraded land, and reduction of shrubland) [33]. Demographic drivers are part of the reason for the expansion of farmland and settlements, reduction of grassland, reduction of shrubs, and expansion of degraded land [33].
Population increase seems to be the main driving force of LULC changes, which is mainly manifested by expanding cultivated land at the expense of vegetation cover [124]. This has led to the expansion of cultivated land, charcoal, fuelwood extraction, and other wood products (such as housing construction products related to population growth) as the main direct drivers of the loss of natural vegetation in the study area [124].

LULC Change Implications for Water Resources Management
The majority of sub-Saharan Africa (SSA) countries are susceptible to LULC changes. LULC change is increasingly burdening the future water bodies. The gradual deterioration of water bodies can lead to landscape component degradation (flora, soil, and fauna), leading to desertification, which in turn leads to health problems, poverty, and loss of biodiversity [125]. Deforestation is one of the changes in LULC, which is considered to be the main cause of changes in hydrological processes such as surface runoff, sediment production, evapotranspiration, groundwater, infiltration, lateral flow, and rainfall interception [126][127][128][129][130]. Changes in LULC are driven by human interaction, which in turn will drive changes that change the availability of water resources. Under various circumstances, natural systems have become agricultural land to feed a growing population [131]. It has been noted that changes in LULC will have an impact on natural resources, especially water resources in semi-arid areas [62]. Land-use and climate change are the two main driving forces that affect the hydrological processes of the watershed [125]. However, the impact of LULC changes on the hydrological system is a complex process, depending on the size of the affected area, the type of land use, and its landscape. The availability of future water resources depends, to a large extent, on land-use planning and management in a constantly changing environment. On the other hand, continued human behavior continues to modify LULC to meet increased demand, especially due to the significant increase in population and the development of better facilities [21].
The intensification of agricultural land and urban areas is related to the reduction of vegetation cover in most areas [124,132], which leads to soil compaction in farmland and an increased resistant surface in urban areas [133,134]. Since the time interval of crop evapotranspiration is relatively short compared with forest or shrubland, it is well known that agricultural land will reduce evapotranspiration However; these LULC classes exhibit low water absorption and release large amounts of water to the outlet of the sub-basin, resulting in high water production. LULC loss also indicates reduced infiltration, which means low base flow and groundwater recharge [135], high surface runoff and peak flow, and low evapotranspiration [136,137]. However, semi-arid areas such as the upper reaches of the Awash River basin have a lot of rainfall and high temperatures, and surface water evaporation is still relatively high. For example, during the period 1970-2008, the surface water evaporation in the Koka basin was estimated to be 1.37 × 109 m 3 [27], which is almost the same as the annual rainfall of 1.97 × 109 m 3 , which leads to a low level of the water available in the basin area. Therefore, high evapotranspiration leads to a decrease in base flow and groundwater, which leads to a decrease in water production and flow; reduced water output may lead to water shortages and extreme water supply anxiety [62].
Assessment of river basins requires LULC change assessment to explore ecological and hydrological conditions for future sustainable water resources management. Several studies have studied the impact of LULC on water resources in many countries around the world. For instance, in eastern and southern Africa, changes in LULC are the main driver of reducing river flow and surface runoff [138], leading to water shortages. Similarly, the study reported by [139] showed that LULC changes are one of the key drivers of hydrological changes in the watershed. LULC change is a complex source of pressure that threatens the sustainability and management of water resources [136,140]. The impact of the LULC change pattern on the sub-basins connecting coastal areas and highland areas has led to economic, social, political, and environmental problems at the national, regional, and local levels in many countries. LULC changes have serious consequences for the natural environment because LULC is directly related to land degradation and leads to environmental changes. Therefore, understanding and paying attention to the existing LULC changes is very important for the policy decision making and regulatory action formulation of future LULC activities.
Changes in LULC not only affect water quality but also affect the quality of water resources through different mechanisms. Changes in LULC, such as the fragmentation and degradation of terrestrial landscape, adversely affect the quality of aquatic ecosystems. Freshwater ecosystems are some of the most threatened ecosystems on the planet, facing human and environmental pressures [141]. For example, due to the increasing pressure and scarcity of freshwater resources, the pollution of freshwater resources is attracting global attention [142]. River water quality is affected by rapid urbanization [143], and agricultural chemicals may leak or seep into the soil due to the expansion of farmland, causing pollution of nearby water sources [62]. In addition, water pollution is common in rapid urbanization, especially in developing countries, because of the lack of proper sanitation and waste treatment facilities (industrial and household). Studies have shown that there is a positive correlation between LULC changes and water quality parameters, indicating that LULC changes contribute to water quality changes, and thus, frequent water quality monitoring and LULC planning and management are recommended to curb watershed pollution [142]. With the predicted changes in LULC in the upper Awash River basin, water resources are expected to continue to be affected; therefore, appropriate actions must be taken to reduce the impact of LULC changes on water resources, especially in these water-scarce sub-basins.
The CA-Markov model showed great efficiency in simulating LULC changes in the upper reaches of the Awash River basin, indicating that the region is vulnerable to changes and requires extensive water resources management and land-use planning in the future. To confirm the capabilities of the model, key informant interviews, household surveys, field observation, and focus group discussions were used to determine the consequences and driving factors of LULC changes in the study area. Thus, the LULC driving factors evaluated through local community surveys showed the applicability of the CA-Markov model in simulating future changes in LULC.
Therefore, the LULC changes caused by urbanization and agricultural intensification have a huge impact on the water resources of the river basin and should be considered in future water resources management.

Conclusions
This study utilized reliable historical data to analyze the dynamics and magnitude of LULC changes in the upstream Awash basin (USAB) and used the CA-Markov model to predict future LULC changes, contributing to a better understanding of possible LULC changes and their driving factors. To confirm the results of remote-sensing image interpretation, the consequences and drivers of change in USAB LULC were assessed through key informant interviews, field observations, household surveys, and focus group discussions. Assessment of historical LULC dynamics is critical for future agricultural and water resource management under agricultural intensification and urbanization. This research plan examined the application of RS and GIS to the dynamic characteristics of LULC in a subwatershed and showed that LULC in the region changed significantly during the reference years of 1984, 2000, and 2019. This study showed that the CA-Markov model can simulate future changes in USAB LULC by providing a reliable 2038 historical LULC model.
The results show that the most significant changes occurred in the spatiotemporal scale of LULC dynamics. From 1984 to 2019, basin-wide agricultural and urban cover increased significantly, while shrub and forest cover decreased slightly. The USAB is expected to continue to experience LULC changes in the future, primarily the expansion of urban land and agricultural land categories. It is estimated that by 2038, the area of shrub and urban sub-watersheds will be about 1524.54 km 2 (16.03%) and 341.11 km 2 (3.59%), respectively, while the coverage area in 2019 will be about 1025.64 km 2 (10.78%) and 163.08 km 2 (1.71%), respectively. However, LULC categories such as agricultural land, forest, grassland, and water bodies are expected to decrease. Population growth and urbanization (mainly from rural-to-urban migration to Addis Ababa) are expected to lead to urban expansion, and the reclamation of agricultural land to meet the food needs of a growing population will lead to expansion of agricultural land; as a result, forest areas continue to reduce. The rapid expansion of agricultural land and urban areas, as well as the reduction of shrubland, forests, and grasslands, is evidence of the human impact of LULC dynamics. Extensive LULC dynamics are human-driven, which in turn affects humans and alters the availability of natural resources, including water, vegetation, soil, and livestock. Therefore, the extent of LULC variation perceived by local community respondents is consistent with the interpretation of remotely sensed LULC images.
Changes in LULC projected in the USAB will have an implication on water resources, such as reductions in water production (surface and groundwater), river flow and water quality, and evaporation in the watershed. The future availability of water resources depends, to a large extent, on the planning and management of land use in this changing environment. Therefore, it is necessary for water resource managers and land-use planners to understand the impact of changes in land use in sub-basins on water resources in order to improve future water resource management. LULC dynamics are caused by population pressures that lead to urbanization and agricultural expansion through unplanned and inappropriate resource management practices to meet the food security needs of a rapidly growing population. Therefore, proper planning needs to be implemented to minimize landuse dynamics; especially in developing countries, whose resource bases are deteriorating and should be strengthened to feed their growing populations.
Agriculture-related LULCs are unavoidable for developing countries such as Ethiopia whose economies depend on agriculture, as most of their arable land is located in large river and lake basin areas. Thus, in the wake of the 2030 Sustainable Development Goals (SDGs), this study draws attention to the need for immediate action to measure sustainability indicators of land and water resources across the upper Awash sub-basin. Evidence of LULC changes and drivers of change are essential for future planning projects, as they provide more information about LULC. More encouraging results on LULC changes and their drivers were obtained using long-term coverage data. The results of this study can be used for future hydrological impact assessments in the sub-basin to contribute to the sustainability of the Koka dam. This study provides water resources managers and landuse planners with valuable information to improve future LULC policies and develop sub-watershed management strategies in the context of sustainable water resources and land-use planning and management. More research using high-resolution multi-temporal satellite imagery is needed to better estimate the dynamics of LULCs and plan for land-use, natural, and environmental resources in sub-watersheds.