Application of SWAT to Assess the Effects of Land Use Change in the Murchison Bay Catchment in Uganda ()
1. Introduction
Land use change can lead to undesirable effects on ecosystems [1] . Land use and land cover changes are significant causes of water, soil and air pollution [2] , which negatively impacts the health of rivers within catchments. The effects of land use changes can be grouped into hydrologic, socio-economic, ecological and environmental [1] . Most river basins have undergone massive change over the past years due to various land use activities [3] . According to [4] , land use and land cover transformations could lead to change in flow pattern of rivers due to temporal changes in the distribution of runoff.
Land use and land cover changes are very common in developing countries whose economies are mainly dependent on agriculture and with rapid human population growth [5] . Vegetation removal to prepare land for agriculture leaves soil susceptible to massive increase in soil erosion by wind and water. This reduces the fertility of the soil rendering it unsuitable for agricultural purposes, as well as transport large volumes of nitrogen, phosphorus, and sediments to streams which can lead to various negative impacts such as increased sedimentation, turbidity, eutrophication and coastal hypoxia of wetlands and rivers. The use of agrochemicals such as herbicides, pesticides, and inorganic fertilizers in modern-day agriculture has hugely contributed to increased levels of pollution in surface water bodies as well as contamination of groundwater through runoff and by way of leaching. This pollution in most cases is toxic to aquatic life [2] and humans. These impacts not only affect the immediate area but their effects can extend to distant regions [3] .
The Murchison Bay catchment in Uganda, which is a major contributor of water for the people in the capital city Kampala, has experienced a lot of land use-land cover transformations over the years. The change of land use from its natural vegetation has caused environmental degradation resulting from inadequate facilities for sewage and sanitation, poor drainage system, and increased industrial pollution and urban agriculture [6] . For instance, different studies [7] - [13] have reported the rising levels of pollution and its increasing impact as a result of expansion of Kampala, wetland encroachment, poor agricultural practices and deforestation in parts of this catchment. Agricultural production in the immediate environs of Kampala has led to rapid devastation of green belts causing deforestation and has led to reduced soil water retention capacity and accelerated erosion [6] . Indiscriminate drainage of wetlands in the city environs for cultivation has also affected the water table. Drainage of wetlands has led them to lose their function to control floods, filter effluents and purification of waters before been discharged into Lake Victoria. As a result of these changes, the cost of treating water has become the major bottleneck to National Water and Sewerage Corporation (NWSC), an institution in charge of treating and supplying water to Kampala City inhabitants [8] .
According to Li et al. [14] , management of water resources and land use patterns are inherently linked. At the catchment level, the type of land use, such as agricultural production, forests, human settlements, industrial and commercial centers, influences the quantity and quality of available water. Continuous monitoring of water quality parameters is expensive, time consuming, and at times difficult [15] . A number of hydrologic and water quality simulations models have been developed to examine the impact of different land use practices on water resources. The results of these models are used to help policy and decision makers determine cost-effective land management and conservation practices to sustain water resources for current and future populations under a changing climate [16] [17] [18] . One of such models is the Soil and Water Assessment Tool (SWAT).
SWAT is a river basin scale and continuous time model that operates on daily time step. This model is designed to quantify and predict the impact of land management practices over long periods on water, sediment and agricultural chemical yields in large complex watersheds with varying soils, land use and management conditions [19] [20] [21] . Expanded routing and pollutant transport abilities including reservoir, point source, pond, wetland, septic tank effects and enhanced sediment routing routines have also been incorporated into the SWAT model. SWAT is a versatile model currently used worldwide to evaluate water quality and hydrological concerns in a number of varying watershed scales and environmental conditions [22] for abreast policy decision making and effective watershed management. Details of SWAT are in [21] [23] .
Assessing the impacts of land use and land cover changes on hydrology forms the foundation for managing watershed and restoring its ecology [24] . In this study, the SWAT model was used to simulate streamflow and also to estimate the amount of sediment yield and nutrients loss from the Murchison Bay Catchment as a result of land use-land cover changes.
2. Materials and Methods
2.1. Study Area
The Murchison Bay is an extension of Lake Victoria situated in the south-east of Kampala which lies between latitudes 00˚10'00"N - 00˚30'00"N and longitudes 32˚35'00'E - 32˚50'00"E with average elevation of 1224 m above sea level. The Bay which is divided into Inner and Outer Bay covers an area of about 60 km2. The Inner Murchison Bay is comparatively a semi-enclosed small water body with an area of 18.4 km2 and length of 5.6 km off the main part of Lake Victoria. The Murchison Bay watershed covers an average combined catchment area of 282 km2 which consists of about 20% swamps and 80% upland [7] [10] with the largest catchment drainage inlet being the Nakivubo channel which passes through Kampala City. This study covers the largest watershed of the Murchison Bay that drains through the Nakivubo Channel (Figure 1). It covers an area of about 40.9 km2. The Nakivubo channel runs from the central populated Kampala district passing through slums, industrial settlements and markets before discharging its water to the Lake Victoria at Murchison Bay. The water in the channel is a mixture of secondary effluents from NWSC sewage treatment at Bugolobi within Kampala city and heavily polluted untreated wastewater from other parts of the city. There used to be a predominant papyrus wetland filtration system but it has been drained and turned into agricultural land or developed for commercial, industrial or residential purposes.
The catchment is within the equatorial belt, and has a moist sub-humid climate. It receives a bi-seasonal rainfall in the periods of March to May and September to November [25] . Mean monthly precipitation ranges from 24 mm in January to about 154
Figure 1. Map of Uganda showing the Murchison Bay Watershed.
mm in October with mean precipitation of about 1293 mm/year.
The soils of the Nakivubo swamp area are alluvial and lacustrine sands, silts, and clays overlying granite gneisses. Alluvial soils within the swamp range from semi-liquid organic material in the very upper layers of the emergent vegetation zones, through reddish ferruginous loams to clays [25] .
2.2. Maintaining the Integrity of the Specifications
The input data used in the study included the Digital Elevation Model (DEM) of Central Uganda, land use-cover map of the Murchison Bay catchment, soil and hydro- meteorological data of the study area.
2.2.1. Digital Elevation Model (DEM)
A DEM for Central Region of Uganda at a resolution of 1 arc second (30 m × 30 m) was retrieved from the United States Geological Surveys (USGS) website (http://gdex.cr.usgs.gov/gdex/). The DEM was used for automatic delineation of the catchment and also to define the stream network and determine sub-basin parameters such as slope [26] .
2.2.2. Land Use and Land Cover Data
Land use cover maps of 1995 and 2003 were used in this study. These were obtained by visual interpretation of remotely sensed images retrieved from USGS Landsat ETM/TM satellites (http://www.earthexplorer.usgs.gov/) in Path 171 and Row 60 at a spatial resolution of 30 m. The images were classified into five land use-cover types in accordance to Anderson et al. [27] Level I generalized classification system using maximum likelihood supervised classification tool in ArcGIS. These are built-up land, wetland, forestland, agricultural land and open water bodies. Table 1 shows the characteristics of land use and land cover in the Murchison Bay catchment.
The major land use and land covers in 1995 were agricultural land (32.42%), forestland (31.15%) and built-up land (26.53%). However, by 2003, built-up land became the most dominated land cover followed by agricultural land while forest land drastically declined. Figure 2 shows the various land use-land covers identified in the Murchison Bay catchment. The main cause of these changes can be attributed to rapid population growth around the catchment which increased by close to 50% between 1991 and 2002 [28] . According to [10] , the rapid growth of Kampala is taking place at the expense of the environment which is highly correlated to the pollution level observed in the Lake Victoria.
2.2.3. Soil Data
The soil map was obtained from Land and Water Resource, FAO soil database [29] and the description of the soils from [30] . The dominant soil associations in the Murchison Bay catchment are acric ferrasols, orthic ferrasols and dystrics gleysols. The soil physical and chemical properties included were texture, soil hydrologic groups, maximum
Table 1. Characteristics of land use and land covers in the Murchison Bay catchment.
Figure 2. Land use and land cover classified map of the Murchison Bay Catchment.
rooting depth, fraction of porosity, moist bulk density, available water capacity, saturated hydraulic conductivity, organic carbon content, electrical conductivity, USLE equation soil erodibility factor, sand, silt and rock fragment contents.
2.2.4. Hydro-Meteorological Data
Hydrological simulation with the SWAT model requires weather data consisting of precipitation, maximum and minimum temperature, relative humidity, solar radiation and wind speed. These parameters were obtained from the National Centres for Environmental Prediction (NCEP) (http://globalweather.tamu.edu/) and NASA agro-cli- matology (http://power.larc.nasa.gov/cgi-bin/cgiwrap/solar/agro.cgi) for the period of 1995 to 2010. The weather data was collected from three stations outside and one station within the study area. Stream discharge data was obtained from the Directorate of Water Resource Management (DWRM) under Uganda’s Ministry of Water and Environment (MWE) at Entebbe. There is only one stream flow gauge in the Murchison Bay catchment located at the Railway Bridge culvert which is very close to the outlet of the delineated catchment at 32.63˚E and 0.30˚N. The observed stream flow data for the period from 1997 to 2007 was used for calibration and validation of the model, however, there were data gaps.
2.3. The SWAT Model Setup
The model uses readily available input data such as Digital Elevation Model (DEM), land use data, soil data and climatic data as described in Section 2.2. In this study, modelling of the hydrological process was carried out using the extension of SWAT for ArcGIS software called ArcSWAT [31] .
ArcSWATv2012.10.1.18 was downloaded from the website for the model (http://swat.tamu.edu/software/arcswat/) and installed in ArcGISv10.1. In the model setup, the first step was to delineate the catchment using the DEM into several connected sub-basins. The sub-basins were further divided into smaller units called hydrologic response units (HRUs). HRUs are lumped land areas within the sub-basin that are comprised of unique land cover, soil, slope and management combinations [21] . A total of 23 sub-basins and 247 HRUs were created. The HRUs were created by defining the thresholds of land use over sub-basin area at 5%, soil class over land use area at 5% and slope class over soil area at 5% using the multiple HRUs definition. The model was run on daily time step for a period of 16 years from 1995 to 2010 with a warm up period of two years.
2.4. Sensitivity Analysis
Sensitivity analysis was performed to choose the most sensitive flow parameters [32] [33] [34] [35] that influence the catchment represented by SWAT to be used for calibration. This was achieved using the global sensitivity approach in semi-automated Sequential Uncertainty Fitting (SUFI2) algorithm. The global sensitivity analysis method takes into consideration, the sensitivity of one parameter relative to the other in order to give their statistical significances [36] . The t-statistics and p-values of the parameters were used to rank to the different parameters considered to influence flow and the final selection done based on the significance of the ranked values. Table 2 shows stream
Table 2. Sensitivity analysis parameters.
flow parameters that were tested for their sensitivity. These are useful in estimating the amount of flow from a catchment [37] .
2.5. Calibration and Validation of the Model
Calibration was accomplished by comparing the output of the SWAT model with the observed data at the same conditions [19] [38] . For calibration and validation, the semi-automated Sequential Uncertainty Fitting (SUFI-2) calibration method within the SWAT Calibration and Uncertainty Procedures (SWAT-CUP) was used. SWAT-CUP version 5.1.6.2 was used. SWAT-CUP [16] is a stand-alone calibration program developed for SWAT that operates on Latin Hypercupe sampling procedures. Due to lack of observed data for sediment and nutrients, the model was calibrated and validated only for streamflow. The model was calibrated with observed daily discharge data for the period of 1997 to 2002 and validated for the period of 2004 to 2008 since there was no observed data for 2003. Though, there were gaps in the observed data, the challenge was addressed by writing it in a format suggested by [36] that SWAT-CUP tool can read.
2.6. Evaluation of Model Performance
To evaluate performance of the model during calibration and validation, statistical measures as well as graphical representations at daily time step were used. This was employed to confirm the relationship between simulated or predicted values and observed values [1] and to verify the robustness of the model [37] . Four statistical measures were employed. They are the coefficient of determination (R2), the Nash-Sutcliffe efficiency (NSE), the percent bias (PBIAS) and the ratio of mean squared error to the standard deviation of the measured data (RSR) [26] [37] [39] . Equations (1)-(4) were used to determine NSE, PBIAS, RSR and R2, respectively. Other details of these measures such as their utility and satisfactory range of values are explained by Moriasi et al. [40] .
(1)
(2)
(3)
(4)
where Oi is the observed daily discharge, Si is the simulated daily discharge, is the average measured discharge, is the average simulated discharge and n is the number of observations.
3. Results and Discussion
3.1. Sensitivity Analysis of Stream Flow Parameters Results
The global sensitivity analysis of 16 flow parameters (Table 2) showed that, only seven were very sensitive to flow. The rankings of the flow parameters are presented in Table 3 while the fitted values for the most sensitive parameters are indicated in Table 4. The most sensitive parameter was the SCS runoff curve number (CN2). The curve number estimates runoff based on the relationship between precipitation, hydrologic soil group and land uses. Other researchers [41] [42] [43] [44] have also found the SCS curve number to be the most sensitive streamflow parameter in modelling hydrology in their studies. The other sensitive parameters included the Manning’s n value for the main channel (CH_N2), effective hydraulic conductivity in the main channel alluvium (CH_K2) and groundwater delay time (GW_DELAY), baseflow alpha factor (ALPHA_BF), groundwater “revap” coefficient and the soil evaporation compensation factor (ESCO). Although, the rest of the parameters were found not to be sensitive to flow in the catchment as their p-values were greater than 5%, other studies [37] [39] [45] [46] [47] have found some of them to be sensitive in their studies. This is expected as conditions such as land use and land covers, soil characteristics and climatic factors vary from one catchment to the other.
3.2. SWAT Model Calibration and Validation Results
The graphical results of the model streamflow simulation performance during the
Table 3. Sensitivity rankings of streamflow parameters in the Murchison Bay catchment.
Table 4. List of sensitive parameters and their calibrated values.
The extension (e.g., .mgt, .rte) refers to the SWAT input file where the parameter occurs. The qualifier (r__) refers to relative change in the parameter where the value from the SWAT database is multiplied by 1 plus a factor in the given range, while the qualifier (v__) refers to the substitution of a parameter by a value from the given range.
calibration and validation periods are shown in Figures 3-6. The year 2003 was not incorporated in the validation period since there were no recorded observed daily discharges in that year. In general, graphical results during calibration (Figure 3 and Figure 5) and validation (Figure 4 and Figure 6) indicated adequate calibration and validation over the range of streamflow discharge, although the calibration results showed a better match than the validation results. Statistical model streamflow discharge simulation performance results during both calibration and validation periods are presented in Table 5. Generally, the statistical results showed a satisfactory performance of the model [40] . Just as with the graphical results, the model simulated streamflow discharge better during the calibration period than the validation period. The relatively low statistical measures during the validation can be attributed to the quality of the observed data
Figure 3. Observed and simulated daily discharge for the Murchison Bay catchment during calibration (1997-2002).
Figure 4. Observed and simulated daily discharge for the Murchison Bay catchment during the validation period (2004-2007).
Figure 5. Scatter plot of observed and simulated discharge for the calibration period (1997-2002).
Figure 6. Scatter plot of observed and simulated discharge for the validation period (2004-2007).
Table 5. Daily time step SWAT calibration and validation statistics.
used as there were many missing days. The results also showed that, the model slightly underestimated flow during both the calibration and validation period as indicated by PBIAS values of approximately 8% and 23%, respectively. According to Moriasi et al. [40] , a PBIAS greater than zero is an indication that the model underestimated flow. During the calibration period (1997 to 2002), the observed and simulated mean daily flows were 3.1 m3/s and 2.91 m3/s, respectively, while during the validation stage (2004- 2007), the observed and simulated mean daily flows were 1.94 m3/s and 1.50 m3/s, respectively.
Also, the underestimations could be attributed to uncertainties that exist in the catchment such as the wetland processes, quality of input data and wastewater dis- charge from point sources which could not be accounted for during the simulation. These plausible reasons including the limitations of using the SCS curve number (CN2) to compute estimate water budget componentshave also been noted by Abbaspour et al. [32] [48] and Betrie et al. [37] to be the cause of over and under prediction of flow. This observation agrees with that of Mutenyo et al. [43] and Qiu et al. [47] who recorded an underestimated streamflow discharge by the SWAT model. However, Qiu et al. [47] suggested that underestimation or overestimation [37] of streamflow discharge by the SWAT model is partly due to the use of curve number,which cannot give accurate pre- diction of runoff for days with several storms. Mutenyo et al. [43] attributed inaccurate simulation of streamflow discharge by the SWAT model to insufficient rain data. It is obvious that, if all these uncertainties are minimised, a well calibrated SWAT model can efficiently predict flow in the catchment for any management purpose.
3.3. Impact of Land Use and Cover Change in the Murchison Bay Catchment
To evaluate the effects of land use and cover change, the SWAT model calibrated and validated for streamflow discharge was used to simulate different land use and cover scenarios on runoff and sediment yield. Land use covers at two different years thus 1995 and 2003 were used as input variables in the SWAT model in order to compare output as a result of their differences in land covers. The simulated results of surface runoff and sediment yield in the catchment of the Murchison Bay under different land use and land cover scenarios using the SWAT model are shown in Table 6. The results showed an increase in runoff and average upland sediment yield, from the catchment except maximum upland sediment yield which reduced very slightly between the two periods. Surface runoff increased from 101 mm/yr in 1995 to 128 mm/yr in 2003 representing 26.7%. This is attributed to the drastic land cover changes that occurred in the catchment. It highlights the effect of changes in land use and cover on hydrology in the catchment. Within the period, there was a sharp decline of forestland from 31.15% in 1995 to 13.91% while built-up land shot from 26.53% to 39.09%.
The observation agrees with the assertion by Nie et al. [24] and Tang et al. [49] that, an increase in urban set-up increases impervious surfaces and declining forest lands, which accelerates surface runoff. Similarly, [50] reported that an increase in annual surface runoff increased when forest lands were converted to urban areas and decreased during when forestland gained a significant expansion. Maalim et al. [51] also made similar observation, areas dominated by urban development has higher runoff and lower runoff in forested areas. On the other hand, different studies [14] [44] [52] have reported that an expansion in forestland and grasslands lead to a decrease in runoff. Although, there is general agreement among researchers that changes in streamflow is due to different types of forestry activities, such as afforestation that may lead to lower runoff generation, others suggest the amount of precipitation [44] [52] and climate change [53] partly play a role. The increase in built-up land within the period might have possibly created impervious layers which reduced infiltration and percolation of water to the shallow aquifers that caused surface runoff to increase. According to Diyer et al. [54] , when industrial or urban development takes place, the new land use changes how water is transported and stored. The combination of constraints connected to compacted lands and impervious surfaces such as parking lots, roads, sidewalks, and roofs leads to reduced infiltration. This causes, for instance, reduced water quality, increased volume and velocity of runoff, increased occurrence and severity of floods, and loss of storage capacity and runoff water in natural vegetation.
Surface Rrunoff is the carrier of all other components such as sediment, nutrient, pesticides, bacteria, agricultural waste, heavy metals, industrial solid and liquid waste that undesirably affects water quality. As surface runoff increases, it is evident that all the aforementioned components will correspondingly increase depending on the activities
Table 6. Estimated runoff and sediment under different land use and covers.
taking place in the catchment. However, the simulated results showed there was not much change in upland sediment yield over the period with increasing runoff. Upland sediment yield rather reduced by at least 0.50% between the two periods. This could still be attributed to the dynamics of land cover that happened at that time. Nevertheless, average upland sediment yield was about 20.50% higher in 2003 than that of 1995. This means that, some sub-basins within the catchment might have been exposed to massive degradation such as construction works and intensive agricultural activities than other sub-basins which resulted in higher sediment yield as runoff increased.
At a rate of 26.54 tons/ha/yr and 26.40 tons/ha/yr which translates to about 108,557 tons/yr and 107,984 tons/yr respectively, the sediment yield estimated in the catchment is greater than the results reported by Kimwaga et al. [55] . In the Simiyu catchment of the Lake Victoria basin, which is much larger than the Murchison Bay catchment, [55] reported a modelled estimate of sediment yield of 98,467 tons/yr though the actual measured sediment yield was 2,075,144 tons/yr. The discrepancy was attributed to the data used to run the model.
The average estimates of sediment yield from upland and runoff from the landscape as well as possible mass loadings from point sources in the Murchison Bay catchment is a clear indication of the deterioration level of water quality that ends up into the Lake Victoria. This is supported by [10] situation analysis report that increasing human population and all their associated economic and social activities have accelerated the rate of delivery of nutrients and caused eutrophication in Lake Victoria. However, it is anticipated that the values of the components estimated have changed tremendously as the catchment is now condensed with more settlements, industries and increased in population of Kampala city from 1.2 million in 2002 to 1.5 million in 2014 [28] .
4. Conclusion and Recommendation
In this study, the SWAT model was used to assess the impact of land use and cover change on the Murchison Bay catchment. The calibrated SWAT model satisfactorily simulated streamflow with satisfactory statistical measures results. Average runoff and sediment yield loss from the catchment was also estimated. The estimations showed that because of rapid changes that occurred over the past, it has led to an increased runoff and average upland yield by 26.73% and 20.50% respectively while maximum upland sediment yield relatively reduced by 0.50% because of increasing impervious surfaces in the catchment. The effects of increasing surface runoff are that, water quality and crop yield in addition to aquatic life are negatively affected as it transports nutrients and other pollutants. This may also explain severe and frequent flooding, which has been experienced in the catchment lately. Therefore, pragmatic strategies for sustainable water and environmental resources in the catchment should be developed for sustainable water and environmental resources in the catchment. This is the first time the SWAT model has been used in the Murchison Bay catchment. The results showed that, if all uncertainties are minimised, a well calibrated SWAT model can generate reasonable hydrologic simulation results in relation to land use, which is useful to water and environmental resources managers and policy and decision makers.
Acknowledgements
The authors wish to thank the METEGA Project and RUFORUM for their financial support.