Document Type : Research Paper

Authors

1 M.Sc. Graduated, Rangeland and Watershed Department, Water and Soil Faculaty, University of Zabol, Zabol, Iran

2 Assistant Professor, Rangeland and Watershed Department, Water and Soil Faculaty, University of Zabol, Zabol, Iran

3 Assistant Professor, Soil and Water Research Department, Golestan Agricultural and Natural Resources Research and Education Center, Agricultural Research, Education and Extension Organization (AREEO), Gorgan, Gorgan, Iran

Abstract

Introduction Spatial variations of soil properties is a natural event, which recognizing these changes is inevitable in order to planning and right management of both agricultural and natural resources. Soil organic carbon (SOC) is the most important factor in soil fertility and quality, climate change and reduction of greenhouse gas emissions. Furthermore, evaluating the rates and spatial distribution of the soil properties, land improvement and restoration can be traced from the carbon sequestration index. According to the above, providing quantitative and qualitative conservation of soil properties such as SOC can be considered an effective way to achieve sustainable development of natural and environmental resources. Digital soil mapping (DSM) can determine the spatial variations of soil organic carbon by exploring the relationship between soil properties and effective environmental variables. Different statistical models such as regression trees and random forest are used in order to communicate between soil characteristics and its spatial distribution. The present study was carried out to investigate the spatial distribution of SOC, as well as, to determine the most important variables affecting their prediction in Ravang watershed in Minab County.
Materials and MethodsRavang watershed with an area of 13821.6 hectares is located in Hormozgan province, Minab city. The maximum and minimum elevations are 357 and 33 meters, respectively. Digital Elevation Model of Ravang watershed was used to extract 17 environmental covariates (such as elevation, aspect, slope, valley depth,…) by SAGA software (http://www.gdem.aster.ersdac.or). Moreover, two environmental covariates related to remote sensing including Normalized Difference Salinity Index (NDSI) and Normalized Difference Vegetation Index (NDVI) were determined in the study area. In addition, the maps of land use, sand, silt, clay and pH were used as covariates in modeling. In order to determining the location of sampling points, the conditioned hyper-cube technique was used. After determining of soil sample location, field sampling was carried out at a depth of 0-30 cm. then, 100 soil samples were taken and the amount of SOC was measured. Random forest model was applied to the relationship between SOC and covariates. The model includes two user-defined parameters, including the number of variables used in the construction of each tree, which expresses the power of each independent tree and the number of trees in each forest. Considering the strength of independent trees, the predictive accuracy of the model increases, conversely, the correlation between them will decrease. The accuracy of the soil organic carbon distribution was also evaluated using root mean square error (RMSE), mean error (ME) and correlation coefficient (R2), which were determined.
Results and Discussion Based on the present study results, elevation, soil silt and sand maps, channel network base level, slope and NDVI are the most important factors on predicting the of SOC variations. The results indicated that RMSE, ME and R2 were 0.36, 0.26 and 0.38, respectively .Results also showed that affecting erosion and sediment, as well as, human effect, have the most impact on the SOC soil spatial distribution in the Ravang watershed. Moreover, result show SOC deficiency in the soil of Ravang watershed due to high salinity, low percent of vegetation cover and land use changes. In addition, drought intensifying and decrease in precipitation have reduced SOC content, which itself causes changes in the texture and chemical properties of the soil and, as a consequence, makes them more susceptible to erosion.
Conclusion The variability of SOC is very high in the study area because of intensive water erosion and land use change. Overall, the results of the present study indicated that the critical condition of soil organic carbon in the Ravang watershed, which requires a comprehensive management of the region's water and soil resources to improve soil conditions and increase the reserves of this important and influential variable in the soil structure. On the other hand, despite of the acceptable performance of the random forest model in estimating of soil properties, due to high variability of some soil properties, model prediction performance may be decreased. 

Keywords

  1. Aksoy, E., and Dirim, M.S. 2009. Soil mapping approach in GIS using Landsat satellite imagery and DEM data. African Journal of Agricultural Research, 4(11): 1295-1302.
  2. Breiman, L. 2001. Random forests. Machine Learning, 45 (1): 5-32. 102
  3. Dai, E.F., Zhai, R.X., Ge, Q.S., and Wu, X. 2014. Detecting the storage and change on topsoil organic carbon in grasslands of Inner Mongolia from 1980s to 2010s. Acta Geographica Sinica, 24: 6. 1035-1046.
  4. Debella-Gilo, M. and B. Etzelmuller. 2009. Spatial prediction of soil classes using digital terrain analysis and multinomial logistic regression modeling integrated in GIS: Examples from Vest fold County, Norway. Catena,77 (1): 8-18.
  5. Farage, P.K., Ardo, J., Olsson, L., Rienzi, E.A., Ball, A.S., and Pretty, J.N. 2007. The Potential of Soil Carbon Sequestration in Three Tropical Dryland Farming Systems of Africa and Latin America: A modelling approach. Soil and Tillage Research, 94 (2): 457- 472.
  6. Franklin, J.1995. Predictive Vegetation Mapping: Geographic Modelling of Bio spatial Patterns in Relation to Environmental Gradients, Progress in Physical Geography, 19(4): 474-499.
  7. Grimm, R., Behrens, T., Marker, M., and Elsenbeer, H. 2008. Soil organic carbon concentrations and stocks on Barro Colorado Island - Digital soil mapping using random forests analysis. Geoderma, 146(1-2):102-113.
  8. Guo, P.T., H.B. Liu., and Wu, W. 2009. Spatial prediction of soil organic matter using terrain attributes in a hilly area. International Conference on Environmental Science and Information Application Technology,Wuhan, China. 3:762-759.
  9. Hengl, T., Heuvelink, G.B.M., Kempen. B.,Leenaars, J.G.B., Walsh, M.G., Shepherd K.D., Sila, A., MacMillan, R.A., Jesus, J.M.D., Tamene, L., and Tondoh, J.E. 2015. Mapping Soil Properties of Africa at 250 m Resolution: Random Forests Significantly Improve Current Predictions. PLoS ONE 10(6): e0125814 doi:10.1371/journal.pone.0125814.
  10. Jiang, P., and Thelen, K.D. 2004. Effect of soil and topographic properties on crop yield in a north-central corn soybean cropping system. Agronomy Journal, 96: 252-258.
  11. Kadkhodapur, M. A. Mirjalili, A.B., and Dashtkian, K. 2009. The study of the permeability status of Miankow aquifer field in natural and planted fields, Watershed Research, (84): 58-53.
  12. Li, X., McCarty, G.W., Karlen, D.L., and Cambardella, C.A. 2018. Topographic metric predictions of soil redistribution and organic carbon in Iowa cropland fields. Catena, 160: 222-232.
  13. Ließ, M., Glaser, B., and Huwe, B. 2012. Uncertainty in the spatial prediction of soil texture comparison of regression tree and random forest models. Geoderma, 170: 70-79.
  14. Malone B.P., McBratney A.B., Minasny B., and Laslett G.M. 2009. Mapping continuous depth functions of soil carbon storage and available water capacity. Geoderma, 154: 138-152.
  15. Mc Bratney, A.B., Santos, M.L.M., and Minasny, B. 2003. On digital soil mapping. Geoderma, 117: 3-52.
  16. Melendez-Pastor, I., Navarro-Pedreño, J., Gómez, I., and Koch, M. 2008. Identifying optimal spectral bands to assess soil properties with VNIR radiometry in semi-arid soils. Geoderma, 147 (3-4):126-132.
  17. Milanesi, P., Holderegger, R., Caniglia, R., Fabbri, E., and Randi, E. 2016. Different habitat suitability models yield different least-cost path distances for landscape genetic analysis. Basic and Applied Ecology, 17(1): 61-7.
  18. Mirakzehi, K.L., Shahriyari, L., Pahlavan Rad, M.R., and Bameri, A. 2015. Application random forest method in prediction of soil classes in low elevation (casa study: Hirmand city). Journal of Soil and Water Conservation Research, 24: 67-84.
  19. Nelson, D.W. and Sommers, L. 1982. Total carbon, organic carbon, and organic matter 1. Methods of soil analysis. Part 2. Chemical and microbiological properties, (methodsofsoilan2), 539-579.
  20. Nicodemus, K.K. 2011. Letter to the Editor: On the stability and ranking of predictors from random forest variable importance measures. Briefings in Bioinformatics, 12: 369-373.
  21. Niknahad Gharmakher, H., and M. Maramaei. 2011. Effects of land use changes on soil properties (Case study: the Kechik catchment). Journal of Soil Management and Sustainable, 1(2): 81-96. (In Persian).
  22. Nourozi, P., Nadiri, A.A., Asghari Moqadam, A., and Qqareh Khani, M. 2015. Projection of aquifer capabilities of Malekan Palin. Knowledge of Water and Soil, 27(2): 61-75. (In Persian).
  23. Naderi, H., Hedayati Zadeh, R., and Doroodi, H. 2007. Physiographic properties (elevation and gradient) on the amount of organic carbon and total nitrogen in soil, 10th Iranian soil science congress, Karaj, College of Agriculture and Natural Resources, University of Tehran.
  24. Pahlavan Rad, M.R., Toumanian, N., and Khormalali, F. 2017. Introduction of digital mapping of soil. Land Management, 4.2 (2): 97-114.
  25. Pahlavan Rad, M.R., Dahmardeh, K.H., and Brungard, C. 2019. Predicting soil organic carbon concentrations in a low relief landscape, eastern Iran. Geoderma Regional, 15. e00195.
  26. Pajand, M.G., Imami, H., and Astaraie, A. 2014. Relationship between topography and some soil features. Soil and Water Journal (Agricultural Sciences and Technology), 29(6): 1699- 1710.
  27. Pandey, C.B., Chaudhari, S.K., Dager, J.C., Singh, G.B., and Singh, R.K. 2010. Soil N mineralization and microbial biomass carbon affected by different tillage levels in a hot humid tropic, Soil and Tillage Research, 11: 33-41.
  28. Peters, J., Verhoest, N., Samaon, R., Boeckx, P., and De Baets, B. 2008. Wetland vegetation distribution modelling for the identification of constraining environmental variables. Landscape Ecology, 23: 1049-1056.
  29. Pham, T.G., Kappas, M., Huynh, C.V., Nguyen, H. 2019. Application of ordinary kriging and regression kriging method for soil properties mapping in hilly region of central Vietnam. International Journal of Geo-Information, 8 (3), 147; doi:10.3390/ijgi8030147.
  30. Post, W.M. 1993. Organic carbon in soil and global carbon cycle. P. 216-302. In: M. Heiman (ed). The Global Carbon Cycle, NATO ASI Series 1: Global Environmental Change. Springer, New York.
  31. Schwender, H. 2008. Identification of SNP interactions using logic regression. Biostatistics 9(1):187-198.
  32. Sidari, M., Ronzello, G., Vecchio, G., and Muscolo, A. 2008. Influence of slope aspects on soil chemical and biochemical properties in a Pinus laricio forest ecosystem of Aspromonte (Southern Italy). European Journal of Soil Biology, 44: 364-372.
  33. Stoorvogel, J.J., Kempen, B., Heuvelink, G.B.M., and de Bruin, S. 2009. Implementation and evaluation of existing knowledge for digital soil mapping in Senegal. Geoderma. 149: 169-170.
  34. Taghizadeh-Mehrjardi, R. 2016. Modern concepts in Soil Science (Pedometric). Ardakan Univ. Press, 311p. (In Persian).
  35. Teng, J., Xiang, T., Huang, Z., Wu, J., Jiang, P., Meng, C., Li, Y., and Fuhrmann, J.J. 2015. Spatial distribution and variability of carbon storage in different sympodial bamboo species in China. Journal of Environmental Management, 168: 46-52.
  36. Vaysse, K. and P. Lagacherie. 2015. Evaluating Digital Soil Mapping approaches for mapping Global Soil Map soil properties from legacy data in Languedoc-Roussillon (France). Geoderma Regional, 4: 20-30.
  37. Zhang, X.C., Zhang, J.L., and Li, S.J. 2009. Spatial variability of soil organic carbon in a watershed on the loess plateau. Pedosphere, 19: 486-495.
  38. Were, K., Bui, D. T., Disk, B., and Singl, B. R. 2015. A comparative assessment of support vector regression, artificial neural networks, and random forest for predicting soil organic carbon stocks across an afromonkane land scape. Ecological Indicator, 394-403.
  39. Wilson, J.P., and Gallant, J.C. 2000. Teain analysis, principle and applications. John Wiley and Sons, Inc. NY.
  40. Zeraatpisheh, M., Ayoubi, S., Jafari, A., Tajik, S., and Finke, P. 2019. Digital mapping of soil properties using multiple machine learning in a semi-arid region, central Iran. Geoderma, 338 (15): 445-452.
  41. Zhao, B., Li, Z., Li, P., Xu, G., Gao, H., Cheng, Y., and Fang, Z. 2017. Spatial distribution of soil organic carbon and its influencing factors under the condition of ecological construction in a hilly-gully watershed of the losse plateau, china. Geoderma, 296:10-17.