ارزیابی حساسیت‌پذیری رخساره فرسایشی هزاردره در حوزه آبخیز فیروزکوه

مقدمهفعالیت‌های انسانی از مهم‌ترین دلایل تغییرات محیطی بوده و انسان‌ها نقش اساسی در گرمایش زمین، تخریب زمین، آلودگی آب و خاک، بالا رفتن سطح آب دریاها، تخریب لایه ازن، جنگل‌زدایی شدید و اسیدی کردن اقیانوس‌ها دارند. در مقیاس جهانی، تخریب و فرسایش یکی از مهم‌ترین پدیده‌های طبیعی است که مناظر را تغییر می‌دهد. فرسایش خاک اثرات درون‌منطقه‌ای و برون‌منطقه‌ای فراوانی را به‌همراه دارد. فرسایش و هدررفت خاک معمولاً در مناطق شیب‌دار اتفاق می‌افتد اما شدت وقوع آن به عامل‌های محیط‌زیستی از قبیل میزان شیب، ویژگی‌های زمین‌شناسی، وضعیت پوشش گیاهی و عوامل اقلیمی بستگی دارد که باعث افزایش آسیب‌پذیری یک منطقه به فرسایش می‌شود. یکی از مهم‌ترین رخساره‌های فرسایشی رخساره هزاردره یا بدلند است. رخساره فرسایشی هزاردره معمولاً به مناطق با رسوبات سست و غیرمتراکم، پوشش گیاهی بسیار ضعیف یا فاقد پوشش گیاهی که به‌دلیل تراکم بالای شیارها و بریدگی‌ها برای کشاورزی غیرقابل استفاده است، اشاره می‌کند. فرسایش هزاردره بیش‌تر در مناطق خشک و نیمه‌خشک مشاهده می‌شود و در واقع برهم‌کنش بارش با مواد زمین‌شناسی سست بسیار مساعد ایجاد و گسترش رخساره هزاردره در مناطق خشک و نیمه‌خشک است. با توجه به اینکه فرسایش خاک یک فرایند بسیار پیچیده و تحت تاثیر ویژگی‌های سطح زمین، نوع خاک و دیگر عامل‌های محیط‌زیستی است، بررسی کمی و تهیه نقشه‌های حساسیت‌پذیری فرسایش امری بحث برانگیز است. هدف اصلی از پژوهش حاضر، تهیه نقشه حساسیت فرسایش هزاردره در حوزه آبخیز فیروزکوه با استاده از مدل نسبت فراوانی است. مواد و روش‌‌هاحوزه آبخیز فیروزکوه به این دلیل به‌عنوان منطقه موردمطالعه انتخاب شد که رخساره فرسایش هزاردره مهم‌ترین نقش را در فرسایش خاک این حوزه به‌دلیل وضعیت اقلیمی، هیدرولوژیکی، توپوگرافی و کاهش پوشش گیاهی و همچنین وجود سازندهای مستعد فرسایش دارد. اولین گام در این پژوهش تهیه نقشه پراکنش رخساره‌های هزاردره و تعیین محل آنها روی نقشه است. این کار با استفاده از تصاویر گوگل ارث و پیمایش زمینی انجام شد. نقشه‌های مربوط به عوامل موثر از منابع مختلف تهیه و وارد محیط نرم‌افزار arc/gis شد. با استفاده از نقاط و خطوط ارتفاعی موجود در نقشه‌های توپوگرافی تهیه شده از سازمان نقشه‌برداری کشور، نقشه مدل رقومی ارتفاع (dem) با اندازه پیکسل 30 متر تهیه شد. نقشه‌های جهت شیب، درجه شیب، انحناء سطح، شاخص رطوبت توپوگرافی (twi) و طبقات ارتفاع با استفاده از نقشه dem و در محیط نرم‌افزارهای arcgis10.3 و saga-gis به‌دست آمد. نقشه زمین‌شناسی منطقه نیز از نقشه زمین‌شناسی ایران با مقیاس 1:100000 استخراج شد. نقشه‌های آبراهه و جاده از نقشه توپوگرافی 1:25000 استخراج شده و فاصله از این عوارض در محیط arc/gis10.3 محاسبه شد. نقشه کاربری اراضی حوزه آبخیز فیروزکوه با استفاده از تصاویر لندست 8 برای سال 2020 و با روش تلفیقی به‌دست آمد. برای تهیه نقشه مشخصات مربوط به خاک 30 نمونه خاک از عمق صفر تا 30 سانتی‌متری برداشت و در آزمایشگاه بررسی شد. نقشه بارش متوسط سالانه با استفاده از داده‌های بارش ایستگاه‌های هواشناسی تهیه شد. بعد از طبقه‌بندی نقشه‌های مربوط به عوامل موثر، وزن هر نقشه با استفاده از مدل نسبت فراوانی محاسبه شد. در مرحله بعد با جمع کردن وزن‌ها، نقشه نهایی حساسیت فرسایش هزاردره تهیه شد. منحنی roc و سطح زیر این منحنی برای ارزیابی دقت مدل نسبت فراوانی استفاده شدند. نتایج و بحثارتباط بین فرسایش هزاردره و عوامل موثر با استفاده از مدل نسبت فراوانی بررسی شد. نتایج نشان داد که بیش‌ترین وزن نسبت فراوانی مربوط به طبقه ارتفاعی 1710 تا 2286 متر، بارش 400 تا 550، شیب بیش‌تر از 35، جهت شمال غرب، فاصله کم‌تر از 1150 متر از آبراهه، سازندهای مارنی، آهکی و شیل، کاربری مرتع، شیب‌های محدب و مقعر، طبقه رس 25 تا 33 درصد، طبقه سیلت 27 تا 35 درصد، گروه هیدرولوژیک c، عمق خاک 57 تا 120 سانتی‌متر، اسیدیته 7.6 تا 8.1، طبقه 6 تا 11 شاخص twi بوده است. بررسی دقت مدل نسبت فراوانی با استفاده از منحنی roc و سطح زیر این منحنی انجام شد. سطح زیر منحنی roc 0.71 به دست آمد که نشان می‌دهد مدل نسبت فراوانی برای تهیه نقشه حساسیت فرسایش هزاردره در حوزه آبخیز فیروزکوه قابل قبول است. مدل نسبت فراوانی با وجود سادگی، به‌دلیل ایجاد ارتباط منطقی بین هزاردره‌ها و عوامل موثر بر آن نتایج قابل قبولی ارائه می‌کند. مطالعات دیگری ازجمله بررسی پتانسیل آب‌های زیرزمینی، تهیه نقشه حساسیت زمین‌لغزش و آسیب‌پذیری یک منطقه به سیلاب نیز با این مدل انجام شده و دقت آن مورد تایید قرار گرفته است. نتیجه‌‌گیریبه‌دلیل شرایط توپوگرافی، اقلیمی و زمین‌شناسی رخساره فرسایشی هزاردره یک پدیده غالب در حوزه آبخیز فیروزکوه است. در این پژوهش نقشه حساسیت فرسایش هزاردره با استفاده از مدل نسبت فراوانی تهیه شد. ارزیابی دقت نشان داد، مدل نسبت فراوانی یک مدل مناسب برای تهیه نقشه حساسیت‌پذیری فرسایش هزاردره در این حوزه آبخیز است. نتایج نشان داد که حدود 50 درصد منطقه حساسیت زیاد و خیلی زیاد به وقوع فرسایش هزاردره دارد، بنابراین توجه به این پدیده و تهیه نقشه حساسیت‌پذیری آن امری ضروری است.

badland erosion susceptibility assessment in the firozkuh watershed

introductionhuman activities are at the core of global environmental change and humans play a key role in global warming, land degradation, air and water pollution, rising sea levels, eroding the ozone layer, extensive deforestation, and acidification of the oceans. soil erosion and degradation is a natural phenomenon altering the relief of the landscape. erosion is often capable of causing several on-site and off-site impacts. erosion and soil loss are common in hilly areas, but their severity will vary depending on the geoenvironmental factors including, steep sloping, geological characteristics, vegetation and climatic factors making it more vulnerable to erosion. one of he most important kind of erosion is badland erosion. the term of badlands currently refers to areas of unconsolidated sediment with little or no vegetation, which are useless for agriculture because of their intensely dissected landscape. badland erosion is observed mostly in arid and semi-arid regions, and the interaction of precipitation with geological materials is responsible for the development of badlands in arid and semi-arid regions. because soil erosion is a complicated process that is influenced by the properties of the land surface and the soil as well as by environmental factors, quantitatively accurate forecasts of soil erosion and susceptibility mapping are challenging. the main goal of this research is to map badland erosion suscepibility in firozkuh watershed using frequency ratio model. materials and methodsfirozkuh watershed was selected as study area because in this watershed, badlands are the most important contributors to soil erosion because of the condition of climatic, hydrologic, topographic, and reduced vegetation conditions, and as well as presence of susceptible soil and geology formations in this region. the first step in this research is to prepare distribution map of the badlands and determine their location on the map. this was done using google earth imagery and field surveys. the maps of condiioning factors were prepared from different sources and entered into the gis environment. digital elevation model (dem) map with the cellsize of 30 meter was prepared using the elevation points and lines in the topographic maps prepared by the national cartography center of iran. slope aspect, slope degree, plan curvature, twi and elevation classes maps was creaed using dem map in arcgis10.3 and saga-gis environment. the geology map of the watershed was extracted from the geologic map of iran with the scale of 1:100000. river and road maps were extracted from 1:25000 topographic map and the distance from these features was calculated in arcgis10.3 environment. the land use of firozkuh watershed was created from landsat 8 images of year 2020  using a synthetic method. to map soil characteristics, 30 samples were taken from depth of 0-30 centimeter and analyzed in the laboratory. aaverage annual rainfall map was developed using rainfall data from meteorological stations. after classifying conditioning factors maps, the weight of each map was calculated using the frequency ration model. in the next step, by combining the weights, the final badland erosion susceptibility map was prepared. the roc curve and the area under the curve were used to assess the accuracy of the frequency ratio models. results and discussionthe relationship between badland erosion and conditioning factors was investigated using the frequency ratio model. the results showed that the highest weight of the frequency ratio is related to the elevation class of 1710 to 2286 meters, rainfall 400 to 550, slope more than 35, northwest aspect, distance less than 1150 meters from drainage network, marl, limestone and shale formations, ranglands, convex and concave slopes, clay 25 to 33%, silt 27 to 35%, hydrological group c, soil depth 57 to 120 cm, ph 7.6 to 1.8, twi class 6 to 11. accuracy assessmen of the freency ratio model was done using roc and area under this curve. the area under the roc curve was 0.71 that showed frequency ratio model is acceptable for badland erosion susceptibility mapping in the firozkuh watershed. despite its simplicity, the freqency ratio model provides acceptable results due to the creation of a logical connection between the badlands and conditioning factors. other studies, including investigating the potential of underground water, landslide susceptibility maping, and the vulnerability to floods, have also been conducted with this model, and its accuracy has been confirmed. conclusionsbecause of topographical, climatic and geological conditions, the badland erosion is a dominant phenomenon in the firozkuh watershed. in this research, badland erosion susceptibility map was prepared using the frequency ratio model. accuracy assessment showed that frequency ratio is a suitable model for badland erosion susceptibility maping in this watershed. the results showed that about 50% of this region has high and very high susceptibility to the badland erosion, so it is necessary to pay attention to this phenomenon and prepare a susceptibility map.