Performance evaluation of varies climate models using observed and regional climate models for the Katar Watershed, Ethiopia

Climate models are fundamental tools to estimates the reliable future climate change and its effects on the water resources and agriculture in basins. However, all climate models are not equally performed for all areas. Therefore, determining the most appropriate climate models for a specific study area is essential. The focus of this study was to evaluate the performance of the regional climate models with regard to simulating precipitation, and temperatures at Katar watershed. This study examines the performance of fourteen CORDEX-AFRICA-220 Regional Climate Models (RCMs) for the period of 1984–2005 using statistical metrics such as Pearson correlation coefficient (R), mean absolute error (MAE), root mean squared error (RMSE), and bias. The findings indicated that GERICS-MPI was better performed in representing Areta, and Bokoji station, GERICS-IPSL was better representing in Assela, Ketergenet, and Sagure station, CCCma-CanESM2-AFR22, and RCA4-ICHEC performed relatively better in representing the mean annual observed rainfall at the Kulumsa, and Ogolcho station respectively. However, RCA4-CSIRO performed weakly in estimation of annual rainfall at all stations. RCM model such as GERICS-MPI was relatively better than the others in replicating the annual pattern of the maximum temperature at Areta, Bokoji, and Ketergenet stations. Similarly, GER-ICS-IPSL were relatively better in replicating the annual maximum temperature at Assela, and Sagure stations, CCCma-CanESM2-AFR22 at Kulumsa station, and RCA4-ICHEC at Ogolcho station performed well in capturing the observed and simulated annual maximum temperature. Better performance was observed on minimum temperature at CCCma-CanESM2-AFR22 at Areta, Assela, and Ketergenet stations, GERICS-MOHE-AFR-22 at Bokoji station, GERICS-MPI at Kulumsa, and Ogolcho stations, RAC4-NOAA-2G at Sagure stations. However, weak performance was observed RCA4-CSIRO at all stations. RCM models of GERICS-MPI, and CCLM4-NCC-AFR-22 performed better than the other RCM models for correction of annual rainfall in Katar watershed. However, poor performance was observed at RCA4-ICHEC model on Katar watershed. The GERICS-MPI model performed well. However, poor performance was observed at RCA4-ICHEC on maximum temperature, and GERICS-NOAA-2M on minimum temperature in Katar watershed.


Introduction
Climate change increased extremely over Africa show the increased extreme weather conditions (very hot days, heat waves, and high fire-danger days) and in other ways it also the cause of flooding (Sidibé et al 2020).This may be due to partly in response to increased human activities such as deforestation and industrial developments which raise the greenhouse gas effect is the main cause changes in climate variables in a variety of ways, across a variety of periods and in various geographical areas (Shivanna 2022).Such changes may have negative impact on climate-sensitive sectors if appropriate mitigation is not undertaken.Climate change consequences are growing more severe in impoverished nations due to their limited adaptive potential (Asfaw et al 2021).Rising temperatures and changes in rainfall patterns endanger Africa's agriculture, water security, and economic development.Sub-Saharan Africa is the most vulnerable to the escalating effects of climate change (Adenuga et al 2021).
Climate models are fundamental tools used in projecting future climate conditions to develop strategies for adaptation and mitigation methods (Luhunga et al 2018).Some advanced countries, such as the United Kingdom and Canada, established nation-specific techniques for analyzing local-scale climate change consequences, but African countries lack such instruments.To solve this issue, proper regional climate models (RCM) must be selected in order to effectively assess the effects of climate change at the local scale (Adugna and Sime 2021;Amin et al. 2022).General circulation models (GCMs) and regional climate models (RCMs) were developed to forecast future climate conditions and predict possible changes in precipitation and temperature trends.Global Climate Models (GCMs) are widely used for climate projections; however, their low spatial resolution can limit their effectiveness in capturing small-scale climate variations influenced by topography (Yimer et al 2022;Matiu et al 2020;Kaini et al. 2019).However, GCM outputs are not used explicitly to analyse the hydrological impact of climate change.The GCM does not provide useful information for hydrological variables at scales smaller than 200 kms (Islama et al. 2023).Therefore, downscaling GCMs to RCMs was introduced to resolve the discrepancy between the low spatial resolution climate models of GCMs and high spatial resolution regional climate models (RCMs) or hydrological point scales (Chokkavarapu and Mandla 2019).
Regional climate models (RCMs) provide a new opportunity for climate change effect analysis since they have a higher spatial resolution and more reliable results on a regional scale than general circulation models (GCMs).Numerous studies have shown that RCM outputs improve the representation of climate change information at the mesoscale by providing spatially and physically coherent outputs with observations (Min et al. 2018;Turco et al. 2017).However, the original RCM outputs still contain considerable bias, which is inherited from the forcing of GCMs or produced by systematic model errors (Min et al. 2018).According to Dibaba et al. (2019); Mutayoba, and Kashaigili, (2017), not all RCMs are equally important in terms of their effectiveness in capturing a localized study area.The inconsistency in the performance of RCMs when capturing different regions and seasons implies an evaluation of the sensitivity of a region using several available RCMs, and it is necessary to select the most effective RCMs for specific regions (Daniel 2023).
Climate change has a greater impact on developing countries like Ethiopia.Nowadays, significant portion of the country is exposed to drought and desertification.Rift Valley River Basin such as Lake Ziway watershed is one of Ethiopia's most vulnerable river basins to climate changes and water variability (Mulugeta et al. 2020).Katar is the one which is affected by climate change in the Lake-Ziway sub-basin which contribute water for lake Ziway that are influenced by climate change (Abraham et al. 2018) which serves as a crucial water source for water supply for irrigation systems in the region.This lake plays a vital role in supporting agricultural activities and meeting the water needs of the aquatic animals such as Fishery which is majorly produces by local peoples.
Varies studies were conducted at varies places on climate change under RCM.Most of the previous studies on Ethiopia have concentrated on climate-change vulnerability and mitigation measures (Zewdu et al. 2020;Dula et al. 2021), some RCM models performed well in simulating air temperature, but poorly in simulating precipitation (Yonas et al. 2023), some models were biased at higher elevations, but work well for low-elevation regions (Hassler, and Lauer 2021;Getaneh et al. 2024;Gebrekidan et al. 2018;Haopeng et al. 2023), and some studies that have focused on the impact of climate change using single RCMs (Dibaba et al. 2019;Mutayoba, and Kashaigili 2017).A study conducted by Adugna and Sime (2021), Yared et al. (2021), Kebede et al. (2023), (Daniel (2023), Adugna and Sime (2021), and Bekan et al. (2023), were conducted using multiple RCM models.However, the performance of RCMs was varied on watershed, and location to location (Dibaba et al 2019;Mutayoba, and Kashaigili 2017;Abraham et al. 2018;Daniel 2023;Yonas et al. 2023;Meaza et al. 2023;Mengistu et al. 2021).Hence, in order to select the appropriate RCMs for a specific location using multiple RCM models decreases the bias of the models for accurate estimation of climate change.Therefore, this study evaluates the performances of RCM outputs using multiple RCM models for precondition of bias correction (BCM) on the Katar watershed for proper projection of climate variables for impact and adaption studies.

Location and topography of the study area
The Katar catchment is located in Ethiopia's Lake Ziway subbasin.The region was defined by geographical coordinates of 7° 45 'to 8° 30' N and 38° 15 'to 3939° 30' E, as shown in Fig. 1.The Katar River is the largest tributary in the subbasin of Lake Ziway and occupies a catchment area of approximately 3580 square kilometers.The river originates at approximately 4000-4250 m.a.s.l. in the eastern volcanic chains of the Arsi Zone and drains south and west of Lake Ziway from the highlands.The elevation of the catchment ranges from 1620 to 4180 m.a.s.l. with a mean elevation of 2266 m.a.s.l.Only one-third of the whole watershed area has an elevation below 1867 m.a.s.l., and more than 56% of the whole watershed has an elevation greater than 2000 m.a.s.l., which magnifies the upland terrain of most parts of the area, as shown in Fig. 1.

Land use/land cover of the study area
The dominant land use in the Katar watershed is agriculture, as shown in Fig. 2. The basin is intensively cultivated overall, and different crops are grown in the catchments using both rain and irrigation.The LULC map was generated by the Ministry of Water Irrigation and Energy (https:// www.mowe.gov.et/).The major land use types of the subbasin are intensively cultivated.Most of the lower portions are moderately cultivated, shrub land, while Afro alpine and some cultivated land cover are found in the upper portions of the basin.

Model climate data
The station distribution and data availability are expressed in Table 1.The IHMS model typically requires elevation, evapotranspiration, and weather (precipitation and temperature) data to determine the impact on flow at the desired location in the watershed.The size and resolution of spatial data are dependent on the availability of sufficient ground data.Meteorological data from the EMI (Ethiopian Meteorology Institute) were obtained for this particular study.The number of meteorological variables and the duration of the data collected also varied from station to station.
Climate data, including rainfall, maximum and minimum temperatures, and other climate parameters for evapotranspiration estimations are needed as shown in Table 2. ARC-GIS 10.4 Used to obtain the physical parameters and spatial    (1) where; Wi is the weighting of individual rainfall stations, Pi is the rainfall data from the gauging station in mm, dpi is the distance from each station to the required point, m is the exponent and the controlling factor fixed by the user, n is the number of gauging stations, and Pp is the required rainfall data in mm.
Double mass curve analysis was also used to determine the accuracy of the data for certain stations (Guoshuai et al. 2023).The cumulative values of one station are compared for the same duration with the aggregate values of the average values of the other stations.The double mass curve allows homogeneity in time series (especially jumps) to be examined when a shift is to be observed, when the rain gauge is to be used, and so on.Consistent data series may be modified to consistent values by proportionality (Fig. 3).It can be seen from the figure above all the stations are consistent.

Regional climate models (RCM) selection
The performance evaluation of regional climate model outputs, nine dynamically downscaled GCM outputs from the CORDEX program archives were used as shown in Table 3.For this study, meteorological, hydrological, and water resource (land use/land cover digital elevation model and elevation) data were acquired.Three dynamically downscaled GCMs were obtained from CORDEX-22, and the remaining GCMs were obtained from CORDEX-44.These data were generated by a global climate model at horizontal resolutions of 25 and 50 km over the African domain for the period 1980-2005.The selected GCM-RCMs are developed by different research institute and universities over different countries.Five RCM models such as Climate Limited-Area Modeling Community Version 4 (CCLM4), Rossby Center Atmospheric Version 4 (RCA4) (Tumsa 2021), Canadian Centre for Climate Modelling and Analysis (CCCMA) (Yonas 2022, Tumsa 2021), Regional Model (GERICS) (Jacob 2001;Meaza et al. 2023) and African climate model completely new model were used.

Statistical analysis of the models
The statistical parameters included bias; coefficient of determination (R), mean absolute error (MAE) and root mean square error (RMSE) are calculated using Eqs.4-7.R, MAE, RMSE, and Bias are commonly used in multiple studies to evaluate the performance of COR-DEX AFRICA-RCMs (Ayugi et al. 2020;Mendez et al. 2022;Yersaw and Mulusew 2024).The Pearson correlation coefficient (R) values showed the linear relationship between observed and simulated by RCMs.Pearson correlation coefficient values can range from − 1 for a perfect negative correlation to 1 for a perfect positive correlation between the simulated and the observed climate variables (Schober et al. 2018).RMSE measures how accurately climate models simulate climate variables.Smaller values of RMSE close to zero had good model performance and vice versa (Chai, and Draxler 2014).The bias measures the systematic error between the observed and simulated climate variables and zero indicates good performance.
where; R = Pearson correlation coefficient, MAE = Mean absolute error, RMSE = Root mean square error, Obsi and Si are the ith observed and simulated datasets, respectively, Bias = Bias (%), and Obs mean and S mean are the means of the observed and simulated variables, respectively.

Performance of RCMs in simulating annual rainfall in Katar sub stations
The performances of RCMs in simulating monthly rainfall in Katar sub catchments were shown in Table 4 and Fig. 4. The correlation (r) of observed and model simulated mean annual rainfall was strong and very strong (0.68-1.00) in all sub stations at all models except CCLM4-ICHEC which showed weak correlation (R ≤ 0.27).This result coincides with a study conducted in Western Africa where the correlation coefficient value indicated that RCMs simulations of the annual rainfall well matched with the observed mean annual rainfall (Ilori, Balogun 2022;Mengistu et al. 2021).The findings indicated that GERICS-MPI was better performed in representing Areta, and Bokoji station, GERICS-IPSL was better representing in Assela, Ketergenet, and Sagure station, CCCma-CanESM2-AFR22, and RCA4-ICHEC performed relatively better in representing the mean annual observed rainfall at the Kulumsa, and Ogolcho station respectively.The performance of RCM models varied at each station.However, RCA4-CSIRO performed weakly in estimation of annual rainfall at all stations.This might be resulting in the numerical model structure and formulation of RCM.Based on the study by Di Luca et al. (2015), there are two categories of factors influencing the RCMs: those that are related to numerical experimental design and those that are related to the particular climate statistics to be evaluated.Another issue is related to the scale difference between observed data (station data) and gridded data (RCM), the lack of enough spatial and temporal observed data, and the presence of orographic features (Rummukainen et al. 2015).
Generally, most of the RCM models overestimate the monthly rainfall at most Katar sub stations.Similarly, Hernández-Díaz et al. (2013) reported that a dynamically downscaled version of CORDEX-AFRICA showed over estimation of the annual rainfall in the Ethiopian highlands and elevated area of Sudan; Ayugi et al. ( 2020) evaluated RCMs in East Africa and found that most of the models and ensembles overestimated the basin-average annual rainfall amount; Yimer et al. (2022), the RCMs tended to overestimate rain fall conducted an assessment of RCMs in the East Africa region and the Ethiopian highlands area.Additionally, Demessie et al. (2023) specifically examined the Guder sub basin of the upper Blue Nile basin, Kebede et al. (2023), Central Rift Valley Lakes Basin of Ethiopia, Tadase and Mamo, (2024) Arsi Zone Southeastern Ethiopia.Their study revealed that the mean annual rainfall in this sub basin was overestimated by the CORDEX-AFRICA-RCMs.
The MAE and RMSE value indicated the error occurred between the observed and the model simulated mean annual rain fall reached up to 276.1 mm, 368.8 mm (Table 3).The lowest MAE value was 2.3 mm, and RMSE was 3.30 mm, which was observed in GER-ICS-IPSL at Sagure station, while the highest MAE was 16.4 mm, and RMSE was 20.80 mm, which was seen in GERICS-MOHE-AFR-22 at Bokoji station.There is no single RCM or GCM output that has better performance in simulating the rainfall of the study area.The results agree with similar studies in different regions (Stefanidis et al. 2020a, b;Demessie et al. 2023;Gebrekidan et al. 2018).Therefore, intensive studies have been conducted to evaluate the skill of global climate models (GCMs) to simulate climate variables, mainly rainfall and temperature, in different parts, and there is a large difference between intermodal performance and simulated climate variables of the same region or station.Similar studies were conducted in Upper Blue Nile Basin of Jemma sub basin by Gebrekidan et al. (2018), and they recommend that the performance of models better in capturing the annual pattern of the rainfall at one station may performed weakly at the other station and had better correlation with observed than any single model (RCM).

Performance of RCMs in simulating annual maximum temperature in Katar sub stations
The performances of RCM models on annual temperature in Katar sub stations were shown in Table 5 and   Most of the models overestimate the annual maximum temperature at all stations (above 50%).Similarly, almost 50% of the CORDEX-RCMs simulations showed overestimation of maximum and minimum temperature (Tumsa 2022).This finding is in line with studies conducted in the Blue Nile basin (Gebrekidan et al. 2018;Demessie et al. 2023).In comparison of the current result, the overestimation was lower at most RCM models compared with the result of Meaza et al. (2023) conducted at Zarima sub-basin north-western Ethiopia.The over estimation was high in highland areas and low in lowland area of the sub basin.Similarly, the study conducted in Tanzania the bias was relatively high in highly elevated area and low in low elevated area (Luhunga et al. 2016).The reason behind to the result was the RCMs performance varies based on locations and topography in upper Blue Nile basin (Dibaba et al. 2019).

Model performance on annual rainfall on Katar watershed
The model performances of annual rainfall in Katar watershed are shown in Table 7 and Fig. 7.The representative RCMs were selected from fourteen RCM outputs for the catchment based on a statistical comparison of the mean annual rainfall.The catchment mean annual rainfall obtained from ground-observed data was compared with the RCM outputs.The smallest value of P BIAS , RMSE, and MAE of the models shows good performances and also vice versa.Most of the RCM models overestimate the annual rainfall at Katar watershed with    RCM models of GERICS-MPI, and CCLM4-NCC-AFR-22 performed better than the other RCM models for correction of annual rainfall in Katar watershed.However, poor performance was observed at RCA4-ICHEC model on Katar watershed.Most of the models overestimate the observed mean, with the highest value of 3094.74 mm for RAC4-ICHEC, followed by GERICS-NOAA-2M, with P BIAS of 257.93 and 201.52%, respectively.The coefficient of variation was also calculated for the observed and model outputs to compare the interannual variability in rainfall in the catchment.The coefficient of variation observed for the 21 years was 6.86%.Except for RAC4-NOAA-2G and GERICS-ICHEC, all the other models had high CVs compared to the observed values.RAC4-NOAA-2G exhibited the best performance according to the CV curve, followed by GERICS-MPI and CCLM4-NCCR-AFR-22.
The outputs from the RCM model show a large bias, with the highest bias of over 100% for four model outputs (CCLM4-CNRM, RCA4-CSIRO, RCA4-ICHEC, and GERICS-NOAA-2M).This indicates that the systemic error in the model accounts for 100% of the annual rainfall, which suggests that the model was unacceptable.GERICS-MPI has a bias of 7.5%, which indicates that the model well captures catchment-wide rainfall.Most models have a relatively high bias, with values greater than + 50%.This value suggested that the observed rainfall was not well captured by the RCM outputs.The GER-ICS-MPI has the smallest values (94.10 and 81.20 mm/ year) for the RMSE and MAE efficiency measurements, respectively.Accordingly, the largest P BIAS (− 257.93%),RMSE (2320.25 mm/year) and MAE (2336.36mm/year) were observed at RAC4-ICHEC which indicates the poor performance of the model.CCLM4-NCCR-AFR-22 has a bias of less than 10%, and RAC4-NOAA-2G has a bias of 11.55%, which indicates very good and good simulation results, respectively.Additionally, they showed the minimum RMSE and MEA next to GERICS-MPI.Furthermore, both climate models exhibited better performances  in terms of the mean, CV, and correlation coefficient than did the other climate models.
The results of the annual rainfall comparison of the RCM outputs are also similar to the findings of Alemseged, and Tom, (2015) on the evaluation of regional climate model simulations of rainfall over the Upper Blue Nile basin.According to the findings of Wakjira et al. (2021), the GERICS model performed best in terms of bias, CV and RMSE, which is similar to the findings of this study, as GERIC-MPI, which exhibited the best performance in annual rainfall comparisons.The result contradicts the findings proposed by Tumsa, (2021), who indicated that RACMO22T and RCA4 were better, performed at upper awash catchments.The reason behind the result was the performance of Regional climate models was different at different catchments along with the specified locations and topographies (Tumsa 2021).The performance of RCM models varies with catchment.According to Dibaba et al. (2019), Adugna and Sime (2021), and Amin et al. ( 2022), all RCMs are not equal when their performance varies within a localized study area.This result also supported by Dibaba et al. (2019), Mutayoba andKashaigili (2017), Abraham et al. (2018), Daniel (2023), Yonas et al. (2023), andMeaza et al. (2023), the performance of RCMs was varied on watershed, and location to location.

Model performance on annual rainfall anomalies
For the 22-year analysis period, 9 normal, 7 dry and 6 wet years for the Katar catchment were obtained during the analysis period, as shown in Fig. 8.In this analysis, historical climate extremes from 1984 to 2005 were used to assess whether the historical climate year was normal, dry, or wet relative to the reduced regional climate model outputs.A comparison of the three model outputs, which showed better performance in terms of annual rainfall, also showed good results in terms of capturing the occurrence of wet and dry years.The performance of RAC4-NOAA-2G was very good at minimum CV (8.70%), P BIAS (11.55%) followed by CCLM4-NCC-AFR-22 (CV of 12.53%, P BIAS of 7.53%) and GERICS-MPI (CV of 9.10%, P BIAS of 7.50%).However, the RCM models such as Fig. 8 Rainfall anomalies in the Katar catchment CCLM4-MOHC, HadGEM2-ES-AFR-22, and CCLM4-ICHEC performed poorly compared to the other RCM models.The RCM models of CCLM4-MOHC and CCLM4-ICHEC over estimates the annual rainfall with PBIAS of -49.97, -10.39% respectively.This result contradicts with the finding of Adugna and Sime (2021) and Girma et al. (2022), who evaluated the performance of different RCMs for temperature and precipitation in a river basin in Ethiopia and found that RACMO22T and HIRHAM5 performed well.This may be due to variation of topography, climate (Daniel 2023;Yonas 2022;Mengistu et al. 2021).

Model performance for monthly rainfall
The monthly rainfall distributions of the observed rain- Based on the above results for CCLM4-CNRM, GER-ICS-IPSL shows poor performance according to the statistical comparison and thorough visual inspection of the monthly data and from the anomaly index graph.From the graph of the annual monthly rainfall cycle, the above climate models were constructed to capture the wet month and the dry month in the catchment, as Fig. 9 Annual monthly rainfall cycle in the Katar catchment they fully underestimate the wet month rainfall and fully overestimate the dry month rainfall across the catchment.Except for GERICS-MPI, RAC4-NOAA-2G, and CCLM4-NCCR-AFR-22, all the other models highly overestimate the wet month rainfall.This result was in line with Meaza et al. (2023), who suggests Some models overestimated rainfall in dry months while underestimating it in wet months.It was better to use the three models for the Katar catchment that were chosen for any studies pertaining to climate change, employing a multimodal climate for climate impact studies to provide more reliable results, while using the above models for impact and adaptation studies yielded erroneous results.Based on the annual rainfall comparison, annual rainfall anomaly, and monthly rainfall cycle analysis, the three models exhibited the best performances compared with the other models; therefore, it is important to apply a multiclimate model to the Katar catchment.Three climate models (GERICS-MPI, RAC4-NOAA-2G, and CCLM4-NCCR-AFR-22) were selected for Katar as representative climate models.

Model performances on annual mean temperature
A statistical summary of all the models' performances is presented in Table 8 and Figs . 10, 11.The ability of fourteen RCM-CORDEX climate models to simulate the observed minimum and maximum temperatures in the Katar watershed at a seasonal time scale was statistically evaluated using historical time series data from 1984 to 2005.All fifteen RCM-simulated models exhibited relatively better performances in simulating minimum and maximum seasonal temperatures than in simulating seasonal rainfall in the catchment.The GERICS-MPI model performed well.However, poor performance was observed at RCA4-ICHEC on maximum temperature, and GERICS-NOAA-2M on minimum temperature in Katar watershed.
The four RCM models exhibited the largest biases, which were greater than 20%.The three selected models (GERICS-MPI, RAC4-NOAA-2G, and CCLM4-NCCR-AFR22) performed better at correlation coefficient (R 2 ) values, indicating the reproduction of the maximum and minimum seasonal temperatures.The three models also exhibited a small bias of less than ± 5%, which indicated that the model exhibited a perfect fit.According to the annual maximum and minimum temperature comparisons, the bias in the minimum temperature was less than the bias in the maximum temperature for all the downscaled climate model outputs.Only two models (GER-ICS-IPSL and GERICS-NOAA-2M) overestimate the catchment annual maximum temperature, and the rest of the models underestimate it slightly.Many studies have also indicated that climate models may overestimate or underestimate the observed maximum temperature (Demissie et al. 2016).This is likely because rainfall and temperature vary spatially with elevation among the considered meteorological stations.This implies that the best-performing RCMs vary within the study basin; thus, site-specific evaluation of RCMs is indispensable (Mengistu et al. 2021).

Model performances on the mean monthly temperature
A graph of the mean maximum monthly temperature of the observed and simulated RCM models is shown in Fig. 12.Most of these simulations have not been able to capture the maximum and minimum mean monthly temperatures.Six models completely underestimate the mean monthly temperature in the catchment area.RAC4-ICHEC reveals the worst trend in the simulation of the average monthly high temperature.Moreover, the GERICS-MPI, RAC4-NOAA-2G, and CCLM4-NCCR-AFR-22 models exhibited better efficiency; thus, these three models could be used for climate effect studies in the Katar catchment area.For the mean monthly temperature, the GERICS-MPI, RAC4-NOAA-2G, and CCLM4-NCCR-AFR-22 models also exhibited better performances than did the other models.

Conclusion
In this study, the performance of the CORDEX Africa RCM in reproducing rainfall and temperature was evaluated using observed data as a reference for the baseline period from 1980 to 2005.The performance of the RCM was found variable at both spatial and temporal scales.The findings indicated that GERICS-MPI was better performed in representing Areta, and Bokoji station, GER-ICS-IPSL was better representing in Assela, Ketergenet, and Sagure station, CCCma-CanESM2-AFR22, and RCA4-ICHEC performed relatively better in representing the mean annual observed rainfall at the Kulumsa, and Ogolcho station respectively.The performance of RCM models varied at each station.However, RCA4-CSIRO performed weakly in estimation of annual rainfall at all stations.Most of the RCM models overestimated the annual rainfall with bias of − 0.23 mm to − 2.21 mm except RAC4-NOAA-2G, GERICS-MPI, and CCCma-CanESM2-AFR22 at Areta station with bias of 1.55, 0.06, and 0.74 mm respectively.Similarly most of the models overestimated the annual rainfall except CCLM4-ICHEC, CCLM4-MOHC, RAC4-NOAA-2G, CCCma-CanESM2-AFR22, GERICS-MOHE-AFR-22, CCLM4-NCC-AFR-22 with bias of 0.53% to 0.83% at Assela station.Models such as RAC4-NOAA-2G, CCCma-CanESM2, GERICS-MOHE-AFR-22, and CCLM4-NCC-AFR-22 with bias of 0.08 mm to 0.63 mm at Bokoji, RCA4-ICHEC, GERICS-IPSL, and CCCma-CanESM2 with   The study examines the effects of climate change on groundwater and surface water resources, agricultural production, soil loss estimation, reservoir operation, and hydropower production.Furthermore, it promotes the development of a reliable climate service evaluation, enabling well-informed decisions for climate adaptation.Ultimately, these endeavours aim to maximize the benefits of addressing climate change impacts.The variability in RCM performance underscores the importance of bias adjustments in RCM simulation results before utilizing them for climate change impact and adaptation studies.
RCMs have their own inherent biases and uncertainties, which could affect the accuracy of the evaluation.Another limitation is the focus on a specific watershed, which may not be representative of broader regional or global climate trends.Moreover, the study may lack consideration of other important factors such as land use changes, socioeconomic dynamics, and adaptation measures, which could influence the performance of climate models.Finally, the study may not account for future climate scenarios and projections, limiting its applicability for long-term planning and adaptation strategies.

Fig. 1
Fig. 1 Location map of the study area

Fig. 3
Fig. 3 Double mass curve of the rainfall data

Fig. 5 .
Fig.5.The result revealed that the individual RCMs were good in representing the yearly pattern of the observed annual maximum and minimum temperature even if some of the simulations were closer to the observed temperature and others were partly different in replicating the annual maximum and minimum temperature.RCM model such as GERICS-MPI was relatively better than the others in replicating the annual pattern of the maximum temperature at Areta, Bokoji, and Ketergenet stations.Similarly, GERICS-IPSL were relatively better in replicating the annual maximum temperature at Assela, and Sagure stations, CCCma-CanESM2-AFR22 at Kulumsa station, and RCA4-ICHEC at Ogolcho station performed well in capturing the observed and simulated annual maximum temperature.Generally, the performance of GERICS-MPI, and GERICS-IPSL were better in capturing annual temprature in most stations.Girma et al. (2022);Yonas et al. (2023) also supported these findings, showing that GERICS-MPI satisfactory captured the maximum and minimum temperature variability.Most models except RAC4-NOAA-2G, GERICS-IPSL, CCCma-CanESM2-AFR22 at Areta, CCLM4-ICHEC, CCLM4-MOHC, RAC4-NOAA-2G, CCCma-CanESM2-AFR22, GERICS-MOHE-AFR-22, and CCLM4-NCC-AFR-22 at Assela, GERICS-MPI at Bokoji, RCA4-ICHEC, RAC4-NOAA-2G, GERICS-MPI, CCCma-CanESM2, GERICS-MOHE-AFR-22, and CCLM4-NCC-AFR-22 at Ketergenet, RCA4-ICHEC, GERICS-IPSL and CCCma-CanESM2 at Kulumsa, RCA4-ICHEC and CCCma-CanESM2 at Ogolcho, RAC4-NOAA-2G, GERICS-IPSL, and GERICS-MPI at Sagure station overestimate the annual maximum temperature.Most of the models overestimate the annual maximum temperature at all stations (above 50%).Similarly, almost 50% of the CORDEX-RCMs simulations showed overestimation of maximum and minimum temperature(Tumsa 2022).This finding is in line with studies conducted in the Blue Nile basin(Gebrekidan et al. 2018;Demessie et al. 2023).In comparison of the current result, the overestimation was lower at most RCM models compared with the result ofMeaza et al. (2023) conducted at Zarima sub-basin north-western Ethiopia.The over estimation was high in highland areas and low in lowland area of the sub basin.Similarly, the study conducted in Tanzania the bias was relatively high in highly elevated area and low in low elevated area(Luhunga et al. 2016).The reason behind to the result was the RCMs performance varies based on locations and topography in upper Blue Nile basin(Dibaba et al. 2019).

Fig. 4
Fig. 4 Performance of RCMs in simulating annual rainfall in Katar sub stations The statistical evaluation of the annual minimum temperature is expressed in Table6and Fig.6.Better performance was observed CCCma-CanESM2-AFR22 at Areta, Assela, and Ketergenet stations, GERICS-MOHE-AFR-22 at Bokoji station, GERICS-MPI at Kulumsa, and Ogolcho stations, RAC4-NOAA-2G at Sagure stations.However, weak performance was observed RCA4-CSIRO at all stations.The result showed that majority of the simulation overestimated the mean annual minimum temperature similar to the annual rainfall and maximum temperature.The annual minimum temperature P BIAS values range was from ( +) 0.02 to (− ) 2.57 ℃.The lowest and highest bias values are observed at CCCma-CanESM2-AFR22, and RCA4-CSIRO simulation at Ketergenet, and Sagure stations respectively.Whereas the underestimation values were ranged from 0.02 to 0.76 ℃ which was observed at CCCma-CanESM2-AFR22 of Ketergenet station, RCA4-ICHEC of Bokoji station.Most of the RCM model overestimates the annual minimum temperature at Katar watershed sub stations.Similarly, almost 50% of the CORDEX-RCMs simulations showed overestimation of minimum temperature at most models season(Tumsa 2022).The minimum MAE value were 0.06 observed at CCCma-CanESM2-AFR22 of Areta station, 0.97 (℃) at CCCma-CanESM2-AFR22 of Assela station, 0.45 (℃) at GERICS-IPSL, and GERICS-MOHE-AFR-22 of Bokoji station, 0.02 (℃) of CCCma-CanESM2-AFR22 of Ketergenet station, 0.13 (℃) of GERICS-MPI of Kulumsa station, 0.08 (℃) of GERICS-MPI at Ogolcho station, and 0.47 (℃) at RAC4-NOAA-2G of Sagure station.The minimum RMSE were also in the respective MAE values.Whereas, the maximum MAE and RMSE were observed at RCA4-CSIRO at all station of the Katar Sub stations.The performance of the RCM models was also varied in stations.This finding is in line with study conducted in south-western Ethiopia where the CORDEX-AFRICA-RCMs showed low variability of the maximum and minimum temperatures on annual and seasonal time scales (Adugna and Sime 2021).

Fig. 5
Fig. 5 Performance of RCMs in simulating annual maximum temperature in Katar sub stations between − 7.5 and − 257.93%,RMSE ranges between 94.1 and 2320.25 mm, MAE between 81.2 and 2336.36 mm.However, CCLM4-ICHEC and CCLM4-MOHC underestimate the annual mean, with the highest value of − 432.58 mm for CCLM4-MOHC at P BIASs of 10.39 and 49.97%, respectively.The result of this study, overestimation percentage was lower compared with the result ofMeaza et al. (2023) conducted at Zarima sub-basin north-western Ethiopia.The overestimation of rainfall in the highland area may be linked with climate model limitations to resolve topography accurately(Qian, and Zubair 2010;Meaza et al. 2023).

Fig. 6
Fig. 6 Performance of RCMs in simulating yearly minimum temperature in Katar sub stations fall and climate model rainfall over the Katar catchment are shown in Fig. 9. Out of the 14 climate model outputs, two models (HadGM2-ES-AFR-22 and CCLM4-MOHC) failed to completely capture the monthly rainfall trend over the catchment area.CCLM4-MOHC and CCLM4-ICHCE were unable to capture the annual rainfall pattern over the area for a rainy month.Three models, two from CORDEX-AFR-44 (GERICS-MPI and RAC4-NOAA-2G) and one from CORDEX-AFR-22 (CCLM4-NCCR-AFR-22), display strong monthly rainfall patterns but overestimate annual monthly rainfall.The observed data indicate a weak annual rainfall cycle.The monthly precipitation ranges across a small range of values (100 and 135 mm) between February and October.The catchment average is 29-50 mm per month for the majority of the months, with the lowest amount of rainfall occurring in December.Most climate models consistently overestimate the amount of monthly rainfall.

Fig. 10
Fig. 10 Performance of RCMs in simulating annual maximum temperature in Katar watershed

Fig. 11
Fig. 11 Performance of RCMs in simulating annual minimum temperature in Katar watershed

Table 2
Data types and their sources https:// asf.alaska.edu/)Used to obtain the physical parameters and spatial information of the watershed, to generate the climate data from CODEX-Africa to the watershed 12.5 × 12.5 m used for watershed delineationClimateNational Meteorological Agency (NMA) (http:// www.ethio met.gov.et/)Rainfall,temperature, and evapotranspiration (1984-2005)are input data for

Table 4
Mean monthly temperature

Table 5
Statically performance of RCM models on annual Maximum temperature

Table 6
Statically performance of RCM models on annual minimum temperature

Table 7
Annual rainfall comparison of dynamically downscaled rainfall from multiple GCM-RCM Fig. 7 Performance of RCMs in simulating annual rainfall in Katar watershed

Table 8
Model performance for annual minimum and maximum temperatures