Climate and extreme rainfall events in the Mono river basin (West Africa): investigating future changes with Regional Climate Models

This study characterizes the future changes in extreme rainfall and air temperature in the Mono river basin where the main economic activity is weather dependent and local populations are highly vulnerable to natural hazards, including flood inundations. Daily precipitation and temperature from observational datasets and Regional Climate Models (RCMs) output from REMO, RegCM, HadRM3, and RCA were used to analyze climatic variations in space and time, and fit a GEV model to investigate the extreme rainfalls and their return periods. The results indicate that the realism of the simulated climate in this domain is mainly controlled by the choice of the RCMs. These RCMs projected a 1 to 1.5 ◦C temperature increase by 2050 while the projected trends for cumulated precipitation are null or very moderate and diverge among models. Contrasting results were obtained for the intense rainfall events, with RegCM and HadRM3 pointing to a significant increase in the intensity of extreme rainfall events. The GEV model is well suited for the prediction of heavy rainfall events although there are uncertainties beyond the 90th percentile. The annual maxima of daily precipitation will also increase by 2050 and could be of benefit to the ecosystem services and socioeconomic activities in the Mono river basin but could also be a threat.


Study Area
The Mono river basin is in West Africa and extends over 560 km from the north to the south [2]. This transboundary watershed covers 15,680 km 2 and is shared by Benin and Togo ( Figure 1). This watershed is home to the Nangbéto hydroelectric dam that has been providing electricity to Benin and Togo since September 1987. As this hydroelectric dam is modifying the hydrology of the river, the main part of this study will focus on the Mono river basin upstream of the dam.
This watershed is patterned in the south by floodplains and plateaus, and higher landforms in the north and north-west, e.g., the Atakora Mountains with a height of 800 m and their southern extensions that are the Togo mountains.
Precipitation in Togo and Benin is constrained by the organization of the atmospheric circulation over West Africa as a whole, i.e., both that of the lower layers (humid monsoon flow from the southwest, and dry Harmattan flow from the northeast) and the flow of air in the lower troposphere and upper atmosphere (respectively African Easterly Jet (AEJ) and Tropical Easterly Jet (TEJ)) [2]. At the scale of the Mono river basin, this circulation is strongly linked to the energy gradients between the coastal plain and the Gulf of Guinea. The precipitations are also modulated by the interaction of the southwesterly humid monsoon flow and the landforms in the northwestern part of the watershed. The average annual precipitation ranges from 900 mm in the southeast, i.e., in the plains along the dry diagonal in the south of Togo, to 1200 mm in the northwestern uplands. There are two climatic domains in the basin: The subequatorial domain in the southern part of the basin with two rainy and two dry seasons and the tropical domain in the northern part with one rainy season and one dry season.
The hydrological regime is of a humid tropical type and consistent with the tropical climate that governs it. The uneven distribution of precipitation combined with heavy rainfall events over very short periods leads to recurring flooding events. These phenomena are aggravated by the increase in surface runoff due to land cover degradation [30]. As a result, the local population or riparian and their properties and livelihoods are more vulnerable than the population living upstream of the basin as evidenced by the severe floods of the past 12 years in the lower valley of the Mono river basin.

Observational Climate Data
We used daily rainfall data from 14 rain gauge stations located in the Mono river basin and its close vicinity, and daily near-surface air temperature from two meteorological stations, namely the Bohicon and Nangbéto stations. Those observed data were used to select the closest model to the observation (reality) for projection. Thus, the model closest to the reality was used. The data cover a 50-year period  and are provided by the Agency for the Air Navigation Safety in Africa

Observational Climate Data
We used daily rainfall data from 14 rain gauge stations located in the Mono river basin and its close vicinity, and daily near-surface air temperature from two meteorological stations, namely the Bohicon and Nangbéto stations. Those observed data were used to select the closest model to the Water 2020, 12, 833 5 of 27 observation (reality) for projection. Thus, the model closest to the reality was used. The data cover a 50-year period (1961-2010) and are provided by the Agency for the Air Navigation Safety in Africa and Madagascar (ASECNA) in Cotonou-Benin, the Direction de la Météorologie Nationale (DMN) of Togo, the Climatology Research Centre (CRC) of the University of Bourgogne in France, and the Global Historical Climate Network (GHCN, [31]).
We compared the gauged rainfall data with the gridded data from the RCM models. The rain-gauge data were spatially interpolated to a 0.5 • resolution upstream of the Nangbéto hydro-electric dam using an ordinary kriging block (Figure 1). A climatological variogram was considered and an adjustment was made to a spherical variogram [22]. The range parameter of the adjusted model is 65 km and indicates the average distance of precipitation decorrelation between two stations. In contrast, the arithmetic mean of the two observed air temperature stations data was made and then used to compare to the model dataset.
In order to describe the climates that prevail in West Africa and compare them with the RCM simulations, gridded monthly data of precipitation and near-surface air temperature from 3 additional datasets were extracted over the region extending over 5-14 • N and 1.2-3 • E. These sets of data covering the years from 1988 to 2009 at a 0.5 • spatial resolution are the GPCP Global Precipitation Climatology Project data [32], the near-surface air temperature data from CRU (Climatic Research Unit, [33]). The CRU data is chosen based on its multidiscipline application, namely in applied climatology, biogeochemical modelling, hydrology, and agricultural meteorology domains [33] and its performance [34].

Data from Regional Climate Models
Precipitation and near-surface air temperature data from the selected four RCMs developed in the framework of the ENSEMBLES European project [35] are listed in Table 1. They are extracted from the AMMA-ENSEMBLES database [17]. These RCMs simulations are centered on West Africa (19.8 • S-35.2 • N and 35.2 • W-31.2 • E), have a horizontal resolution of 50 km, and cover the period 1951-2050 (1980-2050 for the ICTP model). These RCMs simulations are driven by two CMIP3 global climate models (GCMs), namely ECHAM5 for MPI and ICTP models and HadCM3 for HC (METO) and SMHI models, following the IPCC A1B scenario. We used these RCM simulations at a daily time step for a region extending over 5-14 • N and 1.2 • W-3 • E. METO and SMH simulations follow a 360-day calendar. The 18-year period spanning 1988-2005 was used for the comparison of the precipitation, and air temperature RCMs data with observations (GPCP for precipitation and CRU for temperature) were considered for the air temperature. We extracted the CORDEX, CMIP5, and AMMA-ENSEMBLES outputs over the Mono river basin at the outlet of Nangbéto, computed the spatial ensemble averages of each data type, and compared them with the interpolated (observed) data. Since the streamflow gauge data recorded are only available over 1988-2010, we used the period 1988 to 2010 (23 years) as the reference period and the period 2028 to 2050 (23 years) as the scenario period.

Statistical Methods
Observed and simulated climate data were analyzed for the sub-region and the Mono river basin itself. At the sub-regional scale, seasonal average maps were used to compare simulated precipitation and mean air temperature fields with the observational data over 1988-2010, and to assess the projected future changes for 2028-2050. The RCM whose outputs in the reference period is closer to observations can be considered for future projection.
At the catchment scale, the observed and simulated climate data were also spatially averaged and statistically processed. First, monthly averages of rainfall and air temperature were computed to describe the ability of the regional climate models to reproduce the present-day climatic regimes and to assess the projected future changes (2028-2050). The average frequency of the number of rainy days (i.e., a day with a rainfall amount higher than 1 mm) in each grid point in the Mono river basin was also computed. The contribution of the different classes of daily rainfall intensity to the total precipitation recorded on average in the watershed was determined from the calculation of deciles. The 99th percentile is usually used to characterize the extreme events [36][37][38][39]. This is the reason why the 99th percentile is used to assess the importance of heavy rains in the basin.
Observed hydroclimatic series do not always meet the assumption of data stationarity [40,41]. This is often the case for extreme rainfall and flood discharge. In the literature, some scholars [42,43] have suggested describing the distribution of extreme values with three different types of probability distribution for extreme order statistics while others [44][45][46][47] suggested using a generalized extreme value (GEV) distribution, which combines the three types of Fisher-Tippett extreme value distributions [48]. The GEV parameters use the generalized maximum likelihood (GML) method, which includes an additional constraint on the shape parameter to eliminate potentially invalid values of this parameter based on Tramblay et al.'s [4] study.
In this study, we used a GEV distribution to characterize extreme events. In the case of non-stationarity, the GEV parameters are not constant but are dependent on time or other covariates. To avoid having a value of zero for the scale parameter, αt, we considered the parameterization as φt = log(αt) and thus restricted our analysis to linear dependencies [49]. However, many other types of dependencies can be linearized [41,44,48,50]. The return periods of a GEV model are defined by the cumulative distribution function (CDF), F(x), or the probability density function (PDF), f(x), hat are given by the formula: where Pr is the probability of the occurrence of an event and X stands for the random variable, which in this case is a series of annual daily maxima of rainfall. The quantile value, Q T , is a magnitude of the extreme event that has a probability of 1/T and can be exceeded by a single event. A return period, T, associated with a quantile value, Q T , was adapted from Stedinger et al.'s [51] equation and can be expressed as: The selection of the best GEV model adjustment was made by using the deviance test, which is based on the log-likelihood [20]. The spatial distribution of the observed rainfall (1988-2005 annual rainfall average, GPCP data) in the Mono river basin and surrounding areas is presented in Figure 2d. The observed rainfall pattern shows a diagonal axis of a high rainfall amount (over 1320 mm per year) from south-west to north-east along with the Togo and Atakora mountain ranges. A strong decrease in rainfall is also shown in the north towards the Sahelian zone (<600 mm/year). The GPCp reanalysis data captures the observed precipitation gradient over West Africa well [52] and this finding agrees with the study of Ntajal et al. [49] in Mono basin.

Spatial Distribution of the Annual Rainfall and Air Temperature for Three Datasets
The ensemble averages were computed for CMIP5 and RCP8.5 for CORDEX. The CORDEX ( Figure 2c) and CMIP5 ( Figure 2b) display a greater rainfall estimation in the historical period compared to the AMMA-Ensemble ( Figure 2a) and GPCP (Figure 2d), especially at the southern part. The validation made (second column of Figure 2) exhibits that CMIP5 and CORDEX overestimate the rainfall, which is about 340 mm/year. Nevertheless, the AMMA-Ensemble gives an estimation closer to the reanalysis (±20 mm/year).
There is a projected general decrease varying by 16 to 32 mm in the annual total precipitation. For instance, the CMIP5 model presents a decrease of about 24 mm and less than 24 mm, respectively, in the central and southern parts of the basin while the CORDEX estimates a decrease between 24 and 32 mm in the central and southern part of the basin. Nevertheless, in the northern part of the basin, an increase in the annual total precipitation of around 16 mm and in the range of 16 to 24 mm for the CMIP5 and CORDEX models, respectively, is projected. Finally, the AMMA-ENSEMBLES displays a decrease in the projection change ranging from 0 to 32 mm/year. The temperature gradient was well shown by all the datasets, which increase in the northern direction (first column of Figure 3). The observed mean air temperature extracted from Climate Research Unit (CRU) varies slightly, with a maximum of 29.5 • C in the far north (the Sahel in Niger) and a minimum of 26 • C in altitude on the Togolese ridge ( Figure 3d). The AMMA-Ensemble and CORDEX data give an underestimation of the temperature ranging from 0 to 1 • C while the CMIP5 presents an underestimation between 1 and 2 • C (second column of Figure 3) and the lowest rise is given by AMMA-ENSEMBLES. Unlike the rainfall, there is a very good agreement in the projected changes (increase) of air temperature from the AMMA-Ensemble, CMIP5, and CORDEX dataset.
The warming is expected to vary between 1.4 • C in the south to 1.8 • C towards the north. However, AMMA-ENSEMBLES gives the lowest projected changes among all the models (third column).
The performance of CORDEX RCMs in simulating the spatial distribution of the main precipitation and temperature features as well as the occurrence of other physical phenomena over West Africa was demonstrated by Gbobaniyi et al. [53]. Diallo et al. [54] also showed the importance of the ensemble mean of multiple models. They proved that the RCMs-GCMs alone are biased and these biases are reduced by averaging all RCM simulations, suggesting that multi-model RCM ensembles based on different driving GCMs help to compensate systematic errors from both the nested and the driving models [54]. This justified the use of multiple models. According to IPCC assessment reports, the temperature is projected to increase over West Africa for the end of the 21st century from the global climate simulation range between 3 and 6 • C above the late 20th century baseline depending on the emission scenario [55]. However, they also revealed that the onset and length of precipitation could also change.
Water 2020, 12, x FOR PEER REVIEW 8 of 27 century from the global climate simulation range between 3 and 6 °C above the late 20th century baseline depending on the emission scenario [55]. However, they also revealed that the onset and length of precipitation could also change.   AMMA-ENSEMBLES data permits a better understanding of the processes underlying the West-African Monsoon (WAM) [56]. It provides and advances an understanding of the atmospheric processes over West Africa [57]. AMMA data were used in climate change assessment in the region.
The finding of the current study is in accordance with Angelina et al.'s [58] study on the river Niger basin West Africa. They confirmed that all RCMs of AMMA-ENSEMBLES experiment displays an increase in air temperature while there was no consensus among them in precipitation variables.

Spatial Distribution of the Averaged Annual Precipitation and Air Temperature for AMMA-ENSEMBLES Models
An analysis of Figure 4 shows the simulated AMMA-ENSEMBLES model's precipitation dataset. The total annual precipitation data for the reference period (1988-2010 in the first raw) and projected period (2028-2050 in the second raw) and the changes in projection in the third raw are presented in Figure 4 for the MPI, METO, SMH, and ICT models. The total annual precipitation over the basin varies from 1200 to 1800 mm/year whatever the model. The total annual precipitation simulated in the reference period (1988-2010) ranges from 600 to 1800 mm while the observed annual precipitation is between 680 and 1640 mm over the chosen window. The spatial distribution of precipitation is AMMA-ENSEMBLES data permits a better understanding of the processes underlying the West-African Monsoon (WAM) [56]. It provides and advances an understanding of the atmospheric processes over West Africa [57]. AMMA data were used in climate change assessment in the region.
The finding of the current study is in accordance with Angelina et al.'s [58] study on the river Niger basin West Africa. They confirmed that all RCMs of AMMA-ENSEMBLES experiment displays an increase in air temperature while there was no consensus among them in precipitation variables.

Spatial Distribution of the Averaged Annual Precipitation and Air Temperature for AMMA-ENSEMBLES Models
An analysis of Figure 4 shows the simulated AMMA-ENSEMBLES model's precipitation dataset. The total annual precipitation data for the reference period (1988-2010 in the first raw) and projected period (2028-2050 in the second raw) and the changes in projection in the third raw are presented in Figure 4 for the MPI, METO, SMH, and ICT models. The total annual precipitation over the basin varies from 1200 to 1800 mm/year whatever the model. The total annual precipitation simulated in the reference period (1988-2010) ranges from 600 to 1800 mm while the observed annual precipitation is between 680 and 1640 mm over the chosen window. The spatial distribution of precipitation is almost similar in the reference period and the scenario period (2028-2050). In the south, a diagonal band of relatively low rainfall (<1000 mm/year) follows the Ghanaian-Togolese coast in relationto the diffluence of the monsoon winds and the seasonal rise of coastal cold water. However, we noted minor discrepancies, especially in the number of rainy days, as shown in Table 2, for METO and SMH. MPI and ICT simulate the same number of rainy days in the reference period and the scenario or projection period. By comparing it to the reference period (1988-2010) data, the models METO and SMH project a reduction in rainfall amounts during the projection period (2028-2050) by 2% and 4%, respectively [55]. In contrast, the fifth IPCC assessment reports a slight increase in precipitation over West Africa. Nevertheless, they could not be able to show any trend in precipitation in this zone.
almost similar in the reference period and the scenario period (2028-2050). In the south, a diagonal band of relatively low rainfall (<1000 mm/year) follows the Ghanaian-Togolese coast in relationto the diffluence of the monsoon winds and the seasonal rise of coastal cold water. However, we noted minor discrepancies, especially in the number of rainy days, as shown in Table 2, for METO and SMH. MPI and ICT simulate the same number of rainy days in the reference period and the scenario or projection period. By comparing it to the reference period (1988-2010) data, the models METO and SMH project a reduction in rainfall amounts during the projection period (2028-2050) by 2% and 4%, respectively [55]. In contrast, the fifth IPCC assessment reports a slight increase in precipitation over West Africa. Nevertheless, they could not be able to show any trend in precipitation in this zone.    An analysis of the variations between the two periods shows moderate but contrasting changes from one model to the others, and across the domain. METO and ICT have projected a slight increase in rainfall ranging between 2% and 5% by 2050. Meanwhile, MPI projection shows no rainfall change whereas SMH has projected a strong rainfall deficit amounting up to 15% to 20% by 2050 in the south of the Mono river basin. These precipitation changes are associated with an increase in air temperature. This result is in accordance with the IPCC report [55] over the West Africa region. The same conclusion was made by Sylla et al. [59]. Figure 5 shows for the Mono river basin and surrounding areas, the spatial distribution of the simulated mean air temperature. It came out that the spatial distribution of the mean air temperature is well reproduced by the RCMs although most of these RCMs, especially SMH, underestimate the values. MPI simulates temperature variations closer to the observations and a much better spatial distribution. However, MPI still overestimates the average minimum air temperature (e.g., 27.5 • C on average for the whole catchment area) compared to the observed data (26.5 • C).
An analysis of Figure 5 for air temperature changes between the projection and reference periods confirms warming of a 1 • C magnitude order whatever the model. This fairly corroborates the 1.22 • C temperature increase obtained with the historical data and is also in good agreement with the SRES-A1B scenario. METO is the only RCM that projected a warming of more than 1.5 • C by 2050 all over the region, except on the coast of the Gulf of Guinea. Thus, the control of the large-scale driver is key to reproducing an accurate simulation of these dry zone features whose origin is mainly related to ocean-atmosphere interactions along the Ghana-Togo coast [60]. This finding agrees with some previous studies in West Africa. For instance, Sylla et al. [59] reported a gradual warming spatial variable reaching 0.5 • C per decade in recent years over West Africa, with a recovery from drought conditions. However, the amount of total annual precipitation of this recent decade is still lower compared with the decade before the 1970s' and 1980s' drought episodes [61,62]. Moreover, they also confirmed the increasing trend in the projected air temperature, and this could be in the range of 1.5 to 6.5 • C depending on the emission scenario. Furthermore, the projection of precipitation is with a lot of uncertainties. The West Africa region could experience change in the range of ±30% in the future associated with more intense extremes in the future climate but to a lesser extent [55,59].  Figure 6 shows the average rainfall and air temperature regimes in the Mono river basin upstream of Nangbéto as obseved and modeled by the four RCMs.   As shown in Figure 6, the four RCMs reproduce the seasonal cycle of precipitation well. This has been proven earlier by some previous works in the region [54,57,63]. In the middle of the rainy season, MPI, SMH, and ICT overestimated the rainfall amounts while they also underestimated the rainfall amounts in January-February, which corresponds to the dry season. The smallest biases were obtained with ICT, thus making ICT the model that best replicated the rainfall regime of this watershed.

Average Rainfall and Temperature Regimes in the Mono River Basin at the Nangbéto Outlet
In most cases, the RCMs satisfactorily reproduce the temperature seasonality in the watershed. The simulated temperature values are underestimated, except MPI, which overestimates these values from April to July. These biases can significantly influence hydrological impact simulations as they will lead to an under(over)estimation of the high streamflow because of an under(over)estimation of the evaporative loss. Figure 7 reveals the annual cycle of precipitation simulated by the RCMs for the reference period (1988-2010) and the projection period (2028-2050). For SMH, a moderate decrease in the rainfall regime is shown, especially from March to September. SMH fails to reproduce the observed rainfall regime. The other RCMs rarely show a significant decrease, except METO from June to August. A Fisher's test indicates that the changes of the rainfall regime between the reference and projection periods are not statistically significant at a 95% confidence level. Therefore, we assume that only small changes are to be expected in the hydrological regimes by 2050. As shown in Figure 6, the four RCMs reproduce the seasonal cycle of precipitation well. This has been proven earlier by some previous works in the region [54,57,63]. In the middle of the rainy season, MPI, SMH, and ICT overestimated the rainfall amounts while they also underestimated the rainfall amounts in January-February, which corresponds to the dry season. The smallest biases were obtained with ICT, thus making ICT the model that best replicated the rainfall regime of this watershed.
In most cases, the RCMs satisfactorily reproduce the temperature seasonality in the watershed. The simulated temperature values are underestimated, except MPI, which overestimates these values from April to July. These biases can significantly influence hydrological impact simulations as they will lead to an under(over)estimation of the high streamflow because of an under(over)estimation of the evaporative loss. Figure 7 reveals the annual cycle of precipitation simulated by the RCMs for the reference period (1988-2010) and the projection period (2028-2050). For SMH, a moderate decrease in the rainfall regime is shown, especially from March to September. SMH fails to reproduce the observed rainfall regime. The other RCMs rarely show a significant decrease, except METO from June to August. A Fisher's test indicates that the changes of the rainfall regime between the reference and projection periods are not statistically significant at a 95% confidence level. Therefore, we assume that only small changes are to be expected in the hydrological regimes by 2050. Water 2020, 12, x FOR PEER REVIEW 14 of 27  This warming is relatively uniform throughout the year and ranges between +1 and +1.5 °C (a little beyond this range for METO) in accordance with some previous studies [55,59]. Overall, there is a good fit between the seasonal cycles of precipitation and air temperature simulated by the four RCMs, but some of these models still feature biases. The precipitation biases are more dependent on the RCM than the driving GCM. Sylla et al. made the same conclusion [54,59] over West Africa. Indeed, the strongest overestimations of rainfall during the monsoon seasons were obtained with the MPI and SMH experiments. These RCMs were not driven by the same GCM (ECHAM5 for MPI and HadCM3 for SMH). In addition, both ICT and MPI are forced by ECHAM5  This could amplify the evaporative loss in the future and potentially induce a decrease in surface runoff. This warming is relatively uniform throughout the year and ranges between +1 and +1.5 • C (a little beyond this range for METO) in accordance with some previous studies [55,59].  This warming is relatively uniform throughout the year and ranges between +1 and +1.5 °C (a little beyond this range for METO) in accordance with some previous studies [55,59]. Overall, there is a good fit between the seasonal cycles of precipitation and air temperature simulated by the four RCMs, but some of these models still feature biases. The precipitation biases are more dependent on the RCM than the driving GCM. Sylla et al. made the same conclusion [54,59] over West Africa. Indeed, the strongest overestimations of rainfall during the monsoon seasons were obtained with the MPI and SMH experiments. These RCMs were not driven by the same GCM (ECHAM5 for MPI and HadCM3 for SMH). In addition, both ICT and MPI are forced by ECHAM5 Overall, there is a good fit between the seasonal cycles of precipitation and air temperature simulated by the four RCMs, but some of these models still feature biases. The precipitation biases are more dependent on the RCM than the driving GCM. Sylla et al. made the same conclusion [54,59] over West Africa. Indeed, the strongest overestimations of rainfall during the monsoon seasons were obtained with the MPI and SMH experiments. These RCMs were not driven by the same GCM (ECHAM5 for MPI and HadCM3 for SMH). In addition, both ICT and MPI are forced by ECHAM5 but the ICT simulates a more realistic annual precipitation cycle than MPI. Thus, the projection biases originate from both the RCMs and their driven GCMs.

Frequency Analysis of Rainfall Distribution in the Mono River Basin
The cumulative frequency of daily rainfall intensity according to the four RCMs is illustrated in Figure 9. This figure compares the distributions of rainfall intensity between observations and the RCM simulations in the reference period, and between the RCM simulations in the reference and projection periods. For most RCMs, the sill of the cumulative daily distribution is obtained after 200 days, except METO, for which the sill is obtained after 250 days. For SMH and ICT, the distribution of daily rainfall frequencies agreed well with the observations, with many days recording less than 20 mm, and a number of wet days (about 200) similar to the observations. METO shows many wet days with a very high number of medium intensity days (5-10 mm). By contrast, the MPI model underestimates precipitation occurrence, especially for intensity values around 5 to 15 mm/day. The study of Tossou et al. [64] also provides a characterization and analysis of rainfall variability in the Mono river watershed complex.
Water 2020, 12, x FOR PEER REVIEW 15 of 27 but the ICT simulates a more realistic annual precipitation cycle than MPI. Thus, the projection biases originate from both the RCMs and their driven GCMs.

Frequency Analysis of Rainfall Distribution in the Mono River Basin
The cumulative frequency of daily rainfall intensity according to the four RCMs is illustrated in Figure 9. This figure compares the distributions of rainfall intensity between observations and the RCM simulations in the reference period, and between the RCM simulations in the reference and projection periods. For most RCMs, the sill of the cumulative daily distribution is obtained after 200 days, except METO, for which the sill is obtained after 250 days. For SMH and ICT, the distribution of daily rainfall frequencies agreed well with the observations, with many days recording less than 20 mm, and a number of wet days (about 200) similar to the observations. METO shows many wet days with a very high number of medium intensity days (5-10 mm). By contrast, the MPI model underestimates precipitation occurrence, especially for intensity values around 5 to 15 mm/day. The study of Tossou et al. [64] also provides a characterization and analysis of rainfall variability in the Mono river watershed complex. In the projection period, the four RCMs projected a similar distribution of precipitation in the study area as in the reference period, except the SMH model, for which a small change in the distribution of precipitation is displayed, i.e., a slightly reduced number of rain days beyond 10 mm. This close relationship between the distributions from simulations and observations is consistent In the projection period, the four RCMs projected a similar distribution of precipitation in the study area as in the reference period, except the SMH model, for which a small change in the distribution of precipitation is displayed, i.e., a slightly reduced number of rain days beyond 10 mm. This close relationship between the distributions from simulations and observations is consistent with previous findings [3,21,27,46,65,66] and confirms that the RCMs satisfactorily reproduced rainfall at a local scale as stated by Amoussou [2] and other works [4,6,67].
A detailed analysis was performed to overcome the bias challenges regarding the number of rainy days and to better describe the future changes in rainfall distribution. We computed the deciles of the observed daily rainfall and examined for each RCM the contribution of these 10 rainfall classes to the total rainfall as shown by López-Moreno and Beniston [27] in their study on the evolution of rainfall intensity in the Pyrenees. The computation was made for both the reference period (1988-2010) and the projection period (2028-2050). This approach helps to quantify the evolution of the influence of precipitation belonging to different classes of the frequency distribution [27,[68][69][70][71]. Figure 10 displays for the Mono river basin at the outlet of Nangbéto the rainfall distribution from the four RCMs according to the thresholds of the 10 quantiles as obtained from the observations. The threshold values are provided at the bottom of Figure 10. with previous findings [3,21,27,46,65,66] and confirms that the RCMs satisfactorily reproduced rainfall at a local scale as stated by Amoussou [2] and other works [4,6,67]. A detailed analysis was performed to overcome the bias challenges regarding the number of rainy days and to better describe the future changes in rainfall distribution. We computed the deciles of the observed daily rainfall and examined for each RCM the contribution of these 10 rainfall classes to the total rainfall as shown by López-Moreno and Beniston [27] in their study on the evolution of rainfall intensity in the Pyrenees. The computation was made for both the reference period (1988-2010) and the projection period (2028-2050). This approach helps to quantify the evolution of the influence of precipitation belonging to different classes of the frequency distribution [27,[68][69][70][71]. Figure 10 displays for the Mono river basin at the outlet of Nangbéto the rainfall distribution from the four RCMs according to the thresholds of the 10 quantiles as obtained from the observations. The threshold values are provided at the bottom of Figure 10.    Figure 10 indicates that apart from SMH, all RCMs underestimate the distribution of rainfall quantiles, especially the contribution of low rainfall (first four quantiles). This reveals an underestimation of the frequency of low rainfall. For the fifth decile to the eighth decile, the MPI and SMH simulations are close to the observed values. Meanwhile, both ICT and METO overestimate the rainfall in the basin, but this overestimation is stronger with METO. ICT and METO also underestimate the heavy rains (Q9 and Q10), unlike SMH and MPI, which simulate an increase of the daily rainfall (Q10) in the Mono river basin. Between the projection period (2028-2050) and the reference period (1988-2010), the analysis shows a relatively stable distribution of rainfall for most rainfall quantiles from MPI and ICT. For SMH, a strong decrease is found in the contribution of rainfall from the fifth quantile. Meanwhile, an increase of the contribution of the high rainfall intensities to the total rainfall is expected with METO.

Analysis of Extreme Precipitation under Present Day and Future Climates
The hydrological and societal importance of high-intensity rainfall coupled with the disparities between the RCMs regarding the reproduction of the 10th decile of daily rainfall (Q10 > 15 mm) led us to investigate the 99th percentile of daily rainfall in order to better understand the dynamics of heavy rainfall in the Mono river basin. We also analyzed the 5-day cumulative rainfall that likely contributes to flood generation downstream of the basin. Recurrent flood events were experienced in recent years in the Mono river basin [40] because of the high risk of soil saturation that induces a fast runoff. Table 3 displays the 99th percentile of precipitation for the reference and projection periods (1988-2010 and 2028-2050). The values in Table 3 show that in the reference period, about 30 mm/day of rainfall is equivalent to the 99th percentile in this basin, i.e., the highest 1% of the precipitation days according to the historical data. This recorded heavy rainfall intensity followed or preceded by rainy days could trigger catastrophic floods in the Mono river basin.
ICT simulation of heavy rainfall is the closest to observations. METO underestimated the heavy rainfall whereas SMH and MPI largely overestimated these rainfall extremes. A 99th percentile greater than 75 mm/day is too high for a catchment of 15,680 km 2 . However, this is not excluded because the rainfall of a year can precipitate in almost a day, as was the case in Ouagadougou (Burkina-Faso) in 2009 with nearly 290 mm.
In simulation as in projection, the ICT model simulates the heavy rains, but MPI and SMH overestimate them while the low rains are recorded by METO. With the reference data used, the MPI model especially simulates the very heavy rainfall the best, as well as projection periods. In both the reference and projection periods, the evolution is almost identical whatever the climate model [72,73].
In order to comprehend the seasonal distribution of these heavy rains, we also computed per month the 99th percentile of daily rainfall (Figure 11). The distribution of the 99th percentiles of daily rainfall shows that all RCMs underestimated the heavy rainfall from November to February (dry season). However, MPI overestimated the heavy rainfall in February and November while SMH overestimated only that of November in the dry season. The high rainfall observed in December and January is only sporadic because the occurrence of rainfall events during the dry season is extremely rare. SMH and MPI also overestimated the rainfall intensity in the core rainy season months (a maximum of 65 mm in September for MPI). This hinders a realistic estimation of the hydrological risks.
In the projection period, we will focus on ICT simulations and to a lesser extent on METO simulations as these two climate models stand out as being the most realistic for simulating heavy The distribution of the 99th percentiles of daily rainfall shows that all RCMs underestimated the heavy rainfall from November to February (dry season). However, MPI overestimated the heavy rainfall in February and November while SMH overestimated only that of November in the dry season. The high rainfall observed in December and January is only sporadic because the occurrence of rainfall events during the dry season is extremely rare. SMH and MPI also overestimated the rainfall intensity in the core rainy season months (a maximum of 65 mm in September for MPI). This hinders a realistic estimation of the hydrological risks.
In the projection period, we will focus on ICT simulations and to a lesser extent on METO simulations as these two climate models stand out as being the most realistic for simulating heavy rainfall events in the basin. For these two climate models, the heavy rainfall intensities obtained in the projection period are close to those of the reference period even if an increase is noticeable in many months. ICT shows a sharp one-off increase of rainfall intensity in March. However, SMH and MPI predicted a slight decrease in the intensity of the heavy rainfall events from February to August and from August to October, respectively, but these two climate models are the least realistic for the simulation of the heavy rainfall intensities.
The monthly variation of the 99th percentile rainfall maxima shows that the most intense heavy rainfall events occur at the core of the rainy season. However, the so-called dry months sometimes record strong rainfall. Therefore, an analysis of the return periods of the heavy rainfall events can help to better understand the distribution and occurrence of the extreme events in the Mono river basin, and the occurrence of high streamflow in the basin upstream of the dam, and the resulting floods. This result agrees with the study of Kodja [74] in tthe Ouémé river basin, which is located in the same climatic zone as our study basin.

Analysis of the Return Periods of Rainfall Extreme Using a GEV Model
The annual maxima of daily rainfall in the Mono river basin upstream of the outlet of Nangbéto (observed and simulated values) were adjusted to a generalized extreme value (GEV) distribution. The values of the three parameters, namely k, sigma (ς), and mu (µ)) of the GEV distribution (Figure 12), show a similarity between the observed and simulated data for the sigma and mu parameters but strong differences for the k parameter, where K refers to the shape parameter, sigma the scale parameter, and mu stands for the location parameter [75]. Nevertheless, these three parameters were always positive (for observations and simulations in the reference and projection periods) (Figure 12), thus pointing to a Frechet distribution (extreme value distribution type II).
Water 2020, 12, x FOR PEER REVIEW 19 of 27 rainfall events in the basin. For these two climate models, the heavy rainfall intensities obtained in the projection period are close to those of the reference period even if an increase is noticeable in many months. ICT shows a sharp one-off increase of rainfall intensity in March. However, SMH and MPI predicted a slight decrease in the intensity of the heavy rainfall events from February to August and from August to October, respectively, but these two climate models are the least realistic for the simulation of the heavy rainfall intensities. The monthly variation of the 99th percentile rainfall maxima shows that the most intense heavy rainfall events occur at the core of the rainy season. However, the so-called dry months sometimes record strong rainfall. Therefore, an analysis of the return periods of the heavy rainfall events can help to better understand the distribution and occurrence of the extreme events in the Mono river basin, and the occurrence of high streamflow in the basin upstream of the dam, and the resulting floods. This result agrees with the study of Kodja [74] in tthe Ouémé river basin, which is located in the same climatic zone as our study basin.

Analysis of the Return Periods of Rainfall Extreme Using a GEV Model
The annual maxima of daily rainfall in the Mono river basin upstream of the outlet of Nangbéto (observed and simulated values) were adjusted to a generalized extreme value (GEV) distribution. The values of the three parameters, namely k, sigma ( ), and mu (µ)) of the GEV distribution ( Figure  12), show a similarity between the observed and simulated data for the sigma and mu parameters but strong differences for the k parameter, where K refers to the shape parameter, sigma the scale parameter, and mu stands for the location parameter [75]. Nevertheless, these three parameters were always positive (for observations and simulations in the reference and projection periods) (Figure 12), thus pointing to a Frechet distribution (extreme value distribution type II). Furthermore, when k > 0, the GEV distribution is the type II, or Frechet, extreme value distribution. The variance of GEV is not finite when k ≥ 1/2. The GEV distribution has positive density only for values of X, such that k × (X-mu)/sigma > −1 [75]. Figure 13 shows the evolution of the annual maxima of observed daily rainfall. The analysis reveals a strong increase in the annual maximum of daily precipitation that is statistically significant Furthermore, when k > 0, the GEV distribution is the type II, or Frechet, extreme value distribution. The variance of GEV is not finite when k ≥ 1/2. The GEV distribution has positive density only for values of X, such that k × (X-mu)/sigma > −1 [75]. Figure 13 shows the evolution of the annual maxima of observed daily rainfall. The analysis reveals a strong increase in the annual maximum of daily precipitation that is statistically significant at the 95% confidence level. The frequency of high streamflow is not increasing, and this could be explained by the increase in temperature, which induces more evapotranspiration and thus compensates for the impact on high streamflow. This has also been proved by Ntajal et al. [49] downstream of the Mono river basin. However, an in-depth analysis is still needed to prove this assumption because, in principle, very high streamflow is sometimes generated from rapid runoff in response to water excess in soil with a low retention capacity [41,71,76]. This occurs in a very short time period, during which the increase of the evaporative demand is not really a constraint.
Water 2020, 12, x FOR PEER REVIEW 20 of 27 at the 95% confidence level. The frequency of high streamflow is not increasing, and this could be explained by the increase in temperature, which induces more evapotranspiration and thus compensates for the impact on high streamflow. This has also been proved by Ntajal et al. [49] downstream of the Mono river basin. However, an in-depth analysis is still needed to prove this assumption because, in principle, very high streamflow is sometimes generated from rapid runoff in response to water excess in soil with a low retention capacity [41,71,76]. This occurs in a very short time period, during which the increase of the evaporative demand is not really a constraint. This result is consistent with the evolution of discharge maxima flowing into the Nangbéto dam from 1988 to 2010, which pointed to an increase of the annual maxima, which is significant at p = 0.0185 according to a Mann-Kendall test [41,77,78].    This result is consistent with the evolution of discharge maxima flowing into the Nangbéto dam from 1988 to 2010, which pointed to an increase of the annual maxima, which is significant at p = 0.0185 according to a Mann-Kendall test [41,77,78]. Figure 14 compares the distributions of the quantiles of daily rainfall maxima observed over 23 years (1988-2010) to those adjusted by the GEV. at the 95% confidence level. The frequency of high streamflow is not increasing, and this could be explained by the increase in temperature, which induces more evapotranspiration and thus compensates for the impact on high streamflow. This has also been proved by Ntajal et al. [49] downstream of the Mono river basin. However, an in-depth analysis is still needed to prove this assumption because, in principle, very high streamflow is sometimes generated from rapid runoff in response to water excess in soil with a low retention capacity [41,71,76]. This occurs in a very short time period, during which the increase of the evaporative demand is not really a constraint. This result is consistent with the evolution of discharge maxima flowing into the Nangbéto dam from 1988 to 2010, which pointed to an increase of the annual maxima, which is significant at p = 0.0185 according to a Mann-Kendall test [41,77,78].    The annual maxima of daily rainfall are well adjusted to the GEV distribution. The analysis of the return periods of the extreme rainfall events, as displayed in Figure 15, shows that the rainfall events of an intensity higher than 40 mm are likely to occur at a recurrence interval of about 5 years.
Water 2020, 12, x FOR PEER REVIEW 21 of 27 The annual maxima of daily rainfall are well adjusted to the GEV distribution. The analysis of the return periods of the extreme rainfall events, as displayed in Figure 15, shows that the rainfall events of an intensity higher than 40 mm are likely to occur at a recurrence interval of about 5 years. Figure 15. Estimation of the return periods of the annual maxima of rainfall events (mm). The circle represents the observed values over 23 years (1988-2010) and the blue curves refer to the confidence level.
The uncertainty of the rainfall events occurring in a 50-year recurrence interval is higher because of the short length of the available time series and the trend shown in Figure 13. Table 4 presents the return periods for the annual maxima of the daily rainfall simulated with ICT, which stands out as the best of the four models for the reproduction of intense rainfall. The rarest and thus most intense rainfall events remain poorly simulated, as illustrated by the divergence of the 10-day recurrence values (45 mm in the observations, 77 mm in the model). Nevertheless, the projections show a strong increase in rainfall intensities for all return periods described in Table 4. The daily precipitation increases with the return period in the basin, which is in accordance with Ntajal at al.'s [49] and Kodja's [74] studies in the Mono and Ouémé basins, respectively. Ntajal et al. [49] compared the return period at the annual time scale and showed that GEV gives overestimation than the Goodrich exponential distribution model. This could also justify the overestimation of the model. This could induce high magnitude floods, which may pose a critical risk of flood inundation and thus cause damages with severe economic and societal consequences in the lower valley of the Mono river basin. Considering these return periods, the amounts of heavy rains projected are alarming. This confirms the increase in intense rainfall events predicted by the regional climate models. The uncertainty of the rainfall events occurring in a 50-year recurrence interval is higher because of the short length of the available time series and the trend shown in Figure 13. Table 4 presents the return periods for the annual maxima of the daily rainfall simulated with ICT, which stands out as the best of the four models for the reproduction of intense rainfall. The rarest and thus most intense rainfall events remain poorly simulated, as illustrated by the divergence of the 10-day recurrence values (45 mm in the observations, 77 mm in the model). Nevertheless, the projections show a strong increase in rainfall intensities for all return periods described in Table 4. The daily precipitation increases with the return period in the basin, which is in accordance with Ntajal et al.'s [49] and Kodja's [74] studies in the Mono and Ouémé basins, respectively. Ntajal et al. [49] compared the return period at the annual time scale and showed that GEV gives overestimation than the Goodrich exponential distribution model. This could also justify the overestimation of the model. This could induce high magnitude floods, which may pose a critical risk of flood inundation and thus cause damages with severe economic and societal consequences in the lower valley of the Mono river basin. Considering these return periods, the amounts of heavy rains projected are alarming. This confirms the increase in intense rainfall events predicted by the regional climate models.

Conclusions
The comparative analysis of different regional climate models simulating air temperature and precipitation in the Mono river basin and the West African sub-region leads to three major conclusions regarding their ability to reproduce the present-day climate.
First, apart from the few exceptions noted below, the realism of the simulated climate is more controlled by the regional climate model (RCM) than the driven global climate model (GCM). For the same forcing GCM (ECHAM5 or HadCM3 in this study), the precipitation fields are substantially different from one RCM to the others used for the dynamical downscaling. This strong model sensitivity may arise from the RCM structure (model approximations) or from the choice of parameterization techniques made to run each model (e.g., convective scheme).
Second, though the seasonal cycles of air temperature and precipitation are well reproduced by the RCMs, most of them still feature systematic errors related to the magnitude of these variations. In general, the selected RCMs underestimate air temperature in the range 1-3 • C for the Mono river basin. Two RCMs largely overestimate rainfall (by a factor of 1.5 to 2) during the rainy season. This overestimation is due to very high rainfall intensities that could not be compensated by the underestimation of the number of rainy days. The most realistic models regarding the reproduction of the actual seasonal cycle are those that simulate both the frequency of rainfall events and the average intensity of daily rainfall well (ICT and to a lesser extent METO). Finally, almost all RCMs introduced bias with respect to the spatial variation of the climatic variables being analyzed. The north-south gradients of precipitation between the Sahel and the Guinean regions are well reproduced by all RCMs, except METO, although the simulated values for precipitation are not always correct. There are also errors in the southern part of the Guinean region regarding the reproduction of the drought band diagonally stretching from southern Ghana to southern Benin (Dahomey Gap). It was not reproduced by the RCMs forced with ECHAM5 GCM at their boundaries (MPI and ICT).
Third, and this is implicitly derived from the conclusion above, the hierarchy of RCMs in terms of the quality of their simulations differs according to the climate variable considered. For example, MPI tends to offer the most accurate reproduction of air temperature (especially their spatial distribution) and ICT is the most realistic for the precipitation patterns. This conclusion still needs to be taken with caution because though ICT satisfactorily reproduces the rainfall regimes and the frequency distribution of the rainfall events, the spatial distribution of the average annual precipitation is better simulated by the SMH model. Indeed, ICT fails to reproduce the relative dryness of the Ghanaian-Togolese coast in the south of the Guinean region. As well, ICT still features bias regarding the reproduction of the return periods of the annual rainfall maxima, for which it overestimated the frequency despite its good performance for the reproduction of the distribution of daily rainfall. Thus, we conclude that the choice of an appropriate RCM should be based on the research objectives, especially for impact studies.
Though these conclusions call for caution in using climate model simulations to assess future climate change in the sub-region, some lessons can nevertheless be drawn. The four RCMs unanimously project a +1 to +1.5 • C temperature increase in the sub-region by 2028-2050 as compared to the 1988-2010 reference period. The projected trends for cumulated precipitation are null or very moderate and diverge from one RCM to the others. Only SMH projected a significant decrease of precipitation in the Mono river basin. The simulations of the intense rainfall events (99th percentiles) also display contrasted results, but the two regional climate models that better reproduce these extreme rainfall events in the reference period have also unanimously projected a significant increase in their intensity by 2028-2050.
For the intense rainfall events, the Mono river basin at the outlet of Nangbéto recorded between 1988 and 2010 a significant increase in the intensity of the maximum daily rainfall averaged over the watershed, and this increase may be a risk factor. The heaviest rainfall events (99th percentile) occur mainly during the months of July to September, with a peak in September. This fits well with the discharge peaks in the basin. Thus, the intense rainfall events recorded in the tropical domain of the basin could lead to probable risks of flood inundations in the lower valley of the Mono river basin due to an overflow of the river.
The ICT simulations were the closest to the observations in terms of the reproduction of the cumulated precipitation and frequency distribution of daily rainfall. ICT is thus more robust for the analysis of the intense rainfall events and thereby flood forecasting in the Mono river basin. The estimation of extreme rainfall events using ICT showed stability except for the rainfall events occurring in October-November, which slightly increase. The adjusted GEV model satisfactorily predicted the heavy rainfall events, but there are uncertainties beyond the 90th percentile. The adjustment of a Weibull-type GEV distribution is well suited here for the prediction of extreme rainfall events (shape parameter of the GEV distribution = 0), although it is heavy-tailed at the 99th percentile.
The annual maxima of daily precipitation are increasing. Therefore, the impact of heavy rainfall events on discharge will increase under climate change conditions. This knowledge is important to develop an early warning system for the Mono river basin for the sustainability of the ecosystems. This calls for further research to better understand the physical mechanisms underlying the changes in extreme precipitation events and thus improve climate model simulations for short-and medium-term use. Moreover, for future study, a bias correction could be done for those datasets before quantifying the magnitudes of changes due to their coarser resolution.