Trends in climatic and hydrological parameters in the Ajora-Woybo watershed, Omo-Gibe River basin, Ethiopia

The identification of hydro-meteorological trends is essential for analyzing climate change and river discharge at the watershed level. The Ajora-Woybo watershed in Ethiopia was studied for long-term trends in rainfall, temperature, and discharge at the annual, monthly, and seasonal time scales. The rainfall and temperature data extend 1990 to 2020, whereas the discharge data span from1990 to 2015. The Pettitt and Standard Normal Homogeneity Test (SNHT) tests were used to determine homogeneity. The Mann–Kendall and Sen's slope tests, as well as numerous variability measures, were then employed for trend analysis. The degree of relationship between climate variables and river discharge was determined using Pearson correlation coefficients (r). Inhomogeneity was discovered in annual rainfall data from the Angacha and Areka stations. Rainfall and discharge showed insignificant trends over time, with increasing and decreasing variability across stations. Monthly rainfall decreasing trends were observed to be significantly falling in February and March. Rainfall and runoff increase just insignificantly during the Kiremt season. On the other hand, minimum, maximum, and mean annual temperatures showed significant trends with annual increases of 0.01 °C, 0.04 °C, and 0.025 °C, respectively. In this study, the relationships between discharge and temperature and rainfall were found to be moderate and minimal, respectively. Generally, the results of the long-term examination of the hydrological and climate parameters in the watershed show that water resources vary throughout and over time. As a result, designing strategies require due attention. Hydrological and meteorological parameters are essential for analyzing trends in water resource and climate change at the watershed level. Mann-Kendall and Sen's slope tests, along with a number of variability measures, were utilized in conjunction with the time series analysis approach for trend analysis. The analysis of the rainfall, temperature, and discharge in the watershed's data generally demonstrates how the availability of water resource varies over time. Designing suitable plans for water resource management and sustainable development in the watershed is therefore essential. Hydrological and meteorological parameters are essential for analyzing trends in water resource and climate change at the watershed level. Mann-Kendall and Sen's slope tests, along with a number of variability measures, were utilized in conjunction with the time series analysis approach for trend analysis. The analysis of the rainfall, temperature, and discharge in the watershed's data generally demonstrates how the availability of water resource varies over time. Designing suitable plans for water resource management and sustainable development in the watershed is therefore essential.

tion with the time series analysis approach for trend analysis.
• The analysis of the rainfall, temperature, and discharge in the watershed's data generally demonstrates how the availability of water resource varies over time.
Designing suitable plans for water resource management and sustainable development in the watershed is therefore essential.

Introduction
The global hydrological cycle is currently responding to the effects of global warming, which include increased atmospheric water vapor concentration and shifting precipitation patterns [1]. Hydro-meteorological characteristics have considerably improved the most common forms of natural disasters in Sub-Saharan Africa, particularly in drought and humidity [2,3]. Climate change's impact on river basins is a serious problem for Ethiopia in terms of the sustainability of water supplies and agricultural productivity [4]. The country is also sensitive to climate variability, which increases the frequency and magnitude of disasters at the local and national levels [5]. These negative effects of climate change could compound the country's existing social and economic difficulties, especially if people rely on climate-sensitive resources such as rainwater agriculture [6]. As a result, understanding climate variability and change, as well as their relationship to local water resources, is critical for economic growth and the creation of adaptation strategies [7]. Surface water resources are impacted by climate change due to changes in spatial-temporal distribution and magnitude [4], which will therefore affect the water balance in the area and the ecological environment [4,8]. Climate can fluctuate cyclically, monthly, seasonally, annually, and over decades, due to time-based variance (upward or downward trend) [9]. Concerns about the impact of climate change on the behavior of hydrological variables have prompted substantial research into climate and hydrological variable trends [10][11][12]. Furthermore, water management project design is dependent on observable meteorological and hydrological data [13]. This is because river basin integrated water resource management is critical to society's economic well-being and good ecological functioning, and robust analytical methodologies are also necessary to evaluate watershed hydrological behavior [14,15]. The varied terrain of Ethiopia, seasonal fluctuations in the Intertropical Convergence Zone (ITCZ), and accompanying hydrological circulation all have a significant impact on the country's climatic system [2,16]. This demonstrates that the country has a spatially diversified climate system, ranging from semi-arid (locally known as Bereha) in the lowlands to cool climate (Wurch) in mountainous regions with altitudes above 3350 masl. In this regard, three agro-ecologies are identified within the scope of this study: hot climate (kolla) in the lowlands, Midland (woinadega) at 2040 m above sea level in the middle portion of the watershed, and cold temperate climate (dega) in the watershed's highlands [17,18].
Because the hydrological cycle is so closely linked to atmospheric precipitation and temperature, changes in those variables usually result in changes in the hydrological cycle [19]. Furthermore, the seasonality of river flow varies greatly from river to river and is mostly controlled by the local seasonal cycle of precipitation, evaporation demand, and water travel periods from runoff source [20]. As a result, climate change influences all of the processes involved in the water cycle and is one of the primary causes of hydrological changes [21]. Climate change and water resource management have sparked global, regional, and local interest [22]. Many research on the relationship between climate change and water resources have been conducted throughout the world, including Ethiopia. Climate variability and climate change in Ethiopia are mostly represented by variability, i.e. a drop in precipitation and a propensity to increase in temperature example. In addition, precipitation and temperature patterns vary greatly by region. Every ten years, the average annual minimum and maximum temperature rise by around 0.25 °C and 0.10 °C, respectively. Mengistu et al. [23] investigated the geographical and temporal variability, as well as changes in rainfall and maximum and minimum temperatures at seasonal and yearly time scales, in the Upper Blue Nile river basin, Ethiopia, from 1980 to 2011. They discovered that the minimum temperature increased significantly in the basin's northern, central, southern, and southeastern regions during all seasons. On an annual scale, maximum and minimum temperatures climbed dramatically throughout more than one-third of the basin; nevertheless, on an annual and seasonal timeline, the western section of the basin witnessed falling trends. They also found statistically insignificant increases in annual rainfall. Recently, [24], demonstrated annual precipitation variability in the Woleka sub-basin in northern Ethiopia. They discovered intra and inter-annual fluctuation in the basin's climatic parameters. Belihu et al. [25], indicated the annual rainfall pattern showed a considerable decline of around 12 mm each year upstream, which is mostly influenced by the significant decrease in rainfall during the high season. Temperature analyses also found that the catchment, particularly upstream, was warming. The minimum temperature trend indicated a considerable increase of approximately 0.08 °C every year. The stream flow generally shows a downhill trend. Reduced annual and seasonal rainfall, combined with rising temperatures, because increased evaporation and have a direct detrimental impact on stream flow [26]. Discovered that rainfall showed statistically declining tendencies on an annual time scale in the upper reaches of the Omo-Gibe River basin, but seasonal rainfall showed diverse findings in both directions. At the majority of locations, air temperature increased in a statistically significant way. The overall trend of river flow was downward. Woldesenbet  [27] demonstrated that runoff has grown significantly at both the annual and seasonal scales at most monitoring stations in the upper Blue Nile basin. Stream flow trends in the region have increased dramatically on both an annual and seasonal basis at the majority of the gauging stations. The Omo-Gibe basin is one of Ethiopia's most important surface water resources [28]. However, due to changing climatic circumstances and competition in resource demand from the downstream hydroelectric plant, irrigation, and the investment group, the basin's water resources are in high demand. The region's climatic characteristics cause great variability in the temporal and spatial distribution of precipitation over the basin, resulting in a correspondingly high variability in runoff rate. As a result, there is a need to measure and assess past climate change and its relationship with the basin's water resources in general, and the Ajora-Woybo watershed in particular. As a result, analyzing precipitation, temperature, and hydrological behavior is critical for planning and forecasting future climate conditions [29], despite the fact that this analysis only evaluates observable and historical climate conditions over the last three decades. Long-term studies on changes in climatic and hydrological variables in such watersheds can thus aid in our understanding of the watershed's response to climatic change. As a result, more research at the local level is required to assess the nature of the trend and its level of relevance. Consequently, the goal of this article was to look at the trends and changes in climatic and hydrological parameters in the Ajora-Woybo watershed from 1990 to 2020. The objectives of the present study was to examine temporal variations and trends in rainfall, temperature and discharge variables. Policymakers, investors, planners, and even communities (other end-users) require the findings of this study in order to prepare for the predicted trend and change in the area's climate and hydrology. This research will also aid in the development of watershed and water resource-based climate change solutions.
The rest of this study is organized below. In Sect. 2, materials as well as methodologies are presented. In Sect. 3, the results and discussion of the study are presented. In Sect. 4, the conclusions of this study are drawn.

Description of the study area
The Omo-Gibe River basin covers 79,000 km 2 , with a length of 550 km and an average width of 140 km, and includes areas of two national regional states, the Southern Nation Nationalities and Peoples Region (SNNPR) and the Oromia Region. The river basin's entire annual discharge is estimated to be around 16.6 billion cubic meters on average (BMC) [28]. The basin, accounting for 14% of the Ethiopias's yearly runoff. The Walga, Wabe, Derghe and Gilgel-Gibe and Gojeb Rivers in the northeast and the Gilgel-Gibe and Gojeb Rivers in the south-west are the Basin's primary tributaries. The main tributaries from the eastern side of the middle and lower Omo-Gibe watershed are the Ajancho, Soke, Shapa, Woybo, Sana, and Deme Rivers. The basin is a closed river basin with its southern boundary formed by Lake Turkana in Kenya.
The Ajora-Woybo watershed, located on the eastern side of the middle Omo-Gibe River basin, is one of the major watersheds in the Omo-Gibe Basin. It was selected as the case study for this research. The study regions, the Ajora-Woybo watershed, contains four perennial rivers (Ajancho, Soke, Shapa and Woybo Rivers) that are flowing downstream to the lower valley eventually forms two rivers i.e. Ajora and Woybo (from which the name of watershed is given) finally joining into Gibe III dam reservoir ( Fig. 1). In addition, many small intermittent tributaries drain into the rivers. The watershed covers an area of some 1723.7 km 2 with a perimeter of 336.8 km. The watershed lies between 7 o 2′0″N and 7 o 50′0″N latitude and 37 o 30′0″E and 38 o 0′0″E longitude. The enormous irrigation dam project, which includes the proposed dam currently under construction in the watershed, is intended to result in 10,000 hectares of cultivated land, with the Ajancho and Woybo Rivers serving as the primary water sources.
The topography of the watershed, which includes the Omo, Ajancho-Soke, and Gilgel-Gibe Rivers, is exceedingly diversified, with mountainous to hilly terrain cut by sharply incised gorges. The lowest and highest of the watershed lies at an altitude greater than 731 masl and 3036 masl maximum elevation respectively. A gauging site is located at each part of the rivers (i.e. for Ajancho, Soke, Shapa and Woybo) with in the watershed. The climate of the Ajora-Woybo watershed varies in its nature from a hot arid climate in the southern part of the floodplain to a tropical humid one in the highlands that include the extreme north and northeastern part of the watershed. Between these extremes, and for most of the watershed, the climate is tropical sub humid. Each zone has its own weather patterns and agricultural production methods. The annual rainfall varies greatly from the north to the south, ranging from around 1900 mm in the north to less than 300 mm in the south. Furthermore, the rainfall regime in the watershed is bimodal. The mean annual temperature in the watershed varies from 16 °C in the highlands of the north to over 29 °C in the lowlands of the south [18].

Data pre-processing
The daily meteorological (rainfall and temperature) data from 1990 to 2020 and the hydrological data from 1990 to 2015 are collected from the National Meteorological Agency (NMA) and the Ministry of Water and Energy (MWE) office of Ethiopia (Table 1).

Filling of the missed data using multiple imputations (MI)
Multiple imputations are a filling method that provides valid statistical inferences under missing at random condition. The MD (missing data elimination techniques) rates using standard regression procedures and then combine the results of these analyses to obtain the final result, as shown in Eq. (1) [30,31]. MI using chained equations produces m imputations based on sequential imputation regression models of each variable conditioned by   [30,32]. The multiple imputations have been implemented in software such as XLSTAT [32]. MI is a widely used method in hydrology and its advantages were does not overestimate error, simulate missing data multiple times [30,32] and leads to best estimation of missing values [33]. To apply this method, we can proceed as follows: (i) Find k similar databases for each missing value; then the observed values are used to impute MD. (ii) For each MD (Px), we use data imputations Pi by applying single variable regression methods to obtain k different estimate results Ii (Px, Pi). (iii) The result Px will be obtained by combining all imputations results Ii using the average of all the k complete data values, that is: where, k observed values of the respective variable.

Pettitt test
To investigate the existence of the change points or break-point detection in the rainfall and temperature characteristics in time series data different methods can be applied. The Pettittt test was used to identify whether or not there is a point change or jump in the data series [34,35]. It is based on the Mann-Whitney two-sample test, which allows the detection of a single shift at an unknown time and independent of the distribution of data [34]. In addition, many researchers also applied Pettitt test in climatic and hydrological time series (e.g. [24][25][26]34]). It considers a sequence of random variables x 1 ,x 2 ,…, x T , with a change point at τ (x t for t = 1,2,…, τ), and a common distribution function F1(x) and x t for t = τ + 1, … T has a common distribution function and F2(x), and F1(x) ≠ F2(x). 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. To achieve the identification of change point, a statistical index U t , T is defined as follows Then, the approximate significance probability, p (t), of a change point is calculated from an Eq. 3 and the time series changing points of the rainfall and river discharge was checked.

Standard normal homogeneity test (SNHT)
The SNHT was developed by [36] to detect a change in a set of precipitation data. The test is applied to a set of ratios comparing observations from one monitoring station to the average of several stations. The ratios are then standardized. Here, the series of Xi corresponds to the standardized ratios. The null and alternative hypotheses are determined as follows: Ho: The T variables Xi follow an N (0, 1) distribution. Ha: Between times 1 and v, the variables follow an N ( μ 1 , 1) distribution, and between v + 1 and T, they follow an N(μ 2 , 1) distribution. The Alex-Anderson statistic defined by: The statistic T o derives from a calculation comparing the likelihood of the two alternative models. The model corresponding to Ha implies that μ 1 and μ 2 are estimated while determining the v parameter maximizing the likelihood.

Testing for serial correlation
The statistical tests presume that the succeeding data in the series are independent in order to identify a trend in a time series. The presence of serial correlation in the data has a significant impact on the strength of trend tests [37]. When there is a trend, a positive serial correlation causes the null hypothesis of no trend to be incorrectly rejected (type I error). Similar to this, when there is a negative serial correlation, the null hypothesis of no trend is accepted, even when it is wrong (type II error) [37,38]. Lag-1 serial correlation coefficients are computed to examine the data for serial correlation. By determining the lag-1 serial correlation coefficient 1, the time series is examined for serial correlation in a number of trend investigations. The lag-1 serial correlation coefficient (ρ 1 ) is determined as follows for every time series X i = x 1 , x 2 … x n .
(4) D ij = sgn x i − x j , the sgn x i − x j computed using equation 3, Research Article SN Applied Sciences (2023) 5:45 | https://doi.org/10.1007/s42452-022-05270-y where E(x i ) is the mean of the sample and n is the sample size The probability limits for ρ 1 on the correlogram of an independent series is given by Anderson (1941) as By contrasting the ρ 1 value with the essential Student's t-distribution values, the significance of serial correlation was assessed.
Each station contains 16 series (12 monthly, one annual and three seasonal series), and a total of 112 time series generated from seven rain gauge locations (7 × 16) were subjected to serial correlation tests. Two series out of 112 series were determined to be serially associated at a 99% confidence level, according to the study's findings, thus no need of pre-whitening test before applying the MK test.

Mann-Kendall trend detection (MK)
The non-parametric Mann-Kendall [39,40] statistics was chosen to detect the historical trends for rainfall, temperature and stream flow time series data, as it is widely used for water resource planning and management studies (e.g. [25,[41][42][43][44]). It is also used to determine whether a set of data values is increasing or decreasing over time, and whether the trend is statistically significant in either direction. Outliers are well treated by the Seasonal Kendall methods because it uses the rank of the values and the method is usually less affected by the presence of outlier's value and this makes the test not sensitive to outlier's values [25,44]. The observed rainfall, temperature and the discharge trends were analyzed in time series analysis. This is to discern how the historical rainfall and temperature varied over the time 1990 to 2020 with relation to discharge. The daily rainfall, temperature and the discharge data were first calculated as monthly rainfall and temperature (maximum, minimum, and mean) and discharge. Then monthly values were averaged and summed to obtain the seasonal and annual values. In Ethiopia, seasons are classified on the basis of its rainfall patterns. NMA (National Meteorological Agency) [18] states that there are three main seasons: Kiremt (the main rainy season from June to September), Belg (a small rainy season from February to May) and Bega (a dry season from October to January). Trend analyses were carried out on monthly, seasonal and annual basis with reference to other researchers (e.g. [25,26,32]). Yue and Wang [45] pointed out that the significance of the trend depends on the pre-specified significance level, the size of the trend, the sample size, and the amount of variation within a time series. That is, the larger the absolute size of the trend, the more meaningful the tests; the larger the sample size, the more meaningful the tests become.
The test is based on S statistics and each paired observed values x j (j > k) of the random variable will be inspected to find out whether x j > x k or x j < x k where x j and x k are the annual data values in years j and k, j > k respectively and n is number of observation.
S is normally distributed with mean zero, and Variance, Computation for Z score is obtained as: For a given level of significance α, if Z < Z 2 or Z > Z 1 2 , the null hypothesis that trend is absent can be rejected at α% significance level. Where, var(S) is the variance of Mann-Kendall statistic, S. The standardized MK test statistic (Z-score) follows the standard normal distribution with a mean of zero and a standard deviation of one. A positive value of Z indicates an increasing trend (upward trend) while a negative Z value signifies a decreasing trend (downward trend). When the p-value is less than the level of significance ( ), the null hypothesis (Ho) is rejected and a statistically significant trend exists in the rainfall, temperature and river discharge time series. If the p-value is greater than ( ), the null hypothesis (H o ) is accepted. Failing to reject the null hypothesis does not mean that there is no trend. Rather, it is a statement that the available evidence is not sufficient to conclude that there is a trend.

Sen's slope estimator test
It was used for estimating the magnitude of a trend by putting the linear rate of change and the intercept. The slope estimation is given by The median of N values of β i gives the Sen's slope of β, using following Eq. 16: where, xj and xk are the sequential data values, and n is the number of the recorded data. A positive value of β indicates an upward (increasing) trend, and a negative value indicates a downward (decreasing) trend in the time series.

Coefficient of variation
The coefficient of variation measures the overall variability of the rainfall, temperature and river discharge in the area of interest from the mean value. It is calculated as the ratio between the standard deviation and the longterm data mean as expressed in Eq. 17.
where CV is the coefficient of variation, σ is the standard deviation over the period and x is the mean hydro-climatic parameters. Generally, CV is used to classify the degree of variability of events into three: low (CV < 20), moderate (20 < CV < 30), and high (CV > 30). Higher CV shows more variation in parameters and vice versa [42,46].

Correlation analysis between climate and river discharge
Pearson correlation coefficients (r) was calculated for the degree of the linear relation between climate variables (rainfall and temperature) with river discharge. Moreover, the correlation between climate variables and runoff helps us to estimate their relationships. Because it is based on the method of covariance, it is known as the best approach of quantifying the relationship between variables of interest. It indicates the size of the association, or correlation, as well as the direction of the relationship.
= Median if N is even where r is the correlation coefficient, n is the length of the time series, and i is the number of years during the analyzed periods 1990-2020. Xi and Yi are the rainfall, temperature and discharge in the year i, respectively, and X and Y are the mean rainfall and, temperature and discharge, respectively during the studied periods. The obtained value was ranked into a weak correlation (0 <|r|≤ 0.3), a low correlation (0.3 <|r|≤ 0.5), a moderate correlation (0.5 <|r|≤ 0.8), and a strong correlation (0.8 <|r|≤ 1) [43,46]. All MK trend tests, Sen's slope, Pettitt change-point detections and variability analyses were conducted using the XLSTAT (https:// www. xlstat. com) software [24,34,47] (Fig. 2).

Pettitt test and standard normal homogeneity test (SNHT)
Homogeneity tests are used to assess the effects of nonclimatic factors such as changes in instrumentation, observing practices, station relocations, and station environments on climate time series data [48]. For this study, two types of homogeneity test (i.e., Pettitt's test and the SNHT test) were performed at 5% level of significance.
The results of the p-values for the Pettitt test are shown in Table 2. Any p-value less than the significance level of = 0.05 indicates inhomogeneity in the corresponding time series. For five stations (Bele, Bodit, Durame, Gesuba and Sodo), the null hypothesis was satisfied for all months. Two stations (Angacha and Areka) were found to be inhomogeneous as the p values were less than 0.05. The p values obtained from SNHT in the Table 2 show null hypothesis of homogeneous data was rejected for p values less than 0.05. Five stations have showed homogeneous data for all months. The data of two stations were inhomogeneous in the months. The change point detection may be related to occurrence of  Ethiopia as reported by [49]. Durame station is nearby Angacha but there has not been any significant rainfallpattern-shift been observed. This could be connected to the highest elevation area experienced the highest annual rainfall values, and the lowest elevation area saw the lowest annual rainfall values. This discovery is confirmed by the findings of [26,46], who discovered that mean annual rainfall and heights are closely associated in Ethiopia's west Hararge watershed and upper Omo-Gibe basin.

Detection of change points for annual rainfall
The Pettitt test analysis revealed that the rainfall in the upper and central watershed area shows a break point change with a different time series in Angacha and Areka stations (Fig. 3). The change point analysis in the annual rainfall showed that about 28.5% of the stations is inhomogeneous. The year of 2009 was the year of change and the mean annual rainfall was estimated about 1630 mm before shift and 1300 mm after shift at Angacha station. The year of 2013 was year of change and the mean annual rainfall was estimated to be about 1432 mm before shift and 1776 mm after shift at Areka station. Which clearly indicates inconsistency of rainfall data and this may be related to the topography and the high vegetation coverage of the area. The move of the rain gauge to a new location, significant changes in the station's neighborhood, changes in the ecosystem brought on by natural disasters like forest fires and landslides, and errors in observation starting on a specific date could all be contributing factors to the inconsistent rainfall pattern at the Angacha and Areka stations. Generally, except the northern highlands, where a major shift of rainfall was exhibited, the remaining parts, the central and southern parts revealed that there was no significant shift/change point/ of annual rainfall.

Detection of change points in stream flow time series
The change point analysis has shown that there is no change point in river discharge of the watershed at four gauging stations (Fig. 4).

Annual rainfall
The rainfall analysis from the inside and its nearby meteorological stations of Ajora-Woybo watershed was made based on the availability and reliability of the observed gauge stations. The meteorological stations such as Angacha, Areka, Bele, Bodit, Durame, Gesuba and Sodo stations were selected due to the long time availability of data. Accordingly, Angacha and Durame meteorological stations from the upper, Gesuba and Sodo stations from the lower and the remaining from central part of subwatershed have recorded reliable rainfall data for the last 31 years (Table 3 and Fig. 5).
The mean annual rainfall of the watershed during the study period was 1271.8 mm with 140.5 mm standard deviation and CV of 11%. The minimum and maximum ever-recorded rainfalls were 318 mm (in 2017 the driest year) in Gesuba station and 2161.8 mm (in 2019-the wettest year) in Bele station per year. Moreover, June to September is the major rainy season (Kiremt), during this period 48-62% of the annual rainfall has received. As indicated below (Fig. 5) in mean annual rainfall trend, the rainfall is generally homogenous and decreasing but it was insignificant rate. This result agreed with the work done by [26] in southern Ethiopia and [24] in north central Ethiopia. Their finding shows that there is a decreasing RF pattern since about 1981. Viste et al. [50] also showed that the decline in precipitation in the southern part of the country is large enough to produce trends on the national level. CV was used to classify the degree of variability of rainfall events. Based on this, the watershed rainfall is variable in the last two decades, particularly in the Angacha, Bele and Gesuba stations while the remaining stations become less variable.

Mann-Kendall test result of annual rainfall
MK test on annual rainfall data of watershed, the results are obtained in the following manner in Table 4. The MK is based on the calculation of Kendall's tau (measures of connection between two successive annual rainfall years). The analysis results shows that five of the stations successive annual rainfall years are negatively related and decreasing trend (i.e. Angacha, Bele, Bodit, Durame and Gesuba stations). The remaining two stations show an increasing trend but significant only in Areka (0.032) station. Annually, Z-score result was found to be in upward trend in two stations (Areka and Sodo) whereas the remaining five stations has downward trend ( Table 4). This study revealed a combination of upward and downward trends and also [26] reported similar results in upper Omo-Gibe Basin.

Analysis of monthly and seasonal RF
The monthly RF is highest in August and lowest in December in all of the stations. The rainfall trend was increased from February to April and June to August at all stations but its rate is insignificant. Monthly, the RF trend is significantly decreasing in February and March whereas November, December, January, February, March and April months show decreasing trend however it was insignificant. July, August, September and October months show increasing trend but it was insignificant rate (Fig. 6a).
Seasons are classified based on the rainfall pattern in Ethiopia (Fig. 6b). Seasonally, the precipitation was found to be Z = 2.78 in Kiremt, Z = −1.74 in Belg, Z = −0.26 in Bega ( Table 5). All the three seasons show there is no change of rainfall in the watershed at 95% confidence level and increasing for Kiremt and decreasing for Belg and Bega trend in the watershed. The seasonal rainfall trend has varies from Kiremt to Bega season 571.3 mm/season in Angacha station and 270.3 mm/season in Gesuba station. This insignificant rainfall trends may be due to high inter−annual variability [25]. Because rainfall in the studied area is seasonal and primarily falls during the Kiremt and Belg months, such statistical significance have a significant impact on the overall observed annual rainfall. However, if the downward trend in rainfall can last during the dry season months, it could lead to a lengthening of the dry season and a decrease in soil moisture at that time. In addition, the variation of monthly and seasonal RF and declining trend of moisture in time may have a significant impact on different rain−fed agricultural activities. A drought in the agricultural sector can also result from a drop in rainfall during the dry season, which will reduce crop production. As a result, it is crucial to keep an eye on these falling patterns since they may have both direct and indirect effects on the watershed's overall water budget.

Annual temperature
The annual temperature of the maximum (Tmax), minimum (Tmin) and mean (Tmean) along with their standard   Research Article SN Applied Sciences (2023) 5:45 | https://doi.org/10.1007/s42452-022-05270-y deviation (SD) and coefficient of variation (CV) have been statistically computed for each stations (Fig. 7). The temperature variability between the different years and the average annual temperature is 20.1 °C, while the annual mean maximum temperature reached 25.7 °C and the annual mean minimum temperature was 14 °C. The maximum temperature trend was significantly increasing at all station except at Angacha station. The minimum temperature trend was significantly increasing at all station except at Gesuba station. The mean temperature trend was significantly increasing at all station except at Angacha and Gesuba station.
The Mann-Kendall test statistics of the Tmax, Tmin, and Tmean are shown in Table 6. Analysis of temperature using the MK test shows that the climate in the watershed has generally warmed over the past three decades. Based on the trend of minimum and maximum temperatures, the magnitude of change was found to be about 0.01 °C/year to 0.04 °C/year and the mean values also proved that the watershed has warmed by about 0.025 °C/year. In addition, the minimum temperature analysis in southern part of the watershed (Gesuba station) has showed an increasing trend particularly in extremes. The reason for Gesuba station displaying a rising tendency in extremes when compared to the other stations could be attributed to the area's low elevation of 1650 m. This finding is consistent with [25], who found that Ethiopia's climate shifts from its hot and dry lowlands to its cool plateau mostly as a function of elevation. Similar results were also reported by [18] that revealed the mean minimum temperature has been increasing throughout the country even in the cool months by 0.37 °C per decade even though our finding less than the country report. The increasing trend of temperature in the area has significant impact on agricultural activities particularly in soil water demand and it has also led to the loss of more water from the watershed due to evapotranspiration. CV of Tmin highly observed in Bodit station (7.5%) among others.

Monthly and seasonal temperature
From the basic temperature data, the monthly and seasonal, of the maximum, minimum and mean temperatures along with their trend have been statistically computed for each month and for the three seasons: Kiremt, Belg and Bega (Table 7 and Fig. 8).
The monthly maximum temperature trend was increasing at all stations but with insignificant change. In addition, the minimum temperature trend showed increasing trend in most stations but decreasing trend throughout at Bodit station. The maximum monthly temperature is highest in March 31.2 °C at Gesuba station and lowest in July 21.2 °C at Bodit station. On the average, August is the coldest month of the year and March is the warmest (only slightly warmer than February). The lowest mean monthly temperature occurred in July (19.28 °C) (21.88 °C). In the case of seasonal trend analysis, increasing trend of minimum and maximum temperature was detected at all station (Fig. 9). The Mann-Kendall test statistics of the Tmax, Tmin and Tmean are given in Table 7. Towards the north, the temperature was cooler, but to the south, the temperature was warmer. The result of MK test in monthly temperature (maximum and minimum) revealed that the watershed had experienced warming. Mean temperature of Kiremt show decreasing trend which it was insignificant. The slope of the whole months shows that there is a positive value implying an increase in the mean monthly and seasonal temperature. Belay et al.
[1] revealed similar result in southern Ethiopia. These results are consistent with regional studies that revealed an upward trend in Ethiopia's yearly maximum and lowest temperatures as well as the country's average annual temperature. Due to the impact of warming on water availability and various agricultural activities, such

Annual discharge
In the watershed, a long period record has been observed at Woybo, Ajancho, Shapa and Soke hydrometric stations.   [51] and natural causes [9,26]. Based on CV analysis, the watershed annual river flow is highly variable, especially in the Ajancho station (72.5%). The annual trend of the discharge showed a varied trend across stations, which is insignificant. There have been numerous authors, who have reported increasing and decreasing trends in southern and central Ethiopia. Belihu et al. [25] reported increasing and decreasing trend across stations in the Gidabo River basin annually. Gedefaw et al. [43] informed that there was a decrease trend in the Awash River flow. Similarly [26,51] has found an increasing and decreasing trend in stream flow of rivers in upper Omo−Gibe basins from year to year.

Monthly and seasonal discharge
The average monthly and seasonal discharge of the Ajora−Woybo watershed from 1990 to 2015 is depicted in Table 9. The highest monthly discharge was seen in August, which is slightly greater than July and September, and the three months lie in Kiremt season. In January and February, the watershed gets poor amount of discharge. The mean monthly average discharge was low from December to March, and started to increase in the month of April. MK results reveal that no significant trend was found in the monthly discharge of the watershed. Sen's slope shows the months of July, August, September and October have an upward magnitude whereas the remaining months show a downward magnitude of trend. This indicates that monthly discharge generally show an upward and downward trend over several months. Seasonal trend over time at the watershed is shown in Fig. 10. As can be seen from figure, Kiremt season drains the highest amount of water, while Belg and Bega season drains the lowest amount of water with insignificant trend in both cases. The result of trend analysis in discharge across stations indicates that there is an insignificant increasing and decreasing trend and inter seasonal variability of discharge. A similar study conducted at the upper Omo−Gibe River basin has indicated that there is increasing and decreasing trend for seasonal stream flow at Gibe and Abelti gauging station [26]. Generally, the amount of surface water resources of the watershed is decreasing from time to time and this means there will exist shortage of water for irrigation and for the operation of the hydropower dam (Gibe III) which is typically dependent on Woybo, Ajancho, Soke and Shapa rivers. This increasing and decreasing trend of the discharge in this study may be attributed to a decreasing trend in rainfall at some stations and an increasing trend in temperature along with other factors. Different authors suggested the seasonal

Correlation analysis between climate and River discharge
The correlation between precipitation and river discharge was found to be positive and existence of good relation (R 2 = 0.68) in this study. In addition, the correlation between temperature and runoff was found to be negative and very weak (R 2 = −0.04). The decreasing trend of runoff is probably related to the other factors [52]. Therefore, the cause of this change needs further investigation (Fig. 11).

Conclusion
The purpose of this paper was to examine he historical trends and changes in meteorological and hydrological characteristics in the Ajora−Woybo watershed between 1990 and 2020. At each site, the annual, monthly, and seasonal variability of the parameters was investigated. At all locations, inter−annual fluctuation in rainfall, temperature, and discharge was recorded. Mann-Kendall and Sen's slope tests revealed that decreasing and increasing trends in rainfall, temperature, and discharge were seen at the observed stations.
The results show that inhomogeneity was observed in the annual rainfall data of Angacha and Areka stations, while majority of the stations were homogeneous. The rainfall decreases annually but the rate was insignificant trend. Monthly RF is highest in August and lowest in December in all stations. The temperature variability among different years and the average annual temperature is 20.1 °C, while the average annual maximum temperature is 25.7 °C and the average annual minimum temperature is 14 °C. The trend of maximum, minimum and average temperatures was significantly increasing at the watershed level. Based on the trend of minimum and maximum temperatures, the magnitude of change was found to be between 0.01 °C/year and 0.04 °C/year, and the mean values also demonstrate that the watershed is warming by about 0.025 °C/year. The annual trend of the discharge showed an increasing and decreasing tendency from time to time across stations. Seasonally, during kiremt season increasing trend of discharge was observed Research Article SN Applied Sciences (2023) 5:45 | https://doi.org/10.1007/s42452-022-05270-y with insignificant rate. The slope of Sen's discharge shows an upward trend for the months of July, August, September and October, while the remaining months show a downward trend. In the Kiremt and Bega season the highest and lowest amount of water flows in the watershed, respectively. The correlation of runoff with rainfall and temperature was found to be moderate and very weak in this study, respectively. Overall, the findings of the long−term investigation of hydro−meteorological parameters in the Ajora−Woybo watershed show that water resources vary over time. It is conceivable to assume that a decrease in water supplies in this watershed might related to climate change or changes in land use. They are, however, outside the scope of this study; as a result, further study needed in the future.

Funding
The authors did not receive support from any organization for the submitted work.
Data availability Data are available from the corresponding authors upon reasonable request.

Declarations
Competing of interest The authors have no competing interests to declare that are relevant to the content of this article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.