Are the Fouta Djallon Highlands Still the Water Tower of West Africa?

: A large share of surface water resources in Sahelian countries originates from Guinea’s Fouta Djallon highlands, earning the area the name of “the water tower of West Africa”. This paper aims to investigate the recent dynamics of the Fouta Djallon’s hydrological functioning. The evolution of the runo ﬀ and depletion coe ﬃ cients are analyzed as well as their correlations with the rainfall and vegetation cover. The latter is described at three di ﬀ erent space scales and with di ﬀ erent methods. Twenty-ﬁve years after the end of the 1968–1993 major drought, annual discharges continue to slowly increase, nearly reaching a long-term average, as natural reservoirs which emptied to sustain streamﬂows during the drought have been replenishing since the 1990s, explaining the slow increase in discharges. However, another important trend has been detected since the beginning of the drought, i.e., the increase in the depletion coe ﬃ cient of most of the Fouta Djallon upper basins, as a consequence of the reduction in the soil water-holding capacity. After conﬁrming the pertinence and signiﬁcance of this increase and subsequent decrease in the depletion coe ﬃ cient, this paper identiﬁes the factors possibly linked with the basins’ storage capacity trends. The densely populated areas of the summit plateau are also shown to be the ones where vegetation cover is not threatened and where ecological intensiﬁcation of rural activities is ancient.

The Fouta Djallon highlands are unanimously recognized as the 'water tower' of West Africa as they constitute the source of many transboundary rivers and their tributaries in the region. Table 1 gives the dependence rate of some countries on the FD rivers. In this paper, the acronym 'FD' refers to the Fouta Djallon and the broader Guinean Dorsale ( Figure 2) that runs southeastwards from the Fouta Djallon mountains, in Guinea, to the Nimba mountains, on the border with Côte d'Ivoire. By contrast, the acronym 'FDss' refers to the Fouta Djallon stricto sensu, which only constitutes the northwestern side of the Guinean Dorsale. Table 1. Dependence rate (%) from water coming of the Fouta Djallon highlands of the concerned West African countries [5]. The range gathers the sources of the three main West Sudano-Sahelian rivers: the Gambia River, the sources of the Bafing-Falémé complex which form the Senegal River after the confluence of the Bafing and the Bakoye in Mali, and those of the Niger River and its Guinean tributary on the left bank, the Tinkisso [2]. The range also feeds the two main rivers of the Atlantic side of Guinea, the Kouba/Koumba and the Tominé, which become the Koliba and then the Rio Corubal in Guinea Bissau, as well as the Konkouré complex, and the main rivers crossing Sierra Leone and parts of Liberia. The Fouta Djallon highlands are unanimously recognized as the 'water tower' of West Africa as they constitute the source of many transboundary rivers and their tributaries in the region. Table 1 gives the dependence rate of some countries on the FD rivers. In this paper, the acronym 'FD' refers to the Fouta Djallon and the broader Guinean Dorsale ( Figure 2) that runs southeastwards from the Fouta Djallon mountains, in Guinea, to the Nimba mountains, on the border with Côte d'Ivoire. By contrast, the acronym 'FDss' refers to the Fouta Djallon stricto sensu, which only constitutes the northwestern side of the Guinean Dorsale. Table 1. Dependence rate (%) from water coming of the Fouta Djallon highlands of the concerned West African countries [5].

. What Does Not Change
The geological context of the Fouta Djallon differs from a natural water tower. The entire range can be divided into a Precambrian and mostly granitic basement in the south east and a sedimentary ensemble in the north west. The latter is mostly composed of sandstone with doleritic sills and intercalations of argillites ( Figure 2). The separation follows a virtual Conakry-Bamako line [6]. In the FDss, the granite is widely covered by hundreds of meters of sandstone, dolomites, and argillites, with frequent sills of dolerites, none of which constitutes a good natural reservoir. Mamedov et al. (2010) [7] estimated the sandstone-quartzite-argillite thickness as ranging from 1500 to 4000 m. The FDss is itself divided into three geological contexts: a small granitic area southward, the western sandstone area (soft sandstone and shale with some injections of dolerites), and the northwestern doleritic FD, which is also mostly composed of sandstone and shale, but which features a large number of outcropping sills and dolerite veins. Boulvert [6] noticed that these thick sandstone plateaus are separated by the incisions (of more than 700 m) of the so-called 'Atlantic southern rivers' (i.e., the Kokoulo, the Kakrima, and the Koumba), highlighted by spectacular sandstone cornices, gorges, and waterfalls. The Guinean Dorsale is mostly an Archean granitic basement with few sills of dolerites and few outcrops of gneiss and quartzite [2]. All the Upper Niger basin and a part of the Konkouré River basin constitute a granitic basement, locally covered by sandstone and the typical 'lateritic pan'. The common characteristics of the whole mountain area of Guinea is the great extension of the 'bowe', i.e., lateritic plateaus composed of ferruginous and overall aluminous crusts [8].
However, no clear difference can be noticed in their hydrological behavior, most probably due to their similarly high levels of imperviousness. Infiltration exists at the FD scale (faults, cracks, and fractures, very few permeable rocks, local shale, and green rocks), but not at the sample scale. Most of the rocks are either completely impervious (granite, dolerites, quartzite) or mainly impervious (sandstone), except their alterites. In the entire FD, the rocks have a very low water storage capacity [2]. However, Boulvert [8] highlighted the extension of "paleo-crypto karsts" in the Bowe area of Guinea, which covers most of the Fouta Djallon. Therefore, the role of possible occulted streamflows remains to be analyzed.
The main characteristics of the FDss, which are common to the entire FD, can be summarized by the following [9]: The range gathers the sources of the three main West Sudano-Sahelian rivers: the Gambia River, the sources of the Bafing-Falémé complex which form the Senegal River after the confluence of the Bafing and the Bakoye in Mali, and those of the Niger River and its Guinean tributary on the left bank, the Tinkisso [2]. The range also feeds the two main rivers of the Atlantic side of Guinea, the Kouba/Koumba and the Tominé, which become the Koliba and then the Rio Corubal in Guinea Bissau, as well as the Konkouré complex, and the main rivers crossing Sierra Leone and parts of Liberia.

What Does Not Change
The geological context of the Fouta Djallon differs from a natural water tower. The entire range can be divided into a Precambrian and mostly granitic basement in the south east and a sedimentary ensemble in the north west. The latter is mostly composed of sandstone with doleritic sills and intercalations of argillites ( Figure 2). The separation follows a virtual Conakry-Bamako line [6]. In the FDss, the granite is widely covered by hundreds of meters of sandstone, dolomites, and argillites, with frequent sills of dolerites, none of which constitutes a good natural reservoir. Mamedov et al. (2010) [7] estimated the sandstone-quartzite-argillite thickness as ranging from 1500 to 4000 m. The FDss is itself divided into three geological contexts: a small granitic area southward, the western sandstone area (soft sandstone and shale with some injections of dolerites), and the northwestern doleritic FD, which is also mostly composed of sandstone and shale, but which features a large number of outcropping sills and dolerite veins. Boulvert [6] noticed that these thick sandstone plateaus are separated by the incisions (of more than 700 m) of the so-called 'Atlantic southern rivers' (i.e., the Kokoulo, the Kakrima, and the Koumba), highlighted by spectacular sandstone cornices, gorges, and waterfalls. The Guinean Dorsale is mostly an Archean granitic basement with few sills of dolerites and few outcrops of gneiss and quartzite [2]. All the Upper Niger basin and a part of the Konkouré River basin constitute a granitic basement, locally covered by sandstone and the typical 'lateritic pan'. The common characteristics of the whole mountain area of Guinea is the great extension of the 'bowe', i.e., lateritic plateaus composed of ferruginous and overall aluminous crusts [8].
However, no clear difference can be noticed in their hydrological behavior, most probably due to their similarly high levels of imperviousness. Infiltration exists at the FD scale (faults, cracks, Water 2020, 12, 2968 6 of 31 and fractures, very few permeable rocks, local shale, and green rocks), but not at the sample scale. Most of the rocks are either completely impervious (granite, dolerites, quartzite) or mainly impervious (sandstone), except their alterites. In the entire FD, the rocks have a very low water storage capacity [2]. However, Boulvert [8] highlighted the extension of "paleo-crypto karsts" in the Bowe area of Guinea, which covers most of the Fouta Djallon. Therefore, the role of possible occulted streamflows remains to be analyzed.
The main characteristics of the FDss, which are common to the entire FD, can be summarized by the following [9]: -the total absence of a generalized water table and of continuous groundwater reservoir [9].
There are only small and localized aquifers, well known by the local people [10]. -the very low mineral concentration of water (lower than 50 mg L −1 ). The average of dissolved load in total sampling observed by Orange (1990a) [9] is only 35 mg L −1 , due to the sandstone and/or granite predominance and the importance of stable lateritic crusts [9]. This relatively old finding remains overall coherent and in line with a series of small sampling tests made in the FDss by the authors in March and May 2019. In 28 samples taken at an elevation of less than 800 m, the mineral concentration of water was on average 55 mg·L −1 (standard deviation = 74, with values ranging from 6.5 to 414), and in 14 samples taken at an elevation higher than 800 m, the mineral concentration was on average 24.5 mg·L −1 (standard deviation = 14.8, values from 6.5 to 53). -the stationary geochemical regime: the alteration profile deepening (4.6 mm·yr −4 ) is slightly higher than soil erosion (4.2 mm·yr −4 ) [9,10].
The FD's overall impermeability strongly impacts surface water behavior. All of the geological context display the same functioning: very weak discharges during the dry season (from December to June), because they are not supported by groundwater. In contrast, at the occurrence of the first monsoon rainfall events (in June or July), discharges increase suddenly to peak in September, with the rainfall peak [2].
The impervious geology of most of the FD mountains [9,11] leads to high runoff and recession coefficients in the main rivers, due to a very low soil water-holding capacity on the one hand and due to the above-mentioned absence of a generalized water tables on the other.
The monsoon lasts five to six months a year in the core part of the highlands, and annual precipitation is very low during the dry season. The low yet permanent streamflows of the main rivers are then fed by the emptying of the eight thousand small wetlands (commonly settled on alterites), of the thin aquifers, and of the free water of the lowlands alongside rivers and streams of the whole Fouta Djallon and Guinean Dorsale.

What Could Change
D'Aubreville (1947) [12] and then Maignien (1966) [13] were impressed by the great extension of "ferrallitic pans in Guinea". The former thought that the formation of these "bowe" (sing. "bowal"), or "bowalisation" (as the process), is still ongoing, yet most of them have proven to be fossils [8]. These are not evidence of any land degradation processes and results.
The most densely populated area is the summit of the FDss (around Labé and Mali town) and its historical core. This location is very rainy, but it allows very poor water storage. People encounter severe problems in drinking water supply and even greater challenges to obtain water to irrigate a second or third yearly crop and produce enough food for the entire population.
As observed in the Sahel, the decrease in discharge during the prolonged rainfall deficit (1968-1993) of the Sahel's so-called "Great Drought" is approximately twice that of the rainfall [18]. During the drought, the authors highlighted that the severe rainfall deficit provoked for the first time in the Senegal River an exponential increase of the recession (or depletion) coefficient, which was stable for 70 years (a = 0.0186 d −1 ) and which strongly increased to reach 0.038 in 1984 [19]. This long "Great Drought" (Figure 3) led to a lasting degradation of the hydrological system. In Guinean areas, the decrease in discharges lasted 10 to 15 years. The increase in recession coefficients also lasted for years beyond the end of the drought. In 1997, Bricquet et al. [20] showed that despite the observed rainfall recovery, discharges did not begin to increase again at that time. This was due to the necessary refilling of the various natural storages such as surface and sub-surface reservoirs, and streamflows recovered over more than a decade [21]. Figure 4 [1] provides the spatial distribution of rainfall.    However, Gomis observed in 2000 [22] that the rainfall recovery led to a decrease in the recession coefficient of the Gambia River. In the Bafing basin (Upper Senegal basin), this slow and progressive discharge recovery was also characterized by a decrease in the recession coefficient, from 0.025 d −1 in 1972-1994 to 0.022 d −1 during the 1995-2005 period at the Bafing Makana stream gauge station [23].
In some basins, dams have been built or are being planned, and they will strongly modify the annual regime and allow hydroelectric power generation, such as in the Konkouré River basin [24].
This paper aims to study the recent evolution of the FD as a natural water tower through the analysis of runoff and depletion coefficients, along other factors such as rainfall and land cover dynamics.

Material and Methods
This study of the permanency of the Guinean water tower is based on the analysis of the hydrological behavior evolution of 16 basins and sub-basins across these highlands. Basins and their downstream gauging stations are indicated in Figure 5. However, Gomis observed in 2000 [22] that the rainfall recovery led to a decrease in the recession coefficient of the Gambia River. In the Bafing basin (Upper Senegal basin), this slow and progressive discharge recovery was also characterized by a decrease in the recession coefficient, from 0.025 d −1 in 1972-1994 to 0.022 d −1 during the 1995-2005 period at the Bafing Makana stream gauge station [23].
In some basins, dams have been built or are being planned, and they will strongly modify the annual regime and allow hydroelectric power generation, such as in the Konkouré River basin [24].
This paper aims to study the recent evolution of the FD as a natural water tower through the analysis of runoff and depletion coefficients, along other factors such as rainfall and land cover dynamics.

Material and Methods
This study of the permanency of the Guinean water tower is based on the analysis of the hydrological behavior evolution of 16 basins and sub-basins across these highlands. Basins and their downstream gauging stations are indicated in Figure 5. Excel file in supplementary materials gives the list of the 16 considered basins and the availability of their hydrological data, where vertical black lines separate the three climatic periods: 1947-1969, 1970-1993 and 1994-2019. Depletion coefficient (or recession coefficient). The recession coefficient is used after the Maillet law (1905) and as cited by Tallaksen (1995) [25] as follows: Qt = discharge at t time; Q0 = initial discharge at the beginning of the receding waters; α = recession index (= DP). The latter is the inverse of the elapsed time in d −1 and is expressed as The depletion coefficient (DC) used here are the average value of the DC for three durations: one month, two months, and four months (three on the smaller basins). The determined value of the DC is the average value of the DC for each duration around the maximal value of the DC for this duration. The considered period is the one when, on average for the period, the DC is the highest; thus, it is the same for all the years of the period, in the ensuing "Results" part. In a "Supplementary Material" part, we considered, for each station, and for each year, and for three durations (1, 2, and 4 months), the period when the DC reaches its highest value.

Discharges and Runoff Coefficients
Daily discharge values are provided by the basin agencies (ABN/NBA Niger Basin Authority, for the Niger basin; OMVS-Organisation pour la Mise en Valeur du bassin du fleuve Sénégal-for the Senegal basin; OMVG-Organisation pour la Mise en Valeur du bassin du fleuve Gambie-for the Gambia and Corubal basins) and the National Hydraulic Service of Guinea for the Konkouré basin.
Rainfall data of 23 stations are provided by SIEREM (Système d'Informations Environnementales sur les Ressources en Eau et leur Modélisation) [26,27], AMMA Catch (AMMA = African Monsoon Muldisciplinary Analysis) (http://www.amma-catch.org/), ANACIM (Agence Depletion coefficient (or recession coefficient). The recession coefficient is used after the Maillet law (1905) and as cited by Tallaksen (1995) [25] as follows: Q t = discharge at t time; Q 0 = initial discharge at the beginning of the receding waters; α = recession index (= DP). The latter is the inverse of the elapsed time in d −1 and is expressed as The depletion coefficient (DC) used here are the average value of the DC for three durations: one month, two months, and four months (three on the smaller basins). The determined value of the DC is the average value of the DC for each duration around the maximal value of the DC for this duration. The considered period is the one when, on average for the period, the DC is the highest; thus, it is the same for all the years of the period, in the ensuing "Results" part. In a "Supplementary Material" part, we considered, for each station, and for each year, and for three durations (1, 2, and 4 months), the period when the DC reaches its highest value.

Discharges and Runoff Coefficients
Daily discharge values are provided by the basin agencies (ABN/NBA Niger Basin Authority, for the Niger basin; OMVS-Organisation pour la Mise en Valeur du bassin du fleuve Sénégal-for the Senegal basin; OMVG-Organisation pour la Mise en Valeur du bassin du fleuve Gambie-for the Gambia and Corubal basins) and the National Hydraulic Service of Guinea for the Konkouré basin.
Rainfall data of 23 stations are provided by SIEREM (Système d'Informations Environnementales sur les Ressources en Eau et leur Modélisation) [26,27], AMMA Catch (AMMA = African Monsoon Muldisciplinary Analysis) (http://www.amma-catch.org/), ANACIM (Agence Nationale de l'Aviation CIvile et de la Météorologie) (Senegalese data), and from the CHIRPS (Climate Hazards Group InfraRed Precipitation with Station data) data system [28] for recent years (since 1980, 1994, or 2000 according to the basins). The CHIRPS products used here have a spatial resolution of 0.25 • and a daily temporal resolution. The covered area is from 50 • N to 50 • S. They combine satellite data at 0.25 • to soil rainfall data [28]. In order to improve the data accuracy, the CHIRPS data were unskewed again using observed rainfall data for a few stations (see Figure 4)

At the Regional Scale
At the regional scale, the tree cover variations between 2000 and 2019 over the 16 basins were extracted from the Global Forest Change (GFC2019) datasets produced by the University of Maryland, based on Landsat 4, 5, 7, and 8 imagery. These provide 30 m spatial resolution data of tree gains and tree losses compared to land use cover of the year 2000. Cloud-free observations for each pixel, based on quality assessment layers, are used to create composite imagery in bands 3, 4, 5, and 7 during the growing season. Image interpretation relies on a decision tree supervised learning approach and training data, including very high resolution Quickbird imagery as described in Hansen et al. (2013) [29]. Google Earth Engine© (Google, Silicon Valley, CA, USA) geoprocessing capabilities [30] were used here to reduce the time to extract and combine pixel data for the 16 basins.

At the Basin Scale
At the basin scale, a special analysis was conducted in the upper Bafing basin. At the Upper Bafing basin scale, changes in vegetation were searched using a time series analysis (TSA) of the normalized difference vegetation index (NDVI) based on a MODIS (Moderate-Resolution Imaging Spectroradiometer; EOS Earth Observing System, NASA, Washington, D.C., USA) (MODQ13 dataset (250 m resolution; 16 days temporal resolution; from February 2000 to February 2020). Mann-Kendall's tau test is an efficient trend detection method [31] and has been used for this data type for other regions of western Africa [32][33][34]. Mann-Kendall's tau test can be submitted to null hypothesis (P value) testing as the data come from a population with independent realizations and are identically distributed. A significance at 0.001 has been retained for this paper meaning that a sample as extreme as that observed would occur 1 time in 1000 if the null hypothesis is true.

At the Local Scale
At the local scale, some small wetlands of river sources are compared in the central area of the FD. Here, we compare the current status (via Google Earth© images; Google, Silicon Valley, CA, USA) with aerial imagery acquired by the IGN (Institut National de l'Information Géographique et Forestière, Saint Mandé, France) 1953 mission over West Africa. the first one is the 'hyper humid' 1947-1969 period, when rainfall was significantly higher than the long-term average across the whole of West Africa. the second period is the well-known "Great Drought" of West Africa (1970-1993). The rainfall decreased significantly in 1968, but we arbitrarily made the second period begin in 1970 because   Data covering the 1947-2019 timeline was divided into three periods, based on previous segmentation by several authors ( [35][36][37][38]):

Changes of Depletion Coefficients
the first one is the 'hyper humid' 1947-1969 period, when rainfall was significantly higher than the long-term average across the whole of West Africa. -the second period is the well-known "Great Drought" of West Africa . The rainfall decreased significantly in 1968, but we arbitrarily made the second period begin in 1970 because a significant number of stream gauge stations have data starting from 1970. The drought is considered to last until 1993. -the last period spans the years 1994-2019. The year 1994 is when annual rainfall again reached again its long-term mean value. During this period, rainfall is approximately of the long-term trend.
Rainfall predictions in West Africa indicate a regional increase in rainfall during the ongoing 21st century due to the activation of monsoons as a result of the increase in the difference between the continental temperature and SST (sea surface T • ) [39]. However, locally, in the western part of West Africa, rainfall is expected to decrease significantly after 2030 or 2035 (Senegambia, Mauritania, and western Mali), while further south, an increase in ATR (annual total rainfall) across most of the Guinean territory is expected [40,41].
The dynamics of the depletion coefficients (see  (Table 2C) and higher than 20% in five others (Table 2D). -some basins (Niandan, Milo, Gambia, Tominé) show a marked increase in the DC at the end of the last period: in all of these basins, depletion was higher during the 2005-2015 period than during the Great Drought period.
-Extreme values can be observed in the last years (e.g., in the Milo basin), and in 1983-1984 in the Faleme. As described by Olivry (1986), this happens when the higher value of depletion coefficient occurs when the river dries up completely. Consequently, the DC becomes in some cases the drying up coefficient.
General observations in some basins, the depletion coefficient can be different depending on the method selected or on the context. -except the increase in the DC during the drought and decrease after the drought in the large basins, there are no general trends to detect, reflecting great variability across the basins. -overall, contrary to previous observations by Descroix [5] and Descroix et al. [42], which were based on the behavior of the Upper Niger basin tributaries, there is not a general trend of increasing DC, which could demonstrate a general reduction in the soil water-holding capacity. This trend only appears in the Konkouré basin and in the Milo and Niandan basins.
In Table 2, the last column provides the evolution of forest cover by basin based on the Global Forest Change (GFC2019) methodology.
The day when the maximal depletion occurs seems to follow a general behavior, occurring earlier when the depletion coefficient increases (mostly during the drought), later during the first period (yet with numerous gaps in the data), and overall during the last period. This follows an intuitive pattern, with a depletion stage of the river starting later when discharges are abundant and earlier during the dry years. The exception is the case of the Konkouré basin, where the upper basin (Pont de Linsan station) is the one wherein the discharge has decreased the most. The same evolution is shown at the downstream station "Pont de la route de Télimélé", but this station has been influenced since 1999 by the filling and management of the Garafiri dam. This is why its DC decreased strongly after 1999. Table 3 gathers the basins' hydrological characteristics and their changes over the three determined periods: discharge (Q), specific discharge (SD), runoff coefficient (RC), as well as the maximal discharge of the year and its day of occurrence. The general dynamic follows the rainfall trends with a decrease in streamflow between the first period  and the second , followed by an increase in discharge between the second period and the third (1994-2019). Yet two important exceptions have to be highlighted:

Evolution of Runoff Coefficients
- the Konkouré basin, where runoff did not rise again after the dry period, but in this basin, data is not available after 2002, preventing analysis of the long term and obscuring a possible recovery of streamflows. However, since rainfall increased significantly between period 2 and 3, a decrease in the runoff coefficient is noticed. -the Upper Niger basin where runoff increased only in the Tinkisso basin (except for its upper part) over the last period.
This could be explained by a longer time being necessary to refill the natural water reservoir. The maximum daily discharge observed each year shows approximately the same trends (Table 3, last column): the day when this maximum is observed is in most cases later when the discharge increases and earlier when streamflows are low.

Any Explanation from the LULCC?
The analysis of the relationship between, on the one hand, runoff and depletion coefficients and, on the other hand, land use/land cover changes (LULCC) during the period of 2000-2019 is made with the data provided by the 30 m Global Forest Change datasets [29] (Figures 7 and 8). Considering no local field validation has been done and the inconsistencies in some results, the forest evolution shown here must be considered as relative and not as absolute.
Water 2020, 12, x FOR PEER REVIEW 19 of 32

Any Explanation from the LULCC?
The analysis of the relationship between, on the one hand, runoff and depletion coefficients and, on the other hand, land use/land cover changes (LULCC) during the period of 2000-2019 is made with the data provided by the 30 m Global Forest Change datasets [29] (Figures 7 and 8). Considering no local field validation has been done and the inconsistencies in some results, the forest evolution shown here must be considered as relative and not as absolute.
One of the most important trends is the higher deforestation rate measured in the Upper Niger basin (Kouroussa stream gauge) and the Konkouré basin (mostly in the basins of Pont de Linsan and Pont de la Route de Télimélé stations, secondarily Kaba and Nianso stations). These basins are the ones where: -RC paradoxically decreases or remains low during the last period, which is characterized by rainfall recovery, possibly because natural reservoirs continue to be filled. Can this hypothesis however still be valid after nearly 25 years of rainfall recovery? -DC increased strongly (downstream Konkouré basin), moderately (Kaba and Nianso stations), or inversely, and remained unchanged (Upper Niger at Kouroussa station). in Table 3, the basins where DC increased more than 20% are clearly not the ones where the deforestation rate was the highest (except the small Pt de Linsan-Upper Konkouré River basin).

LULCC at the Regional Scale
GFC2019 rasters reveal that the whole south area of the FD has a high value of vegetation cover; only the Upper Bafing basin is extensively deforested, corresponding to the most densely populated area: this is the extension of the 18th and 19th century Peul/Fulani confederation of small theocratic states which dominated West Africa. Deforestation in recent decades on the western side of the FD is not necessarily due to the extension of this high density of the upper plateaus area. It is more likely due to wood-charcoal extraction/forest exploitation. A very low number of blue pixels appear in Figures 7 and 8, reflecting the low vegetation gain detected by GFC2019, which is partly due to methodological reasons of vegetation density and vegetation height. The Bebele basin (Tene River at Bebele stream gauge station), which is presented in Figures 7 and 8, is not considered in this analysis, due to its very limited hydrological data set.
The northern areas on Figure 7 correspond to the driest Sudanian zone, covered with bushes, savannah, and low forest cover. Forest vegetation density falls below 15% here, and this limited vegetation is excluded by the vegetation cover threshold employed here. One of the most important trends is the higher deforestation rate measured in the Upper Niger basin (Kouroussa stream gauge) and the Konkouré basin (mostly in the basins of Pont de Linsan and Pont de la Route de Télimélé stations, secondarily Kaba and Nianso stations). These basins are the ones where: -RC paradoxically decreases or remains low during the last period, which is characterized by rainfall recovery, possibly because natural reservoirs continue to be filled. Can this hypothesis however still be valid after nearly 25 years of rainfall recovery? -DC increased strongly (downstream Konkouré basin), moderately (Kaba and Nianso stations), or inversely, and remained unchanged (Upper Niger at Kouroussa station). -in Table 3, the basins where DC increased more than 20% are clearly not the ones where the deforestation rate was the highest (except the small Pt de Linsan-Upper Konkouré River basin).

Analysis of LULCC
6.1. LULCC at the Regional Scale GFC2019 rasters reveal that the whole south area of the FD has a high value of vegetation cover; only the Upper Bafing basin is extensively deforested, corresponding to the most densely populated area: this is the extension of the 18th and 19th century Peul/Fulani confederation of small theocratic states which dominated West Africa. Deforestation in recent decades on the western side of the FD is not necessarily due to the extension of this high density of the upper plateaus area. It is more likely due to wood-charcoal extraction/forest exploitation. A very low number of blue pixels appear in Figures 7  and 8, reflecting the low vegetation gain detected by GFC2019, which is partly due to methodological reasons of vegetation density and vegetation height. The Bebele basin (Tene River at Bebele stream gauge station), which is presented in Figures 7 and 8, is not considered in this analysis, due to its very limited hydrological data set.
The northern areas on Figure 7 correspond to the driest Sudanian zone, covered with bushes, savannah, and low forest cover. Forest vegetation density falls below 15% here, and this limited vegetation is excluded by the vegetation cover threshold employed here. Figure 9 is a map of the Mann-Kendall tau after applying, for each pixel, the P value test at 0.001 and classifying in three categories: significant positive trend, stable or not significant trend, and significant negative trend.

LULCC at the Upper Bafing Basin Scale
Water 2020, 12, x FOR PEER REVIEW 21 of 32 Figure 9 is a map of the Mann-Kendall tau after applying, for each pixel, the P value test at 0.001 and classifying in three categories: significant positive trend, stable or not significant trend, and significant negative trend. Eighty-six percent of the zone covered by the map (Figure 9) has shown stable values of NDVI from 2000 to 2020 or non-significant trends according to the P value test at 0.001, and 11% of the zone has shown a significant positive trend (greening; biomass increase) for the same period. Only 3% has shown a significant negative trend (browning; biomass decrease).

LULCC at the Upper Bafing Basin Scale
Unlike the whole Fouta Djalon region, the Upper Bafing area is not concerned by important changes in its vegetation over the last 20 years, and the main change is a positive one for vegetation cover, even if it is hard to say without field validation if it corresponds to higher vegetation cover, higher biomass (either for spontaneous vegetation or for crops), or simple betterment of vegetative activities (i.e., phenology with longer green season).
The basins are not very different regarding this analysis. Stability ranges from 84.1% in the Bafing basin to 90.5% in the Niger basin ( Figure 10). Negative trends range from 0.4% in the Gambia basin to 4.4% in the Tinkisso basin. Positive trends range from 6.3% to 14.6% in the Niger and Bafing basins, respectively. Eighty-six percent of the zone covered by the map (Figure 9) has shown stable values of NDVI from 2000 to 2020 or non-significant trends according to the P value test at 0.001, and 11% of the zone has shown a significant positive trend (greening; biomass increase) for the same period. Only 3% has shown a significant negative trend (browning; biomass decrease).
Unlike the whole Fouta Djalon region, the Upper Bafing area is not concerned by important changes in its vegetation over the last 20 years, and the main change is a positive one for vegetation cover, even if it is hard to say without field validation if it corresponds to higher vegetation cover, higher biomass (either for spontaneous vegetation or for crops), or simple betterment of vegetative activities (i.e., phenology with longer green season).
The basins are not very different regarding this analysis. Stability ranges from 84.1% in the Bafing basin to 90.5% in the Niger basin ( Figure 10). Negative trends range from 0.4% in the Gambia basin to 4.4% in the Tinkisso basin. Positive trends range from 6.3% to 14.6% in the Niger and Bafing basins, respectively.

LULCC at the Sources Scale
The comparison between the 1953 aerial picture (IGN-see Section 4.3 aerial mission in Guinea) and the Google Earth satellite scene (2019) does not demonstrate any apparent deforestations ( Figure 11). The quality of both products being very different, it is difficult to do a classical numerical comparison, but the current situation (2019, on the right) seems a priori to be at least as forested as the one in 1953 (Figure 11). This area is named here the "three sources area": effectively, the two creeks (brooks) flowing through the north are the main sources of the Kouba (Koumba) River and the main upper reach of the Koliba River, which becomes the Rio Corubal in Guinea Bissau. The one flowing eastward is the Dima, the upper Gambia River. 'Ore Dima' means "the sources of the Dima River" in Pulaar language (Ore is the head). The other two rivers flowing southwards and the one flowing through the west are the main sources of the Kakrima River, the main tributary of the Konkouré River. This area is mostly composed by a doleritic bedrock.
Changes in vegetation cover evolution over time appear different at these three spatial scales but are difficult to compare considering the different analytical methodologies which are here applied across different areas, periods, and resolutions. However, clearly, some areas are being deforested, though results do not conclude that there is an overall deforestation, contrary to what is commonly described in the literature. Our three studies notably disagree in identifying negative trends in vegetation cover.

LULCC at the Sources Scale
The comparison between the 1953 aerial picture (IGN-see Section 4.3 aerial mission in Guinea) and the Google Earth satellite scene (2019) does not demonstrate any apparent deforestations ( Figure 11). The quality of both products being very different, it is difficult to do a classical numerical comparison, but the current situation (2019, on the right) seems a priori to be at least as forested as the one in 1953 ( Figure 11). This area is named here the "three sources area": effectively, the two creeks (brooks) flowing through the north are the main sources of the Kouba (Koumba) River and the main upper reach of the Koliba River, which becomes the Rio Corubal in Guinea Bissau. The one flowing eastward is the Dima, the upper Gambia River. 'Ore Dima' means "the sources of the Dima River" in Pulaar language (Ore is the head). The other two rivers flowing southwards and the one flowing through the west are the main sources of the Kakrima River, the main tributary of the Konkouré River. This area is mostly composed by a doleritic bedrock.

LULCC at the Sources Scale
The comparison between the 1953 aerial picture (IGN-see Section 4.3 aerial mission in Guinea) and the Google Earth satellite scene (2019) does not demonstrate any apparent deforestations ( Figure 11). The quality of both products being very different, it is difficult to do a classical numerical comparison, but the current situation (2019, on the right) seems a priori to be at least as forested as the one in 1953 ( Figure 11). This area is named here the "three sources area": effectively, the two creeks (brooks) flowing through the north are the main sources of the Kouba (Koumba) River and the main upper reach of the Koliba River, which becomes the Rio Corubal in Guinea Bissau. The one flowing eastward is the Dima, the upper Gambia River. 'Ore Dima' means "the sources of the Dima River" in Pulaar language (Ore is the head). The other two rivers flowing southwards and the one flowing through the west are the main sources of the Kakrima River, the main tributary of the Konkouré River. This area is mostly composed by a doleritic bedrock.
Changes in vegetation cover evolution over time appear different at these three spatial scales but are difficult to compare considering the different analytical methodologies which are here applied across different areas, periods, and resolutions. However, clearly, some areas are being deforested, though results do not conclude that there is an overall deforestation, contrary to what is commonly described in the literature. Our three studies notably disagree in identifying negative trends in vegetation cover.  Changes in vegetation cover evolution over time appear different at these three spatial scales but are difficult to compare considering the different analytical methodologies which are here applied across different areas, periods, and resolutions. However, clearly, some areas are being deforested, though results do not conclude that there is an overall deforestation, contrary to what is commonly described in the literature. Our three studies notably disagree in identifying negative trends in vegetation cover.

Is There Any Synthesis Provided by Data Analysis?
A PCA (Principal Component Analysis) is realized in order to provide clearer information about the variables and observations (basins) sampled in the FD. Table 4 details the variables considered in this analysis. The observations (sampling) are the 16 basins studied. No geological data are accounted for in this study, because it appears that there are no significant differences in the hydrological functioning of the different outcrops in the FD.
Water 2020, 12, x FOR PEER REVIEW 23 of 32

Is There Any Synthesis Provided by Data Analysis?
A PCA (Principal Component Analysis) is realized in order to provide clearer information about the variables and observations (basins) sampled in the FD. Table 4 details the variables considered in this analysis. The observations (sampling) are the 16 basins studied. No geological data are accounted for in this study, because it appears that there are no significant differences in the hydrological functioning of the different outcrops in the FD.
In both Figures 12 and 13: axis 1 (component 1) explains 46.4% of the total variance; it is mostly provided by QS, RAD, RUD, and RC (The definition of all abbreviations is given in Table 4).   Figure 12a shows some redundancy of, on the one hand, hydrological variables (RAD, RUD, QS, RC) and, on the other, the depletion coefficient, separated from their evolution in two different clusters. They are, by contrast, gathered according to axis 3 and axis 4 (Figure 12b). The LULCC is considered here through the regional mapping. Although its results are different (and locally contradictory) with the ones provided by the other LULCC analyses, we considered these to be suitable as relative values rather than as absolute values. The LULCC is also then determined with the same data and methods for the 16 basins. Deforestation as shown in Figure 7 is correlated with  Figure 12a shows some redundancy of, on the one hand, hydrological variables (RAD, RUD, QS, RC) and, on the other, the depletion coefficient, separated from their evolution in two different clusters. They are, by contrast, gathered according to axis 3 and axis 4 (Figure 12b). The LULCC is considered here through the regional mapping. Although its results are different (and locally contradictory) with the ones provided by the other LULCC analyses, we considered these to be suitable as relative values rather than as absolute values. The LULCC is also then determined with the same data and methods for the 16 basins. Deforestation as shown in Figure 7 is correlated with the increase in runoff between period 1 and 2, but negatively correlated with the increase in runoff between period 2 and 3. Relative deforestation is clearly correlated with the increase in the DC (depletion coefficient) in four months during the second period, the one when rainfall recovery should have allowed a significant regreening. However, it is not so clearly correlated with the DC in one month ( Figure 12a); yet following components 3 and 4 (Figure 12b), they seem more clearly correlated. Figure 13 shows firstly a great dispersion of samples (basins) in all the statistical spaces, as evidence of their great geographical and hydrological diversity. Figure 13a seems logically dominated by runoff, and then, the basins of the western slope of the FDss and the ones coming off the Guinean Dorsale are clearly opposed to the ones of the northern part, where streamflow is weaker. No general trend appears in Figure 13b.

Discussion
The following points must be highlighted:

The DC
A strong interannual increase in the DC (see also [5,42]) is observed for the two main tributaries of the Upper Niger River from 1995. However, this is not observed in the Upper Niger (Kouroussa) nor in the Tinkisso River ( Figure 6e). Otherwise, the DC, according to the three durations tested here, can be similar (e.g., in the Niandan River) or very different (e.g., in the Milo River) A great deal was devoted to the DC during the 1990s and the 2000s by Mahé et al. [43,44] (1997 and 2001) and Olivry [18,19,21] (1987, 1993, 2002), as well as Bricquet et al. [20] (1997), Bamba et al. [45] (1996), and Sangaré et al. [46] (Figure 14). The latter showed that the DC increased from the beginning of the 1970s in all the basins, even if more slowly for the Sankarani River (one of the other main tributaries of the Upper Niger River, not investigated here). The exception was then the Milo River, the one whose DC has further increased in recent years. The DC of the Milo River only started its increase at the beginning of the 1980s [39]. A possible role of the basins' water-holding capacity remains to be investigated because the geological context is probably more complex than the schematic regionalization in Figure 2. The Upper Niger River basin is mostly a granitic basement one, however, Brunet Moret et al. (1986) [47] showed that there are sometimes sandstone outcrops in the lower basins of the Milo River. This could have served as a better water reservoir than granite and thus reduced the impact of the first drought spell in the 1970s. Another explanation can be provided by the rainfall's variation coefficient, which is lower in the Upper Milo basin than in most places in West Africa, as was noticed by Mahé et al. [44], most probably due to the specific relief orientation in the area. The Upper Milo basin runs through elevated areas both in the South and in the East, which seem to play a role in the reduction of the interannual variability, as is also observed in Central Ghana [44]. As can be seen ( Figure 14) the DC of the Sankarani river also showed a reduced variability compared to that of other Upper Niger River tributaries, during the Great Drought. In the Senegal River basin, Bodian et al. (2020) [48] observed that the runoff recovery after the Great Drought lasted for five or ten years and is more significant for yearly and maximum discharge and less for low waters, which remain very weak. The same authors concluded that the high value of the DC for the Faleme River at Gourbassi (0.62) at the beginning of the 1980s (Figure 6d) can be attributed to a quasi-total drying up of the streamflow. A similar observation can be made for the Milo River during the 2010s (Figure 6e), but the reason remains to be explained since neither rainfall evolution nor the LULCC are satisfactory explaining factors for this behavior. Water 2020, 12, x FOR PEER REVIEW 26 of 32 Figure 14. Evolution of DC in some Sudano-Guinean tributaries of the Niger River before and during the drought [41].

LULCC (Land Use/Land Cover Change)
As it has been shown above, vegetation cover in West Africa does not decline systematically though the population keeps on increasing at a rapid pace [5]. The Fouta Djallon is one of the best examples of this trend. A negative and Afro-pessimistic belief often links African deforestation with demographic growth. This stereotype anticipated catastrophic scenarios based on extensive peasant practices, which, against the background of strong demographic growth, should threaten the integrity of the Fouta Djallon ecosystems and therefore endanger the whole 'West African water tower'. In fact, local research has opened the debate on the reality and validity of such claims [49]. Brandt et al. (2017) [50] showed localized reductions in woodlands across West Africa but that farmland in Sahelian areas also promotes woody cover, contradicting simplistic ideas of a high negative correlation between population density and woody cover.
Finally, it is worth noting that the 'tapades' traditional intensive gardening and tree cropping practices of the Fouta Djallon, which have played the role of man-made sponges though they have mostly been managed by women, are largely being abandoned. These practices should be rehabilitated in order to promote a sustainable agricultural intensification as the most useful way to infiltrate and store water in a catchment [5].

Some Hydrological Consequences
Except during the drought, the increase in the DC has only been observed in a few basins in recent years (the Konkouré and southern tributaries of the Upper Niger). The partial recovery in the RC (runoff coefficient), as well as the observed increase-at the regional scale-of the occurrence of extreme rainfall events [36,42,56,57], has consequences for flooding frequency and intensity.

LULCC (Land Use/Land Cover Change)
As it has been shown above, vegetation cover in West Africa does not decline systematically though the population keeps on increasing at a rapid pace [5]. The Fouta Djallon is one of the best examples of this trend. A negative and Afro-pessimistic belief often links African deforestation with demographic growth. This stereotype anticipated catastrophic scenarios based on extensive peasant practices, which, against the background of strong demographic growth, should threaten the integrity of the Fouta Djallon ecosystems and therefore endanger the whole 'West African water tower'. In fact, local research has opened the debate on the reality and validity of such claims [49]. Brandt et al. (2017) [50] showed localized reductions in woodlands across West Africa but that farmland in Sahelian areas also promotes woody cover, contradicting simplistic ideas of a high negative correlation between population density and woody cover.
Finally, it is worth noting that the 'tapades' traditional intensive gardening and tree cropping practices of the Fouta Djallon, which have played the role of man-made sponges though they have mostly been managed by women, are largely being abandoned. These practices should be rehabilitated in order to promote a sustainable agricultural intensification as the most useful way to infiltrate and store water in a catchment [5].

Some Hydrological Consequences
Except during the drought, the increase in the DC has only been observed in a few basins in recent years (the Konkouré and southern tributaries of the Upper Niger). The partial recovery in the RC (runoff coefficient), as well as the observed increase-at the regional scale-of the occurrence of extreme rainfall events [36,42,56,57], has consequences for flooding frequency and intensity. This increases streamflow irregularities, accentuating low waters as well as floods. This could have impacts on hydraulic and civil engineering equipment and infrastructures. As illustrated in Figures 15 and 16, the bridges are threatened by some of these extreme flooding events.
This increases streamflow irregularities, accentuating low waters as well as floods. This could have impacts on hydraulic and civil engineering equipment and infrastructures. As illustrated in Figures 15 and 16, the bridges are threatened by some of these extreme flooding events.
However, as this paper attempted to demonstrate, it is not necessarily the sign of a decrease in soil water-holding capacity attributable to LULCC. Depletion coefficients of the main rivers providing freshwater resources to the Sahel are decreasing due to the rainfall recovery. West African agro-ecosystems seem to be developing resilience, which could invert the observed trend and constitute rural resiliency models and sustainable development practices.    This increases streamflow irregularities, accentuating low waters as well as floods. This could have impacts on hydraulic and civil engineering equipment and infrastructures. As illustrated in Figures 15 and 16, the bridges are threatened by some of these extreme flooding events.
However, as this paper attempted to demonstrate, it is not necessarily the sign of a decrease in soil water-holding capacity attributable to LULCC. Depletion coefficients of the main rivers providing freshwater resources to the Sahel are decreasing due to the rainfall recovery. West African agro-ecosystems seem to be developing resilience, which could invert the observed trend and constitute rural resiliency models and sustainable development practices.   However, as this paper attempted to demonstrate, it is not necessarily the sign of a decrease in soil water-holding capacity attributable to LULCC. Depletion coefficients of the main rivers providing freshwater resources to the Sahel are decreasing due to the rainfall recovery. West African agro-ecosystems seem to be developing resilience, which could invert the observed trend and constitute rural resiliency models and sustainable development practices.

Conclusions
The Fouta Djallon seems to have suffered from a reduction in water-holding capacity during the drought (1970s and 1980s). However, after this dry period, in some basins and mostly in those flowing towards Sahelian areas, the expected decrease in the DC due to rainfall recovery is observed (the Niger at Kouroussa, the Falémé at Gourbassi, probably the Gambia and Corubal basin rivers). The Bafing at Daka Saidou did not show any increases in DC during the drought.
In the Upper Niger basin (the Milo and Niandan rivers), the Konkouré, and partially the Koliba/Tominé basins, a significant increase in the DC was noticed during the end of the recent 1994-2019 period. However, in the Upper Niger basin, the DC increased more in the basins where relative deforestation was lower (Milo and Niandan) and decreased slightly in the basin where relative deforestation was stronger (Upper Niger at Kouroussa). Following the PCA analysis, only the increase in the DC in four months between period 2 and period 3 is clearly correlated with relative deforestation. Besides the well-known and well-studied increase in rainfall intensity (increase in extreme events occurrence and increase in rainfall hourly mean intensity), the LULCC could be a discriminant factor of any increase in runoff and depletion coefficients, particularly in the Konkouré basin where deforestation remains an important dynamic. The PCA analysis only discriminates the southwestern basins (the Konkouré and Koliba/Tominé basins).
Contrary to popular belief, there is evidence of successful rural intensification practices in West Africa, which are generally characterized by an increasing vegetation cover and, locally, increasing tree cover. This kind of landscape should logically improve the soil and basins' water-holding capacity and eventually lead to a reduction in the DC (depletion coefficient). In the core of the FD (i.e., the FDss), the 'tapades' intensive system is disappearing, yet no deforestation is observed due to several combined patterns.
Since the sole increase in extreme rainfall events is unlikely to be responsible for the increase in the DC, LULCC seems to be a factor in the reduction in the basin water-holding capacity. However, it is statistically apparent only in the Konkouré and in the Koliba/Tominé basins, where annual rainfall is higher and where the runoff coefficients are significantly higher than in all the other basins, making their susceptibility to changes in runoff conditions more influential than in other basins.
The Fouta Djallon, as the West African water tower, does not appear as threatened. However, there is a stark difference between, on the one hand, the southwestern basins (the Konkouré and Upper Corubal), where LULCC seems to impact the basin water-holding capacity, and, on the other hand, basins such as the Gambia, Senegal, and Niger basins, which have been impacted more by the Great Drought but which are at least partially recovering their natural water storage capacity. This is an optimistic observation for the future of the West African Sahel, whose basins depend the most on the Fouta Djallon for its supply of fresh water resources.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/12/11/2968/s1. Excel 1: 16 considered basins and the availability of their hydrological data. Figure S1a: Evolution of depletion coefficient, Konkouré River basin. Figure S1b: Evolution of rainfall and runoff coefficient, Konkouré River basin. Figure S2a: Evolution of depletion coefficient, Upper Konkouré River basin. Figure S2b: Evolution of rainfall and runoff coefficient, Konkouré River basin. Figure S3a: Evolution of depletion coefficient, Kokoulo River basin. Figure S3b: Evolution of rainfall and runoff coefficient, Kokoulo River basin. Figure S4a: Evolution of depletion coefficient, Kakrima River basin. Figure S4b: Evolution of rainfall and runoff coefficient, Kakrima River basin. Figure S5a: Evolution of depletion coefficient, Tominé River basin. Figure S5b: Evolution of rainfall and runoff coefficient, Tominé River basin. Figure S6a: Evolution of depletion coefficient, Koliba/Koumba River basin. Figure S6b: Evolution of rainfall and runoff coefficient, Koliba/KoumbaRiver basin. Figure S7a: Evolution of depletion coefficient, Koulountou River basin. Figure S7b: Evolution of rainfall and runoff coefficient, Koulountou River basin. Figure S8a: Evolution of depletion coefficient, Gambia River basin. Figure S8b: Evolution of rainfall and runoff coefficient, Gambia River basin. Figure S9a: Evolution of depletion coefficient, Faleme River basin. Figure S9b: Evolution of rainfall and runoff coefficient, Faleme River basin. Figure S10a: Evolution of depletion coefficient, Upper Bafing River basin. Figure S10b: Evolution of rainfall and runoff coefficient, Upper Bafing River basin. Figure S11a: Evolution of depletion coefficient, Bafing River basin Figure S11b: Evolution of rainfall and runoff coefficient, Bafing River basin. Figure S12a: Evolution of depletion coefficient, Upper Tinkisso River basin. Figure S12b: Evolution of rainfall and runoff coefficient, Upper Tinkisso River basin. Figure S13a: Evolution of depletion coefficient, Tinkisso River basin. Figure S13b: Evolution of rainfall and runoff coefficient, Tinkisso River basin. Figure S14a: Evolution of depletion coefficient, Upper Niger River basin. Figure S14b: Evolution of rainfall and runoff coefficient, Upper Niger River basin. Figure S15a: Evolution of depletion coefficient, Niandan River basin. Figure S15b: Evolution of rainfall and runoff coefficient, Niandan River basin. Figure S16a: Evolution of depletion coefficient, Milo River basin. Figure S16b: Evolution of rainfall and runoff coefficient, Milo River basin.