Soil Erosion Assessment Using the RUSLE Model and Geospatial Techniques (Remote Sensing and GIS) in South-Central Niger (Maradi Region)

: A systematic method, incorporating the revised universal soil loss equation model (RUSLE), remote sensing, and the geographic information system (GIS), was used to estimate soil erosion potential and potential area in the Maradi region of south-central Niger. The spatial trend of seasonal soil erosion was obtained by integrating remote sensing environmental variables into a grid-based GIS method. RUSLE is the most commonly used method for estimating soil erosion, and its input variables, such as rainfall erosivity, soil erodibility, slope length and steepness, cover management, and conservation practices, vary greatly over space. These factors were calculated to determine their inﬂuence on average soil erosion in the region. An estimated potential mean annual soil loss of 472.4 t/ac/year, based on RUSLE, was determined for the study area. The potential erosion rates varied from 14.8 to 944.9 t/ac/year. The most eroded areas were identiﬁed in central and west-southern areas, with erosion rates ranging from 237.1 to 944.9 t/ac/year. The spatial erosion maps can serve as a useful reference for deriving land planning and management strategies and provide the opportunity to develop a decision plan for soil erosion prevention and control in south-central Niger.


Introduction
Soil erosion was known to be one of the most important environmental problems in the nineteenth and twentieth centuries, especially for developing country [1,2]. Soil erosion is one of the world's greatest scientific problems. It has various negative consequences, such as land degradation, decreasing water quality, sedimentation of rivers, and destruction of roads [3]. As well as these, it directly affects food production, human health, and the earth's environment. For example, the gradual degradation of habitat is mainly caused by extreme soil erosion [4,5]. Soil erosion is the detachment of soil and its movement and deposition by various methods [6]. The issue has far-reaching political, social, economic, and environmental consequences, as it causes damage, both on-site and off-site [7][8][9][10], in many countries [11,12]. Moreover, soil erosion leads to a decrease in soil fertility, which has become a threat to sustainable agricultural production and water quality in these regions [13,14]. Furthermore, surface water change is generally symbolic of global analysis variance in a research region, where hydrological process might better understand the natural and anthropogenic disruption of meteorological features, as well as projected trends in water management [15,16].
Erosion is the process of the detachment and transport of soil particles by erosion forces. These erosive forces can be wind, ice, or water in the form of raindrop action and surface runoff [17]. On the other hand, erosion by wind, water, ice, or gravity can be defined as the removal of the earth's surface material. The hazards associated with extremely rapid erosion have existed for over a million geological years in most regions of the earth. Nevertheless, the problem has been exacerbated in recent times by increasing human influence on the environment. Furthermore, data collection on soil erosion is typically capital intensive and time-consuming [18,19]. In addition, the effects of soil erosion appear to pose a serious threat to human existence [20,21], As one element of soil degradation, it leads to the deterioration of soil's physical, chemical, and biophysical characteristics.
Furthermore, various research measures of soil erosion have been established, and it is a common opinion that the area affected by erosion can be the best predictor for soil erosion research [22]. Apart from this, soil erosion is the most serious type of land loss that threatens the food security and natural economy of regions, both locally and globally [23,24]. Soil erosion in catchment areas and on the banks of reservoirs and lakes leads to siltation. This may result in increased turbidity, clogging of hydraulic structures, and reduced safety for flood retention [23]. Thus, soil erosion caused by anthropogenic activity in the watershed is a key factor in reservoir sedimentation. The silt is carried downstream by rivers, which eventually lead to reservoirs. The geology, geographical character, climate, plant types, and land use within the drainage basin influence the amount and kind of sediment [24]. According to Kastridis and Kamperidou (2015) [25], the issue of land use is very important and can influence the erosion process in a positive or negative way. Vacca et al. (2000) compared runoff and soil erosion measurements from three hillslopes with varying land uses and concluded that land use has a significant impact on runoff and erosion rate [26]. Additionally, soil erosion is a dynamic phenomenon, with numerous, and various unknown, impacts on the lives of people and the environment [27]. Decreased soil fertility leads to the degradation of natural habitats, such as forests, pastures, and rural ecosystems. Soil erosion is an important index for the evaluation and sustainability of soil conservation strategies [28,29]. Nowadays, poor soil conditions in an area have a negative impact on vegetation development, making the water system and vegetation susceptible [16].
An appropriate model of soil erosion can be applied to measure the soil erosion levels of a specific region. Erosion simulation can be used as a statistical tool to measure soil degradation and soil erosion potential for environmental planning [21]. Therefore, soil erosion has economic and political importance, as well as social and environmental impacts [30]. Brath et al. (2002) [31] Used a regional application of RUSLE to analyze the impact of land use on erosion, based on information about terrain morphology atmosphere, soil, and land use. Several studies have assessed soil erosion losses in different regions of the world using remote sensing and GIS technology [32][33][34][35]. Their analyses showed that changes in land use seriously affected annual soil loss. On the other hand, some studies have attempted to quantitatively evaluate these changes and concluded that the human influence in driving changes in sediment output is more important than climate change [4].
Moreover, determining the economic value of soil erosion requires field knowledge. Many models have been developed to estimate soil loss, such as the universal soil loss equation (USLE) [36] and the revised universal soil loss equation (RUSLE) [37].Several researchers combined erosion models with remote sensing tools and geographic information systems (GIS) [38][39][40] to improve the accuracy of the erosion estimation and to predict spatial maps of soil erosion risk. Moreover, soil erosion is a natural process of soil material removal and transport due to the action of erosive agents such as water, wind, gravity, and human disturbance [41,42]. However, when soil erosion occurs faster because of human disturbance, it has negative ecological and economic impacts. Modelling for erosion in Niger is difficult, due to a lack of data, especially the data related to soils. Key parameters of erosion prediction models, such as soil erodibility, expressed as k-factor, are difficult to obtain [43,44]. While, wind and water erosion are commonly believed to be prevalent in the Sahel, there is no evidence of the rate of soil erosion or the impact on agriculture [43]. To remedy this, indirect approaches are used to estimate these factors to make these studies feasible [45]. Seeking an appropriate model that can accurately measures soil loss from all types of erosion at a given site remains a challenge.
There are several soil erosion models, with varying degrees of complexity. One of the most widely used empirical models for assessing sheet and rill erosion is USLE, developed by Wischmeier and Smith in 1965. USLE was originally developed primarily for estimating soil erosion in cropland or gently sloping topography, and has revised (RUSLE) and modified (MUSLE) versions [46,47].
In this context, the satellite remote sensing can be a convenient way to solve the problem. Remote sensing (RS) provides data over vast regions, with the ability to visit the same land area regularly. Furthermore, it can make an important contribution to regional soil erosion risk assessment through the rapid derivation of appropriate indices [48][49][50]. GIS and remote sensing are widely applied in many research areas, mainly in flood risk assessment, which is directly associated with soil erosion, as extreme hydrological events generate soil erosion and sediment transport. According to Tzioutzios, 2020, by taking advantage of the advanced capabilities offered by GIS, this task can be expedited, so that results and decisions are less time consuming, more efficient, and more accurate [50]. The use of GIS and spatial data plays a leading role in this area; by applying simple spatial techniques to the datasets and combining them with the RUSLE model. GIS and RS allow the use of spatial data of different types. Confirming this, recently made available datasets and technologies (satellite data and GIS) provide an adequate spatial and temporal accuracy, which ultimately leads to better flood risk reduction [51].Additionally, there have been some published studies on the application of RS and GIS technologies for soil erosion modelling in other parts of the world [52], yet little or no work was found on the application of GIS for erosion modelling in south-central Niger (Maradi region).
Therefore, this study aimed to estimate soil erosion potential using RUSLE, RS, and GIS technology in south-central Niger. The specific aims were to: (a) create a method that integrates remote sensing data and GIS with the RUSLE to estimate the spatial distribution of soil erosion of the overall study area, (b) provide a complete map of soil erosion susceptibility and land use/cover changes using RS and GIS, and (c) identify the areas of high and low soil erosion in the study area using the overlay approach.
The spatial erosion maps produced using the RUSLE model, by combining RS and GIS technology, can serve as a useful input for deriving land planning and management strategies and provide an opportunity to develop a decision plan for soil erosion prevention and control in south-central Niger.

Description of the Study Area
The Maradi Region is situated between 13 • and 15 • 26 north latitude and 6 • 16 and 8 • 36 east longitude in the south-central part of Niger ( Figure 1). It is located in central Niger between the Damagram and Ader-Doutchi Rivers, with dune soils (soils on dunes). The area is bordered by Zinder, Tahoua, Agadez, and the Federal Republic of Nigeria. The region covers 41,796 km 2 , or 3.30 percent of the country. These are the land types: agriculture (71.5%), pasture (25%), and woodland (1.463) (3.5 percent).
Maradi has two different climate zones. The northern Sahelian Climate has an average annual rainfall of 200-300 mm. The southern regions have a Sahelian-Sudanese climate with 500-600 mm annual rainfall. It is dry and cold from October to February, with mean annual temperatures of −10 • C; dry and hot from March to May, with temperatures reaching 40 • C; and wet and chilly from June to September. This area also has a rainfall gradient that diminishes from south to north. The average annual rainfall for the previous ten years was 493.3 mm in the south (Madarunfa rainfall station) and 342.7 mm in the north (Fako rainfall station) [53]. for the previous ten years was 493.3 mm in the south (Madarunfa rainfall station) and 342.7 mm in the north (Fako rainfall station) [53].

Data Source Processing
RUSLE has been applied in a GIS approach in a variety of contexts, including mountainous tropical watersheds, large-scale watersheds, agriculturally dominant watersheds, areas with distinct wet and dry seasons, and areas with dynamic changes in land cover patterns, agricultural farmlands, and development. The RUSLE model is composed of three primary databases: (1) Climatic and survey database, including monthly temperature and precipitation data, as well as contours necessary for computing the erosivity factor, slope length, and steepness factors (LS). (2) The crop database stores data essential for calculating the surface cover factor (C). (3) The soil data set comprises soil survey and characterization information, which is used to calculate the soil erodibility factor (K).
The RUSLE model determines the yearly average loss of soil erosion, by taking into account the five elements indicated in Equation (1) [36].
where A is the approximate average soil loss in tonnes per hectare per year (t ha −1 year −1 ), R is the measure of rainfall erosivity factor expressed in (MJ.mm ha −1 h −1 year −1 ), K is the soil erodibility index in (t.hr/MJ.mm), LS is the topographic factor = slope length and slope steepness factor (dimensionless), C is the cover management factor = cropping management factor (dimensionless), and P is the support measures/practice control factor (dimensionless). The RUSLE approach is often used to estimate each element in the model. Previous scholars established a variety of methods for predicting these parameters, including the use of climatic data, soil and geological maps, remotely sensed satellite pictures, empirical formulae, and digital elevation models (DEM) derived from a variety of sources. The

Data Source Processing
RUSLE has been applied in a GIS approach in a variety of contexts, including mountainous tropical watersheds, large-scale watersheds, agriculturally dominant watersheds, areas with distinct wet and dry seasons, and areas with dynamic changes in land cover patterns, agricultural farmlands, and development. The RUSLE model is composed of three primary databases: (1) Climatic and survey database, including monthly temperature and precipitation data, as well as contours necessary for computing the erosivity factor, slope length, and steepness factors (LS). (2) The crop database stores data essential for calculating the surface cover factor (C). (3) The soil data set comprises soil survey and characterization information, which is used to calculate the soil erodibility factor (K).
The RUSLE model determines the yearly average loss of soil erosion, by taking into account the five elements indicated in Equation (1) [36].
where A is the approximate average soil loss in tonnes per hectare per year (t ha −1 year −1 ), R is the measure of rainfall erosivity factor expressed in (MJ·mm ha −1 h −1 year −1 ), K is the soil erodibility index in (t·hr/MJ·mm), LS is the topographic factor = slope length and slope steepness factor (dimensionless), C is the cover management factor = cropping management factor (dimensionless), and P is the support measures/practice control factor (dimensionless). The RUSLE approach is often used to estimate each element in the model. Previous scholars established a variety of methods for predicting these parameters, including the use of climatic data, soil and geological maps, remotely sensed satellite pictures, empirical formulae, and digital elevation models (DEM) derived from a variety of sources. The following parts of this article detail the procedures used to construct the model factors and the findings received.
Each of the factors was derived separately in raster data format, and erosion was determined using the algebra map function. The ArcGIS Spatial Analyst method (10.3) was used to delineate the boundary of the study area.
The RUSLE parameters were determined using a separate equation, whose input was obtained from satellite imagery and DEM. The inputs, their origin, and the equations used are listed in Table 1. Equations available in the literature for calculating the variables were iteratively reviewed and the best equations were selected based on their suitability for use with the available data and their ability to provide estimates comparable to the reported field-based erosion measurements. The estimation of each factor is discussed in more detail in the following sections, see Figure 2. following parts of this article detail the procedures used to construct the model factors and the findings received. Each of the factors was derived separately in raster data format, and erosion was determined using the algebra map function. The ArcGIS Spatial Analyst method (10.3) was used to delineate the boundary of the study area. The RUSLE parameters were determined using a separate equation, whose input was obtained from satellite imagery and DEM. The inputs, their origin, and the equations used are listed in Table 1. Equations available in the literature for calculating the variables were iteratively reviewed and the best equations were selected based on their suitability for use with the available data and their ability to provide estimates comparable to the reported field-based erosion measurements. The estimation of each factor is discussed in more detail in the following sections, see Figure 2.  NASA and the Japan Aerospace Exploration Agency (JAXA) [57] Rainfall Erosivity (R factor): rainfall erosivity is the most important factor in RUSLE and responsible for soil erosion in an area. It is a long-term average of rainfall erosivity. Thus, the R-factor describes the intensity of precipitation in a given location, based on the extent of soil erosion [58,59], which is essential for assessing soil erosion risk under future land use and climate change [59]. The R factor quantifies the influence of precipitation on the amount and rate of runoff factor [50,51]. In the present analysis, the rainfall erosivity Rainfall Erosivity (R factor): rainfall erosivity is the most important factor in RUSLE and responsible for soil erosion in an area. It is a long-term average of rainfall erosivity. Thus, the R-factor describes the intensity of precipitation in a given location, based on the extent of soil erosion [58,59], which is essential for assessing soil erosion risk under future land use and climate change [59]. The R factor quantifies the influence of precipitation on the amount and rate of runoff factor [50,51]. In the present analysis, the rainfall erosivity was calculated using the mean annual rainfall data collected from two monitoring stations in the Maradi region (Madarunfa rainfall station in the south of the region and Fako rainfall station in the north). In this section we used the Renard and Freimund (1994) equation, which establishes a relation between (R) and the average annual rainfall. Furthermore, it can be used for locations with a mean annual precipitation of less than 850 mm. The following formula was used to calculate the R factor: where R is the measure of rainfall erosivity expressed in (MJ·mm ha −1 h −1 year −1 ) in the RUSLE equation and P is the average annual precipitation (mm). Soil Erodibility (K-factor): The K-factor reflects the susceptibility of the soil type to erosion, the transportability of sediments, and the amount and velocity of runoff for an appropriate rainfall input [32]. The soil maps of the region were prepared using FAO Digital Soil Map of the Word Shapefile (DSMW), where the zone of interest was captured over the respective area and the soil data layer was clipped to the study area in ArcGIS environment 10.3. The soil maps and K-factor values obtained from the literature for each soil type were then used to calculate the soil classes for the study area and then entered into ArcGIS and analyzed to obtain a K-factor map. The K-factor values ( Table 2) for the corresponding soil groups were used in the analysis. The K-factor represents the sensitivity of a soil to the isolation and transfer of soil particles. It reflects the combined effect of the soil's properties and shows the general tendency for erosion of a particular soil type [29]. Soils that have different properties have different resistances to erosion. Erodibility, therefore, varies with the physical and biochemical properties of the soil, such as soil texture and structure, aggregate stability, shear strength, infiltration capacity for organic matter, and chemical content. The respective K-factor values for each soil type were assigned based on the existing soil types in the study area [60,61].
Slope length and steepness (LS-factor): slope factor and steepness (LS) are equivalent to topographic factor and relief factor, respectively. This factor is a combined factor of slope length and slope steepness. Furthermore, the slope length factor L calculates the erosion effect of slope length and the slope steepness factor S calculates the erosion effect of slope steepness. The parameters associated with topography were derived from the digital elevation model (DEM). To create the LS-factor map, flow direction and flow accumulation maps were processed after the filling process and created from DEM data in the GIS extension Arc Hydro Tools, to be used as input for the LS-factor calculation. To better understand this, the FA were obtained as flows, by calculating the flow direction with DEM. For each cell, the flow accumulation was determined by the flow that flows through that cell; if the flow accumulation value is greater, the area will form a runoff more easily.
The LS values can be estimated as follows: where LS: slope length and slope steepness factor; FA: flow accumulation Land cover management (C-Factor): According to Molla et al. [29], the C-factor quantifies the cumulative effects on land degradation of trees, crop sequences, and other land cover areas, with values ranging from 0.01 to 1. The C-factor is dimensionless. Although C-factor values ( Table 3) were used in this analysis. The C-factor maps were quantified from the land use/land cover groups of the reservoir catchments. To determine the C-factor it is necessary to rely, as much as possible, on current land use/land cover data, which reflect the present status of the research area [17]. The C-factor is the most important factor in crop management and, perhaps, the most important factor in the RUSLE model, because it is the easiest aspect to manage for erosion mitigation. Although the corresponding proposed C-factor values for different land use/land cover groups were assigned based on the categorized map, these values were taken from previous studies and distributed among the respective land use/land cover forms. Finally, the C-factor map was created using the GIS database. Conservation practice (P-factor): the ratio between the soil loss of a field with a particular conservation method and a field where no conservation is practiced is expressed by the P-factor. It is a measure of the success of land management practices in reducing soil degradation in an area [62,63]; because of the erosion capacity of runoff, by manipulating the runoff conditions, concentration, velocity, and hydraulic forces of runoff, the management activities reduce the rate of soil erosion [32]. Thus, all agricultural activities that can discourage or reduce soil erosion were considered in this review. For further processing, the MODIS MCD12Q1 images were categorized using the maximum probability classifier and the p-value pixels were assigned to each land use and land cover class.

Limitations and Uncertainties
This section discusses a few of the key limitations of our research methodology, including geographic applicability, model uncertainty, input data, and validation.
Many of the limitations and uncertainties in our research could be attributed to existing RUSLE formulations, including uncertainty associated with the model's simple empirical nature and numerous subcomponents; data availability issues; and the model's inability to account for soil loss due to gully erosion, mass sliding events, or the prediction of potential sediment yields in streams.
The RUSLE model ignores gully erosion. Identifying and measuring gullies in the watershed will improve the accuracy of soil loss estimation for future conservation practices, planning, and management. Determining the most sensitive input parameters that guide selective and targeted interventions was not carried out either [63]. Further research on the effect of soil loss on agricultural output and rural lives is needed to establish a relationship between soil loss and unsustainable soil management methods. The study may also not have fully appreciated the quantitative aspects of model outputs, due to subjective community perceptions. In order to make sound decisions, we need more evidence.
The limitations of the USLE make it seem as if we lack a properly informed understanding of the extent of soil erosion in our region, let alone the mechanisms of sediment transport and deposition [64]. Indeed, the validation of the main result of the study was accomplished via the use of the RUSLE model in conjunction with local perspectives. As input parameters, soil, land use/cover, DEM, rainfall, and support practice data were employed. The raster layers were processed to display the relevant input parameters on the ArcGIS platform, and the inputs were then multiplied to calculate the yearly average rate of soil loss and create intensity maps for the research area.
After completing data entry operations and organizing each information layer, the mean yearly soil loss was estimated by multiplying each component layer together according to the RUSLE formula (Equation (1)) in a GIS framework, utilizing the raster calculator of map algebra functions and accompanying packages [65].

Soil Loss Rate Assessment
In this study, we used remote sensing, GIS techniques, and the RUSLE model to estimate the extent and spatial distribution of soil erosion in the study area. Furthermore, five erosion risk factors were calculated, including rainfall erosivity (R-factor), soil erodibility (K-factor), slope duration and steepness (LS-factor), land cover management (C-factor), and soil protection/conservation practices (P-factor).
Rainfall erosivity (R-factor): the calculated R-factor (Figure 3b) ranged from 0.011 to 1644.6 MJ mm/h ha yr. The average rainfall erosivity of the study area was 822.3 MJ mm/h ha yr. It was observed that rainfall is high in the southern part of the study area, as indicated by the results. In addition, the erosivity value of rainfall shows that the R-value is greater in a regions where the rainfall quality and intensity is higher and vice versa. Therefore, the different distributions of rainfall within the study area cannot be satisfactorily explained by a single R-value [66,67]. However, it showed that rainfall erosivity is an important factor in assessing soil erosion for future land-use/land-cover and climate change. Thus, this is due to changes in climate, interannual rainfall distribution, and rainfall intensity.
perspectives. As input parameters, soil, land use/cover, DEM, rainfall, and support practice data were employed. The raster layers were processed to display the relevant input parameters on the ArcGIS platform, and the inputs were then multiplied to calculate the yearly average rate of soil loss and create intensity maps for the research area.
After completing data entry operations and organizing each information layer, the mean yearly soil loss was estimated by multiplying each component layer together according to the RUSLE formula (Equation (1)) in a GIS framework, utilizing the raster calculator of map algebra functions and accompanying packages [65].

Soil Loss Rate Assessment
In this study, we used remote sensing, GIS techniques, and the RUSLE model to estimate the extent and spatial distribution of soil erosion in the study area. Furthermore, five erosion risk factors were calculated, including rainfall erosivity (R-factor), soil erodibility (K-factor), slope duration and steepness (LS-factor), land cover management (C-factor), and soil protection/conservation practices (P-factor).
Rainfall erosivity (R-factor): the calculated R-factor (Figure 3b) ranged from 0.011 to 1644.6 MJ mm/h ha yr. The average rainfall erosivity of the study area was 822.3 MJ mm/h ha yr. It was observed that rainfall is high in the southern part of the study area, as indicated by the results. In addition, the erosivity value of rainfall shows that the R-value is greater in a regions where the rainfall quality and intensity is higher and vice versa. Therefore, the different distributions of rainfall within the study area cannot be satisfactorily explained by a single R-value [66,67]. However, it showed that rainfall erosivity is an important factor in assessing soil erosion for future land-use/land-cover and climate change. Thus, this is due to changes in climate, interannual rainfall distribution, and rainfall intensity. Soil erodibility (K-factor): depending on the soil texture, the K-factor reflects the cumulative effects of the soil properties. However, it shows that a particular soil type is Soil erodibility (K-factor): depending on the soil texture, the K-factor reflects the cumulative effects of the soil properties. However, it shows that a particular soil type is generally resistant to erosion. The soil erodibility factor (K) was determined using the relationship proposed by Wischmeier et al. (1971) [68] (Figure 4b). Moreover, the value of the K element is influenced by the penetration potential of the soil. The spatial distribution of the K factor for the study area showed that soils with higher K values were more susceptible to soil erosion. Although, lower K values were less susceptible to soil erosion [59], ranging from 0.016 t ha h/ha MJ mm to 0.027 t ha h/ha MJ mm, indicating low to moderate soil erodibility.
generally resistant to erosion. The soil erodibility factor (K) was determined using the relationship proposed by Wischmeier et al., (1971) [68] (Figure 4b). Moreover, the value of the K element is influenced by the penetration potential of the soil. The spatial distribution of the K factor for the study area showed that soils with higher K values were more susceptible to soil erosion. Although, lower K values were less susceptible to soil erosion [59], ranging from 0.016 t ha h/ha MJ mm to 0.027 t ha h/ha MJ mm, indicating low to moderate soil erodibility. Slope length and steepness (LS-factor): the LS-factor is used to describe the influence of the topographic component on soil erosion. It considers the topographic aspect of a particular area for soil erosion. Here, the topographic aspect represents the susceptibility of the study area to topographic erosion. It represents the influence of slope length and steepness of the slope on the mechanism of erosion ( Figure 5b). River accumulation and percentage slope were used as inputs for estimating the topographic component. In the central and southern parts of the study area, the slope is high. Moreover, the slope maps identify areas with a high potential runoff and, therefore, high erosion rate; these potentials reduce as the slope degree drops from a high to a low slope degree (Figure 5a). Furthermore, the LS-factor ranged from 0 (lower part) to 111.565 (upper part). Higher LSfactor values prevailed in the central and eastern parts of the study area. The phenomenon indicates that the soil loss will increase with slope length and slope steepness. Soil loss increases faster with slope length than with slope steepness. A holistic analysis of all soil loss factors revealed that LS-factor has a significant effect on the estimation of total soil loss in the study area.
Land cover management (C-factor): the results of available vegetation cover for all land areas in the RUSLE model are numerically expressed by the C-factor, whether in agriculture or soil erosion mitigation measures. For the other main land uses in the LULC, bars were established in the area: shrublands, grasslands, croplands, urban and built-up, and forest, using the algorithm with the highest probability (Figure 6a). In the eastern and northern parts of the study area, the C-factor values were lower and the ground area was characterized by trees. However, in the southern and western parts of the zone, where a large area is characterized by bare soil and built-up areas, the C-factor values were higher. Slope length and steepness (LS-factor): the LS-factor is used to describe the influence of the topographic component on soil erosion. It considers the topographic aspect of a particular area for soil erosion. Here, the topographic aspect represents the susceptibility of the study area to topographic erosion. It represents the influence of slope length and steepness of the slope on the mechanism of erosion ( Figure 5b). River accumulation and percentage slope were used as inputs for estimating the topographic component. In the central and southern parts of the study area, the slope is high. Moreover, the slope maps identify areas with a high potential runoff and, therefore, high erosion rate; these potentials reduce as the slope degree drops from a high to a low slope degree (Figure 5a). Furthermore, the LS-factor ranged from 0 (lower part) to 111.565 (upper part). Higher LS-factor values prevailed in the central and eastern parts of the study area. The phenomenon indicates that the soil loss will increase with slope length and slope steepness. Soil loss increases faster with slope length than with slope steepness. A holistic analysis of all soil loss factors revealed that LS-factor has a significant effect on the estimation of total soil loss in the study area.
Land cover management (C-factor): the results of available vegetation cover for all land areas in the RUSLE model are numerically expressed by the C-factor, whether in agriculture or soil erosion mitigation measures. For the other main land uses in the LULC, bars were established in the area: shrublands, grasslands, croplands, urban and built-up, and forest, using the algorithm with the highest probability (Figure 6a). In the eastern and northern parts of the study area, the C-factor values were lower and the ground area was characterized by trees. However, in the southern and western parts of the zone, where a large area is characterized by bare soil and built-up areas, the C-factor values were higher. Therefore, a C-factor grid map was prepared using the LULC, and the corresponding C-factor values are shown in Figure 6b. Therefore, a C-factor grid map was prepared using the LULC, and the corresponding Cfactor values are shown in Figure 6b.  Conservation practices (P-factor): based on the current land management practices used in the south Central Niger, the P-factor value varied from 0.5 to 1 (Figure 7). Thus, 0.5 represents a very strong facility for protection against human erosion and 1 represents no facility for protection against human erosion [69,70]. Based on the results, lower Pfactor values characterized the core area of the study zone, and higher P-factors characterized the upper part of the study area. In addition, the relevance of conservation  Therefore, a C-factor grid map was prepared using the LULC, and the corresponding Cfactor values are shown in Figure 6b.  Conservation practices (P-factor): based on the current land management practices used in the south Central Niger, the P-factor value varied from 0.5 to 1 (Figure 7). Thus, 0.5 represents a very strong facility for protection against human erosion and 1 represents no facility for protection against human erosion [69,70]. Based on the results, lower Pfactor values characterized the core area of the study zone, and higher P-factors characterized the upper part of the study area. In addition, the relevance of conservation Conservation practices (P-factor): based on the current land management practices used in the south Central Niger, the P-factor value varied from 0.5 to 1 (Figure 7). Thus, 0.5 represents a very strong facility for protection against human erosion and 1 represents no facility for protection against human erosion [69,70]. Based on the results, lower P-factor values characterized the core area of the study zone, and higher P-factors characterized the upper part of the study area. In addition, the relevance of conservation measures increased soil erosion vulnerability. In the study area, the lower P-factor values were concentrated in the southern part of the study area, and the higher P-factor values were distributed in the upper and outer parts of the area. measures increased soil erosion vulnerability. In the study area, the lower P-factor values were concentrated in the southern part of the study area, and the higher P-factor values were distributed in the upper and outer parts of the area.

Estimation of Potential Soil Erosion
For the research region, GIS, RS, and erosion models were merged to assess erosion amplitude and spatial distribution. Rainfall erosivity, soil erodibility, slope length and steepness, land cover management, and soil conservation were all identified as potential erosion factors. Figures 3-7 demonstrate the outcomes of the modelling of these variables.
It is clear from the modelling of rainfall erosivity, that (as shown in Figure 3) the stronger and longer a rainstorm lasts, the larger its erosion potential. Modeling the length and steepness of the slope yielded the results shown in Figure 5. The greater slope there is, the more runoff there will be in the long term. It is clear from Figure 5a that as the slope increases, so do the runoff velocities that cause erosion. This was confirmed by Wischmeier and Smith [69]. This component was also modelled, with the results presented in Figure 6, which shows how this affected land cover. According to the findings, plant cover has an effect on erosion because it intercepts rainwater, lowering its energy and boosting rainfall penetration. Figure 4b shows the soil erodibility factor, which takes into account both soil erosion susceptibility and runoff volume and rate. Therefore, this illustrates how much the soil changes per unit of applied external force of energy, depending on whether water splashes or surface flows detach it.
As a result of the final soil loss model, Figure 8 shows the annual erosion map of the study area, which helps to identify the areas prone to soil erosion. This finding was supported by data from field studies in our regions with a similar environmental condition [71], where there is a high level of erosion that affects the soil and this many detrimental effects due to downstream sedimentation. Table 4 depicts the erosion susceptibility based on LULC and slope categories, in addition to the five erosion factors that were modelled. The RUSLE model estimated yearly soil loss in the study region at 944.9 t/year, which is clearly seen in the findings

Estimation of Potential Soil Erosion
For the research region, GIS, RS, and erosion models were merged to assess erosion amplitude and spatial distribution. Rainfall erosivity, soil erodibility, slope length and steepness, land cover management, and soil conservation were all identified as potential erosion factors. Figures 3-7 demonstrate the outcomes of the modelling of these variables.
It is clear from the modelling of rainfall erosivity, that (as shown in Figure 3) the stronger and longer a rainstorm lasts, the larger its erosion potential. Modeling the length and steepness of the slope yielded the results shown in Figure 5. The greater slope there is, the more runoff there will be in the long term. It is clear from Figure 5a that as the slope increases, so do the runoff velocities that cause erosion. This was confirmed by Wischmeier and Smith [69]. This component was also modelled, with the results presented in Figure 6, which shows how this affected land cover. According to the findings, plant cover has an effect on erosion because it intercepts rainwater, lowering its energy and boosting rainfall penetration. Figure 4b shows the soil erodibility factor, which takes into account both soil erosion susceptibility and runoff volume and rate. Therefore, this illustrates how much the soil changes per unit of applied external force of energy, depending on whether water splashes or surface flows detach it.
As a result of the final soil loss model, Figure 8 shows the annual erosion map of the study area, which helps to identify the areas prone to soil erosion. This finding was supported by data from field studies in our regions with a similar environmental condition [71], where there is a high level of erosion that affects the soil and this many detrimental effects due to downstream sedimentation. Table 4 depicts the erosion susceptibility based on LULC and slope categories, in addition to the five erosion factors that were modelled. The RUSLE model estimated yearly soil loss in the study region at 944.9 t/year, which is clearly seen in the findings (Figure 8). There was potential soil loss in the five classes, ranging from 14.8 to 944.9 t/ac/year. About 12.5% of this region's land has been lost to erosion, or 472.4 t/ac/year. Figure 8 shows the findings, which show that the soil erosion was minimal in all LULC classes. However, even in the somewhat forested areas, there are signs of significant erosion. Visual evaluation of the different components revealed that this was related to the high erodibility of the soil group that crosses the moderate forest LULC Class sections. There is also a considerable probability of excessive erosion in the moderate forest cover class, since it is located in the research area's heavy rainfall zones. Table 4 shows the soil erosion classes and range of the study area.

Conclusions and Recommendations
This paper demonstrates the use of an empirical soil erosion model such as RUSLE in conjunction with RS and GIS technology to estimate the soil erosion potential and potential zones in the Maradi region. Additionally, an effort was made to investigate the effect of changing land use/land cover on erosion rate. The analysis and results indicate that the annual average soil loss in the study area was estimated at 944.9 t/year when the  This explains why the land is so resistant to erosion. There is additional evidence to suggest that steep slopes erode at a considerably greater rate than level areas. As can be seen, the average rate of soil loss and the contribution to total soil loss on steeper slopes are much greater than on mild slopes. Figure 8 shows the findings, which show that the soil erosion was minimal in all LULC classes. However, even in the somewhat forested areas, there are signs of significant erosion. Visual evaluation of the different components revealed that this was related to the high erodibility of the soil group that crosses the moderate forest LULC Class sections. There is also a considerable probability of excessive erosion in the moderate forest cover class, since it is located in the research area's heavy rainfall zones. Table 4 shows the soil erosion classes and range of the study area.

Conclusions and Recommendations
This paper demonstrates the use of an empirical soil erosion model such as RUSLE in conjunction with RS and GIS technology to estimate the soil erosion potential and potential zones in the Maradi region. Additionally, an effort was made to investigate the effect of changing land use/land cover on erosion rate. The analysis and results indicate that the annual average soil loss in the study area was estimated at 944.9 t/year when the RUSLE model was used. Additionally, it was observed that the amount of erosion varies significantly, according to the topography and land use/land cover. Furthermore, the results show that slope length, steepness, and the rainfall erosivity component were the most significant variables in preventing soil erosion. The climate component was often very effective in spatial soil erosion and is a major element in climatic conditions in the RUSLE model of rainfall erosivity. Overall, areas with natural forest cover in the headwater regions had the lowest rate of soil erosion, at 14.82 t/year, whereas places with human intervention had the greatest rate, at 944.9 t/year. The highest probability approach detected shrublands, grasslands, croplands, urban and built-up, and forest areas in the study area. From 14.8 to 944.9 t/ac/year, the five classes had potential soil loss. The mean soil loss in the research region was 472.4 t/ac/year, which represents 12.5 percent of the entire area of this most eroded part. This research highlights the use of RS and GIS in evaluating soil erosion and assisting with erosion loss estimates.
Furthermore, high and severe erosion were mostly identified in high slope gradient zones and intermediate forest LULC areas. The findings showed that sites in the moderate forest LULC class had significant erosion potential, due to an intersecting high erodibility soil group. Since vegetation is a key component in soil erosion, future studies should use MODIS MCD12Q1 images obtained from more recent high-resolution satellite imagery. This will increase the accuracy of the LULC maps and DEMs used to calculate land slopes. Furthermore, validating the soil erosion loss using local data helps improve the quality and accuracy of the findings.
Additional research is required to obtain the right P-factor values within the study region, to accurately assess possible soil erosion. Overall, the findings of this research can assist in conducting proper erosion control measures in the seriously impacted areas. The data collected from this research may help in establishing management scenarios and give alternatives to policy makers for controlling soil erosion threats in the most effective way.

Data Availability Statement:
The data supporting our analyses can be accessed publicly from the links provided in the paper, and we will also share the data soon after the manuscript is published.