Comparing theoretically supported rainfall‐runoff erosivity factors at the Sparacia (South Italy) experimental site

Interpreting rainfall‐runoff erosivity by a process‐oriented scheme allows to conjugate the physical approach to soil loss estimate with the empirical one. Including the effect of runoff in the model permits to distinguish between detachment and transport in the soil erosion process. In this paper, at first, a general definition of the rainfall‐runoff erosivity factor REF_e including the power of both event runoff coefficient Q_R and event rainfall erosivity index EI₃₀ of the Universal Soil Loss Equation (USLE) is proposed. The REF_e factor is applicable to all USLE‐based models (USLE, Modified USLE [USLE‐M] and Modified USLE‐M [USLE‐MM]) and it allows to distinguish between purely empirical models (e.g., Modified USLE‐M [USLE‐MM]) and those supported by applying theoretical dimensional analysis and self‐similarity to Wischmeier and Smith scheme. This last model category includes USLE, USLE‐M, and a new model, named USLE‐M based (USLE‐MB), that uses a rainfall‐runoff erosivity factor in which a power of runoff coefficient multiplies EI₃₀. Using the database of Sparacia experimental site, the USLE‐MB is parameterized and a comparison with soil loss data is carried out. The developed analysis shows that USLE‐MB (characterized by a Nash–Sutcliffe Efficiency Index NSEI equal to 0.73 and a root mean square error RMSE = 11.7 Mg ha^?1) has very similar soil loss estimate performances as compared with the USLE‐M (NSEI = 0.72 and RMSE = 12.0 Mg ha^?1). However, the USLE‐MB yields a maximum discrepancy factor between predicted and measured soil loss values (176) that is much lower than that of USLE‐M (291). In conclusion, the USLE‐MB should be preferred in the context of theoretically supported USLE type models.     

A new version of the USLE‐MM for predicting bare plot soil loss at the Sparacia (South Italy) experimental site

Improving empirical prediction of plot soil erosion at the event temporal scale has both scientific and practical importance. In this investigation, 492 runoff and soil loss data from plots of different lengths, λ (11 ≤ λ ≤ 44 m), and steepness, s (14.9 ≤ s ≤ 26.0%), established at the Sparacia experimental station, in Sicily, South Italy, were used to derive a new version of Universal Soil Loss Equation (USLE)‐MM model, by only assuming a value of one for the topographic length, L, and steepness, S, factors for λ = 22 m and s = 9%, respectively. An erosivity index equal to (Q_REI₃₀)^b1, Q_R and EI₃₀ being the runoff coefficient and the event rainfall erosivity index, respectively, with b₁ > 1 was found to be an appropriate choice for the Sparacia area. The specifically developed functions for L and S did not differ appreciably from other, more widely accepted relationships (maximum differences by a factor of 1.22 for L and 1.09 for S). The new version of the USLE‐MM performed particularly well for highly erosive events, because predicted soil loss differed by not more than a factor of 1.19 from the measured soil loss for measured values of more than 100 Mg ha^?1. The choice of the relationships to predict topographic effects on plot soil loss should not represent a point of particular concern in the application of the USLE‐MM in other environments. However, tests of the empirical approach should be carried out in other experimental areas in an attempt to develop analytical tools, usable at the event temporal scale, reasonably simple and of wide validity. Copyright © 2015 John Wiley & Sons, Ltd.     

Field investigation on rilling in the experimental Sparacia area,South Italy

In this paper the results of a field investigation on rilling carried out in the experimental Sparacia area are reported. The measurements were made on a plot 6 m wide and 22 m long subjected to natural rainfalls. For ten rainfalls the total soil loss (interrill and rill erosion) was collected in a storage system consisting of two tanks arranged in series at the base of the plot. Rill morphology (rill length and cross‐sections) was measured for five rainfall events, while the rill profile was surveyed for three events. First the contribution of each component (rill and interrill erosion) to total soil loss was established. Then the analysis allowed establishment of a power relationship between the rill length and the rill volume. Finally, for three events detailed information on rill erosion and rill morphology allowed verification of the applicability of WEPP and estimation of the rill erodibility constant. Copyright © 2007 John Wiley & Sons, Ltd.     

Testing the use of an image‐based technique to measure gully erosion at Sparacia experimental area

The first part of this investigation was aimed at testing the use of a three‐dimensional (3D) digital terrain model and a quasi‐tridimensional (2.5D) digital elevation model obtained by a large series of oblique images of eroded channels taken from consumer un‐calibrated and non‐metric cameras. For two closed earth channels having a different sinuosity, the ground measurement of some cross sections by a profilometer (P) was carried out and their real volume was also measured. The comparison among the three methods (3D, 2.5D, and P) pointed out that a limited underestimation of the total volume always occurs and that the 3D method is characterized by the minimum difference between measured and real volume. For this reason, 3D model can be used as benchmark. In the subsequent part of the investigation, the three ground measurement methods were applied for surveying of an ephemeral gully (EG) channel at the Sparacia area. The morphological and hydraulic variable values of the 24 surveyed cross sections determined by both 2.5D model and profilometer were compared. This comparison showed that the estimate error is generally less than ±10%. The EG measurements carried out by the three methods supported the applicability both of the empirical relationship between EG length and its eroded volume and the theoretical dimensionless relationship among the morphological variables describing the channelized erosion process. Finally, it was demonstrated that the effect of the distance interval on the EG volume measurement by 3D and 2.5D models is negligible for the investigated EG.     

A comprehensive analysis of Universal Soil Loss Equation-based models at the Sparacia experimental area

Improving Universal Soil Loss Equation (USLE)-based models has large interest because simple and reliable analytical tools are necessary in the perspective of a sustainable land management. At first, in this paper, a general definition of the event rainfall- runoff erosivity factor for the USLE-based models, REF_e = (Q_R)^b1(EI₃₀)^b2, in which Q_R is the event runoff coefficient, EI₃₀ is the single-storm erosion index, and b₁ and b₂ are coefficients, was introduced. The rainfall-runoff erosivity factors of the USLE (b₁ = 0 and b₂ = 1), USLE-M (b₁ = b₂ = 1), USLE-MB (b₁ ≠ 1 and b₂ = 1), USLE-MR (b₁ = 1 and b₂ ≠ 1), USLE-MM (b₁ = b₂ ≠ 1), and USLE-M2 (b₁ ≠ b₂ ≠ 1) can be defined using REF_e. Then the different expressions of REF_e were simultaneously tested against a data set of normalized bare plot soil losses, A_eN, collected at the Sparacia (south Italy) site. As expected, the poorest A_eN predictions were obtained with the USLE. The observed tendency of this model to overestimate small A_eN values and underestimate high A_eN values was reduced by introducing in the soil loss prediction model both Q_R and an exponent for the erosivity term. The fitting to the data was poor with the USLE-MR as compared with the USLE-MB and the USLE-MM. Estimating two distinct exponents (USLE-M2) instead of a single exponent (USLE-MB, USLE-MR, and USLE-MM) did not appreciably improve soil loss prediction. The USLE-MB and the USLE-MM were recognized to be the best performing models among the possible alternatives, and they performed similarly with reference to both the complete data set and different sub-data sets, only including small, intermediate, and severe erosion events. In conclusion, including the runoff coefficient in the soil loss prediction model is important to improve the quality of the predictions, but a great importance has to be paid to the mathematical structure of the model.     

SCALES: a large‐scale assessment model of soil erosion hazard in Basse‐Normandie (northern‐western France)

The cartography of erosion risk is mainly based on the development of models, which evaluate in a qualitative and quantitative manner the physical reproduction of the erosion processes (CORINE, EHU, INRA). These models are mainly semi‐quantitative but can be physically based and spatially distributed (the Pan‐European Soil Erosion Risk Assessment, PESERA). They are characterized by their simplicity and their applicability potential at large temporal and spatial scales. In developing our model SCALES (Spatialisation d'éChelle fine de l'ALéa Erosion des Sols/large‐scale assessment and mapping model of soil erosion hazard), we had in mind several objectives: (1) to map soil erosion at a regional scale with the guarantee of a large accuracy on the local level, (2) to envisage an applicability of the model in European oceanic areas, (3) to focus the erosion hazard estimation on the level of source areas (on‐site erosion), which are the agricultural parcels, (4) to take into account the weight of the temporality of agricultural practices (land‐use concept). Because of these objectives, the nature of variables, which characterize the erosion factors and because of its structure, SCALES differs from other models. Tested in Basse‐Normandie (Calvados 5500 km²) SCALES reveals a strong predisposition of the study area to the soil erosion which should require to be expressed in a wet year. Apart from an internal validation, we tried an intermediate one by comparing our results with those from INRA and PESERA. It appeared that these models under estimate medium erosion levels and differ in the spatial localization of areas with the highest erosion risks. SCALES underlines here the limitations in the use of pedo‐transfer functions and the interpolation of input data with a low resolution. One must not forget however that these models are mainly focused on an interregional comparative approach. Therefore the comparison of SCALES data with those of the INRA and PESERA models cannot result on a convincing validation of our model. For the moment the validation is based on the opinion of local experts, who agree with the qualitative indications delivered by our cartography. An external validation of SCALES is foreseen, which will be based on a thorough inventory of erosion signals in areas with different hazard levels. Copyright © 2010 John Wiley & Sons, Ltd.     

Distance,time and scale in soil erosion processes

This commentary brings together, as a virtual Special Issue, a number of recent papers in Earth Surface Processes and Landforms that are all related to issues of scale in soil erosion. Empirical concepts that were developed in the 1940s are now in need of re‐thinking, and papers are increasingly exploring, through modelling and measurement, appropriate ways to recognize the mechanisms that connect processes across time and space scales. Issues include a more nuanced approach to selective transportation, responses to variability in surface and sub‐surface conditions and the need to analyse measurements in ways that can be transferred between sites and storms. Copyright © 2010 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.     

Spatial and temporal variations in soil erosion and deposition due to land‐levelling in a semi‐arid area of Basilicata (Southern Italy)

This paper examines the changes from 1955 to 2002 in soil erosion and deposition due to changes in land‐use patterns in the semi‐arid territory of Craco, which is characterized by landsliding and badland erosion. The area underwent continuous degradation during the last century due not only to its lithological vulnerability but also to the anthropic pressure favoured by the introduction of Common Agricultural Policy (CAP) measures, which has led to the reclamation of scrub lands and badlands for durum wheat cultivation. Our analysis integrates the Unit Stream Power Erosion Deposition (USPED) model with a geographic information system (GIS) to quantify erosion risk and predict deposition patterns. Soil data, land use inventory, digital elevation data and climatic atlases were used as resource data sets to generate USPED factor values. The obtained results correlate well with field measured erosion data by other researchers. In the investigated 47 years, stable areas decreased by about 280 ha (3·8% of the total surface area), largely attributable to the increase of the low and moderate erosion intensity without significant change in sedimentation. Results from this study have implications related to understanding the geomorphic response of sites that were abandoned following remodelling due to the application of the F measure of Regulation CEE 2078/92. The average annual erosion rates estimated for abandoned and remodelled sites are respectively 15·99 and 10·64 t ha^?1, meaning that the total amount of erosion in 20 years could be estimated at around 100 t ha^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Site‐scale variability of streambank fluvial erodibility parameters as measured with a jet erosion test

The erosion rate of cohesive streambanks is typically modelled using the excess shear stress equation, dependent on two erodibility parameters: critical shear stress and erodibility coefficient. The jet erosion test (JET) has become the most common method for estimating these erodibility parameters in situ. Typically, results from a few JETs are averaged to acquire a single set of parameters for characterizing a streambank layer; however, this may be inadequate for accurately characterizing erodibility. The research objectives were to investigate the variability of JET results from assumed homogeneous streambank layers and to estimate the number of JETs required to accurately characterize erodibility for use in predictive models. On three unique streambanks in Oklahoma and across a range of erodibility, 20 to 30 JETs were conducted over a span of three days at each site. Unique to this research, each JET was analysed using the Blaisdell, scour depth and iterative solutions. The required sample size to accurately estimate the erodibility parameters depended on the JET solution technique, the parameter being estimated, and the degrees of precision and confidence. Conducting three to five JETs per soil layer on a streambank typically provided an order of magnitude estimate of the erodibility parameters. Because the parameters were log‐normally distributed, using empirical equations to predict erosion properties based on soil characteristics will likely contain high uncertainty and thus should be used with caution. This study exemplifies the need to conduct in situ measurements using the JET to accurately characterize streambank resistance to fluvial erosion. Copyright © 2015 John Wiley & Sons, Ltd.     

Implications of model uncertainty for the mapping of hillslope‐scale soil erosion predictions

An integrated modelling approach (MIRSED) which utilizes the process‐based soil erosion model WEPP (Water Erosion Prediction Project) is presented for the assessment of hillslope‐scale soil erosion at five sites throughout England and Wales. The methodology draws upon previous uncertainty analysis of the WEPP hillslope soil erosion model by the authors to qualify model results within an uncertainty framework. A method for incorporating model uncertainty from a range of sources is discussed as a first step towards using and learning from results produced through the GLUE (Generalized Likelihood Uncertainty Estimation) technique. Results are presented and compared to available observed data, which illustrate that levels of uncertainty are significant and must be taken into account if a meaningful understanding of output from models such as WEPP is to be achieved. Furthermore, the collection of quality, observed data is underlined for two reasons: as an essential tool in the development of soil erosion modelling and also to allow further constraint of model uncertainty. Copyright © 2001 John Wiley & Sons, Ltd.     

Supporting USLE‐MM reliability by analyzing soil loss measurement errors

Sampling the collected suspension in a storage tank is a common procedure to obtain soil loss data. A calibration curve of the tank has to be used to obtain actual concentration values from those measured by sampling. However, literature suggests that using a tank calibration curve was not a common procedure in the past. For the clay soil of the Sparacia (Italy) experimental station, this investigation aimed to establish a link between the relative performances of the USLE‐M and USLE‐MM models, usable to predict plot soil loss at the event temporal scale, and soil loss measurement errors. Using all available soil loss data, lower soil loss prediction errors were obtained with the USLE‐MM (exponent of the erosivity term, b₁ > 1) than the USLE‐M (b₁ = 1). A systematic error of the soil loss data is unexpected for the Sparacia soil because the calibration curve does not depend on the water level in the tank. In any case, this type of error does not have any effect on the b₁ exponent. Instead, this exponent decreases as the level of underestimation increases for increasing soil loss values. This type of error can occur at Sparacia if it is assumed that a soil loss measurement can be obtained by a bottle sampler dipped close to the bottom of the tank after mixing the suspension and assuming that the measured concentration coincides with the actual one. In this case, the risk is to obtain a lower b₁ value than the actual one. In conclusion, additional investigations on the factors determining errors in soil loss data collected by a sampling procedure are advisable because these errors can have a noticeable effect on the calibrated empirical models for soil loss prediction.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

Measuring rill erosion at plot scale by a drone‐based technology

The traditional direct method (i.e. metric ruler and rillmeter) of monitoring rill erosion at plot scale is time consuming and invasive because it modifies the surface of the rilled area. Measuring rill features using a drone‐based technology is considered a non‐invasive method allowing a fast field relief. In the experimental Sparacia area a survey by a quadricopter Microdrones md4‐200 was carried out, and this relief allowed the generation of a Digital Elevation Model (DEM), with a mesh size of 1 cm and a resolution elevation equal to 2 mm, for three plots (L, G and C) affected by rill erosion. At first for the experimental L plot, which is 44 m long, the rill features were surveyed by a ‘manual’ method which was carried out by drawing on the PC screen the rill paths obtained by a visual orthophoto interpretation. This manual method was not applicable for the plots in which rills of limited depth occurred and were not detectable. Then, for both L plot and the other experimental plots having a length ranging from 22 to 44 m, an ‘automatic’ extraction method of rills from DEM was applied. Using an appropriate calculation routine, a vector coverage of transects orthogonal to the main flow direction (i.e. the maximum slope steepness path) was generated. The intersection of each plot DEM with the transect coverage allowed to obtain both the cross sections and the main rill morphological features. For the L plot the comparison between the rill morphological features obtained by the two different methods (manual, automatic) was carried out. Finally, the length–volume relationship and a dimensionless relationship proposed in literature were tested for all studied experimental plots. Copyright © 2015 John Wiley & Sons, Ltd.     

Erodibility of soil and organic matter: independence of organic matter resistance to interrill erosion

The enrichment of organic matter in interrill sediment is well documented; however, the respective roles of soil organic matter (SOM) and interrill erosion processes for the enrichment are unclear. In this study, organic matter content of sediment generated on two silts with almost identical textures, but different organic matter contents and aggregations, was tested. Artificial rainfall was applied to the soils in wet, dry and crusted initial conditions to determine the effects of soil moisture and rainfall and drying history on organic matter enrichment in interrill sediment. While erosional response of the soils varied significantly, organic matter enrichment of sediment was not sensitive to initial soil conditions. However, enrichment was higher on the silt with a lower organic matter content and lower interrill erodibility. The results show that enrichment of organic matter in interrill sediment is not directly related to either SOM content or soil interrill erodibility, but is dominated by interrill erosion processes. As a consequence of the complex interaction between soil, organic matter and interrill erosion processes, erodibility of organic matter should be treated as a separate variable in erosion models. Further research on aggregate breakdown, in particular the content and fate of the organic matter in the soil fragments, is required. Copyright © 2007 John Wiley & Sons, Ltd.     

Exploring effects of rainfall intensity and duration on soil erosion at the catchment scale using openLISEM: Prado catchment,SE Spain

In semi‐arid areas, high‐intensity rainfall events are often held responsible for the main part of soil erosion. Long‐term landscape evolution models usually use average annual rainfall as input, making the evaluation of single events impossible. Event‐based soil erosion models are better suited for this purpose but cannot be used to simulate longer timescales and are usually applied to plots or small catchments. In this study, the openLISEM event‐based erosion model was applied to the medium‐sized (～50 km²) Prado catchment in SE Spain. Our aim was to (i) test the model's performance for medium‐sized catchments, (ii) test the ability to simulate four selected typical Mediterranean rainfall events of different magnitude and (iii) explore the relative contribution of these different storms to soil erosion using scenarios of future climate variability. Results show that because of large differences in the hydrologic response between storms of different magnitudes, each event needed to be calibrated separately. The relation between rainfall event characteristics and the calibration factors might help in determining optimal calibration values if event characteristics are known. Calibration of the model features some drawbacks for large catchments due to spatial variability in K_sat values. Scenario calculations show that although ～50% of soil erosion occurs as a result of high frequency, low‐intensity rainfall events, large‐magnitude, low‐frequency events potentially contribute significantly to total soil erosion. The results illustrate the need to incorporate temporal variability in rainfall magnitude–frequency distributions in landscape evolution models. Copyright © 2011 John Wiley & Sons, Ltd.     

Measuring sheet erosion using synthetic color‐contrast aggregates

The bulk of eroded soils measured at the outlets of plots, slopes and watersheds are suspended sediments, splash‐induced sheet erosion. It is depending on rainfall intensity and antecedent soil moisture contents and contributes to a significant proportion of soil loss that usually is ignored in soil erosion and sediment studies. A digital image processing method for tracing and measuring non‐suspended soil particles detached/transported by splash/runoff was therefore used in the present study. Accordingly, fine mineral pumice grains aggregated with white cement and coloured with yellow pigment powder, with the same size, shape and specific gravity as those of natural soil aggregates, called synthetic color‐contrast aggregates, were used as tracers for detecting soil particle movement. Subsequently, the amount of non‐suspended soil particles detached and moved downward the slope was inferred with the help of digital image processing techniques using MATLAB R2010B software (Mathworks, Natick, Massachusetts, USA). The present study was conducted under laboratory conditions with four simulated rainfall intensities between 30–90 mm h^‐1, five antecedent soil moisture contents between 12–44 % v v^‐1 and a slope of 30%, using sandy loam soils taken from a summer rangeland in the Alborz Mountains, Northern Iran. A range of total transported soil between 90.34 and 1360.93 g m^‐2 and net splash erosion between 36.82 and 295.78 g m^‐2were observed. The results also showed the sediment redeposition ratio ranging from 87.27% [sediment delivery ratio (SDR) = 12.73%] to 96.39% (SDR = 3.61%) in various antecedent soil moisture contents of rainfall intensity of 30 mm h^‐1 and from 80.55% (SDR = 19.45%) to 89.42% (SDR = 10.58%) in rainfall intensity of 90 mm h^‐1. Copyright © 2013 John Wiley & Sons, Ltd.     

Relating aerial erosion,soil erosion and sub‐soil erosion to the evolution of Lunan Stone Forest,China

Stone forest (‘Shilin’ in Chinese) is a unique karst landform with a complex evolution process. Based mainly on the characteristics and interrelationships of sub‐soil, soil and sub‐aerial erosion in Lunan karst area, the authors develop a triplex erosion model to describe the evolution of stone forest, and apply it to examine the current development stage and the prospect of the Lunan Stone Forest. The study shows that sub‐soil corrosion, a basic driving force for the vertical scope of a stone forest, usually occurs within 10 m below ground surface but is observed to be most active within the top 2 m, which constitutes the best development zone for stone forest. Under modern climatic conditions, the tip of the stone pillars in Lunan karst area is lowering at a rate of 10·4 mm ka^?1, whereas the base of the stone pillars is deepening at 26·17 mm ka^?1. Therefore, the height of stone pillars is increasing at a rate of 15·77 mm ka^?1. Considering that soil erosion in the study area is as high as 650 mm ka^?1, the visible height of the stone forest is actually increasing at a rate of 639·6 mm ka^?1. However, the best evolution time for Lunan Stone Forest has already passed despite the fact that it is still growing taller at the present time. This is because the soil layer, which plays an extremely significant role in the heightening of stone pillars, is rapidly thinning at a rate of 623·83 mm ka^?1. Copyright © 2006 John Wiley & Sons, Ltd.