Hydro-Meteorological Trends in the Upper Omo-Ghibe River Basin, Ethiopia

This study aimed to assess trends of hydro-meteorological variables in the Upper Omo-Ghibe river basin, Ethiopia. Data records from eleven rainfall, eight air temperature, and five streamflow stations between 1981 to 2008 were investigated. The trends and change points were evaluated for different periods of time depending on data availability. Mann-Kendall and Pettit tests were used to identify trends and change points at a 5% significance level. The tests were applied to mean annual, monthly and seasonal time scales. Rainfall exhibited statistically decreasing trends at a mean annual time scale, while seasonal rainfall depicted heterogeneous results in both directions. For the majority of the stations, air temperature showed statistically significant increasing trends. The magnitude of change in temperature for mean annual, wet and dry season has increased about 0.48, 0.46, and 0.61 ◦C per decade for Jimma station. The Pettit test revealed that the late 1980s and 1990s were the change points. There is generally a decreasing trend in streamflow. The decline in annual rainfall and rise in temperature affected the streamflow negatively. Overall, the results indicate that trend sand change point times varied considerably across the stations and catchments. The identified significant trends can help to support planning decisions for water management.


Introduction
Climate change may lead to alterations in the dynamics of hydro-climatic variables. These changes can have tremendous effects for food, energy and water supply in particular in developing countries. The release of the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report [1] has brought the issue of climate change and climate trends to the forefront of scientific and political communities. Investigations of regional and global climatic changes and variability and their impacts on the society have designed and operated based on the historical pattern of water availability, quality and demand assuming constant climatic behavior.
There are many studies analyzing the trends and the variability of hydroclimatic variables in different parts of the world. For example, [2] studied regional analysis of trend and step changes observed in hydroclimatic variables around the Colorado River Basin. Their results indicated that temperature increased persistent throughout the year across the Colorado River Basin whereas rainfall only notably increased in a few climate divisions during February and remained relatively unchanged. Similarly Abdul Aziz and Burn [3] studied trends and variability of the hydrological regime in the Mackenzie River Basin in northern Canada over different time periods. Their results indicated

Study Area
The study area is the upper Omo-Ghibe River Basin (Gibe III catchment) located in the southwestern highlands of Ethiopia. According to [18], the area receives rainfall for about nine to ten months per year. The basin is located between 6°25' N and 9° 24' N latitude and between 35°36' E and 38°34' E longitude (Figure 1).
The basin is characterized by diverse topographic features with elevations ranging from 746 m asl in the southern part to 3522 m asl in the northern highlands. The rainfall in the basin has a monomodal pattern. The mean annual rainfall during the study period was 1425 mm and 92% of the annual rainfall occurred during the wet season (March-October), while the dry season (November-February) receives only 8% of the annual rainfall [19]. The mean annual temperature is 19.2 °C. The area is drained by some of the major rivers of the country, such as the Omo, Gilgel Gibe, Gojeb, and their numerous tributaries. These have created the dissected terrain. Nearly half of the country's remaining natural forests are found in this region [20].

Data Sources
Daily rainfall data for eleven stations and daily air temperature data for ten stations of the study basin from 1981 to 2008 were provided by the Ethiopian National Meteorological Agency (NMA). After analyzing temperature data with regard to record length and missing data, eight of the ten stations were chosen for further analysis. The daily streamflow data for eight stations in the upper Omo-Ghibe River Basin were collected from the hydrology department, Ministry of Water, Irrigation, and Electricity. However, only five streamflow stations were considered for analysis based on the availability and quality of the data. All data were scanned for outliers and data gaps. A correction was carried out by cross-checking records from neighboring stations. The selection of study catchments and stations was based on the availability of data. Only time series with <10% missing data for daily streamflow and monthly weather data have been used.

Methodology
Daily data were aggregated to monthly data for each dataset. These data were used to detect long-term trends in temperature, rainfall, and streamflow. In this study, the Mann-Kendall test was used to evaluate the existence of monotonic trends in rainfall, temperature, and streamflow time series in the upper Omo-Ghibe River Basin. Pettit tests were employed to detect the step change in the time series [21,22].
The Mann-Kendall test is a rank-based nonparametric test, which has been widely used to analyze trends in rainfall, streamflow and water quality [23][24][25]. It is applicable for non-normally distributed data, which are often encountered in hydrology or climatology [26]. The method is usually less affected by the presence of outliers and missing values [27]. The assumption is based on the null hypothesis meaning the data series is serially independent and identically distributed with no trend and as an alternative hypothesis, a trend exists in the data series. All trend results in this study have been evaluated at the 0.05 significance level to ensure an effective exploration of the trend characteristics in the study area. The null hypothesis of no trend is rejected, when |Z| > 1.96 in Equation (4).
The rank of data values within the time series are compared and the Mann-Kendall test statistic S is computed by using Equation (1).
where x j and x k are the annual data values in years j and k, j > k respectively and n is the number of observations.
The variance of (S) for independent and identically distributed random variables with no tied data computed using Equation (3) and S is a normal distribution with a mean of zero and variance of one. The Z score is computed using Equation (4).
The Z score value was used to evaluate the presence of a statistically significant trend. A positive value of Z indicates an upward trend and a negative value indicates a downward trend. The Mann-Kendall test should be applied to serially independent and uncorrelated data [28]. It is necessary to consider the correlation structure of the data. The correlation structure consists of the serial, or autocorrelation of the data series and the cross-correlation between hydrological variables at different locations. To correct the data for serial correlation, the procedure of trend-free pre-whitening (TFPW) involves estimating a monotonic trend for the series, while removing this trend prior to the pre-whitening the series and finally adding the monotonic trend calculating in the first step to the pre-whitening data series, which is described in Yue et al. [29], was used. The steps used in TFPW starts with detrending the original series by removing the linear trend in the data.
The magnitude of change was analyzed by estimating slope, thus β is estimated using Equation (5) and the detrending series computed using Equation (6).
For all k < j. (j = 2, . . . , n and k = 1, . . . , n−1) where Y t is detrended series, X t is the original data series value at time t and β is the slope. The second step is accomplished checking the presence of auto-correlation with lag-1 in the De-trended series using Equation (7).
where r 1 is the lag-1 autocorrelation coefficient, n is the number of years in the de-trended series, Y t is the trended series, and the over-bar indicates the mean of the trended series. The lag-1 autocorrelation coefficient at 0.05 significance level was evaluated to remove the auto-regressive part in the de-trended series [7]. If the lag-1 autocorrelation coefficient is significant, then the auto-correlation part is removed from the series using Equation (8) and if it is not significant, the Mann-Kendal test can be applied to the original X t data series where Y 1 is a series without an auto-regressive part. Finally, the linear trend is added to the new series using Equation (9) and Mann-Kendal test can be applied to this new series.
where Y 2 is a new series without auto-regressive and linear trend in the original data series.

Detection of Change Points with Pettit Test
To detect the date of a change point within the time series data, the Pettit test was used [22]. The Pettit test is a nonparametric, rank-based and distribution-free test. It is a statistical method to detect abrupt changes in time series. This test has been widely used to detect time change points in climatic and hydrological data [30,31].
The null hypothesis H 0 : No change or τ = T is tested against the alternative hypothesis Ha: Change or 1 ≤ τ < T using the non-parametric statistic K T = Max|U t , T |, 1 ≤ τ ≤ T.
The significance level associated with K T is computed using Equation (12).
P is the probability of the presence of a change point. A p value of less than 0.05 was defined as a significant upward and downward change in the time series considering a 5% significant level. In this study, a spreadsheet model has been set up to detect the change points in the data series.

Temperature Trend Analysis
The Mann-Kendall test analysis of temperature indicated statistically significant increasing trends for minimum, maximum and mean temperatures at the majority of the stations during the study period. This holds true both for mean annual and seasonal data. The spatial map also indicated that the higher the Z score signified that the stronger positive temperature change shown in Figure 2. This positive trend was manifested for the minimum and mean temperature.
The null hypothesis H0: No change or = T is tested against the alternative hypothesis Ha: Change or 1 ≤ < T using the non-parametric statistic KT = Max|Ut,T|, 1 ≤ ≤ T. (11) The significance level associated with KT is computed using Equation (12).
P is the probability of the presence of a change point. A P value of less than 0.05 was defined as a significant upward and downward change in the time series considering a 5% significant level. In this study, a spreadsheet model has been set up to detect the change points in the data series.

Temperature Trend Analysis
The Mann-Kendall test analysis of temperature indicated statistically significant increasing trends for minimum, maximum and mean temperatures at the majority of the stations during the study period. This holds true both for mean annual and seasonal data. The spatial map also indicated that the higher the Z score signified that the stronger positive temperature change shown in Figure  2. This positive trend was manifested for the minimum and mean temperature.

Detection of Change Points for Temperature Time Series
The results from Pettit test indicated that statistically significant increasing change points were identified for the majority of the stations. It is apparent that the change points occurred around the late 1980s, mid-1990s, and early 2000s (Table 1). During these years, the temperature increased about 0.67 °C-0.95 °C on average. This value is the difference between the mean before and after the change. For instance, at Assendabo station the change point for mean annual temperature was in 1995. The mean temperature before the change point was 19.36 °C and the mean temperature after the change point was 20.03 °C, thus, the change was 0.67 °C. Regarding Wolkite station, the change point for the mean annual temperature was in 1996, which means before and after the change point was 19.61 °C and 20.56 °C, the change was therefore 0.95 °C. The selection of the two stations was made due to the range of the change that is statistically significant at 0.01 and 0.05 significance level. Similar patterns were also observed for the rest of the stations.

Detection of Change Points for Temperature Time Series
The results from Pettit test indicated that statistically significant increasing change points were identified for the majority of the stations. It is apparent that the change points occurred around the late 1980s, mid-1990s, and early 2000s (Table 1). During these years, the temperature increased about 0.67 • C-0.95 • C on average. This value is the difference between the mean before and after the change. For instance, at Assendabo station the change point for mean annual temperature was in 1995. The mean temperature before the change point was 19.36 • C and the mean temperature after the change point was 20.03 • C, thus, the change was 0.67 • C. Regarding Wolkite station, the change point for the mean annual temperature was in 1996, which means before and after the change point was 19.61 • C and 20.56 • C, the change was therefore 0.95 • C. The selection of the two stations was made due to the range of the change that is statistically significant at 0.01 and 0.05 significance level. Similar patterns were also observed for the rest of the stations. Table 1. Summarized results of Pettit test for the mean annual and seasonal time scales. The first, second and third rows indicate for the minimum, maximum, and mean temperature, a positive sign in the parentheses indicate an increasing shift, and a negative sign indicates a decreasing shift. Bold numbers indicate p-values for statistically significant change points.

Stations
Mean Annual Wet Season Dry Season

Rainfall
For the Mann-Kendall test, we used all investigated datasets directly, because none of the data series was correlated at a 0.05 significance level, which was tested by a serial correlation test. Table 2 presents variations and trends of the mean annual rainfall at eleven stations. The mean annual rainfall varied between 1146 mm at Hossana station (in the southeastern part of the basin) and 1759 mm at Limugenet station (in the northwestern highlands of the basin). The interannual rainfall variability is shown by the coefficients of variation in (Table 2). They varied from 12.09% at Woliso station (in the northeastern highlands of the basin) and 29% at Bele station (in the southeastern part of the basin). See the location of the stations in Figure 1. The mean annual rainfall showed a decreasing trend in six out of eleven stations. However, only three (Assendabo, Shebe, and Sekoru) stations have experienced statistically changes evaluated at a 0.05 significance level ( Table 3). The rainfall decreased at Shebe and Sekoru, whereas it increased at Assendabo station ( Figure 4).  The mean annual rainfall showed a decreasing trend in six out of eleven stations. However, only three (Assendabo, Shebe, and Sekoru) stations have experienced statistically changes evaluated at a 0.05 significance level ( Table 3). The rainfall decreased at Shebe and Sekoru, whereas it increased at Assendabo station ( Figure 4).

Seasonal Rainfall Trend
The wet season is between March and October and the dry season is between November and February. This classification is also aligned with Korecha [32], in its homogeneity climate classification. The seasonal trend analysis indicated heterogeneous results both for the wet and the dry season. In the wet season, a weak rise of rainfall was observed in six out of eleven stations, but only Assendabo station has experienced a statistically significant increase at 0.05 significance level. The remaining station showed a decreasing but not significant trend (Table 3). However, during the dry season, a decreasing trend was observed in six out of eleven stations but the trend was only statistically significant at Shebe station ( Figure 5).

Seasonal Rainfall Trend
The wet season is between March and October and the dry season is between November and February. This classification is also aligned with Korecha [32], in its homogeneity climate classification. The seasonal trend analysis indicated heterogeneous results both for the wet and the dry season. In the wet season, a weak rise of rainfall was observed in six out of eleven stations, but only Assendabo station has experienced a statistically significant increase at 0.05 significance level. The remaining station showed a decreasing but not significant trend (Table 3). However, during the dry season, a decreasing trend was observed in six out of eleven stations but the trend was only statistically significant at Shebe station ( Figure 5).

Detection of Change Points for Annual Rainfall
The Pettit test indicated that the annual rainfall exhibited an abrupt change in different time series for various stations. The mean annual rainfall at Sekoru, Shebe, and Wolkite stations exhibited a downward change occurred in the years 1998, 1997, and 1997, respectively. Whereas, the stations Assendabo and Bele showed an upward change that occurred in the years 1987 and 1996, respectively. The change point analysis in annual rainfall shown about 45.5% of the total records. The remaining stations indicated no significant change point in rainfall. In this regard, we illustrated four selected stations having the change point of detection in late 1990 and in 1987, shown in Figure 6.
series for various stations. The mean annual rainfall at Sekoru, Shebe, and Wolkite stations exhibited a downward change occurred in the years 1998, 1997, and 1997, respectively. Whereas, the stations Assendabo and Bele showed an upward change that occurred in the years 1987 and 1996, respectively. The change point analysis in annual rainfall shown about 45.5% of the total records. The remaining stations indicated no significant change point in rainfall. In this regard, we illustrated four selected stations having the change point of detection in late 1990 and in 1987, shown in Figure 6.

Detection Change Points of Seasonal Rainfall
The Pettit test for seasonal rainfall exhibited different change points for different seasons. Accordingly, the wet season analysis indicated statistically significant downward change at two stations, Shebe, and Wolkite, occurred in 2000 and 2003, respectively. On the other hand, a significant upward change was observed at Assendabo station occurred in 1992 shown in Figure 7. Overall, the wet season analysis showed about 27.27% of the total records. In the dry season, a statistically significant downward change exhibited at Bonga station in 1997. On the contrary, significant upward change exhibited at Assendabo station in 1987 shown in Figure 7. Since the rainfall is very small in amount and distribution, the changes in the consecutive years of the recording were too much variable. The change point analysis proved that less than 18.18% of the stations showed different change points.

Detection Change Points of Seasonal Rainfall
The Pettit test for seasonal rainfall exhibited different change points for different seasons. Accordingly, the wet season analysis indicated statistically significant downward change at two stations, Shebe, and Wolkite, occurred in 2000 and 2003, respectively. On the other hand, a significant upward change was observed at Assendabo station occurred in 1992 shown in Figure 7. Overall, the wet season analysis showed about 27.27% of the total records. In the dry season, a statistically significant downward change exhibited at Bonga station in 1997. On the contrary, significant upward change exhibited at Assendabo station in 1987 shown in Figure 7. Since the rainfall is very small in amount and distribution, the changes in the consecutive years of the recording were too much variable. The change point analysis proved that less than 18.18% of the stations showed different change points.

Detection Change Points of Monthly Rainfall
The results of change point analysis across different months indicated heterogeneous change points. Overall, change in downward direction was prevalent in the months of the dry and wet season. However, few stations showed change points in an upward direction during the wet season.

Streamflow
Streamflow gauging stations are explained in Table 4, and the results of the Mann-Kendall test for the streamflow time series are summarized in Table 5. The table shows the results of the investigated streamflow trend, which is characterized by different trend directions. Overall, the mean annual and seasonal streamflow analysis indicated a decreasing but statistically insignificant trend. However, only the wet season at Assendabo streamflow station showed a significant increasing trend at the 0.05 significance level.

Detection Change Points of Monthly Rainfall
The results of change point analysis across different months indicated heterogeneous change points. Overall, change in downward direction was prevalent in the months of the dry and wet season. However, few stations showed change points in an upward direction during the wet season.

Streamflow
Streamflow gauging stations are explained in Table 4, and the results of the Mann-Kendall test for the streamflow time series are summarized in Table 5. The table shows the results of the investigated streamflow trend, which is characterized by different trend directions. Overall, the mean annual and seasonal streamflow analysis indicated a decreasing but statistically insignificant trend. However, only the wet season at Assendabo streamflow station showed a significant increasing trend at the 0.05 significance level.   The monthly streamflow analysis in some months indicated upward trends both at the Assendabo and Abelti station at a 0.05 significance level, respectively, shown in Table 5, while Sekoru, Shebe, and Wolkite station showed a mixture of upward and downward trends. However, statistically significant downward trends were observed at Shebe and Wolkite streamflow stations February, March, and December ( Table 5). The remaining months did not exhibit significant streamflow trends. The time series plots of the wet season streamflow for the selected stations are shown in Figure 8.   The monthly streamflow analysis in some months indicated upward trends both at the Assendabo and Abelti station at a 0.05 significance level, respectively, shown in Table 5, while Sekoru, Shebe, and Wolkite station showed a mixture of upward and downward trends. However, statistically significant downward trends were observed at Shebe and Wolkite streamflow stations February, March, and December ( Table 5). The remaining months did not exhibit significant streamflow trends. The time series plots of the wet season streamflow for the selected stations are shown in Figure 8.

Detection of Change Points in Streamflow Time Series
The change point analysis indicated several change points at the majority of the stations. The test result of the mean annual streamflow at Assendabo and Abelti exhibited an upward change that occurred around 1989 and 1987. While the mean annual streamflow at Sekoru showed a downward change in 1998, the remaining two streamflow stations showed a downward change, which was statistically insignificant. Moreover, the seasonal streamflow analysis at Assendabo and Abelti stations exhibited an upward change that occurred around 1987 but for the stations Sekoru, Shebe, andWolkite, a downward change that occurred around 1998, 1998, and 1989, respectively, was detected. The mean annual and seasonal streamflow change point for selected stations is illustrated in Figure 9. change in 1998, the remaining two streamflow stations showed a downward change, which was statistically insignificant. Moreover, the seasonal streamflow analysis at Assendabo and Abelti stations exhibited an upward change that occurred around 1987 but for the stations Sekoru, Shebe, andWolkite, a downward change that occurred around 1998, 1998, and 1989, respectively, was detected. The mean annual and seasonal streamflow change point for selected stations is illustrated in Figure 9.

Discussion
Generally, the test results are characterized by statistically significant increasing trends in temperature, decreasing trends in mean annual rainfall, a mixture of increasing and decreasing trends in seasonal rainfall and a decreasing streamflow trend in the annual and seasonal streamflow. The increases are more pronounced in the dry season than in the wet season ( Figure 2). The map also indicated that the higher the Z score, the stronger the positive temperature change (warming) in the basin. Similar results were reported by the National Meteorological Service Agency of Ethiopia [33]. The magnitude of change in temperature for mean annual and seasonal temperature has increased at a considerable rate (Figure 3). This result coincides with previous studies in other parts in Ethiopia [34,35], and also the Ethiopian panel for climate change [36] report indicated a temperature increment by 1.25 °C per 1 °C change of the globe in the next century. In general, most stations in the catchment

Discussion
Generally, the test results are characterized by statistically significant increasing trends in temperature, decreasing trends in mean annual rainfall, a mixture of increasing and decreasing trends in seasonal rainfall and a decreasing streamflow trend in the annual and seasonal streamflow. The increases are more pronounced in the dry season than in the wet season ( Figure 2). The map also indicated that the higher the Z score, the stronger the positive temperature change (warming) in the basin. Similar results were reported by the National Meteorological Service Agency of Ethiopia [33]. The magnitude of change in temperature for mean annual and seasonal temperature has increased at a considerable rate (Figure 3). This result coincides with previous studies in other parts in Ethiopia [34,35], and also the Ethiopian panel for climate change [36] report indicated a temperature increment by 1.25 • C per 1 • C change of the globe in the next century. In general, most stations in the catchment exhibited a change with continues increment of temperature (Table 1). The observed increase in temperature in the area has a considerable impact on dominant agricultural activities and cascaded hydropower plants in the basin.
The annual rainfall showed a wide variation across the basin and the interannual variability is higher at stations with low elevations and lower at stations with high elevations ( Table 2). The mean annual and dry season rainfall shows generally decreasing trends. This result is consistent with the previous study by Kebede and Bewket [15] that covered southwestern Ethiopia and identified a decreasing trend in annual rainfall. Our study confirmed findings by Seleshi and Zanke [13], who reported a significant declining trend in annual and June to September rainfall in the eastern, southern and southwestern parts of Ethiopia since about 1982. In addition, the Ethiopian Panel on Climate Change [36] revealed a rainfall decline in the southern part of the country.
The Mann-Kendall test analysis for seasonal rainfall trend in a wet season indicated a weak but statistically not significant rise of rainfall in the region except the station Assendabo. The seasonal rainfall has a great contribution to the annual rainfall and therefore affects runoff. Over different periods of investigation on a monthly scale, only February and April showed a significant change across the catchment. The remaining months showed different directions of trends but they were statistically not significant. This result is in agreement with Kebede and Bewket [15] and Demissie et al. [37].
The Pettit test results for the change points indicated different times that were not consistent for temperature, rainfall and streamflow data series. Conversely, a few stations (Sekoru and Assendabo) showed a similar change point for rainfall and streamflow. Accordingly, the change point analysis showed that for most of the streamflow gauging stations, the late 1980s and 1990s are the dominant change times. The rainfall variation and a declining trend of moisture have a significant impact on streamflow in the catchment, large-scale water resources development for hydropower production in the basin and different agricultural production activities since the majority is dependent on the rain-fed agricultural system [7,38].
The analysis of streamflow trends for five streamflow stations in the study area showed different directions of trends. The trend results did not show any spatial coherence, rather, contrasting trends were observed between adjacent streamflow gauging stations, particularly for the annual and wet season. For instance, Great-Ghibe at Abelti showed and increasing streamflow for most of the trend tests, while Wabi at Wolkite showed a decreasing streamflow. The same contrasting trends were observed between Bidru Awana at Sekoru and Gilgle gibe at Assendabo ( Table 5). The observed trends for most of the hydrological variables were consistent with the previous studies [38]. These heterogeneous results suggest human interventions [39] and natural causes [7]. Generally, the annual and seasonal streamflow for the investigated stations discovered a mixture of increasing and decreasing trends. However, these trends were statistically not significant. Nonetheless, a significant increasing trend was found in Gilgelgibe at Assendabo station at 0.05 significance level. In addition, the result of the monthly flow of Gojeb at Shebe, and Bidru Awana at Sekoru stations indicated a significant decreasing trend that is prevalent in the dry season. The monthly flow for Gilgelgibe at Assendabo and Great gibe at Abelti stations showed a significant increasing trend during the wet season.

Conclusions
The trend and change-point analysis was carried out for eight air temperature, eleven rainfall, and five stream gauging stations with a data length from 1981 to 2008 in the Upper Omo-Ghibe basin, Ethiopia. The data were analyzed using Mann-Kendall trend and Pettit test (change point).
The results of temperature analysis revealed that the investigated basin is getting warmer. The spatial distribution showed the minimum and mean temperature, strong trends are prevalent as compared to the maximum temperature both in mean annual and seasonal time series. About 62.5% of the total stations showed a positive trend regarding minimum temperature. The increase in temperature was, on the other hand, more articulated during the dry season. The mean annual rainfall trend analysis indicated decreasing trends for six out of eleven stations. However, only three stations (Assendabo, Shebe, and Sekoru) experienced statistically significant changes at 0.05 significance level. The seasonal trend analysis also revealed heterogeneous results both in the wet and in the dry season. The rainfall in the wet season has shown a weak rise whereas the rainfall in dry season exhibited declining trends. The monthly rainfall trend analysis showed that except February and April which had a significant decreasing and increasing trends respectively, the remaining months had no significant trends at all. The streamflow analysis identified different direction of trends. The flow in Gojeb at Shebe, Bidru Awana at Sekoru, and Wabi at Wolkit in the investigated catchments showed a decreasing trend in mean annual and seasonal time series. On the other hand, Gilgel Gibe at Assendabo and Great Gibe at Abelti exhibited an increasing trend but a statistically significant trend was observed only at Assendabo during the wet season. Monthly flow had the same pattern particularly at January, February, April, and May showed significant change at 0.05 significance level.
The Pettit test regarding change points exhibits different times, which were not consistent for temperature, rainfall and streamflow data series. For most of the streamflow gauging stations, the late 1980s and 1990s are the dominant change points. There is a tendency of decreasing trend and change-point sign for mean annual and dry seasons flow in the investigated catchment. Generally, the finding provides important information on hydro-meteorological changes in the study catchment. It is considered that the changes for the observed hydro-meteorological variables might be related to both changes in climate variability, either local or global, and land use/land cover change over the investigated periods. However, quantification of the causes attributed to the changes in stream flows, either from anthropogenic activities (land use/land cover change) or climate change or climate variability are beyond the focus of this study. Therefore, there is a need for more research relating to land-use change detection using remote sensing and statistical analysis of hydro-meteorological variables. In addition, for proper utilization of the scarce water resources in the region, it is recommended to quantify the effects of land use/land cover and climate change using rainfall-runoff modeling.
Author Contributions: D.J. designed the study, carried out data analyses and wrote the manuscript. B.A. and A.B. reviewed the statistical analysis and substantially revised the text during the review process. B.L. coordinated the study and gave substantial input to the manuscript.