Land use/land cover change impact on streamflow using Arc-SWAT model, in case of Fetam watershed, Abbay Basin, Ethiopia

Changes in Land Use Land Cover (LULC) are currently one of the greatest pressing issues facing the watershed, its hydrological properties of soil, and water management in catchment areas. One of the most important elements impacting streamflow in watersheds is LULC change. The main objective of this study was to evaluate the effect and future predication of LULC change on streamflow of the Fetam watershed by using Cellular Automata (CA)-Markov in IDRISI software. To analyze the impact of land use/cover change on streamflow, the Soil and Water Assessment Tool (SWAT) calibration and validation model was used. LULC map was developed by using Cellular Automata (CA)-Markov in IDRISI software, and the coverage of LULCs was including parameters of cropland, vegetation, grassland, Built-up area/Urban and water body. The findings of this study showed that the major challenges of land use/cove changes were rapid population increase, farming, and industrial activity. During the study period (2000–2020), most portions of the water body, vegetation, and grassland were changed into cropland and constructed by building. Cropland and construction areas increased by 15% and 46.95%, respectively, whereas water bodies, vegetation, and grassland decreased by 62.7%, 70.02%, and 38.1%, respectively. According to the forecasted results for the period of 2030–2040, cropland and built-up areas are increased, while vegetation, grassland, and water bodies were decreased. The SWAT model's calibration and validation performance was evaluated using the streamflow of the most sensitive parameters. For the years 2000–2004, and 2005–2012, the models were calibrated and validated, and the results showed good agreement between observed and simulated streamflow, with NSE and R2 values of 0.88 and 0.72 and 0.9 and 0.85, respectively. The results of this study indicated that the seasonal streamflow was decreased from 2000 to 2010 and 2010–2020 years during the dry and rainy seasons. In general, the impacts of land use/cover change on streamflow are significant considerations for planning and implementing water resource projects. In order to address the risks, effective land-use planning and climate-resilient water management strategies will be improved.


Introduction
The availability and competition for water will vary as a result of land use/land cover change (LULCC) and climate change, which are brought on by regional and global economic development (Dibaba et al. 2020;Disse et al. 2018;Mekonnen et al. 2022). High rainfall variability, rising food and water needs, frequent hydrological extremes, and rapid population growth all contribute to the environment's continued degradation by affecting the availability of various biophysical resources. In addition, to satisfy the demands of an expanding population, agricultural land expansion and intensification, urban area growth, and the requirement to extract wood products and fire fuel are all on the rise (Gashaw et al. 2018;Woldesenbet et al. 2018). Since LULCC is not only a local environmental problem but is also growing into a force of global significance, it will have an impact on important parts of how the earth's system functions when it aggregates globally (Disse et al. 2018). There is a chance that climate change will put further strain on the accessibility and availability of water. Interest in how these changes will affect the hydrological process will grow as a result of the projected negative environmental implications of these 1 3 111 Page 2 of 19 changes (Alemu et al. 2022;Dibaba et al. 2020;Mekonnen et al. 2022). The demand for fresh water is correlated with changes in lifestyle, economic growth, and population size. Water scarcity is a result of both excessive human demand and hydrological unpredictability (Alemu et al. 2022). Natural resources include land, water, air, energy sources, plants, soil, animals, and minerals, all of which are derived from nature. These resources are critical components of a catchment's of environment, providing socio-economic and ecological benefits to society (Chen et al. 2020;Hassaballah et al. 2017;Kidane et al. 2019;Tena et al. 2019;Yihun 2020).
The change in land use/cover in a given area is a result of environmental and socioeconomic issues, as well as human activity over time and space. Increases and decreases in population growth in the system, climate, financial development, and physical elements such as landscape, grade condition, and soil type all have an impact on Land Uses land Cover (LULC) changes : Moges et al. 2017Tewabe and Fentahun 2020).
Land use and land cover change is a major issue when discussing universal dynamics and how they respond to hydrologic features of water and soil administration in a catchment (Chang et al. 2018;Welde and Gebremariam 2017). Land cover/use and climate change has a important effect on the watershed's hydrological processes. The consequences of land use/land cover and climate change are predicted to increase in the future as a result of growing forest area clearing for farming and rising global temperatures (Hu et al. 2019;Hyandye et al. 2018). Increased deforestation, suburbanization, agriculture, and everyday actions of human-caused time-based and spatial changes in land use/ cover, which have influenced flow aquatic path patterns and aquatic stability (Bufebo and Elias 2020;Welde and Gebremariam 2017).
The technique in which the land has been put to use (land use/land cover) is defined, with a focus on the land's useful function in financial events. However, land-cover refers to the physical characteristics of the earth's surface, which are distributed in the form of vegetation, water, soil, and other physical characteristics of the land, including those created entirely by human activities such as continuous cultivation, deforestation, and overgrazing (Bekele 2019;Bufebo and Elias 2020;Saddique et al. 2020).
In Ethiopia, where 85% of the population works in agriculture and is dependent on current water supplies, analyzing and administering existing water resources is a big issue. The surface and rainwater flow model is a useful tool for investigating surface and rainwater systems and administering watersheds (Anmut et al. 2020;Bogale 2020).
Increased population and the constant demand for large amounts of food are the primary factors in Ethiopia, causing natural forests to be converted to urbanization, grasslands, and agricultural lands. As a result, they are the primary contributors to changes in land use and cover (Fentie et al. 2020;Tewabe and Fentahun 2020). The most common variations in land use/cover in our country include the changing of natural forests to farming lands, urban areas, overgrazing land, and other industrial and institutional areas (Fentie et al. 2020;Hyandye et al. 2018;Yohannes et al. 2018). As a consequence, agronomic fields have been enlarged at the expenditure of natural forests in order to deal with infrequent nutrition shortages for growing communities (Abere et al. 2020;Hyandye et al. 2018;Saddique et al. 2020). As a consequence, the landscapes of numerous watersheds must be enhanced, and the extent of the watershed has been increased to create good environments (Fentie et al. 2020;Yohannes et al. 2018).
Streamflow is a critical hydrological variable for water source preparation, design, and advance; this hydrological action is inextricably linked to LULC (Anmut et al. 2020;Bogale 2020). The tendency of overgrazing, farming and deforestation in the Fetam watershed has been expanding from time to time, due to the increase in farming land, which causes the land use/cover change of the catchment spaces. This constant land use/cover shift has changed the water stability of the catchment through altering the quantity and shapes of the mechanisms of streamflow, which include groundwater flow, surface runoff, and lateral flow. These effects have increased the degrees of water management challenges in study areas (Chen et al. 2020;Chimdessa et al. 2019;Mekonnen et al. 2022).
Changes in land use/cover have an impact on evapotranspiration, infiltration, erosion, and an interception in the runoff process (Abidin et al. 2017;Dibaba et al. 2020;Tena et al. 2019). In the current hydrological scenario, this causes a deal of issues. Land use/cover change is a type of increasing the percentage of the waterproof layer, which increases the amount of surface runoff and decreases the amount of water infiltrated into the ground. It also reduces the time of concentration, which causes many interruptions by producing higher volumes of runoff and lowering the amount of water infiltrated into the ground. These and other issues must be thoroughly examined in order to determine the effects of land use on various hydrological procedures in the research area (Bekele 2019;Bufebo and Elias 2020;Hasan et al. 2020;Owar Othow et al. 2017).
At this time, the study area is rapidly developing by infrastructures, such as rapid population growth with no gaps, urbanization, institutional and industrial facilities, increasing agricultural land, and requiring huge amounts of food production, all of which affect land use, land cover, and streamflow changes. The Fetam river watershed, landscape quantification, and economic production of the research area, as well as the country, are all affected by LULC change in some way. The Fetam watershed is one of the most important natural resources; it faces numerous issues as a result of urbanization, population growth, agricultural land, and industrialization, all of which contribute to rapid changes in land use/cover, and land degradation in the Fetam River. Land degradation is defined as the loss of top-fertilized soil due to flooding during the rainy season and a reduction in water supply for irrigation during the dry season; this is the most common concern in the Fetam watershed's streamflow. This study has been started specifically to examine into the watershed's hydrological reaction to the LULC change. LULC prediction was using integrated CA-Markov analytic methods. The evaluation of the isolated and combined LULC change was conducted using the hydrological modeling component of the Soil and Water Assessment Tool (SWAT). The major goal of this study was to determine the impact and future predication of land use and land cover change on streamflow in the Fetam watershed in Ethiopia's upper Blue Nile basin from 2000 to 2020 by using Cellular Automata (CA)-Markov in IDRISI software.

Description of the investigation area
The Fetam basin is one of the southern Upper Blue Nile basins. In the southern Upper Blue Nile basins, the catchment is located between 10o20′51.43″ N and 11o01′23.4″ N latitude and 37o52′02.02″ E and 36o55′2.56″ E longitude (Fig. 1). Fetam Rivers drain Awi Zone Banja, Injibara, Guagusa shikudad woredas and West Gojjam Zone Sekella, Bure Zuria, Bure Ketema, Womberima Woredas in Amhara National Regional State, 432 km north of Addis Ababa. The study area's watershed was estimated to be 757.53 km 2 (Fig. 1).

Climate of Fetam watershed
The year is divided into two seasons based on rainfall patterns: dry and rainy periods. The dry period started from November to April, while the wet period started from June to September. May and October are transitional months with frequent light rainfall. Through the wet period of June to September, the study area receives 70-90% of its yearly precipitation. The research area's mean yearly precipitation is 1896.97 mm, but there is a minor spatial variance within the area. The study area's monthly temperatures fluctuate from month to month, ranging from 9.1 to 29.9 °C. The monthly mean extreme temperature ranges from 23.8 °C in July to 31.1 °C in March, while the monthly mean minimum temperature ranges from 9.9 degrees Celsius in January and December to 14.2 degrees Celsius in May. March and April are the months with the greatest average temperatures. The monthly maximum temperature and monthly mean rainfall of five watershed stations are given in the graph below (Fig. 2).

Methods and data collection
Methods of this study are mainly focusing on assessing the impacts of land use/cover changes on stream flow of Fetam watershed by using Acr-SWAT model and collect the most important data.

SWAT model input data
Landsat imagery from the United States Geological Survey (USGS) commencing in 2000, 2010, and 2020, as well as a shapefile of the Fetam River basin, were used in this study. These data included a DEM with a resolution of 30 m by 30 m, soil types, meteorological, hydrological, and flow data, as well as land use and land cover data.

Soil data
Soil data are one of the main inputs for the SWAT model (Gashaw et al. 2020). The soil map for this study area was provided by Ethiopia's Ministry of Water, Irrigation, and Energy. The Fetam watershed highland is a flat plain watershed dominated by Nitisol and Vertisol, both of which have a main textural class of sandy clay and sandy loam. However, the predominant soil types found in the watershed's hills and lowlands are shallow Leptisols and cambisols.
The Fetam watershed is divided into five main soil classes according to FAO categorization. These are Cambisols, Leptosols, Luvisols Nitisols, and Vertisols (Fig. 3). For distinct layers of each soil type, the SWAT model needed physical and chemical parameters such as organic carbon contented, hydraulic conductivity, soil texture, available water content, and bulk density, some of which will be taken from the Harmonized World Soil Database (HWSD) (Roy and Inamdar 2019).

Hydrological data
Hydrological data, which includes daily streamflow or water discharge data from 2000 to 2020, are collected from the MoWIE's Hydrology Department. These data were used to calibration and validation models of hydrological components. Hydrological data are used to predict the effect of land administration methods such as reservoir administration, land use/cover changes, groundwater withdrawals, and water transmissions in complicated watersheds (Fentie et al. 2020). The model predicts important hydrological variables such as surface stormwater, maximum runoff, evapotranspiration, sediment output, and groundwater inflows for each HRUs unit.

Meteorological data
To simulate the river's hydrological circumstances, the SWAT model required climatological data. The National Meteorological Agency of Ethiopia provided the climatological data required for this investigation, which could be read from a measured data set using a weather generator model. Minimum and Maximum temperatures, rainfall, wind speed, comparative humidity, and sunlight are all essential climatological data, but rainfall was the most important one for this study, which covered four sites from 2000 to 2020. Injibara, Sekela, Tilil, Bure, and Shindi are the stations that are located in and near the study area, shown in (Fig. 4). The missing data are given a value of − 99, which is filled

Sekela
Bure Tilli Injibara Shindi by weather generators and prepared for input into the SWAT model.

Land use/covers data
Another major factor influencing evapotranspiration, runoff, and surface degradation in a watershed is LULC changes. The Fetam watershed land uses/cover changes from time to time due to various issues primarily changing by urbanization, agricultural practice, and irrigation expansion. Before the analysis can begin, LULC change studies frequently necessitate the increase in land cover components. The Landsat data image of the watershed which displays the land use land cover for three different years of 2000, 2010, and 2020 were downloaded from satellite images downloaded from USGS-GLOVIS (www. glovis. usgs. gov) used for land use change detection of the study area.

Methods of data analysis
Land use land cover analysis Image pre-processing and processing The aim of this investigation was to determine changes in LULC delivery in the Fetam watershed during a 20-year period from 2000 to 2020 using Land sat imageries from three groups. For the years 2000-2020, the land sat ETM+ was chosen. Layer stacking is performed on 7 images for landsat5, 8 images for land sat7, and 11 images for land sat8 before being composited into three groups (3, 2, 1) for Land sat 7 and Land sat 5, and a single layer for Land sat 8 (4, 3, 2). To reduce the seasonal change in forest shape and delivery through the year, the developed statistics dates are chosen as closely as likely in the similar yearly period of the developed years. All of the input satellite photos are false-color produced to distinguish and depict the surface landscapes clearly; these images cover the whole Fetam watershed.

LULC classification
This study employed digital image classification to determine land use and cover. The technique of arranging all of the input images into a small number of data classes was well-defined in digital image classification. The CA-Markov model combines the benefits of CA theory with the Markov model's capacity for temporal forecasting. As a result, the CA-Markov model performs better than the other techniques for simulating the temporal and spatial aspects of LULC changes (Yang et al. 2019). Two techniques Markov chain analysis and CA are used in the CA-Markov approach in IDRISI. In this investigation, the methods for projecting LULCCs with CA-Markov described by Pan et al. (2017) were applied within the IDRISI software framework. Creating suitability maps, calculating the transition matrix, and predicting the LULC map are the three steps of the procedure. The Markov chain model in IDRISI developed the LULC transi- tion matrix. Then, using a multi-criteria evaluation (MCE), the suitability maps of the probable driving elements of the LULC translations are developed for each LULC class. To create suitability maps, parameters including cropland, vegetation, grassland, Built-up area/Urban and water body were taken into consideration. Finally, the suitability maps and transition matrix are used to predict the LULC maps. Three LULC maps were developed for the LULCC study for the years 2000, 2010, and 2020. Future land use maps of LULC2030 and LULC2040 were predicted using the CA-Markov algorithm in the IDRISI software using the land use maps of 2010 and 2020 as a baseline. Model calibration and validation are essential for every model prediction. The model's validation outcomes will then determine how useful the model output is. Kappa coefficients derived from Eq. (1) were used to evaluate the CA-Markov model's efficacy in predicting land use maps (Dibaba et al. 2020;Pan et al. 2017). In this study, the 2020 LULC map was simulated using the LULC maps from 2000 and 2020. Then, using the Kappa index, the observed 2020 LULC map was contrasted with the simulated 2020 LULC map. (1) The location of meteorological and hydrological gauge and rainfall distribution systems where P o is the percentage of accurately simulated cells in the equation above, and P c is the anticipated proportional error due to chance between the observed and simulated map. If Kappa ≤ 0.5 indicates sparse agreement, 0.5 ≤ Kappa ≤ 0.75 indicates moderate agreement, 0.75 ≤ Kappa ≤ 1 indicates great agreement, and Kappa = 1 indicates complete agreement.
During the accuracy test of image classification of the years 2000 and 2010, GIS software was used to change the regulator points collected from Google Earth to shapefile, but 2020 collocate on the field. The following diagrams show the classification of land use/cover (Table 1).

Meteorological data analysis
Filling of missing data Data from a certain gauge station are unavailable, or illustrative precipitation is required at interest points. Using the average and normal ratio approach to use the rainfall in this study area, there are several procedures to fill in the missing data from that station. Initiate, calculate the percentage difference, to calculate the filling missing data at necessary station using this equation: where P x is A real average rainfall; P i is Rainfall measured at station 'i'; d i is the distance of sub-region for places 'i'; n is the numeral of neighboring locations. (2) Other missing meteorological data, such as temperature, relative humidity, wind speed, and sunlight hour, were given − 99 values and prepared for SWAT input by weather generators (Gashaw et al. 2020;Chen et al. 2020;Tenagashaw et al. 2022).
Homogeneity and trend test Time series homogeneity tests allow for the identification of a change as well as a time series. It is used to determine whether a series is homogeneous across time. Rainfall records at gauging stations are influenced by the station's location, the tool, and the method used to record and collect data, as well as the observation quality, and the time series may be inhomogeneous. As a result, the data utilized in the modeling of water resources and hydrology procedures would be statistically evaluated for reliability and quality. Rainbow and non-dimensional charts are commonly used to verify the homogeneity of rainfall time series data. As a consequence, while selecting an illustrative meteorological location for the research of areal precipitation approximation, examining the homogeneity of group stations is crucial. The similarity of the selecting gauging locations' monthly rainfall records was computed using the non-dimensional plot method. The following equation is used to test the homogeneity and trend of the average annual rainfall of all stations in and around the Fetam watershed from 2000 to 2020.
where P i.av is the over year's average monthly precipitation for the place I and P ave is the over year's average yearly precipitation for the site i. Because the meteorological and hydrological data are homogeneous. The homogeneity and trend test of rainfall data in five metrological places are shown in Fig. 5.

Checking the consistency of data
A continuous recording is one in which the qualities of the data in the record are retained over the period. To correct for measure consistency, the calculation of a result is employed rather than a missing value. Changes in the monitoring process, changes in the gauge's station, and changes in land use may necessitate modification when maintaining the gauge at its original site is impractical, and when damage occurs regularly. For adjusted inconsistency data, the approach of double mass-curve analysis was employed to validate the consistency of precipitation in a gauge record. Equation 4 should be used to correct it.
(3) PX′ represents corrected precipitation at station x, PX represents initially recorded precipitation at station x, M′ represents the corrected slope of the double mass curve, and M represents the actual gradient of the double mass curve shown in Fig. 6.  Fig. 6 Consistency checking of the five rainfall stations within and around the catchment hydrological sequence that regulates number of nutrients, water, pesticides, and sediment loading into each sub-primary basin's channel. In the land stage of the hydrological rotation, cover evapotranspiration, storage, infiltration, surface runoff, sideways subsurface movement, redistribution, branch stations, and reoccurrence flow were all modeled. The hydrologic cycle's routing stage, defined as the transport of water, nutrients, sediments, and organic compounds through the watershed's channel network to the outflow, was the second division. In the SWAT model, the water balance equation was employed to mimic the hydrological cycle (Assfaw 2020;Saddique et al. 2020). This equation is shown below:

Hydrological component of SWAT
where SW t is the final soil moisture at a time 't' (mm), SW o is the initial soil moisture (mm), t is the time (days), R day is the quantity of rainfall on day I (mm), Q surf is the quantity of surface runoff on day I (mm), E a is the quantity of evaporation and transpiration on day I (mm), W seep is the amount of water seeping into the unsaturated zone from the surface soil on day I (mm).
The weather producer in the SWAT model can be used to simulate or measure precipitation inputs for hydrologic estimation. The delineation of watersheds and the generation of streamflow networks are the first steps in the SWAT model's modeling of the hydrological approach. The SWAT modeling format was completed in time for the first simulation, which was used to evaluate the model. Before assessing the model's performance in the catchment using the SWAT simulation results (Dibaba et al. 2020;Tena et al. 2019).

Sensitivity analysis, calibration and validation
The sensitivity analysis was a crucial element in Arc-SWAT, as it allowed for two different types of performance studies. The first type of study classifies the impact of changing a parameter value on a simulated production quantity, such as normal streamflow, using only one modeled data. The next way of examination uses observed data to calculate the total goodness of fit between the modeled and measured time series (Welde and Gebremariam 2017). The first assessment is used to classify factors that enhance a certain model technique, while the second analysis is used to classify parameters that are influenced by the study watershed's attributes. As a consequence, vulnerability assessment is a strategy for distinguishing the greatest vulnerability characteristics that have a significant impact on the calibration and validation model, as well as defining how the model output differs throughout a range of assumed input variables (Abere et al. 2020;Alemayehu et al. 2019).
Flow simulation for 20 years of recording times is performed after a thorough preparation of the relevant input data into the Arc-SWAT model. The six-year flow data are used as a warm-up period for the simulation, which is used for the evaluation of hydrologic variables and model calibration. The simulation findings cannot be directly used for further study; instead, stream flow must be examined using evaluation, calibration, and validation models (Bufebo and Elias 2020; Tewabe and Fentahun 2020). The method of evaluating model parameters in collaboration with the model forecasted by experimental data for the same scenario was known as hydrological model calibration (Dibaba et al. 2020;Tenagashaw et al. 2022). Data from the time series of discharge at the watershed's outlet are used to calibrate and validate the SWAT model. Measured streamflow data from the years 2000 to 2005 were used to calibrate the model.
Validation is a technique for putting the calibrated model to the test using the self-determining dataset without any additional constraint modifications. The calibrated model was built in opposition to a self-governing set of measured data in order to evaluate the efficiency of future prospective administration procedures (Alemayehu et al. 2019;Yihun 2020). Validation models are based on eight years of streamflow data, covering the years 2006 to 2012.

SWAT CUP
For SWAT alone, the SWAT CUP interface was developed. Using this universal interface, SWAT can easily be linked with any calibration or sensitivity software. The program's connection of the GLUE, Parasol, SUFI2, and MCMC procedures to SWAT is an example of this. In this investigation, sequential uncertainty fits were favored (SUFI2). Automated model calibration involves running the model, altering the unknown model parameters carefully, and then extracting the necessary outputs (corresponding to measured data) from the model output files. An interface's main function is to link a calibration program's input and output to the model.

Model performance evaluation
To compare the model simulation of calibration and validation outputs to the experiential data, model performance estimation was necessary. A variety of strategies could be used to estimate the model performance during the calibration and validation stages. To measure the goodness-of-fit of the model method, two techniques were employed for these investigations: coefficient of determination (R 2 ) and Nash-Sutcliffe Simulation Efficiency (NSE). The R 2 represents how well the model explains the modification in the measured data. The determining factor was the size of the linear relation between the simulated and experimental standards. R 2 is a measure for model quality that goes from 0 to 1, with higher values suggesting lower error variance.
In contrast to the observed data variance, the NSE evaluates the normalized relative amount of remaining variation. The NSE ranges from 1.0 to − ∞, with 1.0 being the best fit. In this technique, the variance between experimental and model results is squared, and low values in time series, such as base flow or lateral subsurface flow, have a minor effect on total NSE. Furthermore, in low flow scenarios, the Nash-Sutcliffe efficiency coefficient is typically unaffected by over or under forecasting (Hassaballah et al. 2017). Within the investigation phase, this issue can be detected by cooperative forecast and experiential values.
In general, model simulations for shorter time phases are weaker than for longer time phases (Guzha 2018). The NSE statistics performance ratings listed above is for monthly time periods and must be improved for a daily time period to be appropriate in this research. Equations (6) and (7) are used to calculate the R 2 and NSE values, respectively.
where R 2 represents the Model approach coefficient of determination, X i represents the measured value (m 3 /s), X av represents the average measured value (m 3 /s), Y i represents the simulated value (m 3 /s), Y av represents the average simulated value (m 3 /s), and NSE represents Nash-Sutcliffe simulation efficiency.

Sensitive parameter
The Ribb watershed was divided into 149 HRUs and 25 sub-basins for this study. With the use of the SWAT CUP (Calibration and Uncertainty Program) model, the sensitive parameters were found (Table 1). In the Upper Ribb watershed CN2, SOL_K, CH_K2 SLP, and SOL_BD were found 1st to 5th sensitive flow parameters, respectively (Table 1). On the other hand, found ESCO, SOIL-AWC, EPCO, CN2, and ALPHA_BF as the first five subtle flow parameters. T-STAT and p values had been used to assess the parameter's level of sensitivity and significance. High absolute T-STAT values for a parameter are regarded as very sensitive, and p values near 0 as highly significant ( Table 2).

Analysis of changes in land use/cover
The simulation of LULC 2020 using the CA-  built-up areas, as well as main components of farmed land or croplands in all parts of the watershed. The overall coverage of cropland, vegetation, grassland, water bodies, and built-up area on the LULC map of 2000 ( Fig. 7 and Table 3) was about 59.936%, 5.32%, 23.19%, 0.069%, and 11.49%, respectively, from the whole part of Fetam watershed. Due to the high population density in the Fetam watershed, this increases the amount of new settlement and farming land available. This is most likely due to deforestation actions that have occurred in the name of farming, urbanization, and rural settlement since the latter was recognized in recent years.
According to land use, land cover changes in 2010 ( Fig. 8 and Table 4) showed that cropland was covered by 79.607%, followed by Grassland, Forest land, building area, and water bodies, which were covered by 15.738%, 2.62%, 1.813%, and 0.222%, respectively. Due to the expansion of farming regions from time to time, the cropland coverage was a very large area during this year.
According to the LULC map of 2020, crops covered roughly 76.918% of the overall area of the Fetam watershed ( Fig. 9 and Table 5). Due to high growth of population in the Fetam watershed, this increases the amount of new settlement and farming land available. Other water bodies, vegetation, grassland, and buildup area, on the other hand, were covered to varying degrees: 0.207%, 4.082%, 15.792%, and 3.001%, respectively. This is due to the annual destruction of forests for farming and the establishment of new urban areas and rural settlements surrounding the Fetam watershed. In the Fetam watershed, cropland increased by 15% over a 20-year period, whereas vegetation land decreased by 70.02%.
Using the 2020 land use map as a baseline map, the LULC 2030 and LULC 2040 future predictions were carried out. Using LULC data between 2010 and 2020, the transition matrix and transfer probability matrix were generated.
According to the forecasted results for the period of 2020-2030, as shown in (Table 6), cropland and built-up areas are increased, while vegetation, grassland, and water bodies were decreased. Cropland, urban built-up and vegetation were expected to continue increasing from 2030 until 2040. In contrast to the outcome of the previous prediction, the rate of cropland and built-up expansion is, however, slower. Due to the farms' reliance on irrigation and the restricted amount of land available for built-up expansion,  this may be feasible. Additionally, the decrease in agricultural expansions may be attributed to the extension of range land through ongoing conservation efforts (soil and water conservation efforts have already started). This might also be influenced by the rapid urbanization. Another explanation might be the expansion of the grasslands brought on by the planned afforestation.

Land use/cover changes detection
The ERDAS Software and GIS were used to detect land use/ cover changes. Table 7 compares the detection of changes in images from 2000 to 2010, 2010 to 2020, and 2000 to 2020. According to above table, the gradual change of LULC, particularly the incremental of cropland by percentages of 13.671% between (2000-2010) and decreased by 3.089% between (2010-2020), and the built-up area was increased by 0. 33% and 4.835% between (2000-2010) and (2010-2020), respectively, while vegetation, grassland, and water bodies were decreased by 4.236%. 9.37%, and 0.315% between (2000-2010) years and 1.462%, 0.146% and 0.138% between (2010-2020) years, respectively, the Fetam river streamflow decreased due to high percentage reductions of water body, vegetation, and grassland by 62.74%, 68.502%, and 38.097%, respectively, as well as increases in cropland and built-up area of 17.655% and 88.5024% in the watershed over 20 years.
According to Fig. 10 water bodies, vegetation and grassland were decreased for both periods of 2000-2010 and 2010-2020 and also 2000-2020 periods. However, cropland is increased and decreased during periods of 2000-2010 and 2010-2020, respectively, but increased starting from periods of 2000-2020. The buildup area was increased throughout three periods, which are 2000-2010, 2010-2020 and 2000-2020 periods.

Calibration and validation model of sensitive analysis
Sensitivity examination of Fetam streamflow for the subbasins was performed using day-to-day observed and

Model of stream flow calibration and validation
The According to Hassaballah et al. (2017), the NSE and R 2 for the daily calibration and validation of the three precipitation products ranged from 0.4 to 0.80 and 0.50 to 0.80, respectively. According to Welde and Gebremariam (2017), the R 2 and NSE for the monthly average streamflow values of calibration and validation are 0.88 and 0.85 and 0.85 and 0.84, respectively. Comparison of monthly average streamflow between actual and simulated values of calibration and validation models is shown in Figs. 11, 12, 13 and 14. During the calibration and validation periods, strong correlation values of 0.72-0.9 were observed and simulated, showing that the NSE and R 2 Values Of The Model, Which Were accurate. However, groundwater percolation is responsible for some of the flow loss. The findings of calibrated and validated streamflow models show that simulated and observed streamflow is in good agreement. As a result of the streamflow data (Table 8), the SWAT model is a good estimate for the Fetam River streamflow (Figs. 12, 14).

Land use/covers change impacts
Satellite images from three separate years showed an effect on the Fetam River's streamflow, which was utilized to construct three different SWAT model simulations, according to the findings of LULC change. The results of the simulated streamflow in the years 2000, 2010, and 2020 were not the same as the results of the whole streamflow computed to its observed streamflow. This indicates that the LULC modification had an impact on the Fetam River's monthly and seasonal streamflow.

Assessments of the effects of changes in land use/ cover on streamflow
Land use/cover change detection of three separate years' satellite image results showed an effect on the streamflow. While calibrating and validating the model using three land use/cover maps, the simulation was run with all   Due to increased agricultural and other land use/cover change, seasonal streamflow in the Fetam watershed decreases in the dry period compared to the wet period, resulting in higher flow during the wet periods. The results of this study indicated that the seasonal streamflow decreased by 5.54 m 3 /s, and 0.38 m 3 /s and 0.23 m 3 /s and 12.28 m 3 /s, respectively, at 2000-2010 and 2010-2020 during the dry and rainy seasons. LULC variations of simulated mean monthly streamflow of Fetam watershed during study periods were presented in 2000, 2010, and 2020 (Fig. 15). For both the calibration and validation phases of land use/ cover, these results showed a relatively increased streamflow from 2000 to 2020.
Changes in land use/cover have an effect on the Fetam River's average monthly and seasonal streamflow throughout the year, as shown below (Table 9).

Impacts of LULC change on groundwater and surface runoff
There are two LULC change scenarios, the first for LULCC from 2020 to 2030 and the second for LULCC from 2020 to 2040 scenarios with the base land use map of 2020, were used to evaluate the hydrological responses of the watershed    to LULC change. Changes in land use/cover had an impact on components of streamflow in the Fetam watershed due to the effects of groundwater discharges and surface runoff throughout the year. The watershed's hydrological responses were taken into account in terms of the processes affecting yearly surface runoff; annual groundwater, total water yield, and potential evapotranspiration were shown in (Table 10). The mean annual surface water runoff and total water yield will rise as a result of the significant anticipated changes in land use between 2020 and 2030 and 2020 to 2040, whereas groundwater recharge and potential evapotranspiration will fall. However, the frequency of LULC changes was correlated with the rates of rise and fall. For instance, it was expected that the first scenario would have a higher rate of agricultural development and a lower rate of forest decrease. As a result, compared to the baseline scenario, the rate of increase in surface flow was slower, rising from 2030 to 2040.
In the first and second scenarios, the surface runoff and water yield were increased by 6.85% and 2.05% and 8.15% and 2.85%, respectively. In the first and second scenarios, the groundwater level and potential evapotranspiration were decreased by 1.5% and 1.45% and 0.85% and 1.76%, respectively, Due to the higher decrease in forest and shrub land, the declining rate is larger for the changes from 2020 to 2030. This is a good example of how changing the land cover can affect groundwater management; for example, if the amount of shrub land grows, groundwater recharge will also increase. Evapotranspiration decreased as a result of the loss of forestland, which caused evapotranspiration to diminish. Furthermore, a lack of soil moisture brought on by increasing runoff may also be a factor in the drop in evapotranspiration. A related study found that as a result of less rainfall, PET decreases when soil moisture is low (Tenagashaw et al. 2022;Woldesenbet et al. 2018).
With the greatest growth of agricultural areas, urban expansion, and decline of forested and shrub land, there is a corresponding greatest increase in surface flow and greatest decrease in groundwater recharge. This shows that converting large areas of LULC to intensive agriculture and settlements will decrease the soil's ability to absorb water, resulting in a significant increase in surface runoff. Groundwater flow declines as a result of decreased soil water infiltration. Additionally, the growing population's increased need for groundwater consumption could worsen the groundwater crisis by increasing groundwater withdrawal. According to recent research by Dibaba et al. (2020) in the Finchaa catchment, the high variations in LULC caused the catchment's springs to run dry and the volume of water to decrease. As a result, flow discharges during the dry season were lower than base flow, while flow discharge during the rainy months was higher than base flow. Changes in land use/cover had a significant impact on runoff output, infiltration rate, and soil water-holding capacity in watersheds, according to this study finding.

Conclusions
This study focuses on the effects of land use/cover change and future prediction of LULC change on streamflow of the Fetam watershed by using Cellular Automata (CA)-Markov in IDRISI software. There are five major parameters were identified for this study of LULC changes, and these are water bodies, vegetation, grassland, cropland, and built-up area. The results of this study indicated that the cropland and built-up area were increased by 15% and 46.95%, respectively, over a 20-years period (2000-2020); however, water bodies, vegetative land, and grassland decreased by 62.74%, 70.02%, and 38.1%, respectively. The results of this study indicated that the future forecasted results of cropland and built-up areas are increased, while vegetation, grassland, and water bodies were decreased for the period of 2020-2030 and 2030-2040. The major causes of this LULC changes are rapid population growths around watersheds, as well as pressures to expand farming land and develop urbanization areas. Performances evaluation of calibration and validation of streamflow were showed good agreement between observed and simulated values of NSE and R 2 , which are 88% and 72% and 90% and 85%, respectively. This study showed that monthly average streamflow decreased from 2000 to 2010, and 2010 to 2020, by 2.885 m 3 /s and 6.33 m 3 /s, respectively. The results of this study indicated that the seasonal streamflow decreased by 5.54 m 3 /s, and 0.38 m 3 /s and 0.23 m 3 /s and 12.28 m 3 /s, respectively, at 2000-2010 and 2010-2020 during the dry and rainy seasons. Impacts of LULC changes in the first and second scenarios, the surface runoff and water yield were increased by 6.85% and 2.05% and 8.15% and 2.85%, respectively. But the groundwater level and potential evapotranspiration were decreased by 1.5% and 1.45% and 0.85% and 1.76%, respectively. Due to high amounts of sediment loads and surface runoff during the rainy season, as well as deforestation and increased cultivated land. Therefore, the effect of land use/cover change on streamflow is a critical issue that affects volumes of streamflow discharges year to year. In general, the impacts of land use/cover change on streamflow are significant considerations for planning and