چکیده:
زمینلغزش، جزء فراوانترین پدیدههای زمینشناسی در حوضههای آبخیز است و از آن به عنوان لندفرمی ژئوموفولوژیگ یاد میشود که سالانه خسارتهای فراوان جانی، مالی و منابع طبیعی در پی دارد. استان مازندران نیز از این قاعده مستثنی نیست. هدف از این تحقیق، پهنهبندی خطر زمینلغزش با استفاده از تئوری بیزین در حوضه آبخیز تالار واقع در استان مازندران است. در این روش، متغیرها از توزیعهای احتمال پیروی میکنند و تصمیمگیری بهینه را میتوان با استدلال بر توزیعهای احتمال و دادههای مشاهده شده انجام داد؛ به این منظور ابتدا با استفاده از نقاط لغزشی موجود، نقشه پراکنش زمینلغزشهای منطقه تهیه شد. سپس نقشه هر یک از عوامل مؤثر بر زمینلغزش شامل درجه شیب، جهت شیب، شکل شیب (تحب و تقعر شیب)، ارتفاع، کاربری اراضی، سنگشناسی، فاصله از جاده، فاصله از آبراهه، فاصله از گسل، شاخص رطوبت توپوگرافی (TWI)، شاخص طول شیب یا حمل رسوب (STI) و پوششگیاهی منطقه رسم شد. درنهایت نیز نقشه پهنهبندی حساسیت زمینلغزش با استفاده از تئوری بیزین برای منطقه مورد مطالعه در پنج کلاس خیلیکم، کم، متوسط، زیاد و خیلی زیاد تهیه شد. نتایج نشان داد که بیشتر زمینلغزشها مربوط به کلاس با شیب 30-15 درجه، با وزن 28/1 و از نظر جهت شیب در غرب و جنوب منطقه به ترتیب با وزن 2/2 و 41/2 رخ دادهاست. همچنین بیشتر زمینلغزشها، در شیب محدب و از نظر ارتفاعی در کلاس 1000-500 متر با وزن 07/4 شکل گرفتهاست. سنگهای ماسهای، شیلی و سیلتی و کاربری مرتع، بیشترین تأثیر را در زمینلغزش داشت. بررسی فاصله زمینلغزشها از جاده، آبراهه و گسل نیز نشان داد که بیشتر لغزشها در فاصله کمتر از 100-0 متری این عوامل رخ دادهاست که رابطه نزدیک بین آنها را نشان میدهد. بیشترین زمینلغزشها از نظر پوششگیاهی، مربوط به طبقه 5/0-3/0 با وزن 56/3 است. تحلیل رطوبت توپوگرافیک و شاخص حمل رسوب (توان آبراهه) نشان داد که بیشترین وزن لغزشی به ترتیب مربوط به طبقات 29/11-39/6 و 6/19-84/11 با وزنهای 009/1 و 17/1 است. علاوه بر این، نتایج نشان داد که زمینلغزشها به ترتیب در کلاسهای با حساسیت کم 40/25 درصد، متوسط 56/31 درصد، زیاد 72/28 درصد و خیلیزیاد 30/14 قرارگرفتهاست. اعتبارسنجی مدل نیز نشان داد که سطح زیرمنحنی ROCبرای مدل بیزین 56/85 است که در طبقه خیلی خوب قرار میگیرد.
1- Introduction
Mass movements are a type of morphodynamical phenomena that are usu ally related to various factors and occur on slopes in mountainous areas. Every year, damages caused by landslides lead to financial loss and death over the world (Fathi et al., 2018). According to previous studies, landslides cause 17% of the world's natural disasters. Mortality rates from 1903 to 2004 for different continents, including Asia, United States, Europe, Africa, and Australia were 29%, 39%, 30%, 1%, and 1%, respectively (Kohorest et al., 2005). Over the last three decades, several research studies have been conducted on landslide susceptibility mapping using different methods for developing their classification. All mapping methods are classified into five different groups including landslide distribution analysis, qualitative, statistical, deterministic, and frequency analyses (Vanvestern, 2003). Iran is prone to landslide phenomena due to natural conditions such as mountain topography, high tectonic and seismicity activity, geological and climatic diversity. Environmental factors affecting the occurrence of landslides are slope degree, aspect, plan curvature, elevation, land use, lithology, distance from the road, river and fault, topographic moisture index, slope length index or sediment transport and vegetation cover (Silakhori et al., 1400).
2- Methodology
Firstly, 134 points were identified using the Iranian landslide database and field survey in order to prepare a landslide susceptibility map using Bayesian theory. In the present study, 12 factors were used including elevation, slope and aspect, plan curvature, distance from the fault, road, and river, land use, geology, sediment transport index (STI), topographic wetness index (TWI), and vegetation cover. Then, topographic (1:50000), geology (1:100000), soil maps, and satellite imagery (Indian Remote Sensing) for 2012, were prepared and classified in ArcMap and ENVI environments. In order to evaluate Bayesian theory in landslide risk analysis, the relative performance curve of the relative efficiency of variables (ROC) was applied. This index is used to determine the accuracy and efficiency of the model (Egan, 1975; Williams et al., 1999). The area under the ROC curve represents the predicted value by describing its ability which accurately estimates the occurred events (landslide occurrence) and its non-occurred events (non-landslide occurrence). Therefore, the area under the curve is used as the model accuracy assessment. In the present study, 134 points of the landslide’s phenomena were used for modeling (70%) and accuracy assessment (30%) (Pourghasemi et al., 2013).
-RESULTS???
3- Discussion & Conclusions
The results of the factors affecting the occurrence of landslides using Bayesian theory in the study area showed that most of the landslides occurred in the class of 15-30 degrees with a weight of 1.28, which is also reported by Eracanoglu and Gokceolu (2004). The reason is that human intervention on these slopes causes more susceptibility (Yalcin et al., 2011). The study of aspect shows that most of the landslides occurred in the west and south directions with a weight of 2.21 and 2.41, respectively, which is confirmed by the results of Shams and Alizadeh (2019). The results of plan curvature showed that most of the landslides happened with a weight of 1.1 in convex slopes, which is close to the results of Pourghasemi (2013). Convex slopes usually have the highest landslides which is also reported in previous studies (Vanvewsten et al., 2003; Jaaferi et al., 2014). The altitude of 500-1000m covered by sandstones, chile and siltstones showed a significant relationship with a high number of landslides (Ayalew & Yamagishi, 2005). Among different land uses, rangeland (weight of 3.48) indicated the most significant relationship in landslide occurrence, which is reported by Shams and Alizadeh (2019).
The results of distance from roads, rivers, and faults showed that most landslides occurred at distances of more than 0-100 meters which confirms the results of previous studies (Pourghasemi & Mohammadi 2016). Also, the relationship between vegetation cover and landslides showed that the highest percentage occurred in class 0.5-0.3 with a weight of 3.56. In addition, the highest slippage for topographic wetness index and sediment transport index (river capacity) occurring were related to 6.39- 11.29 and 11.84-19.6 classes (weights of 1.009 and 1.17), respectively. The area under the curve was calculated at 85.56% in the model validation for the Bayesian model, which was classified as a very good performance. Therefore, the results of our study can play an important role in the management and planning of the Talar watershed.