نوع مقاله : مقاله پژوهشی
نویسندگان
1 دانشجوی دکترای محیط زیست
2 دانشگاه ملایر، گروه محیط زیست
3 هیئت علمی تمام وقت
چکیده
خرس قهوهای (Ursus arctos) گونهای با دامنه پراکنش وسیع در استان کرمانشاه است و این مهم پیششرط بروز تعارض را به وجود آورده است. این مطالعه با هدف شناسایی زیستگاههای کلیدی خرس، نوسان آشیان اکولوژیک و پوشش زیستگاه مطلوب آن توسط مناطق حفاظتشده و شهرستانهای استان انجامگرفته است. در این راستا پس از پایشهای میدانی به مدت 3 سال (از سال 1395 تا 1398) نقاط حضور گونه شناسایی شدند. از 9 متغیر زیستگاهی در غالب متغیرهای اکولوژیک، توپوگرافی و انسانی برای بررسی زیستگاه بالقوه استفاده شد. مدلسازی توزیع با استفاده از الگوریتم آنتروپی بیشینه انجام گرفت. بررسی اهمیت متغیرها در مدلسازی به روش تحلیل جکنایف و اعتبار مدل نیز بهوسیله تابع ROC محاسبه گردید. کمی سازی پهنای آشیان اکولوژیک نیز با استفاده از متریکهای مربوطه در نرمافزار ENMTOOLS محاسبه گردید. بهمنظور محاسبه کمی مساحت زیستگاه از حد آستانه TSS استفاده شد. بر اساس نتایج، پوشش زیستگاه مطلوب گونه توسط مناطق حفاظتشده مناسب نیست و این ضعف در بخشهای جنوبی استان مشهود است؛ بهنحویکه مساحتی معادل 84/265842 هکتار از زیستگاه گونه در خارج از مرز مناطق حفاظتشده و مستعد تعارض با سایر کاربری/پوششهاست. در بین مناطق حفاظتشده استان، مناطق قلاجه و بوزین و مرخیل بیشترین پوشش را از زیستگاه اینگونه چتر داشتند. بر اساس نتایج، حدود نیمی از مساحت شهرستان پاوه معادل 68/46 درصد، مستعد حضور گونه بوده و کمترین احتمال پراکنش نیز در شهرستان قصر شیرین است. نتایج این مطالعه لزوم توجه بیشتر به پیوستگی و پوشش زیستگاه اینگونه در بخشهای جنوبی استان را بیشازپیش آشکار میکند.
کلیدواژهها
موضوعات
عنوان مقاله [English]
Brown bears (Ursus arctos Linnaeus, 1758), an umbrella spices Get involved conflict crisis in Kermanshah provinces
نویسندگان [English]
2 Malayer University, Department of Environmental Sciences
3 Department of Biology, Faculty of Science, Razi University, Kermanshah, Iran
چکیده [English]
Brown bear (Ursus arctos) is a species with high distribution range in Kermanshah province and this has created precondition of conflict. The aim of this study is to identify key bear habitats, ecological Niche fluctuations, and their desirable habitat cover in protected areas and counties of the province. In this regard, after field surveys, the presence of species was identified for 3 years (from 2016 to 2018). Nine habitat variables including ecological, topographic and anthropogenic were used to investigate distribution status. Distribution modeling was performed using Maximum Entropy algorithm. The significance of variables in modeling was calculated by jack knife method and model validity was calculated by ROC function. Width quantification of ecological nests was also calculated using related metrics in ENMTOOLS software. TSS threshold was used to quantify habitat area. Based on the results, appropriate habitat coverage of the species by the protected areas is not suitable and this weakness is evident in the southern parts of the province so that area of 265842.84 hectares of species habitat is outside the protected areas boundary and is prone to conflict with other Land use/land cover. Among the protected areas of the province, Ghaleje and Buzin and Marakhil had the most cover of this an umbrella spices. About half of Paveh county (% 46.68) is suitable for species and the least suitable was located in Qasrshirin county. The results of this study necessitate more attention to the continuity and habitat coverage of this species in the southern parts of the province.
کلیدواژهها [English]
خرس قهوهای (Ursus arctos, Linnaeus, 1758)، گونهی چتر درگیر
بحران تعارض در استان کرمانشاه
پیمان کرمی1، کامران شایسته1*، نصرالله رستگار پویانی2
1 ایران، همدان، دانشگاه ملایر، دانشکده منابع طبیعی و محیطزیست، گروه محیطزیست
2 ایران، کرمانشاه، دانشگاه رازی، دانشکده علوم، گروه زیستشناسی
تاریخ دریافت: 19/07/1399 تاریخ پذیرش: 28/08/1399
چکیده
خرس قهوهای (Ursus arctos) گونهای با دامنه پراکنش وسیع در استان کرمانشاه است و این مهم پیششرط بروز تعارض را به وجود آورده است. این مطالعه با هدف شناسایی زیستگاههای کلیدی خرس، نوسان آشیان اکولوژیک و پوشش زیستگاه مطلوب آن توسط مناطق حفاظتشده و شهرستانهای استان انجامگرفته است. در این راستا پس از پایشهای میدانی به مدت 3 سال (از سال 1395 تا 1398) نقاط حضور گونه شناسایی شدند. از 9 متغیر زیستگاهی در قالب متغیرهای اکولوژیک، توپوگرافی و انسانی برای بررسی زیستگاه بالقوه استفاده شد. مدلسازی توزیع با استفاده از الگوریتم آنتروپی بیشینه انجام گرفت. بررسی اهمیت متغیرها در مدلسازی به روش تحلیل جکنایف و اعتبار مدل نیز بهوسیله تابع ROC محاسبه گردید. کمی سازی پهنای آشیان اکولوژیک نیز با استفاده از متریکهای مربوطه در نرمافزار ENMTOOLS انجام گرفت. بهمنظور محاسبه کمی مساحت زیستگاه از حد آستانه TSS استفاده شد. بر اساس نتایج، پوشش زیستگاه مطلوب گونه توسط مناطق حفاظتشده مناسب نیست و این ضعف در بخشهای جنوبی استان مشهود است؛ بهنحویکه مساحتی معادل 84/265842 هکتار از زیستگاه گونه در خارج از مرز مناطق حفاظتشده و مستعد تعارض با سایر کاربری/پوششهاست. در بین مناطق حفاظتشده استان، مناطق قلاجه و بوزین و مرخیل بیشترین پوشش را از زیستگاه اینگونه چتر داشتند. بر اساس نتایج، حدود نیمی از مساحت شهرستان پاوه معادل 68/46 درصد، مستعد حضور گونه بوده و کمترین احتمال پراکنش نیز در شهرستان قصر شیرین است. نتایج این مطالعه لزوم توجه بیشتر به پیوستگی و پوشش زیستگاه اینگونه در بخشهای جنوبی استان را بیشازپیش آشکار میکند.
واژه های کلیدی: تعارض، کاربری/پوشش، خرس قهوهای، توزیع، کرمانشاه
* نویسنده مسئول، تلفن: 08334260512 ، پست الکترونیکی: ka_shayesteh@yahoo.com
مقدمه
به دنبال افزایش آگاهی از مفهوم تنوع زیستی، زیست شناسان حفاظت به راههای میانبری که بتواند تنوع کل بیوتا (Biota) را حفظ کند علاقه وافری نشان دادند و دراینبین برخی اصطلاحات و مفاهیم نیز کاربردی شدند. یکی از این مفاهیم، گونه چتر بود و در بطن خود این معنا را داشت که از طریق برآورده کردن نیازهای گونه، اکوسیستمها و زیستگاهها نیز حفظ شوند (21). تعارض بین انسان و حیاتوحش به هر نوع تعاملی گفته میشود که نتایج آن بهصورت زیانبار به مردم، حیاتوحش و منابع استفادهشده ظاهر میشود (14). رشد سریع جمعیت و بهتبع آن افزایش مصرف منابع و تخریب زیستگاه طبیعی جانوران، سبب ایجاد درگیری و تعارض بین انسان و حیاتوحش میشود. تعارض بین انسان و حیاتوحش از وجود منابع مشترک و ایجاد رقابت بر سر این منابع محدود است و زمانی که این اثرات در ارتباط با اقتصاد و معیشت مردم نگریسته میشود تبدیل به یک مسئله بحثبرانگیز میشود (1). گونههایی که دارای جثه بزرگ هستند به دلیل نیازهای زیستی بالا آسیبپذیر خواهند بود (17). اندازه جمعیت و زیستگاه گوشتخواران بزرگ دارای روندی کاهشی است و مسبب آن را میتوان اثرات انسانی، تکهتکه شدن زیستگاه و شکار دانست. در میان تمام گونههایی که در معرض انقراض قرار دارند، گوشتخواران بزرگ جثه تهدید پذیری بیشتری دارند. تهدیداتی که با آن مواجه هستند منجر به کاهش در جمعیت، محدوده توزیع جغرافیایی، زیستگاه و غذای آنها شده است که درنهایت ماحصل این تهدیدات، منجر به ارتباط آنها با انسان شده است (20). خرس قهوهای بهعنوان بزرگترین گوشتخوار ایران، گونهای بدون قلمرو گرایی با سیستم تولیدمثلی چندهمسری (Polygamous) است (11). خرسها ازنظر بومشناختی دارای اهمیت زیادی هستند. این گوشتخواران بزرگ جثه در رأس هرم غذایی قرار داشته و به دلیل شبکه گسترده بومشناختی که با اجزای زنده و غیرزنده ایجاد کردهاند، بهعنوان یک گونه کلیدی و چتر برای اکوسیستم محسوب میشوند؛ ازاینرو اطلاعات توزیع و پراکنش، زیستشناسی و اکولوژی گونه برای اجرای برنامههای راهبردی و موقعیت حفاظتی گونه بسیار ضروری است (18).
مواد و روشها
محدوده موردمطالعه: استان کرمانشاه با 14 شهرستان (شکل 1) و مساحتی معادل 24361 کیلومترمربع در مختصات جغرافیایی بین 33 درجه و 36 دقیقه تا 35 درجه و 15 دقیقه عرض شمالی و 45 درجه و 24 دقیقه و 48 درجه و 30 دقیقه طول شرقی از نصفالنهار گرینـویچ در منتهیالیه غربی کشور قرار دارد. از سه طرف با استانهای کردستان، همدان، لرستان و ایلام و از سمت شرق نیز با کشور عراق همجوار است. استان کرمانشاه ناحیهای کوهستانی بوده که بین فلات ایران و جلگه بینالنهرین قرارگرفته و سراسر آن را قلل و ارتفاعات رشتهکوه زاگرس پوشانده است. متوسط میزان بارندگی در مناطق مختلف استان بین 300 تا 800 میلیمتر در نوسان است و متوسط بارش استان در درازمدت بیش از 537 میلیمتر تخمین زده شده است (3).
نقاط حضور: نقاط حضور گونه حاصل پیمایشهای میدانی سالهای 1395 تا 1398 است؛ که در بخشهای مختلف استان انجام گرفت.
شکل 1 - موقعیت استان کرمانشاه در کشور و نقاط حضور گونه در استان
این بازدیدهای میدانی با در نظر گرفتن تمامی لکههای پراکنش گونه در طیفهای مختلف زیستگاهی در سطح استان انجامشده است. به این منظور پس از شناسایی شدن حدود پراکنش، پیمایشهای میدانی از لکهها به عمل آمد و تمام نمایههای بهجامانده از گونه شامل ردپا، سرگین و محلهای استراحت بهوسیله سامانه موقعیتیاب مکانی برداشت شد. بیشترین مشاهده خرس در بخشهای جنوبی و غربی استان، مربوط به منطقه حفاظت شده قلاجه و ارتفاعات دالاهو بود. حد نهایی مناطق مشاهده نقاط حضور گونه، بخشهای شمال غربی استان در حاشیه رودخانه سیروان بود. شکل 1 موقعیت نقاط حضور ثبتشده برای اینگونه را در سطح استان کرمانشاه نمایش میدهد.
متغیرهای زیستگاهی: با بررسی پژوهشهای انجامگرفته بر روی گونه (23،15،7،5) متغیرهای مورداستفاده در فرآیند مدلسازی شناسایی شدند که در سه دسته کلی متغیرهای اکولوژیک، توپوگرافی و انسانی قابلطبقهبندی هستند. متغیر کاربری اراضی با استفاده از تصاویر سنجنده لندست 8 و پس از طی پیشپردازشهای لازم بر روی باندهای تصاویر با استفاده از روش طبقهبندی حداکثر احتمال، تهیه گردید. پس از تهیه نقشه کاربری اراضی، شاخص تکهتکه شدگی (Fragmentation) با استفاده از آنالیز الگو (Pattern Analysis) با فیلتر 3×3 تهیه شد. آمادهسازی متغیرهای ذکرشده در محیط نرمافزار ادریسی انجام گرفت. با استفاده از نقشه کاربری/پوشش، متغیر بی نظمی در پوشش زمین به وسیله متریک میانگین مساحت وزن دادهشده (Area weighted mean shape index) با استفاده از نرمافزار FRAGSTATS محاسبه گردید. پس از الحاق نتایج آن به جدول توصیفی، لایه کاربری /پوشش اراضی بر اساس فیلد اطلاعاتی حاصل از تحلیل مذکور، به رستر تبدیل و در تحلیل وارد شد. متغیر تراکم پوشش گیاهی با استفاده از تصاویر لندست 8 از مورخه 1 فوریه 2015 تا 1 فوریه 2018 محاسبه گردید. ازآنجاکه هدف، استفاده از میانگین شاخص تراکم پوشش گیاهی بود، از سامانه گوگلارثانجین استفاده شد (16). متغیرهای شیب و رطوبت سطح زمین (Compound Topographic Index) با استفاده از نقشه ارتفاع تهیه شدند. شاخص رطوبت سطح که با عنوان آمیختگی توپوگرافی نیز یاد میشود، رطوبت را برای سطح زمین و با در نظر گرفتن آبراههها محاسبه میکند. تمام متغیرهای مورداستفاده با ابعاد 100×100 متر در نرمافزار ArcGIS10.4 آماده شدند. پس از آمادهسازی تمامی متغیرها، میزان همبستگی بین آنها بررسی شد و آن دسته از متغیرهایی که دارای همبستگی بالای 75 درصد بودند وارد تحلیل نشدند. اجرای دستور همبستگی با استفاده از دستور محاسبات باندی در نرمافزار فوق انجام گرفت. جدول 1 تمام متغیرهای مورداستفاده را به تفکیک دامنه تغییر، واحد و منبع نمایش میدهد.
مدل توزیع و کمی سازی آشیان بومشناختی: بهمنظور مدلسازی انتشار گونه از الگوریتم آنتروپی بیشینه در نرمافزار MaxEnt استفاده شد. این نرمافزار با استفاده از رویکرد آنتروپی بیشینه، دادههای حضور را با متغیرهای محیطی مقایسه میکند (2). در این نرمافزار محاسبه اهمیت متغیرها در مدلسازی به روش تحلیل جکنایف و محاسبه اعتبار مدل با استفاده از مساحت سطح زیر منحنی (AUC) در تابع ROC انجام گرفت. 70 درصد داده های حضور برای آموزش و 30 درصد برای آزمون مورد استفاده قرار گرفتند. بهمنظور بررسی سهم هر یک از شهرستانهای استان از دربرگیری زیستگاه مطلوب از حد آستانه (true skill statistics) TSS استفاده گردید. ازاینرو نقاط حضور و نقاط پسزمینه مورد استفاده در اجرا استخراج شدند و از پکیج ROCR در نرمافزار R3.5.2 حد آستانه محاسبه گردید. کمی سازی پهنای آشیان اکولوژیک با استفاده از متریکهای موجود در نرمافزار ENMTools انجام گرفت.
جدول 1- متغیرهای زیستگاهی مورداستفاده در مدلسازی و دامنه تغییر هر متغیر
متغیر زیستگاهی |
دامنه تغییر |
واحد |
منبع |
شاخص رطوبت سطح زمین |
1/19-31/2 |
- |
ارتفاع |
ارتفاع |
3362-108 |
متر |
USGS |
فاصله از زمینهای کشاورزی دیم |
70013-0 |
متر |
نقشه کاربری اراضی |
فاصله از رودخانه |
13366-0 |
متر |
نقشه رودخانههای استان |
فاصله از جاده |
15068-0 |
متر |
نقشه جادههای اصلی استان |
فاصله از چشمه |
54831-0 |
متر |
موقعیت چشمه |
فاصله از مناطق مسکونی |
16105-0 |
متر |
موقعیت مناطق مسکونی |
میانگین شاخص تراکم پوشش گیاهی |
1-01/0 |
- |
گوگلارثانجین |
کاربری اراضی |
7 طبقه |
- |
لندست 8 |
شاخص تکهتکه شدگی سیمای سرزمین |
37/0-0 |
- |
کاربری اراضی |
متریک میانگین مساحت وزن دادهشده |
41/33-17/8 |
0 تا 100 |
کاربری اراضی |
در این راستا از متریکهای B1 و B2 برای بررسی پهنا استفاده شد. شاخص پهنا یا عرض آشیان اکولوژیک با دو شاخص Levin’s B1 و B2 (uncertainty) (13) بیان میشوند. مقادیر مربوط به این شاخصها نیز بین 0 تا 1 در نوسان است؛ که مقادیر نزدیک به 0 به معنای تخصص گرایی و مقادیر نزدیک به 1 به معنای تخصص عمومی است (25). بهمنظور اندازهگیری همپوشانی از 2 شاخص Schoener’s D (22) و Hellinger’s I (26) استفاده شد.
نتایج
مدلسازی: اعتبار سنجی مدل به کمک آماره AUC انجام گرفت. مقدار این آماره برای دادههای آزمون و آموزش در اجرای مدل به ترتیب برابر 90/0 و 87/0 محاسبه گردید؛ که نشان از تائید اجرای مدل است. نتایج آزمون بیونومینال نیز نشان داد که در تمام حد آستانههای موجود تفاوت معناداری با مدل تصادفی وجود دارد (Pvalue<0.001). بر اساس نتایج، متغیرهای کاربری اراضی، فاصله از زمینهای دیم و متریک AWMSI بیشترین تأثیر را در فرآیند مدلسازی زیستگاه داشتهاند. شکل 2 نمودار تغییرات مربوط به اثر متغیرهای محیطی را نمایش میدهد.
مطلوبیت زیستگاه و کمی سازی آشیان: شکل 3 نقشه احتمال حضور گونه یا مطلوبیت زیستگاه را نمایش میدهد. در این شکل بخشهای آبیرنگ نشاندهنده زیستگاههای مطلوب با احتمال حضور بالا و بخشهای قهوهایرنگ نیز بیانکننده زیستگاههای با احتمال حضور پایین است. بر این اساس بیشتر زیستگاههای مطلوب گونه در محدوده منطقه حفاظتشده قلاجه در جنوب استان و در بخشهای شمال غربی استان و در مجاورت کشور عراق قرار دارد. در ادامه شکل، موقعیت قرارگیری زیستگاه مطلوب گونه در شهرستانهای مختلف نمایش دادهشده است. بر اساس این نتایج بخشهای غربی استان در شهرستانهای گیلانغرب، قصر شیرین و سرپل ذهاب فاقد مطلوبیت هستند.
جدول 2 همپوشانی آشیان اکولوژیک خرس را بین مناطق شهرستانهای مختلف استان نمایش میدهد. مقادیر بالای این شاخص به معنای همپوشانی کامل و مقادیر پایین آن به معنای همپوشانی کم است.
شکل 2- اهمیت متغیرهای تأثیرگذار در مدلسازی زیستگاه به روش تحلیل جکنایف
شکل 3- مطلوبیت زیستگاه خرس قهوهای در استان کرمانشاه
بر اساس نتایج این شاخص، بیشترین میزان همپوشانی آشیان اکولوژیک بین شهرستان کرمانشاه و شهرستان هرسین است. کمترین میزان همپوشانی بین شهرستانهای استان نیز بین شهرستانهای کنگاور و گیلانغرب است؛ که نشان از تفاوت در پهنای آشیان گونه در دو شهرستان است.
جدول 2-همپوشانی آشیان اکولوژِیک خرس بر اساس شاخص D
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
1 |
1 |
5/0 |
5/0 |
39/0 |
55/0 |
5/0 |
5/0 |
5/0 |
43/0 |
5/0 |
5/0 |
5/0 |
49/0 |
5/0 |
55/0 |
2 |
× |
1 |
63/0 |
54/0 |
56/0 |
62/0 |
65/0 |
57/0 |
69/0 |
67/0 |
64/0 |
67/0 |
52/0 |
40/0 |
66/0 |
3 |
× |
× |
1 |
66/0 |
49/0 |
59/0 |
66/0 |
61/0 |
69/0 |
65/0 |
61/0 |
69/0 |
62/0 |
44/0 |
74/0 |
4 |
× |
× |
× |
1 |
79/0 |
72/0 |
83/0 |
55/0 |
72/0 |
63/0 |
77/0 |
65/0 |
66/0 |
73/0 |
67/0 |
5 |
× |
× |
× |
× |
1 |
30/0 |
5/0 |
27/0 |
56/0 |
49/0 |
28/0 |
54/0 |
55/0 |
30/0 |
66/0 |
6 |
× |
× |
× |
× |
× |
1 |
74/0 |
66/0 |
84/0 |
75/0 |
70/0 |
73/0 |
71/0 |
64/0 |
77/0 |
7 |
× |
× |
× |
× |
× |
× |
1 |
67/0 |
77/0 |
82/0 |
67/0 |
70/0 |
69/0 |
66/0 |
73/0 |
8 |
× |
× |
× |
× |
× |
× |
× |
1 |
80/0 |
72/0 |
69/0 |
71/0 |
58/0 |
60/0 |
73/0 |
9 |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
54/0 |
49/0 |
67/0 |
57/0 |
35/0 |
64/0 |
10 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
77/0 |
71/0 |
61/0 |
67/0 |
67/0 |
11 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
77/0 |
68/0 |
58/0 |
81/0 |
12 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
63/0 |
56/0 |
67/0 |
13 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
47/0 |
61/0 |
14 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
75/0 |
15 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
* شمارههای ذکرشده مطابق با ردیف و نام شهرستانها در شکل 3 است.
جدول 3 میزان همپوشانی آشیان اکولوژیک خرس در شهرستانهای مختلف استان را نمایش میدهد. مقادیر بالای این شاخص نیز مانند شاخص قبل میزان بالای همپوشانی را نمایش میدهد. بر این اساس بالاترین میزان همپوشانی آشیان اکولوژیک در خصوص زیستگاه گونه مربوط به شهرستان هرسین و شهرستان کرمانشاه و کمترین آن مربوط به گیلانغرب و شهرستان کنگاور است.
حد آستانه و توزیع مطلوبیت: بر اساس حد آستانه TSS مقدار حد آستانه برابر 36/0 شناسایی گردید که در این نقطه برش به ترتیب حساسیت برابر 77 و ویژگی برابر 87 محاسبه شدند. شکل 4 نقطه برش را نمایش میدهد. پس از شناسایی حد آستانه، مقدار مذکور بر روی نقشه پیوسته مطلوبیت زیستگاه اعمال گردید و نقشه گسسته (باینری) مطلوبیت زیستگاه محاسبه شد.
شکل 5 موقعیت زیستگاه مطلوب گونه را به ازای شهرستانهای استان نمایش میدهد. بخشهای قرمزرنگ، زیستگاه نامطلوب و بخشهای آبیرنگ، زیستگاه مطلوب را نمایش میدهند. سهم هر شهرستان از دربرگیری زیستگاه مطلوب گونه در جدول 4 ذکرشده است.
بر اساس نتایج، پهنای آشیان اکولوژیک گونه در محدوده استان شرایط متوسطی را نمایش میدهد.
جدول 3-همپوشانی آشیان اکولوژِیک خرس بر اساس شاخص I
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
1 |
1 |
5/0 |
5/0 |
39/0 |
62/0 |
5/0 |
5/0 |
5/0 |
51/0 |
5/0 |
5/0 |
51/0 |
5/0 |
5/0 |
56/0 |
2 |
× |
1 |
85/0 |
74/0 |
77/0 |
81/0 |
85/0 |
76/0 |
85/0 |
86/0 |
82/0 |
86/0 |
75/0 |
60/0 |
86/0 |
3 |
× |
× |
1 |
85/0 |
71/0 |
75/0 |
82/0 |
80/0 |
85/0 |
85/0 |
77/0 |
88/0 |
81/0 |
61/0 |
89/0 |
4 |
× |
× |
× |
1 |
86/0 |
85/0 |
91/0 |
75/0 |
82/0 |
83/0 |
89/0 |
81/0 |
80/0 |
84/0 |
82/0 |
5 |
× |
× |
× |
× |
1 |
44/0 |
66/0 |
38/0 |
74/0 |
65/0 |
39/0 |
72/0 |
73/0 |
42/0 |
83/0 |
6 |
× |
× |
× |
× |
× |
1 |
90/0 |
85/0 |
93/0 |
90/0 |
89/0 |
88/0 |
84/0 |
79/0 |
90/0 |
7 |
× |
× |
× |
× |
× |
× |
1 |
77/0 |
87/0 |
92/0 |
83/0 |
87/0 |
85/0 |
73/0 |
88/0 |
8 |
× |
× |
× |
× |
× |
× |
× |
1 |
91/0 |
89/0 |
87/0 |
89/0 |
78/0 |
77/0 |
89/0 |
9 |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
64/0 |
59/0 |
80/0 |
72/0 |
45/0 |
84/0 |
10 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
90/0 |
89/0 |
82/0 |
80/0 |
87/0 |
11 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
91/0 |
83/0 |
74/0 |
92/0 |
12 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
84/0 |
64/0 |
86/0 |
13 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
56/0 |
82/0 |
14 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
87/0 |
15 |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
× |
1 |
* شمارههای ذکرشده مطابق با ردیف و نام شهرستانها در شکل 3 است.
شکل 4- حد آستانه TSS بر روی منحنی ROC
از میان شهرستانهای استان بیشترین پهنای آشیان اکولوژیک در محدوده شهرستان پاوه قرار دارد که نشان از وجود شرایط مساعد برای حضور اینگونه است. علاوه بر این شهرستان، شهرستانهای دالاهو و هرسین نیز دارای پهنای آشیان مناسبی برای گونه موردمطالعه هستند. درصد نسبی زیستگاه در سطح شهرستان نشان داد که 68/46 درصد شهرستان پاوه برای گونه دارای مطلوبیت است که بیشترین مقدار با توجه به وسعت شهرستان محسوب میشود. بعد از شهرستان پاوه، شهرستانهای گیلانغرب و دالاهو نیز وسعت بالایی از زیستگاه گونه را شامل میشوند. درمجموع از میان شهرستانهای استان، شهرستانهای کرمانشاه، گیلانغرب، دالاهو و اسلامآباد بیشترین سهم از زیستگاه گونه در سطح استان را دارا میباشند.
شکل 5- توزیع لکههای زیستگاهی در شهرستانهای مختلف استان کرمانشاه
جدول 4- متریکهای پهنای آشیان اکولوژیک، وسعت زیستگاه (به هکتار) و درصد فراوانی کل
ردیف |
نام محدوده |
B1 |
B2 |
وسعت زیستگاه مطلوب |
درصد از شهرستان |
درصد کل استان |
1 |
استان کرمانشاه |
42/0 |
95/0 |
57/332018 |
- |
100 |
2 |
دالاهو |
57/0 |
97/0 |
80/47864 |
30/25 |
41/14 |
3 |
اسلامآباد |
52/0 |
96/0 |
34//40826 |
94/19 |
29/12 |
4 |
قصر شیرین |
02/0 |
81/0 |
50/338 |
17/0 |
10/0 |
5 |
گیلانغرب |
39/0 |
93/0 |
14/51768 |
25/23 |
59/15 |
6 |
هرسین |
55/0 |
96/0 |
05/6751 |
72/6 |
03/2 |
7 |
جوانرود |
60/0 |
97/0 |
27/17067 |
45/22 |
14/5 |
8 |
کنگاور |
50/0 |
96/0 |
26/5793 |
02/13 |
74/1 |
9 |
کرمانشاه |
46/0 |
96/0 |
27/60423 |
49/10 |
19/18 |
10 |
پاوه |
69/0 |
97/0 |
4037161 |
68/46 |
19/11 |
11 |
روانسر |
51/0 |
96/0 |
84/18677 |
61/15 |
62/5 |
12 |
صحنه |
53/0 |
96/0 |
28/7297 |
97/4 |
19/2 |
13 |
ثلاث باباجانی |
40/0 |
94/0 |
19/23363 |
14/14 |
06/7 |
14 |
سرپل ذهاب |
20/0 |
90/0 |
63/1730 |
91/1 |
52/0 |
15 |
سنقر |
50/0 |
96/0 |
19/7342 |
16/3 |
21/2 |
کارایی مناطق حفاظتشده: شکل 6 گستره مطلوبیت زیستگاه را پس از اعمال حد آستانه در سطح استان کرمانشاه مشخص میکند. در این تصویر بخشهای آبیرنگ، زیستگاه مطلوب و مناطق مشکیرنگ، زیستگاههای نامطلوب گونه را نمایش میدهد. مرز مناطق حفاظتشده نیز در این تصویر بر روی لکههای زیستگاهی قرارگرفته است. منطقه شکارممنوع نواکوه، منطقه حفاظتشده قلاجه و بوزین و مرخیل و مجموعه بیستون (منطقه حفاظتشده و پناهگاه حیاتوحش) زیستگاه اینگونه را پوشش دادهاند، اما هنوز بخشهای وسیعی از زیستگاه گونه در خارج از مرزهای مناطق حفاظتشده قرارگرفته است. این خلأهای پوشش در بخشهایی از زیستگاه گونه از منطقه نواکوه تا بوزین و مرخیل، اطراف منطقه حفاظتشده قلاجه به سمت حمیل و هلشی قابلمشاهده است. نتایج نشان میدهد که مساحت زیستگاه مطلوب گونه در خارج از مرز مناطق حفاظتشده استان به 84/265842 هکتار میرسد که مساحت قابلتوجهی است. مناطق شکارممنوع قراویز و زله زرد عملاً پوششی از زیستگاه گونه موردمطالعه ندارند. نتایج توزیع لکههای زیستگاهی نشان داد که پناهگاه حیاتوحش بیستون 64/7278 هکتار، منطقه حفاظتشده بیستون 05/8323 هکتار، امروله و دالاخانی 89/2266 هکتار، منطقه حفاظتشده قلاجه 72/25717 هکتار، منطقه شکارممنوع نواکوه 60/9654 هکتار، منطقه حفاظتشده بدر و پریشان 68/613 هکتار، منطقه حفاظتشده بوزین و مرخیل 46/12319 هکتار و منطقه شکارممنوع هشیلان 64/1 هکتار از وسعت زیستگاههای مطلوب را در برگرفتهاند.
شکل 6- نقشه باینری مطلوبیت زیستگاه خرس در محدوده استان کرمانشاه و مرز مناطق حفاظتشده (1= پناهگاه حیاتوحش بیستون، 2= منطقه حفاظتشده بیستون، 3= منطقه شکارممنوع امروله و دالاخانی، 4= منطقه حفاظتشده قلاجه، 5= منطقه شکارممنوع نواکوه، 6= منطقه شکارممنوع قراویز، 7= منطقه شکارممنوع زله زرد، 8= منطقه حفاظتشده بدر و پریشان، 9= منطقه حفاظتشده بوزین و مرخیل، 10= منطقه شکارممنوع هشیلان)
بحث و نتیجهگیری
نتایج این نوشتار نشان میدهد که توزیع گونه به تفکیک شهرستان و مناطق حفاظت شده استان کرمانشاه به چه شکلی است و هر یک از شهرستانها چه سهمی از پوشش زیستگاه بالقوه گونه دارند؛ عملی که شاید بتواند در آینده استراتژیهای هر شهرستان را تحت تأثیر قرار دهد. در میان شهرستانهای استان، شهرستان پاوه به نسبت، زیستگاه بالقوه بیشتری را برای گونه فراهم کرده و نزدیک به نیمی از وسعت این شهرستان برای حضور گونه شرایط مساعدی دارد. از سوی دیگر، پهنای آشیان اکولوژیک خرس در این شهرستان بالاترین مقدار است (جدول 4). شهرستانهای دالاهو، گیلانغرب و جوانرود نیز به ترتیب در ردههای بعدی ازنظر مساحت زیستگاه در شهرستان قرار دارند؛ بنابراین برنامههای راهبردی برای حفاظت ازاینگونه ارزشمند میتواند با تأکید بر این شهرستانها انجام گیرد.
بر اساس نتایج، متریکهای موردبررسی پهنای آشیان اکولوژیک در شهرستانهای پاوه، جوانرود و کرمانشاه به یکدیگر شباهت دارند که نشان از وجود شرایط مشابه برای این شهرستانها است. بهغیراز سه شهرستان ذکرشده پهنای آشیان گونه در شهرستانهای دالاهو و اسلامآباد غرب نیز مشابه است. وسعت بالای زیستگاه بالقوه (57/332018 هکتار) و عدم پوشش کامل آن توسط مناطق حفاظتشده زمینهساز ایجاد بروز تعارضهایی در خارج از مرز این مناطق شده است. این تعارضها در بخشهای جنوبی استان در محدوده شهرستانهای اسلامآباد غرب، گیلانغرب و دالاهو بهمراتب بیش از سایر بخشهای استان است. عمده تعارضهایی که در این مناطق حاصلشده، ناشی از حضور گونه در زمینهای کشاورزی و باغات بوده و بعضاً صدمات جبرانناپذیری به همراه داشته و حتی در برخی از ارتفاعات قلاجه حمله به چوپانان نیز گزارششده است. در این مطالعه، کاربری اراضی و فاصله از زمینهای دیم بهعنوان مهمترین متغیرهای تأثیرگذار در دامنه انتشار گونه برشمرده شدند. زمینهای دیم منطقه با کشتهای متفاوت تأمینکننده نیازهای غذایی خرس میباشند و شدت حضور اینگونه در زمینهای کشاورزی در فصول برداشت بیشتر است. بیشترین حضور در بخشهای کشاورزی مربوط به شهرستان اسلامآباد غرب، دالاهو و بخش گهواره است. در میان کاربری/پوششهای مختلف، جنگلها بهمراتب نسبت به سایر کاربریها دارای اهمیت بالاتری بودند؛ که مشابه یافته چواردس و همکاران (2013) در مطالعه بر روی اینگونه در کوهستانهای پیرین یونان است (9). در مطالعهای که توسط کانلوپولوس و همکاران (2006) بر روی مطلوبیت زیستگاه خرس در شمال Pindos یونان انجام گرفت نیز نتایج مشابهی حاصل شد که مشابه یافتههای این مطالعه است (12). نفوذ بخشهای روستایی مانند جوامع حاشیه مناطقِ حفاظتشده بیستون و یا جوامع روستایی بخش گهواره و قلاجه به محدوده جنگلهای زاگرس یکی از عوامل مهم ایجاد درگیری و تعارض است.
بیشترین حضور خرس با ماهیت تعارضی در داخل باغهای حاشیه مناطق حفاظتشده و همچنین زمینهای کشاورزی است که پیرامون جنگلهای انبوه بلوط قرارگرفتهاند؛ که همسو با نتایج زایوکفسکا و همکاران (2016) در Carpathians است و نشان داد خرس قهوهای مناطقی با پوشش جنگلی بالا، نزدیک به لبههای جنگل و با فشار انسانی کم را برمیگزیند (27). در مطالعهای که توسط گاستینو و همکاران (2019) در منطقه پرینه فرانسه انجام گرفت نیز نتایج نشان داد که احتمال حمله خرس به دامها در مجاورت مناطق جنگلی بیشتر بوده (10) که مشابه یافتههای این پژوهش است.
متریک میانگین مساحت وزن دادهشده بعد از دو متغیر ذکرشده بیشترین تأثیر را در این مطالعه داشت. این متریک مقدار بینظمی موجود در طبقات مختلف کاربری/پوشش را نمایش میدهد. مقادیر بالای این شاخص نشان از تمایل حضور گونه به مناطقی نامنظم ازنظر کاربری/پوشش را دارد. مقدار مطلوبیت این متریک نشان داد که با افزایش مقدار بینظمی بر مطلوبیت زیستگاه افزوده میشود اما با افزایش مقدار این شاخص از 20 به بالا مقدار مطلوبیت کاسته میشود که نشان از عدم تمایل برای حضور گونه در مناطق بینظمی بالا ازنظر پوشش زمین است؛ بنابراین نتایج این متریک نشان میدهد خرسها با حضور در مناطقی با حدود متوسطی از بینظمی در پوشش زمین تمایل دارند که میتواند تائید کننده تمایل خرس به حضور در جنگل و مناطق حاشیه آن باشد. معیارهای ارتفاع و شاخص تراکم پوشش گیاهی در ردههای بعدی قرار دارند. معیار ارتفاع نشان داد که با افزایش ارتفاع تا حدود 1600 متری بر مطلوبیت زیستگاه گونه افزوده میشود و از این ارتفاع به بعد از مطلوبیت کاسته خواهد شد؛ اما در مطالعهای که توسط سیول (2019) در آنتالیا کشور ترکیه انجام گرفت، ارتفاع مطلوب حضور گونه برابر با 2000 متر ذکرشده (24) که متفاوت از یافتههای این مطالعه است. ارتفاعات بالاتر و کاهش مطلوبیت آن برای گونه موردمطالعه میتواند به دلیل عدم وجود پوششهای جنگلی باشد. در مطالعهای که توسط پسیلکو و همکاران (2004) در کوهستانهای آپنینی ایتالیا انجام گرفت نتایج نشان داد که خرسها به حضور در ارتفاعات بین 800 تا 1000 متری تمایل دارند که این امر میتواند به دلیل ساختار زمین و وجود جنگل باشد (19).
شاخص تراکم پوشش گیاهی نیز نشان داد با افزایش مقدار شاخص تا حدود 6/0 بر مطلوبیت افزوده و از این مقدار به بعد از مطلوبیت زیستگاه کاسته میشود؛ این مقدار تراکم، مشخصکننده مناطق جنگلی است. تمایل خرس به حضور در مناطق جنگلی در سایر مطالعات مربوط به اینگونه نیز ذکرشده است (7، 8). در مطالعه کوچالی و همکاران (1397) در دامنه شمال البرز نیز افزایش مقدار پوشش گیاهی تا رسیدن به مرز جنگل باعث افزایش مطلوبیت زیستگاه گونه میشود (5) که مشابه یافتههای این پژوهش است. ازآنجاییکه مقادیر بالای این شاخص میتواند بیانگر مناطقی با کاربری کشاورزی باشد، بنابراین به نظر میرسد که خرس قهوهای در استان کرمانشاه مناطقی با پوشش جنگلی متراکم را ترجیح میدهد و عمده تعارض رخداده برای گونه نیز در حاشیه این مناطق با کاربریهای باغ و کشاورزی به وجود میآید.
اما میتوان برخی از اقدامات را در راستای کاهش تعارض پیش گرفت. ازآنجاییکه عمده تعارضهای رخداده در محدوده باغات و زمینهای کشاورزی بوده، میتوان از فنسها برای جلوگیری از ورود به این مناطق جلوگیری کرد. ازجملهی این مناطق میتوان به روستاهای شهر گهواره (سرخک سنجابی، دارخور، چنار، نابی لک)، اسلامآباد غرب (روستاهای ترازک، انجیرک، دار بادام، توراب) و کرند غرب (سرخه دیزه، هریر) اشاره کرد. در کنار این عمل، فنس کشی را میتوان برای دامهای اهلی که در فضای بیرون نگهداری میشوند نیز پیاده کرد. آگاهسازی جوامع محلی با نیازهای خرس در فصول مختلف و ویژگیهای رفتاری آن نیز میتواند در راستای تداوم حضور و بقا مفید باشد. اتصال مناطق نوا کوه و قلاجه در جنوب استان میتواند بهعنوان یک رویکرد مدیریتی مناسب برای حفظ یکپارچگی زیستگاه اینگونه در نظر قرار گیرد، این امر میتواند با وضع قوانینی خاص در محدوده تردد وسایل نقلیه راهگشا باشد.