Finding the potential of karst water resources in Birjand city using hierarchical analysis model and geographic information system
Subject Areas : GeohydrologyMohsen Rezaeiarefi 1 , Mohammad Motamedirad 2 , Ibrahim Taghavimoghadam 3 , Mansour Sadeghi 4
1 - Assistant Professor, Geography Education Department, Farhangian University, Tehran, Iran
2 - Assistant Professor, Geography Education Department, Farhangian University, Tehran, Iran
3 - Assistant Professor, Geography Education Department, Farhangian University, Tehran, Iran
4 - PhD in Geomorphology and Environmental Management, Kharazmi University, Tehran, Iran
Keywords: Prospecting, karst water resources, hierarchical analysis model, karst zoning, Birjand city,
Abstract :
Exploration of underground water sources as one of the ways to provide drinking water in the world is necessary and inevitable, considering the ever-increasing need of the world for water. The water crisis in Iran, especially in hot and dry cities, has now become a serious challenge. The high salinity and degree of hardness of water and water shortage in Birjand city will remind the officials related to water resources more and more about the solution of alternative water sources in the near future. Knowing the effective factors in finding the potential of karst water resources has a special place in the field of land management studies. The purpose of this research is to find the potential of karst water resources in Birjand city in South Khorasan province using AHP model. In this research, with the hierarchical analysis model, lithology information layers, distance from fault, slope, direction of slope, height, distance from waterway, distance from fault and land use were considered as factor maps. Also, in order to extract the karst water resource potential model, the shape files were called in the GIS10.3 ARC software environment. Different layers of information were classified as criteria maps by applying expert opinions and field observations. Finally, according to the obtained weight, the potential measurement map of karst water sources in Birjand city was obtained. The results showed that of the total area of Birjand city, 26.40% is in the very low developed class, 27.58% in the low developed class, 23.86% in the middle class, 15.75% in the high developed class and 40.40% 6% is too much in the development class. Therefore, according to the judgment of karst experts and experts, the lithology factor of the area with a value of 0.41 has the highest weight and the most important factor controlling the potential of karst water resource development in the studied area, and the land use factor with a value of 0.017 has the lowest weight. assigned The results show that the lithological factors, the distance from the fault and the distance from the waterway have played the most important role in the current development of karst water resources in this city.
1- اشتوکلین(1368): زمینشناسی و زمینساخت فلات ایران، مترجمان: صادق حداد کاوه و حسن حسنعلی زاده، انتشارات وزارت فرهنگ و آموزش عالی، علمی و فرهنگی، 108.
2- بهنیافر، ابوالفضل و قنبر زاده، هادی (1395): ژئومورفولوژیکارست، مشهد، انتشارات نگاران سبز.
3- پور اکبری، سجاد، کلانتری، نصرا...، مصلح، آرش و عقدکی، یاسر(1399): پتانسیلیابی منابع آب کارستی با استفاده از RS ، GIS و AHP (مطالعه موردی: تاقدیسهای لیلی و کینو در شمال شرق خوزستان). مجله مهندسی منابع آب، سال 13، شماره 2، صص112-99.
4- ثقفی، مهدی، امیر احمدی، ابوالقاسم و ربیعی، مریم(1395): تشخیص واکنشهای ژئومورفیک مخروطه افکنه های جنوب دشت بیرجند به فعالیتهای تکتونیکی و تغییرات آب و هوایی با استفاده از دادههای میدانی، نشریه جغرافیا و برنامهریزی، شماره 55، صص 114-97.
5- دشتی برمکی، مجید، رضایی، محسن و اشجاری، جواد(1394): پتانسیلیابی منابع آب کارست کوههای دوان و شاپور بر اساس تصمیمگیری چندمعیاره، مجله پژوهش آب ایران، شماره1، صص100-89.
6- رضایی عارفی، محسن، زنگنه اسدی، محمدعلی، بهنیافر، ابوالفضل و جوانبخت، محمد(1399): پهنهبندی تحول کارست با استفاده از مدل تحلیل سلسله مراتبی (مطالعه موردی:حوضه کوهستانی کلات، خراسان رضوی)، مجله جغرافیا، سال 18، شماره 64، صص94-79
7- زروش، ناهید، واعظی، عبدالرضا و کریمی، حاجی(1393): ارزیابی پتانسیل توسعه کارست در تاقدیس کبیر کوه استان ایلام با استفاده از تلفیق فازی و روش تحلیل سلسله مراتبی AHPو سنجش از دور و GIS، مجله پژوهشهای ژئومورفولوژی کمی، سال 3، شماره 3، صص 144-157.
8- سپند، ساحل، چیتسازان، منوچهر، رنگزن، کاظم و میرزایی، یحیی(1386): تلفیق سنجش از دور و GIS در پتانسیلیابی منابع آب زیرزمینی محدوده لالی، همایش ژئوماتیک تهران.
9- صفاری، امیر، کیانی، طیبه و زنگنه تبار، ساسان(1398): بررسی عوامل مؤثر در توسعهیافتگی و پهنهبندی کارست کوهستان خورین با استفاده از منطق فازی، نشریه تحقیقات کاربردی علوم جغرافیایی، سال 19، شماره 55، صص36-23.
10- علیزاده، امین(1389): اصول هیدرولوژی کاربردی، انتشارات دانشگاه امام رضا، چاپ سیام، مشهد، ایران.
11- قبادی، محمدحسین(1388): زمینشناسی مهندسی کارست، چاپ دوم، انتشارات دانشگاه بوعلی سینا، همدان، ایران.
12- مزیدی، احمد، کرم، امیر و کوراوند، مژگان(1395): پتانسیلیابی توسعه کارست با استفاده از منطق فازی (مطالعه موردی: حوضه دشت سوسن و دشت ایذه)، مجله پژوهشهای ژئومورفولوژی کمی، سال پنجم، شماره 2، صص 141-130.
13 -مکرم، مرضیه و نگهبان، سعید(1398): بررسی و شناسایی مناطق دارای پتانسیل کارستی شدن با استفاده از روش فازی و مدل تحلیل سلسلهمراتبی، مجله جغرافیا و برنامهریزی محیطی، دوره 30، شماره2، 134-121
14- محمودی، فرجا...، (1385): ژئومورفولوژی ساختمانی، چاپ هشتم، انتشارات پیام نور.
15- نخعی، اعظم، زنگنه اسدی، محمدعلی، بهنیافر، ابوالفضل و گلی مختاری، لیلا(1402): پهنهبندی توسعه کارست در حوضه آبریز بقمچ کشف رود با مدل منطق فازی ANP. مجله پژوهشهای ژئومورفولوژی کمی، سال 11، شماره 4، صص255-230
16- واعظی هیر، عبدالرضا، وفادار، مرتضی و آقایی، واحده(1397): پتانسیلیابی منابع آب زیرزمینی موجود در واحدهای کارستی و سازند سخت کوه مورو-صوفیان با بهکارگیری روشهای تصمیمگیری چند معیاره AHP، SAW و F-AHP، مجله فضای جغرافیایی، سال 18، شماره 64، 234-215
17- یمانی، مجتبی، شمسیپور، علیاکبر، جعفری اقدم، مریم و باقری سید شکری، سجاد (1392): بررسی عوامل مؤثر در توسعهیافتگی و پهنهبندی کارست حوضه چله با استفاده از منطق فازی و تحلیل سلسله مراتبی، استان کرمانشاه، مجله علوم زمین سال، شماره 22، صص 66-57.
18- Abdelkareem, M., & Abdalla, F. (2022): Revealing Potential Areas Of Water Resources Using Integrated Remote-Sensing Data And GIS-Based Analytical Hierarchy Process. Geocarto International, 37(25), 8672-8696.
19- Calic, J., (2011): Karstic Uvula Revisited: Toward A Redefinition Of The Term, Geomorphology, Vol.134, No.1, Pp.32-42
20- Ford, D., & Williams, P. D. (2007): Karst Hydrogeology And Geomorphology. West Sussex, England: John Wiley & Sons
21- Hamdani, N., & Baali, A. (2020). Characterization Of Groundwater Potential Zones Using Analytic Hierarchy Process And Integrated Geomatic Techniques In Central Middle Atlas (Morocco). Applied Geomatics, 12(3), 323-335.
22- Jaiswal, R, K. Ghosh, N, C. Galkate, R, V. Thomas, T. 2015. Multi Criteria Decision Analysis (MCDA) For Watershed Prioritization. International Conference In Water Resources, Coastal And Ocean Engineering (ICWECOE 2015). 1553-1560.
23- Milanovic, P. J. (1981): Karst Hydrogeology. Colorado, CO: Water Resources Ublications
24- Lamoreaux, P.E., 2007, Karst: The Foundation For Concepts In Hydrogeology. Bull Eng Geol Environ, 51. Studies, Vol. 64, No.1, Pp. 23-33.
25- Ozcelik, M., & Sarp, G. (2018): Evaluation Of Sustainable Water Supply Alternatives In Karstified Rock Masses Using GIS And AHP Methodology For Antalya (Turkey) Urban Area. Environmental Earth Sciences, 77, 1-14
26- Palmer, A. N., (2007): Cave Geology: Dayton, Ohio, Cave Books, P. 454.22
27- Saaty, T.L., (1980): The Analytic Hierarchy Process, Mcgraw _ Hill, Inc., Reprinted By RWS Publications, Pittsburgh
28- Waele, J., Gutierrez, F., Audra, P., (2015): Karst Geomorphology: From Hydrological Functioning To Palaeoenvironmental Reconstructions Geomorphology, Vol. 229, Pp. 1-2.
29- White, W. B., (1988): Geomorphology And Hydrology Of Karst ،Oxford University Press. Quinlan, J, Groundwater Monitoring In Karst Terrains, EPA. 600/ X.
فصلنامه جغرافیای طبیعی، سال هفدهم، شماره 64، تابستان 1403 17
صص 38-17
پتانسیلیابی منابع آب کارست در شهرستان بیرجند با استفاده از مدل تحلیل سلسله مراتبی و سیستم اطلاعات جغرافیایی
محسن رضایی عارفی1
استادیار گروه آموزش جغرافیا، دانشگاه فرهنگیان، تهران، ایران
محمد معتمدی راد
استادیار گروه آموزش جغرافیا، دانشگاه فرهنگیان، تهران، ایران
ابراهیم تقوی مقدم
استادیار گروه آموزش جغرافیا، دانشگاه فرهنگیان، تهران، ایران
منصور صادقی
دکتری ژئومورفولوژی و مدیریت محیط، دانشگاه خوارزمی، تهران، ایران
تاریخ دریافت:3/2/1403 تاریخ پذیرش:31/5/1403
چکیده
اکتشاف منابع آب زیرزمینی بهعنوان یکی از راههای تأمین آب شرب در جهان با توجه به نیاز روزافزون جهانیان به آب، امری ضروری و اجتناب ناپذیر است. بحران آب در ایران و بخصوص در شهرهای گرم و خشک در حال حاضر به یک چالش جدی تبدیل شده است. شوری و درجه سختی آب بالا و کمبود آب شهرستان بیرجند راهکار منابع آب جایگزین را بیش از پیش به مسئولان مرتبط با منابع آب در آینده نزدیک گوشزد میکند. شناخت عوامل مؤثر در پتانسیلیابی منابع آب کارست در زمینه مطالعات مربوط به مدیریت سرزمین دارای جایگاه ویژهای است. هدف از این پژوهش پتانسیلیابی منابع آب کارست در شهرستان بیرجند در استان خراسان جنوبی با استفاده از مدل AHP میباشد. در این پژوهش با مدل تحلیل سلسله مراتبی، لایههای اطلاعاتی لیتولوژی، فاصله از گسل، شیب، جهت شیب، ارتفاع، فاصله از آبراهه، فاصله از گسل و کاربری اراضی بهعنوان نقشههای عامل در نظر گرفته شدند. همچنین بهمنظور استخراج مدل پتانسیلیابی منابع آب کارست، شیپ فایلها در محیط نرمافزار GIS10.3 ARCفراخوانی شدند. لایههای اطلاعاتی متفاوت با اعمال نظرات کارشناسی و مشاهدات میدانی بهصورت نقشههای معیار طبقهبندی شدند. درنهایت با توجه به وزن به دست آمده، نقشه پتانسیل سنجی منابع آب کارست در شهرستان بیرجند به دست آمد. نتایج حاصله نشان داد که از کل مساحت شهرستان بیرجند، 40/26 درصد در طبقه خیلی کم توسعهیافته، 58/27 درصد در طبقه کم توسعهیافته، 86/23 درصد در طبقه متوسط، 75/15 درصد در طبقه توسعه زیاد و 40/6 درصد در طبقه توسعه خیلی زیاد گرفته است. از اینرو، طبق قضاوت کارشناسان و متخصصان کارست ، عامل لیتولوژی منطقه با ارزش 41/0 بیشترین وزن و مهمترین عامل کنترل کننده پتانسیل توسعه منابع آب کارست در منطقه مورد مطالعه بوده است و عامل کاربری اراضی با ارزش 017/0 کمترین وزن را به خود اختصاص داده است. نتایج نشانگر آن است که به ترتیب عوامل لیتولوژی، فاصله از گسل و فاصله از آبراهه مهمترین نقش را در توسعه کنونی منابع آب کارست در این شهرستان داشتهاند.
واژگان کلیدی: پتانسیلیابی، منابع آب کارست، مدل تحلیل سلسله مراتبی ، پهنهبندی کارست، شهرستان بیرجند
مقدمه
کارست به مجموعهای از فرایندهای زمینشناسی و پدیدههای حاصل از آنها در پوسته و سطح زمین گفته میشود که عمدتاً پدیدههای ناشی از انحلال سنگها است و با تشکیل بازشدگیها، تخریب و تجزیه ساختمان سنگها، ایجاد نوع ویژهای از رژیم گردش آب، نوع خاص توپوگرافی منطقه و بالاخره تشکیل رژیم خاصی از شبکه زهکشی آشکار میگردد(میلانویچ2،1981). اصطلاح کارست برای توصیف نوع خاصی از چشماندازها به کار میرود که دارای اشکال سطحی، غارها و سیستمهای آب زیرزمینی هستند. این چشماندازها و اشکال در سنگهای قابلحل مانند سنگآهک، دولومیت، مرمر، ژیپس و نمک شکل میگیرند(فورد و ویلیامز3،1:2007). کارست، ناحیهای است که هیدرولوژی و مورفولوژی ویژهای دارد که علت این ویژگیها حلالیت بسیار زیاد در آن و توسعه تخلخل ثانویه و نقش عوامل تکتونیکی در توسعه آن است (وال و همکاران4،2:2015). لند فرمهای کارستی عمدتاً در مناطقی با سنگبستر کربناتی قابل حل بهوسیله انحلال شیمیایی آب با اسیدیته کم شکل میگیرند (پالمر5:454:2007). از نظر ژئومورفولوژیست ها به مناطق شگفتانگیز تشکیلشده در سنگهای کربناته، دولومیته و تبخیری در هر ناحیه، پدیدههای کارستی یا مجموعه زمین شکلهای کارست حقیقی میگویند(لامورکس6، 51:2007). عوامل و فرآیندهای مختلفی، تکامل کارست را تحت تأثیر قرار میدهند، توپوگرافی، لیتولوژی و ویژگیهای زمینساختی بهعنوان عوامل منفعل عمل میکنند و فرآیندهای دیگری مانند: تکتونیک فعال، فرآیندهای هیدرولوژیکی وابسته به اقلیم، انحلال زیستی وابسته به محیطزیست و فرآیندهای ژئومورفولوژی، هوازدگی، شیب، آبراههها و یخبندان بهطور فعال در توسعه اشکال کارست نقش دارند(کالیک7،32:2011). بهمنظور توسعه منابع کارستیک عوامل زیادی مؤثر هستند، روی همرفته 8 عنصر لازم برای ایجاد و توسعه کارست، شرایط اقلیمی، توپوگرافی، سنگشناسی، ستبرای لایههای کربناته، کربن، دمای پایین، فشار و موقعیت زمینساخت هستند(وایت، 1988). مناطق کارستی، اهمیت زیادی از دیدگاههای علمی، محیطی و کاربردی دارند. بهطوری که حدود 20 درصد از مساحت خشکیهای جهان را سنگهای کارستی پوشاندهاند و هرگونه توسعه بر روی این اراضی، مستلزم شناخت کامل از ژئومورفولوژی کارست است(بهنیافر و قنبر زاده، 1395: 16). در سال 2013 میلادی 29درصد از جمعیت جهان از منابع آب کارست استفاده میکردند(فورد و ویلیامز،578،2010). در کشور ما حدود 13 درصد مساحت آن را سازندهای کربناته تشکیل میدهند که 90 درصد آن در ارتفاعات زاگرس واقع شده است و در تأمین آب شرب کشور نقش بسیار مهمی دارند(بهنیافر و قنبر زاده، 1395: 16). در مناطق خشک با توجه به اینکه کارستها مناطق مناسبی برای ذخیره سفرههای آب شیرین محسوب میشوند، شناسایی آنها اهمیت دوچندانی مییابد(مکرم و نگهبان،1398). منابع آبی تودههای کارست کوهستانی اهمیت زیادی برای آب شرب شهرها و سکونتگاهها دارند و از این نظر شناسایی تودههای کربناته کارست ساز برای مجریان برنامهریزی آمایش سرزمین بسیار ضروری است. در منطقه تحت مطالعه که اقلیم آن از نوع گرم و خشک است میانگین درجه حرارت منطقه مورد مطالعه 16.4 درجه سانتیگراد و میانگین بارندگی 50 ساله در محل ایستگاه سینوپتیک شهر بیرجند 177 میلیمتر است(ربیعی، 1395). فعالیتهای تکتونیزه در سازندهای کربناته شدید بوده است و تودههای کارستی را در منطقه مورد مطالعه به وجود آورده است. تاکنون تحقیقات تخصصی در زمینه کارست شهرستان بیرجند صورت نگرفته است. فقط بهصورت اندک در شرق شهرستان در منطقه آهنگران تحقیقات اندکی بر روی کارست منطقه صورت گرفته است که در زمینه ژئوتوریسم و اشکال کارستی بوده است. در نتیجه بررسی و کسب شناخت هر چه بیشتر پتانسیلهای موجود در منطقه هم به لحاظ علمی و هم از نظر اقتصادی اهمیت فراوانی دارد. این مورد بهویژه در دو زمینه ایجاد پتانسیلهای آبی کارستهای شهرستان به دلیل بحران آب در استان خراسان جنوبی و پتانسیل اکوتوریستی در مجاورت شهر بیرجند مورد توجه بیشتر است. با توجه به بحران آب در شهرستان بیرجند و بهرهبرداری بیش از اندازه از منابع آبخانههای آبرفتی موجود، شناخت و بهرهبرداری از منابع جدید یکی از اولویتهای تأمین آب در منطقه است که کارستها برای این منظور بسیار مهماند. در ایـران پهنهبندی اشکال کارستی را محققان مختلفی بررسی کردهاند. دشتی برمکی و همکاران در سال 1394 در پژوهشی با عنوان پتانسیلیابی منابع آب کارست کوههای دوان و شاپور بر اساس تصمیمگیری چندمعیاره به این نتیجه رسیدند که پتانسیل آب کارست دشتک شمالی بیش از جنوبی است. هماهنگی نواحی با پتانسیل بالا با چشمههای با آبدهی زیاد نشان میدهد که مدل ارائه شده میتواند بهطور موفقیتآمیزی برای منطقههای کارستی استفاده شود. مزیدی و همکاران در سال 1395 در مقالهای با عنوان پتانسیلیابی توسعه کارست با استفاده از منطق فازی ( مطالعه موردی: حوضه دشت سوسن و دشت ایذه) به این نتیجه رسیدند که تلفیق لایههای موجود نشان داد که مناطق مستعد به کارستی شدن با مساحت تقریباً 32 درصد در جنوب، شرق، غرب و بخشهایی از شمال شرق منطقه واقعشدهاند. واحد آهک آسماری، گچساران و ایلام ـ سروک شامل این بخش میباشند که از رسوبات کاملاً آهکی تشکیلشدهاند و ازنظر تأمین آب مهمترین سفرههای آبدار کارستی در منطقه مطالعاتی محسوب میشوند. وجود اشکال سطحی کارست از قبیل کارن، ریل کارن و درز و شکافها در سازندهای مزبور نقش مؤثری در پراکندگی چشمهها و آبدهی آنها و کارستی شدن منطقه دارد. همچنین قسمتهای مرکز و بخشهایی از شمال شرق محدوده مطالعاتی با مساحت 50 درصد از پتانسیل بسیار کمی جهت کارستی شدن برخوردار است.
واعظی و همکاران در سال 1397 در پژوهشی با عنوان پتانسیلیابی منابع آب زیرزمینی موجود در واحدهای کارستی و سازند سخت کوه مورو-صوفیان با بهکارگیری روشهای تصمیمگیری چند معیاره AHP، SAW و F-AHP به این نتیجه رسیدند که با روش میانگینگیری وزنی درجهای و با استفاده از تابع حسابگر (Raster Calculator) نقشه پهنهبندی پتانسیل آب زیرزمینی در پنج گروه خیلی کم، کم، متوسط، زیاد و خیلی زیاد طبقهبندی کردند و با موقعیت چشمهها و قناتها اعتبارسنجی کردند. اعتبارسنجی نقشه پهنهبندی به دست آمده با موقعیت چشمههای منطقه نشان داد که مدل F-AHP و SAW نتایج مشابهی را ارائه میدهد. نتیجه این پژوهش نشان دادبخش مرکزی تاقدیس مورو پتانسیل بالایی از آب زیرزمینی را بهویژه در یال شمالی دارد. تراکم شکستگی و نوع لیتولوژی منطقه بیشترین تأثیر را بر تشکیل منابع آب زیرزمینی منطقه داشته است.
مکرم و نگهبان در سال 1398 در مقالهای با عنوان بررسی و شناسایی مناطق دارای پتانسیل کارستی شدن با استفاده از روش فازی و مدل تحلیل سلسلهمراتبی در بخشهایی از غرب استان فارس به این نتیجه رسیدند که در شمال غربی و قسمتی از جنوب شرق استان فارس احتمال کارست با توجه به نقشهها صد درصد است؛ در حالی که مناطق مرکزی کارست کمتری دارند و قسمتهایی از حاشیه منطقه در شرق و غرب بدون کارست هستند. بهمنظور بررسی و تعیین صحت مدل از بعضی نقاط نمونه استفاده و با توجه به نقاط نمونه مشخص شد مدل دقت زیادی برای تعیین مناطق کارستی دارد. بنابراین از روش فازی میتوان بهمثابه روشی دقیق برای بررسی پتانسیل کارستی مناطق مختلف استفاده کرد.
رضایی عارفی و همکاران در سال 1399 در پژوهشی با عنوان پهنهبندی تحول کارست با استفاده از مدل تحلیل سلسله مراتبی مطالعه موردی حوضه کوهستانی کلات در خراسان رضوی به این نتیجه رسیدند که از کل مساحت حوضه کلات، 04/19 درصد در طبقه کمتر توسعهیافته، 57/24 درصد در طبقه توسعهنیافته، 01/42درصد در طبقه متوسط و 38/14 درصد در طبقه توسعهیافته قرار گرفته است. از اینرو، عامل لیتولوژی منطقه با ارزش 53/0 بیشترین وزن و مهمترین عامل کنترلکننده پتانسیل توسعه کارست در منطقه مورد مطالعه بوده است و عامل کاربری اراضی با ارزش 01/0 کمترین وزن را به خود اختصاص داده است که کمترین تأثیر را در کارست زایی کنونی حوضه به خود اختصاص داده است. نتایج نشانگر آن است که به ترتیب عوامل زمینشناسی، ارتفاع، توپوگرافی و فاصله از آبراهه مهمترین نقش را در توسعه کنونی کارست در این منطقه داشتهاند.
پور اکبری و همکاران در سال 1399 در مقالهای با عنوان پتانسیلیابی منابع آب کارستی با استفاده ازRS ، GIS و AHP (مطالعه موردی: تاقدیس-های لیلی و کینو در شمال شرق خوزستان) به این نتیجه رسیدند که بخش مرکزی تاقدیس از بیشترین پتانسیل برخوردار است. از طرف دیگر دماغه شمال غربی (سازند گرو) دارای کمترین پتانسیل است. چشمهها و دولینهای موجود، در محدوده با پتانسیل بالا قرار گرفتهاند.
نخعی و همکاران در سال 1402 در مقالهای با عنوان پهنهبندی توسعه کارست در حوضه آبریز بقمچ کشف رود با مدل منطق فازی ANP به این نتیجه رسیدند که 16.3 درصد از مساحت منطقه در پهنه با پتانسیل توسعه بسیار زیاد، 33.7 درصد با توسعه زیاد، 6.6 درصد با توسعه متوسط ، 3.6 درصد با توسعه کم و 39.8 درصد فاقد کارست است. عامل لیتولوژی منطقه با ارزش 0.233 بیشترین وزن و مهمترین عامل کنترل کننده پتانسیل توسعه کارست در منطقه مورد مطالعه بوده است و عامل فاصله از آبراهه و کاربری اراضی کمترین وزن را به خود اختصاص داده و کمترین تأثیر را در کارست زایی کنونی حوضه به خود اختصاص داده است. نتایج نشانگر آن است که به ترتیب عوامل لیتولوژی، بارش، دما و ارتفاع مهمترین نقش را در توسعه کنونی کارست در این منطقه داشتهاند.
با توجه به اهمیت مناطق کارستی، پژوهشهای به نسبت جامعی در جهان در رابطه با توسعه کارست صورت گرفته است که در این میان میتوان به موارد زیر اشاره کرد:
از کلیک8 و همکاران در سال 2018 در پژوهشی با عنوان ارزیابی جایگزینهای تأمین آب پایدار در تودههای سنگی کارستی شده با استفاده از روش GIS و AHP برای منطقه شهری آنتالیا (ترکیه) به این نتیجه دست یافتند که انتظار میرود در آینده نزدیک تقاضا برای آب با افزایش جمعیت و عملیات کشاورزی در منطقه شهری آنتالیا به میزان قابلتوجهی افزایش یابد. به همین دلیل، منابع آب جایگزین برای آینده آنتالیا مهم هستند. این مطالعه یک روش مفید برای یافتن جایگزینهای تأمین آب برای مناطق شهری مانند آنتالیا ارائه میکند. در این مطالعه، تودههای سنگ کارستی شده برای هیدرو سیستمهای کارستی و منبع آب منافذ سنگ با استفاده از سیستمهای اطلاعات جغرافیایی و فرآیند تحلیل سلسله مراتبی مورد ارزیابی قرار گرفتهاند. علاوه بر این، این مطالعه ارزش عملی تأمین آب از منافذ سنگ و همچنین امکان کاربرد و انتخاب مناطق بالقوه را مورد بحث قرار میدهد.
همدانی و بعالی9 در سال 2020 در مقاله خود با عنوان شناسایی مناطق بالقوه آب زیرزمینی با استفاده از فرآیند تحلیل سلسله مراتبی و تکنیکهای ژئوماتیک یکپارچه در مرکز اطلس میانه (مراکش) به این نتیجه رسیدند که فرآیند تحلیل سلسله مراتبی (AHP) ابزاری ارزشمند است که توسط متخصصان مختلف برای ارتقای ظرفیت مدیریتی خود به کار گرفته شده است. در این پژوهش، یک رویکرد سازگار برای شناسایی مناطق بالقوه آب زیرزمینی (GWP) از طریق استفاده از تکنیکهای AHP، سیستم اطلاعات جغرافیایی (GIS) و (RS) در مراکش استفاده شد. نتایج نشان داد که تنها 3.88 درصد از منطقه مورد مطالعه دارای پتانسیل آب زیرزمینی بسیار خوب و 17.22 درصد دارای پتانسیل خوب، علاوه بر 20.20 درصد با پتانسیل متوسط است، در حالی که در بیشتر مناطق منطقه، 29.89 درصد دارای پتانسیل آب زیرزمینی ضعیف است و 18.60 درصد در پتانسیل آب زیرزمینی بسیار ضعیف و 10.49 درصد با پتانسیل کارستی ندارند.
عبدالکریم و عبدا... در سال 2022 در مقاله با عنوان آشکارسازی مناطق بالقوه منابع آب با استفاده از دادههای سنجش از دور یکپارچه و فرآیند تحلیل سلسله مراتبی مبتنی بر GIS به این نتیجه رسیدند که بخشی از وادی الصحبا عربستان سعودی برای آشکارسازی و مدلسازی مساحت بالقوه آبهای زیرزمینی با استفاده از دادههای یکپارچه سنجش از دور، زمینشناسی و هیدرولوژی از طریق روش فرآیند تحلیل سلسله مراتبی مبتنی بر GIS آزمایش میشود. نقشه اکتشاف خروجی نشان داد که پنج منطقه پیشبینی آب زیرزمینی (GPZs) بسیار کم (11.55%)، کم (28.18%)، متوسط (29.14%)، زیاد (22.41%) و بسیار زیاد (8.73%) هستند. بیشترین مناطق پیش بینی شده بخش وسیعی از الخرج و پاییندست وادی الصهبه را پوشش میدهند. از آنجایی که منطقه تحت مطالعه در نزدیکی شهر بیرجند قرار دارد، لذا بهرهبرداری از آب مناطق کارستی میتواند در تأمین آب شرب شهرستان بیرجند اهمیت زیادی داشته باشد. بنابراین هدف از پژوهش حاضر پتانسیلیابی منابع آب کارست در شهرستان بیرجند در استان خراسان جنوبی با استفاده از مدل AHP است.
موقعیت منطقه مورد مطالعه
محدوده مورد مطالعه، شهرستان بیرجند در استان خراسان جنوبی است. شهرستان بیرجند از نظر موقع ریاضی بین 32 درجه و 35 دقیقه تا 33 درجه و 8 دقیقه عرض جغرافیایی شمالی و 58 درجه و 41 دقیقه تا 59 درجه و 46 دقیقه طول جغرافیایی شرقی قرار گرفته است. این شهرستان با وسعت 3937.5 کیلومترمربع در شرق ایران قرار گرفته است (شکل 1) . بالاترین نقطه ارتفاعی طبق نقشه توپوگرافی 2680 متر در ارتفاعات باغران، حداقل ارتفاع آن 1360 متر در غرب شهرستان و میانگین ارتفاع 1884 متر است. متوسط دمای سالیانه در ارتفاعات 11.4 درجه و در دشت 13.8 درجه سانتیگراد است. متوسط بارندگی سالانه ۱۶۷ میلیمتر است که بیشترین میزان آن مربوط به فصول زمستان و بهار است . بیرجند از شمال به شهرستان قاین، از جنوب به شهرستان نهبندان و کرمان، از شرق به افغانستان و از غرب به شهرستانهای فردوس و طبس محدود میشود. شهرستان بیرجند دارای 4 شهر، 4 بخش و 19 دهستان است. شهرستان بیرجند، از قسمتهای مختلف، شامل جلگه، دشت، و ارتفاعات تشکیلشده است. بیرجند با توجه به وضعیت پستی و بلندیها دارای: 1- اقلیم خشک و گرم شامل مناطق دشت و حاشیه شهرستان 2- اقلیم خشک ملایم که در بخشهای مرتفع شمالی و مرکزی شهرستان یعنی نقاط مرتفع کوههای »باغران« و »مؤمن آباد« مشاهده میشود. لذا بر اساس تقسیمبندی »اقلیمی«، بیرجند جزء مناطق خشک محسوب میشود که دارای زمستانهای سرد و تابستانهای خشک و گرم است. میزان بارش در این شهر با توجه به آب و هوای آن، کم بوده و بیشترین میزان آن، از آذر تا اردیبهشت رخ میدهد که در فصل زمستان اغلب بهصورت بارش برف است. تفاوت نسبتاً زیاد درجه حرارت روز و شب و همچنین اختلاف دمای زمستان و تابستان در این شهرستان به علت نزدیکی آن به مناطق کویری است. فصل گرم بیرجند معمولاً طولانی است و ماههای خرداد و تیر و مرداد را در بر میگیرد و ایام سرد سال شامل ماههای آذر و دی و بهمن است. بیرجند، بهطور میانگین در 72 روز از سال دارای دمای زیر صفر درجه و در ۱۴۲ روز از سال دارای دمای بالای ۳۰ درجه سانتیگراد است.
شکل 1: نقشه موقعیت شهرستان بیرجند
دادهها و روشها
این پژوهش با استفاده از روشهای میدانی، کتابخانهای و اسنادی صورت گرفته است. در ابتدا با استفاده از Dem 30 متر، محدوده مورد مطالعه تعیین گردید. با استفاده از نقشههای توپوگرافی 1:50000 بیرجند و زمینشناسی 1:100000 بیرجند جهت استفاده در شناسایی لند فرمها و سازندهای کربناته استفاده گردید. در این پژوهش با توجه به منابع موجود(فایل رقومی ارتفاع) و تهیه نقشههای مربوطه و بر اساس درجه اهمیت پارامترهای تأثیرگذار در توسعه منابع آب کارست در شهرستان بیرجند فاکتورهای مهم انتخاب شدند. با توجه به نقش عوامل متعدد و توسعه کارستیفیکاسیون درگذشته منطقه و نقش عوامل مختلف در توسـعه کارسـت، نه لایه اطلاعاتی انتخاب شـدند. با استناد به مطالعات میدانی و نظر کارشناسان و مشاوران ، عامل لیتولوژی بهعنوان مهمترین پارامتر در توسعه منابع آب کارست شهرستان بیرجند انتخاب شد و بعد از آن فاصله از گسل با توجه به نقش مهم آن در توسعه درز و شکاف انتخاب گردید. بعد از آن لایه فاصله از آبراهه انتخاب شدند. کمترین اهمیت به لایه کاربری اراضی داده شد. بر اساس منابع مطالعاتی موجود از بین عوامل دخیل در تحول کارســت، در این پژوهش لایههای اطلاعاتی لیتولوژی، فاصله از گسل، اقلیم (بارش و دما)، ارتفاع، شیب، فاصله از آبراهه، جهت شیب و کاربری اراضی بهعنوان نقشههای عامل در نظر گرفته شدند و بهمنظور استخراج پتانسیلیابی منابع آب کارست، لایههای فوق در محیط GIS فراخوانی شد. لایههای اطلاعاتی مختلف با اعمال نظرات کارشناسی و اختصـاص وزن به هر لایه و بازدیدهای میدانی بهصورت نقشههای معیار طبقهبندی شدند و با توجه به وزن به دست آمده، نقشه پتانسیلیابی منابع آب کارست منطقه مورد مطالعه به دست آمد. برای تهیه مدل پتانسیل سنجی منابع آب کارست، لایههای اطلاعاتی فوق به نرمافزار Arc GIS 10.3 معرفی شدند. لایههای اطلاعاتی فوق بهصورت نقشههای معیار طبقهبندی شدهاند و با توجه به درجه اهمیت هرکدام از پارامترها، رتبهای از 1 تا 9 بر اساس روش توماس آل ساعتی10 (ساعتی، 1980:1) به آنها داده شد. انتخاب این ارزشها بر اساس اعمال نظر کارشناسی و شناخت منطقه بوده است. پس از ارزش دهی به هر یک از لایهها بر اساس روش کارشناسی و روش AHP وزن مناسبی که مشخص کننده تأثیر هر عامل در درجه توسعه منابع آب کارست است، به آنها داده شد. لازم به ذکر است که جمع وزنها عدد 1 در نظر گرفته شد. سپس از نرمافزار EXPERT choice که بر اســاس روش تحلیل سلسله مراتبی طراحی شده است، جهت وزن دهی پارامترها استفاده شد. بعد از آن نقشه رستری نه پارامتر مذکور در محیط GIS تهیه و هر یک از این نقشههای رستری طبقهبندی شدند.
مدل تحلیل سلسله مراتبی (AHP)
AHP یک نظریه کلی برای اندازهگیری است. برای استخراج اولویتهای نسبی در مقیاسهای مطلق (غیر متغیر تحت تغییر هویت) از مقایسههای زوجی گسسته و پیوسته در ساختارهای سلسله مراتبی چند سطحی استفاده میشود. این مقایسهها ممکن است از اندازهگیریهای واقعی یا از یک مقیاس بنیادی گرفته شود که قدرت نسبی ترجیحات و احساسات را منعکس میکند. AHP توجه ویژهای به خروج از سازگاری و اندازهگیری این انحراف و وابستگی درون و بین گروههای عناصر ساختار خود دارد. این مدل گستردهترین کاربردهای خود را در تصمیمگیری چند معیاره در برنامهریزی و تخصیص منابع و در حل تعارض یافته است. در شکل کلی آن، AHP یک چارچوب غیرخطی برای انجام تفکر قیاسی و استقرایی بدون استفاده از قیاس است. این امر با در نظر گرفتن چندین عامل بهطور همزمان امکانپذیر میشود و امکان وابستگی و بازخورد فراهم میشود و مبادلات عددی برای رسیدن به یک ترکیب یا نتیجهگیری امکانپذیر میشود ( ساعتی و وارگاس، 2006).
تحلیل سلسله مراتبی AHP
برای انجام تحلیلهای مربوط به روش AHP از نرمافزار اکسپرت چویس استفاده گردید. مراحل مختلف این روش بهصورت زیر است:
- مقایسه زوجی
اولین قدم در این روش مشخص نمودن معیارها برای انجام تحلیل میباشد. بعد از انتخاب معیارها مرحله بعد، مقایسه زوجی معیارها بوده که این امر با استفاده از نظر متخصصان و کارشناسان با توجه به جدول ارائه شده توسط ساعتی انجام گردید. در این مرحله از متخصصان خواسته شد که ارجحیت معیارها را با استفاده از مقیاس ساعتی از عدد 1 (کم اهمیت) تا عدد 9 بیشترین اهمیت مشخص نمایند. بعد از تعیین ارجحیت معیارها برای اطمینان از صحت تحلیلها نرخ سازگاری مقایسات محاسبه گردید. نتیجه محاسبه نرخ سازگاری نشاندهنده عدد کمتر از 1/0 بود که نشان از صحت مقایسات بود(جدول1).
جدول 1: مقایسات زوجی در فرآیند تحلیل سلسله مراتبی به روش (ساعتی، 1980)
مقیاس | قضاوت |
1 | برابر |
2 | بین برابر تا متوسط |
3 | متوسط |
4 | بین متوسط تا زیاد |
5 | زیاد |
6 | بین زیاد تا خیلی زیاد |
7 | خیلی زیاد |
8 | بین خیلی زیاد تا فوقالعاده زیاد |
9 | فوقالعاده زیاد |
بررسی نرخ سازگاری
سازگاری قضاوت ذهنی را میتوان با تخمین نسبت سازگاری که مقایسه بین شاخص ثبات و شاخص سازگاری تصادفی است بررسی کرد (جیسوال11 و همکاران، 2015). برای بررسی نرخ سازگاری از رابطه زیر استفاده گردید.
CR =
در این رابطه CR نرخ سازگاری، CI شاخص ثبات است RI شاخص سازگاری تصادفی است. شاخص سازگاری، معیاری است برای سنجش سازگاری که میتوان با استفاده از معادله زیر تخمین زد:
CI =
در این رابطه ʎmax مقدار ویژه اصلی به دست آمده از ماتریس اولویت است و n اندازه ماتریس مقایسه است.
یافتههای پژوهش
نتایج مدل AHP
در ابتدا پرسشنامه مقایسه زوجی پارامترهای مؤثر بر پتانسیلیابی منابع آب شهرستان بیرجند در اختیار متخصصان کارست قرار گرفت. بعد از انجام مرحله مقایسه زوجی و مقایسه دو به دوی معیارها با یکدیگر ، وزنهای مربوطه وارد نرمافزار اکسپرت چویس (Expert Choice) شد و وزن هر یک از معیارها مشخص گردید (جدول 2)(جدول3). نمودار ضرایب حاصل از مقایسه زوجی معیارهای اصلی مؤثر در پتانسیلیابی منابع آب کارست طبق فرایند تحلیل سلسله مراتبی در نرمافزار اکسپرت چویس با مقدار شاخص ناسازگاری 08/0محاسبه گردید(شکل2). بعد از تعیین وزن هر معیار با استفاده از نرمافزار ArcGis10.3 وزن هر معیار در لایه مربوطه ضرب گردید و لایههای مربوط به هر معیار با یکدیگر ادغام شده و نقشه پتانسیل توسعه منابع آب کارست بر اساس مدل AHP تهیه گردید.
جدول2: وزن دهی عوامل مؤثر بر پتانسیلیابی منابع آب کارست در شهرستان بیرجند
ردیف | نام | لیتولوژی | فاصله از گسل | فاصله از آبراهه | شیب | جهت شیب | ارتفاع | بارش | دما | کاربری اراضی |
1 | لیتولوژی | 1 | 5 | 6 | 7 | 9 | 7 | 6 | 8 | 9 |
2 | فاصله از گسل |
| 1 | 2 | 4 | 6 | 5 | 6 | 7 | 8 |
3 | فاصله از آبراهه |
|
| 1 | 3 | 5 | 4 | 4 | 5 | 5 |
4 | شیب |
|
|
| 1 | 3 | 4 | 2 | 3 | 5 |
5 | جهت شیب |
|
|
|
| 1 | 3/1 | 2/1 | 2/1 | 3 |
6 | ارتفاع |
|
|
|
|
| 1 | 2/1 | 2 | 4 |
7 | بارش |
|
|
|
|
|
| 1 | 3 | 5 |
8 | دما |
|
|
|
|
|
|
| 1 | 3 |
9 | کاربری اراضی |
|
|
|
|
|
|
|
| 1 |
جدول3: وزن لایههای اطلاعاتی به روش AHP و اعمال نظر کارشناسی در منطقه مورد مطالعه
نام لایه | وزن نهایی |
لیتولوژی | 41/. |
فاصله از گسل | 19/. |
فاصله از رودخانه | 13/. |
شیب | 80./. |
بارش | 56./. |
ارتفاع | 44./. |
دما | 32./. |
جهت شیب | 27./. |
کاربری اراضی | 17./. |
جمع | 1 |
شکل2: نمودار ضرایب حاصل از مقایسه زوجی معیارهای اصلی مؤثر در پتانسیلیابی منابع آب کارست طبق فرایند تحلیل سلسله مراتبی
پارامترهای مؤثر در پتانسیلیابی منابع آب کارست
عامل لیتولوژی
به علت تأثیرگذاری شدید نوع سازند بر نفوذ آب و قابلیت متفاوت حفظ آب در سازندهای مختلف این نمایه مد نظر قرار گرفته است(پور اکبری و همکاران، 1399). خصوصیات مختلفی از قبیل میزان رخنمون، ضخامت و خصوصیات لیتولوژیکی سنگهای کارستی و ارتباط با لیتولوژی های دیگر در مقیاس ناحیهای، کنترل کنندههای زمینشناسی توسـعهٔ کارسـت را تشـکیل میدهند (رضایی و همکاران،1399). منطقه مورد مطالعه در بخشی از فلات ایران که از نقطه نظر زمینشناسی منطقه مورد مطالعه در محدوده زیر ناحیه نهبندان-خاش یا فلیش قرار دارد. بهطور کلی در نتیجه تحولات زمینساختی در این منطقه، انواع رخسارههای آذرین رسوبی و دگرگونی دوران دوم و سوم زمینشناسی گسترش یافته است. در دوران کواترنری نیز نهشتههای آبرفتی ناشی از فرایندهای مورفوکلیما، مورفوژنز و مورفوتکتونیک منطقه در قسمتهای مختلف بهویژه در پهنههای مشرف به ارتفاعات گسترش دارد. بهطور عمده سازندهای متعلق با دوران کربونیفر تا ترشیاری در منطقه مورد مطالعه مشاهده میشود. مجموعه سنگهای کوهستان باقران در جنوب شهر بیرجند بهعنوان یک مجموعه آمیزه رنگین معرفی شده است(اشتوکلین، 1368). قدیمیترین رخنمونهای این واحدها به دوره پرمین مربوط به سنگهای پریدوتیت است و وسیعترین گسترش رخنمونها متعلق به رسوبات سطح پایین مخروطه افکنه و رسوبات تراسهای دره با 574.8 کیلومتر مربع و توف آندزیتی داسیتی با 522.6 کیلومتر مربع است. سازندهای آهکی عمده منطقه مورد مطالعه، آهک نومولیتی، سنگآهک، آهک پلاژیک و رسوبات فلیش از نوع آهک میباشد که در شمال غربی و جنوب منطقه مورد مطالعه گسترش دارد. در منطقه مورد مطالعه در مجموع 36 واحد سنگی وجود دارد که بر روی نقشه زمینشناسی(شکل3) مشخص شده است. مطالعات زمینشناسی در مناطق کارستی، با توجه به اهمیت خیلی زیاد نوع لیتولوژی در توسعه منابع آب کارست و کارست زایی اهمیت زیادی دارد. وجود سنگهای انحلال پذیر یکی شرایط اصلی در توسعه پهنههای کارسـتی میباشد. طبق طبقهبندی، آهک نومولیتی، سنگآهک و گل سنگآهکی برای پتانسیل کارست زایی بیشترین وزن را به خود اختصاص داده است. در جدول 4، توزیع مشخصات سنگشناسی منطقه مورد مطالعه همراه با سن، مساحت به درصد و مساحت به کیلومتر مربع آورده شده است.
جدول 4: مشخصات سنگشناسی منطقه مورد مطالعه
ردیف | واحد سنگی | مشخصه به زبان انگلیسی | مشخصه به فارسی | سن | مساحت به درصد | مساحت به کیلومتر مربع |
1 | E1L | Limestone nummolitic | آهک نومولیتی | ائوسن | 1.91 | 75.5 |
2 | E1m | Marl, gypsiferous marl and limestone | مارن، مارن گچی و سنگآهک | اوایل ائوسن | 27/0 | 10.7 |
3 | E2c | Conglomerate and sandstone | کنگلومرا و ماسهسنگ | اواسط ائوسن | 32/0 | 12.7 |
4 | E2s | Sandstone, marl and limestone | ماسهسنگ، مارن، و سنگآهک | ائوسن | 4.39 | 172.8 |
5 | E2sht | Tuffaceous shale and tuff | شیل و توف | ائوسن | 0 | 0.101091 |
6 | Eavt | Andesitic volcanic tuff | توف آتشفشانی آندزیتی | ائوسن | 3.36 | 132.3 |
7 | Ed.avs | Dacitic to Andesitic volcanosediment | رسوبات آتشفشان داسیتی تا آندزیتی | ائوسن | 1.74 | 68.7 |
8 | Ktzl | Thick bedded to massive, white to pinkish orbitolina bearing limestone ( TIZKUH FM ) | سنگآهک دارای بستر ضخیم تا حجیم، سفید تا صورتی مایل به اربیتولین | اوایل کرتاسه | 07/0 | 6.6 |
9 | Kuf | Flysch type sediments including shale, sandstone, limestone and conglomerate | رسوبات نوع فلیش شامل شیل، ماسهسنگ، سنگآهک و کنگلومرا | اواخر کرتاسه | 6.1 | 242.3 |
10 | Kus | Flysch turbidite sandsone with interbedd calcareous mudstone and shale | ماسهسنگ توربیدیت فلیش با گل سنگآهکی و شیل بین لایهای | اواخر کرتاسه | 1.88 | 76.4 |
11 | L.E.Oa.bv | Andesitic to basaltic volcanic | آتشفشانی آندزیتی تا بازالتی | کرتاسه-پالئوسن | 7/0 | 31.4 |
12 | L.E.O.av | Andesitic volcanic | آتشفشان آندزیتی | اواخر ائوسن-الیگوسن | 1.55 | 61.3 |
13 | L.E.Od.avb | Dacitic to andesitic volcanobreccia | برش آتشفشانی داسیتی تا آندزیتی | ائوسن-الیگوسن | 07/0 | 6.7 |
14 | Murm | Ligth - red to brown marl and gypsiferous marl with sandstone intercalations | مارن روشن - قرمز تا قهوهای و مارن گچی با ماسهسنگهای درونآمیزی | میوسن | 0.09 | 3.6 |
15 | OMa.bv | Andesite and andesitic lava flow | جریان گدازه آندزیتی و آندزیتی | الیگوسن-میوسن | 5.38 | 211.9 |
16 | OMat | Andesitic tuff | توف آندزیتی | الیگوسن | 2.31 | 91.3 |
17 | OMd.at | Dacitic Andesitic tuff | توف آندزیتی داسیتی | الیگوسن-میوسن | 5.6 | 522.6 |
18 | OMd.av | Dacitic Andesitic volcanic | داسیتی آندزیتی آتشفشانی | الیگوسن-میوسن | 9.1 | 360.9 |
19 | OMgr-di | Granite to diorite | گرانیت به دیوریت | الیگوسن-میوسن | 0.34 | 13.6 |
20 | Oav | Oligocene andesitic lava flows | جریان گدازه الیگوسن آندزیتی | الیگوسن | 1.03 | 40.9 |
21 | Od.av | Dacitic to andesitic volcanic | آتشفشانی داسیتی تا آندزیتی | الیگوسن | 0.32 | 12.8 |
22 | Pd | Peridotite including harzburgite, dunite, lerzolite and websterite | پریدوتیت شامل هارزبورگیت، دونیت، لرزولیت و وبستریت | ترشیاری-کرتاسه | 1.28 | 50.7 |
23 | PeEf | Flysch turbidite, sandstone and calcareous mudstone | فلیش توربیدیت، ماسهسنگ و گل سنگآهکی | پالئوسن-ائوسن | 17.9 | 704.9 |
24 | PeEm | Marl and gypsiferous marl locally gypsiferous mudstone | مارن و مارن گچی گلسنگ محلی گچی | پالئوسن-ائوسن | 0.40 | 15.9 |
25 | PeEph | Phyllite | فیلیت | پالئوسن-ائوسن | 0.91 | 35.88 |
26 | PlQbv | Basaltic volcanics | آتشفشان های بازالتی | پالئوسن-کواترنری | 86/0 | 33.98 |
27 | Plc | Polymictic conglomerate and sandstone | کنگلومرا پلیمیکتیک و ماسهسنگ | پلیوسن | 8.1 | 319.5 |
28 | Pld.asv | Dacitic to andesitic subvolcanic rocks | سنگهای ساب آتشفشانی داسیتی تا آندزیتی | پلیوسن | 44/0 | 17.6 |
29 | Qbv | Olivine basalt and basalt related to Bazman Volcanism and partly related to Taftan Volcanism | بازالت الیوین و بازالت مربوط به آتشفشان بزمان و تا حدودی مربوط به آتشفشان تفتان | کواترنری | 2/0 | 8.2 |
30 | Qcf | Clay flat | رس | کواترنری | 1.68 | 66.2 |
31 | Qft1 | High level piedmont fan and vally terrace deposits | مخروطه افکنه و تراسهای آبرفتی | کواترنری | 4.1 | 164.4 |
32 | Qft2 | Low level piedment fan and vally terrace deposits | رسوبات سطح پایین مخروطه افکنه و رسوبات تراسهای دره | کواترنری | 14.6 | 574.8 |
33 | db | Diabase | دیاباز | کرتاسه پایانی | 0.02 | 1.1 |
34 | om1 | Tectonized association of peridotites, gabbro, diorite, trondhjemite, diabase and basic volcanics ( Ophiolite Melange ) | ارتباط تکتونیزه پریدوتیت ها، گابرو، دیوریت، ترونجمیت، دیاباز و آتشفشان های پایه (افیولیت ملانژ) | ترشیاری-کرتاسه | 85/0 | 33.6 |
35 | Sp | Spilitic rocks locally with pillow structure | سنگهای اسپیلیتی بهصورت محلی با ساختار بالشی | کرتاسه | 1.44 | 57 |
36 | tm | Tectonic melange - association of ophiolitic components, pelagic limestone, radiolarian chert and shale with or without Eocene sedimentary rocks (Coloured Melange complex ) | ملانژ تکتونیکی - ترکیب اجزای افیولیتی، سنگآهک پلاژیک، چرت رادیولاری و شیل با یا بدون سنگهای رسوبی ائوسن (مجموعه ملانژ رنگی) | تریاس-پالئوسن | 1.5 | 59 |
عامل فاصله از گسل
براثر عملکرد فرآیندهای تکتونیکی، درز و شکافها در جهتهای عمودی، افقی و مایل در تودههای سنگی ایجاد و گسترش مییابند. درزه و شکافهایی که بهصورت یک شبکه گسترش مییابند در شکلگیری و تشکیل اشکال زیرزمینی مانند مجاری به همپیوسته و تشکیل آبخوانهای کارستی نقش بسیارمهمی بازی میکنند. پارامتر فاصله از گسل رابطه معکوس با توسعه کارست دارد. زیرا در مناطق نزدیک گسل به علت تکتونیزه بودن و وجود درزه و شکاف، میزان نفوذپذیری سنگها بیشتر بوده و در نتیجه میزان کارست زایی افزایش مییابد(صفاری و همکاران، 1398). علت اصلی نفوذپذیری سنگهای کربناته، وجود درزهها است که امکان ورود آب و جریان یافتن آن را به داخل توده سنگ فراهم میکند. (قبادی، 1388). شکستگیها و ساختارهای تکتونیکی که اصطلاحاً خطواره نیز نامیده میشوند، به دلیل ایجاد فضاهایی در سازندها و واحدهای زمینشناسی جهت عبور آب و حرکت آن به نقاط پایینتر درون زمین نقاط ضــعیفی تلقی میشوند که اهمیت آنها در سازندهای سخت و آهکی بیشتر است(محمودی، 1385). بنابراین عامل حریم گسل بهعنوان پارامتر مثبتی جهت پتانسیلیابی توسعه منابع آب کارست مد نظر قرار میگیرد. منطقه مورد مطالعه در یکی از نواحی نسبتاً فعال تکتونیکی شرق کشور قرار گرفته که بهتبع آن گسلهایی با امتداد تقریبی شمال غربی- جنوب شرقی به وجود آمدهاند. در وزن دهی به لایه گسل، فرض اصلی بر این است که میزان توسعهیافتگی کارست در مناطق نزدیک به گسل بیشتر است و وزن بیشتری را به خود اختصاص میدهد و مناطق دور از خطوط گسل وزن کمتری را به خود اختصاص داده است. فاصله از گسل در پنج طبقه برای محدوده مورد مطالعه تهیه گردید (شکل4) که تا حریم 500 متری بیشترین وزن را به خود اختصاص داد یعنی بیشتر مستعد توسعه منابع آب کارست است و حریم بیشتر از 200000 متر کمترین وزن را به خود اختصاص داد، یعنی کمترین پتانسیل را در منطقه برای توسعه منابع آب کارست دارا میباشد.
عامل بارش
وجود آب، فاکتور اصلی اقلیمی در توسعهٔ کارست است، این عامل اصلیترین متغیر در کنترل انحلال و فرسایش بوده که بهصورت طبیعی، کارست در مناطقی پیشرفت میکند که میزان بارندگی بالاتر باشد و مناطق خشک یا بسیار سرد مانع از توسعه کارست میگردند(صفاری و همکاران،1397). یکی از مهمترین پارامترهای مؤثر در توسعه کارستیفیکاسیون، اقلیم و در رأس آن بارش است. بارندگی با ایجاد آب زیرزمینـی بهوسیله عواملی مانند توپوگرافی، پوشش گیاهی و لیتولوژی کنترل میشود، این فاکتورها بر مقدار آبی که به درون زمین نفوذ میکند، مؤثر میباشند. با توجه به اینکه حجم بارش در توسعه کارست و پتانسیل آبهای کارستی نقش مهمی ایفا میکند، اقدام به تهیه لایه اطلاعاتی بارش گردید. طبق مبانی نظری هرچه بارندگی بیشتر باشد شرایط برای توسعه کارست مستعدتر است. به عبارتی بارندگی بالای 300 میلیمتر شرایط مساعدی را برای توسعه کارست در یک منطقه فراهم میکند. در منطقه مورد مطالعه لایه بارندگی به چهار طبقه کاملاً مناسب با بارش 300-250 میلیمتر، مناسب با بارش 250-200 میلیمتر، نسبتاً مناسب با بارندگی 200-150 میلیمتر و نامناسب با بارندگی 150-100 میلیمتر تقسسیم شده است. جهت وزن دهی بالاترین امتیاز به بارندگی بین 300-250 میلیمتر و کمترین امتیاز به مناطقی با بارش 150-100 میلیمتر اختصاص داده شد (شکل 5).
عامل دما
یکی دیگر از پارامترهای مهم و تأثیرگذار در پتانسیلیابی منابع آب کارست، عامل دما یا درجه حرارت است. طبق مبانی نظری هر چه دما پایینتر باشد امتیاز بالاتری را در توسعه منابع آب کارست میگیرد. لایه دما در محدوده مورد مطالعه به پنج طبقه تقسیمشده است که بر اساس محاسبه میانگین دمای ایستگاه سینوپتیک شهرستان بیرجند صورت گرفته است. بالاترین امتیاز را دمای بین 13-11 درجه سانتیگراد به خود اختصاص داد و کمترین امتیاز را دمای بین 23-19 درجه سانتیگراد به خود اختصاص داد. (شکل 6).
عامل ارتفاع
در مناطق مرتفعتر، آب و هوا سردتر است و در نتیجه هوازدگی فیزیکی مؤثرتر میباشد. همچنین آب دریافتی از طریق بارشها نیز بیشتر است(صفاری و همکاران،1397). یکی از عوامل فیزیکی که نقش مهمی در میزان بارش، تبخیر و تعرق، درجه حرارت و پوشش گیاهی دارد، ارتفاع میباشد. بهطورکلی عامل ارتفاع یا پستی و بلندی اهمیت زیادی در تغذیه و تخلیه و برونزد چشمههای کارستی در یک منطقه را دارا است. با افزایش ارتفاع در یک منطقه پتانسیل توسعه منابع آب کارست افزایش مییابد. زیرا با افزایش ارتفاع، میزان درجه حرارت کاهش و میزان بارندگی افزایش پیدا میکند(زروش و همکاران، 149:1393). نقشــه رســتری طبقات ارتفاعی حوضــه مورد مطالعه با اســتفاده از دســتور Reclassify در ArcGIS 10.3.1 در8 طبقه به دست آمد (شکل 7). حد ارتفاعی که دارای بیشترین مساحت است در ترازهای ارتفاعی 1700-1500 متر از سطح دریا واقع شده است. کمترین مساحت نیز مربوط به ترازهای ارتفاعی کمتر از 1352 متر است. در وزن دهی به عامل ارتفاع، بیشترین امتیاز در طبقه ارتفاعی 2695-2500 متر و کمترین امتیاز در طبقه ارتفاعی کمتر از 1352 متر را به خود اختصاص داد.
عامل شیب
شیب نقش مهمی در سرعت جریان آب ایفا میکند. این عامل نفوذ آب به درون زمین را کنترل میکند. در مناطقی که شیب ملایم است، رواناب سطحی فرصت بیشتری جهت تماس با محیط و نفوذ را دارد، در حالی که در مناطق با شیب زیاد جریان رواناب راحتتر صورت میگیرد و این امر باعث کاهش نفوذ آب باران میشود (پور اکبری و همکاران،1399). از فاکتورهای توپوگرافی مؤثر که گرادیان هیدرولیکی و جهت حرکت آب زیرزمینی را تعیین میکند، شیب توپوگرافی است. میزان شیب هم در میزان رواناب حاصل از بارش و هم در میزان نفوذ آب به داخل زمین و فرایند انحلال توسط بارندگی نقش مؤثری ایفا میکند(زروش و همکاران، 1393 :150). نقشه نهایی شیب منطقه مورد مطالعه در 8 طبقه از کلاس 3-0 درصد تا بیشتر از 34 درصد به دست آمد (شکل 8). با توجه به مبانی نظری، شیبهای کم و مناطق مسطح بیشترین پتانسیل در کارستی شدن را دارا میباشند. در حوضه مورد مطالعه شیب 15-0 درصد در قلمرو شیب کم میباشند. گسترش پدیده کارست درشیبهای کم به دلیل تماس بیشتر آب با آهک بیشتر بوده و پدیده انحلال در این شیبها تشدید میشود. در حوضه مورد مطالعه بیشترین امتیاز در وزن دهی را شیب بین 0-3 درصد و کمترین امتیاز را شیبهای بالاتر از 34 درصد به خود اختصاص داد.
عامل فاصله از آبراهه
نوع شبکه زهکشی هر منطقه توسط لیتولوژی واحدهای زمینشناسی، توپوگرافی و ساختارهای تکتونیکی و زمینشناسی منطقه کنترل میشود. از عوامل مهم در پتانسیل کارست زایی، هیدرولوژی، تراکم آبراههها و فاصله از آنهاست که هرچه فاصله از آبراهه کمتر باشد میزان توسعه منابع آب کارست بیشتر خواهد بود(رضایی و همکاران،1399). وزن دهی به این صورت است که در مکانهایی که فاصله از آبراهه کم است وزن بیشتری به خود اختصاص میدهد چون فرایند کارستی شدن بیشتر خواهد بود و هرقدر فاصله از آبراهه بیشتر شود، میزان توسعه کارست کمتر خواهد بود. شکل 9 فاصله از آبراهه در منطقه مورد مطالعه را نشان میدهد. فاصله از آبراهه برای منطقه مورد مطالعه به پنج دسته طبقهبندی شده که هرچه فاصله از آبراهه کمتر، وزن بیشتری به خود اختصاص داده است. فاصله از آبراهه تا 500 متری بیشترین وزن را به خود اختصاص داده است یعنی بیشتر مستعد توسعه کارست است زیرا آبراههها دارای نفوذپذیری خوبی برای تغذیه منابع آب کارست میباشند و فاصله بیشتر از 200000 متر کمترین وزن را به خود اختصاص داده اســت یعنی کمترین پتانســیل را در منطقه برای توســعه منابع آب کارســت دارد (شکل 9).
عامل جهت شیب
تأثیر جهت شـیب بر روی تنوع پوشـش گیاهی و بعضـی فرآیندهای هیدرولوژیک مانند ذوب برف اهمیت زیادی دارد)علیزاده، 513:1389). با توجه به اهمیت کیفیت تابش نور خورشـــید در تأمین انرژی مورد نیاز مناطق، جهتگیری دامنهها نقش مهمی در این ارتباط دارد. بهطور کلی، مناطق هموار بیشترین اهمیت را در پتانسیلیابی منابع آب کارست دارا میباشند. مناطق با جهتهای شمالی، شرقی و شمال شرقی به دلیل تبخیر کمتر و ماندگاری بیشتر برف بیشترین امتیاز را به خود اختصاص دادند و بهترین جهتها بهمنظور توسعهیافتگی منابع آب کارست هستند. کمترین امتیاز به جهتهای جنوبی، جنوب شرقی و جنوب غربی داده شد زیرا در این جهتها، تبخیر بیشتر ، نور آفتاب بیشتر و ماندگاری برف هم کمتر است(شکل10).
عامل کاربری اراضی
برای بررسی میزان قابلیت پوشش زمین و کاربری اراضی جهت پتانسیل کارست زایی در محدوده موردنظر، از نقشه کاربری اراضی استفاده شده است. بیشترین وزن به جنگلها، درختزارها و بعد از آن به مراتع و کمترین امتیاز به مناطق شهری و بعد از آن به مناطق کشاورزی داده شد. در مناطق جنگلی به علت وجود برگ گیاهان، میزان هوموس خاک افزایش یافته و نفوذپذیری را افزایش میدهد. همچنین فشار ریشه گیاهان در این مناطق سبب تخریب و شکستگی سنگها شده که عامل مؤثری در افزایش میزان نفوذپذیری است. هر چه تراکم پوشش گیاهی بیشتر باشد، از میزان رواناب ناشی از بارندگیهای شدید جلوگیری شده و زمان لازم جهت نفوذ ریزشهای جوی فراهم میگردد(پور اکبری و همکاران،1399). دلیل اینکه بیشترین وزن به جنگل اختصاص یافته این است که ترشح اسید فولیک از طریق ریشه در جنگل خیلی بیشتر از بقیه پارامترها است (شکل 11).
شکل 3: نقشه زمینشناسی شهرستان بیرجند شکل 4: نقشه فاصله از گسل شهرستان بیرجند
(نقشه زمینشناسی 1:100000 بیرجند)
شکل 5: نقشه هم باران شهرستان بیرجند شکل 6 : نقشه همدمای شهرستان بیرجند
شکل7 : نقشه طبقات ارتفاعی شهرستان بیرجند شکل 8: نقشه شیب شهرستان بیرجند
شکل 9 : نقشه فاصله از آبراهه شهرستان بیرجند شکل 10 : نقشه جهت شیب شهرستان بیرجند
شکل 11: نقشه کاربری اراضی منطقه مورد مطالعه
نقشه نهایی پتانسیلیابی منابع آب کارست با مدل تحلیل سلسله مراتبی
نقشه نهایی پتانسیلیابی منابع آب کارست با استفاده از لایههای رستری لیتولوژی، فاصله از گسل، دما، بارندگی، کاربری اراضی، شیب، جهت شیب، ارتفاع و فاصله از آبراهه با ترکیب وزنهای به دست آمده از روش تحلیل سلسله مراتبی با هر لایه و همپوشانی آنها با استفاده از دستور Raster Calculator دست آمد. بر اساس اصول مبانی کارست و مشاهدات میدانی، نقشه نهایی حاصل از پتانسیلیابی منابع آب کارست در محدوده مورد مطالعه به پنج طبقه: پتانسیل خیلی زیاد، زیاد ، متوسط، کم و خیلی کم طبقهبندی شد. شکل 12 نقشه نهایی پتانسیلیابی منابع آب کارست و جدول 2 مساحت و درصد طبقات توسعه کارست به روش تحلیل سلسله مراتبی را نشان میدهد که از کل مساحت شهرستان بیرجند، 40/26 درصد در طبقه خیلی کم توسعهیافته، 58/27 درصد در طبقه کم توسعهیافته، 86/23 درصد در طبقه متوسط، 75/15 درصد در طبقه توسعه زیاد و 40/6 درصد در طبقه توسعه خیلی زیاد گرفته است که در مناطق جنوبی و شمال غربی محدوده مورد مطالعه قرار گرفته است. بنابراین از کل مساحت شهرستان بیرجند 15/22درصد در طبقه زیاد و خیلی زیاد از نظر توسعه کارست را دارا میباشند(جدول5). از اینرو عامل لیتولوژی منطقه با ارزش 41/0 بیشترین وزن و مهمترین عامل کنترل کننده پتانسیل توسعه منابع آب کارست در منطقه مورد مطالعه بوده است و عامل کاربری اراضی با ارزش 017/0 کمترین وزن را به خود اختصاص داده است که کمترین تأثیر را در کارست زایی کنونی حوضه به خود اختصاص داده است. نتایج نشانگر آن است که به ترتیب عوامل زمینشناسی، فاصله از گسل و فاصله از آبراهه مهمترین نقش را در توسعه کنونی کارست در محدوده مورد مطالعه داشتهاند.
شکل 12 : نقشه پتانسیل منابع آب کارستی در منطقه مورد مطالعه
جدول 5 : مساحت و درصد طبقات توسعه کارست به روش AHP
طبقه | مساحت به کیلومترمربع | درصد مساحت طبقات | |
1 | خیلی کم | 1043 | 40/26 |
2 | کم | 1089 | 58/27 |
3 | متوسط | 942 | 86/23 |
4 | زیاد | 622 | 75/15 |
5 | خیلی زیاد | 253 | 40/6 |
صحت سنجی و اعتبار سنجی نقشه پتانسیل منابع آب کارست
جهت صحت سنجی نقشه نهایی، لایه آبخوانهای شهرستان بیرجند بر روی نقشه نهایی فراخوانی گردید و دقیقاً بر روی مناطق توسعه کارست خیلی زیاد و زیاد قرار گرفت. موقعیت آبخوانها در جنوب شهرستان بیرجند قرار گرفته که بر روی نقشه پتانسیلیابی منابع آب کارست کاملاً نمایان است. (شکل13)
شکل13: تلفیق و اعتبار سنجی لایه آبخوانهای شهرستان بیرجند و لایه پتانسیل کارستی شدن در محدوده مورد مطالعه
نتیجهگیری
هدف از این پژوهش پتانسیل سنجی منابع آب کارست در شهرستان بیرجند در استان خراسان جنوبی با استفاده از مدل AHP است. بر اساس اصول مبانی کارست و مشاهدات میدانی و یافتههای دیگر، هشت پارامتر بهعنوان عوامل مؤثر بر توسعه کارست در منطقه در نظر گرفته شـد. از کل مساحت شهرستان بیرجند، 40/26 درصد در طبقه خیلی کم توسعهیافته، 58/27 درصد در طبقه کم توسعهیافته، 86/23 درصد در طبقه متوسط، 75/15 درصد در طبقه توسعه زیاد و 40/6 درصد در طبقه توسعه خیلی زیاد گرفته است که مناطق با توسعه زیاد و زیاد در مناطق جنوبی و شمال غربی محدوده مورد مطالعه قرار گرفته است. بنابراین از کل مساحت شهرستان بیرجند 15/22درصد در طبقه زیاد و خیلی زیاد از نظر توسعه کارست قرار گرفتند. طبق قضاوت کارشناسان و متخصصان کارست ، عامل لیتولوژی منطقه با ارزش 41/0 بیشترین وزن و مهمترین عامل کنترل کننده پتانسیل توسعه منابع آب کارست در منطقه مورد مطالعه بوده است و عامل کاربری اراضی با ارزش 017/0 کمترین وزن را به خود اختصاص داده است که کمترین تأثیر را در کارست زایی کنونی حوضه به خود اختصاص داده است. نتایج نشانگر آن است که به ترتیب عوامل لیتولوژی، فاصله از گسل و فاصله از آبراهه مهمترین نقش را در توسعه کنونی منابع آب کارست در این شهرستان داشتهاند. جهت صحت سنجی نقشه نهایی، لایه آبخوانهای شهرستان بیرجند بر روی نقشه نهایی فراخوانی گردید و دقیقاً بر روی مناطق توسعه کارست خیلی زیاد و زیاد قرار گرفت. موقعیت آبخوانها در جنوب شهرستان بیرجند قرار گرفته که بر روی نقشه پتانسیلیابی کاملاً نمایان است. با توجه به تهیه نقشه نهایی و اعتبار سنجی نقشه مذکور، میبایست مدیران سازمان آب منطقهای خراسان جنوبی اقدام به اکتشاف منابع آب در مناطق مذکور نمایند. با توجه به اینکه منابع آب شرب شهرستان از شهرستان سربیشه منتقل میگردد، پیشنهاد میگردد جهت جلوگیری از بحران منابع آب در آینده و افزایش کیفیت منابع آب، اکتشاف منابع آب کارست در مناطق جنوبی در ارتفاعات باقران و مناطق شمال غربی شهرستان صورت پذیرد.
منابع
1- اشتوکلین(1368): زمینشناسی و زمینساخت فلات ایران، مترجمان: صادق حداد کاوه و حسن حسنعلی زاده، انتشارات وزارت فرهنگ و آموزش عالی، علمی و فرهنگی، 108.
2- بهنیافر، ابوالفضل و قنبر زاده، هادی (1395): ژئومورفولوژیکارست، مشهد، انتشارات نگاران سبز.
3- پور اکبری، سجاد، کلانتری، نصرا...، مصلح، آرش و عقدکی، یاسر(1399): پتانسیلیابی منابع آب کارستی با استفاده از RS ، GIS و AHP (مطالعه موردی: تاقدیسهای لیلی و کینو در شمال شرق خوزستان). مجله مهندسی منابع آب، سال 13، شماره 2، صص112-99.
4- ثقفی، مهدی، امیر احمدی، ابوالقاسم و ربیعی، مریم(1395): تشخیص واکنشهای ژئومورفیک مخروطه افکنه های جنوب دشت بیرجند به فعالیتهای تکتونیکی و تغییرات آب و هوایی با استفاده از دادههای میدانی، نشریه جغرافیا و برنامهریزی، شماره 55، صص 114-97.
5- دشتی برمکی، مجید، رضایی، محسن و اشجاری، جواد(1394): پتانسیلیابی منابع آب کارست کوههای دوان و شاپور بر اساس تصمیمگیری چندمعیاره، مجله پژوهش آب ایران، شماره1، صص100-89.
6- رضایی عارفی، محسن، زنگنه اسدی، محمدعلی، بهنیافر، ابوالفضل و جوانبخت، محمد(1399): پهنهبندی تحول کارست با استفاده از مدل تحلیل سلسله مراتبی (مطالعه موردی:حوضه کوهستانی کلات، خراسان رضوی)، مجله جغرافیا، سال 18، شماره 64، صص94-79
7- زروش، ناهید، واعظی، عبدالرضا و کریمی، حاجی(1393): ارزیابی پتانسیل توسعه کارست در تاقدیس کبیر کوه استان ایلام با استفاده از تلفیق فازی و روش تحلیل سلسله مراتبی AHPو سنجش از دور و GIS، مجله پژوهشهای ژئومورفولوژی کمی، سال 3، شماره 3، صص 144-157.
8- سپند، ساحل، چیتسازان، منوچهر، رنگزن، کاظم و میرزایی، یحیی(1386): تلفیق سنجش از دور و GIS در پتانسیلیابی منابع آب زیرزمینی محدوده لالی، همایش ژئوماتیک تهران.
9- صفاری، امیر، کیانی، طیبه و زنگنه تبار، ساسان(1398): بررسی عوامل مؤثر در توسعهیافتگی و پهنهبندی کارست کوهستان خورین با استفاده از منطق فازی، نشریه تحقیقات کاربردی علوم جغرافیایی، سال 19، شماره 55، صص36-23.
10- علیزاده، امین(1389): اصول هیدرولوژی کاربردی، انتشارات دانشگاه امام رضا، چاپ سیام، مشهد، ایران.
11- قبادی، محمدحسین(1388): زمینشناسی مهندسی کارست، چاپ دوم، انتشارات دانشگاه بوعلی سینا، همدان، ایران.
12- مزیدی، احمد، کرم، امیر و کوراوند، مژگان(1395): پتانسیلیابی توسعه کارست با استفاده از منطق فازی (مطالعه موردی: حوضه دشت سوسن و دشت ایذه)، مجله پژوهشهای ژئومورفولوژی کمی، سال پنجم، شماره 2، صص 141-130.
13 -مکرم، مرضیه و نگهبان، سعید(1398): بررسی و شناسایی مناطق دارای پتانسیل کارستی شدن با استفاده از روش فازی و مدل تحلیل سلسلهمراتبی، مجله جغرافیا و برنامهریزی محیطی، دوره 30، شماره2، 134-121
14- محمودی، فرجا...، (1385): ژئومورفولوژی ساختمانی، چاپ هشتم، انتشارات پیام نور.
15- نخعی، اعظم، زنگنه اسدی، محمدعلی، بهنیافر، ابوالفضل و گلی مختاری، لیلا(1402): پهنهبندی توسعه کارست در حوضه آبریز بقمچ کشف رود با مدل منطق فازی ANP. مجله پژوهشهای ژئومورفولوژی کمی، سال 11، شماره 4، صص255-230
16- واعظی هیر، عبدالرضا، وفادار، مرتضی و آقایی، واحده(1397): پتانسیلیابی منابع آب زیرزمینی موجود در واحدهای کارستی و سازند سخت کوه مورو-صوفیان با بهکارگیری روشهای تصمیمگیری چند معیاره AHP، SAW و F-AHP، مجله فضای جغرافیایی، سال 18، شماره 64، 234-215
17- یمانی، مجتبی، شمسیپور، علیاکبر، جعفری اقدم، مریم و باقری سید شکری، سجاد (1392): بررسی عوامل مؤثر در توسعهیافتگی و پهنهبندی کارست حوضه چله با استفاده از منطق فازی و تحلیل سلسله مراتبی، استان کرمانشاه، مجله علوم زمین سال، شماره 22، صص 66-57.
18- Abdelkareem, M., & Abdalla, F. (2022): Revealing Potential Areas Of Water Resources Using Integrated Remote-Sensing Data And GIS-Based Analytical Hierarchy Process. Geocarto International, 37(25), 8672-8696.
19- Calic, J., (2011): Karstic Uvula Revisited: Toward A Redefinition Of The Term, Geomorphology, Vol.134, No.1, Pp.32-42
20- Ford, D., & Williams, P. D. (2007): Karst Hydrogeology And Geomorphology. West Sussex, England: John Wiley & Sons
21- Hamdani, N., & Baali, A. (2020). Characterization Of Groundwater Potential Zones Using Analytic Hierarchy Process And Integrated Geomatic Techniques In Central Middle Atlas (Morocco). Applied Geomatics, 12(3), 323-335.
22- Jaiswal, R, K. Ghosh, N, C. Galkate, R, V. Thomas, T. 2015. Multi Criteria Decision Analysis (MCDA) For Watershed Prioritization. International Conference In Water Resources, Coastal And Ocean Engineering (ICWECOE 2015). 1553-1560.
23- Milanovic, P. J. (1981): Karst Hydrogeology. Colorado, CO: Water Resources Ublications
24- Lamoreaux, P.E., 2007, Karst: The Foundation For Concepts In Hydrogeology. Bull Eng Geol Environ, 51. Studies, Vol. 64, No.1, Pp. 23-33.
25- Ozcelik, M., & Sarp, G. (2018): Evaluation Of Sustainable Water Supply Alternatives In Karstified Rock Masses Using GIS And AHP Methodology For Antalya (Turkey) Urban Area. Environmental Earth Sciences, 77, 1-14
26- Palmer, A. N., (2007): Cave Geology: Dayton, Ohio, Cave Books, P. 454.22
27- Saaty, T.L., (1980): The Analytic Hierarchy Process, Mcgraw _ Hill, Inc., Reprinted By RWS Publications, Pittsburgh
28- Waele, J., Gutierrez, F., Audra, P., (2015): Karst Geomorphology: From Hydrological Functioning To Palaeoenvironmental Reconstructions Geomorphology, Vol. 229, Pp. 1-2.
29- White, W. B., (1988): Geomorphology And Hydrology Of Karst ،Oxford University Press. Quinlan, J, Groundwater Monitoring In Karst Terrains, EPA. 600/ X.
[1] * نیسنده مسئول: 09151337352 Email: rezaei.arefi@cfu.ac.ir
[2] 1- Milanovic
[3] - Ford and williams
[5] - Palmer
[6] - Lamoreaux
[7] -Calic
[8] 7- Ozcelik
[9] 8- Hamdani and Baali
[10] 9- saaty