Document Type : مروری
Abstract
Among all natural disasters floods inflict the greatest damages to agriculture, fisheries, housing and infrastructure seriously affecting the social and economic activities of many people around the world. The many floods which have struck the city of Maku, especially northern mountainous sub basins, leaving behind financial losses and damages to human lives indicates a necessity of an integrate study about flooding in the upstream basins of the region. In this study, first, physiographic characteristics of the study area were determined by using satellite imagery, GIS and WMS analytic functions. Then, the flood zone area was identified, with emphasis on the city of Maku, and finally, the stream flooding was analyzed in the city skeleton. The results showed that rough topography, low thickness of the separate elements on bed rock and their low capacity for water storage, poor vegetation, violating the flood boundaries and changing them to artificial land, provides the conditions for the devastating floods especially during the showery precipitation in the dry seasons. As the findings show the runoff from a rainfall with the 120 mm/h intensity was estimated about 2.5 million m3, that 338000 m3 of it was obtained from the small watershed areas overlooking the city.
Keywords
مقدمه
سیل یکی از پیچیده‌ترین و مخرب‌ترین رویدادهای طبیعی است و بیش از هر بلای دیگری، جان و مال انسان و شرایط اجتماعی و اقتصادی جامعه را به مخاطره می‌اندازد، هر چند مدیریت صحیح آن میتواند منبع بزرگی جهت تأمین نیازهای آبی یک منطقه باشد (تلوری، 1376: 9). اگر وقوع این پدیده، ناگهانی و شدید باشد، باعث بهوجود آمدن خسارت‌های زیادی میشود. بررسی مجموعه عوامل محیطی که موجب ایجاد سیل هستند، نشان میدهد که دخالت انسان در چرخه طبیعی آب از طریق تخریب پوشش گیاهی در عرصه‌های آبخیز، کاربری غیراصولی اراضی، توسعه سطوح نفوذناپذیر و امثال آن امکان وقوع سیل را در مناطق گوناگون افزایش میدهد. جمع‌آوری و دفع روانآب‌های ناشی از بارندگی در مسیل رود در واقع از اقدامات ایمنی، بهداشتی و رفاهی به حساب میآیند (قنواتی و همکاران، 1391: 121).
بدیهی است که عوامل اقلیمی و عوامل حوضه‌ای در بروز سیلابها تأثیر دارند، اما برای جلوگیری از آثار زیانبار سیل نمی‌توان در عوامل و عناصر جوی تغییری ایجاد نمود. بنابراین هرگونه راه‌حل اصولی و چاره‌ساز را باید در روی زمین و اختصاصاً در حوضههای آبخیز جستجو کرد. در این ارتباط، اولین اقدامی که برای کاهش خطر سیل مطرح میشود، مهار سیل در سرچشمه آن یعنی زیرحوضههای آبخیز است. مسلماً برای انجام دادن اینکار شناسایی مناطق سیل خیز در داخل حوضه لازم است. تولید روانآب در یک حوضه آبخیز به عوامل متعددی بستگی دارد که از جمله آنها میتوان به خصوصیات حوضه آبخیز، میزان بارش، نفوذ و شرایط پیشین حوضه اشاره کرد (سینگ[1]، 1997: 18).
تجارب مطالعاتی و مدیریتی کشورهای مختلف نشان میدهد که اولین گام در جهت کاهش آثار زیانبار سیل، شناخت مناطق سیل‌گیر و پهنهبندی این مناطق از لحاظ میزان خطر سیل‌گیری است تا بتوان بر اساس نتایج بهدست آمده، درباره نحوه استفاده از اراضی و کاربری‌های مختلف کشاورزی، صنعتی، خدماتی و مکان‌یابی سکونتگاه‌های شهری و روستایی تصمیمگیری اصولی و بهینه نمود و آثار زیانبار سیل را تا حد ممکن به حداقل رسانید. یکی از ابزارهای بسیار مؤثر در مطالعات محیطی و علوم زمین، استفاده از فناوری سنجش از دور و بهرهگیری از دادههای ماهوارهای است (علوی پناه، 1382: 15).
یکی از مخاطرات طبیعی که شهر ماکو در طی سالهای اخیر با آن مواجه بوده، پدیده سیل است. توسعه فیزیکی نامناسب شهر بر روی مسیلهای قدیمی، بهعلت نبود فضاهای مناسب برای ساخت و سازها، هر مسیری را برای جریان روان آب سطحی شهر مسدود کرده است. این امر سبب میشود که هنگام وقوع بارندگی، روان آب حاصل از ارتفاعات اطراف، بهویژه از قسمت شمالی شهر، بهصورت سیل وارد بافت کالبدی شهر شده و خسارت‌های زیادی بهبار آورد.
شهر ماکو در دره‌ای با بستر سنگی سخت و نفوذناپذیر در داخل حوضه آبریز زنگمارچای استقرار یافته است. بیش از 578 کیلومتر مربع از مساحت این حوضه، در بالادست شهر ماکو جا گرفته است و بهدلیل خصوصیات فیزیوگرافی خاص حوضه آبریز و وجود دامنه‌های پرشیب کوهستانی مسلط به شهر در ضلع شمالی و جنوبی تحت تأثیر سیلهای مخرب بهاری و تابستانه قرار دارد. بهعنوان مثال در 17 مرداد ماه 1387، در نتیجه یک بارش رگباری با شدت 20 میلیمتر در 10 دقیقه، در شهر ماکو سیلی جاری شد که باعث ایجاد خسارت بر 350 واحد مسکونی، 700 دستگاه خودرو، 58 واحد صنفی، آب گرفتگی واحدهای مسکونی 150 خانوار و تخریب کلی 22 واحد مسکونی شد. در این جریان 82 نفر مصدوم شدند و 4 نفر نیز جان خود را از دست دادند (هفته نامه آوای ماکو 31 مرداد 1387، شماره 209). مسؤولان اجرایی شهرستان ماکو میزان خسارت ناشی از این سیل را 82 میلیارد و 853 میلیون ریال اعلام کردند. بارش‌های رگباری سالهای 1387 و 1390 و تابستان 1391 که منجر به آبگرفتگی شدید معابر عمومی در ماکو شد، بیانگر شایع بودن وقوع سیل در این شهر و وجود ابهامات گسترده در مسایل مرتبط با مدیریت کمی و کیفی سیلابهای شهری است و برای چاره‌اندیشی آنها و مدیریت یکپارچه شهری و نیز توسعه پایدار، اجرای طرح‌های تحقیقاتی اجتناب ناپذیر است. لذا پژوهش حاضر با هدف شناسایی مناطق مستعد تولید روانآب، پهنهبندی مکان‌های در معرض خطر سیل و یافتن راهکارهای مناسب جهت کاهش آثار زیانبار سیل در شهر ماکو به انجام رسید.
در ارتباط با شناسایی و پهنهبندی مناطق سیل‌گیر با استفاده از روشهای سنجش از دور و GIS و تجزیه و تحلیل میزان خطر، پژوهش‌های زیادی در داخل و خارج انجام گرفته است. در سال 1994 سوان وراکامتون[2] با استفاده از مدل HEC-1 و سیستم اطلاعات جغرافیایی (GIS) آثار تغییر کاربری اراضی بالادست حوضه آبخیز را روی الگوی سیلاب در نواحی پایین دست حوضه مورد ارزیابی قرار داد. نتایج تحقیق نشان داد که با کاهش مساحت جنگل، روانآب حوضه اصلی و زیرحوضه‌ها افزایش یافته و تغییرات کاربری اراضی در بالادست حوضه در بهوجود آمدن سیل در اراضی پایین دست حوضه تأثیر اساسی داشـته است. فرانسیسکو[3] و همکاران (1999) بهمنظور اتصال[4] سیستم GIS با مدلهای هیدرولوژی و هیدرولیکی، مناطق تحت تأثیر سیل را در رودخانهها به نقشه در آوردند. استفان[5] (2002) در تحقیقی در حوضه رودخانه سنگزرد[6] واقع در ایالت مونتانای امریکا به بررسی سیلهای به وقوع پیوسته در سالهای 1996 و 1997 پرداخت و نهایتاً پهنههای سیل 100 و 500 ساله را تعیین کرد و نتیجه گرفت که تغییر کاربری اراضی بیشترین تأثیر را در ایجاد سیلاب دارند. کلدیتز[7] و همکاران (2003) در حوضه پانکو در مکزیک اقدام به تهیه نقشه پهنهبندی سیل را با تلفیق اطلاعات سنجش از دور با روش ژئومورفولوژیکی تهیه کرد. طبق نظر ایشان تلفیق اطلاعات سنجش از دور میتواند نقشههایی با توان تفکیک بالا تهیه کند که برای پهنهبندی سیل در مناطق وسیع میتواند بسیار مفید باشد. فیلیپ[8] و همکارانش (2006: 1180) با مطالعه رکوردهای تاریخی درازمدت از بارندگی، روان آب و دیگر عوامل اقلیمی را برای تحقیق در تمایلات و تغییرات هیدرولوژیکی در حوضه آبریز رودخانه وُلتا[9] (واقع در غرب آفریقا) در دوره زمانی 1901 تا 2002 بررسی کردهاند. یافتههای تحقیق آنها نشان داد که عوامل اقلیمی (شدت بارندگی، نوع بارش و دمای هوا) تأثیر مهمی در چگونگی سیل (شدت و دوام) در ارتباط با ویژگیهای هیدرولوژیکی حوضه (شکل حوضه، جنس زمین، شیب و کاربری اراضی) دارد.
در ایران نیز وهابی (1376) با بهکارگیری روشهای سنجش از دور و GIS، برای پهنهبندی خطر سیل در حوضه آبریز طالقان اقدام کرد. طیبی (1384) با استفاده از مدلهای منطق بولین و فازی در GIS، به مطالعه و پهنهبندی توان سیلخیزی با توجه به ویژگیهای ژئومورفیک حوضه آبریز سبزوار پرداخته است. معیری و انتظاری (1388: 1) علل و عوامل مؤثر بر ایجاد سیلابها و میزان خسارت ناشی از آنها را در شهرستانهای پرخطر استان اصفهان مطالعه کردهاند. کردوانی و همکاران (1388: 1) در تحلیل منطقهای سیلابهای لارستان به این نتیجه رسیدند که مهمترین عوامل وقوع سیلابهای لار شدت و مدت بارش و همچنین توپوگرافی منطقه است و اصلیترین مسیل لار که به بافت قدیم شهر آسیب میرساند، مسیل دربند است مسیل تنگ اسد نیز به بافت جدید شهر آسیب میزند. قهرودی تالی (1388: 167) مدل یکپارچه سیلاب شهری را در بخشهایی از شمال و شمالشرق کلانشهر تهران مطالعه کرده و نشان داد که برای کاهش سیلاب در تهران، میتوان سیلابهای بالادست را کنترل نمود. در این راستا چهار نقطه را برای مدیریت سیلابها مکانگزینی کردهاند. همچنین از نظر ایشان به این دلیل که شبکه ارتباط شهری بدون توجه به جهت زهکشی طبیعی تهران ایجاد شدهاند؛ مسیلهای تهران را از هدایت سیلابها ناتوان ساخته است. لذا این وظیفه را گذرگاههای شهری مانند اتوبانها، خیابانها و کوچهها انجام میدهند. نصری و سبزی (1390) سـیلابها و شبکه مسیلهای تأثیرگذار به شـهر زواره (شهرستان اردستان استان اصفهان) را با توجه به برنامهریزی شهری مورد مطالعه قرار داده و نشان دادند که با شناخت وضعیت شبکه مسیلها، دبی جریان و ویژگیهای مورفومتری حوضه آبریز مشرف به شهر، میتوان برای ایمنسازی شهر را برنامهریزی کرد. در جدیدترین مطالعات نیز بومری و همکاران (1390: 129) پهنههای سیلابی و ویژگیهای فیزیوگرافی و کمی حوضه آبریز دامن (در استان سیستان و بلوچستان) را با استفاده از سیستم اطلاعات جغرافیایی و سنجش از دور مطالعه کرده و نتیجه گرفتند که بهرهگیری از تصاویر ماهوارهای برای مناطقی که فاقد آمار و اطلاعات کافی در خصوص پدیده سیل هستند، بسیار مفید است و با تفسیر تصاویر با وضوح بالا میتوان پهنههای آسیبپذیر از سیل را شناسایی و طبقهبندی کرد.
با عنایت به اهمیت موضوع مدیریت سیلهای سالیانه در شهر ماکو و نبود مطالعات در این راستا، پژوهش حاضر ضمن استفاده از تجربیات مطالعات پیشین در مورد بهکارگیری سنجش از دور و GIS، بهخصوص در مناطقی که کمبود دادههای آماری وجود دارد، به شناسایی مناطق پر خطر و مدیریت بهینه سیلابها میپردازد.
شناخت عرصه پژوهش
ماکو از شهرهای شمالغربی استان آذربایجان غربی است و از نظر موقعیت مطلق جغرافیایی مابین 39 درجه و 16 دقیقه و 37 ثانیه تا 39 درجه و18 دقیقه و 3 ثانیه عرض شمالی و 44 درجه و26 دقیقه و 28 ثانیه تا 44 درجه و33 دقیقه و 29 ثانیه طول شرقی قرار گرفته است. این شهر بر اساس سرشماری سال 1390، دارای 52302 نفر جمعیت و 22842 خانوار بوده است. از نظر موقعیت استقرار، شهر ماکو در تنگه‌ای به طول 14 کیلومتر در بین دو رشته کوه مرتفع «قَیَه داغ» در شمال و «سبد داغی» در جنوب واقع شده است. این دو رشته کوه به موازات هم و با فاصلهای در حدود 500 متر از هم قرار گرفتهاند. رودخانه زنگمار تنها مجرایی است که در داخل بافت کالبدی شهر از غرب به شرق جریان مییابد. طول این رودخانه در داخل کالبد فیزیکی شهر در حدود 5/12 کیلومتر بوده ولی عرض آن متفاوت است و کل آب‌های سطحی شهر و حوضه 57845 هکتاری بالادست خود را زهکشی میکند. بستر و کناره‌های این رودخانه از سنگهای بازالتی حاصل از آتشفشان کوه آرارات شکل یافته است و پیچوخم‌های زیادی دارد. با عنایت به اینکه بیشتر روان آب‌های سطح حوضه در بدو ورود به شهر، وارد مجرای رودخانه زنگمار شده و از آن خارج میشود ولی اغلب سیلهای مخرب شهر ماکو از دامنه‌های جنوبی «قیه داغ» و مسیل «دانالو» وارد شهر میشود و خسارت‌های سنگینی را به بافت شهر تحمیل می‌نماید، این قسمت از حوضه مورد مطالعه و بررسی قرار گرفت (اشکال 1و 8) و (جدول 1).
شکل (1) موقعیت منطقه مورد مطالعه و طبقات ارتفاعی آن
مواد و روشها
برای انجام دادن این تحقیق با توجه به اهداف و ماهیت آن، از روشهای کتابخانهای و بازدیدهای میدانی استفاده کردیم. بدین ترتیب که در مراحل مختلف تحقیق، بازدیدها و برداشتهایی از منطقه مورد مطالعه به عمل آمده است. بهطور کلی برای انجام دادن پژوهش حاضر از دادهها و اطلاعات زیر استفاده شده است:
ـ تصاویر ماهواره SPOT سال 2005؛
ـ نقشههای توپوگرافی منطقه ماکو با مقیاس 50000 :1؛
ـ نقشه زمینشناسی (برگ ماکو) با مقیاس 250000 :1؛
ـ نقشه ارزیابی منابع و قابلیت اراضی استان آذربایجان غربی در مقیاس 250000 :1؛
ـ دادههای ایستگاه‌های هواشناسی، بارانسنجی و دبیسنجی منطقه ماکو؛
ـ یک دستگاه سامانه موقعیتیاب جهانی[10]
تصاویر و نقشههای تهیه شده در بدو امر آمادهسازی و زمین مرجع شدند. تجزیه و تحلیل اطلاعات تهیه شده را با شناسایی و جداسازی محدودهای از حوضه آبریز رودخانه زنگمارچای که در بالادست ساختگاه شهر ماکو واقع شده، آغاز کردیم. در این راستا از قابلیتهای نرمافزارهای Arc/map، spatial analyst plus، Arc/hydro و سامانه مدیریت حوضههای آبریز (WMS)[11] استفاده کردیم. به این منظور ابتدا DEM حوضه را ایجاد و محدوده مورد نظر را در شش واحد هیدرولوژیکی مجزا، تجزیه و تحلیل و ارزیابی کردیم.
یکی از معیارهای اساسی در تولید روان آب حوضه، زمان تمرکز است. که بر اساس فرمولهای تجربی مختلف قابل محاسبه است. در محدوده مورد مطالعه برای محاسبه آن از معادله کرپیچ استفاده شده است:
رابطه (1)
در معادله فوق
tc = زمان تمرکز بر حسب ساعت
L = طول مسیر حرکت آب در داخل حوضه (کیلومتر)
H = اختلاف ارتفاع بین نقطه تمرکز و بلندترین قسمت حوضه (متر)
میزان روان آب در یک حوضه را میتوان متناسب با خصوصیات فیزیکی آن و شدت بارندگی، بر اساس فرمولهای تجربی بر حسب ارتفاع یا حجم روان آب محاسبه کرد. یکی از فرمولهای کاربردی روش SCS است.
رابطه (2)
R = ارتفاع روان آب بر حسب اینچ
P = ارتفاع بارندگی (اینچ)
S = عامل مربوط به نگهداشت آب در سطح زمین است که مقدار آن برابر است با:
رابطه (3)
CS = مقدار نفوذ آب در حوضه است که بر اساس نوع خاک و نوع اراضی از روی جداول تعیین میشود (علیزاده، 1377: 370).
برای بررسی وضعیت پوشش گیاهی منطقه، با بهکارگیری رابطه زیر، شاخص NDVI از تصویر ماهوارهای SPOT استخراج شد:
رابطه (4)
بر اساس این شاخص، تراکم پوشش گیاهی بین 1 و 1- تعیین میشود. هرچه عدد به سمت 1 سوق پیدا کند، نشانگر تراکم زیاد پوشش گیاهی و هر چه به سمت 1- پیش رود، بیانگر کمبود و ضعف پوشش گیاهی است.
برای استخراج نقشه کاربری اراضی منطقه مورد مطالعه از رهیافت طبقهبندی نظارتشده[12] به روش حداکثر احتمال[13] تصاویر ماهوارهای استفاده شد. این روش واریانس و کوواریانس کلاس‌ها را ارزیابی میکند (علوی پناه، 1382 : 311).
خصوصیات فیزیوگرافی منطقه مورد مطالعه
جهت تحلیل ویژگیهای فیزیوگرافی حوضه و تعیین مرز زیرحوضه‌ها، و نیز برای درک بهتر رفتار هیدرولوژیکی آنها در برابر رگبارها، با استفاده از قابلیتهای نرمافزاری Arc/map و spatial analyst plus و Arc/hydro و WMS، ابتدا DEM حوضه ایجاد شد. سپس محدوده مورد نظر در شش واحد هیدرولوژیکی مجزا، مورد تجزیه و تحلیل قرار گرفت (جدول4). نتایج نشان میدهد که بیش از 11 درصد حوضه شیب 25 درجه و بیشتر از آن را دارد و اکثر مناطق پرتگاهی در اطراف شهر ماکو قرار گرفتهاند (جدول 1 و 2 و اشکال2 و 3).
جدول (1) دامنه تغییرات شیب در زیرحوضههای مشرف به شهر ماکو
|
دامنه شیب(درجه) |
5 - 0 |
14 – 6 |
24 – 15 |
37 – 25 |
70 - 38 |
کل |
|
سطح اشغال شده(هکتار) |
6/689 |
8/881 |
751 |
9/367 |
9/125 |
2/2816 |
|
درصد مساحت |
5/24 |
3/31 |
7/26 |
13 |
5/4 |
100 |
شکل (2) نقشه واحدهای هیدرولوژیکی منطقه مورد مطالعه
شکل (3) نقشه شیب حوضه زنگمار در بالا دست شهر ماکو
جدول (2) دامنه تغییرات شیب و میزان سطح اشغال آن در حوضه زنگمار (دست شهر ماکو)
|
درصد مساحت |
سطح اشغال شده (هکتار) |
دامنه شیب (درجه) |
|
22/31 |
18060 |
5 - 0 |
|
12/36 |
20890 |
14 – 6 |
|
56/21 |
12473 |
24 – 15 |
|
61/8 |
4980 |
37 – 25 |
|
5/2 |
1440 |
70 - 38 |
|
100 |
57843 |
کل |
وضعیت پوشش گیاهی
کسب اطلاعات درباره وضعیت پوشش گیاهی، از قبیل میزان و پراکنش آنها، بهلحاظ تأثیری که در مورفوژنز و مورفودینامیک منطقه می‌گذارد، اهمیت زیادی دارد. لذا مطالعه آن میتواند اطلاعات با ارزشی از فرایندهای فرسایشی و تحول منطقه ارائه دهد. تراکم پوشش گیاهی در منطقه یکسان نیست بلکه با توجه به عوامل متعدد مانند ارتفاع، دما، بارش، نوع خاک، توپوگرافی و اقلیم از پراکندگی خاصی دارد. گردآوری اطلاعات در مورد تغییرات پوشش گیاهی با روشهای معمولی بسیار مشکل و پرهزینه است. لذا استفاده از دادههای دورسنجی امکان مطالعه پوشش گیاهی را فراهم می‌آورد (شکل 4).
نوع اراضی
مطالعات میدانی و بررسی نقشههای موجود نشان داد که نوع غالب اراضی در این حوضه، کوهستانی است و در پای دامنه‌های آن ریزش واریزه نیز دیده میشود، اما اراضی پست و هموار دشت‌های سیلابی وسعت کمتری دارند (شکل5).
وضعیت خاک حوضه
خاک عمده حوضه مورد مطالعه با توجه به نوع و کیفیت خاک‌ها در سیستم فائو، بیشتر، از نوع خاکهای رگوسول[14] و لیتوسول[15] است و در اغلب نقاط با محدودیت عمق خاک، شوری و فرسایشپذیری همراه است. نقشه خاک منطقه بر اساس طبقهبندی فائو در شکل(6) نشان داده شده است.
|
|
شکل(4) (الف) تصویر ماهوارهای SPOT-5 و (ب) NDVIمنطقه مورد مطالعه
کاربری اراضی
روشهای مختلفی برای طبقهبندی تصاویر ماهوارهای بهکار می‌روند که از میان آنها طبقهبندی نظارت شده[16] بهروش حداکثر احتمال[17] تاکنون بهعنوان دقیقترین و پراستفادهترین روشها ذکر شده است. این روش واریانس و کوواریانس کلاس‌ها را ارزیابی میکند (علوی پناه، 1382: 311). نقشه کاربری اراضی حوضه مورد مطالعه با استفاده از تصویر ماهوارهای Spot-5 و پس از محاسبات آماری (جدول3 و 4) جهت ایجاد FCC ، در محیط نرمافزاری ENVI با روش حداکثر احتمال استخراج شد. بر اساس این تحلیلها، تقریباً 40 درصد منطقه را انواع مرتع متراکم و نیمهمتراکم تشکیل میدهد، 24 درصد اراضی بدون کاربری است و 17 درصد آن را برونزد سنگی پوشش میدهد (شکل 7).
جدول (3) پارامترهای آماری باندهای طیفی
|
Basic Stats |
Min |
Max |
Mean |
Stdev |
|
Band 1 |
0 |
254 |
01/50 |
87/51 |
|
Band 2 |
0 |
254 |
94/34 |
24/39 |
|
Band 3 |
0 |
254 |
19/39 |
38/42 |
|
Band 4 |
0 |
254 |
37/49 |
93/51 |
جدول (4) ضریب همبستگی بین زوج باندها
|
Correlation |
Band 1 |
Band 2 |
Band 3 |
Band 4 |
|
Band 1 |
1 |
88/0 |
92/0 |
93/0 |
|
Band 2 |
88/0 |
1 |
99/0 |
97/0 |
|
Band 3 |
92/0 |
99/0 |
1 |
97/0 |
|
Band 4 |
93/0 |
97/0 |
97/0 |
1 |
شکل (5) نقشه نوع اراضی منطقه مورد مطالعه (منبع: نقشه زمین شناسی ماکو و نقشه قابلیت اراضی آذربایجان غربی)
شکل (6) نقشه نوع خاک منطقه مورد مطالعه (منبع: نقشه زمینشناسی ماکو و نقشه قابلیت اراضی آذربایجان غربی)
شکل (7) نقشه کاربری اراضی منطقه مورد مطالعه مستخرج از تصویر ماهوارهای SPOT – 2005
یافتههای تحقیق
ارزیابی میزان روانآب سطوح طبیعی بهعلت تأثیرات متقابل بین توپوگرافی، عمق آب زیرزمینی، شدت و تداوم بارش و خصوصیات نفوذپذیری که بهوسیله خاک و پوشش گیاهی تعدیل میشوند، پیچیده‌تر است زیرا تمام سطح حوضه در تولید روانآب نقش یکسانی ندارند. به این ترتیب، نقش شرایط اولیه حوضه در پتانسیل تولید روانآب اهمیت پیدا میکند (Chow and et al; 1988: 532). در خصوص مساحت حوضه و ارتباط آن با مقدار روانآب تولید شده، برخی از محققان، وسعت مساحت را به افزایش حجم روانآب (صدقی؛ 1363: 14-2) و بعضی دیگر افزایش مساحت را به فزونی دبی منوط دانستهاند (Singh, 1997: 18). در منطقه مورد مطالعه وجود گستره وسیعی از سطوح نفوذناپذیر و شیب تند آنها شرایط ایجاد روان آب، و سیلهای مخرب به هنگام رگبارهای زودگذر را فراهم می‌آورد.
جریان سطحی هورتونی بهعنوان یکی از انواع اصلی تولید روانآب و سطوح مشارکتی، زمانی بهوجود می‌آید که شدت بارش بیش از نرخ نفوذ باشد. در این حالت، پس از تأمین ذخایر گودالی جریان سطحی بهوجود می‌آید. معمولاً جریان سطحی هورتونی از سطوحی با خاکهای رسی و کم عمق، سطوح عاری از پوشش گیاهی، سطوح شهری و سطوح با خاکهای فشرده که ظرفیت نفوذ را کاهش میدهند (مانند جادهها و امثال آن) تولید میشود (Stephen, 1999: 1-18). بررسی‌های بهعمل آمده نشان داد که جنس تولید روان آب در شهر ماکو، هورتونی است. قسمت‌هایی از شهر که مابین کوه «قیه داغ» و رودخانه زنگمار و پایین‌تر از «قره کورپی» واقع شده است، به تناسب شرایط نفوذ، شیب توپوگرافیک و نیز وسعت حوضه بالادست، بیش از سایر قسمت‌ها در معرض سیلهای مخرب قرار دارند. اما قسمت‌های دیگر شهر به دلیل ورود جریانها به داخل شاخههای رودخانه زنگمار، چندان آسیبپذیر نیستند (شکل8).
شکل (8) تصاویر google earth از محدوده مورد مطالعه
نمودارهای میانگین بارشها و میانگین حداکثرها، نشان میدهند که اولاً در چند سال اخیر نوسانات بارش سالانه غیر خطی است، ثانیاً بیشترین میزان بارشها در بهار و حداقل آن در تابستان اتفاق می‌افتد و در این فصل پوشش گیاهی آسیبپذیر میشود و نقش حفاظتی خود را از دست میدهد. لذا وقوع یک بارش رگباری و مازاد بر ظرفیت نفوذ، شرایط وقوع سیلابهای مخرب را تسهیل می‌نماید (اشکال 10-9). حجم بارش و روان آب حاصل از این بارش رگباری در هر یک از واحدهای هیدرولوژی در جدول 5 آمده است.
شکل (9) نمودار میانگین بارش سالانه، ایستگاه ماکو87-64
شکل (10) نمودار میانگین بارش‌های حداکثر ماهانه، ایستگاه ماکو87-64
خصوصیات نفوذپذیری حوضه که با نام CN یا شماره منحنی معرفی میشوند (علیزاده، 1377: 363)، ارتفاع روان آب و حجم روانآب در قسمت‌های مختلف حوضه، در جدول 5 و 6 آمده است. برآوردها، حجم روان آب حاصل از یک بارش رگباری 20 میلیمتری در 10 دقیقه را، در حدود 5/2 میلیون مترمکعب نشان میدهد و از این میزان، 338 هزار متر مکعب مستقیم وارد محدوده شهر ماکو میشود و آن را آسیبپذیر می‌سازد. نکته قابل توجه اینکه از حوضه بالادست صخره سنگی معروف شهر در دامنه «قیه داغی» که مساحتی در حدود 200 هکتار دارد، 28 هزار متر مکعب آب وارد شهر میشود (شکل 12).
بررسی‌ها نشان دادند که درحوضه مورد مطالعه 5/21 درصد بارشها تبدیل به روانآب میشود. شکل (11) پهنهبندی مناطق تولید روانآب و در معرض سیل را نشان میدهد.
جدول (5) خصوصیات فیزیکی و هیدرولوژیکی زیرحوضههای رود زنگمار در بالا دست ماکو
|
حجم روانآب (m3) |
ارتفاع روانآب (mm) |
ارتفاع بارندگی (Inch) |
CN متوسط |
حجم بارش (m3) |
مساحت (هکتار) |
زمان تمرکز (Tc) ساعت |
طول آبراهه (L) کیلومتر |
اختلاف ارتفاع (H)متر |
نام زیرحوضه مورد مطالعه |
|
808650 |
2/3 |
78/0 |
80 |
50540600 |
3/25270 |
94/0 |
10 |
1020 |
SB1 |
|
363000 |
7/3 |
78/0 |
85 |
15276000 |
7638 |
76/0 |
9 |
1280 |
SB2 |
|
211514 |
7/3 |
78/0 |
81 |
11433200 |
6/5716 |
7/0 |
7/8 |
1480 |
SB3 |
|
134241 |
3 |
78/0 |
77 |
8949400 |
7/4474 |
75/0 |
5/9 |
1580 |
SB4 |
|
329038 |
4/3 |
78/0 |
78 |
19355200 |
6/9677 |
95/0 |
12 |
1700 |
SB5 |
|
152031 |
3 |
78/0 |
80 |
10135400 |
7/5067 |
5/0 |
7/5 |
960 |
SB6 |
|
2493344 |
3/3 |
78/0 |
80 |
11568600 |
57845 |
5/5 |
55 |
1780 |
کل حوضه* |
* محاسبات به صورت وزنی و توسط نرم افزار صورت گرفته است.
شکل (11) نقشه پهنههای تولید روانآب در محدوده مورد مطالعه
دامنه‌های جنوبی «قیه داغی» که با شیب بسیار تندی به ساختگاه شهر ماکو مشرف است، در مجموع 2/2816 هکتار مساحت دارد. این بخش به چندین زیرحوضه کوچک تقسیمبندی شده است و مقدار CN آن، بهعلت پوشش سنگی و سنگلاخی بودن زمین در اکثر سطوح، رقم 90 تعیین گردید. لذا بر اساس روابط فوق، پارامترهای مربوط محاسبه، و در جدول (6) آورده شده است.
میزان تبدیل بارش به روانآب، برای زیرحوضه مشرف به شهر 60 درصد برآورد شده است.
جدول (6) خصوصیات فیزیکی و هیدرولوژیکی زیرحوضههای مشرف به شهر ماکو
|
حجم روان آب (m3) |
ارتفاع روان آب (mm) |
ارتفاع بارندگی (Inch) |
حجم بارش (m3) |
مساحت (هکتار) |
زمان تمرکز (Tc) ساعت |
طول آبراهه (L)کیلومتر |
اختلاف ارتفاع (H) متر |
زیر حوضه مورد مطالعه |
|
118440 |
12 |
78/0 |
197400 |
987 |
5/0 |
2/5 |
680 |
SB1 |
|
116076 |
12 |
78/0 |
193460 |
3/967 |
48/0 |
05/5 |
740 |
SB2 |
|
26760 |
12 |
78/0 |
44600 |
223 |
31/0 |
2/3 |
580 |
SB3 |
|
24600 |
12 |
78/0 |
41000 |
205 |
16/0 |
95/1 |
740 |
بشگوز |
|
27960 |
12 |
78/0 |
46600 |
233 |
13/0 |
45/1 |
500 |
کوه قیه |
|
337944 |
12 |
78/0 |
563240 |
2/2816 |
48/0 |
2/5 |
800 |
کل* |
* محاسبات به صورت وزنی و توسط نرمافزار صورت گرفته است.
شکل (12) نقشه پهنههای در معرض خطر سیل و آبگرفتگی در شهر ماکو
نواحی آسیبپذیر شهر ماکو در برابر سیل بهشرح زیر مشخص شد و باید عملیات حفاظتی در محلات زیر صورت گیرد (شکل 12):
قسمت مرکزی و بافت قدیمی (هسته اولیه) شهر، در دامنه‌های کوه «قیهداغ»، که جاری شدن سیل به صورت آبشار از بالای کوه به داخل شهر، بافت مسکونی و بافت تجاری اصلی شهر را تهدید میکند (شکل 13).
بافت شهر در دامنه‌های شرقی کوه «قیه داغ» در محله «بشگوز» در روی مسیلهای قدیمی و قسمت شمالی معبر اصلی شهر قرار گرفتهاند. معابر ارتباطی این واحدها، همان مسیلهای قدیمی است که بدون توجه به خطرات عبور سیل، مورد ساخت و ساز واقع شده و تبدیل به معبر گشته‌اند (شکل 14).
قسمت قدیمی شهر در سمت غربی بافت اولیه شهر و مقابل ساختمان شهرداری ماکو موسوم به «محله» که در دل کوه واقع شده است (شکل 15).
بافت تقریباً جدید شهر در سمت غربی کوه «قیه داغ» موسوم به «محله مسکن» که دقیقاً در مسیل کناری راه ارتباطی روستاهای سمت شمالی شهر (روستای دانالو) قرار دارد. این قسمت از بافت شهر در روی مسیل منتهی به رودخانه زنگمار قرار گرفته و آن را مسدود نموده است و تنها محل عبور برای سیلابها، معبر ارتباطی متصل به بلوار اصلی شهر است. سیل جاری شده با ورود به بلوار اصلی شهر بهویژه قسمت شمالی آن، در جهت شرق به طرف بخش مرکزی شهر جریان مییابد. در طول این معبر که تنها معبر اصلی و ارتباطی عمده شهر است، تمامی واحدها، بهویژه واحدهای سمت جنوبی (در جهت شیب منتهی به رودخانه زنگمار) در معرض سیل واقع میشوند.
|
|
شکل (13) جاری شدن سیل از بالای کوه «قیهداغ» به داخل شهر (منبع: آرشیو)
شکل (14) نفوذ سیل به داخل بخش مرکزی و تجاری(منبع: آرشیو)
شکل (15) نفوذ سیل و تخریب پارک در محله «بشگوز» و «محله» (منبع: آرشیو)
نتیجهگیری
پراکنش نامناسب زمانی و مکانی بارش در مناطق خشک و نیمهخشک، علاوه بر ایجاد سیلابهای مخرب، موجب هدر رفتن روانآب سطحی میشود. از این رو مهار روانآب سطحی و بهرهبرداری مناسب از آن میتواند راهکار مناسبی برای جلوگیری از هدر رفتن آب، و تبدیل تهدید‌ها به فرصت‌ها باشد. تعیین عرصه‌های مستعد تولید روانآب و سیل و مناطق در معرض آبگرفتگی یکی از مهمترین عوامل تعیینکننده در موفقیت طرح‌های سیلخیزی است. از سوی دیگر، تعیین مکان‌های مناسب با استفاده از روشهای سنتی و متداول بسیار دشوار بوده، اکثراً باعث بروز خطا میشود. سیستم اطلاعات جغرافیایی به همراه ابزارها و تکنیک‌های آن توانایی آن را دارد تا با تلفیق لایههای مختلف اطلاعاتی در قالب مدلهای مختلف و در حداقل زمان ممکن در مکان‌یابی این عرصهها، مورد استفاده قرار گیرد.
در این مطالعه، علاوه بر پهنهبندی مناطق سیلخیز حوضه، نقاط سیل گیر در داخل شهر ماکو نیز مورد ارزیابی قرار گرفت و مشخص شد که در محدوده مورد مطالعه، بخش کوچکی از حوضه که در شمال شهر با شیب تند و دارای مواد منفصل ناچیز وجود دارد، بیشترین نقش را در تولید سیلهای مخرب و تهدیدکننده شهر ماکو، ایفا میکند که میتوان با احداث کانال‌های انحرافی، جریان آب را در فواصل مختلف در طول شهر، به رودخانه زنگمار هدایت و با یک سد مخزنی در پایین دست، جهت استفاده در مصارف کشاورزی ذخیره نمود. از طرف دیگر نبود زمین لازم برای گسترش شهر و مشکلات مربوط به توسعه فیزیکی آن سبب شده است که تـمامی عرصه‌های موجود به اشغال کاربـری‌های مختلف شهری درآید. وجود دامنه‌های تند و تیز، برونزدهای سنگی و زمین‌های فاقد پوشش سطحی، محدودیت فضا و بارش‌های تند رگباری در فصل خشک، بروز سیلهای مخرب را در این ناحیه مساعدتر می‌سازد و عدم رعایت حریمها و تجاوز به بستر مسیلها و آب گذرها، آن را شدت می‌بخشد و خسارت‌های زیادی را در محله‌های «بشگوز»، «محله»، «مسکن» و بخش قدیمی شهر (هسته اولیه) واقع در پایین دست کوه «قیه داغ» به بار می‌آورد.
از آنجا که نفوذپذیری از عوامل مؤثر در ایجاد روان آب و سیل است و نخستین اثر توسعه کالبد فیزیکی شهر، کاهش نفوذپذیری است؛ باید از هرگونه عملیات توسعه شهری و ساخت و سازها در حریمها جلوگیری شود و با آزادسازی و تبدیل فضاهای اشغالی به فضاهای سبز و مهیا نمودن شرایط نفوذ، جریان روان آب و سیلهای ویرانگر را کاهش داد. همچنین میتوان اندازه مسیر عبور آب‌های سطحی (کانال) را با وسعت حوضه بالادستی آن متناسب ساخت و با احداث کانال‌های انحرافی بزرگتر، جریان سیل را به سمت مجاری اصلی هدایت نمود. لذا یکی از راهکارهای عملی کوتاهمدت، در جهت کاهش آثار زیانبار سیل در شهر آزادسازی حریمها، و در دراز مدت، ایجاد شهرک‌های جدید و انتقال جمعیت، پیشنهاد میشود.
[1]- Singh
[2]- Swan werakamton
[3]- Francisco
[4]- Coupling
[5]- Stephan
[6]- Yellowstone
[7]- Colditz
[8]- philip
[9]- Volta River Basin
[10]- Global Positioning System (GPS)
[11]- Watershed Management System
[12]- Supervised classification
[13]- Maximum likelihood
[14]- Reglesols
[15]- Lithosols
[16]- Supervised classification
[17]- Maximum likelihood
مقدمه
سیل یکی از پیچیده‌ترین و مخرب‌ترین رویدادهای طبیعی است و بیش از هر بلای دیگری، جان و مال انسان و شرایط اجتماعی و اقتصادی جامعه را به مخاطره می‌اندازد، هر چند مدیریت صحیح آن میتواند منبع بزرگی جهت تأمین نیازهای آبی یک منطقه باشد (تلوری، 1376: 9). اگر وقوع این پدیده، ناگهانی و شدید باشد، باعث بهوجود آمدن خسارت‌های زیادی میشود. بررسی مجموعه عوامل محیطی که موجب ایجاد سیل هستند، نشان میدهد که دخالت انسان در چرخه طبیعی آب از طریق تخریب پوشش گیاهی در عرصه‌های آبخیز، کاربری غیراصولی اراضی، توسعه سطوح نفوذناپذیر و امثال آن امکان وقوع سیل را در مناطق گوناگون افزایش میدهد. جمع‌آوری و دفع روانآب‌های ناشی از بارندگی در مسیل رود در واقع از اقدامات ایمنی، بهداشتی و رفاهی به حساب میآیند (قنواتی و همکاران، 1391: 121).
بدیهی است که عوامل اقلیمی و عوامل حوضه‌ای در بروز سیلابها تأثیر دارند، اما برای جلوگیری از آثار زیانبار سیل نمی‌توان در عوامل و عناصر جوی تغییری ایجاد نمود. بنابراین هرگونه راه‌حل اصولی و چاره‌ساز را باید در روی زمین و اختصاصاً در حوضههای آبخیز جستجو کرد. در این ارتباط، اولین اقدامی که برای کاهش خطر سیل مطرح میشود، مهار سیل در سرچشمه آن یعنی زیرحوضههای آبخیز است. مسلماً برای انجام دادن اینکار شناسایی مناطق سیل خیز در داخل حوضه لازم است. تولید روانآب در یک حوضه آبخیز به عوامل متعددی بستگی دارد که از جمله آنها میتوان به خصوصیات حوضه آبخیز، میزان بارش، نفوذ و شرایط پیشین حوضه اشاره کرد (سینگ[1]، 1997: 18).
تجارب مطالعاتی و مدیریتی کشورهای مختلف نشان میدهد که اولین گام در جهت کاهش آثار زیانبار سیل، شناخت مناطق سیل‌گیر و پهنهبندی این مناطق از لحاظ میزان خطر سیل‌گیری است تا بتوان بر اساس نتایج بهدست آمده، درباره نحوه استفاده از اراضی و کاربری‌های مختلف کشاورزی، صنعتی، خدماتی و مکان‌یابی سکونتگاه‌های شهری و روستایی تصمیمگیری اصولی و بهینه نمود و آثار زیانبار سیل را تا حد ممکن به حداقل رسانید. یکی از ابزارهای بسیار مؤثر در مطالعات محیطی و علوم زمین، استفاده از فناوری سنجش از دور و بهرهگیری از دادههای ماهوارهای است (علوی پناه، 1382: 15).
یکی از مخاطرات طبیعی که شهر ماکو در طی سالهای اخیر با آن مواجه بوده، پدیده سیل است. توسعه فیزیکی نامناسب شهر بر روی مسیلهای قدیمی، بهعلت نبود فضاهای مناسب برای ساخت و سازها، هر مسیری را برای جریان روان آب سطحی شهر مسدود کرده است. این امر سبب میشود که هنگام وقوع بارندگی، روان آب حاصل از ارتفاعات اطراف، بهویژه از قسمت شمالی شهر، بهصورت سیل وارد بافت کالبدی شهر شده و خسارت‌های زیادی بهبار آورد.
شهر ماکو در دره‌ای با بستر سنگی سخت و نفوذناپذیر در داخل حوضه آبریز زنگمارچای استقرار یافته است. بیش از 578 کیلومتر مربع از مساحت این حوضه، در بالادست شهر ماکو جا گرفته است و بهدلیل خصوصیات فیزیوگرافی خاص حوضه آبریز و وجود دامنه‌های پرشیب کوهستانی مسلط به شهر در ضلع شمالی و جنوبی تحت تأثیر سیلهای مخرب بهاری و تابستانه قرار دارد. بهعنوان مثال در 17 مرداد ماه 1387، در نتیجه یک بارش رگباری با شدت 20 میلیمتر در 10 دقیقه، در شهر ماکو سیلی جاری شد که باعث ایجاد خسارت بر 350 واحد مسکونی، 700 دستگاه خودرو، 58 واحد صنفی، آب گرفتگی واحدهای مسکونی 150 خانوار و تخریب کلی 22 واحد مسکونی شد. در این جریان 82 نفر مصدوم شدند و 4 نفر نیز جان خود را از دست دادند (هفته نامه آوای ماکو 31 مرداد 1387، شماره 209). مسؤولان اجرایی شهرستان ماکو میزان خسارت ناشی از این سیل را 82 میلیارد و 853 میلیون ریال اعلام کردند. بارش‌های رگباری سالهای 1387 و 1390 و تابستان 1391 که منجر به آبگرفتگی شدید معابر عمومی در ماکو شد، بیانگر شایع بودن وقوع سیل در این شهر و وجود ابهامات گسترده در مسایل مرتبط با مدیریت کمی و کیفی سیلابهای شهری است و برای چاره‌اندیشی آنها و مدیریت یکپارچه شهری و نیز توسعه پایدار، اجرای طرح‌های تحقیقاتی اجتناب ناپذیر است. لذا پژوهش حاضر با هدف شناسایی مناطق مستعد تولید روانآب، پهنهبندی مکان‌های در معرض خطر سیل و یافتن راهکارهای مناسب جهت کاهش آثار زیانبار سیل در شهر ماکو به انجام رسید.
در ارتباط با شناسایی و پهنهبندی مناطق سیل‌گیر با استفاده از روشهای سنجش از دور و GIS و تجزیه و تحلیل میزان خطر، پژوهش‌های زیادی در داخل و خارج انجام گرفته است. در سال 1994 سوان وراکامتون[2] با استفاده از مدل HEC-1 و سیستم اطلاعات جغرافیایی (GIS) آثار تغییر کاربری اراضی بالادست حوضه آبخیز را روی الگوی سیلاب در نواحی پایین دست حوضه مورد ارزیابی قرار داد. نتایج تحقیق نشان داد که با کاهش مساحت جنگل، روانآب حوضه اصلی و زیرحوضه‌ها افزایش یافته و تغییرات کاربری اراضی در بالادست حوضه در بهوجود آمدن سیل در اراضی پایین دست حوضه تأثیر اساسی داشـته است. فرانسیسکو[3] و همکاران (1999) بهمنظور اتصال[4] سیستم GIS با مدلهای هیدرولوژی و هیدرولیکی، مناطق تحت تأثیر سیل را در رودخانهها به نقشه در آوردند. استفان[5] (2002) در تحقیقی در حوضه رودخانه سنگزرد[6] واقع در ایالت مونتانای امریکا به بررسی سیلهای به وقوع پیوسته در سالهای 1996 و 1997 پرداخت و نهایتاً پهنههای سیل 100 و 500 ساله را تعیین کرد و نتیجه گرفت که تغییر کاربری اراضی بیشترین تأثیر را در ایجاد سیلاب دارند. کلدیتز[7] و همکاران (2003) در حوضه پانکو در مکزیک اقدام به تهیه نقشه پهنهبندی سیل را با تلفیق اطلاعات سنجش از دور با روش ژئومورفولوژیکی تهیه کرد. طبق نظر ایشان تلفیق اطلاعات سنجش از دور میتواند نقشههایی با توان تفکیک بالا تهیه کند که برای پهنهبندی سیل در مناطق وسیع میتواند بسیار مفید باشد. فیلیپ[8] و همکارانش (2006: 1180) با مطالعه رکوردهای تاریخی درازمدت از بارندگی، روان آب و دیگر عوامل اقلیمی را برای تحقیق در تمایلات و تغییرات هیدرولوژیکی در حوضه آبریز رودخانه وُلتا[9] (واقع در غرب آفریقا) در دوره زمانی 1901 تا 2002 بررسی کردهاند. یافتههای تحقیق آنها نشان داد که عوامل اقلیمی (شدت بارندگی، نوع بارش و دمای هوا) تأثیر مهمی در چگونگی سیل (شدت و دوام) در ارتباط با ویژگیهای هیدرولوژیکی حوضه (شکل حوضه، جنس زمین، شیب و کاربری اراضی) دارد.
در ایران نیز وهابی (1376) با بهکارگیری روشهای سنجش از دور و GIS، برای پهنهبندی خطر سیل در حوضه آبریز طالقان اقدام کرد. طیبی (1384) با استفاده از مدلهای منطق بولین و فازی در GIS، به مطالعه و پهنهبندی توان سیلخیزی با توجه به ویژگیهای ژئومورفیک حوضه آبریز سبزوار پرداخته است. معیری و انتظاری (1388: 1) علل و عوامل مؤثر بر ایجاد سیلابها و میزان خسارت ناشی از آنها را در شهرستانهای پرخطر استان اصفهان مطالعه کردهاند. کردوانی و همکاران (1388: 1) در تحلیل منطقهای سیلابهای لارستان به این نتیجه رسیدند که مهمترین عوامل وقوع سیلابهای لار شدت و مدت بارش و همچنین توپوگرافی منطقه است و اصلیترین مسیل لار که به بافت قدیم شهر آسیب میرساند، مسیل دربند است مسیل تنگ اسد نیز به بافت جدید شهر آسیب میزند. قهرودی تالی (1388: 167) مدل یکپارچه سیلاب شهری را در بخشهایی از شمال و شمالشرق کلانشهر تهران مطالعه کرده و نشان داد که برای کاهش سیلاب در تهران، میتوان سیلابهای بالادست را کنترل نمود. در این راستا چهار نقطه را برای مدیریت سیلابها مکانگزینی کردهاند. همچنین از نظر ایشان به این دلیل که شبکه ارتباط شهری بدون توجه به جهت زهکشی طبیعی تهران ایجاد شدهاند؛ مسیلهای تهران را از هدایت سیلابها ناتوان ساخته است. لذا این وظیفه را گذرگاههای شهری مانند اتوبانها، خیابانها و کوچهها انجام میدهند. نصری و سبزی (1390) سـیلابها و شبکه مسیلهای تأثیرگذار به شـهر زواره (شهرستان اردستان استان اصفهان) را با توجه به برنامهریزی شهری مورد مطالعه قرار داده و نشان دادند که با شناخت وضعیت شبکه مسیلها، دبی جریان و ویژگیهای مورفومتری حوضه آبریز مشرف به شهر، میتوان برای ایمنسازی شهر را برنامهریزی کرد. در جدیدترین مطالعات نیز بومری و همکاران (1390: 129) پهنههای سیلابی و ویژگیهای فیزیوگرافی و کمی حوضه آبریز دامن (در استان سیستان و بلوچستان) را با استفاده از سیستم اطلاعات جغرافیایی و سنجش از دور مطالعه کرده و نتیجه گرفتند که بهرهگیری از تصاویر ماهوارهای برای مناطقی که فاقد آمار و اطلاعات کافی در خصوص پدیده سیل هستند، بسیار مفید است و با تفسیر تصاویر با وضوح بالا میتوان پهنههای آسیبپذیر از سیل را شناسایی و طبقهبندی کرد.
با عنایت به اهمیت موضوع مدیریت سیلهای سالیانه در شهر ماکو و نبود مطالعات در این راستا، پژوهش حاضر ضمن استفاده از تجربیات مطالعات پیشین در مورد بهکارگیری سنجش از دور و GIS، بهخصوص در مناطقی که کمبود دادههای آماری وجود دارد، به شناسایی مناطق پر خطر و مدیریت بهینه سیلابها میپردازد.
شناخت عرصه پژوهش
ماکو از شهرهای شمالغربی استان آذربایجان غربی است و از نظر موقعیت مطلق جغرافیایی مابین 39 درجه و 16 دقیقه و 37 ثانیه تا 39 درجه و18 دقیقه و 3 ثانیه عرض شمالی و 44 درجه و26 دقیقه و 28 ثانیه تا 44 درجه و33 دقیقه و 29 ثانیه طول شرقی قرار گرفته است. این شهر بر اساس سرشماری سال 1390، دارای 52302 نفر جمعیت و 22842 خانوار بوده است. از نظر موقعیت استقرار، شهر ماکو در تنگه‌ای به طول 14 کیلومتر در بین دو رشته کوه مرتفع «قَیَه داغ» در شمال و «سبد داغی» در جنوب واقع شده است. این دو رشته کوه به موازات هم و با فاصلهای در حدود 500 متر از هم قرار گرفتهاند. رودخانه زنگمار تنها مجرایی است که در داخل بافت کالبدی شهر از غرب به شرق جریان مییابد. طول این رودخانه در داخل کالبد فیزیکی شهر در حدود 5/12 کیلومتر بوده ولی عرض آن متفاوت است و کل آب‌های سطحی شهر و حوضه 57845 هکتاری بالادست خود را زهکشی میکند. بستر و کناره‌های این رودخانه از سنگهای بازالتی حاصل از آتشفشان کوه آرارات شکل یافته است و پیچوخم‌های زیادی دارد. با عنایت به اینکه بیشتر روان آب‌های سطح حوضه در بدو ورود به شهر، وارد مجرای رودخانه زنگمار شده و از آن خارج میشود ولی اغلب سیلهای مخرب شهر ماکو از دامنه‌های جنوبی «قیه داغ» و مسیل «دانالو» وارد شهر میشود و خسارت‌های سنگینی را به بافت شهر تحمیل می‌نماید، این قسمت از حوضه مورد مطالعه و بررسی قرار گرفت (اشکال 1و 8) و (جدول 1).
شکل (1) موقعیت منطقه مورد مطالعه و طبقات ارتفاعی آن
مواد و روشها
برای انجام دادن این تحقیق با توجه به اهداف و ماهیت آن، از روشهای کتابخانهای و بازدیدهای میدانی استفاده کردیم. بدین ترتیب که در مراحل مختلف تحقیق، بازدیدها و برداشتهایی از منطقه مورد مطالعه به عمل آمده است. بهطور کلی برای انجام دادن پژوهش حاضر از دادهها و اطلاعات زیر استفاده شده است:
ـ تصاویر ماهواره SPOT سال 2005؛
ـ نقشههای توپوگرافی منطقه ماکو با مقیاس 50000 :1؛
ـ نقشه زمینشناسی (برگ ماکو) با مقیاس 250000 :1؛
ـ نقشه ارزیابی منابع و قابلیت اراضی استان آذربایجان غربی در مقیاس 250000 :1؛
ـ دادههای ایستگاه‌های هواشناسی، بارانسنجی و دبیسنجی منطقه ماکو؛
ـ یک دستگاه سامانه موقعیتیاب جهانی[10]
تصاویر و نقشههای تهیه شده در بدو امر آمادهسازی و زمین مرجع شدند. تجزیه و تحلیل اطلاعات تهیه شده را با شناسایی و جداسازی محدودهای از حوضه آبریز رودخانه زنگمارچای که در بالادست ساختگاه شهر ماکو واقع شده، آغاز کردیم. در این راستا از قابلیتهای نرمافزارهای Arc/map، spatial analyst plus، Arc/hydro و سامانه مدیریت حوضههای آبریز (WMS)[11] استفاده کردیم. به این منظور ابتدا DEM حوضه را ایجاد و محدوده مورد نظر را در شش واحد هیدرولوژیکی مجزا، تجزیه و تحلیل و ارزیابی کردیم.
یکی از معیارهای اساسی در تولید روان آب حوضه، زمان تمرکز است. که بر اساس فرمولهای تجربی مختلف قابل محاسبه است. در محدوده مورد مطالعه برای محاسبه آن از معادله کرپیچ استفاده شده است:
رابطه (1)
در معادله فوق
tc = زمان تمرکز بر حسب ساعت
L = طول مسیر حرکت آب در داخل حوضه (کیلومتر)
H = اختلاف ارتفاع بین نقطه تمرکز و بلندترین قسمت حوضه (متر)
میزان روان آب در یک حوضه را میتوان متناسب با خصوصیات فیزیکی آن و شدت بارندگی، بر اساس فرمولهای تجربی بر حسب ارتفاع یا حجم روان آب محاسبه کرد. یکی از فرمولهای کاربردی روش SCS است.
رابطه (2)
R = ارتفاع روان آب بر حسب اینچ
P = ارتفاع بارندگی (اینچ)
S = عامل مربوط به نگهداشت آب در سطح زمین است که مقدار آن برابر است با:
رابطه (3)
CS = مقدار نفوذ آب در حوضه است که بر اساس نوع خاک و نوع اراضی از روی جداول تعیین میشود (علیزاده، 1377: 370).
برای بررسی وضعیت پوشش گیاهی منطقه، با بهکارگیری رابطه زیر، شاخص NDVI از تصویر ماهوارهای SPOT استخراج شد:
رابطه (4)
بر اساس این شاخص، تراکم پوشش گیاهی بین 1 و 1- تعیین میشود. هرچه عدد به سمت 1 سوق پیدا کند، نشانگر تراکم زیاد پوشش گیاهی و هر چه به سمت 1- پیش رود، بیانگر کمبود و ضعف پوشش گیاهی است.
برای استخراج نقشه کاربری اراضی منطقه مورد مطالعه از رهیافت طبقهبندی نظارتشده[12] به روش حداکثر احتمال[13] تصاویر ماهوارهای استفاده شد. این روش واریانس و کوواریانس کلاس‌ها را ارزیابی میکند (علوی پناه، 1382 : 311).
خصوصیات فیزیوگرافی منطقه مورد مطالعه
جهت تحلیل ویژگیهای فیزیوگرافی حوضه و تعیین مرز زیرحوضه‌ها، و نیز برای درک بهتر رفتار هیدرولوژیکی آنها در برابر رگبارها، با استفاده از قابلیتهای نرمافزاری Arc/map و spatial analyst plus و Arc/hydro و WMS، ابتدا DEM حوضه ایجاد شد. سپس محدوده مورد نظر در شش واحد هیدرولوژیکی مجزا، مورد تجزیه و تحلیل قرار گرفت (جدول4). نتایج نشان میدهد که بیش از 11 درصد حوضه شیب 25 درجه و بیشتر از آن را دارد و اکثر مناطق پرتگاهی در اطراف شهر ماکو قرار گرفتهاند (جدول 1 و 2 و اشکال2 و 3).
جدول (1) دامنه تغییرات شیب در زیرحوضههای مشرف به شهر ماکو
|
دامنه شیب(درجه) |
5 - 0 |
14 – 6 |
24 – 15 |
37 – 25 |
70 - 38 |
کل |
|
سطح اشغال شده(هکتار) |
6/689 |
8/881 |
751 |
9/367 |
9/125 |
2/2816 |
|
درصد مساحت |
5/24 |
3/31 |
7/26 |
13 |
5/4 |
100 |
شکل (2) نقشه واحدهای هیدرولوژیکی منطقه مورد مطالعه
شکل (3) نقشه شیب حوضه زنگمار در بالا دست شهر ماکو
جدول (2) دامنه تغییرات شیب و میزان سطح اشغال آن در حوضه زنگمار (دست شهر ماکو)
|
درصد مساحت |
سطح اشغال شده (هکتار) |
دامنه شیب (درجه) |
|
22/31 |
18060 |
5 - 0 |
|
12/36 |
20890 |
14 – 6 |
|
56/21 |
12473 |
24 – 15 |
|
61/8 |
4980 |
37 – 25 |
|
5/2 |
1440 |
70 - 38 |
|
100 |
57843 |
کل |
وضعیت پوشش گیاهی
کسب اطلاعات درباره وضعیت پوشش گیاهی، از قبیل میزان و پراکنش آنها، بهلحاظ تأثیری که در مورفوژنز و مورفودینامیک منطقه می‌گذارد، اهمیت زیادی دارد. لذا مطالعه آن میتواند اطلاعات با ارزشی از فرایندهای فرسایشی و تحول منطقه ارائه دهد. تراکم پوشش گیاهی در منطقه یکسان نیست بلکه با توجه به عوامل متعدد مانند ارتفاع، دما، بارش، نوع خاک، توپوگرافی و اقلیم از پراکندگی خاصی دارد. گردآوری اطلاعات در مورد تغییرات پوشش گیاهی با روشهای معمولی بسیار مشکل و پرهزینه است. لذا استفاده از دادههای دورسنجی امکان مطالعه پوشش گیاهی را فراهم می‌آورد (شکل 4).
نوع اراضی
مطالعات میدانی و بررسی نقشههای موجود نشان داد که نوع غالب اراضی در این حوضه، کوهستانی است و در پای دامنه‌های آن ریزش واریزه نیز دیده میشود، اما اراضی پست و هموار دشت‌های سیلابی وسعت کمتری دارند (شکل5).
وضعیت خاک حوضه
خاک عمده حوضه مورد مطالعه با توجه به نوع و کیفیت خاک‌ها در سیستم فائو، بیشتر، از نوع خاکهای رگوسول[14] و لیتوسول[15] است و در اغلب نقاط با محدودیت عمق خاک، شوری و فرسایشپذیری همراه است. نقشه خاک منطقه بر اساس طبقهبندی فائو در شکل(6) نشان داده شده است.
|
|
شکل(4) (الف) تصویر ماهوارهای SPOT-5 و (ب) NDVIمنطقه مورد مطالعه
کاربری اراضی
روشهای مختلفی برای طبقهبندی تصاویر ماهوارهای بهکار می‌روند که از میان آنها طبقهبندی نظارت شده[16] بهروش حداکثر احتمال[17] تاکنون بهعنوان دقیقترین و پراستفادهترین روشها ذکر شده است. این روش واریانس و کوواریانس کلاس‌ها را ارزیابی میکند (علوی پناه، 1382: 311). نقشه کاربری اراضی حوضه مورد مطالعه با استفاده از تصویر ماهوارهای Spot-5 و پس از محاسبات آماری (جدول3 و 4) جهت ایجاد FCC ، در محیط نرمافزاری ENVI با روش حداکثر احتمال استخراج شد. بر اساس این تحلیلها، تقریباً 40 درصد منطقه را انواع مرتع متراکم و نیمهمتراکم تشکیل میدهد، 24 درصد اراضی بدون کاربری است و 17 درصد آن را برونزد سنگی پوشش میدهد (شکل 7).
جدول (3) پارامترهای آماری باندهای طیفی
|
Basic Stats |
Min |
Max |
Mean |
Stdev |
|
Band 1 |
0 |
254 |
01/50 |
87/51 |
|
Band 2 |
0 |
254 |
94/34 |
24/39 |
|
Band 3 |
0 |
254 |
19/39 |
38/42 |
|
Band 4 |
0 |
254 |
37/49 |
93/51 |
جدول (4) ضریب همبستگی بین زوج باندها
|
Correlation |
Band 1 |
Band 2 |
Band 3 |
Band 4 |
|
Band 1 |
1 |
88/0 |
92/0 |
93/0 |
|
Band 2 |
88/0 |
1 |
99/0 |
97/0 |
|
Band 3 |
92/0 |
99/0 |
1 |
97/0 |
|
Band 4 |
93/0 |
97/0 |
97/0 |
1 |
شکل (5) نقشه نوع اراضی منطقه مورد مطالعه (منبع: نقشه زمین شناسی ماکو و نقشه قابلیت اراضی آذربایجان غربی)
شکل (6) نقشه نوع خاک منطقه مورد مطالعه (منبع: نقشه زمینشناسی ماکو و نقشه قابلیت اراضی آذربایجان غربی)
شکل (7) نقشه کاربری اراضی منطقه مورد مطالعه مستخرج از تصویر ماهوارهای SPOT – 2005
یافتههای تحقیق
ارزیابی میزان روانآب سطوح طبیعی بهعلت تأثیرات متقابل بین توپوگرافی، عمق آب زیرزمینی، شدت و تداوم بارش و خصوصیات نفوذپذیری که بهوسیله خاک و پوشش گیاهی تعدیل میشوند، پیچیده‌تر است زیرا تمام سطح حوضه در تولید روانآب نقش یکسانی ندارند. به این ترتیب، نقش شرایط اولیه حوضه در پتانسیل تولید روانآب اهمیت پیدا میکند (Chow and et al; 1988: 532). در خصوص مساحت حوضه و ارتباط آن با مقدار روانآب تولید شده، برخی از محققان، وسعت مساحت را به افزایش حجم روانآب (صدقی؛ 1363: 14-2) و بعضی دیگر افزایش مساحت را به فزونی دبی منوط دانستهاند (Singh, 1997: 18). در منطقه مورد مطالعه وجود گستره وسیعی از سطوح نفوذناپذیر و شیب تند آنها شرایط ایجاد روان آب، و سیلهای مخرب به هنگام رگبارهای زودگذر را فراهم می‌آورد.
جریان سطحی هورتونی بهعنوان یکی از انواع اصلی تولید روانآب و سطوح مشارکتی، زمانی بهوجود می‌آید که شدت بارش بیش از نرخ نفوذ باشد. در این حالت، پس از تأمین ذخایر گودالی جریان سطحی بهوجود می‌آید. معمولاً جریان سطحی هورتونی از سطوحی با خاکهای رسی و کم عمق، سطوح عاری از پوشش گیاهی، سطوح شهری و سطوح با خاکهای فشرده که ظرفیت نفوذ را کاهش میدهند (مانند جادهها و امثال آن) تولید میشود (Stephen, 1999: 1-18). بررسی‌های بهعمل آمده نشان داد که جنس تولید روان آب در شهر ماکو، هورتونی است. قسمت‌هایی از شهر که مابین کوه «قیه داغ» و رودخانه زنگمار و پایین‌تر از «قره کورپی» واقع شده است، به تناسب شرایط نفوذ، شیب توپوگرافیک و نیز وسعت حوضه بالادست، بیش از سایر قسمت‌ها در معرض سیلهای مخرب قرار دارند. اما قسمت‌های دیگر شهر به دلیل ورود جریانها به داخل شاخههای رودخانه زنگمار، چندان آسیبپذیر نیستند (شکل8).
شکل (8) تصاویر google earth از محدوده مورد مطالعه
نمودارهای میانگین بارشها و میانگین حداکثرها، نشان میدهند که اولاً در چند سال اخیر نوسانات بارش سالانه غیر خطی است، ثانیاً بیشترین میزان بارشها در بهار و حداقل آن در تابستان اتفاق می‌افتد و در این فصل پوشش گیاهی آسیبپذیر میشود و نقش حفاظتی خود را از دست میدهد. لذا وقوع یک بارش رگباری و مازاد بر ظرفیت نفوذ، شرایط وقوع سیلابهای مخرب را تسهیل می‌نماید (اشکال 10-9). حجم بارش و روان آب حاصل از این بارش رگباری در هر یک از واحدهای هیدرولوژی در جدول 5 آمده است.
شکل (9) نمودار میانگین بارش سالانه، ایستگاه ماکو87-64
شکل (10) نمودار میانگین بارش‌های حداکثر ماهانه، ایستگاه ماکو87-64
خصوصیات نفوذپذیری حوضه که با نام CN یا شماره منحنی معرفی میشوند (علیزاده، 1377: 363)، ارتفاع روان آب و حجم روانآب در قسمت‌های مختلف حوضه، در جدول 5 و 6 آمده است. برآوردها، حجم روان آب حاصل از یک بارش رگباری 20 میلیمتری در 10 دقیقه را، در حدود 5/2 میلیون مترمکعب نشان میدهد و از این میزان، 338 هزار متر مکعب مستقیم وارد محدوده شهر ماکو میشود و آن را آسیبپذیر می‌سازد. نکته قابل توجه اینکه از حوضه بالادست صخره سنگی معروف شهر در دامنه «قیه داغی» که مساحتی در حدود 200 هکتار دارد، 28 هزار متر مکعب آب وارد شهر میشود (شکل 12).
بررسی‌ها نشان دادند که درحوضه مورد مطالعه 5/21 درصد بارشها تبدیل به روانآب میشود. شکل (11) پهنهبندی مناطق تولید روانآب و در معرض سیل را نشان میدهد.
جدول (5) خصوصیات فیزیکی و هیدرولوژیکی زیرحوضههای رود زنگمار در بالا دست ماکو
|
حجم روانآب (m3) |
ارتفاع روانآب (mm) |
ارتفاع بارندگی (Inch) |
CN متوسط |
حجم بارش (m3) |
مساحت (هکتار) |
زمان تمرکز (Tc) ساعت |
طول آبراهه (L) کیلومتر |
اختلاف ارتفاع (H)متر |
نام زیرحوضه مورد مطالعه |
|
808650 |
2/3 |
78/0 |
80 |
50540600 |
3/25270 |
94/0 |
10 |
1020 |
SB1 |
|
363000 |
7/3 |
78/0 |
85 |
15276000 |
7638 |
76/0 |
9 |
1280 |
SB2 |
|
211514 |
7/3 |
78/0 |
81 |
11433200 |
6/5716 |
7/0 |
7/8 |
1480 |
SB3 |
|
134241 |
3 |
78/0 |
77 |
8949400 |
7/4474 |
75/0 |
5/9 |
1580 |
SB4 |
|
329038 |
4/3 |
78/0 |
78 |
19355200 |
6/9677 |
95/0 |
12 |
1700 |
SB5 |
|
152031 |
3 |
78/0 |
80 |
10135400 |
7/5067 |
5/0 |
7/5 |
960 |
SB6 |
|
2493344 |
3/3 |
78/0 |
80 |
11568600 |
57845 |
5/5 |
55 |
1780 |
کل حوضه* |
* محاسبات به صورت وزنی و توسط نرم افزار صورت گرفته است.
شکل (11) نقشه پهنههای تولید روانآب در محدوده مورد مطالعه
دامنه‌های جنوبی «قیه داغی» که با شیب بسیار تندی به ساختگاه شهر ماکو مشرف است، در مجموع 2/2816 هکتار مساحت دارد. این بخش به چندین زیرحوضه کوچک تقسیمبندی شده است و مقدار CN آن، بهعلت پوشش سنگی و سنگلاخی بودن زمین در اکثر سطوح، رقم 90 تعیین گردید. لذا بر اساس روابط فوق، پارامترهای مربوط محاسبه، و در جدول (6) آورده شده است.
میزان تبدیل بارش به روانآب، برای زیرحوضه مشرف به شهر 60 درصد برآورد شده است.
جدول (6) خصوصیات فیزیکی و هیدرولوژیکی زیرحوضههای مشرف به شهر ماکو
|
حجم روان آب (m3) |
ارتفاع روان آب (mm) |
ارتفاع بارندگی (Inch) |
حجم بارش (m3) |
مساحت (هکتار) |
زمان تمرکز (Tc) ساعت |
طول آبراهه (L)کیلومتر |
اختلاف ارتفاع (H) متر |
زیر حوضه مورد مطالعه |
|
118440 |
12 |
78/0 |
197400 |
987 |
5/0 |
2/5 |
680 |
SB1 |
|
116076 |
12 |
78/0 |
193460 |
3/967 |
48/0 |
05/5 |
740 |
SB2 |
|
26760 |
12 |
78/0 |
44600 |
223 |
31/0 |
2/3 |
580 |
SB3 |
|
24600 |
12 |
78/0 |
41000 |
205 |
16/0 |
95/1 |
740 |
بشگوز |
|
27960 |
12 |
78/0 |
46600 |
233 |
13/0 |
45/1 |
500 |
کوه قیه |
|
337944 |
12 |
78/0 |
563240 |
2/2816 |
48/0 |
2/5 |
800 |
کل* |
* محاسبات به صورت وزنی و توسط نرمافزار صورت گرفته است.
شکل (12) نقشه پهنههای در معرض خطر سیل و آبگرفتگی در شهر ماکو
نواحی آسیبپذیر شهر ماکو در برابر سیل بهشرح زیر مشخص شد و باید عملیات حفاظتی در محلات زیر صورت گیرد (شکل 12):
قسمت مرکزی و بافت قدیمی (هسته اولیه) شهر، در دامنه‌های کوه «قیهداغ»، که جاری شدن سیل به صورت آبشار از بالای کوه به داخل شهر، بافت مسکونی و بافت تجاری اصلی شهر را تهدید میکند (شکل 13).
بافت شهر در دامنه‌های شرقی کوه «قیه داغ» در محله «بشگوز» در روی مسیلهای قدیمی و قسمت شمالی معبر اصلی شهر قرار گرفتهاند. معابر ارتباطی این واحدها، همان مسیلهای قدیمی است که بدون توجه به خطرات عبور سیل، مورد ساخت و ساز واقع شده و تبدیل به معبر گشته‌اند (شکل 14).
قسمت قدیمی شهر در سمت غربی بافت اولیه شهر و مقابل ساختمان شهرداری ماکو موسوم به «محله» که در دل کوه واقع شده است (شکل 15).
بافت تقریباً جدید شهر در سمت غربی کوه «قیه داغ» موسوم به «محله مسکن» که دقیقاً در مسیل کناری راه ارتباطی روستاهای سمت شمالی شهر (روستای دانالو) قرار دارد. این قسمت از بافت شهر در روی مسیل منتهی به رودخانه زنگمار قرار گرفته و آن را مسدود نموده است و تنها محل عبور برای سیلابها، معبر ارتباطی متصل به بلوار اصلی شهر است. سیل جاری شده با ورود به بلوار اصلی شهر بهویژه قسمت شمالی آن، در جهت شرق به طرف بخش مرکزی شهر جریان مییابد. در طول این معبر که تنها معبر اصلی و ارتباطی عمده شهر است، تمامی واحدها، بهویژه واحدهای سمت جنوبی (در جهت شیب منتهی به رودخانه زنگمار) در معرض سیل واقع میشوند.
|
|
شکل (13) جاری شدن سیل از بالای کوه «قیهداغ» به داخل شهر (منبع: آرشیو)
شکل (14) نفوذ سیل به داخل بخش مرکزی و تجاری(منبع: آرشیو)
شکل (15) نفوذ سیل و تخریب پارک در محله «بشگوز» و «محله» (منبع: آرشیو)
نتیجهگیری
پراکنش نامناسب زمانی و مکانی بارش در مناطق خشک و نیمهخشک، علاوه بر ایجاد سیلابهای مخرب، موجب هدر رفتن روانآب سطحی میشود. از این رو مهار روانآب سطحی و بهرهبرداری مناسب از آن میتواند راهکار مناسبی برای جلوگیری از هدر رفتن آب، و تبدیل تهدید‌ها به فرصت‌ها باشد. تعیین عرصه‌های مستعد تولید روانآب و سیل و مناطق در معرض آبگرفتگی یکی از مهمترین عوامل تعیینکننده در موفقیت طرح‌های سیلخیزی است. از سوی دیگر، تعیین مکان‌های مناسب با استفاده از روشهای سنتی و متداول بسیار دشوار بوده، اکثراً باعث بروز خطا میشود. سیستم اطلاعات جغرافیایی به همراه ابزارها و تکنیک‌های آن توانایی آن را دارد تا با تلفیق لایههای مختلف اطلاعاتی در قالب مدلهای مختلف و در حداقل زمان ممکن در مکان‌یابی این عرصهها، مورد استفاده قرار گیرد.
در این مطالعه، علاوه بر پهنهبندی مناطق سیلخیز حوضه، نقاط سیل گیر در داخل شهر ماکو نیز مورد ارزیابی قرار گرفت و مشخص شد که در محدوده مورد مطالعه، بخش کوچکی از حوضه که در شمال شهر با شیب تند و دارای مواد منفصل ناچیز وجود دارد، بیشترین نقش را در تولید سیلهای مخرب و تهدیدکننده شهر ماکو، ایفا میکند که میتوان با احداث کانال‌های انحرافی، جریان آب را در فواصل مختلف در طول شهر، به رودخانه زنگمار هدایت و با یک سد مخزنی در پایین دست، جهت استفاده در مصارف کشاورزی ذخیره نمود. از طرف دیگر نبود زمین لازم برای گسترش شهر و مشکلات مربوط به توسعه فیزیکی آن سبب شده است که تـمامی عرصه‌های موجود به اشغال کاربـری‌های مختلف شهری درآید. وجود دامنه‌های تند و تیز، برونزدهای سنگی و زمین‌های فاقد پوشش سطحی، محدودیت فضا و بارش‌های تند رگباری در فصل خشک، بروز سیلهای مخرب را در این ناحیه مساعدتر می‌سازد و عدم رعایت حریمها و تجاوز به بستر مسیلها و آب گذرها، آن را شدت می‌بخشد و خسارت‌های زیادی را در محله‌های «بشگوز»، «محله»، «مسکن» و بخش قدیمی شهر (هسته اولیه) واقع در پایین دست کوه «قیه داغ» به بار می‌آورد.
از آنجا که نفوذپذیری از عوامل مؤثر در ایجاد روان آب و سیل است و نخستین اثر توسعه کالبد فیزیکی شهر، کاهش نفوذپذیری است؛ باید از هرگونه عملیات توسعه شهری و ساخت و سازها در حریمها جلوگیری شود و با آزادسازی و تبدیل فضاهای اشغالی به فضاهای سبز و مهیا نمودن شرایط نفوذ، جریان روان آب و سیلهای ویرانگر را کاهش داد. همچنین میتوان اندازه مسیر عبور آب‌های سطحی (کانال) را با وسعت حوضه بالادستی آن متناسب ساخت و با احداث کانال‌های انحرافی بزرگتر، جریان سیل را به سمت مجاری اصلی هدایت نمود. لذا یکی از راهکارهای عملی کوتاهمدت، در جهت کاهش آثار زیانبار سیل در شهر آزادسازی حریمها، و در دراز مدت، ایجاد شهرک‌های جدید و انتقال جمعیت، پیشنهاد میشود.
[1]- Singh
[2]- Swan werakamton
[3]- Francisco
[4]- Coupling
[5]- Stephan
[6]- Yellowstone
[7]- Colditz
[8]- philip
[9]- Volta River Basin
[10]- Global Positioning System (GPS)
[11]- Watershed Management System
[12]- Supervised classification
[13]- Maximum likelihood
[14]- Reglesols
[15]- Lithosols
[16]- Supervised classification
[17]- Maximum likelihood
)