The effect of truck traffic and road water content on sediment delivery from unpaved forest roads

A study investigated the effect of truck‐traffic intensity and road water‐content on the quality of runoff water from unsealed forest roads. Three sections of a gravel‐surfaced forest road were instrumented and exposed to low and high levels of truck traffic during wet winter conditions and dry summer conditions between July 2001 and December 2002. Rainfall, runoff, road moisture, and traffic were measured continuously, and suspended and bedload sediments were integrated measurements over 2‐week site‐service intervals. The median suspended sediment concentration from the three road segments under low truck‐traffic conditions (less than nine return truck passes prior to a storm) was 269 mg l^?1, increasing 2·7‐fold to a median of 725 mg l^?1 under high truck‐traffic conditions (greater than or equal to nine return truck passes prior to a storm). These concentrations, and increases due to traffic, are substantially less than most previously reported values. When these data are expressed as modified universal soil loss equation (MUSLE) erodibility values K, accounting for differences in rainfall energy, site characteristics and runoff, high traffic resulted in a road surface that was four times more erodible than the same road under low traffic conditions. Using multiple regression, traffic explained 36% of the variation in MUSLE erodibility, whereas road water content was not significant in the model. There was little difference in the erodibility of the road when trafficked in low water‐content compared with high water‐content conditions (MUSLE K values of 0·0084 versus 0·0080 respectively). This study shows that, for a good quality well‐maintained gravel forest road, the level of truck traffic affects the sediment concentration of water discharging from the road, whereas the water content of the road at the time of that traffic does not (note that traffic is not allowed during runoff events in Victoria). These conclusions are conditional upon the road being adequately maintained so that trafficking does not compromise the lateral drainage of the road profile. Copyright © 2005 John Wiley & Sons, Ltd.     

Landslide scar/soil erodibility mapping using Landsat TM/ETM+ bands 7 and 3 Normalised Difference Index: A case study of central region of Kenya

Landsat series multispectral remote sensing imagery has gained increasing attention in providing solutions to environmental problems such as land degradation which exacerbate soil erosion and landslide disasters in the case of rainfall events. Multispectral data has facilitated the mapping of soils, land-cover and structural geology, all of which are factors affecting landslide occurrence. The main aim of this research was to develop a methodology to visualize and map past landslides as well as identify land degradation effects through soil erosion and land-use using remote sensing techniques in the central region of Kenya. The study area has rugged terrain and rainfall has been the main source of landslide trigger. The methodology comprised visualizing landslide scars using a False Colour Composite (FCC) and mapping soil erodibility using FCC components applying expert based classification. The components of the FCC were: the first independent component (IC1), Principal Component (PC) with most geological information, and a Normalised Difference Index (NDI) involving Landsat TM/ETM+ band 7 and 3.The FCC components formed the inputs for knowledge-based classification with the following 13 classes: runoff, extreme erosions, other erosions, landslide areas, highly erodible, stable, exposed volcanic rocks, agriculture, green forest, new forest regrowth areas, clear, turbid and salty water. Validation of the mapped landslide areas with field GPS locations of landslide affected areas showed that 66% of the points coincided well with landslide areas mapped in the year 2000. The classification maps showed landslide areas on the steep ridge faces, other erosions in agricultural areas, highly erodible zones being already weathered rocks, while runoff were mainly fluvial deposits. Thus, landuse and rainfall processes play a major role in inducing landslides in the study area.     

Complexities of landscape evolution during incision through layered stratigraphy with contrasts in rock strength

Variation in the erodibility of rock units has long been recognized as an important determinant of landscape evolution but has been little studied in landscape evolution models. We use a modified version of the Channel‐Hillslope Integrated Landscape Development (CHILD) model, which explicitly allows for variations in rock strength, to reveal and explore the remarkably rich, complex behavior induced by rock erodibility variations in even very simple geologic settings with invariant climate and tectonics. We study the importance of relative contrasts in erodibility between just two units, the order of these units (whether hard rocks overlie soft or soft rocks overlie hard) and the orientation of the contact between the two units. We emphasize the spatial and temporal evolution of erosion rates, which have important implications for basin analysis, detrital mineral records, and the interpretation of cosmogenic isotope concentrations in detrital samples. Results of the landscape evolution modeling indicate that the stratigraphic order of units in terms of erodibility, the gross orientation of the contact (i.e. dipping away or toward the outlet of the landscape) and the contact dip angle all have measurable effects on landscape evolution, including significant spatial and temporal variations in erosion rates. Steady‐state denudation conditions are unlikely to develop in landscapes with significant contrasts in rock strength in horizontal to moderately tilted rock layers, at least at the scale of the entire landscape. Additionally, our results demonstrate that there is no general relation between rock erodibility and erosion rates in natural settings. Although rock erodibility directly controls the erosion rate constant in our models, it is not uncommon for higher erosion rates to occur in the harder, less erodible rock. Indeed erosion rates may be either greater or less than the rock uplift rate (invariant in time and space in our models) in both hard and soft rocks, depending on the local geology, topography, and the pattern of landscape evolution. Copyright © 2016 John Wiley & Sons, Ltd.     

Laboratory soil piping and internal erosion experiments: evaluation of a soil piping model for low‐compacted soils

Mechanistic models have been proposed for soil piping and internal erosion on well‐compacted levees and dams, but limited research has evaluated these models in less compacted (more erodible) soils typical of hillslopes and streambanks. This study utilized a soil box (50 cm long, 50 cm wide and 20 cm tall) to conduct constant‐head, soil pipe and internal erosion experiments for two soils (clay loam from Dry Creek and sandy loam from Cow Creek streambanks) packed at uniform bulk densities. Initial gravimetric moisture contents prior to packing were 10, 12 and 14% for Dry Creek soil and 8, 12, and 14% for Cow Creek soil. A 1‐cm diameter rod was placed horizontally along the length of the soil bed during packing and carefully removed after packing to create a continuous soil pipe. A constant head was maintained at the inflow end. Flow rates and sediment concentrations were measured from the pipe outlet. Replicate submerged jet erosion tests (JETs) were conducted to derive erodibility parameters for repacked samples at the same moisture contents. Flow rates from the box experiments were used to calibrate the mechanistic model. The influence of the initial moisture content was apparent, with some pipes (8% moisture content) expanding so fast that limited data was collected. The mechanistic model was able to estimate equivalent flow rates to those observed in the experiments, but had difficulty matching observed sediment concentrations when the pipes rapidly expanded. The JETs predicted similar erodibility coefficients compared to the mechanistic model for the more erodible cases but not for the less erodible cases (14% moisture content). Improved models are needed that better define the changing soil pipe cross‐section during supply‐ and transport‐limited internal erosion, especially for piping through lower compacted (more erodible) soils as opposed to more well‐compacted soils resulting from constructing levees and dams. Copyright © 2013 John Wiley & Sons, Ltd.     

Soil erosion resistance effects on rill and gully initiation points and dimensions

Recent research has indicated the large spatial and temporal variation in soil erosion resistance against concentrated flow (SER). This study analyzes this variability in relation to rill and gully initiation locations on slopes and the downslope eroded volumes. The soil erodibility (Kc) and critical flow shear stress (τ_cr), were estimated from topsoil properties and correlated to eroded rill and gully volumes and their initiation points on slopes in the Belgian loess belt. Therefore, concentrated flow paths and topsoil properties were measured in their vicinity. The results show that rill and gully initiation points, and hence the lengths of concentrated flow paths, depend on τ_cr, which is controlled by soil surface conditions and can be predicted from saturated soil shear strength. Soil erosion control measures that increase soil shear strength (e.g. thalweg compaction), can therefore decrease rill and gully lengths. Once a rill or an ephemeral gully is initiated, its cross‐section was found to depend on Kc, which can be estimated from the soil water content, dry bulk density, and the dry density of roots and crop residues incorporated in the topsoil. 74% of the variation in the channel cross‐sectional area measured in the study area could be predicted from the combined effect of flow intensity and these three soil properties, whereas flow intensity alone could only account for 31% of the variation. Soil conservation measures affecting one of the soil properties that control Kc (e.g. double drilling of the thalweg, conservation tillage) can therefore decrease the cross‐sections of the concentrated flow paths. These findings also indicate that rill and gully initiation points are not only topographically controlled but also depend on the SER, which in turn determines the dimensions of these concentrated flow paths. Hence, knowledge of the variability in SER is indispensable. Copyright © 2009 John Wiley & Sons, Ltd.     

Rangeland hydrology and erosion model (RHEM) enhancements for applications on disturbed rangelands

The rangeland hydrology and erosion model (RHEM) is a new process‐based model developed by the USDA Agricultural Research Service. RHEM was initially developed for functionally intact rangelands where concentrated flow erosion is minimal and most soil loss occurs by rain splash and sheet flow erosion processes. Disturbance such as fire or woody plant encroachment can amplify overland flow erosion by increasing the likelihood of concentrated flow formation. In this study, we enhanced RHEM applications on disturbed rangelands by using a new approach for the prediction and parameterization of concentrated flow erosion. The new approach was conceptualized based on observations and results of experimental studies on rangelands disturbed by fire and/or by tree encroachment. The sediment detachment rate for concentrated flow was calculated using soil erodibility and hydraulic (flow width and stream power) parameters. Concentrated flow width was calculated based on flow discharge and slope using an equation developed specifically for disturbed rangelands. Soil detachment was assumed to begin with concentrated flow initiation. A dynamic erodibility concept was applied where concentrated flow erodibility was set to decrease exponentially during a run‐off event because of declining sediment availability. Erodibility was estimated using an empirical parameterization equation as a function of vegetation cover and surface soil texture. A dynamic partial differential sediment continuity equation was used to model the total detachment rate of concentrated flow and rain splash and sheet flow. The enhanced version of the model was evaluated against rainfall simulation data for three different sites that exhibit some degree of disturbance by fire and/or by tree encroachment. The coefficient of determination (R²) and Nash–Sutcliffe efficiency were 0.78 and 0.71, respectively, which indicates the capability of the model using the new approach for predicting soil loss on disturbed rangeland. By using the new concentrated flow modelling approach, the model was enhanced to be a practical tool that utilizes readily available vegetation and soil data for quantifying erosion and assessing erosion risk following rangeland disturbance. Copyright © 2014 John Wiley & Sons, Ltd.     

Bank erosion of legacy sediment at the transition from vertical to lateral stream incision

Streambank erosion is a primary source of suspended sediments in many waterways of the US Atlantic Piedmont. This problem is exacerbated where banks are comprised of fine sediment produced by the intensive land use practices of early European settlers. A stream in this region, Richland Creek incises into banks comprised of three stratigraphic layers associated with historic land use: pre‐European settlement, early European agriculture and development, and water‐powered milldam operation. This study aims to identify the bank processes along a reach of Richland Creek that is eroding towards its pre‐disturbance elevation. The volume of material that has eroded along this stream since the milldam breached was calculated by differencing a reconstructed surface of the pond bed and an aerial lidar digital terrain model (DTM). Immediately downstream from the study reach, the channel is floored by bedrock and immediately upstream the rate of channel erosion approximately doubled along the longitudinal profile of Richland Creek, which indicate that the study reach spans the transition from a channel dominated by vertical incision in the upstream direction to horizontal widening in the downstream direction. The combined hydrometeorological conditions and dominant processes causing reach‐scale cut bank erosion were investigated with analyses of stream stage, precipitation, and streambank volumetric and surfaces change that was measured during nine terrestrial lidar surveys in 2010–2012. The spatial variability of erosion during a simulated precipitation event was examined in a field‐based experiment. Erosion was greatest where mill pond sediment columns detached along vertical desiccation and horizontal seepage cracks. This sediment accumulated on the bank toe throughout the study and was a source of readily‐entrained fine sediment contrary to the upper reaches where depositional accommodation space is more limited. Findings suggest that hotspots of sediment excavation progress upstream, indicating that restoration efforts should focus upon stabilizing banks at these locations. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessment of Soil Erodibility in Taleghan Drainage Basin,Iran, Using Multivariate Statistics

Soil erosion has been recognized as one of the major forms of human-induced soil degradation. Due to land use changes in Iran, erosion has increased 800% between 1951 and 2002, calling for urgent action. But erosion research and policy development are hampered by a lack of information on the underlying factors controlling erosion. Soil types vary in their inherent susceptibility to erosion; but, like most countries, Iran lacks a network of field plots where erodibility is measured. A proxy for erodibility based on existing data and supplemented by an easily measured minimum data set is therefore needed. In this study, we use geological mapping and cluster, principal component, and factor analysis to group soils in the Taleghan Drainage Basin in Iran and subsequently determine their erodibility. First, a geological map of the area was prepared by photogeological methods and on-the-ground verification. Then, three soil profiles were investigated within similar landform units of each geological formation, and soil samples were taken. Physical and chemical properties that might impact soil erodibility (soil texture, pH, electrical conductivity, CaCO₃, and soil organic matter) were used to create a matrix of soil properties and parent material. Application of cluster analysis and factor analysis to the data allowed identification of three geological (parent material) clusters. To investigate the mutual effect of land use and parent material on soil erodibility, a soil erodibility factor was obtained for three land use types in each cluster: rangeland, cropland (irrigated), and dry-land farming (nonirrigated). Geological cluster 1, consisting of marl, gypsum, and gypsiferous mudstone, was the most erodible; geological cluster 2, consisting of recent alluvium, alluvial fan, and landslip deposits, was of intermediate erodibility; and geological cluster 3, consisting of igneous rocks, dolomite, and conglomerate, was the least erodible. Within each geological cluster, dry-land farming was the most erodible, cropland was medium erodible, and rangeland was least erodible. The study suggests that geological and land use maps provide a useful framework for assessing soil erodibility. This work can guide future soil erosion studies and direct soil conservation policy to areas most susceptible to erosion.     

Testing the Universal Soil Loss Equation‐MB equation in plots in Central and South Italy

Planning soil conservation strategies requires predictive techniques at event scale because a large percentage of soil loss over a long‐time period is due to relatively few large storms. Considering runoff is expected to improve soil loss predictions and allows relation of the process‐oriented approach with the empirical one, furthermore, the effects of detachment and transport on soil erosion processes can be distinguished by a runoff component. In this paper, the empirical model USLE‐MB (USLE‐M based), including a rainfall‐runoff erosivity factor in which the event rainfall erosivity index EI₃₀ of the Universal Soil Loss Equation (USLE) multiplies the runoff coefficient Q_R raised to an exponent b₁ > 1 is tested by the measurements carried out for the Masse (10 plots) and Sparacia (22 plots) experimental stations in Italy. For the Masse experimental station, an exponent b₁ > 1 was also estimated by tests carried out by a nozzle‐type rainfall simulator. For each experimental site in fallow conditions, the effect of the sample size of the plot soil loss measurements on the estimate of the b₁ coefficient was also studied by the extraction of a fixed number N of randomly obtained pairs of the normalized soil loss and runoff coefficient. The analysis showed that the variability of b₁ with N is low and that 350 pairs are sufficient to obtain a stable estimate of b₁. A total of 1,262 soil loss data were used to parameterize the model both locally and considering the two sites simultaneously. The b₁ exponent varied between the two sites (1.298–1.520), but using a common exponent (1.386) was possible. Using a common b₁ exponent for the two experimental areas increases the practical interest for the model and allows the estimation of a baseline component of the soil erodibility factor, which is representative of the at‐site soil intrinsic and quasi‐static properties. Development of a single USLE‐MB model appears possible, and sampling other sites is advisable to develop a single USLE‐MB model for general use.     

Effects of biological soil crusts on soil detachment process by overland flow in the Loess Plateau of China

Biological soil crusts (BSCs) cover up to 60 to 70% of the soil surface in grasslands after the ‘Grain for Green’ project was implemented in 1999 to rehabilitate the Loess Plateau. However, few studies exist that quantify the effects of BSCs on the soil detachment process by overland flow in the Loess Plateau. This study investigated the potential effects of BSCs on the soil detachment capacity (D_c), and soil resistance to flowing water erosion reflected by rill erodibility and critical shear stress. Two dominant BSC types that developed in the Loess Plateau (the later successional moss and the early successional cyanobacteria mixed with moss) were tested against natural soil samples collected from two abandoned farmland areas. The samples were subjected to flow scouring under six different shear stresses ranging from 7.15 to 24.08 Pa. The results showed that D_c decreased significantly with crust coverage under both moss and mixed crusts. The mean D_c of bare soil (0.823 kg m^?2 s^?1) was 2.9 to 48.4 times greater than those of moss covered soil (0.017–0.284 kg m^?2 s^?1), while it (3.142 kg m^?2 s^?1) was 4.9 to 149.6 times greater than those of mixed covered soil (0.021–0.641 kg m^?2 s^?1). The relative detachment rate of BSCs compared with bare soils decreased exponentially with increasing BSC coverage for both types of BSCs. The D_c value can be simulated by flow shear stress, cohesion, and BSC coverage using a power function (NSE ≥ 0.59). Rill erodibility also decreased with coverage of both crust types. Rill erodibility of bare soil was 3 to 74 times greater than those of moss covered soil and was 2 to 165 times greater than those of mixed covered soil. Rill erodibility could also be estimated by BSC coverage in the Loess Plateau (NSE ≥ 0.91). The effect of crust coverage on critical shear stress was not significant. Copyright © 2015 John Wiley & Sons, Ltd.