Document Type : Research Paper
Authors
Department of Geomatics Engineering, Faculty of Civil Engineering, University of Tabriz, Tabriz, Iran
Abstract
Highlights
The maximum land subsidence rate of about -79 mm per year in the direction of line of sight of the sensor is seen using SBAS time series method. The maximum subsidence is seen in the southern and southeastern parts of the plain. The obtained land subsidence susceptibility map was classified into 5 classes; namely the areas with the highest sensitivity to subsidence (Very High: 3%), areas with high probability of land subsidence (High: 20%), areas with moderate probability of land subsidence (Medium: 25%) and areas with low and very low susceptibility to land subsidence: (Low: 33%) and (Very low: 19%) classes in the land subsidence susceptibility map.
In this article, two methods have been used to verify the accuracy of the land subsidence sensitivity map that follows: a) Using displacement map produced by InSAR time series analysis method. The area under ROC curve showed the prediction accuracy of 86% b) Sensitivity analysis to the conditioning factors with highest weights. According to the performed validations, the effectiveness of the proposed method for evaluating land subsidence in Shabstar Plain was confirmed.
Conclusion
According to the InSAR time series results, land subsidence is more visible in the southern areas of Shabastar Plain. In the generated land subsidence susceptibility map, these areas were mainly located in the class ‘very high’ sensitivity to land subsidence too. Based on AHP method applied for weighing conditioning factors, among the selected factors, 'underground water level changes', 'geology' and 'land use' were the factors with highest weights; 36.8 %, 22.6% and 18.7% accordingly. The important role of the three factors in causing land subsidence was investigated and the existence of a direct and significant relationship between the important conditioning factors and land subsidence that occurred in the Shabstar plain in the period from 2018 to 2021 was confirmed.
By using the proposed method, it is possible to get a correct understanding of the region in terms of land subsidence occurred and prediction of land subsidence occurrence according to a set of conditioning factors and concerning their weights. This paper showed that the combination of SAR interferometry time series and spatial analysis can be an efficient and low-cost method to evaluate land subsidence in land subsidence prone areas. The method used in this paper can be useful for other areas where there is lack of enough field data about the occurred land subsidence or its preparation is time consuming and expensive. The results of this study can help managers and related authorities plan for necessary actions accurately and timely in the prevention and mitigation phase.
Keywords
Main Subjects
پدیده فرونشست به عنوان جابجایی سطح زمین در جهت عمودی تعریف میشود که میتواند ناشی از فعالیتهای انسانی یا پدیدههای طبیعی باشد. طبق تعریف یونسکو، فرونشست زمین، فرونشست یا نشست سطح زمین است که به دلایل مختلف در مقیاس بزرگ رخ میدهد و میتواند به تدریج یا ناگهانی رخ دهد (لهستان، ۱۹۸۴). فرونشست زمین میتواند خسارات فیزیکی و اجتماعی خاصی را ایجاد کند. بنابراین، ارزیابی به موقع و قابل اعتماد این پدیده در مناطق آسیبپذیر ضروری است. در مطالعه حاضر، دشت شبستر به عنوان یکی از دشتهای پرخطر واقع در استان آذربایجان شرقی به عنوان منطقه مورد مطالعه انتخاب شده است. دادهها و روشها در این مقاله، با استفاده از روش تحلیل سریهای زمانی SBAS، نقشهای از نرخ جابجایی سطح پوسته دشت شبستر بین سالهای ۲۰۱۸ تا ۲۰۲۱ تهیه شد. این دشت تقریباً در طول جغرافیایی ۴۵ درجه و ۱۱ دقیقه تا ۴۶ درجه و ۳ دقیقه و عرض جغرافیایی ۳۸ درجه و ۷ دقیقه تا ۳۸ درجه و ۲۳ دقیقه قرار دارد. اولین گام در تهیه نقشه جابجایی با استفاده از تداخلسنجی راداری (SAR)، انتخاب مجموعه مناسب از تصاویر SAR است. در این مطالعه از تصاویر Sentinel-1A استفاده شده است. از الگوریتم زیرمجموعه خط مبنای کوچک (SBAS) استفاده شد. 143 تداخلنگار برای 46 تصویر Sentinel-1A با استفاده از نرمافزار GMTSAR ایجاد شد، با این معیار که خطوط مبنای زمانی کمتر از 100 روز و خطوط مبنای مکانی کمتر از 100 متر باشند. نقشه جابجایی زمین دشت شبستر که نرخ جابجایی سالانه را نشان میدهد، به دست آمد. برای مدلسازی حساسیت فرونشست زمین منطقه مورد مطالعه، عوامل مؤثر بر فرونشست زمین با توجه به مطالعات قبلی انتخاب شدند. عوامل با استفاده از روش سلسله مراتبی تحلیلی (AHP) وزندهی شدند. نقشههای عوامل طبقهبندی و زیرکلاسها بر اساس نظرات کارشناسان رتبهبندی شدند. از همپوشانی وزنی برای ادغام همه عوامل با در نظر گرفتن وزن هر عامل و رتبه زیرکلاسها در محیط GIS استفاده شد. با ترکیب وزنی معیارها و زیرمعیارهای بهدستآمده در مراحل قبل، نقشه حساسیت فرونشست زمین تهیه شد.