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, conducted by soil scientists from ICAR-National Bureau of Soil Survey and Land Use Planning, examined soils developed on the Archaean rocks of the Shillong Group, a terrain uplifted rapidly during the late Cenozoic roughly 8 to 14 million years ago. Working between Mylliem and Pynursla at elevations above 1200 metres, the team described 55 soil profiles along eleven transects and distilled them into twelve representative soil series. These were grouped by land use: pine forest on summits and side slopes, potato cultivation on the northern slopes, deep grassland soils in the south, and seasonally waterlogged rice paddies in valley bottoms. The soils belong mainly to the orders Ultisols, Inceptisols and Entisols, products of intense weathering under a warm, perhumid monsoonal climate.

The climatic backdrop is extreme. The southern escarpment near Cherrapunji receives an average annual rainfall exceeding 12,000 millimetres, and around 80 percent of the plateau’s precipitation falls during the June-to-September monsoon. Against this assault, the practice of shifting cultivation, or jhum, has proven devastating. Previous studies cited in the paper estimate soil losses of about 40 megagrams per hectare per year on slopes of 21 to 26 degrees under shifting cultivation, rising to a staggering 145 to 170 megagrams per hectare per year on steeper slopes of 31 to 37 degrees during the initial cropping years. By contrast, fallow vegetation cuts losses to about 7 megagrams per hectare per year after one year and 2 after ten, while grass-covered slopes lose only around 2.1 megagrams per hectare per year.

The core of the new work lies in a set of erodibility indices that translate soil chemistry and physics into erosion risk. The researchers calculated the critical shear stress, which measures a soil particle’s resistance to detachment by flowing water; the erosion index of organic matter, known as EIROM; the critical limit of organic matter, or CLOM, the minimum organic matter needed to keep soil structure intact; and the soil structural stability index, SSI. They then estimated the soil erodibility factor K, the key parameter of the Universal Soil Loss Equation, using multiple linear regression models built from these indices. The average K value across the plateau was 0.04 plus or minus 0.02 tonnes hectare hour per hectare megajoule per millimetre, with a coefficient of variation of 45.17 percent, signalling considerable spatial variability in erosion processes across the landscape.

The statistical results were striking. A regression model combining CLOM, critical shear stress and very coarse sand content predicted the erodibility factor with a coefficient of determination of 0.97, an exceptionally tight fit for field soil data. A second model expressed structural stability as a function of EIROM, shear stress, the cation-exchange-capacity-to-clay ratio and K, explaining 78 percent of the variance. Critical shear stress emerged as the dominant control, correlating with the erodibility factor at r = 0.96. In plain terms, the soils’ resistance to being scoured by runoff, together with their organic matter reserves and coarse sand content, almost entirely determines how readily they wash away.

Land use proved decisive. Soils under grass cover were found to be only slightly to moderately erodible, whereas soils under pine plantations and potato cultivation showed markedly higher erosion susceptibility. Potato-growing soils on side slopes are generally shallow to moderately shallow, and intensive cultivation there accelerates the decomposition of organic matter, weakening the aggregates that bind the soil. The structural stability index told a consistent story: surface horizons of some potato soils were relatively stable, but subsurface horizons of potato and grassland soils fell below the critical threshold of 5, indicating structural degradation and heightened vulnerability. Exposed clay-enriched Bt horizons showed only moderate structural stability, with a mean SSI of 5.49, and greater susceptibility to erosion in grassland and potato-growing areas.

The chemistry of these soils compounds their fragility. Nearly all profiles were strongly to extremely acidic, with pH values in surface horizons as low as 3.6 under pine, and exchangeable aluminium dominated the exchange complex while base cations generally remained below 1 centimole per kilogram. Cation exchange capacity was mostly very low, below 15 centimoles per kilogram, reflecting the dominance of low-activity silicate clays formed under relentless leaching. Organic carbon correlated almost perfectly with CEC at r = 0.98, underscoring how much these acidic soils depend on their thin organic blankets for fertility and cohesion. Notably, the erodibility factor showed a significant negative relationship with the exchangeable calcium-to-aluminium ratio, suggesting that calcium helps stabilise aggregates while aluminium dominance under acidity weakens them.

Texture added a geomorphic twist. Elevation correlated negatively with fine sand and total sand but positively with silt, meaning higher ground carries finer-textured soils. Many grassland profiles showed abrupt lithological discontinuities, evidence of buried horizons and episodic colluvial deposition from landslides and slope wash. The researchers attribute the heterogeneous textures, roughly 50 percent silt, 30 percent sand and 20 percent clay in the dominant loamy matrix, to repeated mass-movement events rather than purely in-place soil formation, with coarse sand fractions possibly produced by mechanical crushing of rock fragments at landslide bases.

K-means clustering of eleven erodibility parameters sorted the soils into three groups that mirrored these patterns. The first cluster, dominated by pine-forest horizons, showed high variability in structural stability and organic matter; the second, mostly grassland and potato horizons, displayed moderate structural stability with a mean SSI of 5.49; the third, drawn mainly from the southern plateau, revealed structural degradation, organic matter depletion and the highest variability of all. Across clusters, the depletion of organic carbon and sand fractions in profiles with higher erodibility indices confirmed that land management leaves a measurable fingerprint in the soil’s erosion resistance.

The study’s authors argue that detailed pedological datasets at the soil-series level, such as the benchmark sites established under India’s National Agricultural Technology Project, are essential for monitoring erodibility across agro-ecological regions and for calibrating conservation strategies. With an estimated 5.3 billion tonnes of soil lost to erosion in India each year and nearly two-thirds of Meghalaya’s population dependent on agriculture, the message is urgent: protecting the organic-rich surface layers, maintaining vegetative cover and tailoring interventions to the specific soil series could determine whether the abode of clouds keeps its soil, or watches it slide toward the Bay of Bengal.

Subject of Research: Pedological controls on soil erodibility in high-altitude Ultisols of the Meghalaya Plateau, India

Article Title: Pedological controls on soil erodibility in high-altitude ultisols and associated soils of the Meghalaya plateau, India

Article References: Bhaskar, B. P., Sharma, R. P., Baruah, U., & Butte, P. S. (2026). Pedological controls on soil erodibility in high-altitude ultisols and associated soils of the Meghalaya plateau, India. Discover Soil, 3(1), Article 113. https://doi.org/10.1007/s44378-026-00274-1

Image Credits: AI Generated

DOI: 10.1007/s44378-026-00274-1

Keywords: soil erosion, soil erodibility, Meghalaya Plateau, Ultisols, shifting cultivation, critical shear stress, soil organic carbon, USLE K-factor, soil structural stability, land use, Northeast India, soil conservation

The research focused on completely decomposed granite, a sandy, weakly cohesive soil that is ubiquitous in man-made slopes across Hong Kong and much of subtropical Asia. The team sieved the soil into two fractions, one with a maximum particle size of 5 millimeters and one with a maximum of 2 millimeters, and compacted each to relative densities of 90, 80, and 70 percent. That matrix of six bare-soil treatments formed the baseline. For the biological treatments, the researchers used an elegant compartmental design: pots were split by a 38-micrometer mesh, with bahiia grass (Paspalum notatum) grown on one side and only mycorrhizal fungal hyphae, root exudates, and other dissolved biological agents allowed to cross into the neighboring compartment. This meant one side of each pot contained living roots, while the other received the biochemical products of root activity without any physical roots at all, allowing the team to separate the two great reinforcement pathways that roots provide.

Those two pathways are physical binding and biochemical bonding. Physical binding is the intuitive one: fine roots physically enmesh soil aggregates, the way rebar reinforces concrete. Biochemical bonding is subtler and, until now, much harder to quantify. Root exudates, decomposing organic matter, glomalin-related proteins deposited by mycorrhizal fungi, and microbial polysaccharides all act as glues at the microaggregate scale, stabilizing soil structure without any mechanical entanglement. Earlier studies had estimated that physical binding accounted for roughly 70 to 76 percent of the erosion resistance conferred by herbaceous roots, but those estimates typically relied on short-term growth experiments or dead roots, in which biochemical effects could not be fully eliminated. By constructing bonded soils directly, the Hong Kong team could measure each contribution cleanly for the first time in compacted, engineered slope materials.

The erosion measurements themselves were carried out in a custom-built transparent flume nearly four meters long, in which saturated soil cores were scoured under carefully controlled combinations of flow discharge and bed gradient. Nine to eleven shear stress levels were tested for each treatment, spanning the transition from gentle sheet flow to the concentrated rill flows that carve gullies into bare slopes during intense storms. Detachment rates were then fitted to two competing mathematical models. The first, the classic excess shear stress model, assumes a simple linear relationship between detachment rate and the stress exceeding a critical threshold. It is convenient but, as the authors note, its erodibility coefficient lacks any real physical meaning. The second, Wilson’s model, derives from a two-dimensional representation of individual particles interacting with turbulent flow, and its two parameters, b0 and b1, correspond to the rate at which particles exchange with the flow and the probability that drag forces exceed the stabilizing forces holding them in place.

Across the bare soils, the results confirmed textbook expectations with striking clarity. The coarser 5-millimeter soils were consistently harder to erode than the finer 2-millimeter soils, showing lower erodibility coefficients and higher critical shear stresses, with mean critical shear stress of 2.66 pascals versus just 0.94 pascals for the finer material. Compacting the soil to 90 percent relative density amplified this resistance further, though the difference between 80 and 70 percent densities was modest, suggesting that compaction below 80 percent buys engineers surprisingly little erosion protection. The Wilson’s model generally fit the data better than the excess shear stress model, particularly for the looser soils, reinforcing the argument that a physically based framework is worth its extra mathematical complexity when the goal is mechanistic understanding rather than mere curve fitting.

The biological results were more surprising. Roots and bonding agents consistently reduced erodibility coefficients and raised shear thresholds across every particle size and density combination, but the magnitude of the effect depended strongly on soil conditions. The largest benefits appeared in the finer, moderately dense soils: in the 2-millimeter soil compacted to 80 percent relative density, the presence of roots or bonding agents cut the erodibility coefficient by more than fourfold and the Wilson parameter b0 by nearly sixfold. Paradoxically, the densest compaction, which engineers specify precisely to stabilize slopes, worked against biology. At 90 percent relative density, grass roots grew poorly, producing the smallest shoot biomass, root length density, and root surface area of any treatment, and the resulting reinforcement was correspondingly weak. The soil’s fine fraction, it turns out, is the essential substrate on which roots and their biochemical products build cohesion; without it, even abundant roots have little to grip.

The centerpiece of the study is a new dimensionless indicator, Kbio, derived by algebraically manipulating Wilson’s model. The researchers showed that the combined term sqrt(b0 squared times b1) behaves as a soil-specific erodibility constant that is largely independent of the flow conditions during detachment, and that it is governed by the ratio of particle resistance to the exposure resistance of the surrounding finer particles. Because both roots and bonding agents act primarily by adding cohesion, this term shrinks when biology strengthens the soil, and the ratio of its value before and after biological treatment defines Kbio. A value greater than one signals reduced erodibility. When the team computed Kbio for each treatment, the total biological effect ranged from 1.60 in the densest coarse soil to 4.34 in the finer moderately dense soil, while the physical binding component stayed relatively constant, between 1.11 and 1.63, across all conditions. The biochemical bonding component, by contrast, varied widely and dominated the reinforcement in the finer, moderately dense soils.

This finding upends the conventional wisdom that physical binding is the dominant mechanism. In compacted sandy slope soils, the invisible chemistry of root exudates and fungal hyphae appears to do most of the heavy lifting, at least where fine particles are available to be glued together. The authors validated their modified Wilson’s framework against twenty independent datasets, including laboratory flume tests on peanut-farming soils under different management practices and a field test on an unpaved loess road, and the reparameterized model reproduced the original model’s predictions in eighteen of the twenty cases. They are candid about the limits: the indicator has so far been tested only on sandy soils with low cohesion and high permeability, and extreme turbulence conditions, where instantaneous shear stresses spike far above the mean, remain outside the model’s validated domain.

The practical implications reach well beyond the laboratory. Ecological engineering, the use of living plants as construction materials, is increasingly specified for roadside slopes, dikes, and restored terrain worldwide, yet designers currently rely on empirical root indices that are soil-specific and confounded by multiple factors. A dimensionless, physically meaningful indicator like Kbio could let engineers predict, before planting, how much a given grass species will actually reduce erosion on a given compacted soil, and could guide the choice of both plant species and compaction specifications. The study also sends a caution to geotechnical practice: specifying maximum compaction for slope stability may inadvertently suppress the biological reinforcement that protects the same slope from washing away. The authors call for field validation across different rainfall regimes and soil textures, but the conceptual advance is clear. Erosion resistance, long treated as an empirical black box, can now be decomposed into measurable physical and biological parts, and the invisible half of the story, the chemistry flowing from roots into soil, turns out to be the bigger half after all.

Subject of Research: Effects of plant roots and biochemical bonding on the erodibility of compacted sandy soils under rill flow

Article Title: Effects of plant roots on the erodibility of compacted soils can be explained by a biological indicator derived from Wilson’s model

Article References: Li, M., Leung, A. K., & Chen, X. W. (2026). Effects of plant roots on the erodibility of compacted soils can be explained by a biological indicator derived from Wilson’s model. Plant and Soil. https://doi.org/10.1007/s11104-026-09147-6

Image Credits: AI Generated

DOI: 10.1007/s11104-026-09147-6

Keywords: soil erosion, plant roots, soil erodibility, Wilson's model, rill detachment, compacted soils, biochemical bonding, physical binding, mycorrhizal fungi, aggregate stability, flume tests, ecological engineering