تعداد نشریات | 44 |
تعداد شمارهها | 1,303 |
تعداد مقالات | 16,022 |
تعداد مشاهده مقاله | 52,493,645 |
تعداد دریافت فایل اصل مقاله | 15,219,905 |
برآورد فرسایش و رسوب حوضه ی آبریز رودخانه ی گاوی با استفاده از روش MPSIAC در محیط GIS | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
هیدروژئومورفولوژی | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقاله 1، دوره 3، شماره 6، خرداد 1395، صفحه 1-16 اصل مقاله (923.72 K) | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
نویسندگان | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
داود مختاری1؛ فاطمه محمدزاده گلانی2؛ محمدرضا نیکجو3؛ شمس الله عسگری4 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1استاد گروه ژئو مورفولوژی دانشگاه تبریز | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
2کارشناسی ارشد | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
3استادیار گروه ژئومورفولوژی دانشگاه تبریز | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
4کارشناسی ارشد، مرکز تحقیقات کشاورزی و منابع طبیعی استان ایلام | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
چکیده | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
حوضهی آبریز رودخانهی گاوی واقع در جنوبغرب شهرستان ایلام با مساحت 460 کیلومتر مربع میباشد که وقوع انواع فرسایش در آن قابل ملاحظه بوده، لذا برآورد مقدار فرسایش و رسوب جهت طرحهای حفاظت آب و خاک ضروری مینماید. روش تحقیق در این مطالعه مبتنی بر مطالعات کتابخانهای، میدانی و استفاده از مدلهای تجربی برآورد فرسایش و رسوب می باشد. در این تحقیق مدل تجربی MPSIAC، که بیشترین پارامترهای مؤثر در فرسایش و رسوب را در نظر میگیرد، به عنوان تکنیک منتخب مورد استفاده قرار گرفته و پارامترهای محیطی در قالب این روش مورد بررسی قرار گرفتهاند و در نهایت نتایج حاصل با ایستگاه رسوبسنجی تنگه باجک مقایسه گردید. در این مقاله از نرمافزار GIS به منظور تهیهی لایههای اطلاعاتی استفاده شده است. نـتایج این بررسیها، بیانگر این است که مدل MPSIAC در این حوضه به خوبی پاسخ میدهد، زیرا به عنوان نمونه مقدار رسوب برآورد شده در حوضهی مورد مطالعه با مدل MPSIAC 20/453 تن رسوب در کیلومتر مربع در سال و رسوب برآورد شده در ایستگاه هیدرومتری تنگهی باجک 23/461 تن در سال میباشد، همچنین مشاهده گردید میزان رسوب برآورد شده در این حوضه در کلاس فرسایشی 4 و از نظر طبقهبندی کیفی فرسایش در کلاس درجهی زیاد قرار داشته و طبق تعریف این کلاس میتوان گفت که در این حوضه جابجایی ذرات خاک به میزانی است که اجرای برنامههای حفاظت خاک و آب ضرورت و اولویت داشته و استفاده از اراضی محدودیت زیادی دارد. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
کلیدواژهها | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
حوضه ی آبریز رودخانه ی گاوی؛ فرسایش و رسوب؛ MPSIAC؛ GIS | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
اصل مقاله | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مقدمه زمین از بدو پیدایش، در اثر نیروهای درونی و بیرونی دایماً در حال تحول و دگرگونی بوده است (رجایی، 1373). یکی از عواملی که موجب تغییر و تحول در نقاط مختلف زمین شده پدیدهی فرسایش میباشد که امروزه در تمام جهان به عنوان خطری برای رفاه بشر و بلکه حیات او شناخته شده است. امروزه تعیین بار رسوبی حوضهها برای استفاده در طرحهای توسعهی عمرانی ـ اقتصادی یک ضرورت است. معمولاً سیستم رسوبدهی یک حوضه، بسیار پیچیده و در عین حال دارای تغییرپذیری زمانی و مکانی است و همچنان به منزلهی یک چالش بزرگ پیش روی محققان در این زمینه باقی مانده است (هسو و همکاران[1]، 2012). اگرچه واژهی فرسایش در قرن نوزدهم برای عموم شناخته شده بود، اما واژهی فرسایش خاک بعداً وارد فرهنگ واژگان شد و در قرن بیستم مورد استفاده عمومی قرار گرفت (بیاتیخطیبی، 1382). از لحاظ ژئومورفولوژی حصارهای کوهستانی، محیط شکلزایی در سرزمین ایران را پر انرژی ساخته و در نتیجه به قدرت فرایندهای فرسایشی در این سرزمین افزودهاند (علاییطالقانی، 1388) در ژئومورفولوژی از فرسایش به عنوان یک فرایند نام برده میشود (شایان، 1390). فرآیند فرسایش شامل سه مرحلهی بـرداشت، حمل و رسوبگذاری است (مایر و همکاران[2]، 1969). به طور کلی فرسایش فرایندی است که طی آن ذرات خاک از بستر اصلی خود جابجا شده و به کمک یک عامل انتقال دهنده به مکانی دیگر حمل میشوند. در صورتی که عامل جدا کنندهی ذرات از بستر و انتقال آنها به مکانی دیگر آب باشد به آن فرسایش آبی گفته میشود (علیزاده، 1390). هر ساله بالغ بر 20 تا 52 میلیارد تن رسوب توسط رودخانههای جهان انتقال مییابد و در آبهای ساکن تهنشین میشود (فیضنیا و همکاران، 1381). برآورد رسوب رودخانه نیز از جمله مسایل مهم و کاربردی در مدیریت منابع آب میباشد (احمدی و همکاران، 1391). برآورد صحیح از خطر فرسایش و مشخص کردن مناطق حساس به فرسایش در بهبود توصیههای مدیریت اراضی و مقابله با فرسایش با روشهای مختلف نقش مهمی دارد (ویکتورا و همکاران، 1998). اولین گزارش نسبتاً کامل در مورد فرسایش خاک و لزوم حفاظت آب و خاک در ایران در سال 1327 توسط داون و ریبن کارشناسان فائو به زبان انگلیسی تهیه و منتشر شد. روش پسیاک برای اولین بار در ایران در سال 1352 در حوضهی آبریز سد دز استفاده شد. سپس با توجه به دقت نسبتاً خوب آن در مقایسه با سایر روشها و مدلهای تجربی در برخی از حوضههای آبریز کشور مانند دوخواهران، کهیر، زاینده رود، هلیل رود، دز، سروان، زیردان و اوزن دره مورد استفاده قرار گرفت. مومیپور (1383)، دو مدل فرسایش MPSIAC، RUSLEبا استفاده از دادهای دورسنجی و GIS در زیرحوضهی اوجان چای را مورد مقایسه قرار داده، نتیجهی نهایی مبین این واقعیت است که مدل MPSIACدر شرایط محیطی حوضه، نتایج قابل قبولتری نسبت به دیگر مدل را دارد. عاقلی کهنهشهری و صادقی (1384)، در مقالهای با مرور بر آثار مستقیم و غیرمستقیم فرسایش خاک، هزینههای فرسایش خاک را محاسبه کردند. علیزاده گرجی (1385)، برآورد فرسایش و رسوب با استفاده از مدل MPSIACدر محیط GIS حوضهی کلیجا نستاق را مورد بررسی قرار داده و به این نتیجه رسید که مدل ام پسیاک کاربرد قابل قبولی برای بررسی فرسایش و رسوب حوضهی مورد مطالعه را داشته است. منصوری دانشور و همکار (1391)، خیام و همکاران (1392)، در مقالهای تحت عـنوان مقایسهی کارایی مدلهای MPSIAC و EP، به این نتیجه رسیدند که مدل MPSIAC با دقت 31/98 درصد رسوب حوضه را برآورد نموده است این در حالی است که مدل EPMتنها 4/75 درصد رسوب مشاهدهای حوضه را با متوسط شدت تولید رسوب 5133/111 تن بر هکتار در سال محاسبه نموده است. محاسبهی میزان فرسایش در یک حوضهی آبریز در صورت موجود بودن آمار کافی از دبی و رسوب، با بکارگیری روشهای آماری امکانپذیر میباشد. ولی نبود یا کمبود دادهها در زمینهی فرسایش خاک و تولید رسوب در بسیاری از حوضههای کشور کاربرد روشهای تجربی را برای برآورد فرسایش خاک الزامی میکند. بر این اساس، روش تجربی ام پسیاک با لحاظ کردن نه عامل مؤثر در فرسایش و تولید رسوب شامل زمین شناسی، خاک، اقلیم، رواناب، شیب، پوشش گیاهی زمین، کاربری اراضی، وضعیت فعلی فرسایش و فرسایش رودخانهای و حمل رسوب بیش از سایر مدلها در ایران مورد توجه قرار میگیرد. در این پژوهش، نتایج مدل ام پسایک با نتایج ایستگاه رسوبسنجی تنگهی باجک مقایسه شده است. هدف نهایی از این تحقیق بررسی میزان کارایی و دقت مدل تجربی MPSIAC در برآورد میزان رسوب حوضه و شناسایی عوامل مؤثرتر در فرسایش خاک و تولید رسوب حوضهی مورد مطالعه میباشد. لازم به ذکر میباشد که این پژوهش برای نخستین بار حوضهی مورد مطالعه را مورد بررسی قرار داده است. فرضیاتی که در این ارتباط مطرح شده، لیتولوژی بیشترین نقش را در ایـجاد فرسایش حوضهی آبریز دارد و با توجه به ویـژگیهای طبیعی منطقه مدل MPSIAC میتواند برای برآورد فرسایش و رسوب حوضهی مورد مطالعه استفاده شود. منطقهی مورد مطالعه حوضهی آبریز رودخانهی گاوی یکی از رودخانههای حوضهی آبریز مهران در جنوبغرب شهرستان ایلام با موقعیت جغرافیایی ¢24°46 تا ¢36°46 طول شرقی و ¢12°33 تا ¢28°33 عرض شمالی واقع گردیده است (شکل 1). این حوضه بین حوضهی آبریز رودخانههای چنگوله و کنجانچم واقع گردیده است. کمترین ارتفاع حوضهی مورد مطالعه 409 متر در تنگ باجک محل احداث سد گاوی و بیشترین ارتفاع نیز در کوه ورزرین 2714 متر واقع در شمال حوضه میباشد. با توجه به آمار و اطلاعات 15سال (1377-1392) ایستگاههای بارانسنجی، بیشینه بارش سالانه در محدودهی مورد مطالعه، بین 515 میلیمتر در ایستگاه کشوری تا 376 میلیمتر در ایستگاه پیرمحمد متغیر است. با توجه به تقسیمبندی که نبوی در سال 1355 انجام داده و پهنهی ایران را به 16 زون تقسیم کردهاند (جداریعیوضی، 1376)، حوضهی آبریز مورد مطالعه در زاگرس چین خورده قرار میگیرد. این زون که کوههای زاگرس را در بر میگیرد، ساخت زمینشناسی آن ساده، ملایم و شامل مجموعهای از رشته تاقدیسهای نزدیک و به هم فشرده با سطح محدودی معمولاً قایم و روند شمال غربی، جنوب شرق است. رسوبات چینخوردهی این منطقه به طور متناوب از آهک یا دولومیت همراه با مارن و مارنهای آهکی است که با چینهبندی کم و بیش ظریف مشخصاند. سنگهای مارنی موجود در سطح حوضه بسیار ضعیف در مقابل عوامل فرسایش بوده و سالانه حجم قابل توجهی از رسوبات ناشی از فرسایش این سنگها به همراه سیلابها از این حوضه خارج میشود.
شکل (1) موقعیت حوضهی آبریز رودخانهی گاوی مواد و روش روش مطالعه در این پژوهش شامل برآورد فرسایش و رسوبزایی با استفاده از روش ام پسیاک و سپس مقایسهی نتایج این روش تجربی با برآورد واقعی حاصل شده از طریق دادههای ایستگاه رسوبسنجی منطقهی مورد مطالعه میباشد. دادههای مورد نیاز در این مطالعه شامل، نقشههای توپوگرافی، تصاویر ماهوارهای Land Sat ETM 2014 جهت تهیهی نقشهی پوشش زمین و کاربری اراضی، بازدیدهای صحرایی (جهت بررسی وضعیت فعلی فرسایش و فرسایش رودخانهای) و دادههای آب و هواشناسی سالهای 1377تا 1392 و مطالعات پیشین است. از مطالعات پیشین جهت بررسی پوشش گیاهی و محاسبهی ضریب فرسایشپذیری خاک استفاده شده است. به منظور بالا بردن سرعت و دقت در محاسبات و تهیهی لایههای اطلاعاتی روش MPSIAC، از دانش سنجش از دور و (نرمافزار ArcGIS) جهت ژئورفرنس کردن و رقومی کردن نقشهها و تهیهی نقشههای ثانوی از قبیل نقشههای پوشش زمین، کاربری اراضی، زمینشناسی و شیب مربوط به حوضه استفاده شده است. در روش مورد مطالعه، تأثیر و نقش نه عامل مهم و مؤثر در فرسایش و تولید رسوب بر اساس جدول استاندارد و رابطههای تعریف شده در مدل، محاسبه شده و سپس نتایج به صورت لایههای اطلاعاتی با فرمت رستر وارد محیط سیستم اطلاعات جغرافیایی گردیده است. زمینشناسی (y1) ـ در روش ام پسیاک امتیاز عامل زمینشناسی از رابطهی (1)، جدول (1)، حاصل میشود که بر اساس نوع سنگ، سختی، شکستگی و هوازدگی تعیین میشود. لذا با توجه به مطالعات زمینشناسی حوضهی آبریز رودخانهی گاوی بر اساس سازندهای موجود در این حوضه (شکل2) و همچنین مقاومت آنها در مقابل فرسایش، هر یک از سازندها امتیازدهی شدهاند که کمترین امتیاز در آنها به سازند سروک و بیشترین امتیاز به سازند گورپی نسبت داده شده است (جدول 3). خاکشناسی (Y2)ـ در روش MPSIAC بـرای تعیین عامـل خاک از رابـطهی (2)، جـدول (1)، استفاده میشود: که در آنX2 امتیاز رسوبدهی خاک در روش پسیاک و k عامل فرسایشپذیری خاک در فرمول جهانی فرسایش میباشد. در فرمول جهانی فرسایش[3](USLE) برای تعیین k از مشخصههای درصد سیلت به اضافهی شن خیلی ریز، درصد شن، درصد مادهی آلی، ساختمان خاک و قابلیت نفوذ استفاده به عمل میآید (رفاهی، 1388). به منظور تعیین مشخصههای درصدی جهت تعیین این عامل از گزارش خاکشناسی حوضه (آنالیز فیزیکی شیمیایی پروفیل خاک حوضهی آبریز مهران) استفاده شده است. آب و هوا (Y3)ـ جهت تعیین امتیاز عامل آب و هوا با توجه به رابطهی (3)، جدول (1)، نیاز به مقدار بارش 6 ساعته با دوره بازگشت 2 سال بر حسب میلیمتر میباشد. لذا از روش بارش 24 ساعته با استفاده از فرمول تجربی زیر میتوان حداکثر بارش 6 ساعته با دوره برگشت T یا را تخمین زد. معادلهی (1) در این معادله t برحسب ساعت و T دوره برگشت بر حسب سال و بر حسب میلیمتر میباشد (علیزاده،1390). برای بدست آوردن بارش 24 ساعته با دور برگشت 2 سال جهت تعیین امتیاز این عامل در حوضهی مورد مطالعه از آمار ایستگاه سینوپتیک ایلام و مهران استفاده شده است. رواناب (y4) ـ عامل رواناب به عنوان چهارمین پارامتر مندرج در مدل پسیاک اصلاح شده عبارت است از «بازده یک سیستم آبخیز که از نتیجه، عملکرد ساختمان آبخیز بر روی دادههای آن (نزولات آسمانی) پدیدار میگردد (رفاهی، 1388). برای تعیین امتیاز عامل رواناب در روش MPSIAC از رابطهی (4)، جدول (1)، استفاده شده است که در این رابطه R: ارتفاع رواناب سالانه بر حسب میلیمتر، و q دبی ویژه پیک برحسب متر مکعب بر ثانیه در کیلومتر مربع میباشد. با استفاده از مطالعات هیدرولوژی و منابع آب حوضهی مورد مطالعه پس از مشخص شدن ارتفاع رواناب و دبی حداکثر سالانه در هر یک از واحدهای کاری نسبت به محاسبهی عامل رواناب در روش پسیاک اقدام گردید. جدول (2)، این محاسبات را نشان میدهد. توپوگرافی (Y5)ـ شیب حوضه نقش اساسی در میزان رواناب، مقدار نفوذ، شدت سیلابها و میزان فرسایش دارد (مهدوی، 1384). در این مطالعه جهت استخراج پارامتر شیب، ابتدا نقشههای توپوگرافی رقومی شده و سپس این خطوط در محیط GISبه نقشهی مدل رقومی ارتفاعی DEM تبدیل شدهاند و آنگاه از این نقشه با دستور Slope در محیط نرمافزار، نقشه درصد متوسط شیب استخراج گردید (شکل3)، سپس امتیاز این عامل تعیین گردیده است. پوشش زمین (y6)ـ در مدل پسیاک اصلاح شده از رابطهی (6)، جدول (1)، جهت تعیین امتیاز عامل پوشش زمین استفاده گردیده است. در تحقیق حاضر، به منظور تعیین عامل فوق، ابتدا با استفاده از تصویر ماهوارهای لندست از نوع TM[4] اقدام به تهیهی نقشهی پوشش زمین (شکل4)، از نظر مقدار گیاهی گردید و سپس با مطالعهی میدانی و مشاهدات صحرایی تفاوت واحدهای گیاهی مشخص و مرز تیپها اصلاح شده است و درصد اراضی لخت و فاقد پوشش هر واحد محاسبه شد و در نهایت عدد به دست آمده را در رابطهی مربوط قرار داده و امتیاز عامل پوشش زمین برآورد شده است. نحوهی استفاده از زمین (Y7)ـ در روش ام پسیاک بر مبنای درصد تاج پوشش گیاهی، مشخصهی نحوهی استفاده از زمین محاسبه میگردد. به منظور تعیین امتیاز عامل نحوهی استفاده از اراضی، ابتدا اقدام به تهیهی نقشهی کاربری اراضی (شکل 5) و سپس نقشهی پوشش گیاهی حوضهی آبریز مورد مطالعه با استفاده از تصویر ماهوارهای لندست از نوع TM گردید، سپس مساحت هر تیپ از پوشش گیاهی مشخص شده و در ادامه جهت تعیین وزن برای هر تیپ پوشش گیاهی از نتایج گزارش مطالعات پوشش گیاهی انجام شــده در استان استفاده شده است. در نهایت درصد تاج پوشش گیاهی هر واحد محاسبه و با استفاده از رابطهی (7)، جدول (1)، امتیاز عامل مربوطه تعیین میگردد. جدول (1) عوامل مؤثر در مدل MPSIAC و نحوهی امتیازدهی به آنها
منبع: (رفاهی، 1388) وضعیت فعلی فرسایش در سطح حوضهی (y8)ـ در روش MPSIACبرای تعیین امتیاز عامل وضعیت فعلی فرسایش، از رابطهی (8)، (جدول 1)، استفاده میشود. تعیین این عامل در حوضهی مورد مطالعه با کمک اطلاعات مربوط به دادههای صحرایی، استفاده از نرمافزاز گوگل ارث جهت بررسی مناطق دور از دسترس و با توجه به نوع فرسایش در سطح حوضهی آبریز با استفاده از جدول BLM[5] صورت گرفته است. پس از بررسی 7 عامل، وضعیت سـطح خاک و ارزیابی کـمی، مـجموع نـمرههای داده شـده را با هـم جمع کـرده، سپس بر مجموع حداکثر نمرات هر عامل تقسیم نموده و حاصل را در 100 ضرب میکنیم. جدول (2) مشخصات هیدرولوژیکی حوضهی مورد مطالعه بر حسب زیرحوضه
جدول (3) نحوهی امتیازدهی به عامل زمینشناسی در مدل MPSIAC
منبع: (منابع طبیعی استان ایلام، 1390) فرسایش رودخانهای (Y9)ـ برای تعیین امتیاز این عامل از رابطهی (9)، (جدول 1)، استفاده شده است. به منظور تعیین امتیاز فرسایش رودخانهای، SSf[6] به عنوان یکی از پارامترهای جدول BLM مد نظر قرار گرفته و در نهایت به صورت لایهای اطلاعاتی با فرمت رستری وارد محیط GISگردیده است.
شکل (2) نقشهی زمینشناسی حوضهی آبریز رودخانهی گاوی
شکل (3) نقشهی شیب حوضه، حوضهی آبریز رودخانهی گاوی
شکل (4) نقشهی پوشش زمین حوضهی آبریز رودخانهی گاوی
شکل (5) نقشهی کاربری اراضی حوضهی آبریز رودخانهی گاوی تعیین درجهی رسوبدهی واحدهای هیدرولوِژیکی پس از وزندهی و نهاییسازی لایههای اطلاعاتی که به صورت لایههای وزنی درآمدهاند برای تهیهی نقشهی فرسایش حوضه (شکل 6)، لازم است این عوامل را با هم تلفیق کنیم، لازم به ذکر است لایههای وزنی باید به صورت رستری باشند تا امکان تلفیق آنها وجود داشته باشد. سپس به منظور تعیین میزان فرسایش هر یک از واحدهای هیدرولوژیکی، شناسایی تأثیرگذارترین عامل در تولید رسوب و نیز تشخیص واحدهای بحرانی از نظر تولید رسوب، محاسبهی درجهی رسوبدهی هر یک از این واحدها ضروری است. برای این منظور مجموع امتیاز عوامل نه گانه در هر واحد محاسبه و درجهی رسوبدهی آن واحد مشخص میگردد (جدول 4).
شکل(6) نقشهی شدت فرسایش حوضهی آبریز رودخانهی گاوی به منظور بررسی بیشتر و تعیین میزان فرسایش و رسوبزایی در سطح حوضه از جدول (5)، که در این زمینه تهیه شده است، استفاده میشود. در نهایت برای اعمال دقت بیشتر و پرهیز از اشتباه در عمل درونیابی و برونیابی برای برآورد فرسایش و تولید رسوب طبق مدل ام پسیاک از رابطه بین درجهی رسوبدهی و میزان تولید رسوب، رابطهی (10)، استفاده به عمل میآید (صفامنش و همکاران[7]، 2006): رابطهی (10) 18.6=Qs در این رابطه QSمیزان رسوبدهی سالانه بر حسب متر مکعب در کیلومتر مربع،e عدد نپرین که تقریباً برابر با 718/2 میباشد وR درجهی رسوبدهی یعنی مجموع امتیازات عوامل مختلف در نظر گرفته شده در مدل MPSIACمیباشد. لازم به ذکر میباشد که وزن مخصوص برای حوضهی مورد مطالعه با توجه به مطالعات انجام شده توسط امور آب استان بر روی رودخانه، 5/1 در نظر گرفته میشود. برای محاسبهیSDR از روش لارنس[8] رابطهی (11) استفاده شده است و برای تعیین فرسایش کل بایستی فرسایش ویژه را در کل مساحت هر زیر حوضه ضرب کنیم (افسری و قدوسی 1390). رابطهی (11) SDR = جدول (4) دادههای لازم جهت محاسبه درجهی رسوبدهی واحدهای هیدرولوژیکی حوضهی آبریز رودخانهی گاوی
جدول (5) تعیین میزان تولید رسوب سالانه و کلاس فرسایش خاک در روش MPSIAC
منبع: (رفاهی، 1385) بحث و نتایج این پژوهش با هدف برآورد فرسایش و رسوبزایی حوضهی آبریز رودخانهی گاوی به روش ام پسایک و مقایسهی نتایج با برآورد واقعی انجام شده است. نتایج بررسی روش ام پسیاک و 9 پارامتر این روش برای 4 زیرحوضه تعیین شده (جدول 6)، نشان میدهد که با توجه به مجموع امتیازات پارامترهای مؤثر در فرسایش و رسوب، حوضهی مورد مطالعه بـه صورت کلی با درجهی رسوبدهی 47/77 در کلاس فرسایشی زیـاد قرار میگیرد و میزان تولید رسـوب این کلاس 1429-476 متر مکعب در کیلومتر مربع و 5/2143-714 تن در کیلومتر مربع میباشد. با اعمال رابطهی (10)، میزان رسوبدهی سالانه برابر 13/302 متر مکعب در کیلومتر مربع برآورد گردیده که با دخالت دادن وزن مخصوص (5/1)، معادل 20/453 تن رسوب در کیلومتر مربع است. جدول (6) نتایج مدل MPSIACدر حوضه
نتیجهگیری نتایج حاصله نشان مـیدهد که منطقهی مورد مطالعه از لـحاظ فرسایشپذیری در ردهی فرسایشی زیـاد قـرار میگیرد و طبق تعریف این کلاس میتوان گفت که در این حوضه جابجایی ذرات خاک به میزانی است که اجرای برنامههای حفاظت خاک و آب ضرورت و اولویت داشته و استفاده از اراضی محدودیت زیادی دارد. همچنین مدل MPSIAC در این حوضه به خوبی پاسخ میدهد و فرضیهی اول تأیید میشود، زیرا به عنوان نمونه مقدار رسوب برآورد شده در حوضهی مورد مطالعه با مدل ام پسیاک 20/453 تن رسوب در کیلومتر مربع در سال و رسوب برآورد شده در ایستگاه رسوبسنجی حوضهی مربوطه 23/461 تن رسوب در کیلومتر مربع در سال میباشد، دلیل دقت این روش میتواند این باشد که این روش به عنوان جغرافیاییترین مدل در برگیرندهی 9 پارامتر تأثیرگذار در تولید فرسایش و رسوب میباشد. در مورد فرضیهی دوم قبل از انجام تحقیق با توجه به بازدیدهای به عمل آمده از منطقه مورد مطالعه به نظر میرسید که از بین عوامل مؤثر در فرسایش، لیتولوژی با توجه به سازندهای حساس به فرسایش، بیشترین تأثیر را در فرسایش حوضه داشته باشد، در حالی که نتایج حاصل شده از طریق مدل تجربی عامل کاربری با امتیاز 48/18را به عنوان عوامل مؤثر در فرسایش حوضه در اولویت اول و عامل زمینشناسی با امتیاز 76/5 در اولویت ششم قرار میدهد. بر این اساس فرضیهی فوق رد میشود. در پایان باید اشاره کرد که با توجه به اینکه حوضه از نظر تقسیمبندی واحد اراضی در گروه دشتهای دامنهای با شیب ملایم قرار میگیرد و مناسب کشت غلات به صورت دیم و مرتع آن محل قشلاق عشایر میباشد، در نتیجه این شرایط تأثیر عامل کاربری اراضی را تأیید میکند. [1]- Hsu et al., [2]- Mayer et al., 1- Unirersal Soil Loss Equation 1- Thematic Mapper [5]- Burau of Land Management [6]- Soil Surface Factor [7]- Safamanesh et al., 2006 2- Larens مقدمه زمین از بدو پیدایش، در اثر نیروهای درونی و بیرونی دایماً در حال تحول و دگرگونی بوده است (رجایی، 1373). یکی از عواملی که موجب تغییر و تحول در نقاط مختلف زمین شده پدیدهی فرسایش میباشد که امروزه در تمام جهان به عنوان خطری برای رفاه بشر و بلکه حیات او شناخته شده است. امروزه تعیین بار رسوبی حوضهها برای استفاده در طرحهای توسعهی عمرانی ـ اقتصادی یک ضرورت است. معمولاً سیستم رسوبدهی یک حوضه، بسیار پیچیده و در عین حال دارای تغییرپذیری زمانی و مکانی است و همچنان به منزلهی یک چالش بزرگ پیش روی محققان در این زمینه باقی مانده است (هسو و همکاران[1]، 2012). اگرچه واژهی فرسایش در قرن نوزدهم برای عموم شناخته شده بود، اما واژهی فرسایش خاک بعداً وارد فرهنگ واژگان شد و در قرن بیستم مورد استفاده عمومی قرار گرفت (بیاتیخطیبی، 1382). از لحاظ ژئومورفولوژی حصارهای کوهستانی، محیط شکلزایی در سرزمین ایران را پر انرژی ساخته و در نتیجه به قدرت فرایندهای فرسایشی در این سرزمین افزودهاند (علاییطالقانی، 1388) در ژئومورفولوژی از فرسایش به عنوان یک فرایند نام برده میشود (شایان، 1390). فرآیند فرسایش شامل سه مرحلهی بـرداشت، حمل و رسوبگذاری است (مایر و همکاران[2]، 1969). به طور کلی فرسایش فرایندی است که طی آن ذرات خاک از بستر اصلی خود جابجا شده و به کمک یک عامل انتقال دهنده به مکانی دیگر حمل میشوند. در صورتی که عامل جدا کنندهی ذرات از بستر و انتقال آنها به مکانی دیگر آب باشد به آن فرسایش آبی گفته میشود (علیزاده، 1390). هر ساله بالغ بر 20 تا 52 میلیارد تن رسوب توسط رودخانههای جهان انتقال مییابد و در آبهای ساکن تهنشین میشود (فیضنیا و همکاران، 1381). برآورد رسوب رودخانه نیز از جمله مسایل مهم و کاربردی در مدیریت منابع آب میباشد (احمدی و همکاران، 1391). برآورد صحیح از خطر فرسایش و مشخص کردن مناطق حساس به فرسایش در بهبود توصیههای مدیریت اراضی و مقابله با فرسایش با روشهای مختلف نقش مهمی دارد (ویکتورا و همکاران، 1998). اولین گزارش نسبتاً کامل در مورد فرسایش خاک و لزوم حفاظت آب و خاک در ایران در سال 1327 توسط داون و ریبن کارشناسان فائو به زبان انگلیسی تهیه و منتشر شد. روش پسیاک برای اولین بار در ایران در سال 1352 در حوضهی آبریز سد دز استفاده شد. سپس با توجه به دقت نسبتاً خوب آن در مقایسه با سایر روشها و مدلهای تجربی در برخی از حوضههای آبریز کشور مانند دوخواهران، کهیر، زاینده رود، هلیل رود، دز، سروان، زیردان و اوزن دره مورد استفاده قرار گرفت. مومیپور (1383)، دو مدل فرسایش MPSIAC، RUSLEبا استفاده از دادهای دورسنجی و GIS در زیرحوضهی اوجان چای را مورد مقایسه قرار داده، نتیجهی نهایی مبین این واقعیت است که مدل MPSIACدر شرایط محیطی حوضه، نتایج قابل قبولتری نسبت به دیگر مدل را دارد. عاقلی کهنهشهری و صادقی (1384)، در مقالهای با مرور بر آثار مستقیم و غیرمستقیم فرسایش خاک، هزینههای فرسایش خاک را محاسبه کردند. علیزاده گرجی (1385)، برآورد فرسایش و رسوب با استفاده از مدل MPSIACدر محیط GIS حوضهی کلیجا نستاق را مورد بررسی قرار داده و به این نتیجه رسید که مدل ام پسیاک کاربرد قابل قبولی برای بررسی فرسایش و رسوب حوضهی مورد مطالعه را داشته است. منصوری دانشور و همکار (1391)، خیام و همکاران (1392)، در مقالهای تحت عـنوان مقایسهی کارایی مدلهای MPSIAC و EP، به این نتیجه رسیدند که مدل MPSIAC با دقت 31/98 درصد رسوب حوضه را برآورد نموده است این در حالی است که مدل EPMتنها 4/75 درصد رسوب مشاهدهای حوضه را با متوسط شدت تولید رسوب 5133/111 تن بر هکتار در سال محاسبه نموده است. محاسبهی میزان فرسایش در یک حوضهی آبریز در صورت موجود بودن آمار کافی از دبی و رسوب، با بکارگیری روشهای آماری امکانپذیر میباشد. ولی نبود یا کمبود دادهها در زمینهی فرسایش خاک و تولید رسوب در بسیاری از حوضههای کشور کاربرد روشهای تجربی را برای برآورد فرسایش خاک الزامی میکند. بر این اساس، روش تجربی ام پسیاک با لحاظ کردن نه عامل مؤثر در فرسایش و تولید رسوب شامل زمین شناسی، خاک، اقلیم، رواناب، شیب، پوشش گیاهی زمین، کاربری اراضی، وضعیت فعلی فرسایش و فرسایش رودخانهای و حمل رسوب بیش از سایر مدلها در ایران مورد توجه قرار میگیرد. در این پژوهش، نتایج مدل ام پسایک با نتایج ایستگاه رسوبسنجی تنگهی باجک مقایسه شده است. هدف نهایی از این تحقیق بررسی میزان کارایی و دقت مدل تجربی MPSIAC در برآورد میزان رسوب حوضه و شناسایی عوامل مؤثرتر در فرسایش خاک و تولید رسوب حوضهی مورد مطالعه میباشد. لازم به ذکر میباشد که این پژوهش برای نخستین بار حوضهی مورد مطالعه را مورد بررسی قرار داده است. فرضیاتی که در این ارتباط مطرح شده، لیتولوژی بیشترین نقش را در ایـجاد فرسایش حوضهی آبریز دارد و با توجه به ویـژگیهای طبیعی منطقه مدل MPSIAC میتواند برای برآورد فرسایش و رسوب حوضهی مورد مطالعه استفاده شود. منطقهی مورد مطالعه حوضهی آبریز رودخانهی گاوی یکی از رودخانههای حوضهی آبریز مهران در جنوبغرب شهرستان ایلام با موقعیت جغرافیایی ¢24°46 تا ¢36°46 طول شرقی و ¢12°33 تا ¢28°33 عرض شمالی واقع گردیده است (شکل 1). این حوضه بین حوضهی آبریز رودخانههای چنگوله و کنجانچم واقع گردیده است. کمترین ارتفاع حوضهی مورد مطالعه 409 متر در تنگ باجک محل احداث سد گاوی و بیشترین ارتفاع نیز در کوه ورزرین 2714 متر واقع در شمال حوضه میباشد. با توجه به آمار و اطلاعات 15سال (1377-1392) ایستگاههای بارانسنجی، بیشینه بارش سالانه در محدودهی مورد مطالعه، بین 515 میلیمتر در ایستگاه کشوری تا 376 میلیمتر در ایستگاه پیرمحمد متغیر است. با توجه به تقسیمبندی که نبوی در سال 1355 انجام داده و پهنهی ایران را به 16 زون تقسیم کردهاند (جداریعیوضی، 1376)، حوضهی آبریز مورد مطالعه در زاگرس چین خورده قرار میگیرد. این زون که کوههای زاگرس را در بر میگیرد، ساخت زمینشناسی آن ساده، ملایم و شامل مجموعهای از رشته تاقدیسهای نزدیک و به هم فشرده با سطح محدودی معمولاً قایم و روند شمال غربی، جنوب شرق است. رسوبات چینخوردهی این منطقه به طور متناوب از آهک یا دولومیت همراه با مارن و مارنهای آهکی است که با چینهبندی کم و بیش ظریف مشخصاند. سنگهای مارنی موجود در سطح حوضه بسیار ضعیف در مقابل عوامل فرسایش بوده و سالانه حجم قابل توجهی از رسوبات ناشی از فرسایش این سنگها به همراه سیلابها از این حوضه خارج میشود.
شکل (1) موقعیت حوضهی آبریز رودخانهی گاوی مواد و روش روش مطالعه در این پژوهش شامل برآورد فرسایش و رسوبزایی با استفاده از روش ام پسیاک و سپس مقایسهی نتایج این روش تجربی با برآورد واقعی حاصل شده از طریق دادههای ایستگاه رسوبسنجی منطقهی مورد مطالعه میباشد. دادههای مورد نیاز در این مطالعه شامل، نقشههای توپوگرافی، تصاویر ماهوارهای Land Sat ETM 2014 جهت تهیهی نقشهی پوشش زمین و کاربری اراضی، بازدیدهای صحرایی (جهت بررسی وضعیت فعلی فرسایش و فرسایش رودخانهای) و دادههای آب و هواشناسی سالهای 1377تا 1392 و مطالعات پیشین است. از مطالعات پیشین جهت بررسی پوشش گیاهی و محاسبهی ضریب فرسایشپذیری خاک استفاده شده است. به منظور بالا بردن سرعت و دقت در محاسبات و تهیهی لایههای اطلاعاتی روش MPSIAC، از دانش سنجش از دور و (نرمافزار ArcGIS) جهت ژئورفرنس کردن و رقومی کردن نقشهها و تهیهی نقشههای ثانوی از قبیل نقشههای پوشش زمین، کاربری اراضی، زمینشناسی و شیب مربوط به حوضه استفاده شده است. در روش مورد مطالعه، تأثیر و نقش نه عامل مهم و مؤثر در فرسایش و تولید رسوب بر اساس جدول استاندارد و رابطههای تعریف شده در مدل، محاسبه شده و سپس نتایج به صورت لایههای اطلاعاتی با فرمت رستر وارد محیط سیستم اطلاعات جغرافیایی گردیده است. زمینشناسی (y1) ـ در روش ام پسیاک امتیاز عامل زمینشناسی از رابطهی (1)، جدول (1)، حاصل میشود که بر اساس نوع سنگ، سختی، شکستگی و هوازدگی تعیین میشود. لذا با توجه به مطالعات زمینشناسی حوضهی آبریز رودخانهی گاوی بر اساس سازندهای موجود در این حوضه (شکل2) و همچنین مقاومت آنها در مقابل فرسایش، هر یک از سازندها امتیازدهی شدهاند که کمترین امتیاز در آنها به سازند سروک و بیشترین امتیاز به سازند گورپی نسبت داده شده است (جدول 3). خاکشناسی (Y2)ـ در روش MPSIAC بـرای تعیین عامـل خاک از رابـطهی (2)، جـدول (1)، استفاده میشود: که در آنX2 امتیاز رسوبدهی خاک در روش پسیاک و k عامل فرسایشپذیری خاک در فرمول جهانی فرسایش میباشد. در فرمول جهانی فرسایش[3](USLE) برای تعیین k از مشخصههای درصد سیلت به اضافهی شن خیلی ریز، درصد شن، درصد مادهی آلی، ساختمان خاک و قابلیت نفوذ استفاده به عمل میآید (رفاهی، 1388). به منظور تعیین مشخصههای درصدی جهت تعیین این عامل از گزارش خاکشناسی حوضه (آنالیز فیزیکی شیمیایی پروفیل خاک حوضهی آبریز مهران) استفاده شده است. آب و هوا (Y3)ـ جهت تعیین امتیاز عامل آب و هوا با توجه به رابطهی (3)، جدول (1)، نیاز به مقدار بارش 6 ساعته با دوره بازگشت 2 سال بر حسب میلیمتر میباشد. لذا از روش بارش 24 ساعته با استفاده از فرمول تجربی زیر میتوان حداکثر بارش 6 ساعته با دوره برگشت T یا را تخمین زد. معادلهی (1) در این معادله t برحسب ساعت و T دوره برگشت بر حسب سال و بر حسب میلیمتر میباشد (علیزاده،1390). برای بدست آوردن بارش 24 ساعته با دور برگشت 2 سال جهت تعیین امتیاز این عامل در حوضهی مورد مطالعه از آمار ایستگاه سینوپتیک ایلام و مهران استفاده شده است. رواناب (y4) ـ عامل رواناب به عنوان چهارمین پارامتر مندرج در مدل پسیاک اصلاح شده عبارت است از «بازده یک سیستم آبخیز که از نتیجه، عملکرد ساختمان آبخیز بر روی دادههای آن (نزولات آسمانی) پدیدار میگردد (رفاهی، 1388). برای تعیین امتیاز عامل رواناب در روش MPSIAC از رابطهی (4)، جدول (1)، استفاده شده است که در این رابطه R: ارتفاع رواناب سالانه بر حسب میلیمتر، و q دبی ویژه پیک برحسب متر مکعب بر ثانیه در کیلومتر مربع میباشد. با استفاده از مطالعات هیدرولوژی و منابع آب حوضهی مورد مطالعه پس از مشخص شدن ارتفاع رواناب و دبی حداکثر سالانه در هر یک از واحدهای کاری نسبت به محاسبهی عامل رواناب در روش پسیاک اقدام گردید. جدول (2)، این محاسبات را نشان میدهد. توپوگرافی (Y5)ـ شیب حوضه نقش اساسی در میزان رواناب، مقدار نفوذ، شدت سیلابها و میزان فرسایش دارد (مهدوی، 1384). در این مطالعه جهت استخراج پارامتر شیب، ابتدا نقشههای توپوگرافی رقومی شده و سپس این خطوط در محیط GISبه نقشهی مدل رقومی ارتفاعی DEM تبدیل شدهاند و آنگاه از این نقشه با دستور Slope در محیط نرمافزار، نقشه درصد متوسط شیب استخراج گردید (شکل3)، سپس امتیاز این عامل تعیین گردیده است. پوشش زمین (y6)ـ در مدل پسیاک اصلاح شده از رابطهی (6)، جدول (1)، جهت تعیین امتیاز عامل پوشش زمین استفاده گردیده است. در تحقیق حاضر، به منظور تعیین عامل فوق، ابتدا با استفاده از تصویر ماهوارهای لندست از نوع TM[4] اقدام به تهیهی نقشهی پوشش زمین (شکل4)، از نظر مقدار گیاهی گردید و سپس با مطالعهی میدانی و مشاهدات صحرایی تفاوت واحدهای گیاهی مشخص و مرز تیپها اصلاح شده است و درصد اراضی لخت و فاقد پوشش هر واحد محاسبه شد و در نهایت عدد به دست آمده را در رابطهی مربوط قرار داده و امتیاز عامل پوشش زمین برآورد شده است. نحوهی استفاده از زمین (Y7)ـ در روش ام پسیاک بر مبنای درصد تاج پوشش گیاهی، مشخصهی نحوهی استفاده از زمین محاسبه میگردد. به منظور تعیین امتیاز عامل نحوهی استفاده از اراضی، ابتدا اقدام به تهیهی نقشهی کاربری اراضی (شکل 5) و سپس نقشهی پوشش گیاهی حوضهی آبریز مورد مطالعه با استفاده از تصویر ماهوارهای لندست از نوع TM گردید، سپس مساحت هر تیپ از پوشش گیاهی مشخص شده و در ادامه جهت تعیین وزن برای هر تیپ پوشش گیاهی از نتایج گزارش مطالعات پوشش گیاهی انجام شــده در استان استفاده شده است. در نهایت درصد تاج پوشش گیاهی هر واحد محاسبه و با استفاده از رابطهی (7)، جدول (1)، امتیاز عامل مربوطه تعیین میگردد. جدول (1) عوامل مؤثر در مدل MPSIAC و نحوهی امتیازدهی به آنها
منبع: (رفاهی، 1388) وضعیت فعلی فرسایش در سطح حوضهی (y8)ـ در روش MPSIACبرای تعیین امتیاز عامل وضعیت فعلی فرسایش، از رابطهی (8)، (جدول 1)، استفاده میشود. تعیین این عامل در حوضهی مورد مطالعه با کمک اطلاعات مربوط به دادههای صحرایی، استفاده از نرمافزاز گوگل ارث جهت بررسی مناطق دور از دسترس و با توجه به نوع فرسایش در سطح حوضهی آبریز با استفاده از جدول BLM[5] صورت گرفته است. پس از بررسی 7 عامل، وضعیت سـطح خاک و ارزیابی کـمی، مـجموع نـمرههای داده شـده را با هـم جمع کـرده، سپس بر مجموع حداکثر نمرات هر عامل تقسیم نموده و حاصل را در 100 ضرب میکنیم. جدول (2) مشخصات هیدرولوژیکی حوضهی مورد مطالعه بر حسب زیرحوضه
جدول (3) نحوهی امتیازدهی به عامل زمینشناسی در مدل MPSIAC
منبع: (منابع طبیعی استان ایلام، 1390) فرسایش رودخانهای (Y9)ـ برای تعیین امتیاز این عامل از رابطهی (9)، (جدول 1)، استفاده شده است. به منظور تعیین امتیاز فرسایش رودخانهای، SSf[6] به عنوان یکی از پارامترهای جدول BLM مد نظر قرار گرفته و در نهایت به صورت لایهای اطلاعاتی با فرمت رستری وارد محیط GISگردیده است.
شکل (2) نقشهی زمینشناسی حوضهی آبریز رودخانهی گاوی
شکل (3) نقشهی شیب حوضه، حوضهی آبریز رودخانهی گاوی
شکل (4) نقشهی پوشش زمین حوضهی آبریز رودخانهی گاوی
شکل (5) نقشهی کاربری اراضی حوضهی آبریز رودخانهی گاوی تعیین درجهی رسوبدهی واحدهای هیدرولوِژیکی پس از وزندهی و نهاییسازی لایههای اطلاعاتی که به صورت لایههای وزنی درآمدهاند برای تهیهی نقشهی فرسایش حوضه (شکل 6)، لازم است این عوامل را با هم تلفیق کنیم، لازم به ذکر است لایههای وزنی باید به صورت رستری باشند تا امکان تلفیق آنها وجود داشته باشد. سپس به منظور تعیین میزان فرسایش هر یک از واحدهای هیدرولوژیکی، شناسایی تأثیرگذارترین عامل در تولید رسوب و نیز تشخیص واحدهای بحرانی از نظر تولید رسوب، محاسبهی درجهی رسوبدهی هر یک از این واحدها ضروری است. برای این منظور مجموع امتیاز عوامل نه گانه در هر واحد محاسبه و درجهی رسوبدهی آن واحد مشخص میگردد (جدول 4).
شکل(6) نقشهی شدت فرسایش حوضهی آبریز رودخانهی گاوی به منظور بررسی بیشتر و تعیین میزان فرسایش و رسوبزایی در سطح حوضه از جدول (5)، که در این زمینه تهیه شده است، استفاده میشود. در نهایت برای اعمال دقت بیشتر و پرهیز از اشتباه در عمل درونیابی و برونیابی برای برآورد فرسایش و تولید رسوب طبق مدل ام پسیاک از رابطه بین درجهی رسوبدهی و میزان تولید رسوب، رابطهی (10)، استفاده به عمل میآید (صفامنش و همکاران[7]، 2006): رابطهی (10) 18.6=Qs در این رابطه QSمیزان رسوبدهی سالانه بر حسب متر مکعب در کیلومتر مربع،e عدد نپرین که تقریباً برابر با 718/2 میباشد وR درجهی رسوبدهی یعنی مجموع امتیازات عوامل مختلف در نظر گرفته شده در مدل MPSIACمیباشد. لازم به ذکر میباشد که وزن مخصوص برای حوضهی مورد مطالعه با توجه به مطالعات انجام شده توسط امور آب استان بر روی رودخانه، 5/1 در نظر گرفته میشود. برای محاسبهیSDR از روش لارنس[8] رابطهی (11) استفاده شده است و برای تعیین فرسایش کل بایستی فرسایش ویژه را در کل مساحت هر زیر حوضه ضرب کنیم (افسری و قدوسی 1390). رابطهی (11) SDR = جدول (4) دادههای لازم جهت محاسبه درجهی رسوبدهی واحدهای هیدرولوژیکی حوضهی آبریز رودخانهی گاوی
جدول (5) تعیین میزان تولید رسوب سالانه و کلاس فرسایش خاک در روش MPSIAC
منبع: (رفاهی، 1385) بحث و نتایج این پژوهش با هدف برآورد فرسایش و رسوبزایی حوضهی آبریز رودخانهی گاوی به روش ام پسایک و مقایسهی نتایج با برآورد واقعی انجام شده است. نتایج بررسی روش ام پسیاک و 9 پارامتر این روش برای 4 زیرحوضه تعیین شده (جدول 6)، نشان میدهد که با توجه به مجموع امتیازات پارامترهای مؤثر در فرسایش و رسوب، حوضهی مورد مطالعه بـه صورت کلی با درجهی رسوبدهی 47/77 در کلاس فرسایشی زیـاد قرار میگیرد و میزان تولید رسـوب این کلاس 1429-476 متر مکعب در کیلومتر مربع و 5/2143-714 تن در کیلومتر مربع میباشد. با اعمال رابطهی (10)، میزان رسوبدهی سالانه برابر 13/302 متر مکعب در کیلومتر مربع برآورد گردیده که با دخالت دادن وزن مخصوص (5/1)، معادل 20/453 تن رسوب در کیلومتر مربع است. جدول (6) نتایج مدل MPSIACدر حوضه
نتیجهگیری نتایج حاصله نشان مـیدهد که منطقهی مورد مطالعه از لـحاظ فرسایشپذیری در ردهی فرسایشی زیـاد قـرار میگیرد و طبق تعریف این کلاس میتوان گفت که در این حوضه جابجایی ذرات خاک به میزانی است که اجرای برنامههای حفاظت خاک و آب ضرورت و اولویت داشته و استفاده از اراضی محدودیت زیادی دارد. همچنین مدل MPSIAC در این حوضه به خوبی پاسخ میدهد و فرضیهی اول تأیید میشود، زیرا به عنوان نمونه مقدار رسوب برآورد شده در حوضهی مورد مطالعه با مدل ام پسیاک 20/453 تن رسوب در کیلومتر مربع در سال و رسوب برآورد شده در ایستگاه رسوبسنجی حوضهی مربوطه 23/461 تن رسوب در کیلومتر مربع در سال میباشد، دلیل دقت این روش میتواند این باشد که این روش به عنوان جغرافیاییترین مدل در برگیرندهی 9 پارامتر تأثیرگذار در تولید فرسایش و رسوب میباشد. در مورد فرضیهی دوم قبل از انجام تحقیق با توجه به بازدیدهای به عمل آمده از منطقه مورد مطالعه به نظر میرسید که از بین عوامل مؤثر در فرسایش، لیتولوژی با توجه به سازندهای حساس به فرسایش، بیشترین تأثیر را در فرسایش حوضه داشته باشد، در حالی که نتایج حاصل شده از طریق مدل تجربی عامل کاربری با امتیاز 48/18را به عنوان عوامل مؤثر در فرسایش حوضه در اولویت اول و عامل زمینشناسی با امتیاز 76/5 در اولویت ششم قرار میدهد. بر این اساس فرضیهی فوق رد میشود. در پایان باید اشاره کرد که با توجه به اینکه حوضه از نظر تقسیمبندی واحد اراضی در گروه دشتهای دامنهای با شیب ملایم قرار میگیرد و مناسب کشت غلات به صورت دیم و مرتع آن محل قشلاق عشایر میباشد، در نتیجه این شرایط تأثیر عامل کاربری اراضی را تأیید میکند. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
مراجع | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ـ افسری، روحالله و جمال قدوسی (1390)، «ارزیابی روشهای مختلف تخمین نسبت تحویل رسوب (SDR) تحت شزایط آب و هوایی مختلف (مطالعهی موردی حوضههای آبریز استان مرکزی)»، فصلنامهی جغرافیای طبیعی سال چهارم شماره 12. ـ احمدی، حسن؛ ملکیان، آرش و رقیه عابدی (1391)، «مناسبترین روش آماری برآورد رسوب معلق رودخانهی جاجرود»، فصلنامهی پژوهشی پژوهشهای فرسایش محیطی، شمارهی 2. ـ ارخی، صالح؛ ملکی، سعید و رحیم نظری (1387)، «پهنهبندی شدت فرسایش و تولید رسوب بر اساس مدل MPSIAC با استفاده از GIS (مطالعهی موردی: حوضهی آبریز در استان ایلام)»،همایش ژئوماتیک. ـ ادارهی کل مـنابع طبیعی اسـتان ایلام، «آنالیز فیزیکی شیمیایی پـروفیل خاک حوضهی آبـریز مهران، گـزارش خاکشناسی حوضه». ـ بیاتی خطیبی، مریم (1383)، «فرسایش، فرایندهای فرسایشی و شکلهای ناشی از آنها»، مجلهی رشد آموزش جغرافیا شمارهی 69. ـ پیرمحمدی، زیبا؛ فیضنیا، سادات؛ یونسزاده جلیلی، سهیلا و جهانگیر فقهی (1389)، «مقایسهی مدلهای کمی در برآورد فرسایش و رسوب در زیرحوضهی سامان عرفی چم لرستان»، نخستین کنفرانس پژوهشهای کاربردی منابع آب ایران، دانشگاه صنعتی کرمانشاه. ـ خیام، مسعود؛ غنمی جابر، مصطفی و رسول صمدزاده (1392)، «مقایسه کارایی مدلهای MPSIACوEPM در برآورد فرسایش و رسوبزایی حوضهی آبریز سقز چی چای نمین»، فصلنامهی ژئومورفولوژی کاربردی ایران شمارهی 1. ـ فیضنیا، سادات؛ مجدآبادی فراهانی، فرهاد؛ محسنی ساروی، محسن و محمود عرب خدری (1381)، «طول دورههای آماری مناسب برای برآورد میانگین رسوب سالانه و روابط آن با مساحت، اقلیم، زمینشناسی و پوشش گیاهی حوضهی آبریز»، مجلهی کشاورزی و منابع طبیعی گرگان، سال نهم شمارهی 3. ـ رفاهی حسینقلی (1388)، «فرسایش آبی و کنترل آن»، تهران، انتشارات دانشگاه تهران. ـ شایان، سیاوش و غلامرضا زارع (1390)، «تبیین مفهوم فرسایش از دیدگاه ژئومورفولوژی و مقایسهی آن با دیدگاه منابع طبیعی»، مجلهی پژوهشهای فرسایش محیطی، دانشگاه هرمزگان، سال اول، شمارهی اول. ـ عاقلی کهنه شهری، لطفعلی و حسین صادقی (1384)، «برآورد آثار اقتصادی فرسایش خاک در ایران»، فصلنامهی پژوهشهای اقتصادی، شمارهی 15. ـ علیزاده گرجی، غلامرضا (1385)، «برآورد فرسایش و رسوب با استفاده از مدل MPSIACدر محیط GIS (مطالعه موردی حوضهی کلیجانر ستاق) پایاننامهی کارشناسی ارشد»، دانشکدهی علوم انسانی و اجتماعی، دانشگاه تبریز. ـ عسگری، شمسالله (1386)، «مدلسازی فرسایش خاک وتولید رسوب حوضهی آبخیز گل گل- استان ایلام»، محل انتشار چهارمین همایش ملی علوم و مهندسی آبخیزداری ایران. - مومیپور، مـهدی (1383)، «مقایسه دو مدل MPSIAC، RUSLEبا اسـتفاده از دادههای دورسـنجی و GISمطالعهی موردی زیرحوضهی اوجان چای»،پایاننامهی کارشناسی ارشد، دانشکدهی علوم انسانی علوم اجتماعی، دانشگاه تبریز. ـ مقصودی، مهران و مجتبی یمانی (1388)، «برآورد فرسایش ورسوب از طریق ارزیابی متغیرهای تأثیرگذار در حوضه ی آبریز وزنه با استفاده از GIS»، جغرافیا و توسعه، شمارهی 16 (زمستان). - Dearman,W.H. (1970), “Weathering Classification In The Chavacteriz- Ation of Erosion”Intel.Assoc.Engry Geol”, BuII., No. 13, PP: 123-128. - Hsu, S.M, Wen,H.Y.,Chen, N.C., Hsu, S.Y. and Chi, S. (2012), “Using an Intergrated Method to Estimate Watershed Sediment Yieldduring Heavy Rain Paiod: A Case Study in Hualin County”, Taiwan, Natural Hazards and Earth System Sciences, Vol.12, No. 6. - Li . M,Y.W. (2010),“Effects of Landforms on The Erosion Rate in a Small Watershed by The Cs Tracing Method”, Journal of Environmental Radioactivity, Vol. 101, No. 5. - Morgan, R.P.C, (1995), “Soil erosion and Conservation, Longman Scientific and Technica, John Wiley and Sons”. - Renard, K.G, J.J Ston. (1982), “Sediment Ield from Small Semiarid Rangland Watersheds, USDa-SEA_ARM, Westem Series”, No. 26, PP: 129-144. - Mayer, L.D et al. (1969), “Mathmation Simulation of Processes of Soil Erosion by Water”, Trans.Am..soc.agric.Engric, 12(6), PP: 754-758. - PSIAC, (1968), “Report of the Water Management Subcom-mittee on Factors Affecting Sediment Yild in the Pacific South-west Area and Selection and Evaluation of Measures for Reduc-tion of Erosion and Sediment Yield”, American Society of Civil Engineers,98, Report No. HY12. - Safamanesh, R., AzminSulaiman, W.N., Ramli, M.F., (2006),“Erosion risk Assessment Using an Empirical Model of Pacific South West Inter Agency Committee Method for ZargehWaterhed, Iran”, Journal of spatial Hydrology Vol. 6, No.2. - Victora, C., A. Kacevas,&h.Fiori.(1998), “Soil Credibility Assessments with Simulated Rainfall and with the USLE Monograph in Soil from Uruguay”, Proceeding of 16th World Congress of Soil Science, Montpllier, France. - Mansouri Daneshvar, M.R. and Bagherzadeh, A .(2012), “Evaluation of Sediment Yield in Psiac and Mpsiac Models by Using GIS at Torop Watershed”, Frontiers of Earth Science in China, 83. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
آمار تعداد مشاهده مقاله: 3,127 تعداد دریافت فایل اصل مقاله: 1,559 |