Assessment of Spatio-Temporal Changes of Land Use and Land Cover over South-Western African Basins and Their Relations with Variations of Discharges

West African basins play a vital role in the socio-economic development of the region. They are mostly trans-boundary and sources of different land use practices. This work attempts to assess the spatio-temporal land use and land cover changes over three South Western African basins (Volta, Mono and Sassandra basins) and their influence on discharge. The land use and land cover maps of each basin were developed for 1988, 2002 and 2016. The results show that all the studied basins present an increase in water bodies, built-up, agricultural land and a decline in vegetative areas. These increases in water bodies and land use are as a result of an increase in small reservoirs, of dugouts and of dam constructions. However, the decline in some vegetative clusters could be attributed to the demographic and socio-economic growth as expressed by the expansion of agriculture and urbanization. The basic statistical analysis of precipitation and discharge data reveals that the mean annual discharge varies much more than the total annual precipitation at the three basins. For instance, in the entire Volta basin, the annual precipitation coefficient of variation (CV) is 10% while the annual discharge CV of Nawuni, Saboba and Bui are 43.6%, 36.51% and 47.43%, respectively. In Mono basin, the annual precipitation CV is 11.5% while the Nangbeto and Athieme annual discharge CV are 37.15% and 46.60%, respectively. The annual precipitation CV in Sassandra basin is 7.64% while the annual discharge CV of Soubre and Dakpadou are 29.41% and 37%, respectively. The discharge varies at least three times much more than the precipitation in the studied basins. The same conclusion was found for all months except the driest months (December and January). We showed that this great variation in discharge is mainly due to land use and land cover changes. Beside the hydrological modification of the land use and land cover changes, the climate of the region as well as the water quality and availability and the hydropower generation may be impacted by these changes in land surfaces conditions. Therefore, these impacts should be further assessed to implement appropriate climate services and measures for a sustainable land use and water management.


Introduction
Sub-Saharan Africa has experienced, over the past, some changes in land cover due to drought events [1,2]. After the drought of the 1980s the region has recovered from this loss. Apart from the climatic factors, other drivers of change are biological and socio-political [3] as well as demographical [4] and socio-economical [5] factors. In order to meet the need for food, fiber, water and shelter of an increasing population, the environment has seen a change in forest landscapes, farmlands, waterways and air quality. The cause of land use and land cover changes in West Africa basins was mainly attributed to rapid population growth rate [6,7]. However, this phenomenon consequently has some effects on hydrological dynamics [8,9]. The sustainability of water resources depends mainly on land use planning and management [8]. Some main West African basins, namely, Volta, Niger and Sassandra have already experienced these changes due to a rapid population growth rate and the transboundary condition of the basins.
The West African region has a high population growth rate associated with a high urbanization ratio [10,11]. Apart from urbanization, the West Africa population relies on rain-fed agriculture which requires more and more land to feed its growing population. All of these practices have consequently resulted in land use and land cover changes.
The Volta, Niger, Mono and Sassandra catchments are the main contributors to the socio-economic development of the region and are all shared between at least two countries. These basins are the main sources of the hydroelectricity production [12] of the sub-region. Despite the fact that all the riparian's countries have agriculture as their main activity, the policy for planning and management of land uses and practices may be different and this may vary within the same country. This may result in different land uses and practices within the same basin with impacts on the water availability of the river. For instance, in West Africa, the water consumption and withdrawals have increased [13] and created a strong competition between hydropower and upstream water withdrawal for agriculture [14], especially between Ghana and Burkina Faso [15,16]. Some communities and even governments tend to blame upstream countries for phenomena such as deficits in discharges or rivers, floods and water weeds [14]. Some hydropower facilities located downstream lack water because of too many withdrawals. In fact, in Volta basin, the upstream dams namely Bagre, Kompienga and Bui plants deprive Akosombo hydropower generation. Moreover, the Ghana "One village one dam" project which aims for the construction of at least 570 small reservoirs [17] for irrigation at the Northern part of Ghana of Akosombo dam could also be a real threat to Akosombo hydropower. However, the water losses caused by implementation of small reservoirs are marginal compared to the one due to climate change patterns especially evapotranspiration [18] but the implementation of large reservoirs and/or great number of small reservoirs at upstream of Akosombo could strongly affect Akosombo power generation [19].
The aim of this study is to assess the spatio-temporal changes of land use and land cover of Volta, Mono and Sassandra basins and the relation between these changes and the variations of precipitation and discharge.

Hydro Climatology of the Study Area
The studied basins are located in the southern part of West Africa (Figure 1). Sassandra catchment is located in Ivory Coast and covers an area of 75,479.79 km 2 . Mono catchment spreads across Togo and Benin, and cover an area of 24,282.26 km 2 . The Volta basin with an area of 413,520.15 km 2 is shared among Burkina Faso, Mali, Benin, Togo, Ivory Coast and Ghana. In these basins are constructed some of the largest hydropower plants which are the main sources of electricity generation of the region. The population living in these basins are mainly rural and agriculture is the main activity. The main land cover types in these basins are mainly shrubs, savanna, forest, barren land, agricultural land, built-up and water bodies. The agricultural land covers more than 20 percent of each basin area.  In the Sassandra basin, the Buyo and Soubre dams are constructed with a capacity of 165 Mega Watt (MW) and 275 MW, respectively. In the Mono basin is located the Nangbeto dam with a capacity of 65 MW. In the Volta Basin, there are five dams: the Bagre (16 MW), Kompienga (14 MW), Bui (400 MW), Akosombo (1038 MW) and Kpong (160 MW). The operation of these hydropower plants relies strongly on the hydro-climatic condition of the sub-regions.
The climate of West Africa is controlled by the movement of the Inter Tropical Convergence Zone (ITCZ) [20,21] and influenced by the Monsoon and Harmattan. Their hydroclimatic characteristics vary progressively along a north-south axis. In the southern zone or Guinean zone (transitional equatorial climate), the total annual precipitation is greater than 1100 mm shared in two rainy seasons (March-June and August-November) while the Northern part (transitional tropical climate mainly known as the Sahelian zone) is drier with a total annual rainfall variation ranging between 300 and 500 mm in one rainy season (July to September) [13]. The transitional sub-region between the Sahelian zone and the Guinea zone is called the Sudanian zone with annual precipitation increasing southward between 500-900 mm (in the Sudano-Sahelian area) and between 900 and 1100 mm (in the Sudano-Guinean area) [13]. The pan evaporation is estimated at 2540 mm per year, and an average annual runoff coefficient of about 8.3% [22]. All these factors vary from basin to basin.
For instance, in the Volta basin, the mean monthly runoff in the sub-basin varies from 623 m 3 /s at the peak of the rainy season to about 2 m 3 /s in the dry season. During the rainy season, the rainfall In the Sassandra basin, the Buyo and Soubre dams are constructed with a capacity of 165 Mega Watt (MW) and 275 MW, respectively. In the Mono basin is located the Nangbeto dam with a capacity of 65 MW. In the Volta Basin, there are five dams: the Bagre (16 MW), Kompienga (14 MW), Bui (400 MW), Akosombo (1038 MW) and Kpong (160 MW). The operation of these hydropower plants relies strongly on the hydro-climatic condition of the sub-regions.
The climate of West Africa is controlled by the movement of the Inter Tropical Convergence Zone (ITCZ) [20,21] and influenced by the Monsoon and Harmattan. Their hydroclimatic characteristics vary progressively along a north-south axis. In the southern zone or Guinean zone (transitional equatorial climate), the total annual precipitation is greater than 1100 mm shared in two rainy seasons (March-June and August-November) while the Northern part (transitional tropical climate mainly known as the Sahelian zone) is drier with a total annual rainfall variation ranging between 300 and 500 mm in one rainy season (July to September) [13]. The transitional sub-region between the Sahelian zone and the Guinea zone is called the Sudanian zone with annual precipitation increasing southward between 500-900 mm (in the Sudano-Sahelian area) and between 900 and 1100 mm (in the Sudano-Guinean area) [13]. The pan evaporation is estimated at 2540 mm per year, and an average annual runoff coefficient of about 8.3% [22]. All these factors vary from basin to basin.
For instance, in the Volta basin, the mean monthly runoff in the sub-basin varies from 623 m 3 /s at the peak of the rainy season to about 2 m 3 /s in the dry season. During the rainy season, the rainfall is greater than evapotranspiration, and vice versa, during the dry season. The basin mean, highest and lowest temperatures are, respectively, 26 • C, 44 • C and 15 • C. The relative humidity varies between 20-30% during the Harmattan period in dry season and 70-80% in the rainy season [6].
The total annual precipitation in Mono basin is 1168 mm [23]. The mean monthly discharge at Athieme station varies from 261 m 3 /s (September) at the peak of the rainy season to about 0 m 3 /s in the dry season (February to April) [24]. During the rainy season, the rainfall is greater than evapotranspiration (which varies from 4.5 to 6 mm/day [25]), and vice versa, during the dry season. Mono basin present similar temperature and relative humidity as Volta basin.
The total annual precipitation of Sassandra basin varies from 1250 to 1700 mm (Southward). The mountainous area (Man region) record the highest depth of total annual precipitation which varies from 1600-2500 mm. The potential evapotranspiration varies from 1400 to 1500 mm/year. The mean annual discharge at Soubre station is estimated at 630 m 3 /s for wet year and 475 m 3 /s for dry year [26,27]. The relative humidity and temperature are similarly to the one of Volta basin.
In the three basins, the agricultural sector represents the main socio-economic activity of the riparian following by trading. This makes all the basins exposed to the rapid land use change as the practice remains traditional with an extensive and slash farming system on burned areas.

LULC Data Sources
The data used in this study are from the United State Geological Surveys (USGS). Data of the Landsat Thematic Mapper (TM), the Enhanced Thematic Mapper (ETM) and of the Operational Land Imager (OLI) were assessed for the years 1988, 2002 and 2016 respectively to study the land use and land cover changes. To ensure the best quality, these data were extracted for the driest and cloudless months (January and February). Five (05) scenes are necessary to cover the entire Mono basin corresponding to the following paths/rows : 192/54; 192/55; 192/56; 193/54 and 193/55. To cover the entire Sassandra basin, eight (08) scenes were used corresponding to the following paths/rows: 197/ 55; 197/56; 198/53; 198/54; 198/55; 198/56; 199/54 and 199/55. The largest study basin used twenty-five (25) scenes to cover the entire Volta basin. These scenes belong to the paths from 192 to 197 with the row 52 of the 192th path, the rows from 52 to 56 of path 193, the rows from 52 to 55 of 194th path, the rows from 50 to 55 and 50 to 54 of 195th and 196th paths and finally, the row 51 to 53 of the 197th path. For a given year all these scenes were merged into a single map. The southern part of Mono basin which was cloudy was not considered. The Landsat data level 2 has been ordered for atmospheric correction. After the completion of the ordered data, it has been processing and classified using supervised classification method.

LULC Data Analysis Method
There are mainly three approaches to perform a land cover classification, namely, pixel-based, subpixel-based and object-based approaches divided in different methods each with their advantages and limitations [28]. The state of the art of each approach is well described in [28]. The maximum likelihood pixel-based classification method, which was developed in the 1970s, is the most commonly used method on Landsat images. This method has been used to develop our land cover map.
The land cover types may vary within the same basin as well as from a catchment to another catchment. For instance, for Volta basin, we have nine (09) clusters of land cover namely water bodies, built-up, barren, agricultural land, shrubs, herbaceous savanna, savanna, forest and evergreen forest. Whilst, in Mono basin we have made only six (06) clusters of land cover namely water bodies, built-up, agricultural land, herbaceous savanna, Savanna, and forest. The Sassandra basin has the same number of clusters as Mono basin, but here we have the appearance of evergreen forest cluster and the absence of herbaceous savanna class.
Indeed, the "water bodies" cluster represents the stream line, small reservoirs, lakes and dams. The "built-up" cover type represents urbanized areas and roads as well as land covered by buildings and other man-made structures, residential areas, commercial services, industrial areas, mixed urban and built up lands. The "barren" cluster was defined as land with exposed soil, sand or rocks, desert area or bare land. The "agricultural land" represents the farmland areas or the lands covered with temporary crops followed by harvest period, crop fields and pastures. The "shrubs" cluster represents some woody plants, smaller than a tree, usually having multiple permanent stems branching from or near, the ground or woody plants of relatively low height, having several stems arising from the base and lacking a single trunk. The "savanna" represents grassland with scattered trees, grading into either open plain or woodland while the "herbaceous savanna" represents grassland. The "forest" cluster refers to growth of trees and other plants covering a large area while the "evergreen forests" are forest made up of rainforest trees in tropical zone or evergreen trees.
Land classification can be subjected to errors due to geometric errors, misclassifications, and undefined classes. To statistically quantify these errors, a random selection of pixels of the classified maps was performed to build a confusion matrix [6,7]. In order to calibrate and validate the land cover classification, an accuracy assessment was performed and the kappa coefficient was used as the statistical parameter. This approach consists of the reclassification of land cover data in another way by following the step described by Congalton [29]. The accuracy assessment can be defined as a comparison of a map produced from remotely-sensed data with another map from some other source. A determination is made of how closely the new map produced from the remotely-sensed data matches the source map. It is necessary to present the accuracy assessment as an error matrix as recommended by Congalton [29]. For these assessments, the most common way to represent the classification accuracy of remotely sensed data is in the form of an error matrix. An error matrix or confusion matrix is a square array of numbers set out in rows and columns which express the number of sample units (i.e., pixels, clusters of pixels, or polygons) assigned to a particular category relative to the actual category as verified on the ground. The columns usually represent the reference data while the rows indicate the classification generated from the remotely sensed data. Reference pixels for each category were selected in a stratified random fashion. Ten-pixel times the number of land cover category are randomly selected for each group.
The confusion matrix as presented in Table 1 and consists of: (i) The Total Accuracy (overall accuracy) of the classification which is given as the sum of the truly classified pixel (Diagonal values) over the total number of the selected pixel. The total selected pixel is equal to ten times the number of land cover category. For instance, the total accuracy of Volta land cover of 1988 is calculated as follows: Then the Total Accuracy is 90.37% for Volta land cover classified data of 1988, meaning that for the land cover classification 1988 in Volta basin, 90.37% are correctly classified (Table 1). (ii) The Producer's Accuracy which is the total number of correct pixels in a category divided by the total number of pixels of that category as derived from the reference data. In other words, it is the probability of a reference pixel being correctly classified and the measure of omission error and is expressed as follow: Producers Accuracy f cluster cluster = (Corrected classi f ied pixel) (Total selected Pixel cluster ) × 100 For example, the Producer's Accuracy of water bodies (Table 1) is: Producers Accuracy f cluster water bodies = (87) (80) × 100 This means that among 90 pixels of water bodies 96.67% are truly classified as water bodies by the model and the remainder are omissions. (iii) The omission error in percentage which represents pixels that belong to the true class, but fail to be classified into the proper class. For the water bodies we have: (1 + 1 + 1) × 100/90 = 3.33%. Among the 90 pixels of water bodies 3 are omitted during the classification and regarded as built-up, agricultural land and forest. (iv) The user's accuracy, here the total number of correct pixels in a category is divided by the total number of pixels that were classified in that category as derived from classified data. It is the probability that a pixel classified on the map/image actually represents that category on the ground. The user's accuracy of agricultural land (Table 1) is: (90 × 100)/104 = 86.54%.
Meaning that among 104 pixels classified as agricultural land, it is only 90 pixels which are truly agricultural land. The additional pixels which belong to other classes is the omission error; and (v) The commission error represents the pixels that belong to another class, but are labeled as belonging to the class. The commission error of agricultural land is: (1 + 9 + 1 + 3) × 100/104 = 13.46%. The technical explanation is that 1, 9, 1 and 3 pixels of water bodies, barren land, savanna and forest respectively are classified as agricultural land.
The Table 1 exhibits the confusion matrix of Volta basin 1988 land cover classification constituted of described above factors. The producers' accuracy and the omission errors are presented on the last two raw while the user's accuracy and the commission errors are on the last two columns. As we have nine maps we just presented only Volta basin 1988 confusion matrix and the results of the remains maps was summarized in Table 2.
The kappa coefficient K, a discrete multivariate proxy used in accuracy assessment of thematic maps, is an efficient method to derive information from an image via the confusion matrix. K > 0.80 represents strong agreement and good accuracy; K between 0.40 and 0.80 means a middle accuracy, and K < 0.40 represents a poor accuracy [6,7]. The kappa coefficient is given by the Equation (4).
where r is the number of rows in the matrix, x ij is the number of observations in row i and column i, x i+ and x +i are the marginal totals of row i and column i, respectively, and N is the total number of observations. For this study, the kappa coefficient and the accuracy assessment of the three catchments for the years 1988, 2002, and 2016 are summarized in the Table 2. These factors were obtained from confusion matrix of Volta, Mono and Sassandra basin for the years 1988, 2002 and 2016. Once should recall that the kappa coefficient K > 0.80 represents strong agreement and good accuracy; K between 0.40 and 0.80 means a middle accuracy, and K < 0.40 represents a poor accuracy [6,7]. All the Kappa coefficients of this study are greater than 80% which means the results show strong agreement and good accuracy.

Precipitation and Discharge Data
The Global Precipitation Climatology Project (GPCP) monthly data [30] was used to assess the precipitation variability in Mono, Volta and Sassandra basins due to lack of observed data. The GPCP products are merged data from rain gauge stations, satellites, and sounding observations to estimate monthly rainfall on a 2.5-degree global grid from 1979 to the present [31]. The GPCP data was used to validate some model outputs within African continent and specifically in West Africa region [32][33][34][35][36][37] and one of its weakness is that it gives earlier onset rainfall compare to rainfall product of Tropical Rainfall Measuring Mission (TRMM) satellite. The GPCP data is correlated to observed rainfall in Ghana [38]. The river discharge data was collected from Bui power Authority (Bui 1954-2015), from Volta River Authority (Nawuni and Saboba 1954-2008), from Mono River Authority  and from Sassandra River Authority (Soubre 1954-2004and Dakpadou 1969-2004.
The temporal variation of some hydrological stations data at the downstream of the basins were analyzed following some literatures methods [39,40]. The daily data from the stations are described in Table 2.
These stations were chosen according to the continuity of their data (low missing data) and also base on the fact that there are mostly at the outlet of each basin or sub-basin.
However, some of these stations may be affected by the hydropower plants constructed at upstream. For instance, in the Sassandra basin, the Buyo (1980) was constructed in 1980 at upstream of Soubre with installed capacity of 165 MW. In the White Volta basin, the Bagre and Kompienga hydropower plants were constructed at the upstream of Nawuni and Saboba respectively. The dam construction could have affected the flow regime of the river and in the aim to not attribute any change in studied hydrological data only to the hydropower construction, we also choose some river stations which had not received any of this kind of activity during the considered period. It is in this view that Bui, Athieme and Dapkadou rivers were chosen respectively in Volta, Mono and Sassandra basins.
In order to have the same period study for our analysis, all the discharge data have been taken from 1969 to the end date of the data availability at each station. The consistency, the homogeneity and the stationarity of the data have been checked as well as the continuity in order to assess the missing values. There are several methods to cope with missing data. Each method has its own strengths and weaknesses [41]. The regression-based imputation derived from the single-imputation method and frequently used in hydrology for missing data correction was applied in this study [41,42]. The percentage of missing data is presented in Table 3.
The regression-based imputation method consists firstly of the estimation of regression equations that relates to the variable that contains missing data to a set of variables featuring complete information across all observations in the dataset. Secondly, the equation is used to predict the missing values [41]. This method does not have any bias problem with the mean values. The consistency, the trend, the periodicity and the persistence of the data have also been assessed as recommended by Machiwal and Jha [39]. The trend variation of both data has been checked using Mann-Kendall [43] non-parametric test to determine the significance of the observed trends test, and the Sen's Slope [44] Estimator to assess the direction (upward or downward) of the trend. The trend in both variables (precipitation and discharge) was assessed based on the period of available data, i.e., 1979-2000 for Mono basin; 1979-2004 for Sassandra basin and 1979-2008 for Volta basin. The test was also performed for the whole period of data availability from 1969 to the end data for discharge and from 1979 to 2017 for precipitation data.

Trend Detection
Several methods exist for trend detection in hydrology, namely, (i) Spearman's rho which is a rank-based test for correlation between two variables and can be used to test for a correlation between time and the data series; (ii) another rank-based test called Kendall's tau/Mann-Kendall test which is similar to Spearman's rho, but using a different measure of correlation which has no parametric analogue; (iii) Seasonal Kendall test which is a version of the Mann-Kendall test that allows for seasonality in the data; (iv) The test statistic for linear regression or the regression gradient is one of the most common tests for trend and, in its basic form, assumes that data are normally distributed and (v) other robust regression tests that are a number of robust methods for estimating trend in series and could potentially be used as alternative measures of the change [45].
The non-parametric Mann-Kendall test which is a rank-based non-parametric test for assessing the significance of a trend [43] recommended by the World Meteorological Organization (WMO) for trend detection in hydro-meteorological time series [46] is used in this study. The description of the different factors guiding this test is given below.
The null hypothesis H0 is that a sample of chronologically ordered data is independent and identically distributed. The statistic S is defined as follows: When n ≥ 40, the statistic S is asymptotically normally distributed with mean 0 and variance given by the following equation: where t is the size of a given tied group and ∑ t is the summation over all tied groups in the data sample.
The standardized test statistic K is computed by using the following equation: The standardized statistic K follows the standard normal distribution with mean zero and variance of one. The probability value p of the statistic K of sample data can be estimated using the normal cumulative distribution function as: For independent sample data without trend, the p value should be equal to 0.5. For sample data with large positive trend, the p value should be close to 1.0, whereas a large negative trend should yield a p value close to 0.0. If the sample data are serially correlated, then the data should be pre-whitened and a correction applied to calculate the variance In order to determine the magnitude of the trend, the Sen's slope is performed. We compute at the same times both the slope (i.e., linear rate of change) and the intercept according to Sen's method. The slope of a trend is estimated as follows: where β is the estimate of the slope of the trend and x j is the jth observation. An upward trend is represented by a positive value of β and a downward trend is represented by a negative value of β. These parameters are widely used in climatological and hydrological data trend analysis and the equations to compute them are also well known and well described [7,14,29,30].
In addition, the coefficient of variation of precipitation and discharge were computed in order to assess the relationship between them and land use and land cover changes.

Spatio-Temporal Changes of Land Use and Land Cover
The       For instance, in Volta basin (Table 4a) the land cover clusters such as water bodies, built-up, agricultural land and shrubs have increased between the three periods 1988-2002, 2002-2016 and 1988-2016. However, all the other land cover clusters, namely, the barren land, savanna, herbaceous savanna and evergreen forest have declined except the forest which has increased in the first period before declining in the last two periods. The decline of barren land areas could be associated with the increase of shrubs. In contrast, this could be some limitation of our model to perform the classification. It may consider some barren land as agricultural and/or built-up in the northern part of Volta basin. The increasing of shrubs for all the periods could be explained in one hand by the fact that the region has been recovered from the droughts of 1980s. In second hand, the increase in shrubs For instance, in Volta basin (Table 4a) the land cover clusters such as water bodies, built-up, agricultural land and shrubs have increased between the three periods 1988-2002, 2002-2016 and 1988-2016. However, all the other land cover clusters, namely, the barren land, savanna, herbaceous savanna and evergreen forest have declined except the forest which has increased in the first period before declining in the last two periods. The decline of barren land areas could be associated with the increase of shrubs. In contrast, this could be some limitation of our model to perform the classification. It may consider some barren land as agricultural and/or built-up in the northern part of Volta basin. The increasing of shrubs for all the periods could be explained in one hand by the fact that the region has been recovered from the droughts of 1980s. In second hand, the increase in shrubs and generally some vegetative clusters could be attributed to the tree plantation resulting from agroforestry practice. The same conclusion was made by Vittek et al., in a previous study for the period 1975-1990 for West Africa overall [47].
The result reveals that the Mono basin (Table 4b)  The Sassandra basin (Table 4c) presents an increasing in water bodies, built-up and agricultural land while the evergreen forest has declined. However, some vegetative areas such as the savanna and the forest have decline in the first period and increased for the last two periods.  Overall, the results from the LULC analysis show that water bodies, built-up and agricultural lands have increased rapidly for all studied basins. Such an increase in water bodies can be explained by the construction of small reservoirs, dugouts and dams within the basins. These basins have received more attention in term of small reservoirs and dugouts construction for irrigation, domestic water supplies and for livestock breeding as well as hydropower plants for electricity generation. Also, the observed increase in built-up and agricultural lands is associated with the increase in urbanized and cultivated areas respectively. In addition, the decline in vegetative areas (herbaceous savanna, savanna, forest and evergreen forest) can be explained by the expansion of agricultural and urbanized areas. In fact, the vegetative areas are being reduced to enlarge the land use areas or in another word, the vegetative areas are replaced by agricultural, built up and water bodies.
Over the world, and particularly in the tropics and subtropics areas, deforestation is attributed to large-scale commercial agriculture, local subsistence agriculture, infrastructure, urban expansion and mining [48]. During 2000-2010 period in Africa, 8000 km 2 , 7000 km 2 , 2500 km 2 , 2000 km 2 and 500 km 2 of forest are converted per year respectively into commercial agriculture, local agriculture, infrastructure, mining and urban expansion [48]. According to Food and Agriculture Organization report [48], in Africa, fuelwood charcoal, timber logging, livestock grazing in forest and uncontrolled fires contribute respectively to 50%, 35%, 10% and 5% of forest degradation.
It was highlighted that the main drivers of land cover change are the excessive use of trees for timber, for biomass energy namely charcoal and wood fires as well as the uncontrolled bushfire phenomenon. Wood fuel accounts for 70% of the Sub-Saharan Africa total energy production [49] and is one of the main sources of forest and woodland degradation. Due to the increase in population growth rate, and relative price changes of alternate energy sources for cooking, it is expected that the consumption of charcoal will remain at very high levels or even increase in absolute terms over the next decade [50], thus, woodland and forest degradation could steadily increase. Uncontrolled bushfire, one of the major causes of soil erosion [51] and livestock pasture destruction could also affect seriously the land cover of each basin. Unfortunately, all the basins in West Africa are exposed to this phenomenon. It has also been highlighted that global croplands, pastures, plantations, and urban areas have expanded in recent decades, accompanied by large increases in energy, water, and fertilizer consumption, along with considerable losses of biodiversity [52]. West Africa basin is under vegetative area loss due to land use practices explained above. This is confirmed by Vittek et al., [47] in their study focus on dynamic land cover change in West Africa (1975)(1976)(1977)(1978)(1979)(1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990). They also found out the decreased in all the vegetative clusters.
Lastly, the change in land cover change found is in accord with some previous work in West Africa regions [51,53] and especially in Black Volta basin [6,7].

Statistical Analysis of Precipitation Data
The Mann-Kendall test was applied to assess the trend in the precipitation data variation of each location. The statistic indicators of Mann-Kendall test are the p-value, score (S) or Mann-Kendall statistic, the Z-score, the Kendall tau and their equations were given above. So, the Kendall tau is a measure of correlation; therefore, measures the strength of the relationship between the two variables base on the rank of the data. Whilst, the p-value measures the significance of the trend and the Z-score is an indicator of the statistical significance of the trend at 95% confidence level. Indeed, at the 2-sided p-value test the rejection region of a hypothesis is Z-score ≥ 1.96 or Z-score ≤ −1.96. Whilst, the score (S) expresses the magnitude change of the trend (increasing/upward; decreasing/downward). The periods 1979-2000 for Mono basin stations, 1979-2004 and 1979-2008 for Sassandra and Volta basin stations respectively corresponding to the available discharge data period in each basin are considered for the analysis of trends of precipitation (Table 5). Table 5. Results of Mann-Kendall test of precipitation in the three basins (* represents the statistically significant trends at the 5% significance level (95% confidence level) and ** statistically significant trends at the 10% significance level (90% confidence level). BV: Black Volta sub-basin, WV: White Volta sub-basin, OV: Oti Volta sub-basin, MV: Main Volta sub-basin.). The locations of stations were chosen according to the basin size and available data. Among all the chosen locations, only three stations in the Volta basin (i.e., Goere, Porga and Saboba respectively at the northern part of Black Volta, Northern and southern part of Oti) show a significant upward trend at 95% level of confidence. The p-values of these stations are lower than the (α = 5%) at 95% level of confidence. This is also confirmed by the Z-score values greater than 1.96 which put them in the rejection zone. However, at the 90% level of confidence the Bagre Aval and Nawuni station (White Volta) show a significant upward trend. Though other locations in the Volta basin did not show any trend, the Northern to central locations indicate a positive Mann-Kendall score (S) contrarily to the southern locations (Bui, Kpong and Sunyani).

Basins and
In the Mono and Sassandra basins, the results of Mann-Kendall test reveal no trend in precipitation. Nevertheless, between 1979 and 2000, the central (Nangbeto) and the southern (Athieme and Tetetou) parts of the Mono basin present a negative Mann-Kendall score (S) contrarily to the northern part (Corre-Cope) of this basin. The same conclusion can be made for Sassandra basin  where the northern part (Odienne) shows a positive Mann-Kendall score while the region from the central (Man) to the southern part (Soubre and Dakpadou) shows a negative one.
Overall the northern part of the three basins presents a positive Mann-Kendall score while the southern part shows the contrary, though they did not present any trend. The Sassandra and Mono basins have not experienced any trend, only Volta basin present a statistically significant trend.
Where data are available, the Mann-Kendall test of the precipitation for the extended period until 2017 (1979-2017) was performed (not shown). There is no major significant difference with results of the first period presented in Table 5 apart from the Odienne (upward), Man (downward) in Sassandra basin and Sunyani (downward) in Volta basin which present a statistically significant trend for 1979-2017 period. In addition, the Mono basins do not show again any trend during the 1979 to 2017 period. significant at 90% confidence level for the first period (** in Table 6) and the second period (** in table 7), respectively.   (1979-2008)). NB: * the statistically significant trends at the 5% significance level (95% confidence level) and ** statistically significant trends at the 10% significance level (90% confidence level). The Mann-Kendall test was applied to assess the trend of the hydrological data variation of each station. The same process was followed. We analyzed our data in two times. The first time starts from 1979 (corresponding the start date of precipitation data) to end date of each station and the second time considers the data from 1969 to the end date. For both periods, all the selected hydrological stations of the basin present an increased trend in the discharge. This is confirmed by the positive Mann-Kendall tau and Z-score.

Mann-
The analysis shows that at 95% confidence level (α = 5%) only Saboba station's (* in Tables 6 and 7) trend is found to be statistically significant and this is confirmed by the Z-score (2.10884 > 1.96 for in Table 6 and 2.60479 > 1.96 in Table 7) that belong to the rejection region of 2-sided p-value test. This trend is an upward one. However, Nangbeto and Athieme trend (upward) were statistically significant at 90% confidence level for the first period (** in Table 6) and the second period (** in Table 7), respectively. Table 6. Result of discharge data Mann-Kendall test for studied basins (Mono basin (1979-2000), Sassandra basin  and Volta basin (1979-2008)). NB: * the statistically significant trends at the 5% significance level (95% confidence level) and ** statistically significant trends at the 10% significance level (90% confidence level).   , Sassandra basin  and Volta basin . NB: * the statistically significant trends at the 5% significance level (95% confidence level) and ** statistically significant trends at the 10% significance level (90% confidence level). Overall, all the downstream hydrological stations of the basins show an increase trend (non-statistically significant) in discharge during all the two analysis periods except the Saboba station which present statistically significant trend at 95% level of confidence and the Nangbeto and Athieme hydrological stations which show a statistically significant trend of 90% confidence level for the first period  and the second period . The increase in discharge of all selected stations and the statistically significant upward trend of Saboba could be caused by the increase in spatio-temporal distribution of precipitation and/or another non-climatic factor. The spatial distribution trend of precipitation and discharge for considered period in Table 6 of test are presented in Figure 6.

Stations
Despite the non-statistically significant trend in precipitation observed in Mono basin , the Nangbeto hydrological station shows a statistically significant upward trend at 90% confidence level. Though, the large part of Sassandra basin ) presents a downward (non-statistically significant) trend, the Soubre station at the downstream shows an upward trend (non-statistically significant) discharge. This could be due to land cover and land use change.
Overall, all the downstream hydrological stations of the basins show an increase trend (nonstatistically significant) in discharge during all the two analysis periods except the Saboba station which present statistically significant trend at 95% level of confidence and the Nangbeto and Athieme hydrological stations which show a statistically significant trend of 90% confidence level for the first period  and the second period . The increase in discharge of all selected stations and the statistically significant upward trend of Saboba could be caused by the increase in spatio-temporal distribution of precipitation and/or another non-climatic factor. The spatial distribution trend of precipitation and discharge for considered period in Table 6 of test are presented in Figure 6. . NB: * the statistically significant trends at the 5% significance level (95% confidence level) and ** statistically significant trends at the 10% significance level (90% confidence level).
Despite the non-statistically significant trend in precipitation observed in Mono basin , the Nangbeto hydrological station shows a statistically significant upward trend at 90% confidence level. Though, the large part of Sassandra basin (1979-2004) presents a downward (nonstatistically significant) trend, the Soubre station at the downstream shows an upward trend (nonstatistically significant) discharge. This could be due to land cover and land use change.
The present work shows an increased (non-statistically significant) trend in precipitation at the upstream of Bui and Nawuni hydrological stations of Black and White Volta basins respectively. This leads to an increase in trend of discharge at Bui and Nawuni stations. Even though the results reveal  (1979in Mono, 1979-2004in Sassandra and 1979-2008 in Volta). NB: * the statistically significant trends at the 5% significance level (95% confidence level) and ** statistically significant trends at the 10% significance level (90% confidence level).
The present work shows an increased (non-statistically significant) trend in precipitation at the upstream of Bui and Nawuni hydrological stations of Black and White Volta basins respectively. This leads to an increase in trend of discharge at Bui and Nawuni stations. Even though the results reveal an upward trend in both variables (precipitation and discharge), the magnitude of increased in discharge is greater than in precipitation. For instance, the slope of the regression line is lower for rainfall compare to the discharge. In the Oti Volta sub-basin, the Porga and Saboba stations (upstream) show an upward trend in precipitation and discharge, but the two variables are statistically significant only at Saboba station. Also, the magnitude of change upward trend in discharge (slope of 3.187) is greater the one for precipitation (slope 0.18) in Saboba station . This difference in trend and magnitude of precipitation and discharge could be due to land use and land cover change.
These present results confirm previous studies that found out the same trend in discharge at Akosombo station (Lower Volta basin) and at Bui station (Black Volta) for the period of 1970-2011 [7] [54]. In the Lower Volta basin (Ghana part), most of the rainfall stations data analyzed has experienced an increased trend as well as Akosombo discharge [54] but while in the Black Volta basin [7], several rain-gauges stations present a downward trend in precipitation (except Boura station) while the Bui discharge show an upward trend. Their results also conclude that the land use and land cover change were the probable cause of this difference in trend of precipitation and discharge [7,55]. In order to understand the influence of land use and land cover changes, a basic statistical analysis was further processed.

Basic Statistical Analysis of Precipitation and Discharge
The Table 8 presents the basic statistic namely the standard deviation (SD), the mean and the coefficient of variation (CV) expressed in percentage of precipitation (P in mm) or of discharge (Q in m 3 /s) at monthly and annually scale. The analysis of Table 8 shows that the discharges have a larger variability (greater CV) than the precipitation at annual and monthly level except some driest months namely December and January in all the basins. The annual precipitation coefficient of variation is lower (10%) in Volta basin while the three hydrological station present a greater coefficient of variation (>35%) 43.6%; 36.51% and 47.43% for Nawuni, Saboba and Bui respectively. In Mono basin, the annual precipitation coefficient of variation is lower (11.5%) while the Nangbeto and Athieme annual discharge coefficient of variation are 37.15 and 46.60% respectively. The annual precipitation coefficient of variation in Sassandra basin is 7.64% and the annual discharge coefficient of variation of Soubre and Dakpadou are respectively 29.41% and 37%. In summary, the annual discharge varies at least three times greatly than the annual precipitation. The same conclusion can also be made for all the month from February to November.
Previous studies confirmed that this large variation in discharge despite the lowest variation in precipitation is mainly due to land use and land cover change. Andreini et al. [56] showed that the coefficient of variation of runoff from the Volta sub-watersheds is much more variable than rainfall, due to land use. This has been confirmed recently in Black Volta basin with the same statistical analysis [6] and using Soil and Water Assessment Tool (SWAT) model [7]. The same conclusion has also been made in Chemoga watershed of Blue Nile basin in Ethiopia [57] using the coefficient of variation method. It is apparent from the Volta water balance that land use and land cover changes in the uplands of the basin play a pivotal role in determining the future of the basin water resources. The relative hydrological effects of forest changes and climatic variability are largely dependent on the change magnitudes and watershed characteristics [58]. Some previous work on West Africa basins showed an increase in runoff generation coefficient despite the increase in the number of dams in the basin [59] due to land use practice [60]. It has been demonstrated that such basin has loss their soil water hold capacity [60] and caused of the Sahelian paradox [61,62] which consequently result in an increase of flow of the river.
The analysis of land use and land cover change presented above reveals mostly a decrease in vegetative areas and an increase built-up areas, agricultural land, and water bodies. This change in land cover and land use also will affect the water holding capacity of the soil by decreasing the infiltration rate and increasing the runoff generation capacity of the basin. In the Black Volta basin, this increase in land use and land cover changes leads to an increase the discharge at Bui [6]. This effect combined with climate change can also reduce the underground recharge water capacity as it has been confirmed the decline in underground water recharge of Bandama basin [63]. Water availability could be at risk with changes of land use and land cover.
However, the land use and land cover changes have an important effect on the catchment runoff generation [55,[64][65][66][67][68][69][70] and potential evapotranspiration [6,71,72]. Any change in runoff generation coefficient could affect the hydrodynamic of the river. So, the land use and land cover change could affect the hydrodynamic of the river. This has been found in different basins of the region namely Bandama [63], Senegal [55,64], Mono [25,65], Sahelo-Sudanian [61] and elsewhere [66,67]. The land use practices as the construction of large dams could have an influence on surrounding climate [67] and can lead to microclimatic changes [68]. It has been proved to have high impact on the rainfall [69] and on climate system [73,74] of a given region. The study of Mara river in East Africa region exhibits that the conversion of forests to agriculture and grassland in the basin headwaters is likely to reduce dry season flows and increase peak flows, leading to greater water scarcity at critical times of the year and exacerbating erosion on hillslopes [75]. The influence of land use and land cover on hydrological dynamic of the river is well documented [25,[76][77][78][79][80]. In addition, it is also found that land use and land cover practice in the same basin could be favorable to Bui hydropower plant production under climatic uncertainties [7]. But the resultant effect of this phenomenon in other regional hydropower plants is not known and need to be investigated. So there is a nexus between land use-climate change-energy productions [73].
Beside the hydrological dynamic modification of land use [81,82], it can also be responsible of global carbon cycle and climate changes; can affect regional climate, water and air quality as well as cause biodiversity decline [52] and cause some infectious diseases. The land use change can also impact the global and regional climate [80] through the surface energy budget as well as the carbon cycle [77,83].
The land use and land cover change have strong influence on regional hydroclimate [78], ecosystem and society [77]. The land use/cover classification over Greater Horn of Africa between 1986-2000 exhibits increased conversion of large portions of forests to agriculture, particularly in Kenya, with land (acreage) under crops increasing from 6.55% in 1986 to almost 18% by the year 2000 [84]. This change in land use/cover has resulted into a modest reduction in monthly rainfall totals and also may be contributing to notable shifts in moisture convergence zones and centers of rainfall maxima. A decrease in precipitation and an increase in surface temperature around the lake region was shown when surface vegetation cover is converted to agricultural land cover type [83]. The land use/cover plays important role in climate system of the region.
The changes in land use and land cover vary according to the basin and the land use category. This has a greater impact on the streamflow generation [81], and thus on water availability for different sectors as well as hydropower plant within the basin and could vary from plant to plant within the same region. Therefore, there is a need of some careful scientific study combining the resultant effect of land use and land cover change on hydrological dynamics of each basin as well as its impacts on each hydropower plant of all the basins in West Africa under climatic uncertainties. Due to the divergent impacts of land use on the environment, it is very important to model how this factor is going to change in the future and what are the plausible effects on all the environmental aspects [82]. This factor is also very important to be taking into account while simulating climate change [85].

Conclusions
The West Africa region's basins are under different land use practices due to their trans-boundary feature. All the studied basins are subject to more intense usage of land (built-up and agricultural land) due to population growth rate and socio-economic drivers. Between 1988 and 2016, the annual water bodies (built-up) has increased by 4.8% (13.4%); 8.2% (3.33%) and 5.23% (4.56%) for Volta, Mono and Sassandra basins, respectively, while the agricultural land has increase by 20.27%, 5.57% and 4.19%. The expansion land use area (built-up and agricultural land) has led to the reduction in vegetative areas namely forests and savannas. The savanna (herbaceous savanna) has decreased by 1.51% (2.14%), 0.43% (0.65%) per year for Volta and Mono basins respectively while an increase (0.02%/year) has been observed in savanna within Sassandra basin. The forest has decreased by 2.58% and 3.43% per year for Volta and Mono basin respectively while it has increased by 0.28% per year in Sassandra basin. The evergreen forest has also decreased by 3.45% and 3.36% per year in Volta and Sassandra basins respectively. Nevertheless, the shrubs cluster has increased over time within Volta basin and this could be attributed to the agroforestry practice in the basin. This change could be attributed to the cause of different trend and magnitude in precipitation and discharge observed in the study basin. Thus, the change in land use and land cover could have an influence on runoff generation coefficient of the basin as well as the hydrological dynamic modification. It has been responsible of the decline of the underground water recharge of the basin in West Africa [7,78]. This can also influence surrounding regions [54]. The analysis of the discharge in all selected hydrological stations of the basins reveals an upward trend and this can be due to the spatio-temporal precipitation variability in the basins or the change in land use or land cover. However, the spatio-temporal variation of the precipitation presents a downward trend (non-statistically significant) from the southern to central part of Mono and Sassandra basins and in southern part of Volta basin. Whilst, the Northern part of Mono and Sassandra basins as well from the North to central part of Volta basin present an upward trend (non-statistically significant). Few of the precipitation sites shows a statistically significant trend and this is observed only in Volta basin. Despite this trend all the hydrological station presents an upward trend, only Saboba present a statistically significant at 95% level of confidence and Nangbeto and Athieme stations (Mono basin) at 90% level of confidence for 1979-2000 and 1969-2000 periods respectively. The magnitude of upward trend in discharge is greater than the precipitation one within the same period.
In addition, the coefficient of variation of precipitation and discharge analysis for the three basin reveals that the discharge varies much more (at least three times) than the annual precipitation at the three basins. The same conclusion was made for all the month except the driest months (December and January). This difference in coefficient of variation as well as in trend and magnitude is mainly due to land use and land cover change.
Although the change in land use (conversion of vegetative land into agricultural land) has led to increase in discharge, its resultant effects such as sedimentation and siltation could affect the water quality and availability for all water users' sectors. Thus, all the water user sectors could be affected by this phenomenon under this climatic uncertainty. For instance, it has been demonstrated in Black Volta basin that this change will be favorable for Bui hydroelectric production in the short and medium terms. This impact in the long term is not known and its resultant effects on the existing West Africa hydropower plant needs to be investigated for water-energy security nexus issues.
Finally, there should be a common land use planning and management for each basin for sustainable land use, agricultural practice and energy production.