Watershed Hydrological Responses to Land Cover Changes at Muger Watershed, Upper Blue Nile River Basin, Ethiopia

: Changes in land cover (LC) are the major factors inﬂuencing the hydrological processes within a watershed. Understanding the impacts of LC on watershed hydrology is crucial for planning and predicting land resource utilization, water resources, and sustaining hydrological balance. This study assesses the hydrological responses of LC changes in the Muger watershed located in the Upper Blue Nile River Basin (UBNRB) from 1986 to 2020. We used the Soil and Water Assessment Tool (SWAT) hydrological model to investigate the effects of LC on the hydrological process. The simulations were driven by several datasets, such as watershed elevations, mean climatology, hydrology and soil datasets, and LC satellite maps for three time periods (i.e., satellite imagery taken in 1986, 2003, and 2020). We found that the key LC changes that affected hydrological parameters in the Muger watershed are changes in cultivation land, forest land, and settlement. The expansion of cultivation land and shrinkage of forest and shrub lands triggered surface runoff and a reduction in groundwater between 1986 and 2003. Additionally, settlement was identiﬁed as the primary factor contributing to increases in evapotranspiration (ET) and surface runoff. The LC changes that occurred between 1986 and 2020 reduced the average annual, wet season, and dry season streamﬂow. Between 2003 and 2020, surface runoff decreased by 3.71% due to the effect of land landscape restoration interventions. The outcome of the study can assist decision-makers and planners in preparing adaptable strategies under changing LC conditions within a watershed.


Introduction
Land cover (LC) variations considerably modify a watershed, which, depending on the type of environmental variable, has considerable effects on hydrological parameters by changing the quality, amount, dissemination, and timescale of streamflows, ultimately impacting the water resource management and operations [1]. Generally, LC changes influence the hydrological parameters of the basin by partitioning rainfall routes into the surface and subsurface runoff, as well as through ET. Similarly, LC variations have Longitudes. The basin's yearly rainfall ranges between 800 mm and 2000 mm. The Muger catchment is administratively located in Oromia regional state, approximately 32 km west of Addis Ababa, which is the capital city of the country. The catchment is a sub-watershed of the UBNRB covering 7246 km 2 and joins the UBNRB at 37.93 • E and 9.92 • N. The Muger River is one of the largest tributaries that join the upper Blue Nile from the basin's southeast [30,31].
The catchment is found geographically between 37 • 54 42.57 E and 39 • 01 24.16 E, longitudes and between 9 • 04 55.54 N and 9 • 58 04.64 N latitudes (Figure 1). The altitude of the Muger watershed varies between 953 and 3550 m above sea level (m.a.s.l.). The southern and eastern highlands of the sub-watershed are higher in altitude, ranging from 2600 m.a.s.l. to more than 3550 m.a.s.l. The Muger River's lowlands are lower in elevation, at less than 1700 m.a.s.l. The mean yearly rainfall of the watershed ranges from 833 mm to 1326 mm. In the highlands of the sub-basin, relatively high rainfall has been observed. Around the river and in the lowlands, the mean yearly rainfall is between 833 mm and 1000 mm.

Study Area Description
The UBNRB is the primary stream of the Nile basin, located in western and central Ethiopia between 7°45′ N and 12°45′ N latitudes and between 34°05′ E and 39°45′ E Longitudes. The basin's yearly rainfall ranges between 800 mm and 2000 mm. The Muger catchment is administratively located in Oromia regional state, approximately 32 km west of Addis Ababa, which is the capital city of the country. The catchment is a sub-watershed of the UBNRB covering 7246 km 2 and joins the UBNRB at 37.93°E and 9.92°N. The Muger River is one of the largest tributaries that join the upper Blue Nile from the basin's southeast [30,31].
The catchment is found geographically between 37°54′42.57″ E and 39°01'24.16″ E, longitudes and between 9°04'55.54″ N and 9°58'04.64″ N latitudes (Figure 1). The altitude of the Muger watershed varies between 953 and 3550 m above sea level (m.a.s.l.). The southern and eastern highlands of the sub-watershed are higher in altitude, ranging from 2600 m.a.s.l. to more than 3550 m.a.s.l. The Muger River's lowlands are lower in elevation, at less than 1700 m.a.s.l. The mean yearly rainfall of the watershed ranges from 833 mm to 1326 mm. In the highlands of the sub-basin, relatively high rainfall has been observed. Around the river and in the lowlands, the mean yearly rainfall is between 833 mm and 1000 mm. The watershed's yearly maximum and minimum temperature range between 16 and 31.5 °C and between 3 and 16.5 °C, respectively. Lowlands that are hot to warm and moist are found in the watershed's northwest. Rendzic Leptosols, Eutric Leptosols, Eutric Vertisols, and Chromic Luvisols are the predominant soil types in the basin. In the Muger The watershed's yearly maximum and minimum temperature range between 16 and 31.5 • C and between 3 and 16.5 • C, respectively. Lowlands that are hot to warm and moist are found in the watershed's northwest. Rendzic Leptosols, Eutric Leptosols, Eutric Vertisols, and Chromic Luvisols are the predominant soil types in the basin. In the Muger watershed, small-scale subsistence farming is the main economic activity and the source of sustenance. Cultivation is very intensive in the rainy season.

Input Data
The inputs for the SWAT model are spatial data (i.e., digital elevation model land cover at different periods, soil) and temporal data (weather and streamflow data) from the catchment to simulate the hydrological parameters.

Topographic Factors
The digital elevation model (DEM) is among the crucial data required by the SWAT model. It is used to derive stream networks, watershed boundaries, and terrain slopes, among others. A DEM with a resolution of 30 m × 30 m was obtained from the Ministry of Water, Irrigation, and Energy (MoWIE) of Ethiopia. It was used to assess the hydrological parameters of the catchment. Five classes of slope were used in this study (Table 1 and Figure 2). The majority of the land in the catchment (32.48%) has a slope class between 15 and 30%, and the remaining land has the rest of the slope classes (Table 1).

Land Cover
LC is input data used to determine the characteristics of a catchment and its impact on hydrological responses. The LC data used in this study are based on the study conducted by Teshome et al. [32] (Figure 3, Table 2). It is a crucial element influencing surface

Land Cover
LC is input data used to determine the characteristics of a catchment and its impact on hydrological responses. The LC data used in this study are based on the study conducted by Teshome et al. [32] (Figure 3, Table 2). It is a crucial element influencing surface erosion, runoff, and ET, as well as describing the HRU in a given watershed [33,34].

Soils
For the SWAT model to simulate hydrological components for the different soil types, data on the basic physicochemical soil properties are essential. A soil map was obtained from the MoWIE, Ethiopia (Figure 4). The basic physicochemical characteristics of soils in the study area were obtained from FAO soil categorization. Nine types of soil were found in the watershed (Table 3; Figure 4). Rendzic Leptosols, Eutric Vertisols, and Humic Nitisols are the dominant soil types in the study watershed.   For the SWAT model to simulate hydrological components for the different soil types, data on the basic physicochemical soil properties are essential. A soil map was obtained from the MoWIE, Ethiopia ( Figure 4). The basic physicochemical characteristics of soils in the study area were obtained from FAO soil categorization. Nine types of soil were found in the watershed (Table 3; Figure 4). Rendzic Leptosols, Eutric Vertisols, and Humic Nitisols are the dominant soil types in the study watershed.

Meteorological Data
The daily meteorological data required for the SWAT model include rainfall, tempe ature, relative humidity, solar radiation, and wind speed. Nine meteorological station (Table 4) were selected based on their quantity, quality, duration, consistency, homogen ity, and spatial distribution. The Ethiopian National Meteorological Agency provide data for the years from 1983 through 2020. Furthermore, the missing weather data for s lected stations were filled with data from different nearby stations using both the metho of arithmetic mean and the normal ratio. The double mass curve was used to analyze th homogeneity and consistency of the data, and it was found that the data were consisten

Meteorological Data
The daily meteorological data required for the SWAT model include rainfall, temperature, relative humidity, solar radiation, and wind speed. Nine meteorological stations (Table 4) were selected based on their quantity, quality, duration, consistency, homogeneity, and spatial distribution. The Ethiopian National Meteorological Agency provided data for the years from 1983 through 2020. Furthermore, the missing weather data for selected stations were filled with data from different nearby stations using both the method of arithmetic mean and the normal ratio. The double mass curve was used to analyze the homogeneity and consistency of the data, and it was found that the data were consistent.

Hydrological Data
The daily recorded hydrological (i.e., streamflow) data from the gauging stations for the period from 1983 to 2009 were used for the SWAT model calibration and validation. The streamflow data of the Muger watershed were obtained from MoWIE, Ethiopia. These data were collected from five gauging stations in the catchment, i.e., Muger, Deneba, Sibilu, Gerbi, and Gorfo. These are the only gauging stations in the catchment. Therefore, the streamflow data at the gauging stations were interpolated to predict streamflow values at ungauged locations. The description of the gauging stations is shown in Table 5. The general conceptual framework used for hydrological modeling and assessment of the impacts of land cover change on hydrological processes in the Muger watershed is shown in Figure 5.

Hydrological Modeling
The SWAT model has been used to evaluate the impacts of LC change on hydrological components in the Muger watershed [35]. The model has been used to examine hydrological parameters in large and small catchments around the world, including Ethiopia [12,15,16,35]. The SWAT model is a long-term, semi-distributed, continuous, deterministic, and effective hydrological model [12]. The fundamental attributes of the SWAT model for a given watershed are hydrology, weather, erosion/sedimentation, nutrients, pesticides, plant growth, agricultural management, routing of ponds and reservoirs, and channel routing [36]. In this study, the basin was divided into a number of sub-basins, which were then further divided into smaller areas known as HRUs (hydrological response units). In a subwatershed, HRUs are parts of land that have non-spatial units and a unique combination of homogenous soil, land cover, slope, and management attributes [12,13]. The model uses DEM, soil type, LC, and slope to split the Muger catchment into 33 sub-catchments, which were then divided into 434, 435, and 451 HRUs for 1986, 2003, and 2020 LC data, respectively. The HRU has been defined based on a threshold of 10% for land cover, 20% for soil, and 10% for slope [35]. The general conceptual framework used for hydrological modeling and asse of the impacts of land cover change on hydrological processes in the Muger water shown in Figure 5.

Hydrological Modeling
The SWAT model has been used to evaluate the impacts of LC change on hydro components in the Muger watershed [35]. The model has been used to examine hy ical parameters in large and small catchments around the world, including E

Sensitivity Analysis
The model's statistics, including the t-statistic and p-value, provide an indication of the sensitivity of the model variables [37][38][39][40]. A larger t-stat in absolute values and a lower p-value suggest more sensitive factors [37].

Calibration and Validation
The process of determining model components by comparing simulated parameters with observed data is known as hydrological model calibration [41]. For calibration, the most delicate hydrologic parameters were utilized, and their values were iteratively modified within acceptable upper and lower ranges until satisfactory concordance among measured and simulated streamflow was achieved [12,40]. In this research, calibration and validation were performed using monthly observed streamflow data for 27 years

Performance Evaluation
An evaluation of the model performance is required to investigate the modeling representation process to the real biophysical circumstances [42,43]. The SWAT model performance for the goodness-of-fit test was carried out according to Moriasi et al. [41]. The various indices utilized to assess the model performance between observed and simulated values were NSE, PBIAS, and R 2 . The NSE value ranges from −∞ to 1, with a higher value indicating the model's good performance [38,41,44]. The percent bias (PBIAS) measures how much the observed variable is underestimated or overestimated [41,45]. The coefficient of determination (R 2 ), which ranges between 0 and 1, assesses the consistency of simulated and observed data. Equations (1)-(3) were used to determine R 2 , NSE, PBIAS, and RSR, respectively. Table 6 shows the overall model performance ratings and their key attributes.
where Qsi is the simulated stream, Qoi: is the observed streamflow value, − Qo is the average observed streamflow, − Qs is the mean simulated streamflow.

Model Application
The simulated results were used to analyze the impacts of LC change on hydrologic parameters at the catchment and sub-catchment scales, as well as to analyze the contribution of changes within individual LC categories. Therefore, to analyze the effects of LC on the hydrological parameters of the Muger catchment, the SWAT model was calibrated and validated depending on each of the classified LC inputs (1986,2003, and 2020) while maintaining the DEM, meteorological variables, and soil data constant. A number of studies have used this approach in various parts of the world [19,24,[45][46][47].

Sensitivity Analysis
The analysis of sensitivity of the Muger watershed was assessed using the SWAT-CUP at five gauging stations (Muger, Sibilu, Gorfo, Gerbi, and Deneba). As shown in Table 7, the 16 most influential sensitive parameters for the watershed were identified and graded based on how sensitive they are, then prioritized in the calibration process. In order to obtain good concordance between simulated and observed streamflow at each gauging station, the hydrological component values were iteratively modified within the permitted limitations [37][38][39]48]. However, the remaining parameters vary among stations. As a result, most delicate hydrological components were utilized in the processes of model calibration and validation [49].

SWAT Model Calibration and Validation
Niraula et al. [50] suggested that the hydrological model needs to be calibrated before determining the effects of LC change. Five gauging stations were used in the Muger watershed to calibrate and validate the model. The model calibration was conducted for the measured years from 1986 to 1999 and validated using measured data from 2000 to 2009 using monthly data at each gauging station. A model spin-up period was considered to have occurred during the first three years of simulation (1983)(1984)(1985).
The statistical performance indicators presented in Table 8 demonstrate a good consistency between the simulated and observed streamflow data during the calibrated and validated periods [19,38,41,45]. A graphical comparison of the average monthly simulated and measured streamflow at all gauging stations revealed similar patterns ( Figure 6).   The relative width of the 95PPU was slightly smaller, but it was still within the acceptable range predicted by the model [51]. The model performance is enough to simulate the streamflow in the Muger watershed, as evidenced by the simulated and observed hydrograph, which shows that the model successfully captured the trends of hydrological variations both during the validation and calibration phases. The simulation results revealed that during the calibration and validation stages, the SWAT model successfully replicated the observed streamflow in the investigated watersheds. These results were in line with earlier research conducted in the area [15].

Land Cover Change Impacts on Hydrological Parameters in the Watershed
The impact of LC change on the hydrological processes of the Muger catchment was evaluated using the LC classes of different periods. The effects of LC modification on the watershed's hydrological parameters were simulated using the calibrated SWAT model while taking into account three alternative LC periods (i.e., 1986, 2003, and 2020). The impacts of LC modification on watershed hydrology were independently examined with LC change data while the remaining model calibration variables and other SWAT inputs remained constant. Therefore, LC was the sole factor considered responsible for the The evaluation of the model during the calibration period showed that the SWAT model somewhat underestimated the mean streamflow at Sibilu station by 8.9%. On the other hand, the mean streamflow at Muger, Gorfo, Gerbi, and Deneba stations was overestimated by 4.5%, 20.3%, 13%, and 5.9%, respectively. The model slightly overestimated the mean streamflow at the Sibilu, Gorfo, and Gerbi stations, respectively, by 7.8%, 22.8%, and 17.7% during the validation period, while slightly underestimating it at the Muger and Deneba stations ( Figure 6). This could be due to the input data quality, uncertainties, and a restriction of the SCS in the SWAT simulation because the model takes into account daily average rainfall depth rather than intensity and duration [48,51]. The sample size for streamflow data could contribute to the observed limitations in this study. Similarly, the SWAT model has the weakness of underestimating peak flow and sediment during extreme events.
The relative width of the 95PPU was slightly smaller, but it was still within the acceptable range predicted by the model [52]. The model performance is enough to simulate the streamflow in the Muger watershed, as evidenced by the simulated and observed hydrograph, which shows that the model successfully captured the trends of hydrologi-cal variations both during the validation and calibration phases. The simulation results revealed that during the calibration and validation stages, the SWAT model successfully replicated the observed streamflow in the investigated watersheds. These results were in line with earlier research conducted in the area [15].

Land Cover Change Impacts on Hydrological Parameters in the Watershed
The impact of LC change on the hydrological processes of the Muger catchment was evaluated using the LC classes of different periods. The effects of LC modification on the watershed's hydrological parameters were simulated using the calibrated SWAT model while taking into account three alternative LC periods (i. e., 1986, 2003, and 2020). The impacts of LC modification on watershed hydrology were independently examined with LC change data while the remaining model calibration variables and other SWAT inputs remained constant. Therefore, LC was the sole factor considered responsible for the variations in hydrological parameters. According to the hydrological simulation results, changes in LC had a considerable effect on streamflow and hydrological parameters. The spatial and temporal effects of LC change on the hydrological components of the Muger catchment are presented in Table 9. The simulation of monthly and seasonal hydrological parameters between 1986 and 2020, based on three historical LC change data, demonstrated great volatility. The mean annual surface runoff increased while groundwater recharge and water yield varied as a consequence of the variations in LC from 1986 to 2003 and from 1986 to 2020. The rate of change of hydrological components is depicted in Figure 7. Generally, the effects of the LC change were associated with the increase in the area of cropland together with a decline in vegetation cover. The results indicated that the expansion of cultivation land and the decline of forest cover from 1986 to 2003 were the key contributors to the increment in average annual surface runoff (Tables 2 and 9). The amount of surface runoff changed from 319.91 mm in 1986 to 333.55 mm in 2003 as a consequence of the expansion of cultivated land and settlements and the decline in vegetation.
As the amount of cultivated land increases, soil infiltration capacity decreases due to the compaction of soil, which decreases lateral flow and increases surface runoff. In contrast, the surface runoff declined from 333.55 mm in 2003 to 321.61 mm in 2020 due to the gradual increment of shrubland as a result of the implementation of current soil conservation practices in the watershed. A decrease in lateral flow occurred during the 1986-2020 LC change, which caused a reduction in water yield despite an increase in surface runoff. This could be attributed to the poor infiltration of non-vegetated areas. Similar relationships between an increase in surface runoff and a commensurate decline in lateral flow were found by Baker and Miller [24]. In general, surface runoff increased by 4.09% in the first period (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003) and decreased by 3.71% in the second period (2003-2020). Karamage et al. [22] found out that the increase in runoff could be caused by agricultural land expansion and urban sprawl at the expense of forest areas. Similar findings were presented by Teshome et al. [33], who found that soil loss occurs on cultivated steep slopes.
LC change were associated with the increase in the area of cropland together with a decline in vegetation cover. The results indicated that the expansion of cultivation land and the decline of forest cover from 1986 to 2003 were the key contributors to the increment in average annual surface runoff (Tables 2 and 9). The amount of surface runoff changed from 319.91 mm in 1986 to 333.55 mm in 2003 as a consequence of the expansion of cultivated land and settlements and the decline in vegetation. As the amount of cultivated land increases, soil infiltration capacity decreases due to the compaction of soil, which decreases lateral flow and increases surface runoff. In contrast, the surface runoff declined from 333.55 mm in 2003 to 321.61 mm in 2020 due to the gradual increment of shrubland as a result of the implementation of current soil conservation practices in the watershed. A decrease in lateral flow occurred during the 1986-2020 LC change, which caused a reduction in water yield despite an increase in surface runoff. This could be attributed to the poor infiltration of non-vegetated areas.  (Table 2). This outcome is in line with a study conducted by Dias et al. [53], which found that watersheds with more vegetation tend to have higher ET. A high rate of ET was also reported by Leta et al. [8], which referred to it as the watershed's largest water consumer. A continuous decline in lateral flow was detected in the watershed. Lateral flow decreased by 3.04%, 2.20%, and 5.31% for the LC change scenarios of 1986-2003, 2003-2020, and 1986-2020, respectively (Table 9).
In comparison to the LC baseline year (1986), the watershed's average annual water yield increased in 2003 by 0.014% and decreased in 2020 by 0.05%. Similarly, the water yield gradually declined by 0.064% when the LC changed from 2003 to 2020. The annual average groundwater recharge declined by 1.50% between 1986 and 2003. However, in the second period, groundwater recharge increased by 2.01% as a result of the increase in shrubland in 2020 as compared to the 2003 LC (Table 2). In the first period, the annual streamflow of the watershed declined as urban area and cultivation land increased and forest cover, shrubland, and grassland declined. For all the change analysis periods, the watershed's mean annual streamflow slightly decreased over time.
In addition, the hydrological components were simulated on a sub-basin scale to better comprehend the effect of LC changes. The spatial distribution of hydrological components by LC change from 1986 to 2020 is depicted in Figure 8a,b. The hydrological processes have different characteristics in each sub-watershed. The spatial distribution of increment in ET corresponds to areas that were identified as being covered by forest, grassland, water bodies, and shrubland. Similarly, sub-basins represented by cultivated land generate higher surface runoff. This is consistent with other research studies that indicated higher surface runoff and lower ET in cultivation areas [53,54]. As depicted in Figure 9, the impact of LC change on groundwater flow, surface runoff, lateral flow, water yield, and ET was assessed by comparing monthly mean values. Surface runoff increased from June to September from 1986 to 2020, but water yield increased from June to August and decreased for all other months. Similar conclusions were reached by other researchers [8,45,55]. Flooding may occur as a result of increased surface runoff during rainy seasons, and water scheme operations may be impacted by a decline during the dry season. However, from the graph, the LC change on average monthly water yield and ET was negligible. Similarly, the average monthly streamflow during the wet season decreased from 522.70 m 3  As depicted in Figure 9, the impact of LC change on groundwater flow, surface runoff, lateral flow, water yield, and ET was assessed by comparing monthly mean values. Surface runoff increased from June to September from 1986 to 2020, but water yield increased from June to August and decreased for all other months. Similar conclusions were reached by other researchers [8,45,54]. Flooding may occur as a result of increased surface runoff during rainy seasons, and water scheme operations may be impacted by a decline during the dry season. However, from the graph, the LC change on average monthly water yield and ET was negligible. Similarly, the average monthly streamflow during the wet season decreased from 522.70 m 3   The decline in vegetation cover and the increment in cultivation land and settlement decreased the dry season flow throughout the periods. The findings of this study concur with those of other investigations [26]. The alteration of forests to cultivation land from 1985 to 2011 increased the mean wet monthly flow while it decreased the mean dry monthly flow, according to a research finding in the Angereb watershed [23]. In the Quaternary catchment, South Africa, streamflow increased from 2004 to 2013 due to an increment in cultivation land and a decline in the grass and woodland areas (9.8%) [55]. In general, it is observed that groundwater and surface runoff are more sensitive to changes in LC than lateral flow. The decline in vegetation cover and the increment in cultivation land and settlement decreased the dry season flow throughout the periods. The findings of this study concur with those of other investigations [26]. The alteration of forests to cultivation land from 1985 to 2011 increased the mean wet monthly flow while it decreased the mean dry monthly flow, according to a research finding in the Angereb watershed [23]. In the Quaternary catchment, South Africa, streamflow increased from 2004 to 2013 due to an increment in cultivation land and a decline in the grass and woodland areas (9.8%) [56]. In general, it is observed that groundwater and surface runoff are more sensitive to changes in LC than lateral flow.

Conclusions
This study examined the effects of LC change on hydrological components in the Muger watershed, UBNRB, Ethiopia, from 1986 to 2020. The study used the SWAT process-based hydrological model to identify the impacts of LC change on the hydrological processes of the watershed.
The effect of LC change on hydrological parameters was more highly pronounced at the sub-watershed scale than at the catchment level. The assessment of the impacts of LC change on the hydrologic processes was undertaken at the watershed and sub-watershed levels in order to identify the most severely impacted sections of the basin. The impacts of LC change on hydrological processes could be better studied if more uniform distribution of hydrological gauging stations were available.
The findings of the study demonstrated that from 1986 to 2020, ET increased, lateral flow decreased, and groundwater flow showed fluctuations. The LC change brought about more effects on surface runoff, groundwater, lateral flow, and ET, as indicated by the monthly and annual mean values. However, the impact of LC change is negligible on water yield, particularly when considering monthly average values. Due to its adverse impacts on surface runoff and water production, ET is a key factor in determining the availability of water. A slight change in ET brings about major changes in water availability. Therefore, future efforts should attempt to predict ET from different data types, including remotely sensed data.
The loss of forest directly affects the infiltration rate through an increment in surface runoff during the wet season and a decline during the dry season. The reduction in surface water could also have influenced the availability of water resources in the watershed and aggravated the downstream water shortage, especially between 2003 and 2020. As a result, managing and planning regional and local scale interventions in conservation strategies for water resources and land should be a priority. The findings of this study will be useful for policy-making and planning for the sustainable management of water resources and land in the Muger watershed. The findings from this study provide evidence for improving land and water resources management to ensure sustainable production and flow of water in the basin.
The methods used in this research have identified the impacts of LC change on hydrological parameters and are useful in predicting the hydrological repercussions of LC change in other basins. This approach has produced quantitative data that will help manage land and water resources effectively.

Data Availability Statement:
The data used in this study can be made available by the authors on reasonable request.