Identification of Priority Areas for the Control of Soil Erosion and the Influence of Terrain Factors Using RUSLE and GIS in the Caeté River Basin, Brazilian Amazon

Abstract

Soil erosion poses a significant global environmental challenge, causing land degradation, deforestation, river siltation, and reduced agricultural productivity. Although the Revised Universal Soil Loss Equation (RUSLE) has been widely applied in Brazil, its use in the tropical river basins of the Amazon remains limited. This study aimed to apply a GIS-integrated RUSLE model and compare its soil loss estimates with multiple linear regression (MLR) models based on terrain attributes, aiming to identify priority areas and key geomorphometric drivers of soil erosion in a tropical Amazonian river basin. A digital elevation model based on Shuttle Radar Topography Mission (SRTM) data, land use and land cover (LULC) maps, and rainfall and soil data were applied to the GIS-integrated RUSLE model; we then defined six risk classes—slight (0–2.5 t ha−1 yr−1), slight–moderate (2.5–5), moderate (5–10), moderate–high (10–15), high (15–25), and very high (>25)—and identified priority zones as those in the top two risk classes. The Caeté River Basin (CRB) was classified into six erosion risk categories: low (81.14%), low to moderate (2.97%), moderate (11.88%), moderate to high (0.93%), high (0.03%), and very high (3.05%). The CRB predominantly exhibited a low erosion risk, with higher erosion rates linked to intense rainfall, gentle slopes covered by Arenosols, and human activities. The average annual soil loss was estimated at 2.0 t ha−1 yr−1, with a total loss of 1005.44 t ha−1 yr−1. Additionally, geomorphological and multiple linear regression (MLR) analyses identified seven key variables influencing soil erosion: the convergence index, closed depressions, the topographic wetness index, the channel network distance, and the local curvature, upslope curvature, and local downslope curvature. These variables collectively explained 26% of the variability in soil loss (R2 = 0.26), highlighting the significant role of terrain characteristics in erosion processes. These findings indicate that soil erosion control efforts should focus primarily on areas with Arenosols and regions experiencing increased anthropogenic activity, where the erosion risks are higher. The identification of priority erosion areas enables the development of targeted conservation strategies, particularly for Arenosols and regions under anthropogenic pressure, where the soil losses exceed the tolerance threshold of 10.48 t ha−1 yr−1. These findings directly support the formulation of local environmental policies aimed at mitigating soil degradation by stabilizing vulnerable soils, regulating high-impact land uses, and promoting sustainable practices in critical zones. The GIS-RUSLE framework is supported by consistent rainfall data, as verified by a double mass curve analysis (R2 ranging from 0.64 to 0.77), and offers a replicable methodology for soil conservation planning in tropical basins with similar erosion drivers. This approach offers a science-based foundation to guide soil conservation planning in tropical basins. While effective in identifying erosion-prone areas, it should be complemented in future studies by dynamic models and temporal analyses to better capture the complex erosion processes and land use change impacts in the Amazon.