Intraseasonal Precipitation Variability over West Africa under 1.5 C and 2.0 C Global Warming Scenarios: Results from CORDEX RCMs

This study assessed the performance of 24 simulations, from five regional climate models (RCMs) participating in the Coordinated Regional Climate Downscaling Experiment (CORDEX), in representing spatiotemporal characteristics of precipitation over West Africa, compared to observations. The top five performing RCM simulations were used to assess future precipitation changes over West Africa, under 1.5 ◦C and 2.0 ◦C global warming levels (GWLs), following the representative concentration pathway (RCP) 8.5. The performance evaluation and future change assessment were done using a set of seven ‘descriptors’ of West African precipitation namely the simple precipitation intensity index (SDII), the consecutive wet days (CWD), the number of wet days index (R1MM), the number of wet days with moderate and heavy intensity precipitation (R10MM and R30MM, respectively), and annual and June to September daily mean precipitation (ANN and JJAS, respectively). The performance assessment and future change outlook were done for the CORDEX–Africa subdomains of north West Africa (WA-N), south West Africa (WA-S), and a combination of the two subdomains. While the performance of RCM runs was descriptorand subregionspecific, five model runs emerged as top performers in representing precipitation characteristics over both WA-N and WA-S. The five model runs are CCLM4 forced by ICHEC-EC-EARTH (r12i1p1), RCA4 forced by CCCma-CanESM2 (r1i1p1), RACMO22T forced by MOHC-HadGEM2-ES (r1i1p1), and the ensemble means of simulations made by CCLM4 and RACMO22T. All precipitation descriptors recorded a reduction under the two warming levels, except the SDII which recorded an increase. Unlike the WA-N that showed less frequency and more intense precipitation, the WA-S showed increased frequency and intensity. Given the potential impact that these projected changes may have on West Africa’s socioeconomic activities, adjustments in investment may be required to take advantage of (and enhance system resilience against damage that may result from) the potential changes in precipitation.

. Map of Coordinated Regional climate Downscaling Experiment (CORDEX) West Africa domains (marked WA-N and WA-S) where the analysis was done. Map adapted from https://bit.ly/2BR3qBJ and [25] under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).

RCM Data
Total daily near-surface precipitation data from 24 CORDEX-Africa RCM realisations (bias uncorrected and in native form) resulting from single or multiple ensemble runs (Table 1) were used in this study. The RCMs, described in detail by [26], were forced by CMIP5 GCMs, and have a spatial resolution of about 50 km at the equator [27]. Our study focused on realisations under RCP 8.5 given that it provides the largest set of ensemble members compared to other RCPs. Further, RCP 8.5 represents the most realistic warming scenario considering today's global greenhouse gas emission trajectory [28] and has been widely used for analysis of projections in sub-Saharan Africa [14,[29][30][31].

RCM Data
Total daily near-surface precipitation data from 24 CORDEX-Africa RCM realisations (bias un-corrected and in native form) resulting from single or multiple ensemble runs (Table 1) were used in this study. The RCMs, described in detail by [26], were forced by CMIP5 GCMs, and have a spatial resolution of about 50 km at the equator [27]. Our study focused on realisations under RCP 8.5 given that it provides the largest set of ensemble members compared to other RCPs. Further, RCP 8.5 represents the most realistic warming scenario considering today's global greenhouse gas emission trajectory [28] and has been widely used for analysis of projections in sub-Saharan Africa [14,[29][30][31]. Table 1. List of regional climate models (RCMs) used in our study (data accessed in March 2020 from http://bit.ly/2RoIist).

Observational Data
The daily Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), version 2.0 [32] for the period 1983-2005 was used as the reference observational precipitation data. The CHIRPS dataset ranges from 1981 to near present and integrates satellite imagery (0.05 • resolution) with in situ station data producing a gridded rainfall timeseries. To facilitate comparison, the Tropical Applications of Meteorology using Satellite and ground-based observations, version 3 (TAMSAT3) [33] was used as a second reference dataset for precipitation. TAMSAT3 is daily rainfall dataset based on high-resolution thermal-infrared observations and available from 1983 to the near-present. It is based on the disaggregation of 5-day total TAMSAT rainfall estimates to a daily time-step using daily cold cloud duration. In so doing, the dataset provides a temporally consistent historic and near-real time daily rainfall information for all of Africa. We also considered the daily version of the Global Precipitation Climatology Project dataset (GPCP-Daily). GPCP-Daily is based on the monthly version and produced by optimally merging precipitation estimates from infrared, microwave, and sounder data observed by the international constellation of precipitation-related satellites, and precipitation gauge analyses [34]. However, GPCP-Daily is only available from October 1996 to present, making it too short (9 years) for our baseline . CHIRPS and TAMSAT3 provide 25 and 23 years of analysis, respectively, and have been validated and used widely in studies on West African precipitation [33,35,36].
Computations were done on the native data grids [37], before bilinearly interpolating [38] the processed files to make it possible to compare observations and models. Additionally, computations were done for each year before averaging to get the climatology of descriptors for the periods of interest to minimize the risk posed by simple averaging [39,40].

Data Analysis
This study was done in two parts: (1) evaluating the performance of 24 CORDEX-Africa RCM realisations in reproducing West Africa's historical intra-seasonal precipitation characteristics and (2) estimating future precipitation change projections over the study domain using the top-performing model runs. The choice of top-performing model runs was informed by, in part, the existence of a large set of model runs produced by a growing number of models with differing characteristics and notable interdependence [41]. For instance, the skill of representing observed climate and trends differs from one model to another [42,43]. Using a large group of model simulations can be impracticable and pose a risk of bias given that the performance of models is strongly dependent on the scale and variable being investigated [44]. Additionally, several studies have shown that a large multi-model ensemble performs dismally, compared to an ensemble of few top-performing models in the group [41,[45][46][47].
Nonetheless, our choice of a few more-adequate models was not intended to discredit the performance of the other models. We appreciate the potential risk of biased or overconfident results especially if the number of models is small [48]. However, the risk is minimized when a robust method of model performance assessment is used. In this study, we used nine model runs from a pool of 24 simulations, to estimate future precipitation changes over the study domain.
Specifically, we evaluated the performance of RCM runs based on their skill in reproducing spatial and temporal patterns of West African precipitation (hereinafter 'descriptors'). The descriptors ( Table 2) represent three categories of daily precipitation indices defined by the Expert Team on Climate Change Detection and Indices (ETCCDI) [49]. The three categories are (i) the mean and strong intensity indices; (ii) occurrence indices; and (iii) duration indices, herein represented by (a) simple precipitation intensity index (SDII); (b) R1mm, R10mm, and R30mm; and (c) CWD, respectively. The indices have been widely used to study precipitation over West Africa [10,50] and beyond [47,51,52]. The other two descriptors are the annual (ANN) and seasonal (JJAS) daily mean precipitation.
Performance evaluation was done using simulations from five CORDEX-Africa RCMs forced by GCMs. The period 1983 to 2005 was chosen to correspond to data availability for reference data (CHIRPS and TAMSAT3) available from 1981 and 1983 to the near present, respectively. It also corresponds to historical simulations of the CORDEX-Africa RCMs (up to year 2005). The performance assessment was done in three tiers for a better representation of the study domain. Here, the 24 model runs were evaluated on their ability to represent precipitation characteristics over WA-S, WA-N, and a combination of WA-S and WA-N.
The performance of the RCM runs in reproducing the spatial patterns for all the precipitation descriptors, relative to observations, was done using spatial correlation coefficients (PCC). The temporal part was assessed using an inter-annual variability skill score (IVS) [45], defined as follows: where σ o and σ m represent the standard deviation for observations and models, respectively. An IVS value of 0 implies that σ o is equal to σ m , and the closer an IVS value is to 0, the better its skill in simulating the inter-annual variability [45].
Climate 2020, 8,143 6 of 20 The model performance assessment was done using detrended data to minimize the risk of potentially spurious variable relations [53]. To rank the performance of model runs based on the two criteria (temporal and spatial precipitation characteristics), a comprehensive rating index MR [46]; was used, as follows: where n is the number of descriptors and m the number of RCM runs. The model with an MR value closest to 1 gives the best skill of simulation, compared to other models under consideration [46].
computed for June-September of every year in the series mm/day Results from model ranking (MR) based on IVS and PCC, for the three sub-domains, were plotted in scatter diagrams to facilitate comparison. Here, two sets of MR values were calculated. The first set used PCC values where the model run with the highest correlation coefficient was ranked first, while the one with the lowest coefficient was ranked last. Given that the performance of model runs was ranked against seven descriptors, model run 'X' can be ranked number 1 for descriptor 1 and number 5 for descriptor 2, etc. Hence, a formula (Equation (2)) was used to rank the overall performance of individual model runs in all the seven categories (descriptors). The second set of MR values used model runs that were ranked based on their performance in reproducing year-to-year precipitation characteristics (IVS; Equation (1)). Here, the model run with the lowest IVS was ranked first while the model run with the highest IVS was ranked last. The MR equation was then used to give the overall ranking based on IVS. Lastly, MR values based on PCC were plotted against MR values based on IVS. The top-performing models (relatively) are found in the 1st quadrant (Q1), while the least performing models are found in the Q3 of the scatter plot. Other studies that have used this method include [46,47].
The top-performing model runs for each sub-domain were used to compute future precipitation changes for all the descriptors using projections under the RCP 8.5 scenario. A detailed analysis conducted by [54] identified (using a subset of GCM simulations driving CORDEX-Africa RCMs) years 2022 and 2037 as mid-years for 30-year windows, when GWL1.5 and GWL2.0, respectively, are likely to be first experienced. Their analysis was done with reference to the period 1861-1890 as the pre-industrial period. Our study adapted these mid-years and computed future changes in precipitation over West Africa by subtracting precipitation climatology values in the period 1976-2005 (hereinafter CTL) from the precipitation climatology values in periods 2008-2037, 2023-2052, and 2071-2100, corresponding to GWL1.5, GWL2.0, future (FUT), respectively, under RCP 8.5. The future changes were computed for all the seven descriptors.

RCM Performance Assessment
Scatter plots (Figures S1 to S6 in the Supplementary Materials) for model run rankings based on correlation coefficients (PCCs) and interannual variability scores (IVS) were made, to show the performance of the model runs. A summary of the top-performing model simulations in reproducing spatial and temporal precipitation characteristics over WA-N is shown in Table 3. Here, the top five simulations (agreed upon by both TAMSAT3 and CHIRPS) were CCLM4 forced by MOHC-HadGEM2-ES (r1i1p1), RCA4 forced by NCC-NorESM1-M (r1i1p1), IPSL-IPSL-CM5A-MR (r1i1p1), and ICHEC-EC-EARTH (r12i1p1), and an ensemble mean for model runs produced by RACMO22T.
Additionally, the performance of model simulations in reproducing spatial and temporal precipitation characteristics at both WA-S and WA-N was ranked (Table 5). Here, the model runs agreed upon by both reference data sets (TAMSAT3 and CHIRPS) were CCLM4 forced by ICHEC-EC-EARTH (r12i1p1), RCA4 forced by CCCma-CanESM2 (r1i1p1), RACMO22T forced by MOHC-HadGEM2-ES (r1i1p1), an ensemble mean for simulations done by CCLM4, and an ensemble mean for simulations done by RACMO22T.
With CHIRPS as the reference dataset, RCA4 driven by CCCma-CanESM2 (r1i1p1) and RACMO22T driven by MOHC-HadGEM2-ES (r1i1p1) were the only simulations featuring among the top five performers in reproducing spatial precipitation characteristics for WA-N, WA-S, and both WA-N and WA-S regions. RCA4 driven by CNRM-CERFACS-CNRM-CM5 (r1i1p1) was the only RCM simulation featuring among the top five performers in reproducing interannual precipitation characteristics for WA-N, WA-S, and both WA-N and WA-S regions. Considering the performance of simulations in reproducing both spatial and temporal precipitation characteristics (compared to CHIRPS), none of the model runs featured among the top five performers in all three categories (WA-N, WA-S, and both WA-N and WA-S).
In the case of TAMSAT3 as the reference dataset, no model simulation was featured among the top five performers in reproducing interannual precipitation characteristics in all three categories (WA-N, WA-S, and both WA-N and WA-S). CCLM4 and RACMO22T driven by MOHC-HadGEM2-ES (r1i1p1) and an ensemble mean of model runs driving the CCLM4 featured among the top five performers in reproducing the spatial precipitation characteristics in all three categories.
Overall, five model runs emerged as top performers in reproducing precipitation characteristics (both spatial and temporal) at both WA-N and WA-S, with reference to both TAMSAT3 and CHIRPS. The five model runs were CCLM4 forced by ICHEC-EC-EARTH (r12i1p1), RCA4 forced by CCCma-CanESM2 (r1i1p1), RACMO22T forced by MOHC-HadGEM2-ES (r1i1p1), and the ensemble means of simulations done by CCLM4 and RACMO22T. Assessing the performance of four CORDEX RCMs in reproducing climate characteristics over southern Ghana, Ashaley et al. [55] identified the RCA4 model forced by CCCma-CanESM2 as one of the top-performing models. Our results agree with this finding and add value by providing a list of additional top-performing model runs (and for a bigger domain) over West Africa. Table 4. As in Table 3 Table 3  Notably, the performance of model runs was location-and descriptor-specific. No single model run emerged best in reproducing all the seven descriptors-hence the need for a model ranking formula (Equation (2)). The relatively inconsistent performance of model simulations in reproducing the set precipitation descriptors from one region and criteria to another could be, partly, as a result of the predominantly convective nature of precipitation over the study domain [25,56]. The convective precipitation occurs in the mesoscale which begins at around 5 km-significantly lower than CORDEX-Africa RCMs' resolutions of about 50 km at the equator. Consequently, the relatively coarse resolutions of models make it difficult for them to adequately resolve the convective precipitation systems over the study domain [25,57]. Some of the mesoscale convective systems that affect West Africa include squall-lines [58], organized convective systems [59], and mesoscale convective systems [59,60].

Future Precipitation Changes over West Africa
Using a multi-model ensemble mean of the top-performing model runs determined in Section 3.1, projected changes in the spatial climatology of all seven precipitation descriptors for West Africa were plotted (Figures 2 and 3). During the control period (CTL; 1976-2005), daily mean annual precipitation ranging from about 2 mm at the north of the study domain to about 9 mm at the south was recorded. The JJAS season recorded precipitation ranging from about 4 mm to about 12 mm per day with the highest amount recorded over central (around 7.5 E, 9.1 N) and southern Nigeria (about 3-9.2 E, 5.8-6.2 N). Most of the study domain recorded a simple daily precipitation intensity index of between 8 and 12 mm. The highest intensity (exceeding 12 mm) was recorded over central and southern Nigeria.
Under the GWL1.5, no major significant changes in mean daily precipitation at the annual (ANN) and seasonal (JJAS) were recorded over the study domain. Where recorded, precipitation seemed to reduce by about 1 mm/day for both ANN and JJAS. A GWL2.0 was likely to introduce minor (about 1 mm/day) wet and dry biases at a few places over the study domain. An increase of between 1 and 2 mm/day was recorded for the simple daily precipitation index (SDII) over some parts north of the study domain. The SDII seemed to intensify under GWL2.0. Although the changes under GWL2.0 were more than changes under GWL1.5, most of the study domain recorded no significant changes for ANN, JJAS, and SDII. Under the GWL1.5, no major significant changes in mean daily precipitation at the annual (ANN) and seasonal (JJAS) were recorded over the study domain. Where recorded, precipitation seemed to reduce by about 1 mm/day for both ANN and JJAS. A GWL2.0 was likely to introduce minor (about 1 mm/day) wet and dry biases at a few places over the study domain. An increase of between 1 and 2 mm/day was recorded for the simple daily precipitation index (SDII) over some parts north of the study domain. The SDII seemed to intensify under GWL2.0. Although the changes under GWL2.0 were more than changes under GWL1.5, most of the study domain recorded no significant changes for ANN, JJAS, and SDII. precipitation intensity is likely to increase under a global warming scenario [54,61,62].
Looking at the descriptors that show the frequency of precipitation aspects (Figure 3), consecutive wet days (CWD) over the domain seemed to reduce with latitude with about 15 days/annum at the north and 75 days south of the study region. Minor significant changes (about -5 days) were recorded over a few areas under GWL1.5 and GWL2.0. By the end of the 21st century, most of the study domain recorded a reduction in CWD of between 5 and 15 days/annum. Areas in central Nigeria recorded the largest reduction (about 15 days/annum), while areas south of the study domain recorded no significant changes. Overall, the 1.5 °C and 2.0 °C global warming levels did not show significant effect on precipitation over many parts of the study domain. Specifically, all precipitation descriptors recorded a reduction under the two warming scenarios except SDII which recorded an increase. The northern part of the study domain recorded most changes compared to the south. Considering the future Projections (under RCP 8.5) showed a reduction in daily mean precipitation of up to 2 mm in several areas to the north and east of the study domain for ANN and JJAS, by the end of the 21st century. However, an increase of up to 2 mm in mean daily precipitation for ANN and JJAS was recorded over areas south-west of the study domain. Most of the study domain recorded an increase in SDII by the end of the 21st century. The increase was more pronounced to the north of the domain with changes exceeding 3 mm/day recorded in some areas. No significant changes were recorded over most of Nigeria and southern Togo and Benin (south-east of the study domain). A few areas south-west of Nigeria recorded a decrease in SDII of up to 2 mm.
Our findings on ANN, JJAS, and SDII are consistent with similar works done in the study domain. Evaluating the impacts of GWL1.5 and GWL2.0 over Senegal, Mbaye et al. [61] established a possible marginal decrease in annual and July-September mean precipitation. Their study also pointed to the possibility of a slight increase in heavy precipitation events over Senegal under the two global warming levels for the RCP 8.5 scenario. Looking at Africa in general, Nikulin et al. [54] reported a tendency towards wetter conditions in annual average precipitation over central and eastern Sahel, but with a low model agreement. There seems to be a consensus (consistent with our findings) that the daily precipitation intensity is likely to increase under a global warming scenario [54,61,62].
Looking at the descriptors that show the frequency of precipitation aspects (Figure 3), consecutive wet days (CWD) over the domain seemed to reduce with latitude with about 15 days/annum at the north and 75 days south of the study region. Minor significant changes (about -5 days) were recorded over a few areas under GWL1.5 and GWL2.0. By the end of the 21st century, most of the study domain recorded a reduction in CWD of between 5 and 15 days/annum. Areas in central Nigeria recorded the largest reduction (about 15 days/annum), while areas south of the study domain recorded no significant changes.
An average of about 50 rainy days (R1MM) were recorded in the north of the study domain and at least 150 days/annum in the south. Most of the domain recorded between 10 and 30 days of wet days with moderate intensity precipitation (R10MM). For wet days with heavy intensity precipitation (R30MM), most of the domain recorded less than 15 days/annum. No major significant changes were recorded over the domain for R1MM, R10MM, and R30MM, save for a slight reduction of up to 10 days/annum, especially for R1MM over some areas north of the domain. Over the future (FUT) period, the greatest change (exceeding 15 days/annum in most of the domain) was expected for R1MM. R10MM and R30MM recorded a slight reduction of about 5 days/annum in many parts of the domain except the south-west where an increase of about 5 days/annum was recorded.
Overall, the 1.5 • C and 2.0 • C global warming levels did not show significant effect on precipitation over many parts of the study domain. Specifically, all precipitation descriptors recorded a reduction under the two warming scenarios except SDII which recorded an increase. The northern part of the study domain recorded most changes compared to the south. Considering the future period (FUT, 2071-2100), SDII and R1MM recorded the greatest change, with many areas exceeding 2 mm/day and -15 days/annum, respectively. A few areas to the south-west of the domain recorded more episodes of wet bias than the rest of the study domain.
To enhance the understanding of future precipitation changes over the north and south subdomains of our study area, boxplots were done using a set of nine top-performing model simulations determined in Section 3.1. The boxplots showed the spread of the nine CORDEX RCM runs in reproducing the seven precipitation descriptors over the two subdomains as detailed in the following subsections.

Future Precipitation Changes over the WA-N Sub-Region
First, precipitation changes over the northern part of the study domain (WA-N) were assessed ( Figure 4). Here, the spread of model runs was greatest (relatively) for R30MM, implying a lower model agreement compared to other precipitation descriptors. Compared to the control period (CTL), no noticeable changes were recorded for R1MM and R30MM under GWL1.5. A slight decrease (increase) was recorded for Annual, JJAS, and CWD (SDII and R10MM). period (FUT, 2071-2100), SDII and R1MM recorded the greatest change, with many areas exceeding 2 mm/day and -15 days/annum, respectively. A few areas to the south-west of the domain recorded more episodes of wet bias than the rest of the study domain.
To enhance the understanding of future precipitation changes over the north and south subdomains of our study area, boxplots were done using a set of nine top-performing model simulations determined in Section 3.1. The boxplots showed the spread of the nine CORDEX RCM runs in reproducing the seven precipitation descriptors over the two subdomains as detailed in the following subsections.

Future Precipitation Changes over the WA-N Sub-Region
First, precipitation changes over the northern part of the study domain (WA-N) were assessed ( Figure 4). Here, the spread of model runs was greatest (relatively) for R30MM, implying a lower model agreement compared to other precipitation descriptors. Compared to the control period (CTL), no noticeable changes were recorded for R1MM and R30MM under GWL1.5. A slight decrease (increase) was recorded for Annual, JJAS, and CWD (SDII and R10MM). Under GWL2.0, a further reduction (increase) in annual, JJAS, and CWD (R30MM, R10MM, and SDII) was recorded. No conspicuous changes were recorded for R1MM. Under the future scenario (FUT), discernible negative (positive) changes were recorded for annual, JJAS, R30MM, and CWD (SDII). No major changes were recorded for R1MM and R10MM. As in the spatial plots, the boxplots showed a reducing trend in annual mean daily precipitation (ANN) and CWD, but an increasing trend in SDII. This implies a potential for a decreasing frequency and increasing intensity of Under GWL2.0, a further reduction (increase) in annual, JJAS, and CWD (R30MM, R10MM, and SDII) was recorded. No conspicuous changes were recorded for R1MM. Under the future scenario (FUT), discernible negative (positive) changes were recorded for annual, JJAS, R30MM, and CWD (SDII). No major changes were recorded for R1MM and R10MM. As in the spatial plots, the boxplots showed a reducing trend in annual mean daily precipitation (ANN) and CWD, but an increasing trend in SDII. This implies a potential for a decreasing frequency and increasing intensity of precipitation over the WA-N region.

Future Precipitation Changes over the WA-S Sub-Region
Secondly, an analysis of future change projections over the southern part of the study domain (WA-S; Figure 5) was done. Unlike the case of WA-N, where a decreasing trend was observed, WA-S showed an increasing trend for ANN and JJAS from the GWL1.5 to FUT scenarios. SDII showed an increasing trend while R30MM, which showed no remarkable change over WA-N, showed an increase in the three scenarios relative to the CTL.
Except for ANN, JJAS, and SDII, the difference between the impact of GWL1.5 and GWL2.0 on precipitation over WA-S was minimal. Changes during the FUT period were more pronounced for all the descriptors compared to GWL1.5 and GWL2.0. All but two descriptors (CWD and R1MM) showed a positive change during the FUT period than the CTL period. Overall, unlike the WA-N, where less frequent but more intense precipitation was recorded, the WA-S showed potential for more frequent and intense precipitation under the three scenarios.  Except for ANN, JJAS, and SDII, the difference between the impact of GWL1.5 and GWL2.0 on precipitation over WA-S was minimal. Changes during the FUT period were more pronounced for all the descriptors compared to GWL1.5 and GWL2.0. All but two descriptors (CWD and R1MM) showed a positive change during the FUT period than the CTL period. Overall, unlike the WA-N, where less frequent but more intense precipitation was recorded, the WA-S showed potential for more frequent and intense precipitation under the three scenarios.

Discussion
Our findings agree with similar studies that have shown the possibility of enhanced precipitation over West Africa in the future [63] and add value by (i) delineating the expected future

Discussion
Our findings agree with similar studies that have shown the possibility of enhanced precipitation over West Africa in the future [63] and add value by (i) delineating the expected future changes in two sub-domains of West Africa; and (ii) assessing more precipitation characteristics for a better representation of West African precipitation. For instance, we show that the number of wet days with heavy intensity precipitation (R30MM) is projected to increase (remain unchanged) over WA-S (WA-N) under GWL1.5, GWL2.0, and FUT. Our results also show a reduction in consecutive wet days for both WA-N and WA-S with the reduction over WA-S being more pronounced. These results qualify earlier studies that alluded to the possibility of extreme rainfall events [64], increased precipitation intensity [65,66], and increased dry spells [67] over the study domain.
Over 75% of West Africa's population lives in areas that have been affected at least once in every two years by flood, dust/sandstorm, and drought in the last three decades alone [68]. This comes with great impacts on livelihoods and requires holistic planning and strategies that incorporate different scenarios and early warning development mechanisms. Yet, many countries in West Africa lack adequate meteorological and hydrological capacity to collect, process, and disseminate climate information for early warning to vulnerable communities [68]. Therefore, the results of this study provide useful information to support decision making and investments in water resources management, agriculture, and energy-the most affected sectors of the West African economy.
Our study corroborates previous works that inferred a possibility of more floods and crop water stress over the study domain [69]. Additionally, we show that the southern part of West Africa (WA-S) will potentially experience more intense and frequent precipitation compared to the northern part (WA-N), which records potential for less frequent and more intense precipitation. As West Africa experiences rapid population increase, rapid demographic and socioeconomic transformations, urbanization, and rising incomes, greater stress is placed on food systems. Hence, tough decisions will have to be made to minimize greenhouse gas emissions and adapt agricultural systems in preparation for possible future precipitation changes. The possibility of a drying climate calls for common but differentiated climate change adaptation strategies to address current and future national, regional, and transboundary challenges in the study domain.
As a caveat, our study relied on a study done by [54] to determine periods corresponding to the 1.5 • C and 2.0 • C global warming levels. Further, our study focused on a subset of GCMs participating in CORDEX. Various sub-sets of GCMs and the choice of the reference pre-industrial period may give varying estimates of the future period when the 1.5 • C and 2.0 • C global warming levels are likely to be reached. Hence, we recommend further research to explore various scenarios and enhance the understanding of the impacts of global warming levels on climate systems. Additionally, a large ensemble of climate projections can introduce bias especially when ensembles consist of models with similar components or from the same family [42]. A common issue with CORDEX simulations is that not all RCMs downscale all GCMs. For instance, the RCA4 model contributed 14 out of the 24 simulations considered in our study. The sparse RCM-GCM matrix can potentially introduce a bias towards specific RCMs and hence, potentially lead to overconfidence on future climate change scenarios. Our study does not address these limitations and only provides insights aimed at spurring further research and discussions on future climate changes over the study domain.

Conclusions
Through a relatively low model agreement, this study confirmed the difficulty of climate models (at the current spatial resolutions) to represent West African precipitation adequately due to, in part, mesoscale climate systems [25,57,58]. The performance of the 24 model simulations was descriptor and subdomain specific with no single model run featuring as a top performer in representing the spatial and temporal characteristics of all the seven precipitation descriptors, across all the three sub-domain categories. Nonetheless, five model runs emerged as top performers namely CCLM4 forced by ICHEC-EC-EARTH (r12i1p1), RCA4 forced by CCCma-CanESM2 (r1i1p1), RACMO22T forced by MOHC-HadGEM2-ES (r1i1p1), and the ensemble means of simulations done by CCLM4 and RACMO22T. Assessment of potential future precipitation changes under 1.5 • C and 2.0 • C global warming levels under the RCP 8.5 scenario showed a minimal impact on West African precipitation, relative to the 1976-2005 control period (CTL). All precipitation descriptors recorded a reduction under the two warming levels, except SDII, which recorded an increase. Unlike the WA-N, which showed less frequent and more intense precipitation, the WA-S showed an increase in frequency and intensity. These changes seemed to be more robust in the future period (2071-2100) under the RCP 8.5 scenario.
Given the potentially significant impact these projected changes may have on West Africa's socioeconomic activities, adjustments in investment may be required to take advantage of the heavier and more intense precipitation for enhanced productivity and water availability. The potential decrease in the number of rainy days and the consecutive wet days call for investment in water harvesting initiatives to ensure water availability for domestic and agricultural use during the dry periods. The potential for more heavy and intense precipitation events poses a risk to infrastructure over the study domain, which is vulnerable to floods [70,71].
Supplementary Materials: The following are available online at http://www.mdpi.com/2225-1154/8/12/143/s1, Figure S1: Model run ranking (relative to CHIRPS) based on correlation coefficients (PCCs) and interannual variability scores (IVS), for the WAN region, Figure S2: As in Figure S1 but for the WAS region, Figure S3: Model run ranking (relative to TAMSAT3) based on correlation coefficients (PCCs) and interannual variability scores (IVS), for the WAN region, Figure S4: As in Figure S3 but for WAS region, Figure S5: Model run ranking (relative to TAMSAT3) based on correlation coefficients (PCCs) and interannual variability scores (IVS), for both WAN and WAS region, Figure S6: As in Figure S5