Engineering Geology

Engineering Geology

Gully Erosion Hazard Zoning in the Gharnaveh Watershed, Golestan Province

Authors
Kazem Saber Chenari
Abolareza Bahremand1
Vahed Berdi Sheikh
Chooghi Bairam Komaki
10.18869/acadpub.jeg.13.1.69
Abstract
Introduction

One of the main problems in the Golestan province watersheds is the high degree of erosion and soil degradation, so that the equilibrium between the soil process and the soil erosion is unbalanced, and the erosion rate increases from west to east. Among these, the gully erosion and piping have the highest role. Gully is a canal or stream with the headcut with active erosion, sharpened slope and steep walls that results from the destruction of surface flow (usually during or after the occurrence of precipitation), dissolution movements, and small mass movements. The extent of gully in the eastern parts of Golestan province has caused the land degradation of arable land and landscape and has increased the conservation cost and etc. Because of connecting upstream areas of the basin to the downstream areas, gully has particular importance, which provides the possibility of sediment and pollutant transport, road destruction and financial losses to agricultural lands. In order to prevent and control the development of gully processes from a small scale to large one, it is a versatile utility to identify and extract the areas prone to gully erosion.

Due to the high intensity of gully erosion and its increasing growth in the Gharnaveh watershed, the Garnaveh River has an unstable status and severe eroded gully, and in some areas it has a great depth and vertical lateral walls, as well. Therefore, in this research, the watershed of Garnaveh was selected to prepare the risk areas of gully erosion.

The aim of this research is to determine Gully Erosion Hazard zoning using Frequency Ratio and Gupta & Joshi methods (Gully Nominal Risk Factor-GNRF) in the Garnaveh watershed (Golestan province). Ultimately, the accuracy of the model has been evaluated using quality sum method and Kappa coefficient.

Material and methods

The study area is located in the northern part of Iran, Golestan province. The Garnaveh watershed with an area of about 78430 hectares lies between longitudes 370360 E and 414472 E, and latitudes of 4183819 N and 4155267 N (UTM Zone 40).

At first, gully erosion inventory map with the scale of 1:75,000 (dependent variable) for the Gharnaveh watershed has been prepared using multiple field surveys and satellite images. From total gullies, 70% have been selected randomly for building gully erosion hazard zoning model and the remaining ones (30%) have been used to validate the provided model.

In this research, seven data layers including slope percent, slope aspect, plan curvature, lithology formation, land use types, distance from rivers and distance from roads have been selected as gully erosion controlling factors (covariates/ independent variables) and then they have been digitized in ArcGIS software. The amount of Gully density of each factor class has been calculated from a combination of independent and dependent variables, and the rating of classes have done based on Frequency Ratio and Gully Nominal Risk Factor equations. Finally, the Gully erosion hazard zoning map has been drawn from the summation of weighting maps in ArcGIS. In this map, the value of each pixel is calculated by summing the weights of all the factors in that pixel. The pixel values are categorized based on the natural breaks classifier into very low, low, medium, high and very high hazard zones. Then, an accuracy of zoning map has been evaluated by quality sum method and Kappa coefficient.

Results and discussion

The result of affecting factors classification of the gullies shows that loess deposits formation, rangeland, areas with low distance from road and rivers, northwest aspect, low slope amplitude and concave slopes contain the most susceptibility to gullying. The results of frequency percent comparison of gullies in hazard classes show that from all gully zones in the validation step of the GNRF and frequency ratio models %74.52 and %78.11 of zones are located in the high and very high risk classes, respectively. The result of model validation using the quality sum method and a Kappa coefficient show that the frequency ratio model is a more appropriate model for gully erosion hazard zoning (with the quality sum and a Kappa coefficient of 3 and 0.89, respectively) than the GNRF model (having the quality sum and Kappa coefficient of 1.27 and 0.74, respectively).

Conclusion

In this research, the areas susceptible to gully erosion in the Gharnaveh watershed have been mapped with the frequency ratio and GNRF (for the first time) models. For this purpose, 7 affecting factors (independent variable) and 805 gully zones (dependent variable) were provided to measure the hazard maps of gully erosion. The following results are obtained from this study.

- The geology factors were identified as the most effective factors in the occurrence of gully erosion in the Gharnaveh watershed.

- Based on the gully erosion zoning hazard map of the Gharnaveh watershed, more than 70 percent of gullies are situated in the very high and high hazard classes.

- The produced gully erosion hazard map is useful for planners and engineers to reorganize the areas susceptible to gully erosion hazard, and offers appropriate methods for hazard reduction and management, as well../files/site1/files/131/4Extended_Abstract.pdf
Keywords
Hazard zoning
Gully erosion
Gharnaveh watershed
Frequency ratio
Gopta & Joshi
1. احمدی ح.، " ژئومورفولوژی کاربردی"، جلد 1 فرسایش آبی، انتشارات دانشگاه تهران (1385) 688 . ## 2. دماوندی م. ز.، "بررسی عارضه‌های مورفولوژیکی نهشته‌های لسی در بخش غربی شهر گرگان"، گزارش طرح پژوهشی، دانشگاه علوم کشاورزی و منابع طبیعی گرگان (1385) 45. ## 3. Poesen J., Nachtergaele J., Verstraeten G., Valentin, C., "Gully eosion and environmental change: importance and research needs", Catena. Vol.50 (2003) 91-133. ## 4. تجری ع.، "تأثیر خصوصیات فیزیکی و شیمیایی خاک در فرسایش خندقی و مورفولوژی خندق‌ها (مطالعه موردی: حوضه آبخیز کال‌شور کلاله استان گلستان)"، پایان‌نامۀ کارشناسی ارشد، دانشگاه علوم کشاورزی و منابع طبیعی گرگان. (1391) 122. ## 5. راهی غ. ر.، نظری‌سامانی ع. ا.، احمدی ح.، سلاجقه ع. "تأثیر ویژگی‌های خاک در نوع سازوکار ایجاد خندق و مورفولوژی آن در منطقه گناوه"، مجلۀ منابع طبیعی ایران، دورۀ 62، شمارۀ 4 (1388) 472-459. ## 6. Gómez Gutiérrez A., Schnabel S., Lavado, C., "Use of two nonparametric methods (CART and MARS) to model the potential distribution of gullies in Spanish rangelands", International Congress on Environmental Modeling and Software. (2008) 1398-1405. ## 7. جعفری گرزین ب، کاویان ع.، "ارزیابی وقوع فرسایش خندقی در حوزه آبخیز سرخ‌آباد مازندران با استفاده از سنجش از دور و سامانه اطلاعات جغرافیایی"، مجلۀ علمی پژوهشی علوم و مهندسی آبخیزداری ایران، شمارۀ 17 (1388) 58-55. ## 8. شادفر ص.، "بررسی فرسایش خندقی با استفاده از روش تحلیل سلسله مراتبی در شهرستان رودبار استان گیلان"، پژوهش‌های فرسایش محیطی، شمارۀ 3 (1390) 16-30. ## 9. ملکی ا.، احمدی م.، میلادی ب.، "شبیه‌سازی مناطق مستعد خندق‌زایی با استفاده از روش SPI در رودخانه مرگ"، پژوهش‌های ژئومورفولوژی کمی، دورۀ 1، شمارۀ 3 (1391) 23-38. ## 10. اسفندیاری درآباد ف.، بهشتی جاوید ا.، فتحی م. ح.، "بررسی قابلیت فرسایش خندقی با استفاده از مدل منطق فازی مطالعۀ موردی، حوضۀ سد گلستان- رودخانه قرناوه"، مجموعه مقالات دومین کنفرانس ملی مخاطرات محیطی، دانشگاه خوارزمی (1392). ## 11. شادفر ص.، "کاربرد منطق فازی در بررسی فرسایش خندقی با استفاده ازسیستم اطلاعات جغرافیایی (مطالعه موردی: حوضه طرود)، فصلنامه علمی- پژوهشی اطلاعات جغرافیایی"، شماره 92 (1393) 35-42. ## 12. طهماسبی‌پور ن.، رحمتی الف.، قربانی‌نژاد، س.، "پیش‌‌بینی حساسیت به فرسایش خندقی منطقه سیمره براساس مدل عامل قطعیت و تعیین اهمیت عوامل مؤثر بر آن"، مجلۀ اکوهیدرولوژی، دورۀ 3، شمارۀ 1 (1395) 83-93. ## 13. شهریور ع.، شادفر ص.، خزایی م.، عادلی ب.،"ارزیابی روش‌های پهنه‌بندی فرسایش خندقی (مطالعۀ موردی: حوزۀ آبخیز آبگندی)"، مجلۀ اکوهیدرولوژی، دورۀ 4، شمارۀ 1 (1396) 119-132. ## 14. Rahmati O., Haghizadeh A., Pourghasemi H. A., Noormohamadi F., "Gully erosion susceptibility mapping: the role of GIS-based bivariate statistical models and their comparison", Natural Hazards, Vol.82 (2016) 1231-1258 . ## 15. Lesschen J. P., Kok K., Verburg P. H., Cammeraa L. H., "Identification of Vulnerable Areas for Gully Erosion under Different Scenarios of Land Abandonment in Southeast Spain", Catena, Vol.71 (2007) 110-121. ## 16. Luca F., Conforti M., Robustelli G., "Comparison of GIS-based gullying susceptibility mapping using bivariate and multivariate statistics: Northern Calabria, South Italy", Geomorphology, Vol.134 (2011) 297-308. ## 17. Ogbonna J. U., "Understanding Gully Erosion Vulnerability in Old ImoState Using Geographic Information System and Geostatistics", American Journal of Geographic Information System, Vol.1 (2012) 66-71. ## 18. Conoscenti C., Agnesi V., Angileri S., Cappadonia C., Rotigliano E., Marker M., "A GIS-based approach for gully erosion susceptibility modelling: a test in Sicily, Italy", Environtal Earth Science, Vol.70 (2013) 1179-1195. ## 19. Conoscenti C., Angileri S., Cappadonia C., Rotigliano E., Agnesi V., Marker M., "Gully erosion susceptibility assessment by means of GIS-based logistic regression: a case of Sicily (Italy)", Geomorphology, 204 (2014) 399-411. ## 20. Dube F., Nhapi I., Murwira A., Gumindoga W., Goldin J., Mashauri D.A., "Potential of weight of evidence modeling for gully erosion hazard assessment in Mbire District–Zimbabwe", Physics and Chemistry of the Earth, Vol.67 (2014) 145-152. ## 21. Dewitte O., Daoudi M., Bosco C. Eeckhaut, M.V.D., "Predicting the susceptibility to gully initiation in data-poor regions", Geomorphology, Vol.228 (2015)101-115. ## 22. اداره کل منابع طبیعی و آبخیزداری استان گلستان.، "مطالعات طراحی سازه‌های کنترل سیل و رسوب حوزه آبخیز قرناوه"، شرکت مشاور کاوش پی مشهد، (1390) 2-10. ## 23. Andrew O., Ian P. Prosser., "Gully erosion prediction across a large region: Murray–Darling Basin, Australia"و Soil Research, Vol.4 (2012) 267-277. ## 24. Ozdemir A., Altural T., "A comparative study of frequency ratio, weights of evidence and logistic regression methods for landslide susceptibility mapping: Sultan Mountains, SW Turkey", Journal of Asian Earth Sciences, Vol.64 (2013)180-197. ## 25. Ozdemir A., "GIS-based groundwater spring potential mapping in the Sultan Mountains (Konya, Turkey) using frequency ratio, weights of evidence and logistic regression methods and their comparison", Journal of Hydrology, Vol.411(2011) 290-308. ## 26. Lee S., "Application and verification of fuzzy algebraic operators’ landslide susceptibility mapping", Environmental Geology, Vol.52 (2007) 615-623. ## 27. Pradhan B., Lee S., "Landslide susceptibility assessment and factor effect analysis: backpropagation artificial neural networks and their comparison with frequency ratio and bivariate logistic regression modelling", Environmental Modelling & Software. Vol.25 (2010) 747-759. ## 28. صابرچناری ک.، سلمانی ح.، محمدی م.، "ارزیابی خطر زمین‌لغزش با استفاده از مدل‌های ارزش اطلاعات و LNRF"، نشریۀ اکوهیدرولوژی، شمارۀ 2 (1394) 105-116. ## 29. گرایی پ.، کریمی ح.، "تعیین مناسب‌ترین روش پهنه‌بندی خطر زمین لغزش در حوضۀ آبخیز سد ایلام"، دورۀ 25، شمارۀ 96 (1389) 128-101. ## 30. Gupta R. P., Joshi B. C., "Landslide hazard zonation using the GIS approach- a case study from the Ramganga Catchment, Himalayas", Engineering Geology. Vol.28 (1990) 119-131. ## 31. طولابی س؛ عابدینی م.، علی‌بخشی الف.، "مقایسه و ارزیابی مدل‌های LNRF و AHP در پهنه‌بندی خطر زمین لغزش حوضه آبخیز نوژیان"و فصلنامۀ بین‌المللی پژوهشی تحلیلی زمین پویا. دورۀ 2، شمارۀ 2 (1393) 41-31. ## 32. شیرانی ک.، سیف ع.، "پهنه‌بندی خطر زمین‌لغزش با استفاده از روش‌های آماری (منطقۀ پیشکوه، شهرستان فریدون شهر)"، مجلۀ علوم زمین. دورۀ 22، شمارۀ 8 (1390) 158-149. ## 33. Gee M.D., "Classification of landslide hazard zonation methods and a test of predictive capability", Landslides, Bell (ed.), Balkema, Rotterdam (1991) 947-952. ## 34. Baeza C., Lantada N., Moya J., "Validation and evaluation of two multivariate statistical models for predictive shallow landslide susceptibility mapping of the Eastern Pyrenees (Spain)", Environtal Earth Science. Vol. 61 (2010) 507-523. ## 35. Landis J. R., Koch G. G., "The measurement of observer agreement for categorical data", Biometrics, Vol. 33 (1977) 159-174. ## 36. Mandal S., Maiti R., "Semi-quantitative Approaches for Landslide Assessment and Prediction, Chapter 7, Application of Analytical Hierarchy Process (AHP) and Frequency Ratio (FR) Model in Assessing Landslide Susceptibility and Risk", Natural Hazards (2014) 191-226. ## 37. انتظاری م.، ملکی الف.، مرادی خ.، الفتی س.، "پهنه‌بندی فرسایش خندقی در حوضه آبریز دیره با استفاده از روش تحلیل سلسله مراتبی"، برنامه‌ریزی و آمایش فضا، شمارۀ 17، (1392) 63-86. ## 38. یمانی م.، زمان‌زاده س. م.، احمدی م.، "تحلیل عوامل مؤثر بر شکل‌گیری و توسعه فرسایش خندقی (مطالعه موردی: حوضه دشت کهور در استان فارس)"، مجلۀ کاوش‌های جغرافیایی مناطق بیابانی، شمارۀ 1، (1392) 53-84.. ## 39. Conforti M., Aucell C., Robustelli G., Scarciglia F., "Geomorphology and GIS analysis for mapping gully erosion susceptibility in the Turbolo stream catchment (Northern Calabria, Italy)", Natural Hazards, Vol. 58 (2011) 881-898. ## 40. Dai F. C., Lee C. F., Xu Z. W., "Assessment of landslide susceptibility on the natural terrain of Lantua sland, Hong Kong", Environment Geology, Vol. 40 (2001) 381-391. ## 41. Umar Z., Pradhan B., Ahmad A., Jebur M. N., Tehrany M. S., "Earthquake induced landslide susceptibility mapping using an integrated ensemble frequency ratio and logistic regression models in West Sumatera Province, Indonesia", Catena.Vol. 118 (2014) 124-135. ## 42. Tehrany M. S, Pradhan B, Jebur M. N., "Flood susceptibility mapping using a novel ensemble weightsof-evidence and support vector machine models in GIS", Journal of Hydrology. Vol. 512 (2014) 332-343. ## 43. مخدوم فرخنده م.، "شالوده آمایش سرزمین"، انتشارات دانشگاه تهران (1389) 289. ## 44. قهرودی م.، "مدل پهنه‌بندی خطر ناشی از گسترش فرسایش آبکنـدی در حـوزه آبکنـد کلوچـه بیجار با استفاده از GIS و "RS، طرح تحقیقاتی وزارت نیرو (1382) 52-53. ##