[go: up one dir, main page]
Next Article in Journal
Relative Effect of Location Alternatives on Urban Hydrology. The Case of Greater Port-Harcourt Watershed, Niger Delta
Previous Article in Journal
Flood Hazard Analysis of Proposed Regulator on Shatt Al-Arab River
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Hydrology
Volume 6
Issue 3
10.3390/hydrology6030081
hydrology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Impacts of Climate Change on Sediment Yield from the Logiya Watershed, Lower Awash Basin, Ethiopia

by
Nura Boru Jilo
1,*,
Bogale Gebremariam
2,
Arus Edo Harka
1,
Gezahegn Weldu Woldemariam
3 and
Fiseha Behulu
4
1
Hydraulic and Water Resources Engineering Department, Haramaya Institute of Technology (HiT), Haramaya University, P.O. Box 138 Dire Dawa, Ethiopia
2
Department of Hydraulic and Water Resources Engineering, Arba Minch Water Technology Institute, Arba Minch University, P.O. Box 21 Arba Minch, Ethiopia
3
Geoinformation Science Program, School of Geography and Environmental Studies, Haramaya University, P.O. Box 138 Dire Dawa, Ethiopia
4
Institute of Technology (AAiT), Addis Ababa University, Addis Ababa 385, Ethiopia
*
Author to whom correspondence should be addressed.
Hydrology 2019, 6(3), 81; https://doi.org/10.3390/hydrology6030081
Submission received: 28 June 2019 / Revised: 16 August 2019 / Accepted: 20 August 2019 / Published: 11 September 2019
Browse Figures
Figure 1
<p>Location of the Logiya watershed, meteorological and flow gauging stations.</p> "> Figure 2
<p>Mean monthly precipitation and temperature (Tmax and Tmin) selected stations in the study area (1988–2016).</p> "> Figure 3
<p>Soil types (<b>a</b>) and Land Use/Land Cover (LULC) classes (<b>b</b>) in Logiya watershed.</p> "> Figure 4
<p>Sediment rating curve of the Logiya River at the Logiya gauge station.</p> "> Figure 5
<p>Historical trends of the annual precipitation (<b>a</b>), temperature maximum (Tmax) (<b>b</b>), temperature minimum (Tmin) (<b>c</b>), and sediment yield (<b>d</b>) from 1971–2000 in the Logiya watershed.</p> "> Figure 5 Cont.
<p>Historical trends of the annual precipitation (<b>a</b>), temperature maximum (Tmax) (<b>b</b>), temperature minimum (Tmin) (<b>c</b>), and sediment yield (<b>d</b>) from 1971–2000 in the Logiya watershed.</p> "> Figure 6
<p>Observed and downscaled average monthly maximum temperature (<b>a</b>), minimum temperature (<b>b</b>), and precipitation (<b>c</b>) of observed vs. bias-corrected and un-corrected RCP time series (1988–2016) of the Logiya watershed.</p> "> Figure 7
<p>Soil and Water Assessment Tool (SWAT) model calibration (<b>a</b>) and validation (<b>b</b>) for streamflow at the monthly time scale during 1992–1999 and 2000–2005, respectively.</p> "> Figure 8
<p>SWAT model calibration (<b>a</b>) and validation (<b>b</b>) for sediment yield at the monthly time scale during 1992–1999 and 2000–2005 respectively.</p> "> Figure 9
<p>Mean monthly observed and simulated sediment load (1992–2005).</p> "> Figure 10
<p>Spatial map of the observed sediment yield.</p> "> Figure 11
<p>Spatial map of the annual sediment yield (<b>a</b>) in 2030s RCP4.5 (<b>b</b>) in 2030s RCP8.5 (<b>c</b>) in 2060s RCP4.5 and (<b>d</b>) in 2060s RCP8.5 in the Logiya watershed.</p> "> Figure 12
<p>Projected change of mean monthly sediment yield (<b>a</b>) in the 2030s (<b>b</b>) in the 2060s from the baseline period in the Logiya watershed.</p> "> Figure 13
<p>Mean monthly sediment yield in the 2030s (<b>a</b>) and the 2060s (<b>b</b>) projected at the outlet of the Logiya watershed.</p> ">
Review Reports Versions Notes

Abstract

:
It is anticipated that climate change will impact sediment yield in watersheds. The purpose of this study was to investigate the impacts of climate change on sediment yield from the Logiya watershed in the lower Awash Basin, Ethiopia. Here, we used the coordinated regional climate downscaling experiment (CORDEX)-Africa data outputs of Hadley Global Environment Model 2-Earth System (HadGEM2-ES) under representative concentration pathway (RCP) scenarios (RCP4.5 and RCP8.5). Future scenarios of climate change were analyzed in two-time frames: 2020–2049 (2030s) and 2050–2079 (2060s). Both time frames were analyzed using both RCP scenarios from the baseline period (1971–2000). A Soil and Water Assessment Tool (SWAT) model was constructed to simulate the hydrological and the sedimentological responses to climate change. The model performance was calibrated and validated using the coefficient of determination (R2), Nash–Sutcliffe efficiency (NSE), and percent bias (PBIAS). The results of the calibration and the validation of the sediment yield R2, NSE, and PBIAS were 0.83, 0.79, and −23.4 and 0.85, 0.76, and −25.0, respectively. The results of downscaled precipitation, temperature, and estimated evapotranspiration increased in both emission scenarios. These climate variable increments were expected to result in intensifications in the mean annual sediment yield of 4.42% and 8.08% for RCP4.5 and 7.19% and 10.79% for RCP8.5 by the 2030s and the 2060s, respectively.

1. Introduction

Climate change is becoming a major environmental concern because increasing scientific evidence shows the high concentration of greenhouse gases (GHGs) in the atmosphere, and frequent hydro-meteorological extreme events are becoming the 21st century phenomenon [1,2]. Climate change results in significant impacts on life and natural resources [3]. The consequences of climate change involve adverse impacts on environment, hydrological cycle, water resources, agriculture and food security, human health, terrestrial ecosystems, and biodiversity [3,4,5,6,7,8]. These events are mostly manifested through subsequent influences of climate variables, including precipitation, temperature, and evapotranspiration [9,10,11]. Despite variations across regions, heat waves and extreme precipitation events will become more intense and more frequent. It is also estimated that the mean annual global air temperature induced by GHG emissions will likely increase by 1.4–2.0 °C by the end of the 21st century [10]. The rising GHG concentration due to continued gas emissions from different sources to the atmosphere affects climate variables and consequently alters hydrological cycles [10,12,13,14,15]. The change in climate—mainly due to precipitation and temperature patterns—could significantly influence soil erosion rates, streamflow, and sediment yield, which (directly or indirectly) adversely affects water resource availability and ecosystems [13,16,17,18,19]. In connection to this, numerous studies have indicated that changes in precipitation, temperature, and the interactions of these with land use and land cover change (LULCC) will be the main climate change related stress expected to exacerbate soil loss and sediment transport [20,21,22,23].
Soil erosion is a naturally occurring phenomenon through which the most productive topsoil materials are detached, transported, and deposited downstream by wind, water, and gravitational forces [24]. The displacement of soil layers due to water induced forces is causing continuous erosion, which in turn is changing the Earth’s surface [16,25]. This is increasing over recent years due to climate change and anthropogenic factors, resulting both in on-site and off-site effects [16,24]. Due to all these factors, soil erosion rates and the amount of suspended sediment discharges are being exacerbated across the world.
Over the past decade, the impacts of climate change on sediment yields have been widely studied around the globe [13,19,25,26,27,28]. Also, studies show that sediment transport is highly influenced by extreme precipitation and river discharge due to the non-linear relationship between water discharge and sediment transportation rates [29]. Moreover, Francipane et al. [25] reported that the magnitudes of mean and extreme events of sediment yield are expected to decrease with a high probability in the Walnut Gulch Experimental watershed in southeastern Arizona, USA. However, on their part, Zhang et al. [19] indicated that climate change is expected to increase the annual average sediment yield by 4–32% when compared to the base period for the Zhenjiangguan watershed, China. Azim et al. [12] examined the impact of climate change on the sediment yield of the Naran watershed, Pakistan using global climate model (GCM) predictors from the Special Report on Emission Scenarios (SRES) A2 and B2 and found mean annual sediment yield increases of 5–6% and 9–11% under the A2 and the B2 scenarios for time horizons of 2011–2040 and 2041–2070, respectively. Likewise, Shrestha et al. [12] used the Soil and Water Assessment Tool (SWAT) model in the Nam Ou Basin in northern Laos to predict the impact of climate change on sediment yield and reported changes in an annual sediment yield ranging from a 27% decrease to about a 160% increase. Further, Azari et al., [27] conducted a similar study focusing on the Gorganroud River Basin in northern Iran and found that the annual sediment yield would increase by 47.7%, 44.5%, and 35.9% for the 2040–2069 time series of the A1F1, the A2, and the B1 emission scenarios, respectively.
Farmers in the Horn of Africa, including Ethiopia, are frequently vulnerable to climate change extremes such as floods and droughts, both of which are the most common hydro-meteorological extremes in tropical countries [5,6,7,8,30,31,32]. Ethiopia is a country in which the economy largely depends on rain-fed agriculture, which is highly vulnerable to the impacts of climate change, mainly changes in precipitation patterns [33,34]. Precipitation pattern and character in Ethiopia are predominantly controlled by the oscillation of the Intertropical Convergence Zone. This precipitation has a strong impact on soil erosion. A change in precipitation intensity and seasonal distribution results in a significant impact on the soil erosion rate [20]. According to Zhang and Nearing [35], precipitation and soil loss are directly proportional to each other.
The Logiya watershed is one of the sub-watersheds in the Lower Awash River Basin (LARB). It is located in the arid lowlands of the Afar Region in the northeastern part of Ethiopia. The upstream of the Logiya Watershed extends up to North Wollo Zone of the Amhara Regional State, and the outlet of the watershed is situated in the Afar Regional State. The precipitation and the temperature of the watershed both fluctuate, which subsequently results in a variation of the river flow in the watershed under undue climate change conditions. The Logiya watershed is characterized by severe land degradation and desertification, but high flood also occurs, mostly in the upstream areas during high rain seasons [36].
Several studies have been conducted to assess flooding and droughts caused by climate change in the LARB, including the Logiya watershed [36,37,38,39,40,41,42]. However, none of these studies addressed the climate change impacts on sediment yield both in the LARB and the Logiya watershed. The magnitude of climate change effects on sediment yield varies depending on the region in focus and the climate scenario taken into consideration [43]. Usually, heterogeneity of the climate conditions is wrongly assumed in the previous scenarios by some researchers. However, it is necessary to regionalize the assessment of climate change impacts. In other words, understanding the peculiar effects of climate change on sediment load is critical for watershed management [44].
In order to examine the non-linear interplays between climate-induced processes that influence sediment discharge at the watershed level, hydrological and sediment models can be used in combination with climate projections from global circulation models and regional climate models [25,45]. The present study used representative concentration pathways (RCPs), which supersede the Special Report on Emission Scenario (SRES), a newly developed scenario that provides input for climate models and a distributed hydrological model, SWAT, to estimate potential impacts of climate change on sediment yield in the Lower Awash Logiya watershed. The specific objectives of the study were to (i) assess future precipitation, air temperature, and evapotranspiration; (ii) estimate the spatial variation of sediment yield under different climate scenarios; and (iii) evaluate the impacts of climate change on sediment yield from the Logiya watershed under RCPs scenarios.

2. Materials and Methods

2.1. Description of the Study Area

The Awash River Basin (ARB) is one of the major river basins in Ethiopia. It originates in the highlands of Central Ethiopia and flows northeastwards, where it eventually drains into the Lake Abe. The total length of the river is about 1200 km and has a total drainage area of about 110,000 km2. The ARB is divided into three distinct zones: Upper Basin, Middle Basin, and Lower Basin. The Logiya watershed is among the sub-watersheds of LARB and is situated in the western part of the basin between 11°28′21″ N and 12°04′55″ N and 39°40′30″ E and 40°56′15″ E (Figure 1). The watershed drainage area is about 3151.86 km2 up to the outlet. The elevation of the study area ranges from 3426 to 379 m a.m.s.l. with an average altitude of 890.6 m. Awash River Basin ranges from cold high mountainous zones to semi-desert lowlands with extreme ranges of temperature and precipitation [40]. There are three seasons in the Awash River Basin based on the movement of the Intertropical Convergence Zone (ITCZ) and the amount of precipitation timing. The three seasons are named: Kiremt, which is the main rainy season (June–September), Bega, which is the dry season (October–January), and Belg, the small rainy season (February–May). Therefore, the precipitation was categorized under bi-modal pattern. The mean annual precipitation varies from 1600 mm in the elevated areas to 160 mm in the LARB. In the same way, the mean annual temperature of Awash River Basin ranges from 20.8 °C in the upper part to 29 °C in the lower part. Figure 2 shows observed mean monthly precipitation as well as maximum and minimum temperatures of the Logiya watershed.

2.2. Data Sources

The required input data used to run the SWAT models necessary for this study were drawn from hydrological and meteorological data, soil maps, Land Use/Land Cover (LULC) maps, and digital elevation models (DEM) of the watershed. These were collected from different sources, to be discussed below.

2.2.1. Digital Elevation Model

The DEM is used to describe the elevation of any point in a given area at a specific spatial resolution. The 30 m × 30 m resolution DEM data covering the watershed area was downloaded from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model (GDEM) website [46]. It was required for the SWAT model to delineate the watershed, the sub-watershed, the topography, and the drainage pattern of the study area. The DEM data was loaded into Arc SWAT interface after being re-projected into the World Geodetic System (WGS 84) spheroid, a spatial reference system of the Universal Transverse Mercator (UTM) with Datum Zone 37 N.

2.2.2. Hydrological and Meteorological Data

This includes the daily observed climate data from 1988–2016 comprising data about the precipitation, maximum and minimum temperature, solar radiation, wind speed, and relative humidity. These were collected from the National Meteorological Agency (NMA) of Ethiopia [47]. For this purpose, four stations, namely, Dubti, Mersa, Kobbo, and Bati (Table 1), were selected based on the quality and the availability of long year recorded data. Moreover, hydrological data including streamflow and suspended sediment yield at Logiya gauging station was collected from the Ministry of Water, Irrigation, and Electricity (MoWIE). Daily streamflow data at the outlet of the watershed were used for the calibration and the validation of the SWAT model. The missing hydrological data were filled by linear regression from the nearby Mille river gauged station. Moreover, the meteorological data of each neighboring station were filled by the Normal-Ratio Method, while the homogeneity and the consistency of data were tested and used as inputs for the climate and the hydrological models [48,49,50].

2.2.3. Soil Data

The soil data of the Logiya watershed was obtained from the MoWIE. This data was prepared according to the Food and Agriculture Organization (FAO) classification scheme. Accordingly, there were ten soil types found in the study area (Figure 3a). These are included Calcaric flubisols, Chromic luvisols, Dystric nitisols, Eutric cambisols, Eutric, regosols, Haplic xerosols, Leptosols, Orthic solonchaks, Vertic cambisols, and Vitric andosols, with a proportion of each class contributing 29.04%, 0.15%, 0.90%, 2.72%, 19.16%, 6.57%, 12.22%, 18.76%, 3.27%, and 7.21%, respectively, to the total watershed area.

2.2.4. Land Use/Land Cover

The most dominant LULC classes in the Logiya watershed included grass steppe, shrub, tree steppe and bare soil and very sparse vegetation (Figure 3b). It has been found that LULC has impacts on climate change, which in turn influences surface runoff, soil erosion, sediment load, and evapotranspiration in a watershed as the SWAT model simulation shows [51].

2.2.5. Climate Models and RCPs Scenarios

The coordinated regional climate downscaling experiment (CORDEX) is a program sponsored by the World Climate Research Program (WCRP) to generate historical and future climate projections at the regional scale. A regional climate model (RCM) was used to advance the predictive model of the Earth’s climate in the scientific analysis of the dominant sets of the governing processes. It describes climate change on a regional scale [52]. Thus, the RCM was integrated into the CORDEX-Africa domain with a horizontal grid resolution of 0.44° × 0.44° (50 km × 50 km) [53].
The historical simulations were forced by observed data and anthropogenic atmospheric composition covering the period 1951–2005, whereas the projected scenarios between 2006 and 2100 were derived through RCP4.5 and RCP8.5 scenarios. Also, simulated data of RCP4.5 and RCP8.5 scenarios were generated from Hadley Global Environment Model 2-Earth System (HadGEM2-ES) of the global climate model outputs. These were dynamically downscaled from the CORDEX-Africa database.
Furthermore, future scenario analyses RCP4.5 (medium emission scenarios) and RCP8.5 (high emission scenarios) were used to assess the precipitation, maximum and minimum temperature, evapotranspiration, and sediment yield. The analyses were performed in two-time horizons in a time period of 30 years for both baseline and future periods. The baseline period used for this study was 1971–2000, while future scenario analysis involved 2020–2049 (2030s) and 2050–2079 (2060s).

2.2.6. Accuracy of Precipitation Simulations from the Climate Model

The precipitation simulations from the HadGEM2-ES global climate model outputs were evaluated using observed data. Accuracy of the precipitation was evaluated using statistical methods, including bias, correlation coefficient (CC), and coefficient of variation (CV), computed as follows [54]:
B i a s = R ¯ r c m − R ¯ o b s R ¯ o b s × 100
C V = σ R R ¯ × 100
C o r r e l = ∑ i = 1 N ( R o b s − R ¯ o b s ) ( R r c m − R ¯ r c m ) ∑ i = 1 N ( R o b s − R ¯ o b s ) 2 ( R r c m − R ¯ r c m ) 2
where R ¯ is an average precipitation over the watershed; rcm and obs are subscripts representing the precipitation amount over the watershed either from RCM simulation or observed datasets, respectively; σ indicates the standard deviation of either the RCM or the observed precipitation data, and R represents estimated statistics individually either for RCM or observed precipitation amount.

2.2.7. Bias Correction Method

Using downscaled regional climate data for impact assessment without any bias correction may lead to considerable deviation when a hydrologic model is forced with a biased RCM [10,55,56,57,58]. Hence, a bias correction method was applied for each daily precipitation and temperature dataset derived from CORDEX-Africa to minimize the systematic statistical deviation of climate input data from the observational one before using them for any modeling and future climate change projection. The bias correction method used for this study was non-linear correction to adjust observed and simulated precipitation data. Correction factors were computed from the statistics of the observed and the simulated variables. The principle of this method shows the mean and the standard deviation of the daily precipitation data becoming equal to those of the observed data [59,60,61].
P * = a P b
where P * is the corrected value of precipitation, and a and b are power transformation parameters determined for the period of a year, which were obtained from calibration in the baseline period and subsequently applied to the projection period.
Hence, the bias correction of temperature only involved shifting and scaling to adjust the mean and the variance [62,63]. The corrected daily temperature T c o r r was obtained as:
T c o r r = T ¯ o b s +   σ ( T o b s ) σ ( T r c m ) ( T r c m   −   T ¯ o b s ) + ( T ¯ o b s − T ¯ r c m )
where T c o r r is the corrected daily temperature; T ¯ r c m is the uncorrected daily temperature from the RCM model; T ¯ o b s is observed daily temperature; T ¯ o b s and T ¯ r c m are mean temperatures for observed and simulated datasets, respectively.

2.2.8. Statistical Test for Trend and Variability Analysis

For this study, non-parametric statistical method, the Mann-Kendall (MK) trend test, was used to identify the presence of trends in the time series at the 5% significance level [64,65,66,67]. The MK trend test analysis was recommended by the World Meteorological Organization (WMO) and has been widely used in practice to evaluate the significance of monotonic trends in hydrological and meteorological time series [14,68,69] and to determine the changes (increasing, decreasing, or trendless) in the values of climatic variables.
For a given time series X (x1, x2, …, xn), the statistic S is defined as:
S =   ∑ i   =   1 n ∑ j   =   i   +   1 n − 1 s g n ( x j − x i )
where
S i g n ( X j − X i ) { 1 , x j > x i , 0 , x j = x i , − 1 , x i < x j ,
x i and x j are the data values in the time series, and n is the number of data points in the time series,
The variance S is calculated using the following equation when the value of n ≥ 8;
V a r ( s ) =   n ( n − 1 ) ( 2 n + s ) 18
The standard normal test statistic is calculated using the following equation;
Z =   { ( 1 − S ) V a r ( S ) ,         S > 0 , 0 ,                               S = 0 , ( 1 + S ) V a r ( S ) ,         S < 0 .
The MK statistic is denoted by S follows the standard normal distribution. Positive and negative values of S denote an upward and downward trends, respectively. The null hypothesis of no trend is rejected if |S| > 1.96 at the 5% significance level.

2.3. Sediment Rating Curve Development

Sediment is some of the most important data in the hydrological study for the estimate of sediment transport. Due to the scarcity of continuous sediment data, the sediment rating curve was necessary to develop the relationship between daily streamflow and sediment data measured at the outlet of the Logiya watershed. It displayed the rate of sediment transport as a function of streamflow [26]. The suspended load in milligram per liter was converted to tons/day as follows:
Q s = 0.0864 × Q × C
where Q s is sediment load in (t day−1), Q is streamflow (m3 s−1) and C is suspended sediment concentration (mg L−1). As depicted in Figure 4, the sediment rating curve was developed by using linear regression techniques, and the curve is expressed in the form of a power-law type equation [29,70,71,72].

2.4. SWAT Model Setup

2.4.1. SWAT Model Description

SWAT is one of the most powerful hydrologic models [73]. It is physically semi-distributed, conceptually and computationally efficient and operates on a daily time step according to the soil properties, topography, vegetation and land management practices in the watershed [51]. The model was developed by the United States Department of Agricultural Research Service [74] to predict and qualify the impact of agricultural management practices on water, sediment and chemical yields in different size of catchments with varying soils, land use and management conditions over long periods of time [51,75].
The model has been widely applied to the different sizes of the watersheds and its applications all over the water resources, including the effects of climate and soil erosion [51]. In our study, the SWAT model was used to estimate hydrologic components such as surface runoff, evapotranspiration, groundwater flow, sediment yield and for each hydrologic response unit (HRU). In the hydrological component, surface runoff was estimated separately for each sub-basin of the total watershed area and routed to obtain the total surface runoff for the watershed. Surface runoff volume was estimated from daily rainfall using modified soil conservation service curve number (SCS-CN) and Green and Ampt methods [76]. Potential evapotranspiration can be estimated by three methods: the Hargreaves [77], the Priestley–Taylor [78], and the Penman–Monteith [79]. Based on the availability of observed data, the Hargreaves method was adopted for this study. In each HRU sediment yield was estimated using a modified universal soil loss equation (MUSLE). The modified universal soil loss equation was determined as [74]:
S e d = 11.8 × ( Q s u r f × q p e a k × A r e a H R U ) 2 × K U S L E × C U S L E × P U S L E × L S U S L E × C F R G
where S e d is the sediment yield (metric ton day−1), Q s u r f is the surface runoff volume (mm ha−1), q p e a k is the peak runoff rate in m3 s−1, A r e a H R U is the area of HRU in ha, KUSLE is the soil erodibility factor, CUSLE is the cover and management factor, PUSLE is the support practice factor, LS is topographic factor, and CFRG is the course fragment factor.
The routing phase of SWAT defines the movement of water, nutrients, sediment, and pesticides through the channel network of the watershed into the outlet. The sediment-routing model that simulates sediment transport in the channel network consists of two components operating simultaneously—deposition and degradation [74]. The amounts of deposition and degradation are based on the maximum concentration of sediment in the reach and the concentration of sediment in the reach at the beginning of the time step.
The final amount of sediment in the reach was determined by:
S e d c h = S e d c h , i − S e d d e p + S e d d e g
where S e d c h is the amount of suspended sediment in the reach (metric tons days−1), S e d c h , i is the amount of suspended sediment in the reach at the beginning of the time period (metric tons days−1) S e d d e p is the amount of sediment deposited in the reach segment (metric tons days−1) and S e d d e g is the amount of sediment re-entrained in the reach segment (metric tons days−1). The amount of sediment transported out of the reach was calculated as:
S e d o u t = S e d c h × V o u t V c h
where S e d o u t is the amount of sediment transported out of the reach (metric tons days−1), V o u t is the volume of outflow during the time step (m3), and V c h is the volume of water in the reach segment (m3).

2.4.2. Watershed Delineation

The watershed delineator tool in ArcSWAT 2012 interface allowed us to delineate the watershed and the sub-watersheds using the 30 × 30 m DEM. Flow direction and accumulation are the concepts behind the definition of the stream network of the DEM in SWAT. The monitoring point was added manually, and the numbers of the sub-basin were adjusted accordingly. There were 33 sub-catchments with an area ranging between 239.71 and 47,022.06 hectares.

2.4.3. Hydrologic Response Units (HRU) Analysis

HRU analysis helps to load land use and soil maps as well as to incorporate classification of HRU into different slope classes. HRUs were assigned to each sub-catchment based on a 10% threshold value for LULC, soil, and slope categories, as suggested by the user’s manual of SWAT to increase the chance of inclusion.
Based on the topographical characteristics of the terrain in the DEM, the SWAT model described watershed spatial variability by further splitting the sub-watershed into homogeneous characteristics, lumped land areas, and HRUs based on topography, soil, land use, and slope. Multiple land use/soil/slope methods were used to define the HRU with land use (10%), soil (10%), and slope (15%) thresholds. The number of HRUs was defined by eliminating the percentages of land use, soil, and slope values that covered a percentage of the sub-catchment area less than the threshold level. Accordingly, the Logiya watershed was divided into a total of 234 HRUs of different LULC, soil, and slope combinations.

2.5. SWAT Calibration and Uncertainty Procedures

SWAT Calibration and Uncertainty Procedures (SWAT-CUP) is an interface that was developed for the SWAT model [80]. It is a computer program for integrated various sensitivity analysis, calibration, and validation programs for the SWAT model. Sensitivity analysis is used to check the rate of change in output parameters with respect to the input parameters of the model. For this study, the measured stream flow and sediment yield data obtained from the discharge–sediment rating curve at the outlet of the Logiya watershed were used to calibrate and validate the model by SWAT-CUP using Sequential Uncertainty Fitting version-2 (SUFI-2) algorithm from 1990 to 2005. The first two years were considered as the model warm-up period, the period from 1992 to 1999 was used for calibration, while the remaining data (2000–2005) were used for model validation.
The model performance efficiency was determined by comparing observed against simulated hydrographs. The coefficient of determination (R2), the Nash–Sutcliffe simulation efficiency (NSE), and the percent bias (PBIAS) were used to evaluate the model‘s performance during calibration and validation processes to test the goodness of fit between monthly simulated and observed values. These parameters were important to investigate and determine the accuracy of the SWAT model simulations. The determination coefficient (R2) describes the proportion of variance in measured data from the model. The value of R2 ranged from zero (which indicated the model was poor) to one (which indicated the model was perfect), with higher values indicating less error variance; typical values greater than 0.6 were considered acceptable [53]. The R2 was calculated using the following Equation:
R 2 = [ ∑ i = 0 n ( Q O b s − Q ¯ O b s ) × ( Q S i m − Q ¯ S i m ) ] 2 ∑ i = 0 n ( Q O b s − Q ¯ O b s ) 2 ∑ i = 0 n ( Q S i m − Q ¯ S i m ) 2
If, R2 = 0, none of the variances in the measured data were replicated by the model predictions and if R2 = 1, it indicated that all of the variances in the measured data were replicated by the model predictions. NSE indicated how well the plots of the observed versus the simulated data fit the 1:1 line. NSE was computed as follows (Equation (15)):
N S E = 1 − ∑ i = 0 n ( Q O b s − Q S i m ) 2 ∑ i = 0 n ( Q O b s − Q ¯ O b s ) 2
The value of NSE ranged from negative infinity to one (perfect) where values lower than zero indicated that the mean observed value was a better predictor than the simulated value and revealed unacceptable model performances. NSE values greater than 0.5 showed the simulated value was a better predictor than the mean measured value and was generally viewed as acceptable performance [81].
PBIAS was used to measure the average tendency of the simulated data to be larger or smaller than the observed values. PBIAS is expressed in percentage; a lower absolute value of the PBIAS is better and determines model performance.
P B I A S = ∑ i = 0 n ( Q O b s − Q S i m ) ∑ i = 0 n ( Q O b s ) × 100
where n is the number of observations during the simulation period, QObs is the observed flow data, QSim is the simulated flow value with the respective time, and QObs and QSim are the arithmetic means of the observed and the simulated values, respectively
The optimal value of PBIAS was zero, with low-magnitude values indicating accurate model simulation. Positive values indicated model underestimation bias, and negative values indicated model overestimation bias [82].

3. Results and Discussion

3.1. Statistical Trend Analysis of Historical Climate Variability

The variation trends of annual precipitation, maximum and minimum temperature, and sediment yield in the Logiya watershed from 1971–2000 can be found in Figure 5a–d. In the watershed, the annual precipitation and the sediment yield showed insignificant decreasing trends, and the annual maximum and minimum temperatures showed insignificant increasing trends, which meant the p-value was greater than the significance level α (5%) (Table 2). Therefore, results of mean annual precipitation as well as maximum and minimum temperatures were related to other findings reported for the country. For example, Atlas [83] reported decreasing precipitation and increasing temperature trends in the Afar region during the period 1981–2014. Likewise, Jury and Funk [84] showed annual air temperature increased and precipitation decreased over the time series 1948–2006 in Ethiopia.
The mean, the maximum (Max), the minimum (Min), the standard deviation (SD), and the coefficient of variance (CV) of annual precipitation, maximum and minimum temperature, and sediment yield of the watershed are depicted in Table 2. The mean annual precipitation of the watershed during the baseline period was 125.4 mm with 54.3 mm standard deviation and 43.3% CV. The minimum and the maximum ever-recorded precipitations were 58.0 mm and 304.5 mm, respectively. CV was used to classify the degree of variability in the historical data events as less (CV < 20), moderate (20 < CV < 30), high (CV > 30), or very high (CV > 40%), but CV > 70% indicated extremely high inter-annual variability [69]. The coefficients of variation of annual precipitation and sediment yield were above 40% and 70%. This indicated very high variability of precipitation and sediment yield over the watershed.

3.2. Evaluating Accuracy of Climate Model Simulations

Impact assessment of climate change necessitates evaluation of the climate model’s accuracy. For this reason, in this study, the model’s performance was evaluated using statistical measures such as bias, CV, and CC [54,85]. Observed mean annual precipitation as well as bias-corrected and bias-uncorrected historical RCPs framed for the Logiya watershed were 742.46, 733.65, and 408.86 mm, respectively (Table 3). Despite a slight underestimation in model simulation, the mean annual precipitation of the observed and the model output showed a good relationship. The results revealed the degree of rainfall coefficient variability (CV = 2.60%) for observed and bias-corrected RCP data (CV = 3.12%) and for bias-corrected RCP data, while there was slight underestimation with bias = −1.19%, which suggested that, in the watershed, the precipitation is well captured or represented, and a good linear relationship was established, as evaluated by correlation (CC = 0.56) for the bias-corrected model simulation data. The average result downscaled dynamically with maximum and minimum temperatures, which also virtually indicated a good relationship with the observed data. Therefore, comparisons of observed and corrected RCP data were consistent (Figure 6).

3.3. Future Climate Projection

This study also analyzed future climate patterns such as precipitation, temperature, and evapotranspiration under RCP4.5 and RCP8.5 scenarios. The evaluation was made over two consecutive 30-years period of the 2030s and the 2060s, respectively. Accordingly, the future climate change was estimated for each RCP based on the baseline period (1971–2000). The mean annual precipitations were projected under RCP4.5 to increase by 0.9% and 4.11% for the 2030s and the 2060s, respectively. Likewise, for RCP8.5, the mean annual precipitation increases were found to be 3.55% in the 2030s and 8.69% in the 2060s. In both RCP4.5 and RCP8.5 scenarios, the results showed an increase in precipitation patterns in July–January, whereas they showed decreased precipitation patterns in February–June in the 2030s and the 2060s, respectively. Besides, the results pointed out that the alteration of precipitation patterns depended on the seasonality of the watershed. Supporting the present study’s findings, Intergovernmental Panel on Climate Change (IPCC) [10] reported an increase of precipitation in rainy seasons with heavy precipitation events over the world. The global precipitation pattern is expected to increase up to 20% [10]. Getahun and Lanen [86] reported projected change in main rainy season precipitation with less than ±20% in all GCMs in the Upper Awash River Basin. The predicted increase of precipitation was unevenly distributed in a year, which would further increase the seasonal variation of precipitation in the watershed [19].
In this study, the predicted temperatures in all months might increase in both RCP scenarios. Under RCP4.5, changes in monthly maximum temperatures ranged from 0.1 °C to 0.47 °C with an average annual of 0.31 °C by the 2030s and from 0.07 °C to 0.26 °C with an average annual of 0.17 °C by the 2060s. For the RCP8.5 scenario, the monthly maximum temperature increased from 0.03 °C to 0.73 °C and from 0.09 °C to 0.27 °C with average annuals of 0.33 °C and 0.19 °C for the 2030s and the 2060s, respectively. The results also indicated that all the values of monthly minimum temperature increased from 0.03 °C to 0.45 °C by the 2030s and from 0.03 °C to 0.49 °C by the 2060s under RCP4.5 with average annuals of 0.25 °C and 0.32 °C, respectively, compared to the baseline period. Additionally, in the RCP8.5 scenario, the minimum temperature increased from 0.05 °C to 0.56 °C for the 2030s and from 0.02 °C to 0.57 °C for the 2060s. The average annual minimum temperatures were 0.28 °C and 0.35 °C under RCP8.5 for the 2030s and the 2060s, respectively.
These results agreed with the findings reported by Daba et al. [87] for the Upper Awash Sub-Basin showing an increase in predicted maximum and minimum temperatures, but the changes in temperature were uneven throughout the month. Getahun and Lanen [86] and Kinfe [88] projected maximum and minimum temperature increases in old SRES and in both RCP4.5 and RCP8.5 scenarios in the Upper Awash River Basin. Overall, these findings were in line with previous studies conducted in the Awash Basin [86,87,88,89,90]. The increased temperature was additionally related with an increase in evapotranspiration in all RCP scenarios. Finally, the projected increment in temperature change under RCP4.5 was lower than under RCP8.5 for both maximum and minimum temperatures because RCP8.5 had higher emission of greenhouse gases than the RCP4.5 scenario [10].

3.4. SWAT Model Calibration and Validation

3.4.1. Sensitivity Parameters of Streamflow Analysis

Twelve of the most sensitive parameters were selected for streamflow simulation of the Logiya watershed with their fitted value ranked in Table 4.

3.4.2. Streamflow Calibration and Validation

Streamflow calibration and validation was done by comparing the observed and the simulated flow values for a 16-years period (1990–2005). The first two years (1990–1991) were considered as a warm-up period, the 1992–1999 period was used for the model calibration, and the 2000–2005 period was used for the validation period. Both calibration and validation for streamflow simulation obtained satisfactory results of fit with R2 of 0.8 and 0.77 and NSE coefficients of 0.73 and 0.63, respectively (Table 5 and Figure 7). The results revealed that the variation pattern of simulated streamflow was generally consistent with that of the observed streamflow.

3.4.3. Sediment Yield Calibration and Validation

The most sensitive parameters that affect sediment yields used in the Logiya watershed are the linear parameter for calculating the maximum amount of sediment (SPCON), the USLE equation support particle factor (USLE_P), and the exponent parameter for calculating sediment re-entrained channel sediment routing (SPEXP) (Table 6).
The SWAT-CUP in parallel gives the best simulation of sediment yield using the most sensitive parameter values. Using fitted parameters of calibration, the simulation was undertaken for the validation period, and its performance was checked with the observed data. The best simulation performances of sediment yield results using these fitted parameters and the values for monthly calibration and validation are shown in Table 7 and Figure 8.
Numerous studies (e.g., Setegn et al. [8], Santhi et al. [81], Moriasi et al. [82], Benaman et al. [91]) suggested that the prediction efficiency of the calibrated model can be a good agreement if R2 and NSE values are greater than 0.6 and when the value of PBIAS is between ±15 ≤ and ≤±30. Since, in the present study, the values of R2 and NSE were 0.83 and 0.79, respectively, the calibration and the validation of the observed and the simulated streamflow and sediment yield of the Logiya watershed showed good relations with reasonable accuracy [8,81,82,90].
The results indicated that, for the Logiya watershed, the SWAT model simulated streamflow and sediment yield with reasonable accuracy. Mean annual sediment loads of observed and simulated values in the years 1992–2005 were 19.53 tons ha−1 and 23.66 tons ha−1, respectively. This indicated the estimates of sediment load in the model were greater than the observed values. This may have been due to limited sediment data, quality and constraint of weather data, streamflow, errors during data recording, and uneven distribution of gauging stations in the watershed [75,92,93]. Underestimates and overestimates of mean monthly data have been recorded in the model outputs over the observed high flow months such as February, March, April, July, August, and September (Figure 9).

3.5. Spatial Variation of Sediment Yield under Different Climate Scenarios

It was found that the SWAT model is advantageous in providing quantitative estimates of soil loss and sediment yield and their spatial variations in different sub-watershed by considering HRU [94,95,96]. This result plays a critical role for the detailed understanding and evaluation of the amount of sediment yield as well as for designing an appropriate management scenario at each sub-watershed. The spatially distributed sediment yield in each sub-watershed under different climate scenarios is demonstrated in Figure 10 and Figure 11. During the observed periods, the rate of sediment yield ranged from 3.28 t ha−1 yr−1 to 51.77 t ha−1 yr−1, which corresponded to the sub-watersheds 30 and 27 of the Logiya watershed, respectively, with a mean value of 19.07 t ha−1 yr−1.
Besides, the projected sediment yield estimated under both RCP4.5 and RCP8.5 scenarios showed the spatial variation across the watershed. Severely eroded areas were situated at downstream parts of the study area. These findings agreed with studies conducted on the Lower Awash basin and the Upper Blue Nile River basin [37,94]. For RCP4.5, by the 2030s and the 2060s, the mean values of sediment yields were 2.57 t ha−1 yr−1 and 2.80 t ha−1 yr−1, respectively. Also, under the RCP8.5 scenario, the mean values of sediment yields were 3.31 t ha−1 yr−1 and 4.07 t ha−1 yr−1 for the 2030s and the 2060s, respectively. Results indicated that for RCP4.5 and RCP8.5, the maximum and the minimum sediment loads of the watershed occurred in the sub-watersheds 27 and 19, respectively, in both time series data. In the case of the Logiya watershed, more than two-thirds (76% of the total study area) of the observed soil loss rate was above the tolerable erosion limit suggested by Hurni [97] for various agro-climatic zones of Ethiopia, which ranges from 2 to 18 t ha−1 yr−1.
In addition, the results indicate that more severe erosions occurred in the downstream parts of the watershed, than in the upstream areas of the watershed. This shows that soil erosion was a serious phenomenon in the different sub-basins in this study area during observed and projected time series. During the observed and the projected future RCP scenarios, the maximum sediment load occurred in sub-basin 27 of the watershed. Also, bare land and Eutric Regosols were the most dominant LULC and soil type of sub-basin 27, respectively. All this indicate that LULC was more sensitive to soil erosion in this sub-basin. What’s more, it can be understood that Eutric Regosols soil type mostly occurs in desert regions, and it occupies the main LULC of that area. Kefyalew [37] indicated that the Lower Awash basin was highly affected by flood and soil erosion year-to-year.

3.6. Impact of Climate Change on Sediment Yield and Streamflow in Logiya Watershed

Similar to future precipitation, sediment yield estimates showed an increasing pattern in July–January and a decreasing pattern in February–June in both RCP scenarios (Figure 12). The mean annual sediment yields were expected to increase by 4.42% and 8.08% in the 2030s and the 2060s, respectively, for RCP4.5. They were also expected to increase by 7.19% for the 2030s and the 10.79% for 2060s under the RCP8.5 scenario. The results indicated that precipitation and sediment yield had strong relationships, and this agreed with previous findings [12,23,98,99]. Moreover, the relationship between streamflow and sediment load at the outlet of the Logiya watershed was strongly related. Both streamflow and sediment yield were predicted to increase under both RCP scenarios. Mean annual streamflows were predicted to increase by 1.43% and 3.47% for RCP4.5 in the 2030s and the 2060s, respectively. For the RCP8.5 scenario, the model showed an increase of 2.81% by the 2030s and 5.48% by the 2060s. This implied that sediment yield at the outlet of the watershed was influenced by streamflow. Overall, both RCP scenarios showed fairly similar patterns of increment in projected monthly precipitation, streamflow, and sediment yield in a fair amount in the study area. Similar results were reported by Shrestha et al. [13], Adem et al. [28], and Azari et al. [5], all of whom showed increasing streamflows proportional to increasing sediment yield. Figure 12 displays the projected mean monthly sediment yield at the outlet of the Logiya watershed under all RCP scenarios.
Figure 13 depicts the mean annual sediment load for RCP4.5 which was 23.51 tons ha−1 and 22.82 tons ha−1 for the 2030s and the 2060s, respectively. Also, for the RCP8.5 scenario, they were 24.3 tons ha−1 by 2030s and 23.1 tons ha−1 by the 2060s. The findings indicated that the monthly maximum and minimum sediment loads would happen in August−June. This pattern showed that the observed and projected sediment load in this study area were strongly related.

4. Conclusions

The present study investigated the impacts of climate change on sediment yield from the Logiya watershed in the Lower Awash Basin, Ethiopia. RCMs in CORDEX-Africa derived from the HadGEM2-ES Global Climate Model were used for RCP4.5 and RCP8.5 by the 2030s and by the 2060s as compared to the baseline period (1971–2000). The SWAT hydrological model was used to simulate present and future changes in sediment yield, and the SUFI-2 algorithm in the SWAT-CUP program was used for parameter adjustment. Calibration, validation, and uncertainty analyses for sediment suggested that the SWAT model could be applied to simulate future changes in sediment yields due to ultimate climate change.
Results indicated that mean annual precipitation as well as maximum and minimum temperatures would increase under RCP4.5 and RCP8.5 scenarios. The projected annual precipitations of RCP4.5 in the 2030s and the 2060s would increase by 0.9% and 4.11%, respectively, and in the RCP8.5 scenario would increase from 3.55% and 8.69% for the 2030s and the 2060s, respectively. The RCP8.5 scenario was changed by a higher magnitude than the RCP4.5 scenario because RCP8.5 had a higher greenhouse gas emission scenario with a higher degree of global warming than the RCP4.5 scenario.
Results of climate change impacts on sediment yield clearly showed that sediment loss would increase in all RCP scenarios. The mean annual sediment yields were predicted to increase by 4.42% by the 2030s and by 8.08% by the 2060s under RCP4.5. For RCP8.5, they were found to be 7.19% and 10.79% by the 2030s and the 2060s, respectively. Moreover, the projected sediment yield was found to be a similar pattern of precipitation and streamflow. This indicated that precipitation, streamflow, and sediment load would operate in good relation with this watershed. Finally, the results showed that the amount of sediment load would vary in each sub-watershed due to changes in precipitation and temperature in the 2030s and the 2060s for both RCP4.5 and RCP8.5 scenarios. The maximum and the minimum sediment yields would occur in sub-basins 27 and 19 in the watershed.

Author Contributions

All authors made substantial contribution to the development of this manuscript. N.B.J. was in charge of conceptualization, data analysis and writing of original draft. B.G., A.E.H., G.W.W. and F.B. reviewed, edited and improved the manuscript. All authors read and approved the final manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Seneviratne, S.; Nicholls, N.; Easterling, D.; Goodess, C.; Kanae, S.; Kossin, J.; Marengo, J.; McInnes, K.; Rahimi, M.; Reichstein, M.; et al. Changes in Climate Extremes and Their Impacts on the Natural Physical Environment, Chapter A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change (IPCC); Cambridge University Press: Cambridge, UK; New York, NY, USA, 2012; pp. 109–230. [Google Scholar]
  2. Eldiiar, D.; Xi, C.; Amobichukwu, C.; Amanambu, F.; Rustam, O.; Gulnura, I.; Gulkaiyr, O. Projected Rainfall Erosivity Over Central Asia Based on CMIP5 Climate Models. Water 2019, 11, 897. [Google Scholar] [CrossRef]
  3. Surya, G. Simulating Climate Change Impact on Soil Erosion and Soil Carbon Sequestration; Indian Space Research Organization: Bangalore, India, 2015.
  4. World Bank. Costing Adaptation through Local Institutions Village Survey Results: Ethiopia; World Bank: Washington, DC, USA, 2011. [Google Scholar]
  5. Di Falco, S.; Veronesi, M.; Yesuf, M. Does adaptation to climate change provide food security? A micro-perspective from Ethiopia. Am. J. Agric. Econ. 2011, 93, 829–846. [Google Scholar] [CrossRef]
  6. Wanders, N.; Wada, Y.; Van Lanen, H. Global hydrological droughts in the 21st century under. Earth Syst. Dyn. 2015, 6, 1–15. [Google Scholar] [CrossRef]
  7. Houghton, D. Introduction to Climate Change: Lecture Notes for Meteorologists; Secretariat of the World Meteorological Organization: Geneva, Switzerland, 2002. [Google Scholar]
  8. Setegn, S.; Rayner, D.; Melesse, A.; Dargahi, B.; Srinivasan, R. Impact of climate change on the hydro climatology of Lake Tana basin, Ethiopia. Water Resour. 2011, 47, 1–13. [Google Scholar]
  9. Intergovernmental Panel on Climate Change (IPCC). Climate Change: AR4 Synthesis Report; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  10. Intergovernmental Panel on Climate Change (IPCC). Climate Change: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stoker, T., Qin, D., lattner, K., Tigner, M., Allen, K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013. [Google Scholar]
  11. Pumo, D.; Arnone, E.; Francipane, A.; Caracciolo, D.; Noto, L.V. Potential implications of climate change and urbanization on watershed hydrology. J. Hydrol. 2017, 554, 80–99. [Google Scholar] [CrossRef]
  12. Azim, F.; Shakir, A.; Rehman, H.; Kanwal, A. Impact of climate change on sediment yield for Naran watershed, Pakistan. Int. J. Sediment Res. 2016, 31, 212–219. [Google Scholar] [CrossRef]
  13. Shrestha, B.; Babel, M.; Maskey, S.; Van, A.; Uhlenbrook, S.; Green, A.; Akkharath, I. Impact of climate change on sediment yield in the Mekong River basin: A case study of the Nam Ou basin, Lao PDR. J. Hydrol. Earth Syst. Sci. 2013, 17, 1–20. [Google Scholar] [CrossRef]
  14. Omwoyo, M.A.; Muthama, A.N.J.; Opere, A.; Onwonga, R. Simulating Streamflow in Response to Climate Change in the Upper Ewaso Ngiro Catchment, Kenya. J. Clim. Chang. Sustain. 2017, 1, 12–29. [Google Scholar]
  15. Elena, S.; Chiara, B.; Alexandre, S.G.; John, J.H. Mitigating Climate Change in the Cultural Built Heritage Sector. Climate 2019, 7, 90. [Google Scholar] [CrossRef]
  16. Li, Z.; Fang, H. Impacts of climate change on water erosion: A review. Earth-Sci. Rev. 2016, 163, 94–117. [Google Scholar] [CrossRef]
  17. Pumo, D.; Caracciolo, D.; Viola, F.; Noto, L.V. Climate change effects on the hydrological regime of small non-perennial river basins. Sci. Total Environ. 2016, 542, 76–92. [Google Scholar] [CrossRef] [PubMed]
  18. Mohammad, R.; Ommolbanin, B.; Thomas, P.; Elham, R.S. Modeling the Impact of Climate Change and Land Use Change Scenarios on Soil Erosion at the Minab Dam Watershed. Sustainability 2019, 11, 3353. [Google Scholar] [CrossRef]
  19. Zhang, S.; Li, Z.; Lin, X.; Zhang, C. Assessment of Climate Change and Associated Vegetation Cover Change on Watershed-Scale Runoff and Sediment Yield. Water 2019, 11, 1373. [Google Scholar] [CrossRef]
  20. Nearing, M.; Pruski, F.; O’Neal, M. Expected climate change impacts on soil erosion rates: A review. J. Soil Water Conserv. 2004, 59, 43–50. [Google Scholar]
  21. Peizhen, Z.; Molnar, P.; Downs, W. Increased sedimentation rates and grain sizes 2-4 Myr ago due to the influence of climate change on erosion rates. Nature 2001, 410, 891–897. [Google Scholar] [CrossRef] [PubMed]
  22. Pruski, F.; Nearing, M. Climate-induced changes in erosion during the 21st century for eight US locations. J. Water Resour. 2002, 38, 1–34. [Google Scholar]
  23. Mullan, D.; Mortlock, D.; Fealy, R. Addressing key limitations associated with modelling soil erosion under the impacts of future climate change. J. Agric. For. Meteorol. 2012, 156, 18–30. [Google Scholar] [CrossRef] [Green Version]
  24. Woldemariam, G.W.; Iguala, A.D.; Tekalign, S.; Reddy, R.U. Spatial Modeling of Soil Erosion Risk and Its Implication for Conservation Planning: The Case of the Gobele Watershed, East Hararghe Zone, Ethiopia. Land 2018, 7, 25. [Google Scholar] [CrossRef]
  25. Francipane, A.; Fatichi, S.; Ivanov, V.Y.; Noto, L.V. Stochastic assessment of climate impacts on hydrology and geomorphology of semiarid headwater basins using a physically based model. J. Geophys. Res. Earth Surf. 2015, 120, 507–533. [Google Scholar] [CrossRef] [Green Version]
  26. Julien, P. Erosion and Sedimentation, 2nd ed.; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  27. Lili, W.; Keith, A.; Dennis, C. Impacts of Climate Change on Soil Erosion in the Great Lakes Region. Water 2018, 10, 715. [Google Scholar] [CrossRef]
  28. Azari, M.; Reza, H.; Saghafian, B.; Faramarzi, M. Climate change impacts on streamflow and sediment yield in the North of Iran. Hydrol. Sci. J. 2016, 61, 123–133. [Google Scholar] [CrossRef]
  29. Adem, A.; Tilahun, S.; Ayana, E.; Worqlul, A.; Assefa, T.; Dessu, S.; Melesse, A. Climate Change Impact on Sediment Yield in the Upper Gilgel Abay Catchment, Blue Nile Basin, Ethiopia. Landscape Dynamics, Soils and Hydrological Processes in Varied Climates. Springer Geogr. 2015. [Google Scholar] [CrossRef]
  30. Admassie, A.; Adenew, A. Stakeholders’ Perceptions of Climate Change and Adaptation Strategies in Ethiopia; EEA Research Report; Ethiopian Economic Association: Addis Ababa, Ethiopia, 2008; Available online: http://www.eeaecon.org/Research%20Materials.htm (accessed on 9 June 2017).
  31. Awass, A. Hydrological Drought Analysis-Occurrence, Severity, Risks: The Case of Wabi Shebelle River Basin, Ethiopia. Ph.D. Thesis, University at Siegen, Siegen, Germany, 2009. [Google Scholar]
  32. Hamza, I.; Iyela, A. Land use pattern, Climate change, and its implication for food security in Ethiopia. Ethiop. J. Environ. Stud. Manag. 2012, 5, 26–31. [Google Scholar] [CrossRef]
  33. Deressa, T.; Hassan, R.; Cringler, C. Climate change and agriculture paper perception of and adaptation to climate change by farmers in the Nile basin of Ethiopia. J. Agric. Sci. 2011, 149, 23–31. [Google Scholar] [CrossRef]
  34. Cheung, H.; Gabriel, B.; Ashbindu, S. Trends and spatial distribution of annual and seasonal rainfall in Ethiopia. Int. J. Climatol. 2008, 28, 1723–1734. [Google Scholar] [CrossRef]
  35. Zhang, X.; Nearing, M. Impact of climate change on soil erosion, runoff, and wheat productivity in central Oklahoma. Catena 2005, 61, 185–195. [Google Scholar] [CrossRef]
  36. African Climate Change Resilience Alliance (ACCRA). Climate Trends in Ethiopia: Summary of ACCRA Research in Three Sites. 2011. Available online: https://insights.careinternational.org.uk/media/k2/attachments/ACCRA_Policy_Brief_Ethiopia_Climate_Trends.pdf (accessed on 2 February 2018).
  37. Kefyalew, A. Integrated Flood Management Case Study: Ethiopia Integrated Flood Management; WMO/GWP Associated Programme on Flood Management; World Meteorological Organization: Geneva, Switzerland, 2003. [Google Scholar]
  38. Yidenkachew, A. Flood Frequency Analysis for Lower Awash Sub-Basin (Tributaries from northern Wollo High Lands) Using Swat 2005 Model. Master’s Thesis, Addis Ababa University, Addis Ababa, Ethiopia, 2008. [Google Scholar]
  39. Uqubai, B. Flood Forecasting and Early Warning System in Lower Awash Sub-Basin. Master’s Thesis, Arba Minch University, Arba Minch, Ethiopia, 2010. [Google Scholar]
  40. Getahun, Y.; Gebre, S. Flood Hazard Assessment and Mapping of Flood Inundation Area of the Awash River Basin in Ethiopia using GIS and HEC-GeoRAS/HEC-RAS Model. J. Civ. Environ. Eng. 2015, 5, 179. [Google Scholar]
  41. Simenew, K.; Dejen, T.; Tesfaye, S.; Fekadu, R.; Tesfu, K.; Fufa, D. Characterization of Camel Production System in Afar Pastoralists, North East Ethiopia. Asian J. Agric. Sci. 2013, 5, 16–24. [Google Scholar] [CrossRef]
  42. Simane, B.; Beyene, H.; Deressa, D.; Kumie, A.; Berhane, K.; Samet, J. Review of Climate Change and Health in Ethiopia: Status and Gap Analysis. Ethiop. J. Health Dev. 2016, 30, 28–41. [Google Scholar] [PubMed]
  43. Intergovernmental Panel on Climate Change (IPCC). The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change; Houghton, J.T., Ding, Y., Griggs, D.J., Noguer, M., van der Winden, P.J., Dai, K., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2001; p. 881. [Google Scholar]
  44. Lu, X. Sediment loads response to climate change: A preliminary study of eight large Chinese rivers. Int. J. Sediment Res. 2013, 28, 1–14. [Google Scholar] [CrossRef]
  45. Monte, R.; Nearing, A.; Roel, C.; Jane, S.; Rebecca, A. Climate change impacts on soil erosion in Midwest United States with changes in crop management. Catena 2005, 61, 165–184. [Google Scholar]
  46. Food and Agriculture Organization (FAO). The Digital Soil Map of the World, Land and Water Development Division; FAO: Rome, Italy, 1995; Available online: http://www.fao.org/geonetwork/ (accessed on 15 March 2016).
  47. National Meteorological Agency (NMA). Mean Monthly Rainfall Data; NMA: Addis Ababa, Ethiopia, 2017; Unpublished Document.
  48. McCuen, R. Hydrologic Analysis and Design, 2nd ed.; Pearson: London, UK, 1998. [Google Scholar]
  49. Raes, D.; Mallants, D.; Son, Z. RAINBOW: A software package for analyzing hydrologic data. Trans. Ecol. Environ. 1996, 12. [Google Scholar]
  50. Subramanya, K. Engineering Hydrology; McGraw-Hill: New York, NY, USA, 2008. [Google Scholar]
  51. Neitsch, S.; Arnold, J.; Kiniry, J.; Williams, J.; King, K. Soil and Water Assessment Tool; Theoretical Documentation, Version 2005; Texas A&M University System: College Station, TX, USA, 2005. [Google Scholar]
  52. Lydia, A.; Joseph, O.; Mohammed, Y.; Christopher, O. Projected Climatic and Hydrologic Changes to Lake Victoria Basin Rivers under Three RCP Emission Scenarios for 2015–2100 and Impacts on the Water Sector. Water 2019, 11, 1449. [Google Scholar] [CrossRef]
  53. Endris, H.S.; Lennard, C.; Hewitson, B.; Dosio, A.; Nikulin, G.; Panitz, H.J. Teleconnection responses in multi-GCM driven CORDEX RCMs over Eastern Africa. Clim. Dyn. 2016, 46, 2821–2846. [Google Scholar] [CrossRef]
  54. Haile, A.; Rietjes, T. Evaluation of regional climate model simulations of rainfall over the upper Blue Nile basin. Atmos. Res. 2015, 161, 57–64. [Google Scholar]
  55. Bakker, A.; Hurk, B. Bias Correction and Resampling of RACMO Output for the Hydrological Modelling of the Rhine; Technical Report TR-307; KNMI: De Bil, The Netherlands, 2011. [Google Scholar]
  56. Gadissa, T.; Nyadawa, M.; Behulu, F.; Mutua, B. The Effect of Climate Change on Loss of Lake Volume: Case of Sedimentation in Central Rift Valley Basin, Ethiopia. Hydrology 2018, 5, 67. [Google Scholar] [CrossRef]
  57. Wörner, V.; Kreye, P.; Meon, G. Effects of Bias-Correcting Climate Model Data on the Projection of Future Changes in High Flows. Hydrology 2019, 6, 46. [Google Scholar] [CrossRef]
  58. Graham, L.P.; Haggemann, S.; Jaun, S.; Beniston, M. On interpreting hydrological change from regional climate models. Clim. Chang. 2007, 81, 97–122. [Google Scholar] [CrossRef]
  59. Leander, R.; Buishand, T. Resampling of regional climate model output for the simulation of extreme river flows. J. Hydrol. 2007, 332, 487–496. [Google Scholar] [CrossRef]
  60. Teutschbein, C.; Seibert, J. Bias correction of regional climate model (RCM) simulations possible for non-stationary conditions? J. Hydrol. Earth Syst. Sci. 2013, 17, 5061–5077. [Google Scholar] [CrossRef]
  61. Chaemiso, A.; Abebe, A.; Pingale, S. Assessment of the impact of climate change on surface hydrological processes using SWAT: A case study of Omo-Gibe river basin, Ethiopia. Model. Earth Syst. Environ. 2016, 2, 205. [Google Scholar] [CrossRef]
  62. Terink, W.; Hurkmans, R.; Torfs, P.; Uijlenhoet, R. Evaluation of a bias correction method applied to downscaled precipitation and temperature reanalysis data for Rhine Basins. J. Hydrol. Earth Syst. Sci. 2009, 6, 5377–5413. [Google Scholar] [CrossRef]
  63. Ho, C.; Stephenson, D.; Collins, M.; Ferro, C.; Brown, S. Calibration strategies: A Source of Additional Uncertainty in Climate Change Projections. Bull. Am. Meteorol. Soc. 2012, 93, 21–26. [Google Scholar] [CrossRef]
  64. Mann, H.B. Non-parametric tests against trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  65. Kendall, M.G. Rank Correlation Methods, 4th ed.; Charles Griffin: London, UK, 1975; ISBN 978-0195208375. [Google Scholar]
  66. Amogne, A.; Belay, S.; Ali, H.; Amare, B. Variability and time series trend analysis of rainfall and temperature in northcentral Ethiopia: A case study in Woleka sub-basin. Weather Clim. Extrem. 2018, 19, 29–41. [Google Scholar]
  67. Elsa, C.; Andrés, M.; José, A.; Vincenzo, L. East Africa Rainfall Trends and Variability 1983–2015 Using Three Long-Term Satellite Products. Remote Sens. 2018, 10, 931. [Google Scholar] [CrossRef]
  68. Adeyeri, O.; Lamptey, B.; Lawin, A.; Sanda, I. Spatio-Temporal Precipitation Trend and Homogeneity Analysis in KomaduguYobe Basin, Lake Chad Region. J. Climatol. Weather Forecast. 2017, 5, 214. [Google Scholar]
  69. Birhanu, H.; Yingjun, C.; Zinabu, M.; Miseker, B. Temperature and Precipitation Trend Analysis over the last 30 years in Southern Tigray regional State, Ethiopia. J. Atmos. Pollut. 2017, 5, 18–23. [Google Scholar]
  70. Jie, W.; Hiroshi, I.; Wenchao, S.; Shaowei, N. Development and Interpretation of New Sediment Rating Curve Considering the Effect of Vegetation Cover for Asian Basins. Sci. World J. 2013, 2013, 9. [Google Scholar] [CrossRef]
  71. Nikita, I.T. Fitting sediment rating curves using regression analysis: A case study of Russian Arctic rivers. (In Proceedings of the Sediment Dynamics from the Summit to the Sea, New Orleans, LA, USA, 11–14 December 2014). Proc. IAHS 2014, 367, 193–198. [Google Scholar]
  72. Mamaru, A.M.; Fasikaw, A.Z.; Muluken, L.A.; Getaneh, K.A.; Dessalegn, C.D.; Seifu, A.T.; Tammo, S. Sediment concentration rating curves for a monsoonal climate: Upper Blue Nile. Soil 2016, 2, 337–349. [Google Scholar]
  73. Manfred, K.; Netsanet, C. SWAT-Modeling of the Impact of future Climate Change on the Hydrology and the Water Resources in the Upper Blue Nile River Basin, Ethiopia. In Proceedings of the ICWRER, Koblenz, Germany, 3–7 June 2013; pp. 488–523. [Google Scholar]
  74. Arnold, J.; Williams, J.; Maidment, D. Continuous time water and sediment routing model for large basins. J. Hydraul. Eng. 1995, 121, 171–183. [Google Scholar] [CrossRef]
  75. Abdi, B.; Desalegn, C.; Ekasit, K. Simulation of Sediment Yield using SWAT Model in Fincha Watershed, Ethiopia. Kasetsart J. Nat. Sci. 2012, 46, 283–297. [Google Scholar]
  76. Arnold, J.G.; Fohrer, N. SWAT2000: Current capabilities and research opportunities in applied watershed modeling. Hydrol. Proc. 2005, 19, 563–572. [Google Scholar] [CrossRef]
  77. Hargreaves, G.; Riley, J. Agricultural benefits for Senegal River Basin. J. Irrig. Drain. Eng. 1985, 111, 113–124. [Google Scholar] [CrossRef]
  78. Priestley, C.; Taylor, R. On the assessment of surface heat flux and Evaporation using large-scale parameters. Weather Rev. 1972, 100, 81–92. [Google Scholar] [CrossRef]
  79. Monteith, J. Evaporation and the environment in the state and movement of water in living organisms. In 19th Symposia of the Society for Experimental Biology; Cambridge University Press: London, UK, 1965. [Google Scholar]
  80. Abbaspour, K.; Yang, I.; Maximov, R.; Siber, K.; Bogner, J.; Mieleitner, J.; Zobrist, R. Modeling hydrology and water quality in the pre-alpine/alpine Thur watershed using SWAT. J. Hydrol. 2007, 333, 413–430. [Google Scholar] [CrossRef]
  81. Santhi, C.; Arnold, J.; Williams, W.; Dugas, R.; Srinivasan, M.; Hauck, L. Validation of the SWAT model on a large river basin with point and nonpoint sources. J. Am. Water Resour. Assoc. 2001, 37, 1169–1188. [Google Scholar] [CrossRef]
  82. Moriasi, D.; Arnold, J.; Van, M.; Binger, R.; Harmel, R.; Veith, T. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. TASABE 2007, 50, 885–900. [Google Scholar]
  83. Atlas (Climate Change Adaptation, Thought Leadership and Assessments). Climate Variability and Change in Ethiopia; Technical Report; Ask Order No. AID-OAA-I-14-00013, under the Restoring the Environment through Prosperity, Livelihoods, and Conserving Ecosystems (REPLACE) IDIQ; United States Agency for International Development: Washington, DC, USA, 2015.
  84. Jury, M.; Funk, C. Climatic trends over Ethiopia: Regional signals and drivers. Int. J. Climatol. 2013, 33, 1924–1935. [Google Scholar] [CrossRef]
  85. Guo, D.; Sun, J.; Yu, E. Evaluation of CORDEX regional climate models in simulating temperature and precipitation over the Tibetan Plateau. Atmos. Ocean. Sci. Lett. 2018, 11, 219–227. [Google Scholar] [CrossRef] [Green Version]
  86. Getahun, Y.S. Impact of Climate Change on Hydrology of the Upper Awash River Basin (Ethiopia): Inter-Comparison of Old SRES and New RCP Scenarios. Master’s Thesis, Wageningen University, Wageningen, The Netherlands, 2016. [Google Scholar]
  87. Daba, M.; Tadele, K.; Shimelis, A.; Melesse, A. Hydrological Response to climate change impacts and Adaptation Strategies of Upper Awash Sub-basin, in the Awash Basin, Ethiopia. In Proceedings of the International Conference on Climate-Change Impacts for Scientists & Stakeholders, Postdam, Germany, 11–13 October 2017. [Google Scholar]
  88. Kinfe, H. Impact climate change on the water resources of Awash River Basin. Clim. Res. 1999, 12, 91–96. [Google Scholar]
  89. National Meteorological Agency (NMA). National Meteorology Agency; Final Report on Evaluation Criteria for Identifying High Priority Adaptation Activities prepared by B and M Development Consultants for NMA; NMA: Addis Ababa, Ethiopia, 2007.
  90. McSweeney, C.; New, M.; Lizcano, G. UNDP Climate Change Country Profiles-Ethiopia; UNT Libraries: Denton, TX, USA, 2008. [Google Scholar]
  91. Benaman, J.; Christine, G.; Douglas, A. Calibration and validation of soil and water Assessment Tool on an Agricultural watershed in Upstate New York. J. Hydrol. Eng. ASCE 2005, 10, 363–374. [Google Scholar] [CrossRef]
  92. Sufiyan, A. Sediment Yield Estimation from Upper Awash Watershed: Sedimentation into Koka Reservoir. Master’s Thesis, Arba Minch University, Arba Minch, Ethiopia, 2014. [Google Scholar]
  93. Engidayehu, C. SWAT Based Estimation of Sediment Yield from Guder Watershed, Abbay Basin, Ethiopia. Master’s Thesis, Arba Minch University, Arba Minch, Ethiopia, 2014. [Google Scholar]
  94. Ayele, G.T.; Teshale, E.Z.; Yu, B.; Rutherfurd, I.D.; Jeong, J. Streamflow and Sediment Yield Prediction for Watershed Prioritization in the Upper Blue Nile River Basin, Ethiopia. Water 2017, 9, 782. [Google Scholar] [CrossRef]
  95. Nerantzaki, S.D.; Giannakis, G.V.; Efstathiou, D.; Nikolaidis, N.P.; Sibetheros, I.A.; Karatzas, G.P.; Zacharias, I. Modeling suspended sediment transport and assessing the impacts of climate change in a karstic Mediterranean watershed. Sci. Total. Environ. 2015, 538, 288–297. [Google Scholar] [CrossRef]
  96. Chimdessa, K.; Quraishi, S.; Kebede, A.; Alamirew, T. Effect of Land Use Land Cover and Climate Change on River Flow and Soil Loss in Didessa River Basin, South West Blue Nile, Ethiopia. Hydrology 2019, 6, 2. [Google Scholar] [CrossRef]
  97. Hurni, H. Soil Conservation Manual for Ethiopia; Ministry of Agriculture: Addis Ababa, Ethiopia, 1985.
  98. Dao, N.; Tadashi, S. Assessment of climate change impacts on hydrology and sediment yield in the river catchment, Vietnam. Hydraul. Eng. 2013, 69, 31–36. [Google Scholar]
  99. Cousino, L.K.; Becker, R.H.; Zmijewski, K.A. Modeling the effects of climate change on water, sediment, and nutrient yields from the Maumee River watershed. Journal of Hydrology: Regional Studies. J. Hydrol. Reg. Stud. 2015, 4, 762–775. [Google Scholar]
Hydrology 06 00081 g001 550
Figure 1. Location of the Logiya watershed, meteorological and flow gauging stations.
Figure 1. Location of the Logiya watershed, meteorological and flow gauging stations.
Hydrology 06 00081 g001
Hydrology 06 00081 g002 550
Figure 2. Mean monthly precipitation and temperature (Tmax and Tmin) selected stations in the study area (1988–2016).
Figure 2. Mean monthly precipitation and temperature (Tmax and Tmin) selected stations in the study area (1988–2016).
Hydrology 06 00081 g002
Hydrology 06 00081 g003 550
Figure 3. Soil types (a) and Land Use/Land Cover (LULC) classes (b) in Logiya watershed.
Figure 3. Soil types (a) and Land Use/Land Cover (LULC) classes (b) in Logiya watershed.
Hydrology 06 00081 g003
Hydrology 06 00081 g004 550
Figure 4. Sediment rating curve of the Logiya River at the Logiya gauge station.
Figure 4. Sediment rating curve of the Logiya River at the Logiya gauge station.
Hydrology 06 00081 g004
Hydrology 06 00081 g005a 550Hydrology 06 00081 g005b 550
Figure 5. Historical trends of the annual precipitation (a), temperature maximum (Tmax) (b), temperature minimum (Tmin) (c), and sediment yield (d) from 1971–2000 in the Logiya watershed.
Figure 5. Historical trends of the annual precipitation (a), temperature maximum (Tmax) (b), temperature minimum (Tmin) (c), and sediment yield (d) from 1971–2000 in the Logiya watershed.
Hydrology 06 00081 g005aHydrology 06 00081 g005b
Hydrology 06 00081 g006 550
Figure 6. Observed and downscaled average monthly maximum temperature (a), minimum temperature (b), and precipitation (c) of observed vs. bias-corrected and un-corrected RCP time series (1988–2016) of the Logiya watershed.
Figure 6. Observed and downscaled average monthly maximum temperature (a), minimum temperature (b), and precipitation (c) of observed vs. bias-corrected and un-corrected RCP time series (1988–2016) of the Logiya watershed.
Hydrology 06 00081 g006
Hydrology 06 00081 g007 550
Figure 7. Soil and Water Assessment Tool (SWAT) model calibration (a) and validation (b) for streamflow at the monthly time scale during 1992–1999 and 2000–2005, respectively.
Figure 7. Soil and Water Assessment Tool (SWAT) model calibration (a) and validation (b) for streamflow at the monthly time scale during 1992–1999 and 2000–2005, respectively.
Hydrology 06 00081 g007
Hydrology 06 00081 g008 550
Figure 8. SWAT model calibration (a) and validation (b) for sediment yield at the monthly time scale during 1992–1999 and 2000–2005 respectively.
Figure 8. SWAT model calibration (a) and validation (b) for sediment yield at the monthly time scale during 1992–1999 and 2000–2005 respectively.
Hydrology 06 00081 g008
Hydrology 06 00081 g009 550
Figure 9. Mean monthly observed and simulated sediment load (1992–2005).
Figure 9. Mean monthly observed and simulated sediment load (1992–2005).
Hydrology 06 00081 g009
Hydrology 06 00081 g010 550
Figure 10. Spatial map of the observed sediment yield.
Figure 10. Spatial map of the observed sediment yield.
Hydrology 06 00081 g010
Hydrology 06 00081 g011 550
Figure 11. Spatial map of the annual sediment yield (a) in 2030s RCP4.5 (b) in 2030s RCP8.5 (c) in 2060s RCP4.5 and (d) in 2060s RCP8.5 in the Logiya watershed.
Figure 11. Spatial map of the annual sediment yield (a) in 2030s RCP4.5 (b) in 2030s RCP8.5 (c) in 2060s RCP4.5 and (d) in 2060s RCP8.5 in the Logiya watershed.
Hydrology 06 00081 g011
Hydrology 06 00081 g012 550
Figure 12. Projected change of mean monthly sediment yield (a) in the 2030s (b) in the 2060s from the baseline period in the Logiya watershed.
Figure 12. Projected change of mean monthly sediment yield (a) in the 2030s (b) in the 2060s from the baseline period in the Logiya watershed.
Hydrology 06 00081 g012
Hydrology 06 00081 g013 550
Figure 13. Mean monthly sediment yield in the 2030s (a) and the 2060s (b) projected at the outlet of the Logiya watershed.
Figure 13. Mean monthly sediment yield in the 2030s (a) and the 2060s (b) projected at the outlet of the Logiya watershed.
Hydrology 06 00081 g013
Table
Table 1. Meteorological stations near the Logiya watershed and percentage of missing data.
Table 1. Meteorological stations near the Logiya watershed and percentage of missing data.
StationsLatitudeLongitudeElevation (m)Missing (%)
Dubti11.7241.013766.62
Mersa12.1339.6314707.72
Kobbo11.6639.67157814.75
Bati11.240.0216604.02
Table
Table 2. Statistical summary of annual baseline data (1971–2000).
Table 2. Statistical summary of annual baseline data (1971–2000).
Descriptive Statistics Precipitation (mm)Tmax (°C)Tmin (°C)Sediment Yield (tons)
Max304.531.717.134547.3
Min58.030.215.5381.6
SD54.30.40.47001.4
Mean125.431.116.44528.4
CV (%)43.31.32.3154.6
p-value0.2720.1010.1170.299
Alpha0.050.050.050.05
SD: standard deviation, CV: coefficient of variance.
Table
Table 3. Performance evaluation of observed and downscaled historical representative concentration pathway (RCP) datasets of the average annual precipitation in the Logiya watershed (1988–2016).
Table 3. Performance evaluation of observed and downscaled historical representative concentration pathway (RCP) datasets of the average annual precipitation in the Logiya watershed (1988–2016).
Mean Annual Precipitation (mm)CV (%)Bias (%)Correlation
Observed 742.462.60--
Bias-corrected RCP733.653.12−1.190.56
Un-corrected RCP408.862.32−44.930.34
Table
Table 4. Result of the sensitivity analysis parameters of streamflow [sensitivity parameters maximum, minimum, and fitted values using Sequential Uncertainty Fitting version-2 (SUFI-2)].
Table 4. Result of the sensitivity analysis parameters of streamflow [sensitivity parameters maximum, minimum, and fitted values using Sequential Uncertainty Fitting version-2 (SUFI-2)].
RankParametersDescriptions Fitted ValueMinMax
1ALPHA_BFBase flow alpha factor (days)0.78501
2CH_K2Effective hydraulic conductivity of the main channel77.49−0.05500
3GW_DELAYGroundwater delay (days)362.50500
4CN-2SCS (soil conservation service) runoff curve number for moisture condition—II0.1−0.25+0.25
5REVAPMNThreshold depth of water in the shallow aquifer for “revap” to occur (mm)477.010500
6RCHRG-DPDeep aquifer percolation fraction0.38501
7SOL_KSaturated hydraulic conductivity (mm/hr)1546.7802000
8GWQMNThe threshold depth of water in shallow required for return flow (mm)2134.8705000
9SURLAGSurface runoff lag time16.890.0524
10SOL_AWCSoil available water capacity0.0501
11EPCOSoil evaporation compensation factor0.6901
12SOL_ALBMoist soil albedo0.100.25
Table
Table 5. Streamflow calibration and validation results (monthly) at the Logiya outlet gauging station
Table 5. Streamflow calibration and validation results (monthly) at the Logiya outlet gauging station
VariableCalibrationValidation
R20.80.77
NSE0.730.63
PBIAS+30.5+38.1
NSE: Nash–Sutcliffe efficiency; PBIAS: percent bias.
Table
Table 6. Sensitivity analysis parameters of Sediment yield with fitted values using SUFI-2.
Table 6. Sensitivity analysis parameters of Sediment yield with fitted values using SUFI-2.
RankParametersDescriptionsFitted ValueMinMax
1SPCONLinear factor for channel sediment routing0.0050.00010.01
2SPEXPExponential factor for channel sediment routing1.38712
3USLE-PUSLE support Practice factor0.27501
SPCON: maximum amount of sediment; SPEXP: sediment re-entrained channel sediment routing; USLE-P: universal soil loss equation support particle factor.
Table
Table 7. Monthly sediment yield calibration and validation results at the Logiya outlet gauging
Table 7. Monthly sediment yield calibration and validation results at the Logiya outlet gauging
VariableCalibrationValidation
R20.830.85
NSE0.790.76
PBIAS−23.4−25.0

Share and Cite

MDPI and ACS Style

Jilo, N.B.; Gebremariam, B.; Harka, A.E.; Woldemariam, G.W.; Behulu, F. Evaluation of the Impacts of Climate Change on Sediment Yield from the Logiya Watershed, Lower Awash Basin, Ethiopia. Hydrology 2019, 6, 81. https://doi.org/10.3390/hydrology6030081

AMA Style

Jilo NB, Gebremariam B, Harka AE, Woldemariam GW, Behulu F. Evaluation of the Impacts of Climate Change on Sediment Yield from the Logiya Watershed, Lower Awash Basin, Ethiopia. Hydrology. 2019; 6(3):81. https://doi.org/10.3390/hydrology6030081

Chicago/Turabian Style

Jilo, Nura Boru, Bogale Gebremariam, Arus Edo Harka, Gezahegn Weldu Woldemariam, and Fiseha Behulu. 2019. "Evaluation of the Impacts of Climate Change on Sediment Yield from the Logiya Watershed, Lower Awash Basin, Ethiopia" Hydrology 6, no. 3: 81. https://doi.org/10.3390/hydrology6030081

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop