The highest long-term values were concentrated in northern Kerala, while the smallest erosion-potential domain covered much of southern and parts of central Kerala. Rainfall erosivity describes the capacity of rainfall to detach and transport soil, while erosivity density expresses that erosive potential relative to the amount of rainfall. The two measures help distinguish the effect of how much rain falls from the effect of how intensely it falls.

For the analysis, the researchers converted daily gridded rainfall observations into monthly totals for each year and applied the Arnoldus empirical approach, a monthly-rainfall method widely used for long-term regional assessments. They calculated annual erosivity density as annual rainfall erosivity divided by annual rainfall. Trend significance was evaluated with the modified Mann–Kendall test, which adjusts for serial correlation, and trend magnitude was estimated with Sen’s slope method.

Agglomerative hierarchical clustering with Ward’s linkage then grouped grid cells with similar 120-year mean rainfall erosivity and erosivity density. The analysis identified a pronounced north–south gradient. Northern Kerala contained the highest minimum, maximum and mean rainfall erosivity, as well as the highest erosivity-density values.

The largest values occurred in and around Kasaragod, Kannur, Wayanad, Kozhikode and parts of eastern Palakkad, with additional high-erosivity areas in parts of Thrissur, Ernakulam and Kottayam. Most of central and southern Kerala had low to moderate values. The mapped patterns broadly followed rainfall totals, but rainfall erosivity had a higher coefficient of variation than rainfall, indicating that the erosive potential of Kerala’s rainfall was more variable across space and time than precipitation amount alone.

Four statistically distinct erosivity domains emerged from the clustering. The largest, C1, contained 23 cells and represented comparatively low conditions, mainly across southern and parts of central Kerala. Its mean rainfall erosivity ranged from 1,291 to 7,044 MJ mm ha⁻¹ h⁻¹ y⁻¹, and its erosivity density ranged from 1.34 to 3.38 MJ ha⁻¹ h⁻¹.

The intermediate C2 domain comprised 17 cells across much of central Kerala, with rainfall erosivity from 8,375 to 12,685 MJ mm ha⁻¹ h⁻¹ y⁻¹ and erosivity density from 2.45 to 3.99 MJ ha⁻¹ h⁻¹. Eight northern cells formed the high-erosivity C3 domain, reaching 21,643 MJ mm ha⁻¹ h⁻¹ y⁻¹ for rainfall erosivity and 5.76 MJ ha⁻¹ h⁻¹ for erosivity density. Two northern cells formed C4, the extreme domain, with rainfall erosivity of 26,974 to 29,987 MJ mm ha⁻¹ h⁻¹ y⁻¹ and erosivity density of 6.54 to 7.38 MJ ha⁻¹ h⁻¹.

Long-term averages nevertheless concealed substantial temporal variability. The 1991–2000 decade showed the widest spread in rainfall erosivity, from 1,269 to 38,824 MJ mm ha⁻¹ h⁻¹ y⁻¹, and the widest spread in erosivity density, from 1.26 to 10.53 MJ ha⁻¹ h⁻¹. Those ranges did not coincide with the highest average values: the decade with the greatest mean rainfall erosivity was 1921–1930, at 12,224 MJ mm ha⁻¹ h⁻¹ y⁻¹, and the highest mean erosivity density was 3.74 MJ ha⁻¹ h⁻¹.

Inter-decadal comparisons showed alternating periods of increase and decrease rather than a smooth progression. Rainfall erosivity rose most strongly between the second and third decades, while pronounced decreases occurred in several later transitions, particularly between 1991–2000 and 2001–2010. Comparing the first and last decades likewise showed that declines predominated statewide.

Statewide linear analysis confirmed the general direction of change. Mean annual rainfall decreased at a rate of 3.53 mm per year, with an R² of 0.1119. Mean rainfall erosivity declined by 34.34 MJ mm ha⁻¹ h⁻¹ y⁻¹ per year, with an R² of 0.1302, while erosivity density fell by 0.0074 MJ ha⁻¹ h⁻¹ per year, with an R² of 0.0947.

Grid-level results were more varied. For rainfall erosivity, 39 of 50 cells showed decreasing trends, including 33 statistically significant declines; 10 showed increases, five of them significant, and one showed no trend. Sen’s slopes ranged from −171.72 to 38.79 MJ mm ha⁻¹ h⁻¹ y⁻¹ per year.

For erosivity density, 40 cells declined, including 34 significant decreases, nine increased, including four significant increases, and one showed no trend. The strongest declines occurred in northern districts, whereas increases were confined to a few eastern highland locations. Elevation did not explain the regional pattern by itself.

In an analysis using 300 spatially distributed points sampled from interpolated surfaces, elevation had a moderate negative relationship with mean rainfall, with r = −0.49 and R² = 0.2418. Its relationships with rainfall erosivity and erosivity density were much weaker: r = −0.24 and R² = 0.0589 for erosivity, and r = −0.18 and R² = 0.0325 for erosivity density. Low-elevation locations generally received more than 3,000 mm of annual rainfall and often had greater erosivity, whereas high-elevation sites commonly received 1,000–1,500 mm, but considerable variation occurred within every elevation band.

The results indicate that monsoon circulation, rainfall intensity, seasonal distribution, topographic position and local atmospheric conditions jointly shaped erosivity rather than altitude acting as a simple control. The coexistence of declining long-term averages and episodic high-erosivity periods is an important feature of the record. A reduction in mean rainfall erosivity does not necessarily imply a proportional reduction in erosion hazard, because average conditions and extreme events represent different parts of the rainfall regime.

The monthly-rainfall method also cannot resolve the kinetic energy, duration or within-storm intensity of individual short-duration storms. Although this approach provides a consistent century-scale comparison, it may under- or overestimate local annual erosivity where extreme storms make a disproportionate contribution. The estimates describe climatic potential for soil detachment, not observed soil loss, which also depends on soil properties, vegetation, land use, slope, infiltration and sediment connectivity.

Likewise, the annual erosivity-density calculation describes long-term rainfall aggressiveness per unit rainfall rather than seasonal concentration or event-scale behavior. The spatial divergence also carries implications for erosion modelling. Because Kerala’s rainfall erosivity varies sharply over a comparatively small and topographically diverse area, applying one statewide R-factor could obscure localized high-erosivity conditions.

The four-domain classification provides an alternative framework for parameterizing Revised Universal Soil Loss Equation models, prioritizing watersheds and supporting spatially differentiated conservation planning. It does not by itself establish current erosion rates, however, and the four-cluster solution was selected through dendrogram interpretation rather than quantitative validation with indices such as silhouette width. High-temporal-resolution rainfall observations, climate projections and further analysis of extreme events will be needed to clarify future erosion susceptibility, but the current evidence shows that Kerala’s erosive rainfall has been highly variable despite its broad long-term decline.

Subject of Research: Century-scale analysis of rainfall erosivity and erosivity density across Kerala, India

Article Title: Century-scale spatiotemporal variability and trends of rainfall erosivity and erosivity density in an ecologically sensitive tropical monsoon region of India (1901–2020)

Article References: Valappil, N. K. M., Yusup, Y., Dewi, J. R., Huda, N., & Hamza, V. (2026). Century-scale spatiotemporal variability and trends of rainfall erosivity and erosivity density in an ecologically sensitive tropical monsoon region of India (1901–2020). Environmental Earth Sciences, 85(17), Article 437. https://doi.org/10.1007/s12665-026-13119-5

Image Credits: AI Generated

DOI: 10.1007/s12665-026-13119-5

Keywords: rainfall erosivity, erosivity density, Kerala, monsoon rainfall, soil erosion, Western Ghats, climate variability, Mann–Kendall test, hydroclimatology, RUSLE, spatial analysis, long-term trends

The results are model estimates for the catchment, rather than direct measurements of sediment exported from every slope. RUSLE estimates mean annual soil loss as the product of five factors: rainfall erosivity, soil erodibility, slope length, slope steepness, cover and management, and support practices. Rainfall erosivity, denoted R, was derived from 20 years of daily rainfall data in the NASA database and calculated with an empirical equation calibrated for East African highland conditions.

The resulting R values ranged from 527.507 to 888.624 megajoules per millimetre per hectare per hour per year, with the highest values in the eastern and southern catchment and lower values in central and northwestern areas. Station values were distributed across the landscape through inverse-distance weighting in ArcGIS, producing a continuous raster of rainfall erosivity. Soil erodibility, or K, was calculated from global FAO soil-texture and organic-carbon data.

It ranged from 0.239 to 0.688 tonnes per hectare-hour per hectare-megajoule per millimetre, with the highest values generally found in the south and northeast. These light-textured sandy and silty soils with relatively little organic matter were modeled as more susceptible to detachment and transport than clay-rich or organically enriched soils in the central and northwest. The K layer was constructed pixel by pixel with a raster calculator, while a 30-metre digital elevation model supplied the terrain information needed to estimate the combined slope-length and slope-steepness factor.

The study therefore combined regional climate records, global soil information and digital topography rather than relying on a single erosion gauge. The topographic factor, LS, ranged from zero to 1.678. It increased on long, steep slopes, where runoff has more opportunity to gather speed and transport detached particles, particularly in the northeastern and central catchment.

The cover-management factor, C, was derived from Landsat 8 near-infrared and red reflectance through the Normalized Difference Vegetation Index, or NDVI. Values ranged from 0.0004 in dense forest to 0.5237 in bare or degraded areas. This more than 1,300-fold range made vegetation cover one of the model’s most spatially variable controls.

The support-practice factor was fixed at 1.0 because the study found few formal soil-conservation structures in the catchment, so the mapped losses represent conditions without effective terracing, contour farming or comparable erosion-control practices. Multiplying the five raster layers produced the study’s estimated soil-loss range of zero to 51.4391 tonnes per hectare per year. The highest values occurred where the model combined high rainfall erosivity with high C and LS values, rather than in a single uniform zone.

In the northeastern Akanyaru Hills, annual cropland occupied steep, elongated slopes with locally more than 30-degree gradients. C values of roughly 0.25 to 0.52 represented limited protection during land preparation and early crop development, while R values reached the catchment’s upper range. The most degraded slope units in this zone had estimated losses approaching or exceeding 40 to 51.4 tonnes per hectare per year.

The model thus identifies the intersection of intense rainfall, topographic amplification and incomplete vegetation cover as a defining combination in the northeast. A second major hotspot lay in the southwest, where persistently bare and severely degraded land coincided with C values of approximately 0.45 to 0.524 and K values reaching 0.55 to 0.688. In model terms, this zone lacked even the partial canopy protection present in annual cropland, allowing rainfall and runoff to act directly on exposed mineral soil.

The southern forest margin formed a third vulnerable landscape: a shifting mosaic of remnant forest, newly cleared plots and pioneer grassland. Here, R values approached 888 megajoules per millimetre per hectare per hour per year, and clearing abruptly removed the very low C values associated with intact forest. These are spatial model patterns, not a measured time series of erosion before and after forest clearing.

Not every high-loss location depended on the same combination of conditions. Smallholder fields in the Gisagara and Nyaruguru southern uplands sat on moderately sloping terrain, with LS values around 0.5 to 1.0, but underlain by the catchment’s most erodible sandy and silty soils. The study attributed this pattern to depleted organic carbon and weak aggregate stability, showing why land cover alone cannot provide a complete erosion assessment.

In the Bugesera and Cyohoha transitional zone, modeled per-hectare losses were generally lower, often 5 to 20 tonnes per hectare per year, because slopes were gentler. However, drainage networks connected this mixed agricultural and expanding residential landscape to Cyohoha Lake, increasing the potential delivery of detached sediment to wetlands and a transboundary lake system. The distinction between erosion on a slope and sediment delivery downstream is important: the mapped estimate concerns soil loss across the landscape, while the lake connection identifies a pathway of potential transport rather than a measured sediment load.

The study recommends geographically targeted measures because its model identifies different dominant controls in different parts of the catchment. For steep agricultural slopes, it calls for bench terracing, grass and legume waterways and contour practices that interrupt runoff. Bare and degraded southwestern areas require rapid revegetation, including pioneer grasses, shrubs and agroforestry systems, while agricultural intervals with exposed soil could be protected with fast-establishing cover crops.

In the erodible southern uplands, the authors emphasize sustained organic-matter inputs, residue retention, reduced disturbance of soil structure and cross-slope contour bunds. Around Cyohoha Lake, riparian buffers, sediment-retention structures and controls on bare construction surfaces are proposed to reduce downstream sediment connectivity. These are conservation proposals based on the modeled vulnerability patterns; the research did not test their effectiveness in field experiments.

These findings should be read with the model’s limitations in mind. The analysis also relies on interpolated rainfall erosivity, global soil-texture information, satellite-derived vegetation cover and a constant support-practice factor; the accuracy of individual pixel estimates therefore depends on those inputs. The catchment’s different sub-catchments, including the Akanyaru Hills, Akanyaru Valley Left Bank and Akanyaru Cyohoha, also have different topographic and land-cover conditions, so the relative contribution of each factor varies across the landscape.

The 30-metre output provides a common spatial reference for identifying where modeled risk is concentrated, but it does not establish the amount of soil physically measured at a particular farm, nor does it demonstrate that land-cover change alone caused a measured increase in erosion. It offers a detailed baseline for comparing future conservation outcomes and targeting monitoring where steep slopes, exposed soil, erodible materials and intense rainfall overlap.

Subject of Research: Spatial modeling of soil loss across four land-cover types in Rwanda’s Akanyaru Catchment

Article Title: Quantification of soil loss by using RUSLE model under different land use and land cover conditions in AKANYARU catchment in Rwanda

Article References: Ntirenganya, J. B., Bugenimana, E. D., Niyonkuru, R., Iragena, A., Uwamahoro, P. C., Bosco, T. J., & Uwihirwe, J. (2026). Quantification of soil loss by using RUSLE model under different land use and land cover conditions in AKANYARU catchment in Rwanda. Discover Soil, 3(1), Article 191. https://doi.org/10.1007/s44378-026-00347-1

Image Credits: AI Generated

DOI: 10.1007/s44378-026-00347-1

Keywords: soil erosion, RUSLE, Akanyaru Catchment, land cover, rainfall erosivity, soil erodibility, agricultural slopes, GIS, Rwanda, erosion hotspots, Quantification, soil

The research team, led by Ninu Krishnan Modon Valappil of Universiti Sains Malaysia, together with Yusri Yusup and Vijith Hamza, drew on the India Meteorological Department’s high-resolution daily gridded rainfall dataset, which covers the subcontinent at a quarter-degree grid spacing and extends back to the beginning of the twentieth century. From these daily records, aggregated into monthly totals, the team computed two closely related quantities. The first is rainfall erosivity, often called the R-factor, a term in the Universal Soil Loss Equation family of models that quantifies the kinetic punch delivered by raindrops and the runoff they generate. The second is erosivity density, which normalizes that punch by the amount of rain, effectively asking how destructive each millimeter of rainfall is on average.

To estimate these quantities across the full century, the researchers employed the monthly rainfall-based empirical model introduced by Arnoldus in 1980, a widely used approach when sub-hourly rainfall intensity data are unavailable, as they are for most of the historical record. The resulting values reveal an extraordinary range. Annual rainfall across Kerala varied from as little as 134 millimeters to as much as 5,424 millimeters in individual grid cells and years. Rainfall erosivity ranged from 61 to 58,063 megajoule millimeters per hectare per hour per year, a spread of nearly three orders of magnitude, while erosivity density spanned 0.38 to 17.61 megajoules per hectare per hour. That enormous variability is precisely why the authors argue that averages and totals conceal more than they reveal about erosion risk.

Geographically, the study found a persistent north-south divide. Higher rainfall, higher erosivity, and higher erosivity density were consistently concentrated in northern Kerala, where the Western Ghats force moisture-laden monsoon winds upward and squeeze out intense orographic precipitation. Lower values predominated across much of the southern region. This spatial pattern matters because the Western Ghats are recognized as one of the world’s biodiversity hotspots, and previous work has documented substantial soil loss across the region, including dramatic erosion episodes following the severe Kerala floods of 2018. Knowing where the erosive energy of the climate is concentrated provides a scientific basis for targeting watershed management and soil conservation measures where they will do the most good.

The temporal analysis was where the study broke new ground. Using linear trend analysis alongside seasonal, decadal, and inter-decadal comparisons, the team found that rainfall, erosivity, and erosivity density did not move in lockstep. Instead, the record alternated between phases of increasing and decreasing values on decadal timescales, and seasonal hotspot analysis, performed with the Getis-Ord Gi* statistic, a method for identifying statistically significant spatial clustering, revealed pronounced shifts between monsoon and non-monsoon periods. In other words, the places and times where erosive power concentrates are not fixed features of the landscape but migrate through the decades and across the calendar, responding to the shifting rhythms of the monsoon system.

The conceptual centerpiece of the paper is the decoupling index, a measure borrowed from economics, where decoupling analysis was developed to examine whether economic growth could be separated from environmental damage. Applied here, the index asks a simple question: when rainfall amount changes, does erosivity change proportionally? The answer, across most of Kerala’s twentieth century, was no. Weak coupling predominated throughout the study period, meaning that changes in how much rain fell were only loosely reflected in changes in how erosive that rain was. More strikingly, the results suggest that climatic rainfall erosivity became progressively less dependent on rainfall amount alone as the century wore on, hinting that the character of the rain itself, its intensity, concentration, and timing, has been changing in ways that totals cannot capture.

This decoupling has a physical explanation rooted in how raindrops transfer energy to the ground. Erosivity scales with the kinetic energy of falling drops and with rainfall intensity, not merely with volume. A season that delivers the same total rainfall as another, but in fewer, fiercer bursts, will strip far more soil. Climate change is widely expected to intensify precisely this pattern across the tropics, with warming seas and atmospheres loading more moisture into individual storm events even where total rainfall stagnates or declines. Related studies cited by the authors have documented intensifying erosivity in West Africa, and research along the Western Ghats and the southwest coast of India has documented changes in extreme rainfall, mesoscale convective systems, and moisture transport in recent decades, all consistent with a monsoon regime whose extremes are sharpening.

For Kerala, the practical stakes are considerable. The state’s steep slopes, lateritic soils, dense river networks, and reservoir-dependent agriculture make it acutely sensitive to sediment loss, which chokes reservoirs, degrades farmland, and compounds landslide and flood hazards. The study’s authors frame their results as a foundation for regional soil erosion assessment, watershed management, and climate adaptation in tropical monsoon environments. If planners continue to infer erosion risk from rainfall totals, they may systematically misjudge which decades and districts face the greatest threat. A decade of modest total rainfall punctuated by violent downpours could be more erosive than a wetter, gentler decade, and the decoupling index offers a way to detect exactly that divergence in the historical record.

Methodologically, the study also demonstrates the value of squeezing more from the data that exist. True erosivity calculations ideally require high-temporal-resolution rainfall intensity measurements, which global efforts such as the Global Rainfall Erosivity Database have assembled for recent decades. But century-scale assessment demands the long observational records that only monthly or daily data can provide, and the Arnoldus monthly model, though an approximation, allows researchers to extend erosion-relevant analysis back through periods when no rain gauge recorded intensity. The trade-off is acknowledged in the literature, and the authors’ use of trend analysis, hotspot statistics, and the decoupling index together provides a more robust picture than any single metric could, triangulating on the underlying behavior of the monsoon system.

The broader message extends well beyond Kerala. Rainfall-driven soil erosion is a major cause of land degradation in tropical monsoon regions, where hundreds of millions of people depend on rain-fed agriculture. As global assessments of rainfall erosivity grow more sophisticated, the Kerala study adds a century-scale caution: the relationship between the amount of water falling from the sky and the damage that water does is neither fixed nor guaranteed. In a warming world, that relationship appears to be loosening, and the erosion threat may be growing fastest precisely where rainfall statistics look unremarkable. For the steep, green slopes of the Western Ghats and for monsoon landscapes across Asia, Africa, and South America, the rain that matters most may be the rain that falls hardest, not the rain that falls most.

Subject of Research: Century-scale changes in rainfall erosivity and its decoupling from rainfall amount in Kerala, India

Article Title: Decoupling rainfall amount and rainfall erosivity: century-scale changes in climatic rainfall erosivity across Kerala, India

Article References: Valappil, N. K. M., Yusup, Y., & Hamza, V. (2026). Decoupling rainfall amount and rainfall erosivity: century-scale changes in climatic rainfall erosivity across Kerala, India. Theoretical and Applied Climatology, 157(10), Article 675. https://doi.org/10.1007/s00704-026-06587-z

Image Credits: AI Generated

DOI: 10.1007/s00704-026-06587-z

Keywords: rainfall erosivity, soil erosion, Kerala, monsoon, Western Ghats, erosivity density, decoupling index, climate change, India Meteorological Department, watershed management, RUSLE, land degradation

The study, conducted by Prodip Mandal and Mayank Joshi of the Govind Ballabh Pant National Institute of Himalayan Environment and Siddharth Prizomwala of the Institute of Seismological Research, applied the Revised Universal Soil Loss Equation, widely known as RUSLE, to a 75.52 square kilometre area around the rapidly growing district of Namchi. Their results, published in the journal Discover Geoscience, estimate an average annual soil loss of 39.10 tonnes per hectare per year, with localized hotspots reaching up to 138.62 tonnes per hectare per year. While those figures far exceed the global average erosion rate of roughly 2.4 tonnes per hectare per year, they are modest compared with rates exceeding 800 tonnes per hectare per year recorded in parts of India’s Northeastern Himalaya, where shifting cultivation and deforested wastelands leave slopes fully exposed.

RUSLE is an empirical model originally developed by the United States Department of Agriculture that computes average annual soil loss as the product of five factors: rainfall erosivity, soil erodibility, slope length and steepness, land-cover management, and conservation practices. Its enduring appeal lies in its adaptability. Because each factor can be derived from remotely sensed data and processed within a geographic information system, the model works well in data-scarce mountain regions where rain gauges with sub-hourly intensity measurements and detailed field soil surveys simply do not exist. That is precisely the situation in much of the Sikkim Himalaya, a tectonically young fold mountain belt built of weakly consolidated rocks from the Daling Group, threaded with active faults and hammered by prolonged monsoonal rainfall that makes the region a natural hotspot for erosion and landslides alike.

To build each factor layer, the researchers assembled a mix of datasets. Rainfall erosivity, which captures the capacity of falling rain to detach and transport soil particles, was estimated from daily rainfall records at six stations in and around the study area using an empirical relationship suited to Indian conditions, yielding values between 717.47 and 826.03 megajoule millimetres per hectare per hour per year, highest in the wetter southern portion of the map. Soil erodibility was computed with the EPIC equation from sand, silt, clay and organic carbon contents drawn from the ISRIC Soil Grids database, producing values between roughly 0.0137 and 0.0156, with the more erodible soils at higher elevations where organic carbon is abundant. The topographic factor, combining slope length and steepness, was derived from a 12.5-metre ALOS PALSAR digital elevation model and reached values as high as 19.05 in dissected terrain near drainage channels.

The remaining two factors encode human influence on the land surface. The cover management factor, ranging from 0.004 to 1.0, was assigned to five land-use classes mapped from European Space Agency satellite imagery acquired in October 2024 and refined with high-resolution imagery from Google Earth: dense forest, moderately dense forest, sparse forest, agricultural land and built-up areas. Forested slopes, where canopy intercepts raindrops and roots bind the soil, received the lowest values, while bare and open land received the highest. The conservation support practice factor, spanning 0.55 to 1.0, reflects measures such as terracing, contour farming and retaining walls that locally slow runoff. When the five layers were multiplied together in raster format and the resulting erosion values classified with the Jenks natural breaks method, the picture that emerged was strikingly uneven: 89.52 percent of the study area falls in the low-risk class, 6.52 percent in the moderate class and only 3.96 percent in the high-risk class.

Where does the sediment actually come from? The answer, overwhelmingly, is the steepest farmland. Among land-use classes, agricultural land and open land showed the highest erosion rates, peaking at 138.62 tonnes per hectare per year, whereas dense forest lost only 21.53 tonnes and moderately dense forest 20.23 tonnes per hectare per year. Together, dense and moderately dense forest blanket more than 55 percent of the study area, and the analysis makes clear that this vegetation is the single most important brake on soil loss, absorbing raindrop kinetic energy and reducing runoff velocity before it can carve rills into hillsides. Interestingly, the urban core of Namchi itself recorded relatively low erosion of 85.76 tonnes per hectare per year, an unexpected finding the researchers attribute to impervious concrete surfaces and the retaining walls residents build to protect hillside properties, which shield the soil from direct raindrop impact even as construction churns the surrounding periphery.

Slope proved equally decisive. The highest erosion rates occurred on slopes between 35 and 50 degrees, where gravity accelerates runoff and shallow, weakly stable soils sit atop fragile Himalayan bedrock. Somewhat counterintuitively, the steepest slopes of all, above 50 degrees, showed the lowest measured rates at just 6.33 tonnes per hectare per year, because exposed hard rock and difficult access limit both soil availability and human interference there. The geographic pattern also follows the rain: the southern part of the study area, closer to the moisture-laden winds arriving from the Bay of Bengal, receives more rainfall and more sunshine on its south-facing slopes, and consequently loses more soil than the northern segment.

Seeking to understand what drives these patterns, the team extracted six topographic and hydrological factors from the digital elevation model, including elevation, slope, aspect, relative relief, topographic wetness index and drainage density, and examined their pairwise relationships with Pearson correlation coefficients after checking for multicollinearity. Elevation showed weak but statistically significant positive correlations with slope, relative relief and drainage density, and a weak negative correlation with the topographic wetness index, which is higher in low-lying zones where moisture accumulates and deposition rather than erosion dominates. Aspect correlated positively with slope and relief, and the strongest relationship in the matrix linked rainfall erosivity with drainage density, underscoring how climatic forcing and landscape dissection interact to concentrate erosion along densely channelled terrain. Overall, relative relief and slope emerged as the most influential contributors to soil loss in the region.

The implications reach well beyond geomorphology. The authors frame their findings against the United Nations Sustainable Development Goals, noting that erosion degrades spring catchments that mountain communities depend on for drinking water, undermining SDG 6 on clean water; strips fertile topsoil from farms, deepening rural inequality under SDG 10; and threatens the sustainability of fast-growing Himalayan towns under SDG 11. Because erosion here is most severe during the monsoon, targeted interventions such as contour farming, terracing and slope protection on the moderate-to-steep agricultural slopes of the southern study area could deliver outsized benefits. Prior studies cited in the work suggest that well-implemented conservation practices can cut erosion by a factor of 2.5 to 7, and the Namchi results hint that even informal measures by local residents already mute erosion around settlements.

The study is candid about its limits. Empirical erosivity equations built on daily rainfall cannot fully capture the short, violent cloudbursts that characterize early-monsoon events in the Eastern Himalaya and frequently trigger landslides, and globally generalized soil databases may smooth over local heterogeneity in the rugged terrain. Climate change, intensifying extreme precipitation and accelerating infrastructure development are all expected to shift erosion rates in ways the current model does not project. Even so, the team argues that RUSLE, paired with freely available satellite data, offers a practical blueprint for other data-scarce Himalayan districts: by identifying exactly which slopes, land uses and aspect zones bleed soil fastest, it gives planners a prioritized map for conservation investment in one of the most erosion-prone mountain belts on the planet.

Subject of Research: Quantitative assessment of monsoon-driven soil erosion using the RUSLE model in the Sikkim Himalaya

Article Title: RUSLE-based quantitative soil erosion assessment in the monsoon-dominated region of Sikkim Himalaya

Article References: Mandal, P., Joshi, M., & Prizomwala, S. (2026). RUSLE-based quantitative soil erosion assessment in the monsoon-dominated region of Sikkim Himalaya. Discover Geoscience, 4(1), Article 348. https://doi.org/10.1007/s44288-026-00706-z

Image Credits: AI Generated

DOI: 10.1007/s44288-026-00706-z

Keywords: soil erosion, RUSLE, Sikkim Himalaya, monsoon rainfall, land use, slope classes, soil conservation, GIS and remote sensing, Namchi, rainfall erosivity, forest cover, sustainable development goals