Document Type : Research Paper

Authors

1 Department of Soil Science, Faculty of Agriculture, Shahid Chamran University of Ahavz

2 Department of soil science, Faculty of Agriculture, Shahid Chamran University of Ahvaz

3 Department of Soil science, Faculty of Agriculture, Shahid Chamran University of Ahvaz

Abstract

Introduction
In recent decades, the potential of soil erosion as serious threat against sustainable land management (SLM) and soil sustainability has been recognized. Furthermore, human-induced soil erosion lead to harmful environmental effects and transportation of sediment to water bodies is accompanied by loss of nutrients and eutrophication. Therefore, there is a need to focus on soil erosion outcomes to prevent its environmental impacts and mitigate the negative feedbacks of soil erosion. One of the most important factors that affecting the amount of erosion and sediment yield is soil erodibility (K-factor in USLE), which is an index of soil susceptibility. Moreover, K-factor is one of the 6 effective factors on soil erosion in universal soil loss equation and mainly is representative of soil properties. Regarding the difficulty of soil erodibility ´s measurement, therefore we need to an easily measurable parameter. In this study soil aggregate is used as an index of soil erodibility factor. Generally, topography alters from one slope to another one, which conduces to dramatically changes of soil erosion intensity, therefore topography in terms of geomorphometric parameters is effective on soil erosion processes; morphometric analysis and soil erosion modeling are robustly interconnected. Also, digital elevation model (DEM) is a base to extract geomorphometric parameters and the quality of DEM is effective on geomorphometric characteristics.
Methods and material
This study conducted with the aim of modelling the aggregate stability using multivariate linear regression and K-means clustering algorithm as an index of soil erodibility factor. For this purpose a study area with an approximate area of 20000 hectares selected. This area is located in Dehdez which is in the north east of Khuzestan province. Based on the climatologic data this area has xeric and thermic regimes. The main land uses in this area including rangeland and forest and in the rangeland the crop cover is mostly grasses. According to the purpose of this study and degree of homogeneity the soil sampling accomplished using completely randomized method. From each land use 25 soil samples collected (the depth of soil sampling was 0-30 cm). Then the soil samples prepared and the mean weight diameter (MWD) as an index of soil aggregate stability using wet sieving, soil texture with hydrometer method and organic carbon content measured in at least 3 replications. In another part of this project the digital elevation model (DEM) of study area extracted from SRTM (Shuttle radar topographic mission) with the spatial resolution of
30 m. After DEM preparation all preprocessing operations performed on DEM using Arc Hydro extension in ArcGIS 10.3 software in order to remove sinks and enhance the quality of DEM. In the next step the primary and secondary geomorphometric parameters extracted from this DEM using ArcGIS 10.3 and SAGA 6.3. software. Finally, multivariate linear regression and K-means clustering models developed between geomorphometric parameters and soil properties as dependent variables and aggregate stability as independent variable to achieve statistical indices in order to evaluate developed models.
Results and discussion
Our results showed that there was a significant correlation between some geomorphometric parameters and soil aggregate stability; therefore there is a possibility to apply these parameters in order to parameterize aggregate stability. Moreover, the results of modelling indicate that regression models using geomorphometric parameters and soil properties was able to cover approximately 75 and 70 percentage of spatial variation of aggregate stability in the rangeland and forest land use of the study area, respectively. While clustering-regression models were able to explain 77 and 82 percentage of the spatial variation of the aggregates stability in the first and second cluster respectively. Also, the results of validation of developed models in this study showed that the root mean square error (RMSE) of regression models for rangeland and forest land use was 0.33 and 0.26 respectively and RMSE of clustering-regression models for first and second clusters was 0.93 and 0.62 respectively. According to importance of soil aggregate stability on soil erodibility factor (K-factor) and the difficulty of soil erodibility ´s measurement therefore these developed models are useful tools in order to predict soil erosion and based on the distribution of soil erosion and deposition using STI (Sediment transport index) and TWI (Topographical wetness index) could be able to select and apply the best management practices in the critical areas. Indeed STI map and TWI map are indices of spatial distribution of soil erosion and deposition in the studied area therefore using these geomorphometric indices we able to control soil erosion and its negative feedbacks.
Conclusion
Regarding the difficulty of soil erodibility ´s measurement therefore we used some indices to simplify this process and our results illustrated that it is possible to develop some regression models in order to estimate aggregate stability as an index of soil erodibility. In summary our results confirm that geomorphometric parameters are easily available parameters based on the DEM to predict soil erodibility.
Keywords: Digital elevation model (DEM), geomorphometry, modelling, soil erosion, topographical indices

Keywords

  1. References

    1. Barthes, B. and Roose, E. 2002. Aggregate stability as an indicator of soil susceptibility to runoff and erosion validation at several levels. Catena, 47(2002):133-149.
    2. Besalatpour, A.A., Ayoubi, S., Hajabbasi, M.A., Mosaddeghi, M., and Schulin, R. 2013. Estimating wet soil aggregate stability from easily available properties in a highly mountainous watershed. Catena, 111: 72-79.‏
    3. Bronick, C.J. and Lal, R. 2005. Soil structure and management: a review. Geoderma, 124(1): 3-22.
    4. Canasveras, J.C., Barrón, V., Del Campillo, M.C., Torrent, J., and Gómez, J.A. 2010. Estimation of aggregate stability indices in Mediterranean soils by diffuse reflectance spectroscopy. Geoderma, 158(1-2): 78-84.‏
    5. Caron, J., Espindola, C.R., and Angers, D.A. 1996. Soil structural stability during rapid wetting: influence of land use on some aggregate properties. Soil Science Society of America Journal, 60(3): 901-908.‏
    6. Carter, M.R. 1992. Influence of reduced tillage systems on organic matter, microbial biomass, macro-aggregate distribution and structural stability of the surface soil in a humid climate. Soil and Tillage Research, 23(4): 361-372.‏
    7. Deng, Z.Q., de Lima, J.L. and Jung, H.S. 2008. Sediment transport rate-based model for rainfall-induced soil erosion. Catena, 76(1): 54-62.‏
    8. Florinsky, I. 2016. Digital terrain analysis in soil science and geology. Second ed., Academic Press,‏ Amsterdam, PP: 7-68.‏
    9. Gee, G.W. and Bauder, J.W. 1986. Particle-size Analysis. In Klute A. (Ed). Methods of Soil Analysis: Part 1—Physical and Mineralogical Methods. American Society of Agronomy, Madison, PP: 383–441.
    10. Gobin, A., Campling, P., and Feyen, J. 2001. Soil-landscape modelling to quantify spatial variability of soil texture. Physics and Chemistry of the Earth, Part B: Hydrology, Oceans and Atmosphere, 26(1): 41-45.‏
    11. Golchin, a. 2017. Soil organic matter. Jahad Daneshgahi press, Zanjan. 98 pp.
    12. Gopp, N.V., Nechaeva, T.V., Savenkov, O.A., Smirnova, N.V., and Smirnov, V.V. 2017. The methods of geomorphometry and digital soil mapping for assessing spatial variability in the properties of agrogray soils on a slope. Eurasian Soil Science, 50(1): 20-29.‏
    13. Jakšík, O., Kodešová, R., Kubiš, A., Stehlíková, I., Drábek, O., and Kapička, A. 2015. Soil aggregate stability within morphologically diverse areas. Catena, 127: 287-299.‏
    14. Odeha, I.O.A., McBratney, A.B., and Chittleborough, D.J. 1994. Spatial prediction of soil properties from landform attributes derived from a digital elevation model. Geoderma, 63(3-4): 197-214.‏
    15. Pierson, F.B. and Mulla, D.J. 1990. Aggregate stability in the Palouse region of Washington: effect of landscape position. Soil Science Society of America Journal, 54(5): 1407-1412.‏
    16. Renschler, C.S. and Harbor, J. 2002. Soil erosion assessment tools from point to regional scales—the role of geomorphologists in land management research and implementation. Geomorphology, 47(2): 189-209.
    17. Walkley, A. and Black, I.A. 1934. An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Science, 37(1): 29-38.‏
    18. Wang, G.G., Gartner, X., Liu, H., and Anderson, A. 2001. Uncertainty assessment of soil erodibility factor for revised universal soil loss equation. Catena, 46: 1-14.
    19. Wilson, J.P. and Gallant, J.C. 2000. Digital terrain analysis. In Wilson, J. P. and Gallant, J. C. (Eds). Terrain analysis: principles and applications. John Wileys and Sons, New York, PP: 1-22.
    20. Yoder, R.E. 1936. A direct method of aggregate analyses and a study of the physical nature of erosion losses. Agronomy Journal, 28:337–351.
    21. Zadorova, T., Jakšík, O., Kodešová, R. and Penížek, V. 2011. Influence of terrain attributes and soil properties on soil aggregate stability. Soil Water Research, 6(3): 111-119.
    22. Zhang, H.Y., Shi, Z.H., Fang, N.F., and Guo, M.H. 2015. Linking watershed geomorphic characteristics to sediment yield: Evidence from the Loess Plateau of China. Geomorphology, 234: 19-27.