Document Type : Original Article
Authors
Abstract
Keywords
1- مقدمه
خشکسالی به عنوان یکی از مهمترین بلایای طبیعی بوده که با وقوع تدریجی آرام و خزنده اش بر ابعاد مختلف زندگی بشر تاثیر می گذارد. این بلا به عنوان یک پدیده ناگوار اقلیمی که بطور مستقیم جوامع را از طریق محدودیت دردسترسی به منابع آب تحت تاثیر قرار می دهد، هزینه های اقتصادی، اجتماعی و محیطی زیادی را به همراه دارد (گودارد[1] و همکاران، 2003)
در گذشته، رویکرد مرسوم پایش خشکسالی، بر پایه مشاهدات اقلیمی هواشناسی و استفاده از شاخص های اقلیمی بوده است که این روش ها به علت استفاده از آمار ایستگاههای هواشناسی به صورت نقطه ای می باشند و از طرفی پراکندگی ایستگاهها بخصوص در مناطق خشک محدودیت هایی را در بررسی خشکسالی بوجود آورده است. در چهار دهه اخیر، سنجش از دور ابزار پایش خشکسالی را به صورت وسیع و گسترده فراهم آورده و مدل های پایش خشکسالی زیادی ارائه شده است که عموما برپایه شاخص های گیاهی، دمای سطحی زمین (LST)[2] ، رطوبت و بازتاب در ناحیه مرئی و مادون قرمز است که از آن جمله می توان به شاخص محتوای آب برگ(LWCI)[3] (هانت[4] و همکاران، 1987)، شاخص آنومالی پوشش گیاهی (AVI)[5](چن[6] و همکاران، 1994)، شاخص فراهم بودن آب (VSWI)[7] (کلارسون[8] و همکاران، 1994) و شاخص شرایط پوشش گیاهی (VCI)[9] (کوگان[10]، 1995) اشاره نمود.
به منظور محاسبه شاخصهای فوق الذکر نیاز به بازیابی LST و آلبیدو از تصاویر ماهواره ای، تصحیحات آتمسفری داده های ماهواره ای،تجزیه پیکسل های آمیخته و تبدیل باندهای طیفی باریک به باندهای طیفی پهن می باشد (لیانگ[11]، 2003؛ ژائو[12] و همکاران، 2000). بنابراین خطای نهایی استخراج و کمی کردن اطلاعات خشکسالی می تواند زیاد باشد. از طرفی دمای سطحی زمین و شاخص اختلاف پوشش گیاهی نرمال شده (NDVI)[13] پارامترهایی هستند که خشکسالی پس از وقوع با تاخیر زمانی روی آنها تاثیر می گذارد و این مانع پایش لحظه ای[14] خشکسالی می شود. از سوی دیگر شاخص های ذکر شده به منظور پایش خشکسالی نیاز به سری تصاویر چند زمانه ماهواره ای دارند که این خود باعث افزایش حجم پردازش داده ها، هزینه و زمان است.
به منظور رفع این محدودیت ها و مشکلات، غلام[15] در سال 2006 شاخص خشکسالی عمودی (PDI)[16] را بر پایه بازتاب در ناحیه قرمز و مادون قرمز نزدیک ارائه نمود. نتایج حاصل از این شاخص همبستگی بالای 75 درصد با شاخص های خشکسالی مبتنی بر داده های اقلیمی و داده های رطوبتی و همچنین سایر شاخص های مبتنی بر تصاویر ماهواره ای نشان داد.
شاخص PDI نسبت به سایر شاخص های خشکسالی مبتنی بر سنجش از دور، ساده تر و کارآتر است و تنها براساس رفتار رطوبتی پدیده ها در باند قرمز و مادون قرمز محاسبه می شود.
از آنجایی که رطوبت و پوشش گیاهی از جمله عوامل مهم در بررسی خشکسالی می باشند، غلام و همکاران در سال 2007 شاخص PDI با عنوان شاخص خشکسالی عمودی اصلاح شده (MPDI)[17] بهبود بخشیدند. این شاخص علاوه بر رفتار طیفی رطوبت در باند قرمز و مادون قرمز ،درصد پوشش گیاهی(FV)[18] را نیز در نظر می گیرد. نتایج نشان داد که این شاخص نسبت به شاخص PDI، از درصد همبستگی بیشتری با شاخص خشکسالی محاسبه شده از طریق داده های زمینی برخوردار است.
در ایران مطالعات متعددی به منظور بررسی و ضعیت خشکسالی با استفاده از داده های ماهواره ای صورت گرفته است که از آن جمله می توان به مطالعه چنار در سال 1380 در استان های شمال غرب کشور، صفری در سال 1382 در منطقه زابل و باعقیده در سال 1386 اشاره نمود.
در مطالعات انجام شده مذکور عمدتاً روند تغییرات پوشش گیاهی را با استفاده از شاخص NDVI و بوسیله داده های ماهواره ای AVHRR [19] که دارای دقت مکانی پایین (یک کیلومتر) هستند، بررسی نموده اند. این داده ها علی رغم اینکه قدرت تفکیک زمانی بالایی دارند ولی می توانند به علت تفکیک مکانی پایین محدودیت هایی را ایجاد نمایند.
نکته قابل توجه دیگر این است که در تحقیقات انجام گرفته در زمینه خشکسالی با استفاده از سنجش از دور، عمدتا به بررسی خشکسالی با استفاده از شاخص های خشکسالی مبتنی بر تصویر پرداخته شده است که این شاخص ها تنها درصد پوشش گیاهی و دمای سطحی زمین را در نظر می گیرند در حالیکه علاوه بر موارد ذکر شده عوامل متعدد دیگری از جمله عوامل محیطی (شیب، ارتفاع، جهت، درصد پوشش گیاهی و ...) و عوامل اقلیمی (بارندگی سالیانه، دما و تبخیر) نیز بر خشکسالی تاثیرگذار است.
بر این اساس هدف از تحقیق حاضر پهنه بندی شدت خشکسالی حوضه شیطور در استان یزد با تلفیق داده های محیطی، ماهواره ای و هواشناسی با استفاده از روش های دانش مبنا در محیط GIS و انتخاب روش مناسب تلفیق می باشد.
2- محدوده مورد مطالعه
ناحیه مطالعاتی از نظر جغرافیایی بین عرضهای ً25 َ31 ْ31 و ً15 َ45 ْ31 شمالی و طولهای ً46 َ45 ْ55 و ً22 َ03 ْ56 شرقی واقع شده است. این منطقه در استان یزد، در بخش جنوبی شهرستان بافق در فاصله حدود 25 کیلومتری شهر بافق از سمت شرق، در حد فاصل مرز استان یزد و کرمان قرار گرفته است (شکل شماره 1).
شکل (1)- موقعیت محدوده ی مورد مطالعه
مساحت کل محدوده شیطور برابر 416 کیلومترمربع می باشد. بلندترین نقطه ارتفاعی با 2989 متر و پست ترین آن 1844 متر می باشد. متوسط بارندگی محدوده ی مورد مطالعه 171 میلیمتر، دمای متوسط سالانه 14.1 درجه سانتی گراد و اقلیم منطقه به روش دومارتن اصلاح شده از نوع خشک سرد می باشد و جزو مناطق خشک مرکزی ایران محسوب می شود.
پوشش گیاهی حوزه از 4 تیپ گیاهی اصلی شامل درمنه زار دشتی ،درمنه زارکوهی ، ،هزار خار دشتی و گون تشکیل شده است.
3- داده های مورد استفاده
به منظور انجام تحقیق حاضر، از تصاویر سنجنده AVNIR2[20] ماهواره ALOS [21] مربوط به تاریخ 27 تیر 1388 (با قدرت تفکیک مکانی 10 متر و 4 باند طیفی)، نقشه های توپوگرافی با مقیاس 1:25000 جهت زمین مرجع نمودن تصاویر و همچنین داده های اقلیمی و هواشناسی ایستگاههای موجود در محدوده منطقه مورد مطالعه استفاده شده است.
4- عوامل موثر بر خشکسالی
از عوامل موثر بر خشکسالی می توان به موارد زیر اشاره نمود (شایگان، 1381)؛
ارتفاع: تغییرات ارتفاعی بر درجه حرارت و مقدار رطوبت تاثیر مستقیم دارد. لذا هر چه ارتفاع یک منطقه بیشتر باشد شدت خشکسالی کمتر خواهد بود. لازم به ذکر است این حالت زمانی درست است که سایر پارامترهای موثر که در ادامه به آن اشاره خواهد شد ثابت باقی بماند.
شیب: شیب زمین بر مقدار آب نفوذی و پایداری سطح خاک تاثیر گذار است. هرچه شیب بیشتر باشد، آب با سرعت بیشتری از منطقه خارج شده و کمتر نفوذ می کند. لذا با توجه به اینکه با افزایش شیب، آب در دسترس گیاهان کاهش پیدا می کند و خاک رطوبت کمتری دارد، نواحی با شیب بیشتر (با ثابت در نظر گرفتن سایر عوامل) نسبت به خشکسالی حساسیت بیشتری دارند.
جهت: جهت قرار گرفتن عوارض سطح زمین بر میزان دریافت نور تاثیر گذار است. در نیمکره شمالی شیبهای رو به جنوب گرمتر و خشکتر از شیبهای رو به شمال است. جهت و زاویه ی شیب تغییر محسوسی در شدت نور و مدت آفتاب روزانه ایجاد می کند که این امر به اختلافات فاحش دمای موثر بر پوشش گیاهی می انجامد. بنابراین به ترتیب شیب های جنوبی، غربی، شرقی و شمالی به خشکسالی حساس می باشند.
تراکم شبکه آّبراهه ها: تراکم شبکه ی آبراهه ها بر میزان در دسترس بودن آب در یک منطقه تاثیر گذار است، هر چه تراکم آبراهه ها بیشتر باشد، سرعت تخلیه ی آب بیشتر بوده و آب زمان کمتری برای نفوذ دارد و در نتیجه میزان رطوبت خاک کمتر خواهد بود. بنابراین با افزایش تراکم زهکشی میزان حساسیت به خشکسالی بیشتر خواهد بود.
بارندگی: با توجه به شرایط وقوع خشکسالی، باران مهمترین پارامتر موثر در ایجاد خشکسالی است. مناطق با بارندگی بیشتر حساسیت کمتری به خشکسالی دارند.
دما: دما به صورت غیر مستقیم بر مقدار آب قابل دسترس تاثیر گذار است. لذا هر چه دما بیشتر باشد حساسیت به خشکسالی بیشتر خواهد بود.
تبخیر: تبخیر باعث از دسترس خارج شدن آب قابل دسترس در یک منطقه می شود. این پدیده باعث کاهش آب سطحی، رطوبت خاک و حتی آب زیرزمینی می شود. بنابراین مناطق با تبخیر بیشتر حساسیت بیشتری نسبت به خشکسالی دارند.
پوشش گیاهی: خشکسالی اثر مستقیم بر پوشش گیاهی دارد یعنی هرچه شدت خشکسالی بیشتر باشد، کاهش پوشش گیاهی بیشتر است. همچنین می توان گفت هرچه تراکم پوشش گیاهی بیشتر باشد نیاز آبی منطقه بیشتر و در نتیجه حساسیت به خشکسالی بیشتر است. بنابراین پوشش گیاهی متراکم تر آسیب پذیری بیشتری به شرایط کم آبی خواهد داشت.
کاربری اراضی: بطور کلی اکوسیستم های مصنوعی نسبت به اکوسیستم های طبیعی از حساسیت بیشتری نسبت به خشکسالی برخوردارند. بنابراین اراضی باغی- زراعی حساسیت بیشتری به خشکسالی نسبت به اراضی مرتعی دارند.
5- روش تحقیق
5-1- پیش پردازش تصاویر
با توجه به اینکه هدف از این تحقیق محاسبه درصد پوشش گیاهی با استفاده از میزان بازتاب[22] می باشد، بنابراین عدد رقومی ثبت شده توسط سنجنده، ابتدا به رادیانس[23] تبدیل شد. به منظور محاسبه رادیانس برای تصاویر سنجنده AVNIR از رابطه زیر استفاده می شود:
که Li رادیانس در باند i، Gain عدد ثابتی است که مقدار آن برای هر باند تعریف شده است و DNi عدد رقومی ثبت شده توسط سنجنده می باشد. جهت دستیابی به بازتاب واقعی زمینی و همچنین تصحیح اتمسفری تصاویر، از الگوریتم فلش[24] در نرم افزار ENVI4.7 استفاده شد. پارامترهای مورد نیاز تصحیح تصاویر AVNIR2 و مقادیر تنظیم شده در جدول شماره (1) آمده است.
جدول 1) پارامترهای تنظیم شده به منظور تصحیح اتمسفری تصویر مورد استفاده
مقدار تنظیم شده |
پارامتر مورد نیاز |
مقدار تنظیم شده |
پارامتر مورد نیاز |
18/7/2009 |
تاریخ تهیه ی تصویر |
81/31 |
عرض جغرافیایی مرکز تصویر |
2 |
متوسط ارتفاع زمینی (کیلومتر) |
10/56 |
طول جغرافیایی مرکز تصویر |
روستایی |
مدل آئروسل |
692 |
ارتفاع ماهواره (کیلومتر) |
خیر |
بازیابی آئروسل و آب |
تابستان عرضهای میانی |
مدل اتمسفری |
به منظور تصحیح هندسی تصاویر، تعداد 17 نقطه به صورت پراکنده در تمام تصویر (942 کیلومتر مربع) و از روی نقشه های 25000/1 انتخاب شد، سپس با استفاده از معادلات پلی نومیال درجه دو با دقت 0.65 پیکسل و روش نمونه برداری نزدیکترین همسایه، بازنویسی [25] شد.
5-2- تهیه لایه های اطلاعاتی پارامترهای موثر بر خشکسالی
با توجه به اینکه هدف این تحقیق، بررسی میزان تاثیر عوامل موثر بر خشکسالی و در نهایت رسیدن به یک نقشه ی حساسیت خشکسالی است، بنابراین لازم است تا هر یک از پارامترهای موثر برخشکسالی وارد تحلیل گردند و سپس با استفاده از این تحلیل ها به نقشه ی نهایی دست پیدا کرد.
بدین منظور مدل رقومی ارتفاع با استفاده از لایه های اطلاعاتی منحنی میزان، نقاط ارتفاعی و آبراهه نقشه های توپوگرافی با مقیاس 25000/1 با دقت 10 متر تهیه گردید. لازم به ذکر است که دقت اسمی مسطحاتی نقشه های 1:25000، 5 متر می باشد ولی از آنجایی که قدرت تفکیک مکانی تصاویر استفاده شده 10 متر می باشد، بنابراین ابعاد 10 متر به عنوان اندازه ی پایه برای تمامی لایه های اطلاعاتی در نظر گرفته شد.
در مرحله بعد، لایه های اطلاعاتی شیب و جهت از مدل رقومی ارتفاع تهیه گردید. پس از محاسبه ی شیب برحسب درصد، نقشه ی حاصله با توجه به طبقه بندی کارپنتر به 6 کلاس و نقشه ی جهت حاصله نیز، به نقشه ی با چهار جهت اصلی طبقه بندی گردید. در نهایت به منظور حذف مجموعه پیکسل های آمیخته، نواحی با مساحت کمتر از چهار هکتار با طبقات بزرگتر ترکیب گردید.
به منظور استخراج اتوماتیک شبکه ی آبراهه ها و زیرحوضه ها، ابتدا چاله های هیدرولوژیکی[26]، با استفاده از اعمال فیلتر روی مدل رقومی ارتفاع حذف گردید. سپس مدل رستری جهت جریان[27] و مدل جریان تجمعی[28] تهیه گردید. در نهایت به منظور استخراج اتوماتیک آبراهه ها، پیکسل های با جریان تجمعی بیشتر از 1000 به عنوان آبراهه درنظر گرفته شد و بر اساس آن زیرحوضه های موجود در حوضه استخراج گردید.
به منظور محاسبه ی تراکم شبکه ی آبراهه ها، مجموع طول آبراهه های هر زیرحوضه بر مساحت آن زیرحوضه تقسیم و برای هر زیر حوضه تراکم شبکه ی آبراهه ها محاسبه شد.
به منظور تهیه لایه بارندگی از رابطه رگرسیون بدست آمده بین ارتفاع ایستگاههای هواشناسی و میزان بارندگی سالیانه بدست آمد (شکل 2).
شکل (2). رابطه ارتفاع – بارندگی ایستگاههای موجود در محدوده مورد مطالعه
به منظور تهیه لایه اطلاعاتی متوسط دمای سالانه، پس از بازسازی آمار ناقص، جهت تبدیل داده های نقطه ای به سطح پیوسته از روش درونیابی عکس فاصله با توان 2 استفاده شد به منظور تهیه ی لایه ی تبخیر نیز مانند لایه ی دما عمل شد و از روش IDW استفاده شد.
ابراهیمی خوسفی در سال 1389 در تحقیقی که با استفاده از شاخص های گیاهی انجام داد به این نتیجه رسید که شاخص MSAVI1 در مقایسه با سایر شاخص ها با دقت بالاتری بیانگر درصد پوشش گیاهی منطقه می باشد. لذا به منظور محاسبه درصد پوشش گیاهی در این تحقیق از شاخص MSAVI1 استفاده شد. این شاخص با استفاده از رابطه زیر محاسبه می شود:
RRED بازتاب در ناحیه قرمز، RNIR بازتاب در ناحیه مادون قرمز نزدیک، NDVI و WDVI دو شاخص رایج پوشش گیاهی و L: ضریب کاهش دهنده اثر خاک می باشد که بسته به نوع پوشش بین 0 تا 1 تغییر می کند.
پس از محاسبه ی شاخص فوق، با استفاده از رابطه ی زیر، درصد پوشش گیاهی در محدوده ی مورد مطالعه محاسبه شد.
تهیه نقشه کاربری اراضی با استفاده از تصاویر ماهواره ای تهیه و محدوده مورد مطالعه به دو بخش مرتع و اراضی زراعی- باغی تقسیم گردید.
5-3- استاندارد سازی لایه های اطلاعاتی
پس از تعیین نوع تاثیر هر یک از پارامترهای موثر بر خشکسالی، به منظور قابل مقایسه کردن لایه های مختلف اطلاعاتی، بایستی لایه ها استاندارد سازی شوند. در تحقیق حاضر به منظور استاندارد سازی لایه ها از روش منطق فازی خطی استفاده شد. نحوه ی استاندارد سازی برای هر لایه در ادامه آمده است:
5-4- تعیین وزن پارامترهای موثر برخشکسالی
در تحقیق حاضر، به منظور تعیین وزن عوامل موثر برخشکسالی، از روش سلسله مراتبی (AHP)[29] استفاده شد. در این روش با تشکیل ماتریس مقایسه ی زوجی، وزن عوامل مختلف نسبت به یکدیگر با نظر کارشناسان تعیین می گردد. با تشکیل ماتریس مقایسه ی زوجی، ضرایب برداری هر یک از عوامل محاسبه و اوزان مربوط به عوامل موثر تعیین گردیدند. در نهایت شاخص سازگاری به منظور تایید درستی اوزان محاسبه شده، محاسبه گردید.
5-5- ترکیب لایه ها، تهیه نقشه نهایی و ارزیابی نتایج
به منظور ترکیب لایه های مختلف، از روش شاخص وزنی[30] و اپراتورهای مختلف منطق فازی (AND، OR، ضرب و جمع) استفاده گردید و نقشه ی خطر پذیری به ازای هر روش تهیه گردید. پس از تهیه ی نقشه ی نهایی خطر پذیر خشکسالی، به منظور ارزیابی نتایج حاصله از شاخص عمودی خشکسالی اصلاح شده (MPDI) استفاده گردید. این شاخص از رابطه زیر محاسبه می شود:
که a: شیب خط خاک، RNIR: بازتاب در باند مادون قرمز نزدیک، RRED: بازتاب در باند قرمز، Fv: درصد پوشش گیاهی، RV,NIR و RV,RED بازتاب پوشش گیاهی در باند مادون قرمز نزدیک و باند قرمز می باشند که به ترتیب ضرایب ثابت 5/0 و 05/0 در نظر گرفته می شوند (غلام و همکاران، 2007).
پس از ارزیابی نقشه های حاصل از دور روش شاخص وزنی و اپراتورهای مختلف فازی، روشی که بیشتری همبستگی را با شاخص MPDI داشت انتخاب گردید و براساس آن نقشه ی نهایی خطرپذیری خشکسالی تهیه گردید.
6- نتایج و بحث
جدول شماره (2) ماتریس مقایسه زوجی تشکیل شده و جدول شماره (3) وزنهای اختصاص یافته به هر پارامتر را نشان می دهد.
جدول 2). ماتریس مقایسه ی زوجی پارامترهای موثر بر خشکسالی
تبخیر |
ارتفاع |
تراکم زهکشی |
پوشش گیاهی |
کاربری اراضی |
شیب |
جهت |
باران |
دما |
|
تبخیر |
1 |
1 |
3 |
33/0 |
2/0 |
1 |
3 |
2/0 |
1 |
ارتفاع |
1 |
1 |
3 |
33/0 |
33/0 |
33/0 |
3 |
33/0 |
3 |
تراکم زهکشی |
33/0 |
33/0 |
1 |
14/0 |
11/0 |
33/0 |
1 |
11/0 |
33/0 |
پوشش گیاهی |
3 |
3 |
7 |
1 |
14/0 |
3 |
5 |
33/0 |
5 |
کاربری اراضی |
5 |
3 |
9 |
7 |
1 |
3 |
7 |
1 |
7 |
شیب |
1 |
3 |
3 |
33/0 |
33/0 |
1 |
3 |
2/0 |
1 |
جهت |
33/0 |
33/0 |
1 |
2/0 |
14/0 |
33/0 |
1 |
14/0 |
33/0 |
باران |
5 |
3 |
9 |
3 |
1 |
5 |
7 |
1 |
7 |
دما |
1 |
33/0 |
3 |
2/0 |
14/0 |
1 |
3 |
14/0 |
1 |
جدول 3). وزن های تعیین شده به روش سلسله مراتبی برای پارامترهای موثر بر خشکسالی
تبخیر |
ارتفاع |
تراکم زهکشی |
پوشش گیاهی |
کاربری اراضی |
شیب |
جهت |
باران |
دما |
|
0561/0 |
0683/0 |
0218/0 |
1418/0 |
3051/0 |
0744/0 |
0245/0 |
2619/0 |
0461/0 |
وزن |
همانگونه که در جدول شماره ی (3) مشاهده می شود، بیشترین وزن به پارامتر کاربری اراضی اختصاص داده شده است چرا که انسان با دخالت در اکوسیستم های طبیعی بیشترین تاثیر را بر روی خشکسالی می گذارد. در واقع اراضی باغی و زراعی در مواقع خشکسالی نسبت به اراضی مرتعی حساسیت بیشتری دارند. کمترین وزن به تراکم زهکشی و جهت داده شده است. این دو پارامتر اگرچه می توانند بر روی خشکسالی تاثیر داشته باشند ولی نسبت به پارامترهایی مانند باران و کاربری اراضی تاثیر بسیار کمتری دارند.در نهایت نرخ سازگاری وزن های محاسبه شده برابر 0.082 محاسبه شد و از آنجا که این مقدار کمتر از آستانه ی قابل قبول (0.1) می باشد بنابراین صحت وزن های حاصله تایید می گردد.
نتایج حاصل از همبستگی روشهای مختلف ترکیب با شاخص MPDI در جدول شماره (4) آمده است.
جدول شماره 4- ضریب همبستگی بین روش های مختلف تلفیق و شاخص MPDI
ردیف |
روش تلفیق |
ضریب همبستگی با شاخص PDI |
1 |
شاخص وزنی |
0.81 |
2 |
Fuzzy And |
0.17 |
3 |
Fuzzy Or |
0.42 |
4 |
Fuzzy Sum(جمع) |
0.32 |
5 |
Fuzzy Product(ضرب) |
0.21 |
نتایج حاصله نشان داد که از بین روش شاخص وزنی و اپراتورهای مختلف منطق فازی، روش شاخص وزنی دارای بیشترین همبستگی(R2=0.81) با شاخص خشکسالی محاسبه شده دارد. با توجه به این جدول، اپراتورهای ضرب فازی و اپراتور حداقل، دارای کمترین همبستگی با شاخص MPDI می باشند چرا که این اپراتورها دارای بیشترین حساسیت و کمترین ریسک می باشند و حداقل ها را در نظر می گیرند ولی در دنیای واقعی، پدیده ی خشکسالی پدیده ای بسیار حساس به یک پارامتر خاص نمی باشد. به عنوان مثال اگر در منطقه ای بارندگی از حداقل آستانه ی وقوع خشکسالی کمتر باشد ولی سایر پارامترها در حد مناسبی باشند، خشکسالی اتفاق نمی افتد. اپراتورهای حداقل و ضرب، براساس پایین ترین مقدار پارامترهای وارد شده در مدل عمل می نمایند و این باعث شده است تا نتوانند به خوبی جهت تعیین مناطق در معرض خشکسالی عمل نمایند. اپراتورهای جمع و حداکثر نیز نسبت به روش شاخص وزنی از همبستگی خوبی با شاخص عمودی خشکسالی برخوردار نیستند زیرا این اپراتورها حداکثرها را در نظر می گیرند. به عبارت دیگر این اپراتورها دارای بالاترین ریسک می باشد و تنها براساس حداکثر حساسیت مربوط به یک پارامتر خاص عمل می نمایند. لذا شاخص وزنی به عنوان مناسبترین شاخص جهت تلفیق لایه های اطلاعاتی شناخته شد که نقشه ی نهایی حاصل از آن در شکل شماره (3) آمده است.
شکل 3). نقشه ی هایی خطر پذیری حوضه ی شیطور به روش شاخص وزنی
با توجه به شکل شماره (3)، اراضی باغی– زراعی که در مناطق مرکزی و جنوبی حوضه واقع شده اند (منطقه شماره 1) دارای بیشترین خطر در برابر وقوع خشکسالی می باشند. در واقع این مناطق از نظر پارامترهای دما، بارندگی، تبخیر، کاربری اراضی و ارتفاع حساسیت نسبتا بالایی به پدیده ی خشکسالی دارند و به همین دلیل پس از ترکیب عوامل فوق، این مناطق به عنوان مناطق با ریسک بالا معرفی شده اند.
سایر مناطق مرکزی و جنوبی حوضه، دارای خطر پذیری بین 40 تا 60 درصد نسبت به سایر مناطق حوضه می باشند. در این مناطق از حوضه، بارندگی حداقل، دما حداکثر، تبخیر حداکثر، ارتفاع حداقل و تراکم زهکشی حداکثر می باشد که باعث شده است تا این منطقه نسبت به مناطق شمالی و شرقی حوضه حساسیت بیشتری به خشکسالی نشان دهد. پارامترهایی تعدیل کننده ی این منطقه نسبت به نواحی با حساسیت بیش از 60 درصد، عبارتند از:
درصد پوشش گیاهی: که نسبتا کم است و حساسیت کمتری نسبت به نواحی با پوشش گیاهی بیشتر دارد.
کاربری مرتعی: که نسبت به کاربری باغی – زراعی حساسیت کمتری به خشکسالی دارد.
مناطق شمالی حوضه که در برگیرنده ی ارتفاعات می باشد و با رنگ تیره (منطقه شماره 2) نمایش داده شده اند، دارای کمترین حساسیت به پدیده ی خشکسالی می باشند زیرا در این مناطق ارتفاع و بارندگی نسبتا زیاد و دما و تبخیر نسبتا کم می باشد. ولی پارامترهایی نظیر شیب و پوشش گیاهی باعث شده اند تا این بخش از حوضه نیز تا حدی حساسیت به پدیده خشکسالی داشته باشد و در وضعیت ایده آل و بدون ریسک قرار نگیرد.
7- نتیجه گیری
در عمده تحقیقات صورت گرفته در زمینه خشکسالی در ایران از شاخص های خشکسالی مبتنی برتصاویر ماهواره ای به منظور تهیه نقشه شدت خشکسالی استفاده شده است. چنار در سال 1380 در مطالعه ای در استان های شمال غرب کشور به بررسی خشکسالی با استفاده از تصاویر ماهواره ای NOAA و شاخص NDVI در طول سال های 1998 تا 2000 پرداخت و نتایج حاصله بیانگر وقوع خشکسالی در سال های مورد مطالعه می باشد.
صفری در سال 1382 نیز به بررسی خشکسالی با استفاده از تصاویر AVHRR در منطقه زابل پرداخت که نتایج حاصله بیانگر کاهش سطح آب هامون، دبی رودخانه هیرمند و کاهش شاخص های گیاهی در طی شرایط خشکسالی می باشد.
حاصل نتایج تحقیقات فوق یک نقشه خشکسالی مربوط به یک زمان خاص می باشد که می تواند در زمینه مدیریت بحران خشکسالی کمک شایانی به مدیران و برنامه ریزان نماید. همچنین در مطالعات مذکور، تنها از تصاویر ماهواره ای به منظور بررسی خشکسالی استفاده شده و تاثیر عوامل اقلیمی و محیطی در تهیه نقشه خروجی لحاظ نشده است.
شایگان در سال 1381، در مطالعه ای به بررسی عوامل موثر بر خشکسالی و نحوه وزن دهی به این عوامل به منظور تجزیه و تحلیل ریسک پذیری خشکسالی پرداخت. ایشان در تحقیق خود پس از تعیین عوامل موثر بر خشکسالی، به مقایسه وزنهای اختصاص یافته به هر پارامتر به دو روش کریتیک و AHP پرداخت و در نهایت به این نتیجه رسید که روش AHP دقت بیشتری در تعیین وزن ها دارد.
با توجه به بررسی های نگارنده، تاکنون مطالعه ای در زمینه بررسی دقت روش های مختلف تلفیق پارامترهای موثر بر خشکسالی صورت نگرفته است. لذا در این تحقیق سعی گردید تا دقت روش های مختلف دانش مبنا در محیط GIS به منظور ترکیب عوامل موثر بر خشکسالی بررسی گردد.
نتایج تحقیق حاضر در مقایسه با نتایج حاصل از شاخص های مبتنی بر تصاویر ماهواره ای، نشان می دهد که نقشه حاصل از این تحقیق تنها بیانگر شدت خشکسالی در یک زمان خاص نمی باشد بلکه از آنجا که عوامل متعدد محیطی و اقلیمی در تهیه این نقشه دخیل داده شده است، لذا نواحی ای که در معرض خطر بیشتر خشکسالی قرار دارند در زمان وقوع خشکسالی آسیب بیشتری می بینند و لازم است با توجه به نقشه حاصله برنامه ریزی های مناسبی به منظور کاهش اثرات خشکسالی بخصوص در مناطق با ریسک بالا صورت پذیرد. به عبارت دیگر، نقشه حاصله می تواند به مدیریت ریسک خشکسالی کمک نماید و آمادگی مقابله با آن قبل از وقوع و گسترش آن را فراهم نماید، بطوریکه با آگاهی از شرایط فعلی منطقه، تصمیمات لازم و از قبل طراحی شده برای مقابله با آن اعمال گردد.
[1] . Goddard
[2] . Land Surface Temperature
[3] . Leaf water Content Index
[4] . Hunt
[5] . Anomaly Vegetation Index
[6] . Chen
[7] . Vegetation Supply Water Index
[8] . Clarson
[9] . Vegetation Condition Index
[10] . Kogan
[11] . Liang
[12] . Zhao
[13] . Normalized Difference Vegetation Index
[14] . Real Time
[15] Ghulam
[16] . Perpendicular Drought Index
[17] . Modified Perpendicular Drought Index
[18] . Fraction of Vegetation
[19] . Advanced Very High Resolution Radiometer
[20] . Advanced Visible and Near Infrared Radiometer type 2
[21] . Advanced Land Observing Satellite
[22] . Reflectance
[23] . Radiance
[24] . FLAASH
[25] . Resample
[26] . Sink
[27] . Flow direction
[28] . Flow accumulation
[29] . Analytic Hierarchy Process
[30] . Index Overlay