مدل‌سازی مبتنی بر احتمال برای تحلیل کمی ریسک آتش‌سوزی در مناطق حفاظت شده استان گیلان

آتش‌سوزی‌ یکی از آشفتگی‌های اصلی بوم‌سازگان‌های جنگلی است که آثار تخریبی زیادی روی پوشش گیاهی، خاک و آب داشته است. در تحقیق حاضر با استفاده از دادۀ تاریخی آتش‌سوزی (2022-1992) در 20 منطقه حفاظت شده (256488 هکتار) در استان گیلان، طبقه‌بندی‌های اندازه و تکرار آتش‌سوزی با استفاده از مدل‌های فرایند نقطه‌ای (ppms)، مدل‌سازی‌های احتمال احتراق و احتمال سوختن مبتنی بر شبیه‌ساز flammap انجام شد. نتایج توزیع کلاس‌های اندازۀ آتش‌سوزی نشان داد که آتش‌سوزی‌های کمتر از 10 هکتار به‌ترتیب 62.4 و 30.8 درصد از تعداد آتش‌سوزی‌ها و مساحت سوخته‌شده را به‌خود اختصاص داده‌اند. آتش‌سوزی‌های بین 10 تا 50 هکتار به‌ترتیب 32.8 و 46.2 درصد از تعداد آتش‌سوزی‌ها و مساحت سوخته‌شده را شامل شده‌اند. آتش‌سوزی‌های 100-50 هکتار به‌ترتیب 2.7 و 14.5 درصد از تعداد آتش‌سوزی‌ها و مساحت سوخته‌شده را پوشش داده‌اند. در نهایت، یک آتش‌سوزی بیش از 100 هکتار کمتر از 0.5 درصد از تعداد آتش‌سوزی‌ها را شامل شده است، اما به تنهایی 8.4 درصد از مساحت سوخته‌شده را به‌خود اختصاص داده است. توزیع تکرار آتش‌سوزی در منطقه مورد مطالعه از 0 تا 6 متغیر است. 40 درصد از مناطق حفاظت شده هیچ آتش‌سوزی را تجربه نکرده است. 13 درصد از این مناطق 2-1 تکرار آتش‌سوزی داشته‌ند. هم‌چنین، 48 درصد از این مناطق بیش از 2 تکرار آتش‌سوزی داشته‌ند. از نظر احتمال احتراق، 35 درصد از منطقه مورد مطالعه دارای مقادیر خیلی کم و کم، 36 درصد دارای مقادیر متوسط، 18 و 11 درصد به‌ترتیب دارای مقادیر زیاد و خیلی زیاد بودند. از نظر احتمال سوختن، 88 درصد از منطقه مورد مطالعه دارای ارزش‌های کم تا متوسط و 12 درصد دارای مقادیر زیاد و بسیار زیاد بودند. به‌طور کلی، بیشترین رخداد آتش‌سوزی در جنوب و مرکز منطقه مورد مطالعه که پوشیده از جنگل‌های پهن‌برگ متراکم با انباشت زیاد بار ماده‌سوختنی آماده اشتعال در فصل آتش‌سوزی (از اواخر خرداد تا اسفند) است، مشاهده شد که نیازمند امکان‌سنجی برای توسعۀ گزینه‌های مدیریت آتش‌سوزی است. تحقیق حاضر بر ضرورت کاربرد نتایج مدل‌سازی ریسک آتش‌سوزی در پیش‌گیری و مدیریت آتش‌سوزی با هدف حفظ خدمات بوم‌شناختی حاصل از پوشش گیاهی، آب و خاک در استان گیلان و سایر مناطق با شرایط مشابه، تاکید می‌کند.

probability-based modeling for a quantitative wildfire risk analysis in the protected areas of guilan province

introductionwildfires reveal evidence of forest soil, water, and vegetation disturbances resulting from various interacting natural and human factors that create patterns that vary spatially and temporally. fire risk assessment allows for identifying these factors and estimating their area of influence, thereby determining locations at high fire risk. fire risk assessment typically involves the ignition probability (ip) and burn probability (bp) modeling of natural and human-made resources, as well as identifying resource responses to fires of varying severity. in recent decades, wildfires have caused significant damage to the hyrcanian forests in northern iran, even in the protected areas. therefore, this study focuses on the spatial distributions of fire size, fire frequency, ip, and bp as essential components of the fire risk framework. first, a historical fire database (including ignition points, burned area, etc.) was prepared using available resources and field surveying. second, a modeling approach using a limited number of auxiliary variables representing the fire environment (fuel, topography, and weather) and the historical fires (1992-2022) was implemented to calculate ip and bp. the spatial distribution of these parameters helps improve decision-making in fire prevention and control strategies.materials and methodsguilan province is located in northern iran and has an area of 14,044 km2 with an average elevation of 741 m above sea level. this province has 20 natural protected areas, which cover a total of 256,488 ha. arcgis 10.8 was used to create a historical fire database in the study area by digitizing fires between 1992 and 2022 from maps or by importing information directly from the previous gis datasets. point process models (ppms) were used to analyze the spatial distribution of fire frequency. ppms are a regression approach to model point data (i.e., geographic coordinates) for the number of times a 100-m pixel burns between 1992 and 2022. an existing raster map of the study area was converted to points by calculating the center of each pixel, and each point was assigned a frequency. furthermore, ip was calculated as the average ignition probability occurring over a year in a raster pixel. to help use the fire ignition density to plan preventive activities, the output values were classified into five classes reflecting ignition occurrences (from very low to very high). finally, the fire risk using bp was assessed by considering topography, fuel loads, and weather using flammap. to calculate bp, 1000 random ignition points were created based on the distribution of historical ignition points in the study area. the maximum fire simulation time was set to 6 hr (the average fire duration in the area).results and discussion there are 186 recorded fires (total burned area of 2232 ha) with an average annual number of 6 fires (average burned area of 12 ha) in the protected areas. fires <10 ha accounted for 62.4% and 30.8% of the fire number and the burned area, respectively. fires (10-50 ha) accounted for 32.8% and 46.2% of the fire number and the burned area, respectively. fires (50-100 ha) accounted for 2.7% and 14.5% of the fire number and the burned area, respectively. finally, fires >100 ha accounted for <0.5% of the fire number but alone accounted for 8.4% of the burned area. the distribution of fire frequency ranged from 0 to 6. the largest protected areas (40%) experienced no fires. 13% of these areas had 1-2 fire frequencies. furthermore, 48% of this area had more than two fire frequencies. about 35% of the study area had very low and low ip values, 36% had medium ip, and 18 and 11% had high and very high ip, respectively. 88% of the study area had low and moderate bp values, and 12% had high and very high predicted values. two fire regimes can be distinguished in the area, one with relatively high fire frequency and bp (mainly at higher elevations) and the other with low fire frequency and bp (at lower elevations). high fire frequency and bp is very limited in extent and occurs in the patches in the southern area. in contrast, low fire frequency and bp regime is the most widespread regime in the area (except for the southern part). conclusionaccording to the simulated patterns of fire frequency, ip, and bp in the study area, a clear distinction between the actual historical fire perimeters and the predicted burn pattern is that there are areas of moderate to high ip where fires have not occurred in the past 30 years. this is particularly evident in the southern and central parts of the area, where fires have either not occurred or have been very limited in extent. therefore, a justifiable assessment could be that the likelihood of fire spread and vegetation communities undergoing extensive and long-term changes following the fire is high shortly. although this study focuses on the protected forest areas, this approach can be applied to fire risk modeling at larger scales. this allows for broader application in natural resources management and planning at regional and national levels. it also provides a comprehensive tool for assessing and managing forest vegetation, soil, and water vulnerability.