Impact of Land-Use Changes on Sediment Load and Capacity Reduction of Lake Ziway, Ethiopia ()
1. Introduction
Lakes are an essential part of ecosystems that can sustain a healthy aquatic life, and provide human socio-economic needs. It is the responsibility of human being to continue to practice stewardship in our lakes by keeping them healthy for all. Among the parameters, land use and land cover change has been a factor that alters the hydrologic response of the watersheds influencing stream flow and sediment yield changes on lakes [1]. Similarly, agricultural practices and Land use patterns have changed over time due to economic benefits, and these changes are also contributing to runoff [2]. [3] sate that Land use is one of the most critical factors that affect runoff, Evapotranspiration, and surface erosion in a watershed. Land use change occurs because of a natural process and human actions. The central highlands of Ethiopia are characterized as a region of high rates of land degradation and soil erosion [4]. The effects of land use on groundwater recharge and surface runoff and how these are affected by land use changes are of interest for sustainable water resources management [5].
The land is an indispensable natural resource that living and non-living things are prevailing together. Land use is how human beings employ the land and its resources, like for agriculture, urban development, grazing, industrialization, logging, mining. In contrast, Land cover refers to the biophysical characteristics of earth’s surface, including the distribution of vegetation, water, soil, and other physical features of the land [6]. Expansion and intensification of agriculture, the growth of urban areas and extraction of timber, and others are likely accelerated the changes in the land-use over the last decades to satisfy the demands of an increasing population [7].
Lake sedimentation is one of the factors reducing life in the lake with the sediments being transferred from the surrounding catchment by rivers and streams draining from different land-uses. Upstream sediment flows are accelerated significantly beyond natural conditions due to unsuitable agricultural practices in some areas and the rapid conversion of rural lands into urban and suburban Land uses in other areas. According to [8], clearing of forests enhanced surface runoff loaded with suspended sediments into the lake.
Lake Ziway, which is the largest freshwater lake in the Ethiopian Rift Valley, is fed principally by the Meki & Katar Rivers sub-watershed. There have been many activities, over the last century, have modified the land-use in the region [9]. Irrigated agriculture is a common practice around Lake Ziway, by pumping and diverting water from the lake and from the rivers that flow into the lake. Due to this, the LAND-USE change in the basin accelerates the impacts on the streamflow and sediment yield. In turn, the sediment which transports from the watershed affects water depth, level, volume, surface area, and it degrades the water quality of the Lake Ziway. Therefore, this research aimed to quantify the sediment load inter to the Lake and the reduction of lake water depth due to the changes of land-uses at the watershed area of the lake.
2. Materials and Methods
2.1. Description of Study Area
Lake Ziway watershed located in the central rift valley (CRV) of Ethiopia encompasses an area of 7148 km2, falling between gradients 7˚22'36"N and 8˚18'21"N latitude and 37˚53'40"E and 39˚28'9"E longitude (Figure 1). It is at 160 km south of the capital-Addis Ababa. The watersheds of the two principal rivers, Meki and katar Rivers, feeding Lake Ziway lay in two administrative regions of Ethiopia, 73.6% in Oromia National Regional State (ONRS) and the remaining part in Southern Nation Nationalities and People Regions (SNNPRS). The watershed stretches from edges of the Gurage Mountains in the west to the Arsi Mountains in the east [8]. The watershed includes the rift floor, two escarpment areas, two major river inlets the Katar and Meki Rivers, and one river outlet the Bulbula River. The Meki River originates in the highlands of Guraghe and travels a distance of about 100 km from the highlands at an altitude of 3600 m to 1636 m before draining into the Ziway Lake. The total catchment area of the river near Meki town is 2433 km2. The Katar River is originating in the highlands of Arsi has a drainage area of about 3400 km2. The water level of Ziway Lake influences by the outflow to the Bulbula River.
2.2. Land Use/Cover (LAND-USE) Classification
The changes of Land use/cover in the watersheds captured from the Landsat image data from USGS/Earth-Explore (https://earthexplorer.usgs.gov) public sources data for three different decadal as Land use/cover map of 1988 for the period (1988-1998), Land use/cover map of 2002 for the period (1998-2008), and
Land use/cover map of 2015 for the period (2008-2018) (Table 1). Maximum likelihood supervised image classification with ERDAS and ArcGIS interface used to classify images in different Land use/cover classes. LAND-USE maps for different scenarios years were considered to model to evaluate the impact of landuse dynamics on sediment loading to the Lake Ziway.
2.3. Model Setup
Soil and Water Assessment Tool (SWAT), as a hydrological model, was set with its ArcSWAT interface to evaluate the land-use changes on the sediment load Lake Ziway, in Ethiopia. The hydrological responses unit (HRU) of the model established using the topographic information extracted from 30 m resolution digital elevation model (DEM) the project area, and soil information extracted from the national geospatial soil data [10] as static data and Land use/cover information as a change parameter characterized with three decades. The climate data, including rainfall, temperature, sunshine hours, wind speed, and relative humidity for thirteen meteorological stations, which distributed inside or around the basin, collected from the National Meteorological Agency (NMA), Ethiopia. Besides, the stream flow data for three hydrological gauging stations on the three rivers, Meki, Katar and Bulbula collected from Ministry of Water, Irrigation, and Electricity (MoWIE) and used for calibration and validation of the model.
2.4. Sediment Analysis
Soil erosion involves detachment soil particles and transportation by runoff with subsequent deposition of the sediments in river courses and land depressions [11]. Sediment is detached from soil surface both by the raindrop impact and the shearing force of flowing water. SWAT model calculates the soil erosion and sediment yield caused by rainfall and runoff within each HRU with the Modified Universal Soil Loss Equation (MUSLE). It is determined using Equation [12]. Sediment is one of the most critical data in the hydrological study that is an estimation of sediment transport. The rate of sediment transport is given in terms of sediment discharge or as a flux-averaged concentration [13], which mean the mass of sedimentary material passes with the flow of water at a given cross section in unit time. It can be expressed in terms of weight or volume per time.
Table 1. Summary of spatial datasets used in this study (Source: https://www.usgs.gov).
2.5. Sediment Rating Curve
Due to the lack of long term observed sediment data in the study area, flow-sediment rating curve established by Alemu et al. [14] for Meki at Meki station, Katar at Abura station and Bulbula River station used to generate the long-term sediment yield to estimate the sediment yield discharged to the lake.
3. Results
3.1. Land Use Land Cover Changes
Spatial and temporal analyses carried out to assess the land-use changes over the catchment of Lake Ziway using ERDAS imagine 14 and ArcGIS 10.4 student license. The land-use changes were evaluated based on three representative Landsat images as 1988 for the period of (1988-1998), 2002 for the periods of (1998-2008), and 2015 for the periods of (2008-2018). Each year image classified into 7 basic Land use classes to detect the changes and its classification accuracy evaluated within confusion matrix and kappa coefficient using available ground truth data (Table 2). Therefore the bareland, rangeland, wetlands, and water bodies showed a decreasing rate in rank order and urban area, and agricultural lands showed an increasing rate (Table 3, Figure 2).
3.2. Sediment Yield Modeling
The amount of sediment yield from Lake Ziway watersheds simulated using a spatially semi-distributed SWAT model with five sediment parameters, which identified as the most sensitive parameters, significantly affect sediment yield (Table 4).
Table 2. Land cover classification accuracy over the confusion matrix.
Table 3. Land use/cover change over the last 3 decade in Lake ziway catchment.
Table 4. Results of sensitivity analysis for sediment parameters.
Figure 2. LAND-USE map of the study area.
3.3. Sediment Yield Calibration
The three years period from (1993-1996) for 1988 LAND-USE and the five years from (1999-2003) for 2002 LAND-USE used for the model calibration at both stations, and adjustment of the parameter values carried out within the allowable ranges recommended by SWAT model developers in [15] [16]. As shown in Table 5 and Figure 3, the results of the calibration analysis values indicate that the model is a well-performing trend of stream sediment flow and good correlation agreement with the observed mean in the watershed since R2 > 0.60 and NSE > 0.50 [17].
3.4. Sediment Yield Validation
The validation of the sediment yield of the watershed was carried out for the years January 1, 1996, to December 31, 1999, for 1988LAND-USE and January 1, 2003, to December 2007 for 2002 LAND-USE without further adjustment of the parameters. As shown below in Table 6 and Figure 4 the results of the validation, the model is a good performing trend of sediment yield. The land-use change effects on sediment yield of the watershed evaluated using the validated
Table 5. The calibration period performances and sediment yield variability.
Table 6. The validation period performances and sediment yield variability.
Figure 3. Monthly calibrated sediment yield A-Meki and B-Katar for 1988LAND-USE and C-Meki and D-Katar for 2002LAND-USE.
sediment yield results for the three different Land use/cover changes. The annual sediment yield of Lake Ziway watershed was increased from year to year due to Land use cover changes.
Figure 4. Monthly validated sediment yield A—Meki and B—Katar for 1988LAND-USE and C—Meki and D—Katar for 2002LAND-USE.
3.5. Sediment Yield over the Two Gauged Watersheds
The validated Soil and Water Assessment Tools (SWAT) model was run annually for twenty-eight years (1990-2018), to get annual sediment yield spatially within a sub-basins level using three reference Land use maps changing with decadal periods. As resulted in an increase of sediment yield by 72.759 t/km2 for the first period (1988-2002) and by 45.434 t/km2 for the second period (2002-2015) in the Meki sub-watershed, and also an increase of sediment yield by 82.575 t/km2 and 58.623 t/km2 for the first and second period, respectively, in the Katar sub-watershed observed. Moreover, for over the study period (1990-2018), there was an increase of sediment yield by 118.192 t/km2 and 141.198 t/km2 in the Meki and Katar sub watershed, respectively (Table 7). Average annual sediment yield for each sub-basin predicted by the model used to map and identify the significant sediment source area within the two watersheds (Figure 5). The total average annual sediment yield rate from both sub-basins is 3.59 t/ha/yr (2.57 × 106 tons/yr) for the period (1988-1998), 4.36 t/ha/yr (3.12 × 106 tons/yr) for the period (1998-2002), and 4.89 t/ha/yr (3.49 × 106 tons/yr) for the period (2008-2018), which indicates the increasing trend of sediment yield in the Lake Ziway due to Land use/cover change (Table 8).
Table 7. Mean annual sediment yield in both sub basins of Lake Ziway Watershed.
Table 8. Summary of sediment yield (tons/ha/yr) for different land use/cover changes.
Figure 5. Spatial variability of sediment yield in Meki and Katar by three decades LAND-USES.
3.6. Sediment Deposition in Lake Ziway Due to LAND-USE Change
Soil erosion decreases soil fertility of the uplands and causes siltation of lakes and reservoirs [18]. The sediment deposition in the lake has to computed based on the sediment budget which is the net sediment volume entering into the lake as a combination of suspended and bed-load sediment yield entering into the lake through the gauged (Meki and Katar) and ungauged river system mines the suspended sediment yield leaving the Lake through Bulbula River [19] [20]. Therefore, net annual sediment deposition in the lake due to Land use/cover of 1988, 2002, and 2015 is computed as 1.83 × 106 tons/yr, 2.2 × 106 tons/yr, and 2.48 × 106 tons/yr, respectively. Dividing this sediment mass by calculated bulk density of 1.22 ton/m3 for the Lake Ziway Watershed, the volumetric sediment deposition rate is 1.5 × 106 m3 per year, 1.81 × 106 m3 per year, and 2.033 × 106 m3 per year for the Land use/cover of 1988, 2002, and 2015, respectively. The incremental changes in the rate of the deposit volume clearly show the contribution of the land-use changes on the sediment deposition rate in the Lake Ziway.
4. Discussion
The growth of population and human activities as the primary factor transforms the natural land-use system [21], which dominated by vegetation coverage according to the climatic regimes, to the functional role of land for economic activities [22]. It is justified in this study as the land-use of Lake Ziway watershed changes through the increase of economic land-use types; agricultural land, urban and built area increases, while the decease of natural land-use types like rangeland, water bodies, and forest lands over the last thirteen years. Land use and cover features play a significant role in the erosion and sedimentation process of a catchment. Therefore, land-use changes from one period to another have shown an increasing trend in the sediment yield of the Lake Ziway watershed. Changes in human Land use have increased soil erosion rates globally [23]. Increasing sediment deposits in the lakes is one of the most essential and prevalent anthropogenic impacts on aquatic ecosystems. Examination of the sediment deposition is part of the approaches in understanding ecological and climatic changes and other anthropogenic impacts in lakes [24]. The gradual increase of sediment thickness through time may threaten the ecosystem and very existence of a Lake. The net deposition rate of sediment in the lake Ziway, were 1.5 mcm/yr, 1.81 mcm/yr, and 2.033 mcm/yr for the land-use of 1988, 2002, and 2015, respectively. The incremental changes in the rate of the deposit volume clearly show the contribution of the land-use changes on the sediment load in the Lake Ziway. The cumulative sediment deposition over the study period (1988-2018) also estimated as 55 mcm. The bathymetric survey of Lake Ziway, which did in 2006 [25], fixed the lower depth at 1627 m.a.s.l., with zero area and volume. The cumulative sediment volume deposited over the last 12 years after the bathymetric survey was 25.76 mcm, which able to cover 14.33 km2 lake areas and raise the lower depth of the lake to the level of 1631.3 m.a.s.l. It means the cumulative sediment depth change over the last 12 years became 4.3 m, which gave 36 cm/yr an average incremental depth changes on the lake.
5. Conclusions
The impacts of the land-use change on the spatial distribution of sediment yield evaluated using ERDAS and ArcGIS integrated SWAT Model. The accuracy of the three land-use maps developed as representative of three periods (1988-1998), (1998-2008), and (2008-2018) from TM, ETM+, and Ls8 satellite images, checked using the Confusion Matrix with a kappa coefficient of 0.79, 0.82, and 0.85, respectively. From the Land use and land cover change analysis in the study period, Lake Ziway watersheds were changed from its dominant rangeland and forest land coverage to agricultural land and urban/built area coverage over the last three decades.
Following the land-use changes, the changes in sediment load over the lake also evaluated in three decades. The total average annual sediment yields were 3.59 t/ha/yr for the period (1988-1998), 4.36 t/ha/yr for the period (1998-2008), and 4.89 t/ha/yr for the period (2008-2018). The increasing trend of sediment yield in the Lake Ziway watershed indicated the impact of Land use/cover change in the basin. Similarly, the net volume of deposited sediment in the lake also showed an incremental trend with the land-use changes as 1.5 mcm/yr, 1.81 mcm/yr, and 2.033 mcm/yr for the period of 1988-1998, 1999-2008, and 2009-2018, respectively. The sediment deposition over the lake induced reduction of depth and water holding capacity of the lake also by 4.3 m and 25.76 MCM for the last 12 years.
Besides, the study concluded that there is an annual average depth reduction of 36 cm over the lake area due to land-use changes over its catchment area, which has significant impacts on the capacity and existence of the lake unless it mitigate with appropriate watershed management measures.
Acknowledgements
The authors wish to thank MoWIE and NMA of Ethiopia for support regarding hydro-meteorology data. We are also grateful to the school of Civil and Environmental Engineering in Addis Ababa University Institute of Technology to facilitate all things to conduct this research.