2.1. Study framework
Covering areas of 61810 km2 in Garoua and 67012 km2 in Bongor respectively (Table 1), the Benoue and Logone basins are located in the northern tropical zone, between shared by Cameroon, Chad and the Central African Republic (roughly between 6 °5'–11°N and 13°–17°E). The largest parts of the Benoue and Logone basins are located in Cameroon and Chad respectively (Fig. 1). These basins are subject to a tropical Sudanese climate recording an average of precipitation of above 1000 mm/year (Table 2). This climate is characterized by a dry and hot season extending on average from November to April, and a wet and less hot season from May to October during which most of the precipitation occurs. The relief of these basins are rugged. There are mountain ranges including the Adamaoua and Poli with altitudes varying between 1000 m and 1500 m, but also plateaus and several areas of undulations of altitudes between 200 m and 800 m. These reliefs are essentially underlain by granites and crystalline schists, in vertical beds and generally oriented South-West, North-East (Rodier 1953). Savannas and gallery forests are the main forms of vegetation in these basins (Rodier 1953).
Table 1
Gauging stations and corresponding characteristics in the studied basins
Station | Watercourse | Geographical coordinates | Drainage area (km2) | Data series lenght |
| | Latitude N | Longitude E | | Rainfall | Discharge |
Garoua | Benoue | 9.29 | 13.4 | 61810 | 1950-51 to 2015-16 | 1950-51 to 2015-16 |
Bongor | Logone | 10.26 | 15.41 | 67012 | 1950-51 to 2017-18 | 1952-53-2017-18 |
Table 2
Statistics of precipitation and average flows (annual and seasonal) of the basins over their respective study periods, but also on both sides of the identified ruptures
Periods | Interannual | Year of rupture | Mean | CV (%) |
| Mean | CV (%) | | Before rupture | After rupture | Before rupture | After rupture |
Precipitation |
Benoue |
Annual | 1050 | 0.1 | 1970-71 | 1114 | 1020 | 0.06 | 0.1 |
Wet season (May-October) | 989 | 0.09 | 1970-71 | 1040 | 966 | 0.06 | 0.1 |
Dry season (November-April) | 61 | 0.43 | 1970-71 | 74 | 54 | 0.33 | 0.47 |
Logone |
Annual | 1221 | 0.09 | 1969-70 | 1298 | 1185 | 0.06 | 0.08 |
Wet season (May-October) | 1117 | 0.09 | 1970-71 | 1175 | 1090 | 0.08 | 0.08 |
Dry season (November-April) | 104 | 0.31 | 1969-70 | 125 | 92 | 0.25 | 0.3 |
Discharges |
Benoue |
Annual | 320.4 | 0.39 | 1970-71 | 388 | 275 | 0.22 | 0.46 |
Wet season (May-October) | 563 | 0.4 | 1970-71 | 724 | 454 | 0.22 | 0.45 |
Dry season (November-April) | 77 | 0.63 | 1984-85 | 41 | 151 | 0.43 | 0.3 |
Logone |
Annual | 496 | 0.33 | 1970-71 | 616 | 441 | 0.1 | 0.38 |
Wet season (May-October) | 786 | 0.32 | 1970-71 | 932 | 720 | 0.08 | 0.37 |
Dry season (November-April) | 197 | 0.47 | 1970-71 | 269 | 163 | 0.23 | 0.57 |
CV : Coefficient of variation |
2.2. Data source
Several datasets (hydrometeorological and land cover) were used in this work. Available for the periods 1950–2015 and 1952–2017, the flow series for Benoue and Logone come from the SIEREM database (Environmental Information System on Water Resources and Modeling). That of Benoue has significant gaps between 1998 and 2009. The gaps contained in that of Logone are mainly located during the intervals 1980–1983, 1994–1995 and 2008–2012.
The rainfall data used in this study are from CRU (Climate Research Unit) of the University of East Anglia in the United Kingdom. These data have been available since 1901 via the site https://climexp.knmi.nl/selectfield_obs2.cgi?id=2833fad3fef1bedc6761d5cba64775f0/ in NetCDF format, on a monthly time step and at a spatial resolution of 0.25° x 0.25 °. Precipitation and temperature data from CRU have been used to validate CMIP models in Logone basin which is a sub-basin of the Lake Chad basin (Nkiaka et al. 2018).
The spatial data used for the study of land cover in the Benoue basin are Landsat 8 satellite images from January to March 2018 and Landsat TM from January to March 1987, corresponding to path/row 185/53, 185 /54, 185/55, 184/53, 184/54, 184/55, 183/53 and 183/54. All of these images are made available to the general public free of charge by the National Aeronautics and Space Administration (NASA), via the US Geological Survey website (https://earthexplorer.usgs.gov/), in GeoTIFF format. The downloaded images taken during the dry season were preferred to those of the rainy seasons since they are less affected by cloud disturbances.
The demographic data used in this work are those of the population censuses in Cameroon (1976, 1987, 2005 and 2018) and Chad (1993 and 2009). These data were collected from the Central Bureau of Censuses and Population Studies.
2.3. Data analysis
The analysis of precipitation, average river discharge and runoff coefficients (at annual and seasonal time steps) was carried out using Pettitt (Pettitt 1979) and Mann-Kendall (Mann 1945; Kendall 1975; Yue et al. 2002) tests, at the 95% significance level.
The Pettitt test seems to be the most suitable for the analysis of incomplete series such as ours because it can separate the series into two periods with an overall distinct behavior, which avoids the detection of false discontinuities which can sometimes be observed with other tests such as Hubert segmentation. Its principle consists of dividing the studied series (of N size) into two sub-samples of sizes m and n, respectively. We then calculated the sum of the ranks of the elements of each sub-sample in the total sample. A statistical study is then carried out based on the two sums thus determined, then it is tested according to the hypothesis that the two subsamples do not belong to the same population. The Pettitt test is non-parametric and derives from that of Mann Whitney. The absence of a discontinuity in the series (Xi) of size N constitutes the null hypothesis. Its implementation supposes that for any instant T between 1 and N, the time series (Xi) 1 to t and t + 1 to N belong to the same population. The variable to be tested is the maximum in the absolute value of the variable Ut, N defined by:
$$\mathbf{U}\mathbf{t},\mathbf{N}= {\mathbf{Ʃ}}_{\mathbf{i}=1 }^{\mathbf{t}}{\mathbf{Ʃ}}_{\mathbf{i}=\mathbf{t}+1 }^{\mathbf{N}}\mathbf{D}\mathbf{i}\mathbf{j}$$
where Dij = Sign (Xi - Xj) with: sign (x) = 1 if x > 0, 0 if x = 0 and − 1 if x < 0. If the null hypothesis is rejected, an estimate of the date of discontinuity is given by defining the maximum in the absolute value of the variable Ut, N.
In addition to the Pettitt test, the Mann-Kendall test was also used to analyze precipitation, average river discharge and runoff coefficients (at annual and seasonal time steps). This test is based on the test statistic “S” defined as follow:
$$\mathbf{S}=\sum _{\mathbf{i}=1}^{\mathbf{n}-1}\sum _{\mathbf{j}=\mathbf{i}+1}^{\mathbf{n}}\mathbf{s}\mathbf{g}\mathbf{n}(\mathbf{x}\mathbf{j}-\mathbf{x}\mathbf{i})$$
Where the xj are the sequential data values, n is the length of the data set, and sgn = (θ) if θ > 1, 0 if θ = 0 and − 1 if θ < 0. There is no significant trend in the series analyzed when the calculated p-value is above the chosen significance level.
To assess the behavior of extreme flows, the Indicators of Hydrologic Alteration (IHA) tool, version 7.1, developed by The Nature Conservancy was used. This tool offers the possibility of comparing the parameters characterizing the flow regimes under different conditions (Richter et al. 1998). It uses daily discharge values and produces several important statistics. Only four of them were considered essential for this study, among which are the average, the coefficient of variation (CV) of extreme discharge and the Julian date of the annual minimum and maximum. By dividing the series of values in the period before and after the discontinuity, the tool calculates the change that occurred in the evolution of each of these parameters after the discontinuity. We can thus analyze not only the sign of change between the two periods but also the magnitude of the difference.
Landsat images were classified using the supervised maximum likelihood classification, using SNAP software (open access). This enabled us to perform a diachronic analysis of the evolution of land use in the basins studied. This operation was preceded by operations of preprocessing and recognition of objects in the field by photography and GPS (Global Positioning System). Satellite image preprocessing refers to all the process applied to raw data to correct geometric and radiometric errors that characterize certain satellite images. These errors are generally due to technical problems with the satellites and interactions between outgoing electromagnetic radiation and atmospheric aerosols, also called “atmospheric noise”. The atmospheric disturbances are influenced by various factors that are present on the day of acquisition, including weather, fires, and other human activities. They affect all the images acquired by passive satellites including Landsat 4-5-7 and 8. The downloaded Landsat images being orthorectified, the preprocessing involved atmospheric correction of the images and reprojection into the local system (WGS_84_UTM_Zone_32N). For this, neo-channels were created, to increase the readability of the data by enhancing certain properties less obvious in the original image, thus showing more clearly the elements of the scene. Three indices are therefore created, namely: the Normalized Difference Vegetation Index (NDVI, Eq. 1), the brightness index (BI, Eq. 2) and the Normalized Difference Water index (NDWI, Eq. 3). These indices respectively highlight vegetated surfaces, sterile (non-chlorophyllin) elements such as urban areas and water bodies. The formulae used in creating these indices are:
\(\mathbf{N}\mathbf{D}\mathbf{V}\mathbf{I}=\frac{\mathbf{N}\mathbf{I}\mathbf{R}-\mathbf{R}}{\mathbf{N}\mathbf{I}\mathbf{R}+\mathbf{R}}\) Eq. 1
\(\mathbf{B}\mathbf{I}={({\mathbf{R}}^{2}+{\mathbf{N}\mathbf{I}\mathbf{R}}^{2})}^{0.5}\) Eq. 2
\(\mathbf{N}\mathbf{D}\mathbf{W}\mathbf{I}=\frac{\mathbf{N}\mathbf{I}\mathbf{R}-\mathbf{M}\mathbf{W}\mathbf{I}\mathbf{R}}{\mathbf{N}\mathbf{I}\mathbf{R}+\mathbf{M}\mathbf{W}\mathbf{I}\mathbf{R}}\) Eq. 3
where NIR: ground reflectance of the surface in the near-infrared channel; R: ground reflectance of the surface in the red channel and MWIR: ground reflectance of the surface in the mid-wave infrared channel. Due to the fact that the study area extends over several scenes, the enhancement operations were followed by the mosaic of the different scenes used on each date. The use of Google Earth, as well as the spaces sampled from the GPS, made it possible to identify with certainty the impervious areas (buildings, savannas, bare soils, and crops), water bodies (large rivers, lakes and ponds) and forest (secondary, degraded, non-degraded and swampy) of each mosaic. Before the classification, the separability of the spectral signatures of the sampled objects to avoid interclass confusion was assessed by calculating the “transformed divergence” index. The value of this index is between 0 and 2. A value > 1.8 indicates a good separability between two given classes. The different classes used in this study show good separability between them, irrespective of the image considered, with indices > 1.9. The validation of the classifications obtained was carried out using the confusion matrix, making it possible to obtain treatment details to validate the choice of training plots. After validating the land use/land cover maps, the statistical and spatial differences of each class between 1973 and 2018 were evaluated.