3.2. Data Analysis Methods
The RUSLE model, widely used [
17] by researchers, which is the revised version of the universal soil loss equation, was chosen to quantify mean annual water erosion in the Mékrou watershed. This choice was based above all on the availability of data for this watershed, on climatic, topographical and soil parameters, and on land use and conservation methods. It is presented in the form of equation (1):
where A - is the estimated annual soil loss (t/ha/yr); R - is the rainfall erosivity factor (MJ.mm/ha.h.y); K - is the soil erodibility factor (t.ha.h/ha.MJ.mm); LS - is the topographical factor (dimensionless); C - is the vegetation factor (dimensionless); and P - is the soil conservation and protection practices factor (dimensionless).
Each of the components of formula (1) is estimated independently according to the following criteria and formulas:
- The precipitation erosivity factor (R):
This factor is an indicator of the potential of precipitation to induce soil loss through the detachment and transport of particles. It is a function of several precipitation-related parameters, including quantity, energy, duration, intensity, drop size, type of precipitation and the resulting runoff rate [
1,
18] . However, in order to enable its estimation in view of the lack of data, in particular rainfall intensity [
19], following his research in some countries in Africa (Côte d'Ivoire, Burkina Faso, Senegal, Niger, Chad, Cameroon and Madagascar), proposed the following linear equation, as a function of annual rainfall alone:
Specifically for the West Africa sub-region, it is suggested taking the dividend α = 0.5 (i.e. R = 0.5* Pl; where Pl is the annual rainfall in mm).
- The soil erodibility index (K):
This index takes values between 0 and 1. The less sensitive the soil is to water erosion, the more K tends towards 0 [
1]. Originally, the K factor was determined experimentally as a function of soil texture and structure, organic matter content and permeability [
20,
21,
22]. In this study, however, the Williams [
23] equation was used (Formula 3), which relates the K factor to soil properties (clay, sand and silt content, and organic carbon fraction) available in the FAO database for the Mékrou watershed (
Table 1).
where the components Fcsand, Fcl-si, Forgc, and Fhinsand characterize the composition of the soil in terms of sand, clay/silt, organic and carbon respectively. They are calculated by the following formulae 4 to 7:
where m
s is the percentage of sand (diameter 0.05-2.00 mm); m
silt, the percentage of silt (diameter 0.002-0.05 mm); m
c, the percentage of clay (diameter < 0.002 mm); and orgC,the organic carbon content (%).
- The topographical factor (LS) :
The topographic factor LS indicates the gradient responsible for runoff velocity. The steeper the slope of the land, the higher the velocity and erosive potential of the runoff [
24,
25]. The Formula (8) of Moore and Burch [
26] below was used to calculate the LS factor:
where flow accumulation represents the cumulative upstream contributing area for a given cell, Cell size is the size of the pixel which is 30 m resolution, 22.13 is a uniform length (m) to avoid uncertainties related to the influence of the length of the slope on erosion, Sin Slope is the sine of the slope and 0.0896 is a constant.
- The vegetation factor (C):
This factor represents the effect of cultivation and land management practices on erosion. It is closely linked to land use types and therefore depends mainly on vegetation type, growth stage and vegetation cover rate. The C factor was estimated using the equation of Durigon [
27] , which will subsequently be improved by Colman [
28] in certain regions as a function of the NDVI value (Formula 9).
where
, NIR is the near infrared part of the electromagnetic spectrum and RED is the red part of the electromagnetic spectrum.
The Normalised Difference Vegetation Index (NDVI), also known as relative biomass, is a standardized index that allows users to generate images displaying greenness. NDVI uses the properties of red (RED) and near infrared (NIR) light to quantify vegetation, which can be useful for classifying high and low vegetation and indicating drought.
In the current study, NDVI is established using Landsat images introduced in the ARCGIS 10.8 tool and assigned band values according to the Landsat image type.
For the Landsat 4-5 and 7 images, band 4 is used for the near infrared (NIR) and band 3 for (RED), which is the red part of the electromagnetic spectrum.
As for the Landsat 8 images, band 5 is employed for the near infrared (NIR) and band 4 for (RED), which is the red part of the electromagnetic spectrum.
- The soil conservation and protection practices (P):
The values of the P factor, ranging from 0 to 1, depend on the soil conservation and protection practices implemented to mitigate soil erosion (
Table 2), such as contour farming, terraces and strip cropping [
29]. The closer P is to 0, the greater the soil's resistance to human-induced erosion; conversely, the closer it is to 1, the lower the soil's resistance.
By combining all the information related to the quantification of water erosion in the Mékrou watershed, the latter was calculated for the study period (1981-2020) and for the decades 1981-1990, 1991-2000, 2001-2010 and 2011-2020. Water erosion was mapped for each of the 10-year sub-periods, and a cross-comparison was made to analyze the spatiotemporal dynamics of erosion at the scale of the Mekrou watershed.