Search

The Universal Soil Loss Equation (USLE) and the Revised Universal Soil Loss Equation (RUSLE) have been widely used for predicting average soil loss. Slope length is an important topographical parameter of the L factor in USLE/RUSLE. Among the widely studied GIS procedures for extracting slope length, the distributed watershed erosion slope length (DWESL) based on the unit contributing area estimation method, which considers two-dimensional runoff process and cutoff factors, is a relatively complete model for calculating slope length. Slope length in the DWESL model is primarily calculated using conventional flow direction algorithms such as D8, Dinf, MS and MFD-md. However, DWESL outputs require further improvement due to the errors in the usual estimates of the uphill contributing area and the effective contour length of discrete elements. Combined with a theoretical differential equation of specific catchment area on hillsides, the calculation of the DWESL model was optimized without estimating the uphill contributing area or the effective contour length for each cell. The proposed integration method based on the topographical features slope line, contour curvature and cutoff factors (ITF method) was used to extract slope length from the raster digital elevation. Slope length extracted using the ITF method had the smallest error in verification of mathematical surfaces (average RRMSE = 0.0573), and its spatial distribution was more consistent with the structure of the terrain surface for all test data, relative to the conventional flow direction algorithms in the original DWESL model. The proposed ITF method could provide a reference for predicting soil erosion using the USLE/RUSLE model.

The Universal Soil Loss Equation (USLE) and its revised version (RUSLE) are often used to estimate soil erosion at regional landscape scales. USLE/RUSLE contain parameters for slope length factor (L) and slope steepness factor (S), usually combined as LS. However a major limitation is the difficulty in extracting the LS factor. Methods to estimate LS based on geographic information systems have been developed in the last two decades. L can be calculated for large watersheds using the unit contributing area (UCA) or the slope length (λ) as input parameters. Due to the absence of an estimation of slope length, the UCA method is insufficiently accurate. Improvement of the spatial accuracy of slope length and LS factor is still necessary for estimating soil erosion. The purpose of this study was to develop an improved method to estimate the slope length and LS factor. We combined the algorithm for multiple-flow direction (MFD) used in the UCA method with the LS-TOOL (LS-TOOLSFD) algorithms, taking into account the calculation errors and cutoff conditions for distance, to obtain slope length (λ) and the LS factor. The new method, LS-TOOLMFD, was applied and validated in a catchment with complexly variable slopes. The slope length and LS calculated by LS-TOOLMFD both agreed better with field data than with the calculations using the LS-TOOLSFD and UCA methods, respectively. We then integrated the LS-TOOLMFD algorithm into LS-TOOL developed in Microsoft's.NET environment using C# with a user-friendly interface. The method can automatically calculate slope length, slope steepness, L, S, and LS factor, providing the results as ASCII files that can be easily used in GIS software and erosion models. This study is an important step forward in conducting accurate large-scale erosion evaluation.

A hydrologically correct digital elevation model (DEM) forms a basis for realistic environmental modelling, especially in complex terrain. We have performed a study in the Coarse Sandy Hilly Catchments (CSHC) of the Loess Plateau, China, which demonstrates pragmatic, yet effective methods for improving the quality of the DEM by: (1) identifying and correcting source topographic data errors and (2) optimising ANUDEM algorithm parameters. Improvement in the DEM based on fixing over 1100 errors in the input topographic data, and optimising key ANUDEM parameters was assessed using higher accuracy independent validation of 32 contributing areas and 1474 spot heights, and by semi-quantitative analysis of DEM derivatives produced from ANUDEM and Triangular Irregular Network (TIN) algorithms. Improvement in the ANUDEM DEM over the original TIN DEM was shown where the percentage of the total absolute difference in contributing areas reduced from 10.43 to 3.51%, and the bias between the spot heights and DEM elevations reduced from 45 to 32 m. Large improvement in DEM quality was gained by using ANUDEM instead of TIN, with smaller improvement gained by fixing source data errors, and optimising ANUDEM parameters.