Evaluating and Predicting the Effects of Land Use Changes on Hydrology in Wami River Basin, Tanzania

Understanding the variation in the hydrological response of a basin associated with land use changes is essential for developing management strategies for water resources. The impact of hydrological changes caused by expected land use changes may be severe for the Wami river system, given its role as a crucial area for water, providing food and livelihoods. The objective of this study is to examine the influence of land use changes on various elements of the hydrological processes of the basin. Hybrid classification, which includes unsupervised and supervised classification techniques, is used to process the images (2000 and 2016), while CA–Markov chain analysis is used to forecast and simulate the 2032 land use state. In the current study, a combined approach—including a Soil and Water Assessment Tool (SWAT) model and Partial Least Squares Regression (PLSR)—is used to explore the influences of individual land use classes on fluctuations in the hydrological components. From the study, it is evident that land use has changed across the basin since 2000 (which is expected to continue in 2032), as well as that the hydrological effects caused by land use changes were observed. It has been found that the major land use changes that affected hydrology components in the basin were expansion of cultivation land, built-up area and grassland, and decline in natural forests and woodland during the study period. These findings provide baseline information for decision-makers and stakeholders concerning land and water resources for better planning and management decisions in the basin resources’ use.


Introduction
Due to a global increase in population, a substantial rise in fluctuations in many types of land use (e.g., urbanization, deforestation, and agriculture) has been documented, resulting in challenges to water resource availability [1]. Globally, basin resources have been facing severe pressure to provide for the needs of an increasing population, causing intensified agricultural practices and deforestation [2]. As a result, the future accessibility of adequate supplies of water for agricultural and human needs has become critical in many areas [3,4]. Several areas face water scarcity due to variations in land use types, pollution, agriculture, climate change, and human activities [5]. Land use changes cause adverse effects on catchments by affecting infiltration, flood peaks, evapotranspiration, groundwater recharge, sediment transportation, surface runoff, and water quality, among others [6]. In several parts of the earth, in the last few years there has been a severe change in the land use and consequently change in the hydrological and ecological systems [7]. At the watershed scale, land cover changes pose risks of water resources due to loss of vegetative cover, which control the movement of run-off into river systems [8,9].
The influence of land use fluctuations on catchment hydrological components has been growing, where semi-arid and arid areas are particularly susceptible [10]. The Wami River Basin is a significant area, due to its different benefits to a diverse range of shareholders [11]. The effects of hydrological changes caused by expected changes in land use may be severe for the Wami river system, given its primary role for water, providing food and livelihoods [12]. This alteration towards increasing land use has caused several fluctuations in services and roles, subsequently causing an inclusive degradation in provisioning services from the natural resources of the Wami River Basin [13]. Thus, changes in freshwater resource availability is likely to be one of the most critical consequences affecting the sustainable growth of livelihoods (and life) in the basin [14,15]. Relative change in streamflow is the most common estimate for studying the effect of land-use change on water resources and is considered to be useful in decision-making processes for water resource management [16][17][18][19][20][21][22][23].
Many studies have reported the impact of land use changes on hydrology at different spatial and temporal scales [24][25][26][27][28][29][30]. Techniques for evaluating the hydrological impacts of land use changes in basins include hydrological modeling, multivariate statistics, and paired catchments (e.g., [31][32][33]). Physically based models, such as the Soil and Water Assessment Tool (SWAT), are well recognized models for analyzing the impact of land management practices on water, sediment, and agricultural chemical yields in large and complex catchments [34]. SWAT is a semi-distributed physically based simulation model, and can predict the impacts of land use and management practices on hydrological regimes in watersheds with varying soil, land use and management conditions over long periods, primarily as a strategic planning tool [35]. The SWAT model application has gained reliability and consistency while being widely used to assess the effect of land use dynamics on hydrology [36][37][38][39][40][41][42].
Assessing changes in hydrological processes caused by land use change provides a demanding challenge for ecological researchers [43,44]. Studying the inconsistencies in the hydrology feedback of a catchment linked to land use change is crucial for developing an action plan and management of the water resources [45,46]. Therefore, in this study, we investigate the effects of land use change on various hydrological elements of the basin [47]. Initially, we calibrate and validate the SWAT model and then assess the effects of land use changes on the hydrological processes in the basin. These results are intended to assist in managing and monitoring land practices and active water management plans in basins similar to the Wami River basin, which are subject to agricultural expansion, population growth, and urbanization.

Study Area
The study location of the Wami River Basin (Figure 1) ranges between 5-7 • S and 36-39 • E covering the semi-arid parts in central Tanzania. It covers 41,167 km 2 and is characterized by humid inland swamps linked to the Indian Ocean. The basin topographically ranges from 2-2370 m above sea level. Kinyasungwe, Mkondoa, and Wami are three major sub-catchments forming the Wami River Basin [48]. Climatic conditions in the Wami River Basin are both spatially and seasonally variable with an average annual rainfall estimated to be 550-750 mm in the highlands near Dodoma, 900-1000 mm in the middle parts of the basin near Dakawa and at the river's estuary in Bagamoyo [49]. The lowlands are warm, whereas the highlands are cold. The Wami River Basin is comprised of five basin draining zones-(i) upland plains, (ii) mountain torrents, (iii) inland plains, (iv) rejuvenated cascades, and (v) coastal plains-which characterize the topography of the study area. The river's path varies: a relatively straight and confined pattern is observed in elevated and mountainous areas, but a meandering system dominates the coastal and inland plains, interfered with in a few areas by minor cascades at certain gradient levels [50]. On the other hand, a cascade zone can be seen near Wami at Mandera, upstream from the bridge crossing on the Chalinze-Segera highway [51].

Land Use Change Analysis and Prediction
Land use assessment was carried out based on images from Landsat-5, Landsat-7, and Landsat-8. The USGS Global Visualization Viewer website provided the images, from which those with 30 m resolution and less cloud cover were chosen using Path/Row 168/64, 168/65, 166/64, 167/64, and 167/65 ( Table 1). The image classification was hybrid, comprised of both unsupervised and supervised classification techniques [52,53]. The Iterative Self-Organizing Data Analysis (ISODATA) clustering algorithm carried out unsupervised classification, while the Maximum Likelihood Classification (MLC) algorithm was engaged for supervised classification [54,55]. The delineated types were woodland, bushland, wetland, cultivated land, grassland, built-up area, water, and natural forest ( Table 2). To improve classification accuracy, post-classification refinement was used [56]. The mixed pixels problem was later addressed by visual interpretation. Later, the stratified random method was used to assess the accuracy for both of two land use maps in the study area representing different land use classes. Moreover, the non-parametric Kappa test was used for classification accuracy, in order to account for the diagonal and confusion matrix elements [57].

Land Use Change Analysis and Prediction
Land use assessment was carried out based on images from Landsat-5, Landsat-7, and Landsat-8. The USGS Global Visualization Viewer (https://glovis.usgs.gov) provided the images, from which those with 30 m resolution and less cloud cover were chosen using Path/Row 168/64, 168/65, 166/64, 167/64, and 167/65 ( Table 1). The image classification was hybrid, comprised of both unsupervised and supervised classification techniques [52,53]. The Iterative Self-Organizing Data Analysis (ISODATA) clustering algorithm carried out unsupervised classification, while the Maximum Likelihood Classification (MLC) algorithm was engaged for supervised classification [54,55]. The delineated types were woodland, bushland, wetland, cultivated land, grassland, built-up area, water, and natural forest ( Table 2). To improve classification accuracy, post-classification refinement was used [56]. The mixed pixels problem was later addressed by visual interpretation. Later, the stratified random method was used to assess the accuracy for both of two land use maps in the study area representing different land use classes. Moreover, the non-parametric Kappa test was used for classification accuracy, in order to account for the diagonal and confusion matrix elements [57].

Bushland
Mainly comprised of plants that are multi-stemmed from a single root base.

Woodland
An assemblage of trees with canopy ranging from 20-80% but which may, on rare occasions, be closed entirely. Wetland Low-lying, uncultivated ground where water collects; a bog or marsh. Cultivated land Crop fields and fallow lands.
Built-up area Residential, commercial, industrial, transportation, roads, and mixed urban. Grassland Mainly composed of grass.
Natural forest A continuous stand of trees, many of which may attain a height of 50 m; includes natural forest, mangroves, and plantation forests.

Water
River, open water, lakes, ponds, and reservoirs.
Simulation of future land use change was based on CA-Markov chain analysis. The CA-Markov chain is a statistical tool which can define the chance of the land use changing from one date to another by establishing a transitional probability matrix between the first period and second period, based on the neighborhood effects [58][59][60]. The CA-Markov model was developed using the IDRISI Selva 17.0 software following two significant steps: first, to calculate the conversion probability, a conversion area matrix followed by layers of conditional probability is constructed by employing Markov chain analysis; the second stage is spatial description of the land use analysis, simulated by a CA spatial operator and multi-criteria evaluation (MCE) to predict and simulate the future land use change [61,62]. The classified land use maps for 2000 (which represents the past) and 2016 (which represents the present) were used as input data into the model, in order to predict the expected change for the year 2032. We performed model validation by comparison of the 2016 simulated land use map with the actual satellite land use map based on Kappa statistics [63]. The VALIDATE tool was used to compute the Kappa statistics for the projected land use and Kappa index >70% indicates a valid model [64].

Soil and Water Assessment Tool (SWAT) Model
The SWAT model was applied in the Wami River Basin to assess the effects of land use changes on hydrological components. A significant number of SWAT models have been applied to study hydrology in watersheds in different areas of the world. More details about the SWAT model can be found in Arnold et al. [65] and at http://swat-model.tamu.edu/. The SWAT model input data included a Digital Elevation Model (DEM), climate data, soil data, and a land use map. The DEM was obtained from the Shuttle Radar Topography Mission (SRTM) data set with a resolution of 30 m × 30 m. The soil data was obtained from the FAO with 1 km resolution. Land use data for two periods (2000 and 2016) and the projected land use (2032) maps were used to assess the effects of land use changes on the hydrology of the Wami River Basin. The observed climatic data, including daily rainfall and minimum and maximum temperatures values from 1990-2012 for the Wami River Basin, were obtained from the Tanzania Meteorological Agency (TMA). Daily runoff data from the Mandera gauge station (at the outlet of the catchment) from 1990-2012 was collected from the Ministry of Water, Tanzania. The SWAT 2012 version was used, based on the slope, soil data, and land use. The study area was divided into 23 sub-basins, which were later subdivided into 71 HRUs. The SWAT model simulates the land stage of the hydrological cycle established on the water balance equation [65]. The Soil Conservation Service Curve Number (SCS CN) method was employed in the study to compute the runoff, using Equation 2 [66][67][68]: where SW t is the final water content (mm), SW 0 is the initial water content (mm), t is the time (days), R day is the amount of precipitation on a specific day (mm), Q sur f is the runoff amount on a specific day (mm), E a is the evapotranspiration amount on specific day (mm), W sweep is the amount of water percolation into the vadose zones on a specific day (mm), and W gw is the return amount of flow on a specific day (mm).
where Q sur f is the daily surface runoff (mm), R day is the precipitation depth for the day (mm), and S is the retention parameter (mm). The retention parameter (S) is given by the equation where S is the drainable volume of soil water per unit area of the saturated thickness (mm/day) and CN is a curve number. I a is usually given as 0.2S. Therefore, Equation 2 becomes Nineteen flow parameters were selected to detect the SWAT sensitive parameters, based on the literature [69][70][71][72][73][74]. Global sensitivity analysis [73,[75][76][77] was employed in SWAT-CUP 2012 version 5.1.4 by allowing each parameter to be individually varied at a time [78]. To analyze the significance sensitivity, the p-value and t-Stat were used [73,[75][76][77]: a high t-test value means high sensitivity, whereas a lower p-value indicates higher significance [73,75,77]. Model calibration is a process of adjusting model parameters, within the suggested ranges, based on observed data to approve the same response over time [78][79][80], and validation is the process of proving the representation of the parameters by simulating the observed data with an independent data set exclusive of changing model parameters [78][79][80]. The Sequential Uncertainty Fitting (SUFI-2) algorithm was used to carry out calibration and validation, as established in the SWAT-CUP user manual [75]. SUFI-2 is a semi-automated calibration and uncertainty analysis algorithm [81] that accounts for all sources of uncertainty, including uncertainty in the driving variables (e.g., rainfall), conceptual model, parameters and measured data [69,79,80]. SUFI-2 is the most popular calibration and uncertainty analysis program and, has been used in many studies [69,77,80,82]. The data set was split into a warm-up period (1990-1999), a calibration period (2000)(2001)(2002)(2003)(2004)(2005)(2006), and a validation period (2007)(2008)(2009)(2010)(2011)(2012). There are several indices available to check the performance of the SWAT model; the following were used to evaluate the performance of the SWAT model, as recommended by Moriasi et al. [83]: The detailed land use maps (2000 and 2016) and the predicted land use map (2032) were used individually, whereas the other SWAT model inputs were the same as those used to analyze the basin's response to land use change. The simulated results were used to assess the effects of historical and future land use changes on water hydrology and to measure the input of changes in separate land use types to changes in hydrological components at the catchment.

Pearson Correlation and Partial Least Squares Regression
The relationship between land use change and hydrological components was assessed using the pairwise Pearson correlation [39,41,74]. Pearson correlation analysis was adopted to analyze the linear correlations between land use change (dependent variables) and hydrological components (independent variables). Significant values were determined by a 95% confidence interval from the sample size, n equal to four (independent variables), and the Pearson correlation coefficient. Thus, land use change influences the hydrological components more firmly as the Pearson correlation coefficient value approaches 1. Moreover, the hydrological effects of separate land use changes were investigated using correlation analysis by the Partial Least Squares Regression (PLSR) model. PLSR is a multivariate analysis method (Equation 9) which combines and generalizes features from primary element analysis and multiple regressions [84]. PLSR forecasts a dependent variable(s) set from an independent variables set [84][85][86][87][88]. One of the benefit features of PLSR is that the relations between the dependent and independent variables can be concluded from the weights and regression coefficients of single dependent variable in the furthermost independent components [85,86]. Thus, it is possible to control which dependent variable strongly relates with independent variables [88].
Using multiple regression, changes in hydrological components between two land use maps in 2000 and projected 2032, respectively, were related to changes of land use to quantify the effect of individual land use at the basin scale. Dependent variables were the changes for eight land use classes (i.e., natural forest, woodland, bushland, grassland, water, wetland, cultivated land and built-up area) while independent variables (responses) were changes for four hydrological components (i.e., surface runoff, groundwater flow, evapotranspiration and water yield). The PLSR models were created to recognize the key land use types that control hydrological components. To reduce the drawback of over-fitting, the suitable number of components of individual PLSR model was determined by cross-validation to attain an optimum balance among the described variation in the response and the projecting ability of the model (goodness of prediction: Q 2 ). In PLSR modelling, the significance of a predictor for both the dependent and the independent variables is set by the variable importance for the projection (VIP), with large VIP values are the best significant for describing the dependent variable.
The percentage of variance and cross-validated goodness of prediction (Q 2 ) described for the dependent variables (hydrological components), together with the cross-validated root mean squared error as the change among the predicted and observed values of each individual pass, were analyzed for each model. The regression coefficients of the PLSR models were used to express the direction of the correlation between the changes in individual land use types and hydrological components. XLSTAT was used to perform the statistical analysis of the dataset [41].
where y is the dependent variable, b 0 is the intercept, x is the independent variables from 1 to i, and b is the coefficients of the x variables.

Accuracy Assessment and CA-Markov Validation
The accuracy results, including accuracy per specific land use class (i.e., producer accuracy (PA) and user accuracy (UA)), are presented in Table 3. The overall accuracy was 91.87% and 99.01% for the 2000 and 2016 maps, respectively. The Kappa coefficient was above 0.90 for the two classified images.
According to Anderson et al. [89], the level of accuracy in the classification of land use classes from remote sensing data must be at a minimum of 85%, while the Kappa coefficient must be higher than 80% [90]. Hence, the data set showed an excellent agreement with the image and the ground truths. CA-Markov validation was achieved, with a Receiver Operating Characteristic (ROC) value of 87.5%. The Kappa statistics values such as Kno (86.20%), Klocation (83.3%), KlocationStrata (89.8%) and Kstandard (81.9%) were also above 80%, which shows the excellent capacity of the model to simulate the 2032 land use [91,92].

Land Use Change Analysis
The land use situations of the Wami River Basin between 2000 and 2032 are shown in Table 4, Table 5 and Figure 2, while the spatial distribution of all land use is represented in Figure 3. During 2000-2016, a small amount of land use change happened in the basin; further, these changes are anticipated to continue in 2032. For example, grassland increased from 5.74% in 2000 to 14.25% in 2016 and is predicted to continue to increase to 19.30% by 2032 (Table 4). Cultivated land area expanded from 19.05% in 2000 to 26.34% in 2016 and continued increasing to 28.02% in 2032. Similarly, built-up areas increased between 2000 and 2016 (0.02% to 0.16%), and are expected to increase to 0.28% by 2032. Moreover, bushland increased from 44.48% in 2000 to 47.17% in 2016, but is expected to decrease to 44.42% in 2032. Nevertheless, the yearly expansion percentage of grassland until 2032 is higher than cultivated land, followed by bushland and built-up areas. Net loss was observed in natural forests, woodland, wetlands, and water. Natural forest decreased from 6.99% in 2000 to 2.09% in 2016 and continued to decrease to 1.02% in 2032 (Table 4). Woodland decreased from 19.93% in 2000 to 9.43% in 2016 and continued to decrease to 6.46% in 2032. Similarly, the extent of wetland decreased between 2000 and 2016 (3.50% to 0.54%), and is expected to decrease in to 0.47% in 2032. Moreover, water decreased from 0.19% in 2000 to 0.02% in 2016; it is expected to remain constant in the period up to 2032. The annual reduction rate of woodland until 2032 was greater than that of natural forests, followed by wetlands and water. These results are in agreement with a previous study in the basin [53].

Sensitive Parameters.
The sensitivity analysis results identified six sensitive parameters controlling the output variable ( Table 6). The parameters were adjusted to suit the ranges, until a satisfactory agreement was obtained between observed and simulated flows. Parameter values (upper and lower bound, and fitted value) used for model sensitivity analysis were found within the rage of the suggested SWAT user's manual [93]. The three most sensitive parameters were R_CN2.mgt, SOL_AWC.sol, and V_GW_DELAY.gw.

Sensitive Parameters.
The sensitivity analysis results identified six sensitive parameters controlling the output variable ( Table 6). The parameters were adjusted to suit the ranges, until a satisfactory agreement was

Sensitive Parameters.
The sensitivity analysis results identified six sensitive parameters controlling the output variable ( Table 6). The parameters were adjusted to suit the ranges, until a satisfactory agreement was obtained between observed and simulated flows. Parameter values (upper and lower bound, and fitted value) used for model sensitivity analysis were found within the rage of the suggested SWAT user's manual [93]. The three most sensitive parameters were R_CN2.mgt, SOL_AWC.sol, and V_GW_DELAY.gw.

Impacts of Land Use Changes in Hydrological Components
Investigation of the variations in the hydrological components influenced by land use changes showed effects on the components ( Table 8). The analysis focused on hydrological components including surface runoff, groundwater flow, evapotranspiration and water yield. The results are The achieved R 2 (0.93 for calibration and 0.83 for validation) values express good evenness between the observed and simulated data [72,74,76,83]. NSE above 0.60, PBIAS < 10% and RSR < 0.6 were also achieved while positive values of PBIAS indicate the under estimation of the model.

Impacts of Land Use Changes in Hydrological Components
Investigation of the variations in the hydrological components influenced by land use changes showed effects on the components ( Table 8). The analysis focused on hydrological components including surface runoff, groundwater flow, evapotranspiration and water yield. The results are presented in Table 7, from which it can be seen that the hydrological components were affected differently, while surface runoff was the most significant contributor to streamflow in the Wami River Basin. With land use 2000, the results for the hydrological components were as follows: surface runoff, 64.61 mm; groundwater flow, 87.45 mm; evapotranspiration, 511.8 mm; and water yield, 168.38 mm. Moreover, the results showed that, with land use 2016, the surface runoff, groundwater flow, evapotranspiration, and water yield were 68.84 mm, 86.76 mm, 508.7 mm, and 171.63 mm, respectively, while those for 2032 (projected) were 71.26 mm, 85.24 mm, 507.6 mm, and 171.77 mm, respectively. Furthermore, between 2000 and 2016, the surface runoff increased by 4.23 mm while the water yield increased by 3.25 mm, while groundwater decreased by 0.69 mm and evapotranspiration decreased by 3.1 mm over the same period. During the study period (2016-2032), an increasing trend was observed in surface runoff (by 2.42 mm) and water yield (by 1.14 mm), while a decreasing trend was observed in groundwater flow (by 1.52 mm) and evapotranspiration (by 1.1 mm).

Impacts of Individual Land Use Changes in Hydrological Components
The pairwise relationships of the eight land use types and hydrological components are presented in Table 9. Based on the results, nearly all land use types showed a strong relationship with the individual hydrological components. For example, natural forest and woodland both showed negative (but significant) correlations with water yield and surface runoff, and positive correlations with groundwater flow and ET. Bushland had a positive correlation with water yield, surface runoff, and groundwater, but a negative correlation with ET. Furthermore, it was observed that grassland, cultivated land, and built-up areas had significant positive correlations with water yield and surface runoff, while their correlations with groundwater and ET were negative. However, water and wetland had significant negative correlations with surface runoff and water yield, and significant positive correlations with groundwater and ET.  The results of the PLSR model for hydrological components in the Wami River Basin are given in Table 10 and a summary of the PLSR weights and variable importance of the projected values (VIP) of the hydrological components is given in Table 11. The PLSR model Q 2 , cumulative, and variation in the response values were above 0.60, signifying that the model was correctly predicted. The minimum Root Mean Press (predicted residual sum of squares) was used to obtain three components of water balance in the Wami River Basin. The variance of changes in water balance for the 1 st , 2 nd , and 3 rd components was explained in the cumulative model by 95.1%, 98.2%, and 99.8%, respectively (Table 10). The explained variance did not significantly increase with additional components (Table 11).  The first component was dominated by natural forest, woodland, water and wetland on the negative side, and by built-up area and cultivated land on the positive side. Moreover, the second component was dominated by natural forest, woodland, bushland, water and wetland on the positive side, and cultivated land and built-up area on the negative side. The last component was dominated by all land classes, where the negative side included grassland, built-up area and cultivated land, while the positive side included natural forest, woodland, bushland, water and wetland. The results of the relative importance of predictors of the model showed that bushland (VIP < 1) had less significance and that bushland and grassland were relatively less important in affecting water balance components, compared to other land use types.
Furthermore, the PLSR coefficients showed the effects of individual land use types on hydrological components in the periods from 2000 to 2032 in the Wami River Basin, which are given in Table 12. The regression coefficients indicated that natural forest, woodland, bushland, grassland, water and wetland influenced water yield and surface runoff negatively, while the effects of built-up area and cultivated land on these components were positive. Moreover, the regression coefficient results revealed that cultivated land and built-up area influenced groundwater flow and evapotranspiration negatively, while the effects of natural forest, woodland, bushland, grassland, water and wetland on these components were positive.

Discussion
The main objectives of this study were to investigate the effects of land use change on various elements of the hydrological processes in the Wami River Basin. The change of different land use types occurred in different directions and degrees during the study period. For example, natural resource decreased except for grassland while increasing changes were observed in cultivated land and built-up area. By analyzing land use changes, this paper indicated that the Wami River Basin was experiencing cultivated land and built-up area development during the study period. The land use pattern observed agrees with findings from Dickens [94] and Anderson [95] in which it was found that challenges for river conservation and management include expansion of irrigated agriculture and increasing domestic water supply withdrawal due to population increase. Furthermore, the basin was affected by the Southern Agricultural Growth Corridor of Tanzania Project with the aim of accelerated agricultural expansion and agricultural intensification [96]. Thus, irrigation is desirable and comes at the cost of already stressed fresh water supplies [97].
Simulating a land use change in Wami River Basin revealed the changes in basin hydrological processes. From the study, it is evident that land use has changed across the basin since 2000 and is expected to continue in 2032, while the effects of land use changes on the hydrological components were observed. The increase of bushland, grassland, built-up areas and cultivated land, accompanied by a decrease of the natural forests, woodland, water and wetlands during the study period increased water yield and surface runoff while reducing the groundwater flow and evapotranspiration. The ET is the key element in understanding the effect of change in land use on water production, and in this study a significant decrease was shown, which could be caused by increase in grassland [98]. The results agree with Neitsch et al.'s [93] statements, showing that the SWAT model computes the ET based on the maximum plant transpiration rate, canopy water evaporation intercepted and the maximum soil evaporation rate.
The PLSR results demonstrated that all land uses are significant and were relatively important in affecting water balance components, except for bushland and grassland. Therefore, in terms of water stability, the natural forest, woodland, wetland, cultivated land and built-up area should be regulated to a rational level to decrease the related impact on hydrological components. Despite relatively small changes in runoff generation due to the land use change, the findings suggest that increasing cultivated land and the built-up area in the Wami River Basin can cause severe threats in the future [55]. This will be more complicated with the current rapid population growth, climate change and economic growth within the basin [99]. The results can be used to identify and prioritize the implementation of best management practices in the most vulnerable areas to protect the natural resource in the basin [100]. Additionally, this study showed the potential of using a hydrologic modelling framework for a data-scarce basin in Tanzania. The use of Landsat images produced precise land use types in this analysis and showed the land use change in the basin. These results can be used to identify and highlight the application of best management practices in the most exposed areas.

Conclusions
Effects of changes in land use classes due to hydrology components in the Wami River Basin were assessed using hydrological modeling and PLSR. It was found that the major land use changes that affected hydrology components in the basin were expansion of cultivation land, built-up area and grassland, and decline in natural forests and woodland during the study period. The VIP values for these land-use types are greater than 1, and hence they are considered of great importance for the prediction of changes in the hydrology components. It can be concluded from the analysis that cultivation land, built-up area and grassland is directly proportional to the surface run-off and water yield, but inversely proportional to the groundwater flow and evapotranspiration. Furthermore, the decline in other land use classes is directly proportional to the groundwater flow and evapotranspiration, but inversely proportional to the surface run-off and water yield. These findings provide baseline information for decision-makers and stakeholders concerning land and water resources for better planning and management decisions in the basin. The approach used in this study has attributed contributions of changes in land use to hydrological components and might be useful to predict hydrological consequences of land use changes to another basins. Therefore, more studies should be undertaken to examine the effects of land use changes with different temporal and spatial scale in future research.
Author Contributions: The four authors contributed in a substantial way to the manuscript. All authors discussed the structure of the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research was funded by the Ministry of Water; the United Republic of Tanzania and the APC was funded by Dresden University of Technology (TU Dresden).