Insects as zoogeomorphic agents: an extended review

Insects are the largest and most diverse group of living organisms on Earth, playing a critical but underestimated role as agents of geomorphic change. Burrowing insects create micro-scale landforms such as subterranean tunnels and surface mounds and, by this way, exert an influence on hydrology, soil erosion and sediment transfer at a wider landscape scale. However, social insects represented by ants and termites were the main taxa studied as geomorphic agents and ecosystem engineers. This article proposes an extended and critical literature review of insects as zoogeomorphic agents, with reference to various taxonomic orders and families of insects having a burrowing behaviour. It provides a large overview of their primary and secondary impacts on Earth surface systems, both supported by naturalistic evidence and available quantitative data. Some evolutionary insights are discussed based on fossil evidence of geomorphic work by insects and, at finer temporal scale, on recent advances in radiometric and luminescence dating of insect mounds. Finally, this article explores the fruitful links between geomorphology and entomology, and suggests several research perspectives in order to develop an integrated understanding of the importance of insects in Earth surface processes and landforms. © 2020 John Wiley & Sons, Ltd.     

Effect of subsurface water level on gully headcut retreat in tropical highlands of Ethiopia

Gully erosion is a major cause of soil loss and severe land degradation in sub-humid Ethiopia. The objective of this study was to investigate the role and the effect of subsurface water level change on gully headcut retreat, gully formation and expansion in high rainfall tropical regions in the Ethiopian highlands. During the rainy seasons of 2017–2019, the expansion rate of 16 fixed gullies was measured and subsurface water levels were measured by piezometers installed near gully heads. During the study period, headcut retreats ranged from 0.70 to 2.35 m, with a mean value of 1.49 ± 0.56 m year⁻¹, and average depth of the surface water level varied between 1.12 and 2.82 m, with a mean value of 2.62 m. Gully cross-section areas ranged from 2.90 to 20.90 m², with an average of 9.31 ± 4.80 m². Volumetric retreat of gully headcuts ranged from 4.49 to 40.55 m³ and averaged 13.34 ± 9.10 m³. Soil loss from individual gullies ranged from 5.79 to 52.31 t year⁻¹ and averaged 17.21 ± 11.74 t year⁻¹. The headcut retreat rate and sediment yield were closely related over the three study seasons. Elevated subsurface water levels facilitated the slumping of gully banks and heads, causing high sediment yield. When the soil was saturated, bank collapse and headcut retreat were favoured by the combination of elevated subsurface water and high rainfall. This study indicates that area exclosures are effective in controlling subsurface water level, thus reducing gully headcut retreat and associated soil loss.     

Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

Convergent development of experimental landforms with differing rainfall and substrate properties

Four runs of experimental landform development, with the same uplift rate, different rainfall intensity, and the same material of different permeability adjusted by the degree of compaction, showed complicated effects of rainfall and mound-forming material. In the run with more rainfall on less permeable material, low separated ridges developed in the uplifted area, because abundant overland flow promoted valley erosion and slope processes from early stages. In the run with less rainfall on less permeable material, valley incision proceeded mostly in major valleys where surface water converges. Canyons developed during early stages and later a high massive mountain emerged. The effect of rainfall difference, however, appeared completely opposite on more permeable material accompanied by lower shear strength. In the run with more rainfall on more permeable material, a massive mountain similar to that with less rainfall on less permeable material appeared, and low separated ridges appeared in the run with less rainfall on more permeable material as in the run with more rainfall on less permeable material. In the former case, similar amount of water available for Hortonian overland flow in early stages estimated from rainfall rate and permeability can explain the development of similar landforms. In the latter case, while abundant surface water with more rainfall on less permeable material made fluvial erosion active from early stages, the deficiency in surface water with less rainfall on more permeable material apparently attenuated fluvial erosion but possibly accentuated slope processes and slope failures by seepage water flow through more permeable material of low shear strength. The active erosion from early stages apparently resulted in the development of enduring similar low landforms later in the dynamic equilibrium stage. These experimental results indicate that similar landforms can emerge from different environmental and lithologic controls, and that process does not necessarily follow from form.     

Influence of thawed soil depth on rainfall erosion of frozen bare meadow soil in the Qinghai–Tibet Plateau

The Qinghai–Tibet Plateau has a vast area of approximately 70×10⁴ km² of alpine meadow under the impacts of soil freezing and thawing, thereby inducing intensive water erosion. Quantifying the rainfall erosion process of partially thawed soil provides the basis for model simulation of soil erosion on cold-region hillslopes. In this study, we conducted a laboratory experiment on rainfall-induced erosion of partially thawed soil slope under four slope gradients (5, 10, 15, and 20°), three rainfall intensities (30, 60, and 90 mm h⁻¹), and three thawed soil depths (1, 2, and 10 cm). The results indicated that shallow thawed soil depth aggravated soil erosion of partially thawed soil slopes under low hydrodynamic conditions (rainfall intensity of 30 mm h⁻¹ and slope gradient ≤ 15°), whereas it inhibited erosion under high hydrodynamic conditions (rainfall intensity ≥ 60 mm h⁻¹ or slope gradient > 15°). Soil erosion was controlled by the thawed soil depth and runoff hydrodynamic conditions. When the sediment supply was sufficient, the shallow thawed soil depth had a higher erosion potential and a larger sediment concentration. On the contrary, when the sediment supply was insufficient, the shallow thawed soil depth resulted in lower sediment erosion and a smaller sediment concentration. The hydrodynamic runoff conditions determined whether the sediment supply was sufficient. We propose a model to predict sediment delivery under different slope gradients, rainfall intensities, and thawed soil depths. The model, with a Nash–Sutcliffe efficiency of 0.95, accurately predicted the sediment delivery under different conditions, which was helpful for quantification of the complex feedback of sediment delivery to the factors influencing rainfall erosion of partially thawed soil. This study provides valuable insights into the rainfall erosion mechanism of partially thawed soil slopes in the Qinghai–Tibet Plateau and provides a basis for further studies on soil erosion under different hydrodynamic conditions.     

Soil erosion and sediment transport in Tanzania: Part I – sediment source tracing in three neighbouring river catchments

Water bodies in Tanzania are experiencing increased siltation, which is threatening water quality, ecosystem health, and livelihood security in the region. This phenomenon is caused by increasing rates of upstream soil erosion and downstream sediment transport. However, a lack of knowledge on the contributions from different catchment zones, land-use types, and dominant erosion processes, to the transported sediment is undermining the mitigation of soil degradation at the source of the problem. In this context, complementary sediment source tracing techniques were applied in three Tanzanian river systems to further the understanding of the complex dynamics of soil erosion and sediment transport in the region. Analysis of the geochemical and biochemical fingerprints revealed a highly complex and variable soil system that could be grouped in distinct classes. These soil classes were unmixed against riverine sediment fingerprints using the Bayesian MixSIAR model, yielding proportionate source contributions for each catchment. This sediment source tracing indicated that hillslope erosion on the open rangelands and maize croplands in the mid-zone contributed over 75% of the transported sediment load in all three river systems during the sampling time-period. By integrating geochemical and biochemical fingerprints in sediment source tracing techniques, this study demonstrated links between land use, soil erosion and downstream sediment transport in Tanzania. This evidence can guide land managers in designing targeted interventions that safeguard both soil health and water quality.     

A modified Raupach's model applicable for shear-stress partitioning on surfaces covered with dense and flat-shaped gravel roughness elements

A commonly used measure to prevent soil wind erosion is to cover the surface with gravel. Gravel can inhibit soil erosion by covering the surface directly, changing the airflow field near the surface and sharing the shear stress of wind. Similar to other roughness elements, the protective effect of gravel on soil is usually expressed in terms of the ratio of the shear stress on the exposed soil surface to the total shear stress on the rough surface due to wind, i.e. through a shear-stress partitioning model. However, the existing shear-stress partitioning models, represented by Raupach's model (RM93), are only applicable when the lateral coverage of the roughness elements, λ < 0.10, and the applicability of the models to flat-shaped roughness elements is unclear. The purpose of this study is to verify the applicability of RM93 for dense and flat-shaped gravel roughness elements by using shear-stress data from wind-tunnel measurements pertaining to roughness elements with different densities (0.013 ≤ λ ≤ 0.318) and flat shapes (height-to-width ratios in the range 0.20 ≤ H/W ≤ 0.63), and to modify RM93 to enhance its predictive ability. The results indicate that RM93 cannot accurately predict the shear-stress partitioning for surfaces covered by densely distributed and flat-shaped gravel roughness elements. This phenomenon occurs because, when roughness elements are distributed densely or are flat-shaped, the proportion of the shear stress on the top surface of the roughness elements (τ_c) to the total shear stress (τ) is large; in this case, τ_c plays a dominant role and serves as an essential component in the shear-stress partitioning model. Consequently, RM93 is modified by incorporating τ_c into the calculation of τ. Under conditions of λ < 0.32 and H/W > 0.2, the modified RM93 can yield satisfactory predictions regarding the shear-stress partitioning.     

Signal crayfish burrowing,bank retreat and sediment supply to rivers: A biophysical sediment budget

Burrowing into riverbanks by animals transfers sediment directly into river channels and has been hypothesised to accelerate bank erosion and promote mass failure. A field monitoring study on two UK rivers invaded by signal crayfish (Pacifastacus leniusculus) assessed the impact of burrowing on bank erosion processes. Erosion pins were installed in 17 riverbanks across a gradient of crayfish burrow densities and monitored for 22 months. Bank retreat increased significantly with crayfish burrow density. At the bank scale (<6 m river length), high crayfish burrow densities were associated with accelerated bank retreat of up to 253% and more than a doubling of the area of bank collapse compared with banks without burrows. Direct sediment supply by burrowing activity contributed 0.2% and 0.6% of total sediment at the reach (1.1 km) and local bank (<6 m) scales. However, accelerated bank retreat caused by burrows contributed 12.2% and 29.8% of the total sediment supply at the reach and bank scales. Together, burrowing and the associated acceleration of retreat and collapse supplied an additional 25.4 t km⁻¹ a⁻¹ of floodplain sediments at one site, demonstrating the substantial impact that signal crayfish can have on fine sediment supply. For the first time, an empirical relation linking animal burrow characteristics to riverbank retreat is presented. The study adds to a small number of sediment budget studies that compare sediment fluxes driven by biotic and abiotic energy but is unique in isolating and measuring the substantial interactive effect of the acceleration of abiotic bank erosion facilitated by biotic activity. Biotic energy expended through burrowing represents an energy surcharge to the river system that can augment sediment erosion by geophysical mechanisms.     

Hillslope erosion in a grassland environment: Calibration and evaluation of the SIBERIA landscape evolution model

Field measurement and modelling of soil erosion provides insights into landscape systems as well as the potential for enhanced landscape management. There are a number of field and numerical methods by which soil erosion and deposition can be quantified. Here we examine the capability of the SIBERIA landscape evolution model to quantify short-term erosion and deposition on a well-managed cattle grazing landscape on the east coast of Australia. The model is calibrated by two methods (1) a geomorphological approach using a site digital elevation model (DEM) and soil data and (2) a laboratory-scale flume. The two calibration processes resulted in similar model input parameters and estimated erosion rates of 3.1 t ha⁻¹ year⁻¹ and 4.4 t ha⁻¹ year⁻¹, respectively. These were found to closely match erosion rates estimated using the environmental tracer ¹³⁷Cs (2.7–4.8 t ha⁻¹ year⁻¹). However, erosion and deposition estimated at individual points along the hillslope was not well correlated with ¹³⁷Cs at the same position due to the temporal averaging of the model and microtopography. Sensitivity analysis showed the model was more sensitive to parameterisation than sub-DEM-scale topography. This places confidence in the model's ability to estimate erosion and deposition across an entire hillslope and catchment on decadal time scales. We also highlight the robustness and flexibility of the calibration methods.