A New Statistical Approach to Assess Climate Variability in the White Bandama Watershed, Northern Côte d’Ivoire ()
1. Introduction
Climate change is becoming nowadays an important center of interest for the scientific community and also for the world population. Results from numerous researches confirm that no part of the planet is spared the hydro-climatic changes. In fact, this variability sometimes unpredictable of climate by the earth surface becomes recurring since the 1950s through rainfall deregulation, temperature increase and the occurrence of floods and droughts [1] [2] . In West Africa, drought conditions persistence has been observed during the last decades. The most recognized effect of climate interannual variability in this region is the severe drought that occurred during the 1970s and 1980s, with harmful consequences to population well-being and environment [3] [4] . Particularly, for the Sudano-Sahelian Africa, the rainfall deficit has been around 13% and 27% respectively during the 1970s and 1980s decades [5] .
Some authors have examined rainfall trend in Côte d’Ivoire in order to understand the effects of climate variability. Results showed decreasing trend in precipitation and streamflow [6] [7] [8] [9] . The white Bandama watershed, in northern Côte d’Ivoire, belongs to the large Sudano-Sahelian area in West Africa. The watershed follows the general climatic trend observed in West Africa. Several studies have assessed the variability of climate conditions during the second part of the 20th century [3] [10] [11] . They reached opposite findings. [12] suggested that the dryness has persisted through the 1990s with a few rainy years. In contrast, other studies [10] [13] showed that the Sahelian region has recorded a wet period since the beginning of the 1990s and consequences of this phenomenon are floods recorded in some cities [14] [15] . Unfortunately, most of these studies, both ancient and recent, are limited to the 2000s, because the availability of reliable data from measurement network has gradually decreased and ceased with the advent of political crisis occurred in the country in 2002. Moreover, statistical methods used, generally, are Pettitt, Buishand, Lee and Heghinian and Hubert statistical tests, implemented inside the Khronostat software [16] . In many areas around the world, studies undertaken about detecting trends and breaks in climate variables were made by using the Mann-Kendall and the Student’s t test [17] - [31] .
The main goal of this paper is to (re)assess climate and hydrological trends which occurred in the white Bandama basin during the period 1950-2000 and to identify their magnitudes. A methodological framework including several statistical tests is implemented, using the standardized index method to detect trends in long term hydro-climatic data, the Mann-Kendall statistical test to know the level of significance of these trends, and the Kendall-Theil Robust Line statistical test to highlight their magnitude. Finally, the Student’s t statistical test is used to determine the break years.
2. Study Area
The white Bandama basin in northern Côte d’Ivoire, is part of the Bandama river basin (Figure 1). It is located between 9˚22' and 10˚26' North latitude and 5˚00' et 6˚30' West longitude. It covers Ferkessédougou, Ouangolodougou, Dikodougou, Korhogo, M’bengué, Sinematiali, Boundiali and Niakaramadougou cities. Its area is estimated to 12,750 Km2 and population is around 1,607,497 people. The population predominantly relies on agriculture and livestock breeding. Because of its strategic geographical position (borders with neighbouring countries), the watershed is the center of several economic and agricultural activities and has a high demographic growth [32] .
The watershed river network includes the Bandama river and its tributaries (Solomougou, Lokpoho, Lafigué, Badenou). Its climate is classified as tropical regime of transition [33] , with one rainy season from May to October and one drought season from November to April. The annual rainfall average is evaluated to 1230 mm (period 1950-2015). The runoff regime depends on the contrasted climatic regime. There is one period of low runoff from December to May. During that period, runoff generally is from subsurface return flow. The second runoff period is from June to November and represents the period of high flow. The most important floods occur in September with 300 m3/s on the 1962-1997 period [11] . Average relative humidity varies between 35% and 79%. Insolation values are spread over 160.6 hours (July) to 273.8 hours (January). Mean annual temperature is 26.6˚C. The highest values are obtained during the drought season with a peak in March (29.5˚C) and the lowest values are during rainy season with minimum in August (24.7˚C). All these climate information are from Korhogo synoptic station (SODEXAM, 2000) [34] .
3. Data and Methods
3.1. Hydro-Climatic Data Set
Hydro-climatic data used in this work cover the entire study area at the monthly
scale and result from several entities. One part of climate data over the period 1950-2000 was provided by the Society of Development and Exploitation, Aeoroportuary, Aeronautic and Meteorology (SODEXAM): these are monthly rainfall at Boundiali, Ouangolodougou, Niofoin, Niellé, M’bengué, Sinématiali, and Korhogo stations (see Figure 1 for location, Table 1) and monthly temperature at Korhogo synoptic station. Monthly rainfall and temperature of Ferkessédougou station are from the meteorological station of the African Sugar Company of Côte d’Ivoire (SUCAF CI) and were recorded from 1944 to 2015. Korhogo synoptic station was used for temperature data and to estimate potential evapotranspiration by the Thornthwaite formula [35] . Taking into account, temperature and evaporation low spatial variability, information provided from Korhogo synoptic station are considered to be sufficient for the entire studied watershed [33] .
Discharge data are from the network of measuring stations managed by the Ministry of Water and Forest. That data set is from Badikaha station located at the basin outlet and several surroundings stations: Tortiya, Tawara, Seguekiélé, and Sirasso (see Figure 1 for location, Table 2).
3.2. Rainfall and Discharge Data Correction
For this study rainfall data from Boundiali, Ferkessédougou, Korhogo and Ouangolodougou stations were used because they present a long time serie and fewgaps. When it was necessary to fill the gaps, data from other ten neighbour stations which operated during the missing months were used. The method of
Table 1. Rainfall stations used for this study.
Table 2. Discharge stations used for this study.
replacing the missing value by a weighted average of station trends monthly mean was also used. To fill gaps in discharge data got from Badikaha hydrometric station, the test of similar proportionality between hydrometric stations located on the same river was used. In fact, we applied the linear regression and correlation test between the main hydrometric Badikaha station and the neighbour discharge stations from Tortiya, Tawara, Seguekiélé, and Sirasso [36] .
3.3. Statistical Tests for the Analysis of Climate Trend
The methodological framework is summarized in Figure 2. At first, in order to detect change in of database, the standardized variable statistical test was selected. Secondly, in order to understand trends in long time data, Mann-Kendall statistical test was used. Thirdly, the magnitude of trends was detected by calculating the slope of the line. Fourthly, change-point in hydro-climatic data were detected by the Student’s t test. The Mann-Kendall and Student’s t tests were processed with the computer software TREND [37] . To determine the magnitude of trends, we employed the line slope statistical test, by using the Kendall-Theil Robust Line (KTR Line-version 1.0) software developed by the US Geological survey [38] . These visual basic programs-available via www.toolkit.net.au (accessed on the 12/07/2018) are designed to facilitate statistical testing for trend, change and randomness in hydrological and other time series data. TREND has 12 statistical tests, based on the WMO/UNESCO Expert Workshop on Trend/Change detection and on the CRC for Catchment Hydrology publication Hydrological Recipes [39] .
3.3.1. Standardized Variable
The standard or standardized variable is a positive or negative number of standard deviations by which the value of an observation, data point or variable is above or below the mean value of what is being observed or measured. Variables above the mean are positive, while variables below the mean are negative. The standardized variable of any climate variable can be defined by [40] [41] :
where
is the standardized variable, X is the annual mean of the variable,
is the annual mean of the serie X and
is its standard deviation.
Figure 2. Study methodological framework.
3.3.2. The Mann-Kendall Statistical Test
The non-parametric Mann-Kendall test [37] [42] [43] [44] [45] [46] is commonly employed to quantify the significance of trends in series of environmental, climate or hydrological data. The null hypothesis, H0, is that there is no trend in data. The alternative hypothesis, H1, is that there is an increasing or decreasing trend. This method tests whether there is a trend in the time series data. It is a non-parametric test. The n time series values
are replaced by their ranks
(starting at 1 for the lowest up to n). The test statistic S is:
where sgn (x) = 1 for x > 0, sgn (x) = 0 for x = 0, sgn (x) = ?1 for x < 0.
Positive and negative values of Z indicate increasing and decreasing trends respectively. Testing trends is done at the specific a significance level. When p-value < α, the null hypothesis is rejected and a significant trend exists in the time series. P-value is obtained from the standard normal distribution table. In this study, significance levels
and
were used. At the 5% significance level, the null hypothesis of no trend is accepted if p-value > 0.05 and rejected if p-value < 0.05. At the 10% significance level, the null hypothesis of no trend is accepted if p-value > 0.1 and rejected if p-value < 0.1 [21] .
3.3.3. The Kendall-Theil Robust Line Statistical Test
This test [38] allowed us to calculate parameters for robust, non-parametric estimates of linear-regression coefficients between two continuous variables. This method is based on a linear regression model which is written as:
where
is the explanatory (independent, predictor, or X) variable for each data point (i);
is the response (dependent, predicted, or Y) variable for each data point (i);
is the residual error or uncertainty in the predicted Y value for each data point (i); m is the estimated slope; b is the estimated intercept; n is the number of XY data points in the sample.
The slope of the line (m) is estimated as the median of all pairwise slopes between each pair of points in the data set. Each individual slope estimate
for the line connecting the ith and jth data point is calculated by use of the equation [47] - [52] :
3.3.4. The Student’s t Statistical Test
This method allows us to underline a change-point by checking whether the means in two different periods are different. The test assumes that the data are normally distributed. The Student’s t test statistic t is (critical test statistic values for various significance levels can be obtained from Student’s t statistic tables) [53] [54] [55] :
where:
and
are the means of the first and second periods respectively; m and n are the number of observations in the first and second periods respectively; S is the sample standard deviation (of the entire m and n observations).
4. Results
4.1. Interannual Variation in Hydro-Climatic Data
By calculating the standardized index of climatic variables at different time scales, the evolution of annual rainfall within the basin was determined and a wet period, a dry period and a normal period could be highlighted (Table 3).
A wet period is considered when the indexes are negative and located below the x-axis shown by zero and vice versa for a dry period. Thus, normal period is when indexes are distributed in a balanced way on either side of the x-axis.
Results allowed us to show three main trends more or less similar at the annual scale. From 1950 to 1970, climate was characterized by wet conditions. For Ouangolodougou station, this wet period stopped in 1960. This wet trend is confirmed by rainfall index around +2 and mean annual discharge (135 m3/s) greater than the long term serie mean annual discharge (83 m3/s). Also the mean annual rainfall (1500 mm) is greater than the long-term mean annual rainfall (1200
Table 3. Wet, dry and normal periods characterization from 1950 to 2000.
mm). Then, exceptionally, from 1961 to 1974, Ouangolodougou rainfall station recorded a normal period with rainfall index around 0. During the period with normal conditions, some dry years were observed (1958, 1961, and 1967) over the study basin. Finally, from 1970 to 2000, climate was characterized by dry conditions, with low water, demonstrated by negative rainfall and discharge indexes. During that dry period, the mean annual rainfall was around 1000 mm and the mean annual discharge (46 m3/s) was lower than the long-term annual mean discharge (83 m3/s). The year 1990 was the rainiest at Korhogo station (Figure 3).
Temperature and potential evapotranspiration indexes are under 0 from 1950 to 1982, demonstrating a colder period (26.4˚C) than the region annual mean (26.7˚C). Then, these two climate variables got an upward trend until 2000 and the mean annual temperature was 27˚C (Figure 4).
The results obtained by the standardized index method underline that in the whole basin there has been a dry period since the 1970s, marked by rainfall and river discharge downward trend. At the opposite, an upward trend for temperature and potential evapotranspiration is observed since that period.
4.2. Trend Significance and Magnitude in Hydro-Climatic Data
The Mann-Kendall and Kendall-Theil Robust Line tests were applied to rainfall, temperature, potential evapotranspiration and discharge long-term data at annual, seasonal and monthly time scale (Table 4 and Table 5). The Mann-Kendall test showed that there is no significant statistical trend in annual rainfall for all the watershed’s stations at level of significance α = 5%. At this same time scale, only Boundiali and Ferkessédougou stations have recorded a downward trend at level of significance α = 10% with a weak magnitude around 0. By contrast, the rainy season has declined across the basin, except at Ferkessédougou station, where there has not been a significant trend at the level of α = 5% and a trend magnitude around 1. For the rainy season (from May to October), different
Figure 3. Rainfall interannual variation in the white Bandama basin at (a) Boundiali; (b) Ferkessédougou; (c) Korhogo and (d) Ouangolodougou station.
(a)(b)(c)
Figure 4. (a) Temperature; (b) Potential evapotranspiration; (c) Discharge interannual variations, respectively at Korhogo and Badikaha station.
Table 4. Results of rainfall trends at annual, seasonal and monthly scale (1950-2000).
Table 5. Results of temperature and potential evapotranspiration trends at annual, seasonal and monthly time scales (1950-2000).
trends were observed across the basin’s stations. Rainfall in August, the rainiest month, has decreased at the level of significance α = 5% and α = 10%. Nevertheless, at Boundiali station, there has been no trend for August rainfall and the magnitude was −1.95. In Ferkessédougou station, by considering α = 5% significance level, there is only a significant trend decline in August rainfall. At 10% significance level, the decline trend concerns annual, rainy season, and August, September and October rainfall. The Sen’s slope for August rainfall is the highest in the long term series (−2.65). Stations of Korhogo and Ouangolodougou present the same trends in the annual precipitation, rainy season and August. These two stations have no significant trend in annual precipitation at 5% and 10% significance level but they record a same significant decline trend in the rainy season at both significance levels and almost the same Sen’s slope value. It means −0.13 and −0.11 for Korhogo and Ouangolodougou respectively. At Korhogo synoptic station, temperature and potential evapotranspiration have an upward trend, except for November, in which, there is no significant trend at the both significance level (5% and 10%).
Discharge is the only parameter which presents a unique trend at all time scales across the basin. It has a significant decline trend at 5% and 10% at the annual scale, for the rainy season and the rainy months. September and October discharge have the highest values for the Kendall-Theil Robust Line test, respectively −10.85 and −12.02. According to the geographical location through the basin (Figure 5), rainfall has a trend or not at 5% and 10% significance level for any time scale (annual, seasonal and monthly), while temperature and potential evapotranspiration have increased over the basin. At the basin outlet, located in the south, discharge showed a unique decline tendency at 5 and 10% significance level.
4.3. Change-Year Detection in Hydrometeorological Series
Student’s t test showed several break points. A change point was revealed in the long-term precipitation data in 1970 for Ouangolodougou station, in 1971 for both Ferkessédougou and Korhogo stations, and in 1975 for Boundiali station (Table 6). Student’s t test statistic method also showed a break year in 1982 for both temperature and evaporation in the basin (Table 7).
Discharge data in the basin exhibit alternation from abundance to paucity.
(a) (b) (c)
Figure 5. Spatial trends of hydro-climatic variables at annual (a); seasonal (b); monthly (c) time scale by the Mann-Kendall statistical test from 1950 to 2000.
Table 6. Detection of change-year in precipitation data at each station by Student’s t test over the period 1950-2000.
Table 7. Detection of change year of temperature, evaporation and discharge across the catchment by Student’s t test.
River flow tended to decrease over the period 1962-1997. The change point was recorded, by using Student’s t test statistic method, in 1977 (Table 7).
5. Discussion
The study firstly applied the standardized index method in order to determine the interannual variability of hydrometeorological parameters over the period 1950-2000. This test showed a decline trend for rainfall and discharge data, whereas temperature and potential evapotranspiration had an upward trend. A wet period was highlighted before the 1970s, a break/change period from the 1970s to the 1980s, and a dry period from the break to the 1990s. Our findings are in agreement with several studies conducted in Côte d’Ivoire [11] [34] [56] [57] [58] [59] [60] and into Sahelian part of West Africa, that show a drought from the 1970s to 1993 [1] [3] [12] [61] [62] [63] . According to these authors, precipitation variations in the Sudano-sahelian region of West Africa, based on the rainfall standardized index, indicate that the 1970s and 1980s were dry due to fluctuations in the Intertropical Convergence Zone (ITCZ) shifting. The ITCZ is defined as a convergence zone of the northeasterly Harmattan winds that originate in the Sahara and the southwest monsoon flow that emanates from the Atlantic [64] .
The second step of the analysis was to understand how the hydrometeorological trend could be considered in the study basin. Both the Mann-Kendall and Kendall-Theil Robust Line statistical tests were applied over the period from 1950 to 2000. Despite the severe drought during the 1970s and the 1980s, the annual rainfall downward trend, for the period 1950-2000, is not statistically significant. According to some authors, this situation could be explained as an annual rainfall recovery in West Africa during the last decade largely attributed to a climate shift in the mid-1990s to prevailing positive phases of the Atlantic Multidecadal Oscillation [13] [65] [66] [67] . To corroborate this assumption, they relate to floods that occurred in some urban areas of the Sahel. In fact, Korhogo and Bondoukou, two mains cities of the study basin, were devasted by torrential rainfalls in 2003 and 2006 [15] . Also, Niamey was deluged by heavy rainfalls in 1998, Dakar, Saint-Louis and Kaolack in Senegal were invaded by high floods in 1999 and again in 2000 [13] [14] .
As the annual rainfall trend is not significant, on the contrary, trends in the rainy season and monthly rainfall during the rainy season are strongly decreasing at the level of significance α = 5% and α = 10%. These results were also found by several authors, particularly, in Bangladesh [51] , in stream flow distributions in Pacific Northwest of the United States [68] , in East Africa [69] , and in Nigeria [70] . They explain, on the one hand, that discrepancies between annual rainfall trends and rainy season could mean that the main factors of trends in annual and seasonal rainfall may be different. Also, effects of seasonal rainfall trend can be hidden by the overall trend of annual rainfall [21] . Different trends between annual rainfall and seasonal time scale could also be explained by the quality of data. Rainfall data from stations in the watershed contain gaps. The filling of these shortcomings can spoil the homogeneity of records.
Temperature and potential evapotranspiration significant rising trend at 5% and 10% significance level for annual, seasonal and monthly scale highlight that the study area has got higher during the last decades. People in rural zones of West Africa support themselves mainly from subsistence level farming; hence, climate warming threatens them over the long term. Temperature rising across West Africa is described in many researches. According to [71] , by using remote sensing data, temperature has increased of 0.5˚C to 0.8˚C from 1970 to 2010, with a greater magnitude of change in the 1990s. Several reasons explain this temperature warming in the Sudano-sahelian area. At first, generally around the world, the concentration of anthropogenic greenhouse gas has increased in the atmosphere [72] [73] [74] . Secondly, anthropogenic pressure is becoming higher. Especially, across sub-Sahara Africa, [75] showed that from 1975 to 2000, there has been a 57% increase in agricultural areas and 15% increase in barren (largely desert) areas was accompanied by a 16% decrease in total forest cover.
Discharge trend was estimated only at the watershed outlet because of homogenous and long-term data lacks. Results display a significant decline trend over the study period at annual, seasonal time scale and during the rainy season. The statistical tests were significant at the level of α = 5% and 10%. Moreover, a break year was detected in 1977, which separates discharge long-term data according to two different periods by annual module: 100 m3/s over the wet period before 1977 and 35 m3/s over the dry period after that year. Research conducted by [11] at Tortiya hydrometric station (see Figure 1 for location) also showed a discharge decline, with a change year in 1974 [4] . Also found a decline in the Baniriver discharge in Mali for the period 1950-2000. Break years found with the Student’s t statistical test in hydrometeorological data at annual scale are in compliance with other authors’ studies about climate variability across Côte d’Ivoire [36] [76] [77] [78] . These breaks coincide with the period from 1970s to 1980s, marked by a severe drought in West Africa. Spatial distribution of hydro-climatic stations displaying upward, downward and no trend at annual, seasonal and monthly time scale from 1950 to 2000 is presented in Figure 5. The synoptic station in the center zone of the studied area presents an upward trend at the significance level of 5% and 10% for temperature and potential evapotranspiration. Also, the flow station, in the southern zone (outlet) of the studied watershed records a downward trend at the significance level of 5% and 10% for river discharge. Lastly, the very small number of rainfall stations and their weak geographical distribution do not allow concluding about specific regional rainfall trends across the watershed.
The discharge declining trend is stronger than the rainfall trend. Consequently, we conclude that discharge tendency is only partly influenced by precipitation trend. Hence, human activities and high potential evapotranspiration could be the main causes of discharge decline. Similar results were found in other areas of the world [79] [80] [81] , for instance in China. Discharge decline due to anthropogenic pressure could be explained by the large number of projects carried out across the study area in the 1980s and 1990s for agricultural activities, water supply, fight against waterborne diseases, and development of livestock [7] [82] . These development programmes allowed slowing down the exodus of young people and rural inhabitants towards big and developed cities of Côte d’Ivoire by improving their life level but the disadvantage was the increase of human pressure on water resources and environment.
6. Conclusion and Perspectives
The objective of this study about the white Bandama basin in the northern Côte d’Ivoire was focused on hydro-climatic trends and their magnitude. A methodological framework including several statistical tests was implemented: 1) the standardized index allowed us to identify trends into hydro-climatic data; 2) the non-parametric Mann-Kendall test was used to quantify the significance of trends in climate variable data; 3) the Kendall-Theil Robust Line statistical test permitted us to know the magnitude of trends; 4) and the Student’s t test allowed us to detect break/change years in climate and hydrological data. The downward trend in annual rainfall for the period 1950-2000 is not statistically significant, while temperature and potential evapotranspiration have registered a significant increasing trend. Spatially, the Badikaha hydrometric station in the basin outlet in the southern zone was submitted to significant decreasing trend in discharge. Korhogo synoptic station in the center of the study area has displayed significant increasing trend for temperature and potential evapotranspiration. All these changes in climate and hydrological data started from the 1970s to the 1980s. We conclude that anthropogenic pressure and high evapotranspiration were the main driving factors of discharge downward trend during the period 1950-2000.
Acknowledgements
This research was supported by the Swiss Confederation through the excellence scholarship for foreign students obtained by Franck Zokou YAO.