Journal of Earth Science Research
Journal of Earth Science Research(JESR)
Frequency: Annually
Watershed-level Rainfall Erosivity Mapping Using GIS-based Geostatistical Modeling
Due to the significant role of rainfall on soil erosion in one hand and a limited number of rain gauge stations to record daily rainfall, on the other hand, using modeling techniques are inevitable. The main objective of this research is to map rainfall erosivity index for Meymeh watershed of Ilam Province, Iran. Input data were collected at 26 rain gauge stations for 23 years (1987-2010). After selecting representative stations in the study area, monthly and annual rainfall was statistically reconstructed for the period. Then, using Fournier equation, the index was calculated for all stations. Next, GIS-based geostatistical methods were used to produce the maps of the rainfall erosivity factor. An intensive variography analysis was used to indicate a spatial correlation between the observed and the estimated values of rainfall erosivity. Weighed interpolation methods have been used for erosivity index mapping, including linear and nonlinear interpolations such as IDW and GPI. Besides, RBF, LPI and Kriging method were tested. The results indicate that the LPI method with RMSE = 51.28 had the least error. It seems that LPI is an efficient model for rainfall erosivity factor mapping. The simplicity of the method makes it easily retrievable.
Keywords:Watershed; Geo-statistical Modeling; Rainfall Erosivity; GIS; Iran
Author: Hassan Fathizad,Mohammad Ali Hakimzadeh,Seyed Rashid Fallah Shamsi,Soraya Yaghobi


  1. W. H. Wischmeier and D. D. Smith, “Predicting rainfall erosion. Losses: a guide to conservation planning,” Agriculture Handbook, US Department of Agriculture, Washington, DC, vol. 537, pp. 58, 1978.
  2. W. H. Wischmeier, “A soil erodibility nomo-graph for farmland and construction sites,” J Soil and Water Conserv, vol. 26, pp. 189-193, 1971.
  3. D. L. Corwin, M. Sorensen, and J. D. Rhoades, “Using GIS to locate salinity on irrigated soils,” Proc. 8th Conf. Computing in Civil Engineering in Conjunction with A.E.C System ‘92, TCCP. ASCE-Dallas, TX, June, vol. 7(9), pp. 468-485, 1992.
  4. E. Hosseini, J. Gallichand, and J. Caron, “Comparison of several interpolators for smoothing hydraulic conductivity data in southwest Iran,” ASAE, vol. 36(6), pp. 1687-1693, 1993.
  5. L. Z. Gitas, K. Douros, C. Minakou, G. N. Silleos, and C. G. Karydas, “Multi-temporal soil erosion risk assessment in N. Chalkidiki using a modified USLE raster model,” EARSesl eProccedings, vol. 8, no. 1, pp. 45-52, 2009.
  6. M. Karimi Nazar, A. Fakhireh, S. Feiznia, A. R. Rashki, and S. J. Miri Soleiman, “Evaluation of some geostatistical methods for estimating wind erosion threshold velocity in Sistan plain,” Range and Watershed Management- Natural Resources of Iran, vol. 62(3), pp. 405-417, 2009.
  7. S. Arekhi, Y. Niazi, and A. M. Kalteh, “Soil erosion and sediment yield modeling using RS and GIS techniques: a case study, Iran,” Arabian Journal of Geosciences, vol. 5(2), pp. 285-296, 2010.
  8. H. R. Matinfar, A. Shabani, and S. Aziz Qalati, “Investigate the spatial changes of some nutrients of soil using geostatistical methods (case study: Silakhor plain),” The Second National Conference on Agriculture and Sustainable Development, Opportunities and Challenges That Lie Ahead, no. 25, pp. 119-121, 2010.
  9. A. Shabani, H. R. Matinfar, S. Arkhay, and S. Rahimi Herabadi, “Modeling the rainfall erosivity factor using geostatistical method (case study: Ilam dam watershed),” Journal of Applied RS and GIS Techniques in natural resources Sciences, vol. 2(2), pp. 55-67, 2011.
  10. W. W. Hargrove, “Interpolation of Rainfall in Switzerland Using Regularized Splines with Tension, Geographic Information and Spatial Technologies Group Oak Ridge National Laboratory.” http: 2001.
  11. D. Kastelec and K. Kosmelj, “Spatial interpolation of mean yearly precipitation using universal Kriging,” Developments in Statistics, vol. 17, pp. 149-162, 2002.
  12. P. Goovaerts, “Geostatistical approach for incorporating elevation into spatial interpolation rainfall,” Journal of Hydrology, Amsterdam, vol. 228(1-2), pp. 133-129, 2000.
  13. M. E. N. Mingxin, Y. U. Zhenrong, and X. U. Hao, “Study on the spatial pattern of rainfall erosivity based on geo-statistics in Hebei Province, China, Front,” Agric. China, vol. 2(3), pp. 281-289, 2008.
  14. F. Onori, P. D. Bonis, and S. Grauso, “Soil erosion prediction at the basin scale using the revised universal soil loss equation (RUSLE) in a catchment of Sicily (southern Italy),” Environ Geol, vol. 50, pp. 1129-1140, 2006.
  15. J. W. Osborne, “Improving your data transformations: applying the Box-Cox transformation, practical assessment,” Research & Evaluation, vol. 15(12), pp. 1-9, 2010.
  16. K. G. Renard and J. R. Freimund, “Using monthly precipitation data to estimate the R factor in the revised USLE,” J Hydrol, vol. 157, pp. 287-306, 1994.
  17. B. Yu and C. J. Rosewell, “A robust estimator of the R factor for the universal soil loss equation,” Trans Am Soc Agric Eng, vol. 39(2), pp. 559-561, 1996.
  18. V. Ferro, G. Giordano, and M. Lovino, “Iso-erosivity and erosion risk map for Sicily,” Hydrol Sci J, vol. 36(6), pp. 549-564, 1991.
  19. R. P. C. Morgan, Soil Erosion and Conservation, Addison-Wesley, London, p. 198, 1995.
  20. A. M. Silva, “Rainfall erosivity map for Brazil,” Catena, vol. 57, pp. 251-259, 2004.
  21. H. M. J. Arnoldus, “Methodology used to determine the maximum potential average annual soils loss due to sheet and rill erosion in Morocco,” FAO Soils Bulletin, vol. 34, pp. 39-51, 1977.
  22. H. M. J. Arnoldus, “An approximation of the rainfall factor in the Universal Soil Loss Equation,” M. De Boodt, D. Gabriels, Eds., Assessment of Erosion, Chichester, New York, pp. 127-132, 1980.
  23. S. H. Hakimkhani, M. H. Mahdian, and M. Arabkhedri, “Mapping rainfall erosivity for Salt Lake Basin,” Journal of Natural Resources of Iran, vol. 60(6), pp. 712-726, 2008.
  24. M. E. Hohn, Geostatistics and Petroleum Geology, Kluwer Academic Publisher, Netherlands, p. 235, 1998.
  25. G. Y. Lu and D. W. Wong, “An adaptive inverse-distance weighting spatial interpolation technique,” Comp. Geosci, vol. 34, pp. 1044-1055, 2008.
  26. K. Hirsche, S. Boerner, C. Kalkomey, and C. Gastaldi, “Avoiding pitfalls in geostatistical reservoir characterization: A survival guide: The leading edge,” vol. 17, pp. 493-504, 1998.
  27. K. Johnston, J. M. Ver Hoef, K. Krivoruchko, and N. Lucas, “Using geostatistical analyst,” Environmental Systems Research Institute, Inc (ESRI), p. 335, 2001.
  28. G. Bohling, “Introduction to Geo-Statistics and Variogram Analysis, Assistant Scientist Kansas Geological Survey Ferro, V., Giordano, G., and Lovino, M., 1991. Iso-erosivity and erosion risk map for Sicily,” Hydrol Sci J, vol. 36(6), pp. 549-564, 2005.
  29. R. Webster and M. A. Oliver, Geostatistics for Environmental Scientists, Wiley Press, p. 271, 2000.
  30. A. A. Hassanipak, Geo-statistics, Tehran University, p. 380, 2007.
  31. A. R. Zabihi, K. Soleimani, M. Shabani, and S. Abravesh, “Examining the spatial distribution of annual rainfall using geostatistical methods (case study: Qom Province),” Physical Geography Research, vol. 78, pp. 112-101, 2011.
  32. K. J. Johnston, “Using ArcGIS geo-statistics analyst. New York: ESRI. KOULI, equation (RUSLE) in a catchment of Sicily (southern Italy),” Environ Geol, vol. 50, pp. 1129-1140, 2000.
  33. M. Vafakhah, M. Mohseni Saravi, M. Mahdavi, and S. K. Alavi Panah, “The use of geo-statistics for estimation of snow depth and density in Orazan watershed),” Iranian Journal of Watershed Management Science and Engineering, vol. 2(4), pp. 49-55, 2008.
  34. H. Habashi, S. M. Hosseini, Jahangard Mohammadi, and R. Rahmani, “The use of geo-statistics in the study of the soils of forest areas,” Journal of Agricultural Sciences and Natural Resources, Gorgan, vol. 14(1), pp. 11-20, 2006.
  35. S. Tuominen, S. Fish, and S. Poso, “Combining remote sensing, data from earlier inventories and geostatistical interpolation in multi-source forest inventory,” Can. J. For. Res, vol. 33, pp. 624-634, 2003.
  36. G. H. Zahedi Amiri, “Relation between ground vegetation and soil characteristic in a mixed hardwood stand,” Ph.D. thesis, University of Gent, Belgium, Academic Press, p. 319, 1998.
  37. K. Meusburger, A. Steel, P. Panagos, L. Montanarella, and C. Alewell, “Spatial and temporal variability of rainfall erosivity factor for Switzerland,” Hydrology and Earth System Sciences, vol. 16, pp. 167-177, 2012.
  38. N. Hudson, Soil Conservation, 3rd ed., B. T. Bats ford Limited, London, p. 391, 1995.
  39. E. A. Mikhailova, R. B. Bryant, S. J. Schwager, and S. D. Smith, “Predicting rainfall erosivity in Honduras,” Soil Science Society of America Journal, vol. 61(1), pp. 273-279, 1997.
  40. A. Shamshad, M. N. Azhari, M. H. Isa, W. M. A. Wan Hussin, and B. P. Parida, “Development of an appropriate procedure for estimation of RUSLE EI30 index and preparation of erosivity maps for Pulau Penang in Peninsular Malaysia,” CATENA, vol. 72, pp. 423-432, 2008.