Next Article in Journal
Mapping Open Data and Big Data to Address Climate Resilience of Urban Informal Settlements in Sub-Saharan Africa
Previous Article in Journal
Temporal and Spatial Distribution of Lightning Activity over Bulgaria during the Period 2012–2021 Based on ATDnet Lightning Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Understanding Future Climate in the Upper Awash Basin (UASB) with Selected Climate Model Outputs under CMIP6

by
Yonas Abebe Balcha
1,*,
Andreas Malcherek
1 and
Tena Alamirew
2
1
Department of Civil Engineering and Environmental Sciences, Universität der Bundeswehr München, Werner Heisenberg Weg-39, 85577 Neubiberg, Germany
2
Ethiopian Water Resources Institute, Water and Land Resource Center Addis Abab University, Addis Ababa P.O. Box 3880, Ethiopia
*
Author to whom correspondence should be addressed.
Climate 2022, 10(12), 185; https://doi.org/10.3390/cli10120185
Submission received: 7 September 2022 / Revised: 12 November 2022 / Accepted: 15 November 2022 / Published: 22 November 2022

Abstract

:
Climate change makes the climate system of a given region unpredictable and increases the risk of water-related problems. GCMs (global climate models) help in understanding future climate conditions over a given region. In this study, 12 GCMs from the CMIP6 (coupled model intercomparison project six) were evaluated and ranked based on their abilities to describe the historical observed series. The ensemble mean of bias-adjusted best five models of average annual precipitation showed an increment with an uncertainty range of (2.0–11.9) and change in the mean of 6.4% for SSP2-4.5 and (6.1–16.1) 10.6% for SSP5-8.5 in 2040–2069 relative to the historical period. Similarly, for 2070–2099, increments of (2.2–15.0) 7.9% and (11.8–29.4) 19.7% were predicted for the two scenarios, respectively. The average annual maximum temperature series showed increments of (1.3–2.0) 1.6 °C for SSP2-4.5 and (1.7–2.3) 2.0 °C for SSP5-8.5 in 2040–2069. At the same time, increments of (1.7–2.3) 2.0 °C and (2.8–3.2) 3.0 °C were predicted for 2070–2099. Furthermore, it was predicted that the average annual minimum temperature series will have increments of (1.6–2.3) 2.0 °C and (2.2–2.9) 2.5 °C for 2040–2069 and (2.1–2.7) 2.4 °C and (3.7–4.2) 4.0 °C for 2070–2099 for the two scenarios, respectively. An increase in precipitation with increased land degradation in the sub-basin results in a higher risk of flood events in the future. Improved soil and water conservation practices may minimize the adverse impacts of future climate change on the loss of agricultural productivity.

1. Introduction

Securing water that satisfies every human need is the main priority and a challenge faced by many countries throughout the world. Socioeconomic growth, changing consumption patterns, and population growth are believed to be the main driving factors for the ever-increasing water demands [1]. The global freshwater composition that is accessible to humans is less than 1%, which makes it a limited resource [2]. In addition, the naturally uneven distribution of rainfall globally further adds to the complexity of making this resource accessible in adequate amounts to all.
Aside from the natural variation in rainfall distribution, in the past few decades, climate change has further contributed to this variability [3]. Anthropogenically-driven climate variability manifests itself via droughts and floods. Two recent events that support this point are the flood events that occurred in Mecca, Saudi Arabia, in April 2021 [4], and the drought event in southeastern Alaska, USA, in July 2019 [5].
According to a study by the NMA (National Meteorological Agency) on climate change adaptation, the major impacts of climate variability in Ethiopia are food insecurity, outbreaks of water-borne diseases, land degradation, and damage to infrastructure [6]. In particular, the Awash basin is a basin that supports a large population in several major cities and towns, including the capital city Addis Ababa [7]. The basin is highly utilized in the country in terms of water use (mainly for irrigation purposes). Most of the irrigation is practiced in the middle and lower parts of the basin. In addition, intensive rain-fed agriculture dominates the upper basin areas [8,9]. This part of the basin is mostly highland and is densely populated as compared to the middle and lower parts of the basin.
The Awash River is a major river that drains the Awash basin and a large proportion of its average annual flow comes from the upper basin areas. This is due to a relatively high average annual precipitation that falls on these parts of the basin [10]. In recent years, water security issues have become more evident. This is due to the occurrences of more frequent flood and drought events, especially in the middle and lower parts of the basin [11,12].
One of the widely used techniques in studying the impacts of climate change on the water resources of a given place is through the use of outputs from the Coupled Model Intercomparison Project of global climate models (CMIP-GCMs). CMIP-GCMs are mathematical equations that describe the global climate system with three-dimensional grids and are used to simulate the effects of greenhouse gas emissions on climate [13]. Since their introduction, the outputs from CMIP-GCMs have been essential in understanding the Earth’s past, present, and future climate. The use of CMIP-GCMs is often limited to global or regional scales mainly due to their coarse resolutions. However, with the proper application involving a combination of downscaling and bias correction techniques, CMIP-GCMs can be useful in understanding the local climate and play a vital role in the decision-making process for water resource planning and management [14,15].
In recent years, the challenge of predicting future climate has not only focused on climate modeling but also on how to use it. This is due to the sheer number of climate model outputs currently available worldwide. As a result, models often have to be evaluated for their performances in simulating the climate characteristics of the areas in which they are to be applied. The criteria to use for climate model evaluations depend on the climate model selection goals [16].
There are different versions and experiments of climate models generated at different institutions and locations throughout the world. This makes the inter-comparison between models highly cumbersome. This is the reason behind the establishment of the CMIP (climate model intercomparison project), i.e., to place the various model development efforts throughout the world (by different institutions) into one framework [17]. Currently, new sets of models belonging to the sixth CMIP framework have been released. The main differences between models released at each stage of the framework are mainly improvements in the model setups, resolutions, scenarios, and parametrization factors [18].
The CMIP6 archive currently holds model output results from more than 30 climate research centers found on different continents throughout the globe [17]. Each model is different; thus, the outputs generated for the same experiments and scenarios are different [19]. As a result, not all models perform well in certain locations. This is why there is a need to evaluate the capabilities of models and select those that have better performances in describing the local climate [16,20].
In practice, one climate model (or a small ensemble of climate models) is selected for climate change impact studies. In most cases, the selection could be based on a single criterion or a whole set of criteria. Climate models are often selected based on their skills to simulate the present and (near) past climate. This approach is known as the past performance approach. Another approach is called the envelop approach, where models covering a wide range of projections for one or more climate variables of interest are selected from the pool of available models [16,21].
To date, there is no clearly defined methodology to be used for evaluating and selecting a single or a group of climate models [16]. The selection method and process mainly depends on the aim of the study and the variables involved [16]. Most of the works to date have used statistical or data mining techniques. Moreover, there is a difference in the output (single or ensemble model), analysis period (daily, monthly, seasonal, and annual), and level (station, grid, regional, or spatially averaged).
The most common statistical techniques applied for the evaluation of climate models are performance indicators, descriptive statistics, and indices [22,23,24,25]. Data mining methods, such as singular value decomposition (SVD), principal component analysis (PCA), hierarchical clustering, symmetrical uncertainty (SU), and canonical correlation analysis (CCA) are also applied [20,21,26,27,28,29]. The main difference between these techniques is that data mining is a technique that is used to detect patterns and relationships in data, whereas statistical techniques are used to quantify data [30]. Using a combination of these various techniques instead of a single method helps to reduce the uncertainty of the wrong selection and improve it by including various characteristics of the time series. Techniques that are often used and help in incorporating multiple selection methods are skill score (SS) and multi-criteria decision analysis (MCDA) [16,31].
In the past, few attempts aimed to understand the climate of the Awash basin at large (and the upper Awash basin in particular). One of the earliest works on the Awash basin was by Hailemariam [32], where the outputs from three randomly selected CMIP3 GCMs (CCCM, GFD3, and GFDL) together with two scenarios were used in predicting a future runoff condition. After evaluating the potential of GCMs for their annual cycles, seasonal biases, variability, and trends, three CMIP5 GCMs were selected and used to study the impact of climate change by Taye et al. [33]. The impact of climate change on the river basin was studied with help of a few selected models from the CMIP5 archive for different scenarios identified from previous works [34]. Similarly, an ensemble of two GCMs [35] and three selected GCMs [36] from the CMIP5 archive was used to characterize the river flow in the first case and to estimate the river nutrient load in the second one. The ensemble mean of five randomly selected GCMs from the CMIP5 archive was used in characterizing the hydrometeorological situation for the upper Awash basin by Emiru et al. [37]. Most recently, similar research that was published also suggested a group of four climate models, which came from different ensembles of the CMIP6 as an output [38]. Two approaches that combined envelop and past performance approaches were put to use in evaluating the outputs of CMIP6 after [16].
Previous climate change projection efforts in the Awash or upper Awash sub-basin focused on putting into use a few randomly selected climate models without thoroughly evaluating their capability to simulate the climate system over UASB. Therefore, the novelty of this study was through addressing this shortcoming through the execution of more robust evaluation criteria that mainly relied on the application of statistical and data mining techniques. The main outcome of this study will be in identifying a model/s that has an overall good performance over the sub-basin. In addition, a prediction of future climate conditions with the newly selected model will be performed later and compared with similar earlier findings in the sub-basin.

2. Materials and Methods

2.1. Study Area and Data Used

The study focused on the upper Awash sub-basin, which is found in central Ethiopia (Figure 1), with a longitude of 37 ° 57 4 E–39 ° 17 28 E and a latitude of 8 ° 4 52 N–9 ° 19 47 N. The basin has a total drainage area of 12,043 km 2 and contributes the majority of the annual flow to the larger Awash basin. The mean annual precipitation ranges from 861 to 1223 mm at Boneya and Addis Ababa, respectively. The mean annual temperature is in the range of 17 °C at Addis Ababa and 21 °C in Melkasa. The topography of the drainage basin ranges from 3561 m in altitude near ArbGebya and as low as 1547 m at Koka lake.
For the climate model evaluation, a total of 15 rainfall and 10 temperature stations were used. The stations used are shown in Table 1.
Overall, 12 climate models were preliminarily identified from the WCRP CMIP6 archive (https://esgf-node.llnl.gov/search/cmip6/ (accessed on 17 August 2021)) to be used for the study based on model availability for all variables (precipitation, Tmax, and Tmin) for the historical time period. The criteria used to filter the models were daily data, 100 km nominal resolution, and source type: AOGCM and from the r1i1p1f1 variant. The models were accessed from this site on 17 August 2021. The list of models used for analysis is shown in Table 2.

2.2. Methodology

2.2.1. Climate Model Selection

The process begins with extracting the outputs of GCMs for the given station and time period with the help of the CDO (climate data operator) [39]. GCM outputs were extracted to each climate station using the inverse distance weighting (IDW) interpolation technique under CDO. The comparison between the outputs of GCMs and observed data is based on spatially averaged values for the sub-basin. For spatial averaging of station point data, the Thiessen polygon technique was applied. All historical climate series were bias adjusted using the corresponding observed series at each station.
The intentions of the methodologies applied here mainly depend on displaying the important characters of the climate time series that play a critical role in climate impact studies. Thus, based on the intensive literature review conducted on the current methodologies and the objective of this paper, a combination of methodologies was used. The selection or evaluation was based on historical/observed data and was performed at four-time scales or levels—monthly, seasonal (JJAS and MAM), and annual. The two seasons were selected because they are important rainy seasons that determine annual water availability in the sub-basin. It is assumed that the analysis of these time scales will help researchers grasp all statistical characteristics of the time series.
In total, six methodologies were applied, which highlight the various characteristics of the time series. The first technique attempts to fit probability distributions to all climatic variables and checks which climatic models have similar distributions to the observed climate variable. This was performed using the gamlss package [40] under R programming software [41]. The result of the distribution analysis was evaluated between the observed and climate model series using the Akaike information criterion (AIC). Then the climate model or models that have similar distribution(s) is/are assumed to be the model(s) that is/are capable of better describing the observed series.
The second technique observes the trend using the Mann–Kendall (MK) test [42,43]. In detecting the trend, the MK test calculates statistic ‘S’ by ranking the data and calculating the sign as described in Equation (1) [42].
S = i = 1 n 1 j = i + 1 n s g n ( x j x i )
where xj are the data ranked from i = 1, 2, 3, , n−1, and xi are the data ranked from j = i + 1, 2, , n
The sign is the difference between the original data xi, shortened by one data point, and the data itself without the first data point xj.
s g n ( x j x i ) = + 1 i f ( x j x i ) > 0 0 i f ( x j x i ) = 0 1 i f ( x j x i ) < 0
It is observed that when the data point (observation) is more (n ≥ 10), statistic ‘S’ becomes normally distributed with the mean (E(S)) equal to zero and variance calculated as follows:
V a r ( S ) = n ( n 1 ) ( 2 n + 5 ) t = 1 m t 1 ( t 1 1 ) ( 2 t 1 + 5 ) 18
where n is the number of data points and t i are the ties of the sample data series. The test statistic ( Z c ) is calculated as:
Z c = S 1 σ i f S > 0 0 i f S = 0 S + 1 σ i f S < 0
where σ is the standard deviation of the statistic ‘S’. From this test statistic, (Zc) was concluded, a positive value indicates an upward trend, and a negative value is the contrary. According to [43], this approach works well for data where there is no significant correlation at lag 1, if there is a significant correlation at lag 1, the modified MK test is applied. The method calculates the significance of the trend by modifying the variance of the MK test statistic (‘S’) by the ESS (effective sample size). The significance is determined based on the p-value at a significance level of 0.05. Therefore, if the p-value is greater than 0.05, then we accept the null hypothesis that there is no significant trend in the data, and if it is less than or equal, the reverse is true. The results of the MK test can further be supported by Sen’s slope method [44], which calculates the magnitude of the trend (slope ( T i )) as shown in Equation (5) below:
T i = X j X i j i
The median of the N value of T i represents Sen’s slope estimator, which is calculated as follows;
If N appears to be odd:
Q m e d = T ( N + 1 ) 2
If N is even, then:
Q m e d = ( T N 2 + T ( N + 2 ) 2 ) 2
Here, a positive Q m e d value in the time series indicates an upward trend and a negative value the reverse.
In the third approach, four performance metrics were applied to identify which climate model/s better simulated the observed series. The selected measures were coefficient of determination ( R 2 ), root mean squared error (RMSE), mean absolute error (MAE), and BIAS. The performance measures were each applied for all time steps. Moreover, a rank was assigned to each model based on the magnitude of the performance measure. Models that had overall good performances at all four-time steps were identified by summing the ranks at the corresponding time steps and ranking again.
Another technique applied a time series clustering technique to identify which climate model could capture the stochastic process of the observed series much better. For this, the integrated periodogram algorithm developed by [45] was used in the TSclust package [46] under R. This data mining approach attempts to group time series into clusters based on the integrated periodogram ( d I P ) as a distance (dissimilarity) measure. The periodogram technique helps with setting the comparison to be in the frequency domain and this enables characterizing each time series in terms of its underlying stochastic behavior. The dissimilarity measure was calculated using Equation (8) [45],
d I P ( X T , Y T ) = π π | F X T ( λ ) F Y T ( λ ) | d λ
where the normalized cumulative periodograms at each data point j are given by:
F X T ( λ j ) = 1 C X T i = 1 j I X T ( λ i )
F Y T ( λ j ) = 1 C Y T i = 1 j I Y T ( λ i )
The weights used for normalizing the periodograms are as shown below, where m is the number of data points in the original data series:
C X T = i = 1 m I X T ( λ i )
C Y T = i = 1 m I Y T ( λ i )
The periodograms calculated at each time stamp (k) of the observed series are given by:
I X T ( λ k ) = 1 T | t = 1 T X t e i λ k t | 2
I Y T ( λ k ) = 1 T | t = 1 T Y t e i λ k t | 2
Here, λ K = 2 π k T is the frequency component corresponding to the input data sequence k ( k = 1 , 2 , . . . , n ); n = T 1 2 depends on the total data length of the observed series (T). In addition, X t and Y t correspond to each pair of time series for which the dissimilarity measure is to be calculated. Moreover, the term i in Equations (13) and (14) indicates the imaginary term from the Fourier transform of each pair of series.
A pairwise matrix of dissimilarity measure can be produced from the above procedure at which clustering is performed using the agglomerative hierarchical clustering technique under TSclust package in R that applies the hclust() function from the stats package. The classification into a cluster group or merge between clusters is through complete linkage criteria [45].
The fifth technique used is principal component analysis (PCA), which is a dimensionality reduction technique that transforms a large number of correlated variables into much smaller uncorrelated variables called principal components (PCs) [47]. Here, the 11 variables (1 observed + 10 climate models) were displayed into two-dimensional axes called principal components (PCs). The original data points for each variable can be plotted into a two-dimensional space called a score plot. This is done by projecting all data points into those two PCs using the loading vectors obtained from the covariance matrix. The score plot indicates that points closer to the origin are closer to average, points near to each other are similar, and points further outward are outliers. The correlation between each variable can be better viewed using a loading plot. This plot is obtained by plotting the eigenvectors and indicates the contribution of each loading to the PC. The relative length of the vector indicates its contribution to each PC and the angle between the vectors indicates the similarity between the variables. That means that if the angles between the two adjacent loading vectors of variables are smaller, they are more correlated, if orthogonal, they are not related, it they are in the reverse direction, then they are negatively correlated. A combined plot of the score and loading plot is known as a biplot [48].
In the end, the spatial performances of the models were evaluated via correlation analyses and the RMSE as model fit criteria. The two techniques were applied at the station level and later interpreted using a rank-based methodology.

2.2.2. Downscaling, Bias Adjustment, and Future Scenarios

The outputs of GCMs have to be downscaled and bias-adjusted before applying them to real-world situations. This is mainly because the outputs of climate models have biases due to imperfect conceptualization and parametrization, insufficient length of data records, quality of reference data sets, and insufficient spatial resolution [49,50]. Therefore, often the downscaling of GCM outputs is performed by applying either dynamical or statistical downscaling approaches. Dynamic downscaling produces RCMs, which are at finer scales (<50 Km) using GCM boundary conditions and are computationally expensive procedures [51]. More often, statistical downscaling is preferred compared to the former one since there is no drawback (because of the high computational requirements and the ability to bias-adjust). Here, the bias adjustment was performed by creating a statistical link between the large-scale predictors and a finer-scale predictor.
Furthermore, statistical downscaling approaches can be classified into two types—perfect prognosis (PP) and model output statistics (MOS) [52]. PP methods work by creating day-to-day statistical relationships between large-scale predictors and local-scale predictors during the calibration phase. Moreover, the calibrated statistical relationship will be used later for downscaling a predictor variable. Examples of PP could be regression, weather generators, and analog methods. In MOS, there is no day-to-day correspondence, rather, a statistical transfer function is applied during the calibration phase and later used for the downscaling operation [53]. Additive correction/scaling, variance correction, and quantile mapping are some of the techniques under MOS.
Here, both the precipitation and temperature stations were first downscaled to each ground station shown in Table 1 using the IDW technique (15) and later bias-corrected using the quantile matching (QM) approach.
V S t . = ( V G P 1 D 1 + V G P 2 D 2 + V G P 3 D 3 + V G P 4 D 4 ) i = 1 4 1 D i
where V S t . is the variable at the station level, V G P i is the variable at the grid point level, and D i is the distance from the grid point to the specific station.
The QM method uses empirical CDF obtained from the actual observations and requires no assumption of the underlying distribution; it was applied here to the climate model with a top rank [54]. The bias adjustments in this study were performed under R with the help of the qmap package [50] on the daily climate series. The theoretical assumption of QM was to adjust the bias and is shown in Equation (16) as follows:
x f = F f 1 ( F o ( x ) )
where: x f is the value of the bias-adjusted future climate variable; F o ( x ) is the ECDF of the observed variable; F f 1 is the inverse ECDF of the future GCM; and x is the value of the future GCM variable before the bias adjustment.
Three future scenarios (SSP1-2.6, SSP2-4.5, and SSP5-8.5) and three time periods (near-century (2022–2039), mid-century (2040–2069), and end of the century (2070–2099)) were used to understand the future climate of the upper Awash basin. The three scenarios each consider sustainable, middle, and worst future socioeconomic conditions. The worst-case scenario is the upper boundary in terms of the range of scenarios available and it can be considered an update of RCP8.5 of CMIP5. The second one, which is SSP2-4.5, indicates the medium pathway for the future increase in greenhouse gas emissions, which is similar to the RCP4.5 of CMIP5.
In addition, uncertainty limits were added to the calculated expected mean changes (for precipitation ( x ¯ p r e d i c t e d x ¯ o b s e r v e d ) / ( x ¯ o b s e r v e d ) and the Temp ( x ¯ p r e d i c t e d x ¯ o b s e r v e d ) ) of the mean annuals by using the standard deviation across each series (predicted and observed). For example, for the predicted series, the limits were = x ¯ p r e d i c t e d ± x s t a n d a r d d e v i a t i o n , p r e d i c t e d . The same technique was applied when calculating the limits of the observed.

3. Results

3.1. Selection of Climate Model

All the methods used were selected to highlight the very important characteristics of a climate time series, which could play a critical role in addressing the future climate impact in the study area. The evaluation criteria were distribution fitting, trend analysis, performance measures, DTW-based hierarchical clustering, PCA analysis, spatial correlation, and RMSE (root mean squared error). A total of twelve climate models as shown in Table 2 were identified initially based on the common availability of climate models for the three climatic variables (precipitation, maximum temperature, and minimum temperature). The models were evaluated for their statistical characteristics based on seven different evaluation criteria and different time levels (monthly, monthly average, seasonal, and annual). The monthly series was obtained by summing/averaging precipitation/temperature daily series in each month. The monthly average was conducted by taking those monthly values and further averaging them across the total period of analysis for each month.

3.1.1. Identification of Distribution

For the observed and climate model series, the best possible distributions were identified with the help of gamlss package in R [40]. For example, for the precipitation data of the monthly average, JJAS, MAM, and annual series, the observed data were found to have normal (NO), logistic (LO), gamma (GA), and normal (NO) distributions, respectively. The definitions for each distribution type used in fitting the three climate variables are given in Appendix A. The PDF and CDF plots for only the monthly average precipitation series of the observed and climate models are shown in Figure 2 below.
The results of the distribution analysis for all three climatic variables and different time steps of the analysis of both the observed and climate model series are summarized in Table 3. For the precipitation series, it can be deduced that the MRI-ESM2-0 model was capable of having similar distributions with the observed series at all time steps. Moreover, TaiESM1 had similar distributions at the monthly average and annual time steps. NorESM2-MM had similar distributions only at seasonal time steps and ECEARTH3-CC only at MAM.
Regarding the maximum temperature series, none of the climate models showed consistent similarities across different time steps, as shown in Table 3. Here, three models had similar distribution types at monthly average time steps. The GFDL-CM4 model had the same distributions with the observed series twice across those four time steps. Moreover, none of the models had similar distributions at the annual level. Few models shared a similar distribution type with the observed series in the case of minimum temperature. Here, similar distributions were observed only for model CMCC-ESM2 at the monthly average, and MPI-ESM1-2-HR and GFDL-ESM4 at the MAM time steps, respectively. There were no models that were capable of having similar distributions at the annual time step, as shown in Table 3 below.

3.1.2. Trend Analysis

The trend analysis was performed similarly at four different levels (the same as in the distribution analysis). The precipitation series indicates that almost all of the models had similar trends to the observed series across all time steps. The exceptions are ECEARTH3 and ECEARTH3-CC at JJAS and INM-CM4-8 at the monthly average and annual time level, as shown in Table 4. Even though there was no significant trend, the Sen slope values of model TaiESM1 indicated the same trend direction as observed in the series but with different magnitudes.
For the maximum temperature series, ECEARTH3-CC did not have a similar trend to the observed series at any of the time steps. Four models (CMCC-ESM2, NorESM2-MM, TaiESM1, and INM-CM5-0) had similar trends at three-time steps out of four. Moreover, GFDL-ESM4 was the only model that had similar trend characteristics as the observed series at all time steps, as shown in Table 4.
Furthermore, among the ten models, CMCC-ESM2, ECEARTH3-CC, GFDL-CM4, and TaiESM1 had similar trends to the observed minimum temperature series across all time steps, whereas models ECEARTH3, GFDL-ESM4, MRI-ESM2-0, and INM-CM5-0 had the least. Two models, MPI-ESM1-2-HR and NorESM2-MM, showed similarities at all time steps except for the JJAS season, as shown in Table 4.

3.1.3. Performance Metrics

Four different performance indicators were implemented to identify models with good prediction capabilities across all time steps. Table 5 shows the summarized results of each model across those four time steps against four performance measures. The ranking shown in the summary table for each performance measure was obtained by ranking the sum of ranks across each time step.
In the precipitation series, ECEARTH3-CC is the model that had a relatively poor overall performance in predicting the observed series. Moreover, the two other ECEARTH3 families were the top performers. For the maximum temperature series, the NorESM2-MM model performed poorly and TaiESM1 was the top-ranked model in terms of simulating the observed series, as shown in the summary of Table 5. From Table 5, it can also be observed that GFDL-CM4 and ECEARTH3-CC are models that performed well and poorly when simulating the observed minimum temperature series, respectively.

3.1.4. Time Series Clustering

A time-series clustering technique was used in order to identify which climate model was more capable of simulating the observed series in terms of capturing the seasonal variation. This technique uses a distance measure calculated from the cumulative periodograms of the observed and climate models. Using the distance matrix calculated from the integrated periodograms, a complete link agglomerative hierarchical clustering was performed on the monthly precipitation and maximum and minimum temperature series.
From Figure 3, it can be observed that for the precipitation series, the twelve models were clustered into four groups. Moreover, as expected, models that came from the same modeling institution, such as ECEARTH3 and GFDL families, were clustered together. It can be seen that two models, INM-CM4-8 and MPI-ESM1-2-HR, were able to be in the same cluster as the observed series, which means the two are good at simulating seasonal characteristics. Next in line, in terms of having closer seasonal characteristics to the observed series, were INM-CM5-0 and CMCC-ESM2. There was a clear sub-grouping with the rest of the models, which indicated less potential in simulating the seasonal behaviors of the observed series.
Similarly, for the monthly maximum temperature series, a clear clustering between models can be observed. Here, the same models, which are seen in the precipitation series, are the ones that have similar seasonal characteristics as shown in Figure 4.
The clustering of the monthly minimum temperature series is different from the two previous cases, as shown in Figure 5. Here, MRI-ESM2-0 and MPI-ESM1-2-HR performed well in capturing the seasonal variation. Additionally, the dissimilarities between the observed series, the ECEARTH3 family, and NorESM2-MM are relatively higher as compared to the rest of the climate models.

3.1.5. Principal Component Analysis (PCA)

From the PCA analysis of monthly precipitation data, most of the climate models were aligned along the PC1 with the exceptions of TaiESM1, CMCC-ESM2, and MRI-ESM2. Nearly all of the variances of models ECEARTH3-Veg and ECAERTH3-CC are explained by PC1. The models MPI-ESM1-2-HR and ECEARTH3 consist of one group that is closely correlated. As shown in Figure 6, the observed series is strongly correlated to MRI-ESM2.
The ECEARTH3 families and GFDL-CM4 showed good correlations with the observed maximum temperature series, as shown in Figure 7. The MRI-ESM2-0 and CMCC-ESM2 showed the least correlation.
Contrarily to the monthly maximum series, in the monthly minimum series, the ECEARTH3 families showed poor correlation. The most highly correlated model was GFDL-ESM4 followed by the groups of models that included GFDL-CM4, NorESM2-MM, and MRI-ESM2-0, as shown in Figure 8.

3.1.6. Spatial Performances of the Models

The spatial performances of the models were evaluated first by correlating the observed values and climate model outputs at each station for all three climatic variables. For example, Figure 9 shows the box plot of the monthly precipitation correlation coefficient values for each climate model at all stations. From this figure, we can infer that models MPI-ESM1-2-HR, ECEARTH3, and MRI-ESM2 had relatively better correlation to the observed monthly precipitations across the basin. Models CMCC-ESM2 and TaiESM1 were the poorly performing ones.
Similarly, the correlation box plots for the maximum temperature series (Figure 10) show that two models, MPI-ESM1-2-HR and ECEARTH3, were better correlated; MRI-ESM2 performed the poorest.
The correlation coefficient values for the minimum temperature (Figure 11) are not consistent with the precipitation and maximum temperature series. Here, models belonging to the GFDL family showed relatively good performances against the others, while the ECEARTH families had the lowest correlations.
In addition to the analysis of the spatial distribution of the correlation coefficients, RMSE was used to observe which climate model had the better model fit to the observed series. As can be observed from Figure 12, models GFDL-CM4, MRI-ESM2, and GFDL-ESM4 had less spatial RMSE values for the monthly precipitation when compared to the rest. INM-CM4-8 and TaiESM1 had the least spatial performances. In all models, there were fewer errors in simulating the rainfall of stations located at higher elevations.

3.2. Model Ranking and Interpretation

In the interpretation and summary of the results, a simple rank-based system was used. A higher rank (lowest in terms of magnitude) was given to a specific climate model for the desired effect. For example, in the interpretation of the results of the distribution fitting—if a specific model was fit with the same model as the observed series, then a value of one (unit) was assigned, and if otherwise, then zero. These values would be summed across the four analysis periods for each model and given a rank. The ranking was from one to twelve similar to the total number of climate models used in the study. A model had a higher rank in this case if the sum of ranks was higher.
A similar approach is also used in the interpretation of the trend analysis. That means, if a similar trend is obtained for a specific model, then a value of one is assigned, and zero if the contrary. The rest of the ranking procedure is the same as the one used in the interpretation of the distribution analysis. For the case of the performance measure, a higher rank was given to a specific model that scored well according to the criteria of each measure. Since four performance measures were used for each model, the final was based on ranking the sum of ranks.
To interpret the cluster analysis, the result of the distance (dissimilarity) matrix was used. A model with a minimum distance from the observed series was ranked higher than the one with a large distance. To summarize the PCA analysis, the angle between the loading vectors of the observed and climate models was used in the ranking criteria. Here, similarly among the ten models, a model with fewer angles was much closer to the observed series, so a higher-ranked was assigned to it. Finally, the areal-averaged values of two spatial measures (correlation and RMSE) were given a rank.
From the result of the precipitation time-series analysis across each method, the TaiESM1 model performed poorly, as shown in Table 6. The top three models that showed overall good performances were MRI-ESM2, GFDL-ESM4, and MPI-ESM1-2-HR.
There is no common model that has a good performance across the three climate variables as seen in Table 6. The poorest performance for the maximum temperature series was observed in MRI-ESM2. In the case of the minimum temperature series, the ECEARTH3 family performed poorly; contrarily, the GFDL family was one of the top performers, as shown in Table 6.
Finally, the ranking methodology was applied for the identification of possible models with overall good performances for all three climatic variables of concern, as shown in Table 7. It can be observed that the MPI-ESM1-2-HR model is one of the models with the highest overall ranking when compared across the three climatic variables.

3.3. Future Climate Projections

In order to comment on the future climate conditions in the study area, the QM technique was applied as a tool to bias-adjust the top five climate models (MPI-ESM1-2-HR, GFDL-ESM4, INM-CM5-0, MRI-ESM2, and ECEARTH3). Future climate projections were done using an ensemble average of those five models. Because of the absence of the required scenarios, the GFDL-CM4 model was not used in future climate projection. In terms of the precipitation variables, for three time periods of analysis (near, middle, and end of the century), a general increase in the precipitation amount was observed (Table 8) as compared to the base period (1980–2009). The magnitude of increment in the mean annual areal precipitation for the near, middle, and end of the century were 1.3%, 10.6%, and 19.7%, respectively, for the worst-case socioeconomic scenario. Moreover, for the middle socioeconomic scenario, increments of 2.7%, 6.4%, and 7.9% were seen for the near, middle, and end of the century, respectively. Here, except for the sustainable scenario, a higher magnitude of increments was expected by the end of this century. When the seasonal patterns of change were investigated, the magnitude of the increment for the spring rainfall was slightly higher than the summer, as shown in Table 9. Here, we can see a predicted decrement in the summer rainfall, especially for the near-time period.
The mean annual areal maximum temperature in the middle scenario showed increments of 1.1 °C, 1.6 °C, and 2.0 °C for the near, middle, and end of the century, respectively. With the same pattern, increases of 1.2 °C, 2.0 °C, and 3.0 °C were seen for the worst case, respectively. Similarly, for the minimum temperature series in the middle scenario, increases of 1.3 °C, 2.0 °C, and 2.4 °C were observed for near, middle, and end of the century, respectively. The worst-case scenario showed increments of 1.4 °C, 2.5 °C, and 4.0 °C for the near, middle, and end of the century, respectively. It can be concluded that—similar to the precipitation variable—higher magnitudes of temperature increments are observed at the end of the century. Moreover, the increment amounts are relatively higher for minimum temperature than the maximum temperature.
The results of the future climate analysis discussed in the previous two paragraphs indicate the areal-averaged outputs only. In order to understand how the changing pattern may look spatially, four representative stations were selected across the basin. The relative geographic locations (north, east, south, west) in the sub-basin were the main criteria for selections; the intention was to observe if there was any spatial pattern in future climate across the sub-basin.
Analysis of future precipitation (Table 10) shows that there seems to be a relatively higher precipitation mean change and larger uncertainty in the mean change prediction for stations Ejere and Hombole. In the case of future temperature changes, stations Deberzeit and Bui had higher changes when compared to the other two stations. Moreover, stations Bui and Tulubolo had larger uncertainties in predicting the mean change for both maximum and minimum temperature series. There was no clear spatial pattern to observe in the change of future precipitation across the sub-basin. However, the future changes in temperature values show that stations located at lower elevations had higher relative changes.

4. Discussion

Models that are incorporated in the new CMIP6 archive showed comparable or improved capabilities in terms of simulating the climate of the globe as compared to CMIP5 [55]. This is because the models in the CMIP6 archive showed improved model resolution, setups, scenarios, and parametrization compared to all previous CMIP versions [18]. This study mainly focused on the application of statistical and data mining techniques to identify possible climate models from the new CMIP6 archive, which can simulate the climate system of UASB. This was made possible after evaluating all models against the observed historical climate series using the mentioned evaluation criteria. All comparisons were performed on spatially averaged data for the 1980–2009 period of analysis.
No specific model had a consistent performance across the evaluation criteria and climatic variables as shown in Table 6. Model MRI-ESM2-0 represented the distributional properties of the precipitation series. GFDL-CM4 and MPI-ESM1-2-HR performed well for the minimum and maximum temperature series, respectively. In terms of capturing the trend behavior the GFDL-ESM4 model performed equally well for the precipitation and maximum temperature series; models GFDL-CM4, TaiESM1, CMCC-ESM2, and ECEARTH3-CC had equal potential for the minimum temperature series. The outputs of the model fit performance measures indicated that no specific models had similar performances across the three climate variables. The seasonal patterns were better captured by MPI-ESM1-HR for all three climate variables. Not all three climate variables were correlated to the same climate models in the PCA analysis. MRI-ESM2, EC-EARTH3-CC, and GFDL-ESM4 had better correlations with the precipitation, maximum temperature, and minimum temperature series, respectively.
The results of the two spatial evaluation criteria indicate that MRI-ESM2 had better spatial correlation but GFDL-CM4 was a good model fit for the precipitation series. MPI-ESM1-2-HR and ECEARTH3 ranked the highest for the maximum temperature serious based on those criteria, respectively. In addition, the GFDL-CM4 model ranked number one in the minimum temperature series for both criteria.
According to a recent study by [56], the summer rainfall over Ethiopia is influenced by SST conditions over the Gulf of Guinea and the southern Pacific Ocean. The spring rainfall was influenced by SST over the North Atlantic Ocean. A study by [57] indicated that the summer rainfall season (JAS) climate system over the UASB is highly influenced by the equatorial pacific ocean temperature. ENSO and ITCZ that arise from the ocean temperature variations over the equatorial pacific and low-pressure zone near the equator are also two important climate conditions that shape the climate over UASB.
Most of the climate models evaluated in this study have different performances in terms of simulating the global climate process. A recent work presented by [55] showed that models EC-EARTH3, MPI-ESM1-2-HR, and TaiESM1 have good capabilities in terms of capturing ENSO teleconnections. The MPI-ESM1-2-HR model is also capable of producing teleconnections over the Indian and North Atlantic oceans reasonably well [58]. As shown in the results section, the models are among the top-ranked models in this study.
Another important reason why the top-ranked MPI-ESM1-2-HR model performed well may be due to the fact that it has higher atmospheric and ocean model resolutions relative to the ten models. This model has 0.94 ° atmospheric resolution with 95 verticals and the model top at 0.01 hPa. Moreover, the ocean model is having 0.4 ° horizontal resolution with 40 levels [58,59]. The poorest performing model, which was CMCC-ESM2, has an atmospheric resolution of 0.9 ° × 1.25 ° with 30 verticals and a model top at 2 hPa. Moreover, the ocean model has a resolution of 0.33 ° [60]. Since the topography of the study area is complex, it is expected that the model resolution plays a vital role in capturing all important climate characteristics of the study area.
Past climate change studies on the basin that applied CMIP5 outputs have all predicted an increase in temperatures for all time periods and scenarios [33,34,35]. However, there is no clear agreement on future precipitation changes. The works by Taye et al. [33] and Tadese et al. [34] have indicated an increase in the summer (JJAS) and a decrease in the spring (MAM) rainfalls for both the 2050s and 2070s and two scenarios (RCP4.5 and RCP8.5). Daba and You [35] predicted an increase in summer rainfall for the 2020s and a decrease in both the summer and spring rainfalls for the 2050s and 2080s for those two same scenarios. However, the previous AR5 for CMIP5 and the new AR6 for the CMIP6 IPCC reports predict that there will be increases in the annual mean precipitation and maximum temperature values for both the mid and end of the century over East Africa [61,62]. So, it can be seen that the findings from these previous studies are in disagreement with the IPCC report. These fallacies may be due to the poor performances of these previous climate models in simulating the climate conditions over East Africa or improper model selections. In this study, it was predicted that there are increases in both precipitation and temperature Table 8 presents different magnitudes for three scenarios, especially for the mid- and far-time periods. Contrary to CMIP5 outputs, it was found that there is a slightly higher increment in the spring rainfall than in the summer. In addition, the uncertainty bands of future precipitation changes indicate a relatively larger uncertainty toward the end-of-century, while for the temperature series, the reverse is true. This might be due to the characteristics of precipitation prediction, which makes them more likely to be influenced by local factors than temperature. It can be observed that all of the findings of this study on future climate projections align with the new AR6 IPCC report on future climate conditions over Northeast Africa [62].

5. Conclusions

The outputs of the selection process have indicated that it is difficult to conclude the performance of a single model since the performance of each model is different for each evaluation criterion. However, the MPI-ESM1-2-HR model showed an overall good performance when evaluated for all time steps and evaluation criteria.
When observing the outputs of each evaluation criterion across the three climatic variables, the GFDL-CM4 model seems to be the best model across the three measures (PDF, PCA, and spatial). The second evaluation criteria for trend behavior shows that the CMCC-ESM2 model could capture the trend well. Outputs of model fit performance measures ( R 2 , RMSE, MAE, and BIAS) indicate that INM-CM4-8 outperformed the others. MPI-ESM1-2-HR is the best model in terms of capturing the stochastic behavior (underlying periodic patterns) of the observed series.
Downscaling and bias adjustment to the scenarios obtained from the top five high-ranking models (ECAERTH3, GFDL-ESM4, MPI-ESM1-2-HR, MRI-ESM2, and INM-CM5-0) were performed using the IDW and QM approaches for all three climate variables (precipitation, maximum temperature, and minimum temperature). Based on the results of downscaling and bias adjustment, predictions of future climate conditions for three socioeconomic scenarios (sustainable, middle, and worst) across three time periods (near, middle, and end of the century) were conducted. The results of future climate prediction showed relative increments compared to the base period (1980–2009) for all three climate variables. The magnitude of the increment was relatively higher for the end than near the middle of this century. In addition, the analysis of seasonal projection indicated that the magnitude of the increment for spring rainfall was higher than for the summer. It can be observed that outputs of this study are in agreement with the new AR6 IPCC report about expected future climate changes in Northeast Africa.
The intensifying agricultural activity and urbanization in the study area [7] with the projected increase in precipitation amount could lead to higher risks of flooding over UASB in the future.
The results of this study could be verified through the inclusion of other ensemble variants in addition to r1i1p1f1.

Author Contributions

Conceptualization, Y.A.B. and T.A.; methodology, A.M.; software, Y.A.B.; validation, A.M., T.A. and Y.A.B.; formal analysis, Y.A.B.; investigation, Y.A.B.; resources, Y.A.B.; data curation, Y.A.B.; writing—original draft preparation, Y.A.B.; writing—review and editing, T.A.; visualization, Y.A.B.; supervision, A.M.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the DAAD Home Grown Program 2019 and Universität der Bundeswehr München.

Data Availability Statement

All the data and analysis done is available upon request.

Acknowledgments

We would like to thank Ethiopian Meteorological Agency and WCRP (World Climate Research Program) for making the data available for this research. We also thank the editors, reviewers, and colleagues for their constructive comments when developing this manuscript.

Conflicts of Interest

The authors declare no conflict of interest in the development of this manuscript.

Appendix A

A description of the distribution types under the gamlss package [63] used in Section 3.1.1 is listed as follows:
Normal (or Gaussian) distribution (NO( μ , σ )): the normal distribution is parameterized; its PDF is defined by
f Y ( y | μ , σ ) = 1 2 π σ exp ( y μ ) 2 2 σ 2
For < y< , where < μ < and σ > 0; where μ is the mean and σ is the standard deviation of Y.
Logistic distribution (LO( μ , σ )): The Logistic distribution is parameterized and its PDF is defined by
f Y ( y | μ , σ ) = 1 σ e x p y μ σ 1 + e x p y μ σ 2
For < y < , where < μ < and σ > 0; where μ is the mean and σ is the standard deviation of Y.
Gumbel distribution (GU( μ , σ )): the PDF is defined by
f Y ( y | μ , σ ) = 1 σ e x p y μ σ e x p y μ σ
For < y < , where < μ < and σ > 0.
Reverse Gumbel distribution (RG( μ , σ )): PDF is defined by
f Y ( y | μ , σ ) = 1 σ e x p y μ σ e x p ( y μ ) σ
For < y < , where < μ < and σ > 0.
Power exponential distributions (PE):
PE( μ , σ , ν ) PDF is defined by
f Y ( y | μ , σ , ν ) = ν e x p | z c | ν 2 c σ Γ 1 ν
PE2( μ , σ , ν ) PDF is defined by:
f Y ( y | μ , σ , ν ) = ν e x p | z | ν 2 c σ Γ 1 ν
For < y < , where < μ < and ν > 0. Moreover, c 2 = Γ ( 1 / ν ) Γ ( 3 / ν ) 1 .
Normal exponential t-distribution (NET( μ , σ , ν , τ )) is a four-parameter distribution and PDF is defined by
f Y ( y | μ , σ , ν , τ ) = c σ e x p z 2 2 , w h e n | z | ν e x p ν | z | + ν 2 2 , w h e n ν < | z | τ e x p ν τ log | z | τ ν τ + ν 2 2 , w h e n | z | > τ
For < y < , where < μ < , ν > 1, and τ > ν . Moreover, c = ( c 1 + c 2 + c 3 ) 1 , where, C 1 = 2 π 1 2 Φ ( ν ) , C 2 = 2 ν e x p ν 2 2 , and C 3 = 2 ( ν τ 1 ) ν e x p ν τ + ν 2 2 , where Φ (.) is the cumulative distribution function of the standard normal distribution.
Sinh-Arcsinh (SHASH( μ , σ , ν , τ )): its PDF is defined by
f Y ( y | μ , σ , ν , τ ) = c 2 π σ ( 1 + r 2 ) 1 / 2 e z 2 / 2
where z = 1 2 e x p τ s i n h 1 ( r ) e x p ν s i n h 1 ( r ) , and c = 1 2 τ e x p τ s i n h 1 ( r ) + ν e x p ν s i n h 1 ( r ) , and r = ( y μ ) / σ for < y < , where < μ < and σ > 0, ν > 0, and τ > 0.
Skew exponential power type I distribution (SEP1( μ , σ , ν , τ )): its PDF is defined by
f Y ( y | μ , σ , ν , τ ) = 2 σ f Z 1 ( z ) F Z 1 ( ν z )
For < y < , where < μ < and σ > 0, < ν < , and τ > 0, and where z = ( y μ ) / σ and f Z 1 and F Z 1 are the PDF and CDF of Z 1 P E 2 ( 0 , τ 1 / τ , τ ) a power exponential type 2 distribution with f Z 1 = σ 1 e x p [ | z | τ / τ ] , where α = 2 τ ( 1 / τ ) 1 Γ ( 1 / τ ) .
Skew exponential power type II distribution (SEP2( μ , σ , ν , τ )): its PDF is defined by
f Y ( y | μ , σ , ν , τ ) = 2 σ f Z 1 ( z ) Φ ( ω )
For < y < , where < μ < and σ > 0, < ν < , and τ > 0, and where z = ( y μ ) / σ and ω = s i g n ( z ) | z | τ / 2 and f Z 1 is the PDF of Z 1 P E 2 ( 0 , τ 1 / τ , τ ) and Φ ( ω ) is the CDF of the standard normal variable evaluated at ω .
Skew exponential power type III distribution (SEP3( μ , σ , ν , τ )): its PDF is defined by
f Y ( y | μ , σ , ν , τ ) = c σ e x p [ 1 2 | ν z | τ ] I ( y < μ ) + e x p [ 1 2 | z ν | τ ] I ( y μ )
For < y < , where < μ < , σ > 0, ν > 0, and τ > 0 and where z = ( y μ ) / σ and c = ν τ / [ ( 1 + ν 2 ) 2 1 / τ Γ ( 1 τ ) ] .
Gamma distribution (GA( μ , σ )): its PDF is defined by
f Y ( y | μ , σ ) = 1 2 π σ 2 1 y e x p [ l o g ( y ) μ ] 2 2 σ 2
For y > 0, where μ > 0 and σ > 0.
Inverse Gaussian distribution (IG( μ , σ )): its PDF is defined by
f Y ( y | μ , σ ) = 1 2 π σ 2 y 3 e x p 1 2 μ 2 σ 2 y ( y μ ) 2
For y > 0, where μ > 0 and σ > 0.
Weibull distributions (WEI1, WEI2, WEI3): the parametrization and PDF are defined by
The first parametrization (WEI1 (WEI( μ , ν ))) is as follows,
f Y ( y | μ , σ ) = σ y σ 1 μ σ e x p y μ σ
For y > 0, where μ > 0 and σ > 0.
Second parametrization (WEI2)
f Y ( y | μ , σ ) = σ μ y σ 1 e μ y σ
For y > 0, where μ > 0 and σ > 0.
Second parametrization (WEI3):
f Y ( y | μ , σ ) = σ β y β σ 1 e x p y β σ
For y > 0, where μ > 0, σ > 0 and where β = μ / Γ ( 1 σ + 1 ) .
Box–Cox power exponential distribution (BCPE( μ , σ , ν , τ )): its parametrization and PDF are defined by
f Y ( y | μ , σ , ν , τ ) = y ν 1 f T ( z ) μ ν σ F T ( 1 σ | ν | )
For y > 0, μ > 0, σ > 0, τ > 0, < ν < , z = ( y μ ) / σ , and where f T and F T are the PDF and CDF of a random variable T.
Generalized t-distribution (GT( μ , σ , ν , τ )): its parametrization and PDF are defined by
f Y ( y | μ , σ , ν , τ ) = τ 2 σ ν 1 / T B ( 1 / τ , ν ) 1 + | z | τ / ν ν + ( 1 / τ ) 1
For < y < , where < μ < , σ > 0, ν > 0, and τ > 0 and where z = ( y μ ) / σ .

References

  1. Water, C. Leaving No One Behind. In The United Nations World Water Development Report Paris; UNESCO: Paris, France, 2019. [Google Scholar]
  2. Gleick, P.H. Water in Crisis; Oxford University Press: New York, NY, USA, 1993; Volume 100. [Google Scholar]
  3. Thakural, L.; Kumar, S.; Jain, S.K.; Ahmad, T. The impact of climate change on rainfall variability: A study in central himalayas. In Climate Change Impacts; Springer: Singapore, 2018; pp. 181–192. [Google Scholar]
  4. Shaaban, F.; Othman, A.; Habeebullah, T.M.; El-Saoud, W.A. An integrated GPR and geoinformatics approach for assessing potential risks of flash floods on high-voltage towers, Makkah, Saudi Arabia. Environ. Earth Sci. 2021, 80, 1–15. [Google Scholar] [CrossRef]
  5. Claire, S.; Denise, G. DroughtScape. 2019, pp. 3–6. Available online: https://digitalcommons.unl.edu/cgi/viewcontent.cgi?article=1052&context=droughtscape (accessed on 14 November 2022).
  6. NMA. Climate Change National Adaptation Programme of Action (NAPA) of Ethiopia; National Meteorological Services Agency (NMA), Ministry of Water Resources, Federal Democratic Republic of Ethiopia: Addis Ababa, Ethiopia, 2007. [Google Scholar]
  7. Shawul, A.A.; Chakma, S. Spatiotemporal detection of land use/land cover change in the large basin using integrated approaches of remote sensing and GIS in the Upper Awash basin, Ethiopia. Environ. Earth Sci. 2019, 78, 1–13. [Google Scholar] [CrossRef]
  8. Nanesa, K. Awash River’s the Ongoing Irrigation Practices, Future Projects and its Impacts on the Environment of Awash River Basin. Irrig. Drain. Syst. Eng. 2021, 10. [Google Scholar] [CrossRef]
  9. Aregahegn, Z.; Zerihun, M. Study on Irrigation Water Quality in the Rift Valley Areas of Awash River Basin, Ethiopia. Appl. Environ. Soil Sci. 2021, 2021, 8844745. [Google Scholar] [CrossRef]
  10. Mersha, A.N.; de Fraiture, C.; Masih, I.; Alamirew, T. Dilemmas of integrated water resources management implementation in the Awash River Basin, Ethiopia: Irrigation development versus environmental flows. Water Environ. J. 2021, 35, 402–416. [Google Scholar] [CrossRef]
  11. Edossa, D.C.; Babel, M.S.; Gupta, A.D. Drought analysis in the Awash river basin, Ethiopia. Water Resour. Manag. 2010, 24, 1441–1460. [Google Scholar] [CrossRef]
  12. Shenduli, P.R.; Van Andel, S.J.; Mamo, S.; Masih, I. Improving hydrological prediction with global datasets: Experiences with Brahmaputra, Upper Awash and Kaap catchments. In Proceedings of the E-Proceedings of the 37th IAHR World Congress, Kuala Lumpur, Malaysia, 13–18 August 2017. [Google Scholar]
  13. Edwards, P.N. History of climate modeling. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 128–139. [Google Scholar] [CrossRef] [Green Version]
  14. Bhuvandas, N.; Timbadiya, P.V.; Patel, P.L.; Porey, P.D. Review of downscaling methods in climate change and their role in hydrological studies. Int. J. Environ. Ecol. Geol. Mar. Eng. 2014, 8, 713–718. [Google Scholar]
  15. Soriano, E.; Mediero, L.; Garijo, C. Selection of bias correction methods to assess the impact of climate change on flood frequency curves. Water 2019, 11, 2266. [Google Scholar] [CrossRef] [Green Version]
  16. Lutz, A.F.; ter Maat, H.W.; Biemans, H.; Shrestha, A.B.; Wester, P.; Immerzeel, W.W. Selecting representative climate models for climate change impact studies: An advanced envelope-based selection approach. Int. J. Climatol. 2016, 36, 3988–4005. [Google Scholar] [CrossRef] [Green Version]
  17. Brief, C. How Do Climate Models Work? 2018. Available online: https://www.carbonbrief.org/qa-how-do-climate-models-work/ (accessed on 17 August 2021).
  18. Ayugi, B.; Zhihong, J.; Zhu, H.; Ngoma, H.; Babaousmail, H.; Rizwan, K.; Dike, V. Comparison of CMIP6 and CMIP5 models in simulating mean and extreme precipitation over East Africa. Int. J. Climatol. 2021, 41, 6474–6496. [Google Scholar] [CrossRef]
  19. Wilby, R. Evaluating climate model outputs for hydrological applications. Hydrol. Sci. J. J. Des. Sci. Hydrol. 2010, 55, 1090–1093. [Google Scholar] [CrossRef]
  20. Khayyun, T.S.; Alwan, I.A.; Hayder, A.M. Selection of suitable precipitation CMIP-5 sets of GCMs for Iraq using a symmetrical uncertainty filter. In Proceedings of the IOP Conference Series: Materials Science and Engineering; IOP Publishing: Bristol, UK, 2020; Volume 671, p. 012013. [Google Scholar]
  21. Homsi, R.; Shiru, M.S.; Shahid, S.; Ismail, T.; Harun, S.B.; Al-Ansari, N.; Chau, K.W.; Yaseen, Z.M. Precipitation projection using a CMIP5 GCM ensemble model: A regional investigation of Syria. Eng. Appl. Comput. Fluid Mech. 2020, 14, 90–106. [Google Scholar] [CrossRef]
  22. Samadi, S.Z.; Sagareswar, G.; Tajiki, M. Comparison of general circulation models: Methodology for selecting the best GCM in Kermanshah Synoptic Station, Iran. Int. J. Glob. Warm. 2010, 2, 347–365. [Google Scholar] [CrossRef]
  23. Pitman, A.; Perkins, S. Reducing uncertainty in selecting climate models for hydrological impact assessments. IAHS Publ. 2007, 313, 3. [Google Scholar]
  24. Alemseged, T.H.; Tom, R. Evaluation of regional climate model simulations of rainfall over the Upper Blue Nile basin. Atmos. Res. 2015, 161, 57–64. [Google Scholar] [CrossRef]
  25. Ongoma, V.; Chen, H.; Gao, C. Evaluation of CMIP5 twentieth century rainfall simulation over the equatorial East Africa. Theor. Appl. Climatol. 2019, 135, 893–910. [Google Scholar] [CrossRef]
  26. Nashwan, M.S.; Shahid, S. Symmetrical uncertainty and random forest for the evaluation of gridded precipitation and temperature data. Atmos. Res. 2019, 230, 104632. [Google Scholar] [CrossRef]
  27. Ahmadalipour, A.; Rana, A.; Moradkhani, H.; Sharma, A. Multi-criteria evaluation of CMIP5 GCMs for climate change impact analysis. Theor. Appl. Climatol. 2017, 128, 71–87. [Google Scholar] [CrossRef] [Green Version]
  28. Rana, A.; Madan, S.; Bengtsson, L. Performance evaluation of regional climate models (RCMs) in determining precipitation characteristics for Gothenburg, Sweden. Hydrol. Res. 2014, 45, 703–714. [Google Scholar] [CrossRef]
  29. Wilcke, R.A.; Bärring, L. Selecting regional climate scenarios for impact modelling studies. Environ. Model. Softw. 2016, 78, 191–201. [Google Scholar] [CrossRef] [Green Version]
  30. Friedman, J.H. Data Mining and Statistics: What’s the connection? Comput. Sci. Stat. 1998, 29, 3–9. [Google Scholar]
  31. Srinivasa Raju, K.; Nagesh Kumar, D. Ranking general circulation models for India using TOPSIS. J. Water Clim. Chang. 2015, 6, 288–299. [Google Scholar] [CrossRef]
  32. Hailemariam, K. Impact of climate change on the water resources of Awash River Basin, Ethiopia. Clim. Res. 1999, 12, 91–96. [Google Scholar] [CrossRef] [Green Version]
  33. Taye, M.T.; Dyer, E.; Hirpa, F.A.; Charles, K. Climate change impact on water resources in the Awash basin, Ethiopia. Water 2018, 10, 1560. [Google Scholar] [CrossRef] [Green Version]
  34. Tadese, M.T.; Kumar, L.; Koech, R. Climate change projections in the Awash River Basin of Ethiopia using Global and Regional climate models. Int. J. Climatol. 2020, 40, 3649–3666. [Google Scholar] [CrossRef]
  35. Daba, M.H.; You, S. Assessment of climate change impacts on river flow regimes in the upstream of Awash Basin, Ethiopia: Based on IPCC fifth assessment report (AR5) climate change scenarios. Hydrology 2020, 7, 98. [Google Scholar] [CrossRef]
  36. Bussi, G.; Whitehead, P.G.; Jin, L.; Taye, M.T.; Dyer, E.; Hirpa, F.A.; Yimer, Y.A.; Charles, K.J. Impacts of Climate Change and Population Growth on River Nutrient Loads in a Data Scarce Region: The Upper Awash River (Ethiopia). Sustainability 2021, 13, 1254. [Google Scholar] [CrossRef]
  37. Emiru, N.C.; Recha, J.W.; Thompson, J.R.; Belay, A.; Aynekulu, E.; Manyevere, A.; Demissie, T.D.; Osano, P.M.; Hussein, J.; Molla, M.B.; et al. Impact of Climate Change on the Hydrology of the Upper Awash River Basin, Ethiopia. Hydrology 2022, 9, 3. [Google Scholar] [CrossRef]
  38. Gebresellase, S.H.; Wu, Z.; Xu, H.; Muhammad, W.I. Evaluation and selection of CMIP6 climate models in Upper Awash Basin (UBA), Ethiopia. Theor. Appl. Climatol. 2022, 149, 1521–1547. [Google Scholar] [CrossRef]
  39. Schulzweida, U. CDO User Guide. 2021. Available online: https://zenodo.org/record/5614769#.Y3wxceRBxPY (accessed on 14 November 2022). [CrossRef]
  40. Rigby, R.A.; Stasinopoulos, D.M. Generalized additive models for location, scale and shape,(with discussion). Appl. Stat. 2005, 54, 507–554. [Google Scholar] [CrossRef] [Green Version]
  41. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  42. Mann, H.B. Nonparametric Tests against Trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
  43. Yue, S.; Wang, C. The Mann-Kendall test modified by effective sample size to detect trend in serially correlated hydrological series. Water Resour. Manag. 2004, 18, 201–218. [Google Scholar] [CrossRef]
  44. Sen, P.K. Estimates of the regression coefficient based on Kendall’s tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
  45. De Lucas, D.C. Classification Techniques for Time Series and Functional Data. Ph.D. Thesis, Universidad Carlos III de Madrid, Madrid, Spain, 2010. [Google Scholar]
  46. Montero, P.; Vilar, J.A. TSclust: An R Package for Time Series Clustering. J. Stat. Softw. 2014, 62, 1–43. [Google Scholar] [CrossRef] [Green Version]
  47. Kabacoff, R.I. R in Action: Data Analysis and Graphics with R; Manning Publishings Co.: New York, NY, USA, 2015. [Google Scholar]
  48. Hartmann, K.; Krois, J.; Waske, B. E-Learning Project SOGA: Statistics and Geospatial Data Analysis; Department of Earth Sciences, Freie Universitaet Berlin: Berlin, Germany, 2018; Volume 33. [Google Scholar]
  49. Sarvina, Y.; Pluntke, T.; Bernhofer, C. Comparing bias correction methods to improve modelled precipitation extremes. J. Meteorol. Dan Geofis. 2019, 19, 103–110. [Google Scholar] [CrossRef]
  50. Gudmundsson, L.; Bremnes, J.B.; Haugen, J.E.; Engen-Skaugen, T. Downscaling RCM precipitation to the station scale using statistical transformations—A comparison of methods. Hydrol. Earth Syst. Sci. 2012, 16, 3383–3390. [Google Scholar] [CrossRef] [Green Version]
  51. Tabari, H.; Paz, S.M.; Buekenhout, D.; Willems, P. Comparison of statistical downscaling methods for climate change impact analysis on precipitation-driven drought. Hydrol. Earth Syst. Sci. 2021, 25, 3493–3517. [Google Scholar] [CrossRef]
  52. Maraun, D. Bias correcting climate change simulations-a critical review. Curr. Clim. Chang. Rep. 2016, 2, 211–220. [Google Scholar] [CrossRef] [Green Version]
  53. Iturbide, M.; Bedia, J.; Herrera, S.; Baño-Medina, J.; Fernández, J.; Frías, M.; Manzanas, R.; San-Martín, D.; Cimadevilla, E.; Cofiño, A.; et al. climate4r: An r-based open framework for reproducible climate data access and post-processing. Environ. Modell. Softw. 2018. [Google Scholar] [CrossRef]
  54. Boé, J.; Terray, L.; Habets, F.; Martin, E. Statistical and dynamical downscaling of the Seine basin climate for hydro-meteorological studies. Int. J. Climatol. J. R. Meteorol. Soc. 2007, 27, 1643–1655. [Google Scholar] [CrossRef]
  55. Planton, Y.Y.; Guilyardi, E.; Wittenberg, A.T.; Lee, J.; Gleckler, P.J.; Bayr, T.; McGregor, S.; McPhaden, M.J.; Power, S.; Roehrig, R.; et al. Evaluating climate models with the CLIVAR 2020 ENSO metrics package. Bull. Am. Meteorol. Soc. 2021, 102, E193–E217. [Google Scholar] [CrossRef]
  56. Alhamshry, A.; Fenta, A.A.; Yasuda, H.; Kimura, R.; Shimizu, K. Seasonal rainfall variability in Ethiopia and its long-term link to global sea surface temperatures. Water 2020, 12, 55. [Google Scholar] [CrossRef] [Green Version]
  57. Taye, M.T.; Dyer, E.; Charles, K.J.; Hirons, L.C. Potential predictability of the Ethiopian summer rains: Understanding local variations and their implications for water management decisions. Sci. Total Environ. 2021, 755, 142604. [Google Scholar] [CrossRef]
  58. Jungclaus, J.; Fischer, N.; Haak, H.; Lohmann, K.; Marotzke, J.; Matei, D.; Mikolajewicz, U.; Notz, D.; Von Storch, J. Characteristics of the ocean simulations in the Max Planck Institute Ocean Model (MPIOM) the ocean component of the MPI-Earth system model. J. Adv. Model. Earth Syst. 2013, 5, 422–446. [Google Scholar] [CrossRef]
  59. Gutjahr, O.; Putrasahan, D.; Lohmann, K.; Jungclaus, J.H.; von Storch, J.S.; Brüggemann, N.; Haak, H.; Stössel, A. Max planck institute earth system model (MPI-ESM1. 2) for the high-resolution model intercomparison project (HighResMIP). Geosci. Model Dev. 2019, 12, 3241–3281. [Google Scholar] [CrossRef] [Green Version]
  60. Cherchi, A.; Fogli, P.G.; Lovato, T.; Peano, D.; Iovino, D.; Gualdi, S.; Masina, S.; Scoccimarro, E.; Materia, S.; Bellucci, A.; et al. Global mean climate and main patterns of variability in the CMCC-CM2 coupled model. J. Adv. Model. Earth Syst. 2019, 11, 185–209. [Google Scholar] [CrossRef] [Green Version]
  61. Stocker, T.F.; Qin, D.; Plattner, G.K.; Tignor, M.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, B.; Midgley, B. IPCC 2013: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 1535. [Google Scholar]
  62. Masson-Delmotte, V.; Zhai, P.; Pirani, A.; Connors, S.; Péan, C.; Berger, S.; Caud, N.; Chen, Y.; Goldfarb, L.; Gomis, M.; et al. IPCC 2021: The Physical Science Basis, Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2021; p. 2391. [Google Scholar] [CrossRef]
  63. Stasinopoulos, M.; Rigby, B.; Akantziliotou, C. Instructions on How to Use the Gamlss Package in R, 2nd ed.; Gamlss: London, UK, 2008. [Google Scholar]
Figure 1. Map of the upper Awash sub-basin with major rivers and locations of metrological stations.
Figure 1. Map of the upper Awash sub-basin with major rivers and locations of metrological stations.
Climate 10 00185 g001
Figure 2. Fitted distribution for the monthly average precipitation (black line) series of the observed and 10 climate models (indicated with different colors).
Figure 2. Fitted distribution for the monthly average precipitation (black line) series of the observed and 10 climate models (indicated with different colors).
Climate 10 00185 g002
Figure 3. Clusters for the monthly precipitation series for the observed and climate model series (the left figure shows the cluster separation with blue color contrast, and the right figure shows the clustering of climate models and observed data in terms of time series similarities based on seasonal patterns).
Figure 3. Clusters for the monthly precipitation series for the observed and climate model series (the left figure shows the cluster separation with blue color contrast, and the right figure shows the clustering of climate models and observed data in terms of time series similarities based on seasonal patterns).
Climate 10 00185 g003
Figure 4. Clusters for the monthly maximum temperature series for the observed and climate model series (the left figure shows the cluster separation with the blue color contrast, and the right figure shows the clustering of climate models and observed data in terms of time series similarity based on the seasonal pattern).
Figure 4. Clusters for the monthly maximum temperature series for the observed and climate model series (the left figure shows the cluster separation with the blue color contrast, and the right figure shows the clustering of climate models and observed data in terms of time series similarity based on the seasonal pattern).
Climate 10 00185 g004
Figure 5. Clusters for the monthly minimum temperature series for the observed and climate model series (the left figure shows the cluster separation with the blue color contrast, and the right figure shows the clustering of climate models and observed data in terms of the time series similarity based on the seasonal pattern).
Figure 5. Clusters for the monthly minimum temperature series for the observed and climate model series (the left figure shows the cluster separation with the blue color contrast, and the right figure shows the clustering of climate models and observed data in terms of the time series similarity based on the seasonal pattern).
Climate 10 00185 g005
Figure 6. PCA biplot for the monthly precipitation series of the observed and climate model series.
Figure 6. PCA biplot for the monthly precipitation series of the observed and climate model series.
Climate 10 00185 g006
Figure 7. PCA biplot for the monthly maximum temperature series of the observed and climate model series. The Figure on the left is a zoomed-out part of the main biplot on the right.
Figure 7. PCA biplot for the monthly maximum temperature series of the observed and climate model series. The Figure on the left is a zoomed-out part of the main biplot on the right.
Climate 10 00185 g007
Figure 8. PCA biplot for the monthly minimum temperature series of the observed and climate model series. The figure on the left is a zoomed-out part of the main biplot on the right.
Figure 8. PCA biplot for the monthly minimum temperature series of the observed and climate model series. The figure on the left is a zoomed-out part of the main biplot on the right.
Climate 10 00185 g008
Figure 9. Box plots of the monthly precipitation correlation coefficients for each climate mode obtained across all climate stations in UASB.
Figure 9. Box plots of the monthly precipitation correlation coefficients for each climate mode obtained across all climate stations in UASB.
Climate 10 00185 g009
Figure 10. Box plots of monthly Tmax correlation coefficients for each climate mode obtained across all climate stations in UASB.
Figure 10. Box plots of monthly Tmax correlation coefficients for each climate mode obtained across all climate stations in UASB.
Climate 10 00185 g010
Figure 11. Box plots of monthly Tmin correlation coefficients for each climate mode obtained across all climate stations.
Figure 11. Box plots of monthly Tmin correlation coefficients for each climate mode obtained across all climate stations.
Climate 10 00185 g011
Figure 12. Maps of IDW interpolated the monthly precipitation RMSE values of the 12 climate models.
Figure 12. Maps of IDW interpolated the monthly precipitation RMSE values of the 12 climate models.
Climate 10 00185 g012
Table 1. List of names of meteorological stations in the upper Awash sub-basin with their locations (UTM (Universal Transverse Mercator)) and altitude (m) were used for the analysis.
Table 1. List of names of meteorological stations in the upper Awash sub-basin with their locations (UTM (Universal Transverse Mercator)) and altitude (m) were used for the analysis.
No.Station NameLocation (y)Location (x)Altitude
1Addis Ababa Obs *472,248.08996,952.42386
2Addis Alem **432,225.93999,552.912372
3Aleltu516,771.631,016,119.332648
4Ambo **372,449.79993,358.732068
5Asgori *426,775.97971,700.312072
6Boneya460,591.24971,046.062251
7Bui **450,940.01920,899.892054
8ChefeDonsa513,542.53991,537.732392
9Debrezeit *494,500.33965,370.791900
10Ejere528,246.66969,798.672254
11Enselale435,870.79987,532.372000
12Ginchi404,738.43996,808.092132
13Hombole475,209.16925,006.741743
14Huruta *537,697.21900,012.172044
15Melkasa *534,861.76928,532.981540
16Mojo511,901.68951,220.71763
17Nazret531,179.91945,113.491622
18Sebeta459,322.66986,027.982220
19Teji430,354.91976,481.482091
20Tulu Bolo *414,188.26958,456.672100
21Welenchiti547,305.49958,395.381458
22Woliso **388,113.76945,249.622058
N.B: The ones without and with a single star (*) are used for precipitation and the others with one star (*) and a double star (**) are used to differentiate stations used for temperature analysis.
Table 2. List of climate models used in the selection process.
Table 2. List of climate models used in the selection process.
ModelInstituteCountry
MRI-ESM2-0Meteorological Research InstituteJapan
ECEARTH3-CCEC-Earth consortiumSweden
NorESM2-MMNorwegian Climate CenterNorway
TaiESM1Academia SinicaTaiwan
ECEARTH3.VegEC-Earth consortiumSweden
MPI-ESM1.2.HRMax Planck Institute for MeteorologyGermany
ECEARTH3EC-Earth consortiumSweden
CMCC-ESM2Centro Euro-Mediterraneo sui Cambiamenti ClimaticiItaly
GFDL-CM4NOAAUSA
GFDL-ESM4NOAAUSA
INM-CM4-8Institute for Numerical MathematicsRussia
INM-CM5-0Institute for Numerical MathematicsRussia
Table 3. Fitted distributions for all three climatic variables across the four time steps of both the observed and climate model series.
Table 3. Fitted distributions for all three climatic variables across the four time steps of both the observed and climate model series.
Monthly AverageJJASMAMAnnual
Data TypeAICDTAICDTAICDTAICDT
Precipitation
Observed221.62NO359.84LO348.147GA370.72NO
CMCC-ESM2210.09SN2229.72SN2272.3SEP1359.19SN2
ECEARTH3213.9SN2268.16WEI277.63WEI3362.99SN2
ECEARTH3-CC229.69RG257.03GA285.2GA378.79RG
ECEARTH3-Veg228.64SEP1249.17SEP3274.55RG381.82IGAMMA
GFDL-CM4196.06WEI3254.98WEI266.21NO345.16WEI
GFDL-ESM4214.78SEP3233.74SEP2259.64NO364.56IGAMMA
MPI-ESM1-2-HR217.1SEP1263.67WEI3275.74SEP1366.19SEP1
MRI-ESM2232.6NO261.15LO271.72GA381.7NO
NorESM2-MM229.68LO277.26LO295.01GA378.77LO
TaiESM1209.12NO237.44IGAMMA272.38LO358.22NO
INM-CM4-8233.88WEI281.42SEP2238.5WEI3382.97WEI3
INM-CM5-0226.62SN2280.22NO261.28WEI2377.17SN2
Maximum Temperature
Observed45.97SN228.93RG51.31WEI13.19PE2
CMCC-ESM232.12WEI352.57SN237.6LO32.12WEI3
ECEARTH330.55NO44.69WEI353.16GT30.55NO
ECEARTH3-CC30.06SN239.55GG42.85SN230.06SN2
ECEARTH3-Veg26.53NO42.93SN245.23NO26.53NO
GFDL-CM48.59SN237.81RG45.1IGAMMA8.59SN2
GFDL-ESM427.57SN242.36WEI349.34LO27.57SN2
MPI-ESM1-2-HR14.35NO34.3SN257.13NO14.35NO
MRI-ESM210.66NO42.33WEI42.31PE210.66NO
NorESM2-MM37.52SHASH58.02SN255.42WEI337.52SHASH
TaiESM12.78WEI331.94WEI25.06NET2.78WEI3
INM-CM4-814.07WEI314.57IGAMMA31.78PE14.07WEI3
INM-CM5-011.02WEI27.42LO35.8NET11.02WEI
Minimum Temperature
Observed42.35SEP130.9NO34.09NO16.98IGAMMA
CMCC-ESM233.67SEP124.05SN245.87GA33.67SEP1
ECEARTH344.68WEI347.64WEI23.2NET44.68WEI3
ECEARTH3-CC40.14SN252.83WEI31.7GU40.14SN2
ECEARTH3-Veg45.88WEI316.64BCPEo28.47PE45.88WEI3
GFDL-CM43.73IGAMMA-6.15RG28.03IGAMMA3.73IGAMMA
GFDL-ESM431WEI15.08WEI338.77NO31WEI
MPI-ESM1-2-HR21.35IGAMMA10.89SEP33.36NO21.35IGAMMA
MRI-ESM229.13RG9.45GT30.08RG29.13RG
NorESM2-MM50SHASH37.67SEP145.7WEI50SHASH
TaiESM116.68NET12.72LO34.01LO16.68NET
INM-CM4-816.43WEI-6.83WEI33.46SN216.43WEI
INM-CM5-012.78SEP3-7.68SN240IG12.78SEP3
NB: highlighted in grey are models that shared the same distribution as the observed time series across each time period of the analysis. DT = distribution type.
Table 4. Results of the trend analysis for all three climatic variables across the four time steps for the observed and climate model series.
Table 4. Results of the trend analysis for all three climatic variables across the four time steps for the observed and climate model series.
Data TypeAnnualMonthly Ave.JJASMAM
Z-valSen S.Z-valSen S.Z-valSen S.Z-valSen S.
Precipitation
Observed−0.43−0.91−0.43−0.08−0.07−0.060.070.05
CMCC-ESM20.040.090.040.01−1.39−0.321.820.89
ECEARTH31.393.621.390.302.641.06−0.21−0.05
ECEARTH3-CC1.785.611.780.473.141.00−0.39−0.29
ECEARTH3-Veg0.391.510.390.130.640.21−0.36−0.27
GFDL-CM41.683.141.680.261.040.33−0.18−0.05
GFDL-ESM4−0.25−1.02−0.25−0.09−1.43−0.25−0.86−0.45
MPI-ESM1-2-HR0.210.550.210.05−0.25−0.080.610.17
MRI-ESM21.574.521.570.380.50.250.460.31
NorESM2-MM1.213.291.210.270.140.10.140.08
TaiESM1−0.86−1.94−0.86−0.16−1.61−0.50.640.24
INM-CM4-82.874.252.870.351.321.041.040.29
INM-CM5-00.712.710.710.230.00−0.011.140.52
Maximum Temperature
Observed7.360.037.360.031.270.022.380.05
CMCC-ESM22.130.022.130.021.000.030.340.01
ECEARTH34.320.044.320.042.540.051.80.05
ECEARTH3-CC1.290.011.290.013.230.031.000.02
ECEARTH3-Veg3.40.033.40.032.970.040.370.01
GFDL-CM41.320.021.320.020.40.010.790.02
GFDL-ESM42.480.022.480.020.320.012.380.05
MPI-ESM1-2-HR0.080.000.080.000.260.010.000.00
MRI-ESM20.000.000.000.000.690.02−1.47−0.01
NorESM2-MM1.990.021.990.020.260.010.630.01
TaiESM16.630.026.630.021.950.031.690.02
INM-CM4-82.480.032.480.032.850.041.740.02
INM-CM5-02.990.022.990.021.560.021.480.02
Minimum Temperature
Observed1.80.021.80.022.320.031.110.02
CMCC-ESM21.570.021.570.022.370.020.640.01
ECEARTH32.80.062.80.063.380.073.010.03
ECEARTH3-CC1.410.031.410.034.020.051.690.03
ECEARTH3-Veg4.020.044.020.042.640.051.890.01
GFDL-CM41.220.011.220.012.480.02−0.79−0.01
GFDL-ESM44.760.044.760.044.130.022.460.04
MPI-ESM1-2-HR1.000.011.000.010.790.01−0.87−0.01
MRI-ESM25.90.045.90.045.580.042.090.03
NorESM2-MM0.870.010.870.011.760.020.480.02
TaiESM10.790.010.790.016.240.021.690.02
INM-CM4-88.950.028.950.026.830.021.060.01
INM-CM5-02.380.022.380.023.730.012.010.03
NB: the ones shaded in gray show that the models have similar trend test outcomes, i.e, if both the observed and the model are significant/insignificant, then the z-values are highlighted.
Table 5. Summary of model performances for all three climatic variables across four performance measures for the observed and climate model series.
Table 5. Summary of model performances for all three climatic variables across four performance measures for the observed and climate model series.
Ranks
Data Type R 2 RMSEMAEBIASSum of Rank
Precipitation
ECEARTH3622414
ECEARTH3-Veg2611019
GFDL-ESM43181224
CMCC-ESM21059125
INM-CM4-81107826
NorESM2-MM5116527
MPI-ESM1-2-HR4851128
MRI-ESM2993728
TaiESM17310929
GFDL-CM48412630
INM-CM5.012124230
ECEARTH3-CC11711332
Maximum Temperature
TaiESM1611311
INM-CM5-0422614
ECEARTH3149115
INM-CM4-8733417
ECEARTH3-CC2741023
MPI-ESM1-2-HR865524
GFDL-ESM43561226
ECEARTH3-Veg988227
GFDL-CM41097733
MRI-ESM2512111139
CMCC-ESM2121010840
NorESM2-MM111112943
Minimum Temperature
GFDL-CM441117
MPI-ESM1-2-HR1231117
TaiESM1344617
INM-CM4-85321222
ECEARTH3-Veg288725
INM-CM5-01067427
MRI-ESM21275327
GFDL-ESM49561030
CMCC-ESM2699832
NorESM2-MM81212234
ECEARTH371110937
ECEARTH3-CC111011537
Table 6. Summary rank based on seven criteria used for evaluating the three climate model time series.
Table 6. Summary rank based on seven criteria used for evaluating the three climate model time series.
ModelPDF RankTrend RankPM RankCluster RankPCA RankSpatCorrSpaRMSESum of RankFinal Rank
Precipitation
MRI-ESM21186113211
GFDL-ESM45134352232
MPI-ESM1-2-HR5171544273
GFDL-CM451108231304
ECEARTH351019425365
NorESM2-MM2165888386
ECEARTH3-Veg51212777417
INM-CM5-051113101010508
CMCC-ESM25141011119519
ECEARTH3-CC410127666519
INM-CM4-85125299125411
TaiESM1219111212115812
Maximum Temperature
ECEARTH34633431241
INM-CM5-04225843282
INM-CM4-84644524293
MPI-ESM1-2-HR4961612293
TaiESM14212978335
ECEARTH3-CC21258155386
ECEARTH3-Veg4686266386
CMCC-ESM2421171087498
GFDL-CM419911399519
GFDL-ESM421710111010519
NorESM2-MM42129711125711
MRI-ESM24910121212117012
Minimum Temperature
GFDL-CM42114211121
MPI-ESM1-2-HR1521667282
TaiESM15137743303
GFDL-ESM42983135314
MRI-ESM25972422314
CMCC-ESM22195588386
INM-CM5-05966854437
INM-CM4-85748979498
NorESM2-MM551011396498
ECEARTH3-Veg57591010105610
ECEARTH3-CC5112101111116111
ECEARTH35911121212127312
NB: PDFR (PDF Rank), TrendR (Trend Rank), PMR (performance measure rank), ClusterR (cluster rank), PCAR (PCA rank), SpatCorr (spatial Correlation), SpaRMSE (spatial RMSE), SRank (sum of rank), FRank (final rank).
Table 7. Summary results for all three climatic variables.
Table 7. Summary results for all three climatic variables.
Data TypePrecip RankTmax RankTmin RankSum of RankFinal Rank
MPI-ESM1-2-HR33281
GFDL-CM4491142
GFDL-ESM4294153
INM-CM5-0827174
MRI-ESM21124174
ECEARTH35112186
TaiESM11253207
INM-CN4-81138228
CMCC-ESM2986239
ECEARTH3-Veg7610239
NorESM2-MM61182511
ECEARTH3-CC96112612
Table 8. Mean annual areal-averaged future climate change over UASB for three time periods and two scenarios.
Table 8. Mean annual areal-averaged future climate change over UASB for three time periods and two scenarios.
ScenariosNearMidEnd
Precipitation (% Δ )SSP5-8.5(−2.1–5.3) 1.3(6.1–16.1) 10.6(11.8–29.4) 19.7
SSP2-4.5(−0.7–6.8) 2.7(2.0–11.9) 6.4(2.2–15.0) 7.9
SSP1-2.6(−1.2–9.0) 3.4(1.4–16.6) 8.2(−0.3–14.2) 6.2
Tmax (°C)SSP5-8.5(0.8–1.5) 1.2(1.7–2.3) 2.0(2.8–3.2) 3.0
SSP2-4.5(0.7–1.5) 1.1(1.3–2.0) 1.6(1.7–2.3) 2.0
SSP1-2.6(0.7–1.4) 1.1(0.9–1.7) 1.3(0.9–1.8) 1.3
Tmin (°C)SSP5-8.5(1.0–1.9) 1.4(2.2–2.9) 2.5(3.7–4.2) 4.0
SSP2-4.5(1.0–1.7) 1.3(1.6–2.3) 2.0(2.1–2.7) 2.4
SSP1-2.6(0.9–1.7) 1.3(1.2–2.0) 1.6(1.1–2.0) 1.6
NB: numbers in the brackets indicate uncertainty limits of mean changes of the mean annuals calculated by using the standard deviation across the years of each series (predicted and observed).
Table 9. Seasonal areal-averaged precipitation future climate change (% Δ ) over UASB for three time periods and two scenarios.
Table 9. Seasonal areal-averaged precipitation future climate change (% Δ ) over UASB for three time periods and two scenarios.
MAM JJAS
NearMidEndNearMidEnd
SSP5-8.5(−0.15–0.46) 0.03(−0.10–0.47) 0.08(−0.04–0.58) 0.15(−0.12–0.02) −0.06(−0.06–0.11) 0.01(0.02–0.16) 0.08
SSP2-4.5(−0.13 −0.6) 0.08(−0.08–0.45) 0.09(−0.11–0.48) 0.08(−0.11–0.04) −0.04(−0.10–0.07) −0.02(−0.06–0.08) 0.00
SSP1-2.6(−0.14–0.6) 0.07(−0.06–0.56) 0.14(−0.07–0.43) 0.09(−0.11–0.06) −0.04(−0.09–0.09) −0.01(−0.08–0.07) −0.02
Table 10. Mean annual future climate change at the selected four meteorological stations in the UASB for the three time periods and two scenarios.
Table 10. Mean annual future climate change at the selected four meteorological stations in the UASB for the three time periods and two scenarios.
StationsScenariosNearMidEnd
Precip (% Δ )
AASSP5-8.5(4.3–9.7) 7.4(6.0–25.9) 17.5(28.6–33.2) 31.3
SSP2-4.5(24.6–26.3) 25.3(10.5–13.3) 11.7(10.0–10.0) 10.0
EjereSSP5-8.5(−1.6–13.2) 3.8(21.1–23.9) 22.8(29.7–46.9) 36.2
SSP2-4.5(11.1–37.9) 20.9(5.2–28.7) 14.0(6.9–24.1) 13.4
HomboleSSP5-8.5(−23.4–63.7) −2.6(3.4–77.7) 22.0(7.5–123.3) 36.6
SSP2-4.5(−14.5–93.9) 11.4(−12.5–93.7) 14.1(−12.0–86.1) 12.6
GinchiSSP5-8.5(−3.3–5.1) −0.1(14.2–15.1) 14.6(18.0–34.6) 24.5
SSP2-4.5(6.3–26.5) 14.1(2.3–17.8) 8.4(0.2–21.3) 8.5
Tmax (°C)
AASSP5-8.5(1.0–1.3) 1.2(2.4–2.4) 2.4(3.3–3.6) 3.5
SSP2-4.5(0.8–0.9) 0.8(2.1–2.3) 2.2(2.6–2.8) 2.7
BuiSSP5-8.5(0.6–2.1) 1.3(2.4–3.4) 2.9(3.8–4.6) 4.2
SSP2-4.5(0.5–1.6) 1.1(2.1–3.3) 2.7(2.7–4.0) 3.3
DebrezeitSSP5-8.5(1.3–1.5) 1.4(2.9–3.2) 3.0(4.0–4.6) 4.3
SSP2-4.5(1.0 -1.1) 1.1(2.7–2.8) 2.7(3.5–3.5) 3.5
TuluboloSSP5-8.5(−0.1–2.8) 1.4(1.3–3.8) 2.6(2.5–4.9) 3.7
SSP2-4.5(−0.3–2.5) 1.1(1.0–3.8) 2.4(1.5–4.4) 3.0
Tmin (°C)
AASSP5-8.5(1.5–1.7) 1.6(2.6–2.7) 2.7(3.8–4.4) 4.1
SSP2-4.5(1.1–1.2) 1.2(2.2–2.3) 2.2(2.6–3.0) 2.8
BuiSSP5-8.5(1.3–3.8) 2.5(3.5–5.4) 4.5(5.4–7.0) 6.2
SSP2-4.5(1.0–3.0) 2.0(2.9–5.1) 4.0(3.3–5.8) 4.6
DebrezeitSSP5-8.5(1.5–1.9) 1.7(3.2–3.3) 3.2(4.4–4.9) 4.7
SSP2-4.5(1.2–1.3) 1.2(2.7–2.8) 2.8(3.1–3.6) 3.3
TuluboloSSP5-8.5(0.9–1.8) 1.4(1.9–4.1) 3.0(3.2–5.2) 4.2
SSP2-4.5(0.8–1.4) 1.1(1.4–3.9) 2.7(1.8–4.5) 3.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Balcha, Y.A.; Malcherek, A.; Alamirew, T. Understanding Future Climate in the Upper Awash Basin (UASB) with Selected Climate Model Outputs under CMIP6. Climate 2022, 10, 185. https://doi.org/10.3390/cli10120185

AMA Style

Balcha YA, Malcherek A, Alamirew T. Understanding Future Climate in the Upper Awash Basin (UASB) with Selected Climate Model Outputs under CMIP6. Climate. 2022; 10(12):185. https://doi.org/10.3390/cli10120185

Chicago/Turabian Style

Balcha, Yonas Abebe, Andreas Malcherek, and Tena Alamirew. 2022. "Understanding Future Climate in the Upper Awash Basin (UASB) with Selected Climate Model Outputs under CMIP6" Climate 10, no. 12: 185. https://doi.org/10.3390/cli10120185

APA Style

Balcha, Y. A., Malcherek, A., & Alamirew, T. (2022). Understanding Future Climate in the Upper Awash Basin (UASB) with Selected Climate Model Outputs under CMIP6. Climate, 10(12), 185. https://doi.org/10.3390/cli10120185

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

Article Metrics

Back to TopTop