Assessment of Climate Change Impact on Water Resources in the Upper Senegal Basin (West Africa) ()
1. Introduction
Long-term observational records and climate projections provide abundant evidence that freshwater resources are vulnerable and have the potential to be strongly impacted by climate change, with wide-ranging cones- quences for human societies and ecosystems [1] . These consequences are more severe in regions dominated by arid and semiarid climate such as West Africa (WA). Then, Africa as a whole is one of the most vulnerable con- tinents due to its high exposure and low adaptive capacity [2] . Several impact studies over WA have shown that water resources are significantly impacted by climate change [3] -[9] . For much of Africa, knowledge about recent climate change is limited, due to weak climate monitoring, and gaps in coverage that continue to exist [10] . The well known droughts in the 1970s have led to a decline of water flows in many African river basins. Also, the Senegal River Basin, which is situated in WA, has faced those droughts. Its annual average flow at Bakel (reference station) fell from 1374 m3/sec over the period 1903-1950 to 840 m3/sec in the period 1950-1972 and further decreased to only 419 m3/sec in the period 1973-2002 [11] . The drought has reduced rain-fed agriculture, decreased the seasonal flooding of wetlands, limited economic development, and in the overall, enhanced poverty [12] . Such water shortage has obviously affected the main activities in the basin (agriculture, fishery, hydropower generation, etc.). As water is the agent that delivers many of the impacts of climate change to society, for example to the energy, agriculture, and transport sectors [10] , it is fundamental to investigate how climate change will affect this valuable resource in the future. A series of studies has already been carried out over the Senegal River basin in order to analyze the ability of models to simulate the basin’s climate, but also to quantify projected climate changes and to identify potential climate change impacts. [3] evaluated the ability of four General Circulation Models (GCMs) to reproduce the basin’s climate and they found that all four models failed to reproduce the rainfall volumes in the Sahelian zone, and the seasonal dynamics of rainfall in the Guinean zones. The biases were mainly introduced by the limited quality of the GCM input data, which can also be related to the coarse resolution of the GCMs, which does not allow representing the important small scale processes. Analyzing projected changes, an increase in dry spell length by 2050 was found by [4] , and [13] by analyzing regional climate models (RCMs). The higher resolution of RCMs compared to GCMs makes it possible to simulate regional climate features more realistically such as orographic precipitation, extreme events, and regional scale climate anomalies, or non-linear effects [14] . But the climate models (e.g. RCMs) are subject to systematic biases when comparing simulated meteorological variables for the current climate to observations and these biases can affect hydrological simulations considerably [15] . In addition, [16] stated that models will not be able to simulate accurate predictions if the areal precipitation is not adequately represented. To date, there is a lack of impact studies on the seasonal changes of hydrological variables (such as soil moisture, evapotranspiration, and runoff), extreme streamflow changes, and their variability using high resolution climate data as input for a hydrological model over Upper Senegal Basin (USB). However, the better understanding of potential future changes in the spatial and temporal variability of the hydrological cycle is fundamental to inform local societies and water resources managers in order to increase awareness and to support the development of adaptation strategies, and also agricultural productivity which accounts for more than 80% of the water withdrawals in the basin [12] .
In the present study, bias corrected data are used as input for a hydrological model to overcome the well known biases in GCM/RCM output. Taking into account the concern of the potential alteration of the climate change signal due to bias correction [17] , both uncorrected and bias corrected data from the regional climate model REMO are used to drive the Max Planck Institute for Meteorology-Hydrology model (MPI-HM). The main purposes of this study are to investigate if there is an added value of applying bias correction, to assess the potential changes of hydrological variables and water resources availability, and to explore whether bias correction may affect the climate change signal of projected hydrological variables over the USB.
The study is structured as following: data and methods are described in Section 2; the results and discussion are presented in Section 3; and the conclusion is given in Section 4.
2. Data and Methods
2.1. Upper Basin Senegal
The current study focuses on the USB in WA (Figure 1). From its source in Guinea the Senegal River flows through the western Sahel region in Mali, Mauritania and Senegal and has a catchment size of about 218,000 km2. Its three main tributaries (Bakoye, Bafing and Faleme) provide over 80% of its flows and are within the upper basin. The basin is subject to a large north-south precipitation gradient ranging from 200 mm/year in the north to more than 1800 mm/year in the south. The predominantly natural vegetation of the region follows this rainfall gradient, ranging from semi-arid savannah in the North to sub-humid forest in the South [18] .
Bakel is considered to be the reference station of the Senegal River due to its location below the confluence with the last major tributary (Faleme). At this station, the average annual discharge is about 690 m3/sec, which corresponds to an annual input of around 22 × 109 m3. The annual discharge ranges between a minimum of 6.9 bm3 and a maximum of 41.5 bm3 [19] .
2.2. Climate Data
Daily data are taken from a set of simulations (historical and scenario) conducted with the regional climate model REMO [20] . At the model boundaries, the REMO simulations are forced with data from the global climate model MPI-ESM-LR [21] following the RCP4.5 and RCP8.5 emission scenario pathways. REMO data is available in the context of the Coordinated Regional Climate Downscaling Experiment [22] over Africa at 0.44 resolution for the period 1950 to 2100 and it has already been used over Africa [23] -[26] . The data analyses focus on the periods (1971-2000) as reference period and (2071-2100) as scenario period. As precipitation and temperature are the key drivers for the hydrological regime of the river and climate change has its main impact through changes in these two variables [9] both were bias corrected following the method of [27] . Together with radiation, wind and air moisture data, the bias corrected and uncorrected temperatures are then used in the Penman-Monteith equation to estimate the potential evapotranspiration. Evapotranspiration is generally the most substantial component of the water budget, being the pathway for local rainfall recycling, in particular, over the Sahel region [28] . These three variables (precipitation, temperature, and potential evapotranspiration) were used as climate input of the Max Planck Institute for Meteorology-Hydrological Model.
2.3. The Max Planck Institute for Meteorology-Hydrological Model
The MPI-HM [29] is a global hydrological model which consists of the Simplified Land Surface (SL) scheme [30] and the Hydrological Discharge (HD) model [31] -[33] . This model is chosen because it has been widely applied in impact studies over many different rivers around the world, including West African rivers [34] [17] -[35] . It uses a half degree resolution (0.5˚) and daily time steps. The SL scheme is designed to compute the vertical water fluxes while the HD model computes the lateral water fluxes and river routing on the land surface (Figure 2). The SL scheme uses a soil bucket scheme for the calculation of the vertical water balance. Its output comprises surface runoff and drainage that are used directly by the HD model which is a state of the art river routing model. In the HD model, the lateral water flow is separated into the three flow process of overland flow, base flow and river flow. Overland flow uses surface runoff as input, base flow is fed by drainage from the soil and the inflow from upstream grid boxes contributes to river flow. The sum of the three flow processes equals to the total outflow from a grid box [30] . Initially, the two sub-models were being run separately, this coupled ver- sion (MPI-HM) has been constructed by [29] who also implemented a dynamical wetland extent scheme. In this study, neither wetland dynamics are taken into account nor human activities (dams, water withdraw, etc.) nor
Figure 1. Location of upper senegal basin (green area).
land use changes. The MPI-HM model requires as climate input temperature, precipitation, and potential evapotranspiration fields; it uses land sea mask, glacier mask, total field capacity and plant-available soil water capacity from LSP2 dataset [36] that has recently been updated [37] . The other boundary data are sub-grid slope, elevation, river routing direction, and storage retention times. In this study, the global standard version of the model has been modified to meet the local conditions in USB. The overland flow retention constant, the river flow retention constant, minimum soil moisture threshold for drainage and critical soil moisture for drainage, and the evaporation factor coefficients have been optimized by manual adjustments.
The simulated variables of interest are river discharge, surface runoff, drainage, evapotranspiration, and soil moisture. All simulations were carried out from 1965 to 2100; the first six years of simulations were considered as spin-up and removed from the analysis time series in order to avoid the influence of initial conditions in the water storages. All fluxes are computed at every land surface grid cell for every time step. More details about the functioning of the model can be found in [29] [30] -[33] .
Model Performance
In order to evaluate the performance of the optimized MPI-HM model over the USB, the following statistics are used:
Nash Sutcliffe Efficiency [38] which is commonly used to assess hydrological models goodness fit:
Figure 2. Structure of MPI-HM. Orange boxes indicate input variables, green boxes are water storages, and black arrows indicate water fluxes.
(1)
Percent Bias (PBIAS): it is used to estimate the model error
(2)
denotes the observed discharge at the month j, is the simulated one, and is the averaged observed streamflow. NSE ranges from to 1.0, with high values indicating better agreement. NSE = 1 indicates a perfect match of simulated and observed data; NSE = 0 corresponds to the model predictions matching the mean of the measured data, and NSE < 0 shows that the measured mean is a better predictor than the model [39] .
The PBIAS measures the tendency of the simulated variable to be larger or smaller than the corresponding observed value, its optimal value is 0.0, and positive values indicate a tendency for over-estimation whereas negative values indicate a tendency for under-estimation [40] .
2.4. Simulation Analyses
The analyses have been mainly focused on the seasonal cycles of river discharge, runoff, actual and potential evapotranspiration, soil moisture and precipitation. In addition, the seasonal spatial variability was investigated from July to October (JASO) which represents the high flow period. The changes represent the differences between the scenario period (2071-2100) and the reference period (1971-2000). Then, the change of available water resources which depends mainly on mean annual runoff is estimated as following:
(3)
is the annual mean runoff in the scenario period denotes the annual mean runoff in the reference period, and EWR is the Environmental Water Requirements that represents 30% of the mean annual runoff in the basin [41] . Thus, Equation (3) can be simplified to:
(4)
3. Results and Discussion
3.1. Precipitation Biases and Changes
This section assesses mainly the biases and changes of precipitation that is one the most important climate forcing data of the MPI-HM model. Then, bias corrected and uncorrected simulations are compared in order to see how their different variations may affect the hydrological simulations. Figure 3 shows the precipitation biases for uncorrected and bias corrected data during the reference period.
The wet biases during the main rain season (Figure 3(b)) and dry biases during the post and pre-monsoon in the uncorrected data have been removed in the corrected data (Figure 3(a)). Thus, the bias correction has well improved the variations of precipitation. This result is similar to that found also in the other forcing data, i.e. temperature and potential evapotranspiration (not shown). As for the projected precipitation changes (Figure 4), no substantial change in rainfall occurs from November to March for all datasets in all scenarios.
However, a general decrease of precipitation is found in the summer period (particularly in RCP8.5). This de- crease is higher in the uncorrected data than in the bias corrected data.
The climate change signal exhibits the same tendency of decrease and increase of precipitation, except in September where a noticeable increase is found in the uncorrected data with RCP4.5. The bias correction tends to reduce the magnitude of the climate change signal of precipitation over the USB. These differences in magnitude may also affect the hydrological simulations by using these data as climate forcing. In addition, the PET is projected to increase in all months for all data (not shown) by the end of 21st century. Then, the decrease of precipitation combined with the increase of PET is likely to considerably affect the water resources in the USB.
3.2. Evaluation of River Discharge Simulations during the Present Day Climate
The MPI-HM river discharge simulations are evaluated by comparing simulated and observed hydrographs. The
(a) (b)
Figure 3. Precipitation biases of bias corrected (a) and uncorrected (b) precipitation.
Figure 4. Precipitation changes of bias corrected and uncorrected simulations.
main goal within this section is to assess whether the bias corrected input enables a better reproduction of observed river regimes, and how the biases of climate forcing data affect simulated river discharge. In Figure 5(a), the simulated mean river discharge is compared to observed mean discharge during the reference period (1971- 2000) at the outlet of the USB (Bakel). This figure demonstrates a clear improvement of simulated discharge when using bias corrected input. The uncorrected data highly overestimate the monthly discharge. During the dry season (January to June), the corrected data tend to overestimate the low flow and then slightly underestimate the beginning of high flow (July-August). The overestimation during the dry season may be due to the limitation of the hydrology model in properly handling extreme dry climate. As for the low flows (10th percentile) and high flows (90th percentile) presented in Figure 5(b) and Figure 5(c), respectively, the model largely overestimates the river flows by uncorrected data forcing as well. Using bias corrected input, the peak flow in September is well captured although the model overestimates the low flows in the dry season (Figure 5(b)). The huge biases from the uncorrected data input are related to the wet biases of uncorrected precipitation (as seen in Section 3.1). In addition, an underestimation of ET could also contribute to the overestimation of river discharge during the main rainy season.
(a) (b)(c)
Figure 5. Simulated and observed mean (a) 10th percentile (b) and 90th (c) percentile of river discharge at Bakel.
Moreover, the observations are affected by flow regulation at the Manantali dam (upstream), which is not considered in the MPI-HM and might explain a part of the mismatch. The results show that the major differences of simulations occur during the high flow period with very huge values in the uncorrected data.
3.3. Model Performance
The model performance is evaluated by using the simulated mean monthly discharge. A Nash-Sutcliffe Efficiency coefficient of 0.92 is found between the simulated and the observed streamflow with bias corrected input. This high efficiency shows a good agreement between both datasets, indicating the chosen model combination (REMO + MPI-HM) is able to reproduce well the timing and the magnitude of observed streamflows. However, with the uncorrected input the efficiency is −21.98 that shows how far the simulation deviates to the observation. Furthermore, the PBIAS of bias corrected input is 3.88%, which shows that the simulated discharge slightly deviate from those observed; by contrast to the uncorrected input where the PBIAS is up to 300% which depicts the huge overestimation of the model.
Hence, in this evaluation part, the bias corrected input provides realistic river discharge simulations. Therefore, bias correction is useful for a better representation of present day river discharge, so that it can help to reconstruct historical mean streamflows which are fundamental bases for future adaptations strategies planning.
3.4. Potential Changes of Hydrological Variables with and without Bias Correction
This section focuses mainly on the potential changes of river discharge, runoff, evapotranspiration, soil moisture by the end of 21st century (2071-2100) under RCPs 4.5 and 8.5. In addition, the change of the projected signal from corrected and uncorrected input is considered in order to see if the bias correction alters the climate change signal.
3.4.1. Monthly River Discharge Changes
The changes in river discharge at the upper basin outlet (Bakel) and its 10th and 90th percentiles changes are presented in Figure 6. For all simulations, a general decrease is found, particularly in RCP8.5, except for August and September where an increase is found in RCP4.5. The decrease is less pronounced for the low flows than the high flows, and the RCP4.5 low flows show even a slightly increase during the peak flows period. This increase can be explained by lower evapotranspiration and small increase of precipitation. The bias corrected input neither substantially changes the signal in the mean discharge nor do the signals for low and high flows, even though slight differences exist in term of magnitude.
3.4.2. Runoff Seasonal Change
Figure 7 displays the seasonal change of runoff from July to October. Considerable decrease of runoff is found in the majority of the basin where the decrease reaches −80% in the eastern basin with RCP8.5. With this later scenario, the runoff decline is higher in the uncorrected data than in the bias corrected simulation. As with the moderate scenario (RCP4.5), it is noticeable also that runoff is likely to decrease in the USB, except localized parts such as the mountainous areas where slight increase is depicted. The runoff decrease found over the basin in all scenarios and data is more related to the decline of precipitation and higher evaporation. Additionally, the bias corrected and the uncorrected simulations present similar spatial pattern of the climate change signal in both scenarios. However, small differences exist in the magnitude of the projected signals of both simulations.
3.4.3. Evapotranspiration Seasonal Change
The seasonal change of actual evapotranspiration presented in Figure 8 shows a decrease of ET from 12N to 17.5N for uncorrected and bias corrected data with RCP8.5 where the peak decrease reaches −40% in the northern basin in both data. As with RCP4.5, a small decrease (−10% to 0%) is found south of latitude 13N for both datasets where the northern basin may face a decrease of −10%. The decrease of ET can be explained by the decline of soil moisture content. Furthermore, no change is found in the southern basin for all data and scenarios. From these results, it is noticeable that the bias corrected data have the same spatial pattern of the projected climate change signal as the uncorrected data. However, slight differences occur in the magnitude of the signal that seems to be reduced by the bias correction.
Figure 6. Mean monthly discharge change and percentiles (10th and 90th) river discharge changes.
Figure 7. Seasonal changes of total runoff [UC: uncorrected, BC: bias corrected].
Figure 8. Seasonal changes of evapotranspiration [UC: uncorrected, BC: bias corrected].
3.4.4. Soil Moisture Seasonal Change
Figure 9 displays the seasonal change of soil moisture. By the end of 21st century, the soil water content is projected to drastically decrease over the majority of the basin from the central to the northern part (−40% to −10%) in RCP8.5 for bias corrected and uncorrected data. The soil moisture decrease ranges from −20% to 10% in the moderate scenario (RCP4.5) for both simulations. This decrease is due to the decrease of precipitation and the increase of temperature. As seen previously with ET changes, similar spatial pattern is found in the climate change signal of soil moisture between the raw and corrected simulations. This shows that the projected signal is not altered by the bias correction.
3.4.5. Changes in the Seasonal Cycle of Evapotranspiration, Soil Moisture, and Runoff
The seasonal cycle changes of soil moisture, evapotranspiration and runoff are shown in Figure 10. Soil moisture and evapotranspiration are generally decreasing in all data and scenarios. The decrease is more pronounced with RCP8.5 than with RCP4.5. The maximum ET changes are reached in June and November (−0.8 mm/day) for all data. This is mainly related to the increasing temperatures whose changes also peak during these months. The rising temperatures increase the atmospheric water demand and, thus, enhance evaporation. In case of soil moisture decrease, the peaks are reached in September with RCP4.5 (−50 mm) and RCP8.5 (−135 mm) for uncorrected and bias corrected simulations. Hence, the climate change signal within the same scenario has a similar tendency in both data. As for runoff, the main changes (decreases) occur during the period of high flows (JASO) for all simulations and scenarios. In the uncorrected data, the changes are higher in RCP8.5 (−0.8 mm/day) than in the corrected data (−0.2 mm/day). With RCP4.5 scenario, the potential change of runoff signal with bias corrected input is very small by comparison to the runoff signal with uncorrected climate input. It is
Figure 9. Seasonal changes of soil moisture [UC: uncorrected, BC: bias corrected].
Figure 10. Climatological changes of evapotranspiration, runoff, and soil moisture [UC: uncorrected, BC: bias corrected].
noticeable with this later scenario that the uncorrected data depict two noticeable runoff peaks in August (increase) and September (decrease); these peaks are not found with corrected input. This later finding shows that the bias correction may change the projected runoff signal in some months.
3.4.6. Annual Changes of Available Water Resources
Figure 11 clearly shows that by the end of 21st century, available water resources are likely to decrease drastically in the majority of the basin, particularly under RCP8.5 (where the peak reaches −100%). This decrease is enhanced northward and becomes more acute in the most northeast part. Consistent with results mentioned above, no change of available water resources is found around the Guinean highlands, except for the uncorrected input under RCP4.5. In some localized parts of the basin, slight increases are found particularly in the uncorrected data (RCP4.5). All datasets and both scenarios exhibit similar spatial pattern. The decrease of available water follows the change pattern of runoff which is affected by the decrease of precipitation and the increase of potential evapotranspiration. Thus, the bias correction can alter the climate change signal of available water resources in some specific localities, but the spatial pattern is generally preserved.
4. Conclusions
An assessment of climate change impacts on water resources in the USB has been carried out. Despite considerable progress in recent years, output of both global and regional circulation models is still afflicted with bi- ases to a degree that precludes its direct use, especially in climate change impact studies [42] . Then, to cope with this issue, bias correction is often used in hydrological impact studies although for specific locations and
Figure 11. Available water resources changes [UC: uncorrected, BC: bias corrected].
seasons it may affect the climate change signal [17] . For these reasons, bias corrected and uncorrected climate data from the regional climate model REMO were used as input for the MPI-HM model to simulate river discharge, runoff, soil moisture and evapotranspiration. This is in agreement with the work of [43] who recommended presenting both results with and without bias correction in situations where only a few climate simulations are used in case that the climate change signals of outliers can be modified by bias correction. The global MPI-HM model parameters were optimized for local conditions of the USB in order to better capture the river flow variations. The results during the historical period (1971-2000) show that the bias corrected input give better representation of the river flow regimes at the outlet of the upper basin (Bakel) in the seasonal cycle. The same applies to the simulation of low flows (10th percentile) and high flows (90th percentile), thereby showing a slight overestimation of stream flow during the dry season. This can be related to the influence of other global parameters for which a tuning to local conditions might be necessary, or by missing processes such as the representation of soil crusts or human impacts. The results from uncorrected climate input are largely overestimated, which indicates the necessity of the bias correction. In this primary part, we conclude that the bias correction of precipitation and temperature has a great added value in reproducing present day river flows.
Then, by the end of 21st century (2071-2100) under the RCPs 4.5 and 8.5, a general decrease of river discharge, runoff, actual evapotranspiration, soil moisture is projected even though there are some localized increases in some parts of the basin (particularly in Guinean highlands) with the uncorrected simulations. This decrease is mainly related to the decline of precipitation. The most extremes changes of soil moisture, ET, and runoff are likely to occur in the northern basin which is the driest and hottest part of the USB. Additionally, the available water resources exhibit substantial decrease (from −100% to −25%) in the majority of the basin for all data, except the Guinean highlands where an increase (50%) is found under RCP4.5 in the uncorrected data. The comparison of results from uncorrected and bias corrected input (in all variables) demonstrates that the bias correction does not substantially change the signal of future changes of hydrological variables for both scenarios over the Upper Senegal Basin even though there are slight differences in term of magnitude of the projected signal in some part of the basin.
The projected changes over the USB are associated with various sources of uncertainties, especially those arising from the choice of the various models used in the climate model―hydrology model modeling chain. These comprise the driving GCM, MPI-ESM-LR, the RCM REMO used for downscaling and the hydrological model MPI-HM. In addition, the effect of the bias correction technique is another source of uncertainty. These uncertainties need to be taken into account when interpreting the results. To reduce these uncertainties, it would be desirable to conduct analogous hydrological simulations with different components of the modeling chain, i.e. using different RCM outputs, different bias correction techniques, and different hydrological models, which would enable uncertainty analyses in future investigations.
Furthermore, general improvements will be needed such as implementation of soil crusting and a deeper soil water layer into the MPI-HM which might result in better representation of the hydrological regimes. The implementation of a deeper layer will allow more water to infiltrate into the soil, and increase the bucket from which the more slow process of drainage occurs. This will reduce the surface runoff that dominates the flows during the high flow period where the overestimation mainly occurs. Implementing the effect of soil crusting can take into account the quick runoff commonly caused by rainstorms over crusted areas in very dry conditions. We recommend also including the impacts of human activities and land uses changes on river flow for further investigations. The implementation of human activities such as the Manantali dam where a lot of water evaporates from its reservoir, will take into account the water losses from this dam reservoir, the denaturalized flow and reduce also the overestimation of river flow because this evaporation is currently not included in MPI-HM. In this semi-arid region, the degradation of the natural vegetation for farming, agricultural, and urban purposes, reduce infiltration and soil water holding capacity. These situations are likely to be exacerbated in the future and can enhance runoff, so that taking those processes into account will allow for a better representation of current and potential future flow regimes.
Acknowledgements
The authors thank the German Federal Ministry of Education and Research (BMBF) for funding the West African Science Service Center on Climate Change and Adapted Land Use (WASCAL) that has funded this work, and the University of Abomey Calavi-Benin which hosts the Graduate Research Program of Climate Change and Water Resources. We are further thankful to the THY group (Max Planck Institute, Hamburg, Germany) and the Climate Service Center (Hamburg, Germany) where this work was mainly conducted and for their various supports and to Christopher Moseley (Max Planck Institute, Hamburg, Germany) for his great technical assistance.
NOTES
*Corresponding author.