Assessment of Regional Climate Models over Côte D ' Ivoire and Analysis of Future Projections over West Africa

The ability of six Regional Climate Models (RCMs) used in AMMA-ENSEMBLES project is assessed over six meteorological stations in Côte d’Ivoire. The ensemble mean of the models is also used for the prediction of climate change over West Africa. The study focused on two periods: the period 1995-2005, the present-day simulations, is used to evaluate the skills of the models over the country and the years 2010-2013, for assessment of the future climate change scenario used. The results show that the skills of the models vary from one station to another and from one season to another. None of the models considered, presents an excellent performance over the entire country and in all the seasons. Generally, the ensemble mean of all the models presents better results when compared with the observation. These results suggest that the choice of any model for study over the country may depend on the focus of interest: intensity or variability of the rain and also on area of interest. The projection for 2020-2040, future climate change over West Africa shows that the Sahel exhibits a tendency to be drier while wetter Guinean coast is observed.


Introduction
West African Monsoon is a coupled atmospheric, oceanic and land system.These coupled systems include mesoscale convective systems (MCSs); modulated by the westward moving African Easterly Waves (AEW), The Mid-tropospheric African Easterly Jet and Inter-hemispheric Tropical Easterly Jet.The monsoon is also strongly dependent on regional and global Sea Surface Temperature (SST) variations [1] [2].Lamb [3] reports that the Atlantic SST anomaly influences the sub-Saharan summer rainfall.He further showed that there is negative SST anomaly during spring (April-May) and a positive one during the summer (June-September) which is associated with dry year and vice versa.Apart from SST anomaly many studies has related rainfall to other changes in surface conditions.For instance, many studies have been linked soil moisture to precipitation over West Africa [4]- [6].It has a positive feedback with the precipitation and may influence the position of the storm system [7] [8].The Inter-tropical Discontinuity (ITD) has been mostly pointed out for its influence on the rain distribution over this region mainly on the northward shift and the southward retreat of the WAM [9] [10].
Regional Climate Models (RCMs) have been used in acceptable manners to understand the West African Monsoon (WAM) at different spatial and temporal scales [11]- [18].Recently, in the framework of the West African Monsoon Modeling and Evaluating (WAMME) project, [16] observed a relative difference between the simulations made at the Universities of Cocody (Côte d'Ivoire) and Abdou Moumouni (Niger) while executing an inter-comparison between five RCMs.This bias was attributed to the different kind of Sea Surface Temperatures (SST) used for the forcing of different WAM simulations.During The African Monsoon Multidisciplinary Analysis (AMMA) program, many efforts have been made for better understanding the WAM through a long-term collaboration to study the climatic and environmental feedback mechanisms involved in the African monsoon, and in some of its consequences on society and human health.
The major setback in the study of WAM is the availability of continuous and reliable datasets.The AMMA-ENSEMBLES project funded by the European Commission's 6th Framework Program is an opportunity to have deeper insight into the behavior of West African Monsoon system.This platform had already been used to undertake study over West Africa.For instance, [19] focusing in an interannual variability of rainfall over Sahel report that main of the models used provide a good estimation of the rainfall variability because of their skills in the representation of the WAM features.Other efforts have been done using data from similar project: The Coordinated Regional Downscaling Experiment (CORDEX) program recently established by the World Climate Research Program (WCRP).Among them, [20] assessed the ability of 10 Regional Climate Models (RCM) to replicates the observed features of WAM.They found out the capability of these models in the replication of temporal and spatial distribution of the precipitation and temperature over West Africa but with biases.These biases vary from one model to another.They also reported like [21] that using the ensemble mean of these models enhanced the representation of these features.This result was corroborated by the work of [22], focusing only on precipitation [16] [23].More recently, [21], discovered significant differences firstly between the observation and the models and secondly among the RCMs themselves in term of daily rain amount, intensity, extremes and the duration of rainfall events inside the WAM.According to them, this underlined the problematic of the sweeping uncertainties in the simulation of daily precipitation features.These results reinforced the belief that precipitation is the most difficult parameter to simulate accurately in RCMs.
This study will assess the performance of six models used in the AMMA-ENSEMBLES project in their basic representation of rainfall.In analysis of the result, we will consider two time windows.The first one, 1990 to 2005 which is related to the historical and current period, and the short and second one is likely an opportunity to evaluate these models in the projection of the future.So that this result may indicate a degree of reliability in the use of the dataset in Côte d'Ivoire particularly and West Africa in general.Secondly we will also perform a projection over the period 2020-2040 taking 1990-2005 period as the reference.

Domains, Data and Methods
Under the AMMA-ENSEMBLES project, an ensemble of simulations [24] has been made over Africa domain.For these simulations, all of the RCMs are driven for the period 1990-2007 (first window) by the same lateral boundary conditions taken from the European Centre for Medium-Range Weather Forecasts reanalysis ERAinterim ( [25] [26]).For the projection 2010-2013 (second window), the simulations have been done under the scenario A1B and the lateral boundary conditions are provided by ECHAM5-r3 or HadCM3Q0.The objective is to build an international coordinated framework for generating improved regional climate change projections worldwide.In this work, we will firstly assess the basic performance over Côte d'Ivoire of six models (METO-HC_HadRNM3 (HadRM), ICTP-REGCM3 (RegCM3), METNO-HIRHAM, SMHI-RCA (RCA), DMI-HIRHAM5, KNMI-RACMO2 (RACMO)).Secondly we will make a projection over West Africa.All the simulations had been made over the same domain and havethe same horizontal resolution (0.44˚ × 0.44˚, ~50 km).Simulation outputs are freely available from http://ensemblesrt3.dmi.dk/.Table 1 gives an overview of the RCMs configurations for these simulations [27].
Two time windows are considered in this work.The first time window was between 1990 and 2005, the lateral boundary was forced by ERA-interim data and while for second time slice (2010-2013), was forced by ECHAM5-r3 or HadCM3Q0 data (see Table 1).In the first windows, the ability of the models to reproduce the past and current rainfall pattern is assessed while the later one is to assess the performance of those models in the projection of the future scenarios.The models outputs are compared to six (6) weather stations data from the National Meteorology Office of Côte d'Ivoire.Non-availability of continuous dataset due to socio-economic crises in Côte d'Ivoire limited the studies to six meteorological stations, of Abidjan, Bondoukou, Daloa, Dimbokro, Gagnoa and Tabou have been used.This lack of data impresses to focus only in the southern climatic zone (Abidjan, Tabou) and the center climatic zone (Bondoukou, Dimbokro, Gagnoa and Daloa) part of the country instead of the whole area as described in [28] [29].In fact these works suggested three climatic zones (south, center, north) in Côte d'Ivoire which may be related to the vegetation cover; forest in the south, wooded savanna in the center and grassy savanna in the north.Figure 1 gives an overview of the climatic zone of Côte d'Ivoire adapted from [28] [29] and the location of the stations used over the country.
Taylor diagram is used as statistical tool in order to check how close are the model estimation to the observation [30], root mean square (RMS) is also used to test the deviation of model results from observations.The standard deviation (mm of rain) indicates the amplitude of variation.So that a simulation that is in good agreement with the observed will be located near the point marked "obs" (observation) on the x-axis.This simulation must have relative high correlation and low RSM errors (RMSE).It should have furthermore the correct standard deviation which will indicate that the pattern variation has the right amplitude.This work focuses specially on the monthly evolution, the seasonal cycle and also the estimation of the maximum of rain amount.
The regional analysis over West Africa is made with a special focus over three basins; the Volta, Senegal and Bandama basins.The Volta basin covers an area of 400,000 km 2 , stretches from 5˚30'N in Ghana to 14˚39'N in Mali and its widest stretch is from 5˚30'W to 2˚E.It is spread over six West African countries (43% in Burkina Faso, 42% in Ghana, and 15% in Togo, Benin, Côte d'Ivoire and Mali).The Volta basin is divided into two parts (north and south) according to the local climate.The Senegal basin covers around 300,000 km 2 .The Senegal River is 1800 kilometers long.This basin stretches from 8˚N in Guinea to 14˚39'N in Mauritania and its widest stretch is from 18˚W to 5˚W.Its riparian countries are Guinea, Mali, Mauritania, and Senegal.The main part of the basin is located in the Sahelian and desert zone.The Bandama River is the longest river in Côte d'Ivoire with a length around 800 kilometers.In this study the term Bandama refers to all the Rivers (Bandama, Comoe, Sassandra and Cavally) in Côte d'Ivoire.Then we will have a general idea of the basins in the country.
In the following section the term ensemble refers to the average of the estimation of the different models used.In order to check the performance of each model over Côte d'Ivoire, the considered six stations were averaged for the observation, the model outputs and the ensemble mean.This may represent a summary of the performance of each model and the ensemble over the entire domain.Figure 3 shows that DMI-HIRHAM model result is quite different to the other models.In fact DMI-HIRHAM standard deviation is closest to the observed while the other models and the ensemble present less standard deviation than the observation.This model DMI-HIRHAM has a high correlation coefficient (~0.6) however the RMSE is high (110) that means a simulated amount is different from the observation.The models RegCM3, RCA, METNO-HIRHAM, HadRM have RMSE value ~95, and correlation coefficient between 0.4 and 0.6 and the difference between their standard deviations and the observed is about 30.The model RACMO present the smallest value of RMSE (~80) and the highest correlation coefficient (~0.63) however the difference between its standard deviation and the observation one is about 40.On the following sections, the annual cycle of the precipitation and the time series at each selected station (Figure 4) and its associate statistics in Taylor Diagram (Figure 5) are presented and discussed.This is to further analyze the model for their capability in replicating the mean variation of the precipitation over the year and how close the models representation are to the observation.

Simulation of Annual Cycle of Rainfall
Figure 4 shows the seasonal variability of precipitation in the six stations considered.The bimodal nature of the rainy season is clearly shown by all the models even though they did not accurately replicate the magnitude.DMI-HIRHAM presents the closest variation compared to the observation in Abidjan, Daloa and Gagnoa.In Abidjan, DMI-HIRHAM replicates the variation in the first peak of rainfall season and a wet bias in the second peak; this was corroborated by the Taylor diagram presented (Figure 5).In that station, all other models were at various degrees of dry biases.In Daloa, Dimbokro and Gagnoa, DMI-HIRHAM demonstrates highest wet bias while other models were at varying degree of accuracy as compared to the observation.The result of the seasonal variability corroborates the result of Taylor plot presented in earlier section.The ensemble mean shows a realistic variation of precipitation but pain to replication of the magnitude.In Abidjan, Dimbokro, and Gagnoa, ensemble replicates the variation while it underestimates the first maximum.There is a lag in the indication of the month of maxima rainfall between the ensemble mean and the observation.This is clearly seen at the station of Daloa where this ensemble indicates the first maximum in June while the observation shows the same event in April-June.At the station of Bondoukou, the second maximum of rain is in October in the observation while the ensemble simulates it one month earlier.Overall the second rainy season is relatively well captured in comparison to the first season.In all the stations, the ensemble shows high value of correlation coefficient, low value of the centered root mean square error and has confident standard deviation.For example at Bondoukou and Daloa the correlation coefficients are above 0.95, with a standard deviation similar to the observed.Low RMSE makes ensemble closest to the observation in the simulation of the annual cycle as compared to the result of individual models.

Seasonal Analysis
The season May-June corresponds to the intensification of the first and major rainy season in the southern part of the country followed by a shift toward the northern part and beyond [9].In fact, [9] focusing in the monsoon dynamics found out that July-August is the Little Dry Season (LDS) in the southern part and the period of intensification of the monsoon rain in the northern part.September-October correspond to the retreat of the rain belt toward the south.It is then the second and short rainy season in southern part of the country [28] [31].Therefore a good performance of a model at any of these seasons maybe very important in the seasonal forecast of climate in Côte d'Ivoire and generally over West Africa since rainfall is one of the most important indicators of the monsoon variability [32].
Interannual evolution of the seasons; May-June, July-August, and September-October from 1990 to 2005 is analyzed over the selected six stations in Côte d'Ivoire.Time series computed on these seasons show that all the models pain to simulate the interannual evolution of these seasons (figures not shown).The amount of precipitation is somehow not well estimated, the computed correlation coefficients between models and observation are very low (generally <0.4) with a high RMSE and a standard deviation different for the observation one.For this reason, the ensemble which is the average of all the models, presents less satisfactory results.However for the period May-June the Taylor diagram shows that at Bondoukou, there is a good representation of the season by the model METNO-HIRHAM with a high correlation coefficient (~0.75) and a relative low RMSE.In September-October only the ensemble shows a result sensibly close to the observation especially at Bondoukou while the models HadRM, RCA and METNO-HIRHAM seem outliers.In July-August at the station of Dimbokro, the model RACMO reproduces fairly well the pattern of the rainfall; it has a standard deviation close to the observation with a correlation coefficient of about 0.6 and relatively low RMSE.In spite of these lackluster performances, the simulation of the mean values of these peaks is reasonably better estimated.In Figure 6 and Figure 7 we computed the mean over the years of the amount of precipitation at the different peaks of the rainy seasons and at the different stations.Figure 6 shows the bias in percentage of each model rain amount simulated compared to the observation for the first peak.The bias is computed in June which, represents the period when this peak occurs at Abidjan, Dimbokro, Gagnoa and Tabou while for Bondoukou and Daloa it occurs in May.Almost all models underestimate the amount of rainfall during the first maximum at all the stations.For example, at Dimbokro, Gagnoa and Tabou all the models show negative bias except DMI-HIRHAM that shows positive bias in Abidjan and Bondoukou.DMI-HIRHAM shows along the coastal areas, the best estimation of the first maximum even though the bias is ~−25% at the station of Tabou.Indeed it exhibits 1.86% as bias in Abidjan and the bias of 25% at Tabou is the best.This underlines that DMI-HIRHAM may perform well in the littoral.For the other stations, the models pain to proper estimate the first peak accordingly the ensemble exhibits bias beyond 26% except at Bondoukou (18%) and Daloa (5%).
Figure 7 shows the bias in percentage of the rain amount simulated compared to observation for the second rainfall peak.The bias is computed in October which is the period of this peak in Abidjan, Dimbokro, Gagnoa and Tabou while it occurs in September for Daloa and Bondoukou.
Contrary to the result of the first maxima, there is a good estimation of the second peak.For example, at Abidjan, RegCM3 and RACMO are in good agreement with the observation with respective biases of −0.66% and 2.28%.DMI-HIRHAM overestimates the rainfall at that peak, the rests perform well.DMI-HIRHAM shows a bias of ~100% at Abidjan and overestimates in other stations except in Bondoukou, where it underestimates the amount of rainfall (−23.81%).At Bondoukou, METNO-HIRHAM performed best with a bias of 0.92%.For this station, three other models are under 20% of bias (RegCM3, HadRM and RACMO).In summary, there is a good estimation by the models except DMI-HIRHAM of the mean rainfall amount during the second peak of the rainy season at the different stations.

Projection Assessment (2010-2013)
In this section, the accuracy of the scenario used to run the projection of the near future is assessed.Since the different simulations have been run before 2010, this period 2010-2013 is somewhat a great occasion to ascertain the accuracy of the projection over those years before using it for prediction.

Inter-Seasonal Variability
Figure 8 shows the time series of the annual cycle of rainfall over the period 2010-2013.The models did not replicate well the inter-seasonal pattern of rainfall.Contrary to bimodal structure of rainfall especially along the coast and the center regions [28], all the models in the present simulation show a unimodal structure.This result may be due to an unsuitable climate change scenario.Indeed this lower performance may be related to error in the input GCMs (ECHAM5-r3, HadCM3Q0) as the models replicated correctly the bimodal structure when driven by ERA-Interim presented in earlier section.Observations show bimodal structure of rainfall in five of the six stations, in Abidjan, Only RegCM replicated this nature of the seasonal cycle while other models show unimodal rainfall in this station contrary to Observation.The models also show dry bias in the first peak of rainfall season and wet bias in October.In all other stations, there was a dry bias in Tabou and in Bondoukou while wet bias in Daloa and Gagnoa.To some reasonable extent however, RegCM and METNO-HIRHAM replicated rainfall structure in Abidjan and Daloa and not well in other areas.The ensemble only performed fairly well at the station Daloa.

Rainfall Monthly Analysis
Monthly evolution of the rain at each station is presented in Figure 9.In the coastal stations (Abidjan and Tabou), the ensemble mean of all the models did not accurately pick the daily cycle of the rainfall.In June 2010, the ensemble shows an amount of about 190 mm while it is observed to above 400 mm in observations.The ensemble, did not only underestimated the amount of rain but also shift the position of the rainfall maxima, there was lag or lead of one to two months.The more realistic estimation of rain in general is achieved at the station of Bondoukou, Dimbokro and Gagnoa in term of amount and variability even it still have a one-month shift which corroborated by Taylor diagram (figure not shown).At Bondoukou, ensemble did not show the LDS centered in August in the years 2010, 2012 and 2013.For the station of Dimbokro, the ensemble showed unimodal structure of the seasonal cycle over the entire period.Individually, along the littoral, the best estimation is given by METNO-HIRHAM (Figure 9).This result is also supported by the associated Taylor diagram (figure not  In order to have a clear idea of the ability of each model in the representation of the seasonal peaks we extracted and plotted these different peaks.The selection of the month is based on the correspondent seasonal cycle.Thus the month of June is selected for first peak for Abidjan, Bondoukou, Dimbokro and Tabou while May for Daloa and Gagnoa.In the second peak, the month of September is only selected for Tabou and while October for the rest (Abidjan, Bondoukou, Daloa, Dimbokro, Gagnoa and Tabou).
Figure 10 shows the bias of the simulated rain amount in first maxima for the period 2010-2013.It is ob-served that the models underestimated the rain amount by about 20% to 80%.In the coastal area, all the models underestimated the amount of rain and only RCA presented a bias under 30% in Abidjan.Because of the wide gap between the models estimations and the observations, the computed ensemble in the first peak shows a bias above 30% excepting at Bondoukou (the most northern of the stations used) and Dimbokro (in the center of the country).In this center region, (Dimbokro) HadRM underestimated rainfall by about 8.63% which is the best estimation, while for the same location it was RCA that performed best for 1990-2005 simulations.At Daloa, DMI-HIRHAM showed a bias under 20% while in Abidjan and Tabou none of the models showed a bias less than 20%; the lowest is provided by RCA with −20.98% of bias.The best estimations by ensemble are observed at Bondoukou and Dimbokro with a bias of −6.25% and −7.46% respectively.Thus ensemble may be used to have realistic and good projection of the first maximum in the coastal region of the country.The bias of the simulated rain amount for the second maximum is shown in Figure 11 for the period 2010-  2013.At Tabou, all the models underestimated the maxima of precipitation like for the first peak with a bias of above 50%.Unlike the simulation of the period 1990-2005 presented in the earlier section, none of the models performed well in the estimation of this event.However for Abidjan, there is no clear tendency for all; some models underestimated while others overestimate.The best estimation is given by RCA (−0.02% of bias) followed by RegCM3 (9% of bias).The computed ensemble displayed a good estimation of this event with 3% of bias unlike in the first maximum.At Bondoukou the best estimation is given by the ensemble, when among the models, HadRM (7%) presented the best estimation of the second peak followed by RACMO (9.5%).METNO-HIRHAM (3% of bias) gave the closest estimation of the second peak preceded by the computed ensemble (0.04% of bias) and followed by the RACMO (8%).At Daloa in the center of the country, HadRM performed best with 2% of bias.In this part ensemble is able to represent the second peak but the bias is little bit high (28%).In first peak of rainfall season, all the models failed in the estimation of the rainfall amount, while RegCM (−8.70), and RCA (−0.02%) performed well in the second peak.Ensemble showed satisfactory result except at Tabou which high bias of about −70.15%.This may indicate that in Tabou, there are some local physical processes that models failed to take into account.In summary, the ensemble mean of all the models outperformed the skill of individual models and gives better estimation of the precipitation amount.It is then considered more suitable for the projection study.This part of the work is the projection of the West African climate over 2020-2040.Because of the result in the previous sections, we used the ensemble mean of the models.In fact the previous sessions of this work emphasized the improvement of representation of rainfall pattern when using ensemble.In addition, [19] [20] [22] [23] independently showed that ensemble mean outperformed the skill of almost all the individual models.Attention is given to the seasons June-September (JJAS) which is the rainy season in the Sahel, May-July (MJJ) and September-November (SON).MJJ and SON is the first and second rainy season respectively in the Guinean coast area.Figure 12 shows the difference between mean rainfall of 1990-2005 and 2020-2040 for the seasons JJAS, MJJ and SON.There is for all the considered seasons, a decrease in the amount of rainfall over the north and an increase in the south in 2020-2040 compare to the reference period (1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005).During the season JJAS, the neutral (stationary) area (limit between area of decrease and increase) is located around 10˚N.The stationary zone is located around 5˚N during MJJ and SON.It's seemed that this neutral area follows the mean position of the ITD.[9] [34] showed that the mean positions of the ITD are located around these latitudes respectively during these periods.The maximum of precipitation between May and beginning of June is over the Gulf of Guinea (0˚N -6˚N).This rain belt shifts around 10˚N between mid-June and early July then move back to Gulf of Guinea in September.Furthermore the northern part (Sahelian zone) is well known to have less precipitation than the southern part.Thus this led to conclude that the area which has less precipitation will be drier and the zone which has more precipitation will be wetter.This suggests an increase climatic stresses for the local population who depend exclusively on monsoon rainfall for livelihood.There may also be an increase in the probability of drought and flood events in the north and the south respectively.This is in agreement with the work done of [33] that supposed a prevailing drying tendency in sub-Saharan Africa.However there are divergent results over the coastal region.Indeed this southern zone may have more precipitation contrary to [33] who report likely deficit of precipitation over this zone.

Changes in Precipitation Interannual Variability in Regard to 1990-2005 of Three Different Basins over West Africa
The effect of changing climate on three basins over the region is presented in this section.These basins are Volta basin (south and north), Senegal basin and the Bandama basin.The first two basins were selected because of their impact on socio-economic activities of the surrounding countries.The third one gives idea of evolution of rivers over Côte d'Ivoire and their associated impacts.This section focuses on the yearly evolution of precipitation over the season May-October (MJJASO) which is rainfall season in these areas.Figure 13 shows an expected that the basins located in the south received more precipitation than those in the north.For instance, the Bandama and the Volta south basins have their mean value around 900 mm for the period 1990-2005 and 1100 mm over the period 2020-2040 while the northern basins obtained less than 900 mm of rain.In fact, the northern

Changes in Precipitation Annual Cycle over Volta, Senegal and Bandama Basins
The annual cycle of the period 2020-2040 has been compared to 1990-2005 in order to have an insight into the inter-seasonal variability in the near future (Figure 14).It was observed that for the basins located in the Sahel, the projection present the same structure while over the Guinean coast this projection shows a modified clima- tology.The simulation shows that in the period 2020-2040, there will be one rainy season for both Sahelian and coastal area regions, which implies the disappearance of LDS and bimodal rainfall regime in the Guinea coast as observed in 1990-2005 (Figure 14(a)).Two possibilities can be deduced from this result.The first one suggests that there may be modification of the climate over the Guinean coast which is less plausible and it may be probably inaccuracy of the projection of the models used in this study.The period 2010-2013 has been considered to check the first surmise which is related to the modification of variability of the annual cycle over theses basins.Thus the annual cycle of the rainfall from the ensemble mean is compared to the TRMM data (figure not shown).As we expected, this ensemble mean showed once again a unimodal rainfall pattern while TRMM exhibited a bimodal rainy season over the coastal area.Therefore it seems that the scenario consider to run the simulations was not too realistic to be able to capture the climatic features of the Guinean coast region.In agreement of this, many studies [33] [35] pointed out the inherent problems of Regional Climate Models (RCMs) on future precipitation change in West Africa.

Conclusions
The study consisted of the assessment of the ability of RCMs from AMMA-ENSEMBLES in their basic representation of the rain pattern.So the most suitable RCMs would be identified for climate study and policy maker in Côte d'Ivoire in particularly and generally over West Africa.The evaluation, first time period (1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005), consisted of time series and statistical (Taylor diagram) analyses over six stations while the second set (2010-  The analysis of the validation experiment indicates a performance of RCMs varying with the focus of interest.Some models are relatively good in retrieving the first rainfall peak, case of DMI-HIRHAM in the coastal area, while RegCM3, HadRNM3, RACMO, METNO-HIRHAM and KNMI_RACMO show good performance in reproducing the second peak.However, the computed ensemble mean presents a representation close to the observation.This emphasizes that the use of ensemble mean outperforms the skills of almost all the single RCM either in magnitude or spatial distribution.Thus, it illustrates the advantage of multi-models assessment of the WAM rainfall as presented in previous studies [20] [22] [33] [36].
The projection for period 2020-2040 over West Africa reveals decreasing amount of rainfall at the north of 10˚N and an increase at the south of 10˚N.This result suggests a rise of climatic stresses on local population, namely, an augmentation of probability of drought in the northern part and flood in the southern part.There is however limitation in the accuracy of this projection especially over the Guinea coast in the replication of the annual cycle.It seems that the scenario used in the simulations is not able to capture all the features of the variability of the rainfall in the West African coastal area.The elaboration of future scenarios that take into account the specificity of the southern part of West African climate might be necessary to make better and realistic projections.Furthermore, simulation of the climate at finer resolution may contribute to a better representation of physical processes of this region.

Figure 2 .
Figure 2. Comparison of monthly precipitation statistics (standard deviation, correlation coefficient, and root mean square) of AMMA-ENSEMBLES data each other and with observation over selected stations in Côte d'Ivoire.aconvergence between the statistics of RegCM3, RACMO, HadRM, METNO-HIRHAM, RCA at the stations of Abidjan and Daloa since they are confined in a small area.In all the models, DMI-HIRHAM precipitation result is the closest to the observation except at Bondoukou and Daloa.This suggests that DMI-HIRHAM has the best skill out of all the models in replicating rainfall with the correlation coefficient of about 0.6 even though the centered root mean square error (RMSE) is to some extent high.The plot of the monthly time series (not shown) for the period(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005) shows that DMI-HIRHAM generally underestimates the rainfall amount.The skill of the model is better over the coastal areas of Tabou and Abidjan than in other areas.The difference between the values of the standard deviation displayed by DMI-HIRHAM and the observation is very low (≤ 5) in Tabou and about 20 in Abidjan while the rest of the models are up to about 40.In this littoral area (Abidjan and Tabou), all the models except DMI-HIRHAM have a standard deviation around 80 while the standard deviation of the observation is about 140.This high difference expresses low amplitude of variation on the simulated rainfall compare to the observed.In addition, the high RMSE suggests over or under estimation of precipitation by models.DMI-HIRHAM shows a variation relatively close to the observed with high correlation coefficient, but it has a high RMSE.Beyond the coastal zone, the models display closer results to the observation.For instance, at Bondoukou while the standard deviation of the observation is 70, the standard deviation of the simulated is between 50 and 70 with a correlation coefficient above 0.6.In these stations the RMSE has its smallest values (less than 80) while along the coast the value is above 100.The ensemble mean presents the highest correlation coefficients of >0.7 except at Daloa.The variation of this ensemble seems not in good agreement with the observation.The difference between the standard deviation of the observation and ensemble is above 70.At Tabou and Abidjan, at the station of Gagnoa and Dimbokro this difference is ~30.The most right variation display with ensemble compare to the observation is obtained at the stations of Daloa and Bondoukou.In order to check the performance of each model over Côte d'Ivoire, the considered six stations were averaged for the observation, the model outputs and the ensemble mean.This may represent a summary of the performance of each model and the ensemble over the entire domain.Figure3shows that DMI-HIRHAM model result is quite different to the other models.In fact DMI-HIRHAM standard deviation is closest to the observed while the other models and the ensemble present less standard deviation than the observation.This model DMI-HIRHAM has a high correlation coefficient (~0.6) however the RMSE is high (110) that means a simulated amount is different from the observation.The models RegCM3, RCA, METNO-HIRHAM, HadRM have RMSE value ~95, and correlation coefficient between 0.4 and 0.6 and the difference between their standard deviations and the observed is about 30.The model RACMO present the smallest value of RMSE (~80) and the highest correlation coefficient (~0.63) however the difference between its standard deviation and the observation

Figure 3 .
Figure 3.Comparison of monthly precipitation statistics (standard deviation, correlation coefficient, and root mean square) of AMMA-ENSEMBLES data each other and with observation over all the entire stations combined in a single one.

Figure 4 .
Figure 4. Annual cycle of precipitation computed over the period 1990-2005.

Figure 5 .
Figure 5.Comparison of precipitation annual cycle statistics (standard deviation, correlation coefficient, and root mean square) of AMMA-ENSEMBLES data each other and with observation over selected stations in Côte d'Ivoire.

Figure 6 .
Figure 6.Bias of the mean rainfall amount computed on the month of the first peak of the rainy season over the period 1990 to 2005.

Figure 7 .
Figure 7. Bias of the mean rainfall amount computed on the month of the second peak of the rainy season over the period 1990 to 2005.

Figure 8 .
Figure 8. Monthly annual cycle of rainfall computed over the period 2010-2013.The observation data (blue) and the ensemble (orange) have been displayed in bold line.The models output as for them are displayed in fine line.

Figure 9 .
Figure 9. Monthly evolution of rainfall over the period 2010-2013.The observation data (blue) and the ensemble (orange) have been displayed in bold line.The models output as for them are displayed in fine line.shown), While RACMO and HadRM are in good agreement with the observation in Bondoukou.At Daloa, METNO-HIRHAM, DMI-HIRHAM and HadRM exhibitthe best simulations while only HadRM presents realistic estimation for Gagnoa and Dimbokro.In order to have a clear idea of the ability of each model in the representation of the seasonal peaks we extracted and plotted these different peaks.The selection of the month is based on the correspondent seasonal cycle.Thus the month of June is selected for first peak for Abidjan, Bondoukou, Dimbokro and Tabou while May for Daloa and Gagnoa.In the second peak, the month of September is only selected for Tabou and while October for the rest (Abidjan, Bondoukou, Daloa, Dimbokro, Gagnoa and Tabou).Figure10shows the bias of the simulated rain amount in first maxima for the period 2010-2013.It is ob-

Figure 10 .
Figure 10.Bias of the mean rainfall amount computed on the month of the first peak of the rainy season over the period 2010-2013.

Figure 11 .
Figure 11.Bias of the mean rainfall amount computed on the month of the second peak of the rainy season over the period 2010-2013.

4. 1 .
Analysis of Rainy Seasons over West Africa for the Period 1990-2005 and 2020-2040

Figure 12 .
Figure 12.Ensemble mean difference of the periods 1990-2005 and 2020-2040 computed over the seasons JJAS, MJJ and SON.part of Volta basin has a mean value around 800 mm for both 1990-2005 and 2020-2040 periods.Over the Senegal basin, the mean value fluctuated around 650 mm for the period 1990-2005 and 600 mm for the period 2020-2040.This result is in agreement with the result obtained in Section 4.1 which suggested an increase in the south (coastal area) and a decrease in the north in the amount of rain even the northern part of the Volta basin showed a steady evolution.

Figure 14 .
Figure 14.Annual cycle average rainfall computed for the periods 1990-2005 and 2040-2040 over the Volta, Senegal and Bandama basins.

Table 1 .
Summary of the different Regional Climate Models used in this work and their respective main physical schemes and latest references.