Hydrological Response of the Wami–Ruvu Basin to Land-Use and Land-Cover Changes and Its Impacts for the Future

The evaluation of the hydrological responses of river basins to land-use and land-cover (LULC) changes is crucial for sustaining water resources. We assessed the impact of LULC changes (1990–2018) on three hydrological components (water yield (WYLD), evapotranspiration (ET), and sediment yield (SYLD)) of the Wami–Ruvu Basin (WRB) in Tanzania, using the Soil and Water Assessment Tool (SWAT). The 1990 LULC imagery was used for SWAT simulation, and imagery from 2000, 2010, and 2018 was used for comparison with modelled hydrological parameters. The model was calibrated (1993–2008) and validated (2009–2018) in the SWAT-CUP after allowing three years (1990–1992) for the warm-up period. The results showed a decrease in WYLD (3.11 mm) and an increase in ET (29.71 mm) and SYLD (from 0.12 t/h to 1.5 t/h). The impact of LULC changes on WYLD, ET, and SYLD showed that the increase in agriculture and built-up areas and bushland, and the contraction of forest led to the hydrological instability of the WRB. These results were further assessed with climatic factors, which revealed a decrease in precipitation and an increase in temperature by 1 °C. This situation seems to look more adverse in the future, based on the LULC of the year 2036 as predicted by the CA–Markov model. Our study calls for urgent intervention by re-planning LULC and re-assessing hydrological changes timely.


Introduction
In recent decades, changes in the hydrological processes of a wide range of river basins have been observed with increasing frequency due to global trends in human population, climate change, and underlying surface features at global, regional, and local scales [1][2][3][4]. In addition to population, climate, and underlying surface features, landuse and land-cover changes (LULCC) have a potentially large impact on hydrological processes [5,6]. These include changes in water and sediment yields [7], surface runoff, baseflow, streamflow [8][9][10][11], soil water content [12], groundwater recharge [13], and annual river discharge [14]. In most cases, the influence of LULCCs is highly visible due to the modification of dense vegetation into agriculture that, by 2017, occupied 37% of the global land surface and employed about 50% of the world population [15,16]. Despite this, in developing countries (e.g., Tanzania) where arable lands are found in river basins, such statistics are lower due to lack of land-use planning; agriculture is critical for the livelihood of more than 70% of the general population in such countries [17,18]. Vegetation can have a

Study Area
The WRB is located between 5 • -7 • S and 36 • -39 • E in Tanzania, East Africa, (Figure 1) and encompasses a total area of approximately 66,820 km 2 , with an altitudinal gradient of approximately 2267 m above sea level. The WRB consists of two sub-basins, the Wami sub-basin (43,946 km 2 ) and the Ruvu sub-basin (22, Figure A3) [45,46].
The basin climate is characterised by semi-arid conditions in the western areas and a humid climate with high rainfall along the eastern areas towards the Indian Ocean. The average rainfall is approximately 500-760 mm per year in the west of the semi-arid highlands, 900-1000 mm in the central and coastal regions of the Wami sub-basin, and 800-2000 mm in the Ruvu sub-basin. The climate of the Ruvu sub-basin varies according to the topography. Most of its areas receive rainfall ranging between 800 and 1200 mm, except for the eastern slopes of the Uluguru Mountains, which have a mean annual rainfall in excess of 2000 mm, and the Ukaguru Mountains, which receive between 1000-1800 mm annually, as revealed by Kashaigili [47]. The WRB annual mean temperature ranges from 12 • C to 32 • C [48]. The basin population was estimated to be 10 million in 2018 [49]. Furthermore, according to the project conducted by the Japan International Cooperation Agency (JICA) [50], the total population in the WRB is forecasted to increase to 12.58 million in 2035.

SWAT Model Description
The SWAT hydrological model is a continuous time-and spatially-distributed basinscale model in which hydrological processes on the surface and subsurface are coupled in order to simulate the impacts of changes in management practices [35]. It was introduced in the 1990s to assist water resource managers in assessing hydrological responses for proper management [51,52]. As a built-in program, SWAT stands for the Soil and Water Assessment Tool in computing environments. It takes different forms of naming/extension depending on the key GIS software installed with it, such as ArcSWAT for the ArcGIS software, MSWAT for Map window software, and QSWAT for QGIS software. Despite the difference in its extension, this model used the same basic algorithms guided by the water balance equation (i.e., Equation (1)). Further details can be found in [53].
Notes: SW t is the final soil water content (mm), SW o is the initial soil water content on day i (mm), t is the time in days, R i is the amount of precipitation on day i (mm), Q surf,i is the amount of surface runoff on day i (mm), ET i is the amount of evapotranspiration on day i (mm), P i is the amount of water entering the vandose zone from the soil profile on day i (mm), and Q gw,i is the amount of return flow on day i (mm). In this study, we used ArcSWAT extension in the ArcGIS software.

Data Requirements for the Model Input
Four types of data were required to build up the SWAT model: topography; soils; LULC data; and basic climatic data. These data were acquired from different sources and were fit for the appropriate model inputs; each was pre-processed as explained below.

Topography and Soil Data
The topographical data were extracted from the digital elevation model, mapped by the Advanced Land Observing Satellite (ALOS) using the Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM) with a horizontal resolution of 30 m [54]. The choice of the ALOS DEM was made because of its highly accuracy results that have been tested in many countries such as Russia [55], Argentina [56], the Philippines [57], and Cameroon [58], which is close to the case study area. The topographic data were used to generate the stream networks of the WRB, including river length, slope, and width.
The soil data were retrieved from the Agricultural Organization (FAO) and used to determine soil parameters, such as texture, hydrologic soil group, and the available water content of the soil, as needed to run the SWAT model [59]. These data were available in the GIS shapefile form covering all over the world, and its soil codes matched the MSWAT extension. Thus, we added a new column to the ArcSWAT database for the soil data with names similar to those from the downloaded shapefile and used it when defining the hydrological response units (HRUs). In general, the soil data used in this study were chosen due to their high accuracy and wide use in many studies, such as in the Karnali River Basin of Nepal's Himalaya [60], Giba catchment, Tigray, Ethiopia [61], and in the Wami sub-basin-A part of the WRB in Tanzania by [62]. Figure 2 shows the topographical and soil data used for the model input. Table 1 provides a detailed description of the soil data types in the WRB.    The data for LULCs were obtained from Landsat satellite imagery, shared by the United States Geological Survey [63]. The data were downloaded at https://glovis.usgs.gov (accessed on 18 July 2021) from six Landsat scenes along WRS-2 paths 166 to 168 and rows 64 and 65 for 1990, 2000, 2010, and 2018. These datasets were retrieved from three different sensors of TM, ETM, and OLI during the dry season and with the threshold of less than 10% for cloud cover. The choice of this source of imagery was based on their availability that covers the full length of the studied time, the same spatial resolution of 30 m, and common Red-Green-Blue bands, enhancing the visual interpretation during classification. The 1990 imagery was used as a direct input for the simulation of the SWAT model for 29 years and hence, to provide a better understanding, we used the rest of the imagery for the comparative analysis of the modelled hydrological parameters, as detailed in the subsequent sections. The details of the dates, Landsat, and sensors of the downloaded data are presented in Table 2. All Landsat scenes were atmospherically corrected using FLAASH in ENVI 5.3 software through radiometric calibration and dark subtraction processes, and their digital numbers were converted into the top of atmosphere reflectance before being mosaicked into a single composite image cropped to the study area for LULC classification. To determine the LULC types, we identified training datasets and used a classifier selection tool to encode the sampling data for the classification process. The training datasets were randomly divided into two subsets: 60% was used for the training of the classifier and 40% was used for the accuracy assessment. Finally, the images were classified using artificial neural networks (ANNs). An ANN is a supervised classification method composed of several small processing units imitated with artificial neurons, and it is widely applied in recognition processes, such as LULC classification [64]. Eight LULC types were obtained from this process, as shown in Table 3. The classified imagery was further assessed for accuracy using Kappa coefficient analysis by comparing the pixel information to Google Earth and sentinel−2 at 10 m spatial resolution. The overall accuracy assessment was 80.8% for the LULC images in 1990, 78.9% in 2000, 88.4% in 2010, and 86% in 2018. These results are similar to those obtained in other studies, such as those by Abineh and Bogale [65] and Rwanga and Ndambuki [66]. For the use of the LULC imagery in the model, the LULC types were required to match the SWAT code database. Thus, we encoded them accordingly, as shown in Figure 3.

Hydroclimatic Data
The observed climatic data used in this study include the daily precipitation, temperature, wind speed, relative humidity, and solar radiation, which were obtained from the Tanzania Meteorological Agency (TMA), the University of Dar es Salaam-Water Resources Engineering Department database, and the Ministry of Water, Dodoma. All these data were sorted and processed in Excel for use in the SWAT model and, in order to obtain prior knowledge of the basin hydrology, the basic data of precipitation and temperature were assessed to determine their trends using OriginLab. For calibration and validation purposes, the flow data at the outlets of the sub-basin (1H8 and IG2) with 37% of missing data [70] were used in this study. These data cover a period of 29 years (1990-2018) and were collected from the WRB office in the Morogoro region. The choice of the use of the IG2 and 1H8 outlets was due to the quality of their data as indicated in many previous studies, such as [41,70], and to topography, as they are both found downstream. Thus, they are the main collectors of the entire streamflow of the basin. Detailed information on the outlets is presented in Table 4. Figure 4 is a flowchart used for processing all data used in this study.

SWAT Model Input
The SWAT modelling process began with watershed delineation based on digital elevation models (DEMs) using ArcSWAT in ArcGIS software. The watershed delineation involved five stages: DEM setup; stream definition; outlet and inlet definition; watershed outlet (s) selection; and the definition and calculation of the sub-basin parameters. Based on the topographical nature of the basin, both delineated stream and sub-stream networks were composed within the Wami and Ruvu sub-basins. Their outlets were also the same, coded as FLOW_OUT_13 and FLOW_OUT_24 for IG2 and 1H8, respectively ( Figure 5). Subsequently, the slope information obtained from the delineation process and the other soil and land-use data allowed us to define and compute the HRUs of the WRB. This step yielded 40 HRUs, which were, in principle, the result of the overlays of the identified data characteristics [71]. It is worth noting that in each HRU, water balance was the driving force behind all processes, such as surface and groundwater flow, the interception of precipitation, and the distribution of water in different soil profiles, and could impact the movement of sediments. Moreover, as defined by the HRUs, the SWAT model required critical threshold values for land-use, soil type, and slope, and the default values for most applications were 20% for land-use, 10% for soil, and 20% for slope [72]. The default values were applied sequentially during the definition of the HRUs to limit the number of HRUs in the WRB. Finally, before we ran the model, the SWAT model required the input of the climatic data that we imported to work with the HRU information in order to obtain the final results.

SWAT Model Input
The SWAT modelling process began with watershed delineation based on digital elevation models (DEMs) using ArcSWAT in ArcGIS software. The watershed delineation involved five stages: DEM setup; stream definition; outlet and inlet definition; watershed outlet (s) selection; and the definition and calculation of the sub-basin parameters. Based on the topographical nature of the basin, both delineated stream and sub-stream networks were composed within the Wami and Ruvu sub-basins. Their outlets were also the same, coded as FLOW_OUT_13 and FLOW_OUT_24 for IG2 and 1H8, respectively ( Figure 5). Subsequently, the slope information obtained from the delineation process and the other soil and land-use data allowed us to define and compute the HRUs of the WRB. This step yielded 40 HRUs, which were, in principle, the result of the overlays of the identified data characteristics [71]. It is worth noting that in each HRU, water balance was the driving force behind all processes, such as surface and groundwater flow, the interception of precipitation, and the distribution of water in different soil profiles, and could impact the movement of sediments. Moreover, as defined by the HRUs, the SWAT model required critical threshold values for land-use, soil type, and slope, and the default values for most applications were 20% for land-use, 10% for soil, and 20% for slope [72]. The default values were applied sequentially during the definition of the HRUs to limit the number of HRUs in the WRB. Finally, before we ran the model, the SWAT model required the input of the climatic data that we imported to work with the HRU information in order to obtain the final results.

Model Performance and Evaluation
After the simulation, the SWAT model was calibrated and validated using monthly flow data for the two gauging stations IG2 and 1H8 for FLOW_OUT_13 and

Model Performance and Evaluation
After the simulation, the SWAT model was calibrated and validated using monthly flow data for the two gauging stations IG2 and 1H8 for FLOW_OUT_13 and FLOW_OUT_24, respectively. The flow data at the outlets of sub-basins IG2 and 1H8 were considered because this was the only hydrological data available from the WRB from 1990 to 2018. The model calibration (1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008) and validation (2009-2018) periods were carried out by comparing the simulated and measured flow data, with the period of (1990-1992) left out as a warm-up period. Manual and auto-calibration routines [73] were performed. Manual calibration was performed to reduce the acceptable parameter ranges at each station [4]. The optimised parameter ranges suggested by the SWAT user's manual and other studies, such as those conducted by [74,75], were incorporated into a recently developed semiauto-calibration approach, SUFI-2, in the SWAT-CUP [76]. The sequential uncertainty fitting II (SUFI-2) [77] method calibration and uncertainty program [73] was used, with the Nash-Sutcliffe efficiency (NSE) [78], as the objective function to identify the most sensitive parameters to be calibrated in the WRB. The details of the statistical evaluation methods that include the NSE, the ratio of the root, mean square error to the standard deviation of measured data (RSR), and the percent bias (PBIAS) are presented in Table 5, and further details can be obtained from Moriasi et al. [79] and Moriasi et al. [80]. Table 5. The SWAT statistical equations used for calibration and validation processes and their level of performance.

Statistical Equation Value
Rating Performance

Results
The results of this research are presented in three phases: one for the LULCCs that are real and simulated based on the CA-Markov model, and two for the identified trends of the basic climatic data and the hydrological components modelled in the SWAT model. The results of the first and second phases are presented in a straightforward manner, as described in the methodology. The SWAT model results for each hydrological component were generated annually for the entirety of the study period; thus, for a better understanding, we chose the years 1990, 2000, 2010, and 2018 to compare the detected land-use and land-cover changes.

Trends of LULCCs
Tables 6 and 7 summarise the temporal trends of the LULCC analysis from 1990 to 2018 based on eight targeted LULC types. The spatial representation of each LULC from 1990 to 2018 is shown in Figure 6. Initially, the temporal trend of LULC in 1990, as indicated by the percentage of the total area of the WRB, was dominated by forest (58.1), followed by bushland (16.7), grassland (13.6), agriculture (10.6), water and wetland (0.6), bare soil (0.4), and built-up areas (0.1). This trend varied significantly for all LULC types in the years 2000, 2010, and 2018, as highlighted by the land covered by grassland and bushland areas ( Table 6). Grassland areas increased tremendously in 2000 and declined tremendously in 2018. In contrast, bushland was the only natural area that increased from 16.7% in 1990 to 24.9% in 2018, despite a decrease in 2000 (Table 7).
Moreover, the trend of LULCC indicates a general decrease in natural areas, including forests (15.4%), grasslands (6.7%), and water and wetland cover (0.3%) ( Table 7). On the other hand, the land uses supporting economic activities and population development increased, as indicated by agriculture and built-up areas ( Figure 6). Taken altogether, agriculture is the LULC with the highest gain (11.6%) and forest is the LULC with the highest loss (15.4%) over the yearly intervals from 1990 to 2018 (Table 7). These LULC results agree with the results of a previous study conducted by Twisa and Buchroithner [81] in the Wami sub-basin.
On the other hand, the trends of LULC simulated by the CA-Markov model show that the main land-use classes of agriculture, built-up areas, and forest agree with real land-use changes observed previously (Table 8)   that the main land-use classes of agriculture, built-up areas, and forest agree with re land-use changes observed previously (Table 8). These results were confirmed based o the relative accuracy (Kno = 0.8425, Klocation = 0.8357, Kstandard = 0.7987) between th classified map of 2018 and the modelled maps of 2036, as shown in in Figure 7.

Precipitation and Temperature Trends
The results of the precipitation and temperature studies from the Wami and Ruvu sub-basins showed an increase in temperature and a decrease in precipitation (Figure 8). The Wami sub-basin temperature varied between 30.5 • C and 32 • C from the 1990s to 2018. However, the Ruvu sub-basin increased from 28.5 • C to 29.5 • C. The variations in the temperature in these sub-basins demonstrate that the WRB temperature increased by approximately 1 • C from the 1990s to 2018 (Figure 8). Similarly, there was a sharp decrease in precipitation for both sub-basins during the study period. These trends were punctuated with concurrent events, such as in 1997, when both sub-basins recorded the highest precipitation, and in 2003 and 2010, when the same sub-basins recorded the lowest precipitation over the entire study period.

Precipitation and Temperature Trends
The results of the precipitation and temperature studies from the Wami and Ruvu sub-basins showed an increase in temperature and a decrease in precipitation (Figure 8). The Wami sub-basin temperature varied between 30.5 °C and 32 °C from the 1990s to 2018. However, the Ruvu sub-basin increased from 28.5 °C to 29.5 °C. The variations in the temperature in these sub-basins demonstrate that the WRB temperature increased by approximately 1 °C from the 1990s to 2018 (Figure 8). Similarly, there was a sharp decrease in precipitation for both sub-basins during the study period. These trends were punctuated with concurrent events, such as in 1997, when both sub-basins recorded the highest precipitation, and in 2003 and 2010, when the same sub-basins recorded the lowest precipitation over the entire study period.

Spatio-Temporal Water Yield (WYLD) Distribution
Spatially, the results of the WYLD distribution show an annual variation between the four years specified by the study methods ( Figure 9). These results also comprise the temporal trends for each sub-basin, which indicates a reduction in WYLD from 169.38 mm to 166.27 mm, and from 173.59 mm to 170.54 mm for the Wami and Ruvu sub-basins, respectively ( Table 9). The SWAT-defined WYLD was the aggregate sum of water leaving the HRU and entering the principle channel during the time step, i.e., including surface runoff and ground water flow [82]. The WYLD results also comprise the major components of surface flow and groundwater recharge, which revealed the increase in surface flow and decrease in groundwater flows for both sub-basins (Table 9). Similar results were obtained in previous studies by Nobert and Jeremiah [41] and Twisa [75], who used a section of the Wami sub-basin to study the effects of LULCCs on hydrology. In general, the results for the spatial pattern of WYLD in Figure 9 show that most of the agricultural areas in the western areas of the WRB have low water yields. These results indicate that the significant increase in LULCCs associated with expanding agricultural areas that draw water for irrigation has an impact on the changes in WYLD in the WRB (see Figures 6 and 9 for comparison). Correspondingly, our study associated the WYLD reductions with the semi-arid climatic conditions in western areas of the WRB.

Spatio-Temporal Water Yield (WYLD) Distribution
Spatially, the results of the WYLD distribution show an annual variation between the four years specified by the study methods ( Figure 9). These results also comprise the temporal trends for each sub-basin, which indicates a reduction in WYLD from 169.38 mm to 166.27 mm, and from 173.59 mm to 170.54 mm for the Wami and Ruvu sub-basins, respectively ( Table 9). The SWAT-defined WYLD was the aggregate sum of water leaving the HRU and entering the principle channel during the time step, i.e., including surface runoff and ground water flow [82]. The WYLD results also comprise the major components of surface flow and groundwater recharge, which revealed the increase in surface flow and decrease in groundwater flows for both sub-basins (Table 9). Similar results were obtained in previous studies by Nobert and Jeremiah [41] and Twisa [75], who used a section of the Wami sub-basin to study the effects of LULCCs on hydrology. In general, the results for the spatial pattern of WYLD in Figure 9 show that most of the agricultural areas in the western areas of the WRB have low water yields. These results indicate that the significant increase in LULCCs associated with expanding agricultural areas that draw water for irrigation has an impact on the changes in WYLD in the WRB (see Figures 6 and 9 for comparison). Correspondingly, our study associated the WYLD reductions with the semi-arid climatic conditions in western areas of the WRB.

Simulated Evapotranspiration Trend (ET)
The SWAT-simulated results for ET show fluctuations from 1990, with a gradual increase in 2018 ( Figure 10). From 1990 to 2018, the trends of ET increased by 29.71 mm, from 95.01 mm to 124.72 mm. This increase indicates the impact of LULC on ET, mainly resulting from the physical changes of the WRB land surface, which affects the efficiency of ET as a process (also see Figures 6 and 10). Furthermore, Figure 11, which indicates the spatial distribution of ET in the basin, shows that ET in the far western areas is substantially lower than the eastern parts and towards the coastal areas of the WRB. This variation could be attributed to the bushland, which has expanded toward the coastal areas ( Figures  6 and 11). Similarly, the climatic conditions of the river basin, especially in the coastal and estuarine environments, are also likely to have contributed towards this variation. In general, the variation in ET, as shown in Figure 10, is largely consistent with previous studies, such as Wambura [83] who used a Moderate Resolution Imaging Spectroradiometer

Simulated Evapotranspiration Trend (ET)
The SWAT-simulated results for ET show fluctuations from 1990, with a gradual increase in 2018 ( Figure 10). From 1990 to 2018, the trends of ET increased by 29.71 mm, from 95.01 mm to 124.72 mm. This increase indicates the impact of LULC on ET, mainly resulting from the physical changes of the WRB land surface, which affects the efficiency of ET as a process (also see Figures 6 and 10). Furthermore, Figure 11, which indicates the spatial distribution of ET in the basin, shows that ET in the far western areas is substantially lower than the eastern parts and towards the coastal areas of the WRB. This variation could be attributed to the bushland, which has expanded toward the coastal areas (Figures 6 and 11). Similarly, the climatic conditions of the river basin, especially in the coastal and estuarine environments, are also likely to have contributed towards this variation. In general, the variation in ET, as shown in Figure 10, is largely consistent with previous studies, such as Wambura [83] who used a Moderate Resolution Imaging Spectro-radiometer (MODIS) imagery program to study ET and potential evapotranspiration (PET) in the Ruvu sub-basin. (MODIS) imagery program to study ET and potential evapotranspiration (PET) in the Ruvu sub-basin.

Spatio-Temporal Sediment Yield (SYLD) Distribution
Considerable changes in sediment yield in tons per hectare (SYLD t/h) occurred in the WRB from 1990 to 2018. The main change was a significant increase from 0.12 t/h to 1.5 t/h from 1990 to 2018 ( Figure 12). Figure 13 shows the spatial distribution in the SYLD simulated from 1990 to 2018. In the starting year of the simulation, 1990, SYLD was more exposed in HRU 34 of the Ruvu sub-basin. In 2000, most SYLDs were imitated in HRUs 30, 33, and 37 of the same Ruvu sub-basin. In contrast, for the year 2010, the simulated results indicated that both sub-basins of the WRB were affected by the SYLD, with the Wami sub-basin being affected most at the downstream area (HRU 13). The Ruvu sub-basin affected both upstream and downstream, mostly in HRUs 30 and 37. Subsequently, for the last simulated year of 2018, the Wami sub-basin was the only area in the WRB affected by the SYLD, particularly in the upstream areas (HRUs 10 and 2). In general, the spatial distribution of the SYLD from 1990 to 2018 is characterised by notable features, such as in 1990, 2000, and 2010, where the Ruvu sub-basin was the only area simulated with the significant SYLD; in 2018, the situation changed to the Wami sub-basin, a situation that corresponds to the expansion of agricultural land cover (see Figures 6 and 13 for comparison). Considerable changes in sediment yield in tons per hectare (SYLD t/h) occurred in the WRB from 1990 to 2018. The main change was a significant increase from 0.12 t/h to 1.5 t/h from 1990 to 2018 ( Figure 12). Figure 13 shows the spatial distribution in the SYLD simulated from 1990 to 2018. In the starting year of the simulation, 1990, SYLD was more exposed in HRU 34 of the Ruvu sub-basin. In 2000, most SYLDs were imitated in HRUs 30, 33, and 37 of the same Ruvu sub-basin. In contrast, for the year 2010, the simulated results indicated that both sub-basins of the WRB were affected by the SYLD, with the Wami sub-basin being affected most at the downstream area (HRU 13). The Ruvu subbasin affected both upstream and downstream, mostly in HRUs 30 and 37. Subsequently, for the last simulated year of 2018, the Wami sub-basin was the only area in the WRB affected by the SYLD, particularly in the upstream areas (HRUs 10 and 2). In general, the spatial distribution of the SYLD from 1990 to 2018 is characterised by notable features, such as in 1990, 2000, and 2010, where the Ruvu sub-basin was the only area simulated with the significant SYLD; in 2018, the situation changed to the Wami sub-basin, a situation that corresponds to the expansion of agricultural land cover (see Figures 6 and 13 for comparison).

Calibration and Validation of the SWAT Model
According to [84], a hydrological model should be calibrated to estimate the model parameter values that provide the best possible hydrological results of interest. Hence, in this study, the SWAT model was calibrated using two gauging stations IG2 and 1H8 for FLOW_OUT_13 and FLOW_OUT_24. The calibration (1993-2008) and validation (2009-2018) were performed using the monthly flow data and the process started with a sensitivity analysis that identified seven key sensitive parameters governing the calibrated and validated outputs (Table 10). These parameters were manually adjusted iteratively to suit the satisfactory range before the auto-calibration was ran. The parameter ranges (minimum, maximum, and the fitted value) used for the model sensitivity analysis were found within the suggested ranges of the SWAT user's manual and based on other studies, such as those conducted by [74,75]. Of the seven parameters, the four most sensitive parameters in our study were R_CN2.mgt, SURLAG.bsn, SOL.AWC.sol, and V_ALPHA-BF.gw.

Calibration and Validation of the SWAT Model
According to [84], a hydrological model should be calibrated to estimate the model parameter values that provide the best possible hydrological results of interest. Hence, in this study, the SWAT model was calibrated using two gauging stations IG2 and 1H8 for FLOW_OUT_13 and FLOW_OUT_24. The calibration (1993-2008) and validation (2009-2018) were performed using the monthly flow data and the process started with a sensitivity analysis that identified seven key sensitive parameters governing the calibrated and validated outputs (Table 10). These parameters were manually adjusted iteratively to suit the satisfactory range before the auto-calibration was ran. The parameter ranges (minimum, maximum, and the fitted value) used for the model sensitivity analysis were found within the suggested ranges of the SWAT user's manual and based on other studies, such as those conducted by [74,75]. Of the seven parameters, the four most sensitive parameters in our study were R_CN2.mgt, SURLAG.bsn, SOL.AWC.sol, and V_ALPHA-BF.gw.   Figure 14 shows the comparative results of the observed and simulated flows for the outlets FLOW_OUT_13 and /FLOW_OUT_24 of IG2 and 1H8 during the calibration (1993-2008) and validation periods (2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017)(2018). The graphical comparison between the observed and simulated annual river flow exhibits good agreement between the simulated and observed flows for both the calibration and validation periods ( Figure 14). Table 11 presents the statistical values that govern the SWAT-WRB model, the NSE, RSR, and PBIAS, which were found to be satisfactory for both the calibrated and simulated periods; these results also comply with others in the literature, such as in [79,80]. Therefore, our results indicate that the SWAT model performed well in simulating the hydrological conditions of the WRB.

Impacts of LULCC on River Basin Hydrology over 29 Years (1990 to 2018)
LULC changes in the WRB have already been reported. The results for the studied years (1990,2000,2010, and 2018) on the different types of LULC indicate that most of the WRB natural areas, including forests and grasslands, were intensively converted into agricultural land from 1990 to 2018; further, these changes are anticipated to continue towards 2036 (Table 8). For example, agricultural land doubled from 10.6% in 1990 to 22.2% in 2018 and is predicted to continue to increase to 31.0% by 2036 (Tables 6 and 8). Similarly, built-up areas increased between 1990 and 2018 (0.1% to 0.9%) and are expected to increase to 1.4% by 2036. Moreover, bushland increased from 16.7% in 1990 to 24.9% in 2018, despite a decrease in 2000, but is expected to increase to 27.1% by 2036. This land-use pattern reflects the influence of the population increase from 3.5 to 10 million people in the basin area from 1988 to 2018 [49]. This LULC pattern also highlights the common influence of economic forces that reinforce anthropogenic stimulus on LULC change in the river basin [85,86], and this is the main reason that the hydrological conditions of the wide-ranging river basins, including the WRB, have changed due to the transformation of natural vegetation, as substantiated by the SWAT hydrological model.
The observed LULC pattern is well established in the scientific literature-the deforestation of native vegetation influences several hydrological processes, resulting in changes in water and sediment yields [87], surface runoff [88], base and streamflow [33,89], river discharge [90], and latent heat fluxes to the hydrological cycle, including transpiration, interception loss, and evapotranspiration [91]. In the WRB, an increased conversion of the largest areas of natural vegetation into agriculture cover has played a great role in the increase in sediment yields generation from 0.12 t/h to 1.5 t/h, and the changes in WRB hydrological processes including WYLD, which decreased by 3.11 mm from 1990 to 2018. These results are in agreement with a previous study in the basin, particularly in the Wami sub-basin [41]. Spatially, the changes in WYLD and SYLD were consistent with the land-use pattern associated with agricultural expansion, which seems to intensively affect the western areas of the WRB, especially in recent years ( Figure 15).
Other studies, including [92][93][94][95], have also noted similar LULC impacts on sediment loads. When sediments are eroded from their original places, some are deposited on the way while others continue to river basin streams [96]. Both transported and deposited sediments alter the basin's hydrological components, such as the total WYLD, which was found to decrease in the WRB, particularly in the western part ( Figure 14). Recalling statistically, the change in WYLD was observed from its components that indicated a decrease in groundwater flow with an increase in surface flow by 3.23 mm and 4.11 mm for Wami and Ruvu sub-basins, respectively ( Table 9). The results suggest that the decrease of about 1,028,091 ha in forest cover within the basin area from 1990 to 2018 might also have caused a decrease in WYLD due to the decline in interception, which leads to more surface runoff in most cases. In addition, we associated the decrease in WYLD with the decrease in rainfall that was studied over the sub-basins of the WRB from 1990 to 2018 (Figure 8).
On the other hand, from the SWAT-simulated hydrograph for the observed and simulated water flow for both calibrated (1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008) and validated periods (2009-2018), the LULCC impacts that occurred from the 1990s to 2018 led to a decrease in annual flow, especially in 2018 ( Figure 14). This is despite the good agreement between the simulated and observed discharges for both the calibration and validation periods, as revealed by the model performance statistics of FLOW_OUT_13 for the Wami sub  Table 11. Apart from the clear impacts of LULCCs associated with agriculture in the two sub-basins, the difference in the annual average flow observed in our study might have been caused by the difference in the climatic conditions of the two sub-basins: the semi-arid conditions of the far western areas of the FLOW_OUT_13 and the humidity and high rainfall of most areas covered by the FLOW_OUT_24. In addition, the high sensitivity of the R_CN2 parameter, which ranked first in the sensitivity analysis and governed the curve number of surface flow in the river basin (Table 10), could also contribute to the observed flow difference for the two sub-basin outlets. The high sensitivity of CN2 to streamflow has been reported in other studies, such as those conducted by Ndomba et al. [97], Mulungu and Munishi [98], Birhanu and Zemadim [99] in Tanzania river basins, and Baker and Miller [52] in East African watersheds.
Parallel to the sensitivity of the streamflow parameters and the LULCCs and their impacts on river flow, the rise and substantial variation in the spatial pattern of ET found in our study from 1990 to 2018 (Figures 10 and 11) could be attributed to the same source. As noted in Figures 6 and 11, the variation in ET seems to be mostly attributed to the increase in bushland. These results strongly agree with the SWAT model, which computes the ET based on the plant transpiration rate, maximum soil evaporation rate, and canopy water evaporation interception [74]. Moreover, as our study revealed an increase in temperature in the basin by 1 • C from 1990 to 2018, this could also be the reason for the observed increase in ET. Based on the comparison with the spatial LULC distribution (Figures 6 and 11), the ET decreased more from the far west to the central and lower land areas. A study by FIU-GLOWS [46] reported similar findings of an increase in evaporative water loss in the lower streams of the WRB. It is worth noting that, as the lower land areas were also highly populated as indicated by the rapid increase in built-up areas (Figure 6), the increase in ET to the lower areas of the WRB might further lead to hydrological complications in the WRB rivers, tributaries, and reservoirs and affect the amount of water resources needed by the populations residing in those areas.
Therefore, based on the results of this study and those drawn from other research, including Wambura [42] and Wambura [35], it is clear that the LULC changes contributed to modifications in the basin land and hydrological parameters. The trends of the increase in all tested parameters signify the instability of the WRB hydrology. It is worth noting that the future trends of the modelled parameters depend on the good practice of the LULC in the river basin. For this reason, the researchers investigated and identified the land-use plan (LUP) currently being implemented in the basin (Appendix A, Figure A4) and predicted the LULC change of the river basin up to 2036 in order to determine the future of the LULC and its anticipated impacts on the identified hydrological parameters. In the first case, the existing LUP has six main land uses. In comparison with the current LULC of 2018, land-use with high hydrological effects (i.e., agriculture and built-up areas) is allocated more space (>50%). This implies that the current plan favours the hydrological instability of the WRB. The predicted LULC change up to 2036 is consistent with current patterns (Figure 7). The anticipated increase in agricultural land, bushland, and built-up areas and reduction in natural forest and water and wetlands follows the current LULC trend in the region; if this trend continues without remedial measures, its impacts on the hydrological components will increase considerably in the future, as simulated in our study. Furthermore, this will become more complicated with climate change, rapid population growth, and the need to achieve economic growth within the basin [49].

Hydrological Stability of the WRB
The SWAT results demonstrate that the LULC types considered in this study are significant and, therefore, important in affecting the hydrological components of the WRB. As statistically shown, the decrease in WYLD and the increase in SYLD and ET clearly define the instability of the WRB. On this basis, the detected changes in natural forests, built-up areas, and agricultural areas should be regulated in order to decrease the related impacts on the hydrological components and to obtain the hydrological stability of the basin. In this course, the land-use plan, which is in practice in the basin, is among the key areas for intervention since it has awarded a high priority to agricultural development. Generally, agriculture is one of the major income sources [45,95,100] and is supported with several strategies, such as Kilimo Kwanza (Agricultural First) and the Southern Agriculture corridor of Tanzania [101,102]. Such initiatives are implemented more in the arable lands of river basins, including the WRB [102]. Despite its significance, if not controlled by rationing as in the current LUP, its adverse impacts on basin hydrology are expected, as evidenced by the predicted LULC changes up to the year 2036.
Parallel to these observations, the specific hydrological response results showed that the far western areas of the WRB significantly reduced WYLD and increased SYLD. In contrast, the north-eastern and low-land areas also showed the same hydrological complexities due to indicators such as the rapid increase in built-up land and the significant increase in ET. Thus, to reverse this situation, our study proposes re-planning the river basin land-use and re-assessing the hydrologic changes to sustain the hydrological stability of the WRB. For the first case, we think that performing suitability analyses before allocating land-use could provide a better way of maintaining the hydrological conditions of the WRB. The latter should consider assessing the hydrological responses of the WRB using more parameters and with different models, including SWAT, in order to better understand the ecological conditions of the basin. Taken together, proper LUP and the regular monitoring of LULCCs and climatic factors are essential for the stability of the WRB.

Conclusions
This study assessed the impact of LULC changes on the hydrology of the WRB over 29 years (1990 to 2018) using SWAT in ArcGIS software and predicted its future trends using the CA-Markov model up to the year 2036. Landsat data were used to acquire land-use and land-cover changes and, together with hydro-meteorological, topological, and global soil data, to model the hydrology of the WRB. The substantial LULCCs are mainly due to the conversion of natural forests into agricultural areas. The rapid decrease in the forest areas, from 58.1% in 1990 to 42.7% in 2018, is primarily caused by the inefficient land-use plans drawing substantial land for agricultural use, which grew from 10.6% in 1990 to 22.2% in 2018 and is expected to increase to 31.0% by 2036. Parallel to these findings, the hydrological response to the LULC changes included a decrease in WYLD by 3.11 mm, an increase in ET by 29.71 mm, and an increase in SYLD from 0.12 t/h to 1.5 t/h, which define the hydrological instability of the WRB. Following this observation and examination of the causes, we advise possible mitigation measures that may be used to improve the current conditions of the river basin. Generally, the findings of this study provide baseline information for hydrologists, land-use planners, policymakers, resource managers, and decision-makers to elicit better planning and management decisions in the Wami-Ruvu Basin. Future studies should focus on investigating the hydrological conditions of the WRB by considering other key drivers of change besides LULCCs. They should be further extended to evaluate other hydrological parameters in order to better understand the hydrology of the WRB for the optimal use of water resources.