Next Article in Journal
Using Geographically Weighted Regression to Study the Seasonal Influence of Potential Risk Factors on the Incidence of HFMD on the Chinese Mainland
Next Article in Special Issue
GIS-Based Spatial and Multi-Criteria Assessment of Riverine Flood Potential: A Case Study of the Nitra River Basin, Slovakia
Previous Article in Journal
Multi-Scale Coal Fire Detection Based on an Improved Active Contour Model from Landsat-8 Satellite and UAV Images
Previous Article in Special Issue
Time-Series Clustering for Home Dwell Time during COVID-19: What Can We Learn from It?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forest Fire Hazards Vulnerability and Risk Assessment in Sirmaur District Forest of Himachal Pradesh (India): A Geospatial Approach

1
Centre for Climate Change and Water Research, Suresh Gyan Vihar University, Jaipur 302017, India
2
Faculty of Geotechnical Engineering, University of Zagreb, Hallerova aleja 7, 42000 Varaždin, Croatia
3
Department of Civil Engineering, University North, 42000 Varaždin, Croatia
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2021, 10(7), 447; https://doi.org/10.3390/ijgi10070447
Submission received: 23 April 2021 / Revised: 3 June 2021 / Accepted: 25 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue GIScience for Risk Management in Big Data Era)

Abstract

:
The Himachal Pradesh district’s biggest natural disaster is the forest fire. Forest fire threat evaluation, model construction, and forest management using geographic information system techniques will be important in this proposed report. A simulation was conducted to evaluate the driving forces of fires and their movement, and a hybrid strategy for wildfire control and geostatistics was developed to evaluate the impact on forests. The various methods we included herein are those based on information, such as knowledge-based AHP-crisp for figuring out forest-fire risk, using such variables as forest type, topography, land-use and land cover, geology, geomorphology, settlement, drainage, and road. The models for forest-fire ignition, progression, and action are built on various spatial scales, which are three-dimensional layers. To create a forest fire risk model using three different methods, a study was made to find out how much could be lost in a certain amount of time using three samples. Precedent fire mapping validation was used to produce the risk maps, and ground truths were used to verify them. The accuracy was highest in the form of using “knowledge base” methods, and the predictive value was lowest in the use of an analytic hierarchy process or AHP (crisp). Half of the area, about 53.92%, was in the low-risk to no-risk zones. Very-high- to high-risk zones cover about 24.66% of the area of the Sirmaur district. The middle to northwest regions are in very-high- to high-risk zones for forest fires. These effects have been studied for forest fire suppression and management. Management, planning, and abatement steps for the future were offered as suitable solutions.

1. Introduction

India is one of the world’s mega-rich areas in biodiversity, both in terms of fauna and flora. The incredible richness of India’s biodiversity includes 6 million square kilometers of forest cover. The forests of the country are both environmentally and economically valuable. In India, there are only around 1.7 million hectares of well-defined conifer forest, consisting of important species such as fir, teak, sal, and chir pine. The estimated timber growth in these forests is close to 2 billion cubic meters, which is worth around 60 million rupees (100 million USD) [1]. Large forest fires impact the climate, human health, and property, as well as resulting in danger to life. Global attention to forest fires has increased recently due to their major long-term threat to forest habitats and public safety, as well as shorter-term risks to property and human lives [2]. Fire plays an important role in ecological processes, including altering the composition of plant populations, conserving water, enhancing soil quality, and promoting biodiversity. Due to their role in natural plant succession, forest fires are an important mechanism for initiating new growth in the ecosystem. A region equivalent to 6 million km2 of land in India has been destroyed over the past two centuries due to forest fires [3]. Throughout that time, the Forest Survey of India (FSI) has been carrying out fieldwork in the country’s different forest plots since 1965 in order to record fires. For statistical significance, long-term observations of forest conditions were analyzed, and the results were released in [1]. A total of 95% of the country’s forest fires are of human origin, according to FSI research. As per data from the Indian Forest Service, 50% of the forests are likely to be vulnerable to burning. A few national and foreign NGOs (non-governmental organizations), including the Survey of India and state agencies, are also conducting research on forest fires. “India has carried out several forest fire incident case studies in various locations and this booklet summarizes it all their findings” [4].
Despite all of these interventions, the forest-fire-danger data bank remains inadequate. New methods will be needed to circumvent this problem. Satellite data can be very valuable in terms of new fire management strategies. Some new innovations, such as the use of satellites, should be applied in conjunction with the country’s extensive field-based fire data-collection programmed. According to analysts, the research data and facts that have been published show that the country is in a very serious crisis that requires immediate action to be resolved. In India, there are no details on the area and value of the forests burnt, as well as the amount of forest devastated by fire, which makes it impossible to determine the losses. They are incomplete because the number and extent of fires is not available. The explanation given for this is the fear of having to bear responsibility. There is thought to be about 1 million hectares of forest loss due to fire annually. According to an assessment conducted by the Indian Forest Service, each year, there are 3.73 million hectares of wildfires in India that have an effect on forests [5]. There are very few fire outbreaks in India that are caused by natural factors. There is widespread (99%) agreement among observers that the majority these fires are of the people’s own doing, and there is a connection to their socioeconomic status. The chief sources of forest fires in India are grazing, shifting cultivation, and the collection of non-wood products.
Forest-fire risk modeling uses remote sensing, and GIS, in 1990 [6], used an advanced very high-resolution radiometer to detect the fire risk for Rajaji National Park in Uttar Pradesh in India, which is considered to have a more than 50% higher fire risk. This paper also presents one of the earliest papers on this topic for areas in India. The authors of [7] used remote sensing to map vegetation types and using slope, proximity to settlements, and distance from roads to predict fire risk in Madhya Pradesh in India. The authors of [8] used neural networks, knowledge-intensive systems, in order to predict forest fires based on temperature, humidity, rainfall, and fire history. Similar to [8], the author in [9] used fuzzy sets and semi-triangular membership function in order to perform long-term prediction of forest-fire risk based on fire history and drought cycles. The authors in [10,11] also considered different areas for fire risk prediction. For areas in India, multiple authors considered estimating forest-fire risks [12,13,14,15]. The mentioned authors provide an overview of forest fire risk for different areas in India. The authors of [16] used multicriteria decision analysis and an analytical hierarchy process to estimate Bhajji forest fire risk in Himachal Pradesh in India. This paper presents only one forest part of Himachal Pradesh, while our paper presents the whole area of the Sirmaur district.

2. Study Area

The Sirmaur district is in the outer Himalayas and is also commonly known as the Shivalik range. The district shown in Figure 1 and is located between the latitudes of 31°00′ North and 77°00′ South, and 77°00′77′ longitude, between the base latitudes of longitude and the Surveys of India sheet ‘A’ and ‘L’ (Figure 1). It is located at an altitude ranging from 300 to 3647 m, with four major rivers. The forested area is around 48,682 hectares [17].
It is connected to the capital Shimla to the west and to the east. In the south is the state of Haryana and in the east is the Indian state of Uttar Pradesh.
The Monsoon that generally starts in the Sirmaur district is from the third week of June and lasts up to the second week of September. April and May are also the rainiest months of the spring season in the Sirmaur district. In July and August, rainfall is more active. The early rainy season provides 80% of the rainfall and the average rainfall of the Sirmaur district is about 1324.3 mm. Figure 2 shows the average rainfall of the district from temporal variation between 2001 and 2020.

3. Data and Methodology

The following resources and equipment are needed to complete the analysis. Below are the following various data used for the forest fire analysis and methodology with their flow chart for analysis. The land cover layer was created by combining digital image classification and visual interpretation techniques to view satellite images. As supplementary data, the region’s current land use/cover map was used. The key remotely sensed data source in this study was Sentinel 2 imagery from 20 May 2019, resized to a spatial resolution of 10 m [18]. Forest and wild land fire risks, including elevation, slope and aspect, population density and precipitation, were incorporated into ArcGIS for further analysis. The reasons for fires range greatly. The full data used for forest fire risk incorporates slope, elevation, aspect, land use and land cover, geology, geomorphology, buffer zone of road, drainage buffer, built up area buffer, and forest type. Slope, elevation, and aspect data were calculated from the Cartosat-1 digital elevation model [19]. Based on [20], with the increase of elevation, the probability of fire decreases and aspect is only significant on the warmer sides. Land use and land cover is defined with Sentinel-2 imagery [21]. Geology and geomorphology data were downloaded from Geological survey of India [22]. Geology and geomorphology data do not directly affect the risks of forest fires but wildfires can cause problems such as landslides and erosions [23,24]. The roads were downloaded from Geological Survey of India and the road and built-up area proximity of the 200-m buffer zone can significantly increase forest fire risk [25,26]. Therefore, a 200 m buffer zone is included around roads and built-up areas.
The moderate resolution imaging spectroradiometer (MODIS) has been widely used for wildfire incidence and distribution detecting and fire risk assessments. To check the data distribution and accessibility under the Universal Transverse Mercator (UTM) projection, 1:50,000 topographic paper maps with the MODIS hotspot data were mapped to the 1:50,000 topographic paper maps of the WGS84 datum. The comparatively higher temperature on the day of the fire was due to the burned and ash-covered pixels nearby, releasing the latent heat more quickly than usual, resulting in a higher temperature for ID11.
The full workflow is presented in Figure 3.

3.1. Analytical Hierarchy Process (AHP)

As a multi-criteria approach, the analytical hierarchy process (AHP) is one of the best for defining forest fire locations. Many studies have explored the proof of AHP’s claimed forest fire effectiveness [27]. In AHP strategy, the enormous choices are made comprehensible by knowledge from math and experts. We have used the nine-point non-stop scale to review the two parameters. Both odd numbers (1, 3, 5, and 7) and even numbers (2, 4, and 6) refer to equal and humble characteristics, but the odd ones stand out (Table 1). The AHP is a methodology focused on objective objectives branching into multiple criteria (objective branches) that is used for multiple-criteria problems. Table 2 shows the random consistency index (RI).

3.2. Deriving the Weights Using AHP/ANP

Analytic network processing (ANP) is a way of evaluating different parameters by prioritizing and compiling them [29]. Each layer is multilayered, indicating that the relationships between these interdependent groups are too complicated to understand without some effort. This results in their relationship to the different classes being calculated by ANP, so the relationship among these nine has been figured out thusly. The thematic layer weights using ANP and AHP can be derived by following these steps:
  • Construction of a model: many models have been found for the forest reserve map ability, based on a literature review. It is imperative that the problem is identified at both an abstract and thematic level before putting it into model layers.
  • Generation of pairwise comparison matrices: with this handy table—aside from arbitrary designations of importance as found in arbitrary scales—the relative values are assigned as follows: 1 means the same importance of the two themes, while a score of 9 represents an extreme priority of one of the themes [30]. In Table 3, the order of the classes indicates the way we want to implement the priority. Saaty’s nine-scale levels for the delineation of groundwater capacity were used in a grid. It aims to capture what is unknown in opinions with AHP views [29]. Referring to the previous entry, the consistency index (CI) is a measure of the consistency by using Equation (1):
CI = λmax − n/n − 1
As a result, n, Lambda max or λmax is the most frequently occurring Eigenvalue of the pairwise comparison matrix. The measure of consistency of pairwise consistency (MC) is expressed as Equation (2):
CR = CI/RI
In addition to the Empirical Tables, we have a Ranking Tables section that allows RIs to be obtained by reference to different measures and countries that apply the indexes. We have Ranking Tables that include the Empirical Table, ratios, and countries that use those indexes as measures. Table 2 gives the value of RI for various values of n [31]. As long as the variability does not exceed 0.1, it is permissible. In cases where CR is higher than 10%, our judgments must be refined. The 10 pairwise consistencies are presented in Table 3. CR: 0.089; count value: 10.00; λmax: 11.179; CI: 0.131; CR: 0.09; constant: 1.49.
The physical multi-criteria model includes the following variables: slope, elevation, aspect, land use and land cover, geology, geomorphology, buffer zone of road, drainage buffer, built up area buffer, and forest type. Electronic clinometers were used to measure the slope of the lines on the map. After that, it was cross-checked with the digital elevation model of Cartosat 1. The data source for the topography maps is the Cartosat 1 digital elevation model with a spectral resolution of the panchromatic band of 2.5 m [19]. Table 4 shows the weighted value assigned to each class according to the fire hazards.
According to the AHP, participation is encouraged for a community decision to be ensured. The theory of AHP is one of hierarchical structure. The hierarchy makes it possible to evaluate the quality of an individual’s contribution at lower levels of the hierarchy. The role of the AHP is to help the decision-makers to choose the alternative with the greatest overall impact on the different priorities when they have to choose among many alternatives. Nonetheless, multi-criteria spatial issues such as land suitability are not mentioned in the land analysis. Currently, GIS-based multi-criteria analysis can be used to pick sustainable sites for new development. Many alternate decisions may be present, so the AHP helps to focus on the best of them. A literature survey demonstrates that the use of GIS for spatial problems has increased steadily over the last three decades. It provides validation for researchers for the realization of the value of a GIS-based multiple-criteria decision analysis or MCDA. Nowadays, computer technology is being increasingly applied to various decision-making tasks. As a result, it is important for geographic information systems (GIS) to be incorporated into sustainable development decisions.

3.2.1. Forest Type

Spatial data is the most important since it must be understood using a map. The radiometric maps are of the highest resolution for points far away from major landmasses due to reduced ground checking. The knowledge derived from remote sensing techniques has already been employed in India and abroad. The photo interpretation and georeferencing was performed. The forest types were as depicted in Figure 4, according to the visualization and translation of the data presented. Later on in the process, the layout was adjusted using topology.

3.2.2. Aspect Index Map

Due to the greater sunshine in the states to the southwest and southeast being particularly fire-prone, as well as direct sunlight in the south and southeast (Figure 5), those regions were given higher indices. The south is played up to full (Appendix A Table A1). Fires erupt and move faster in the southern regions where sunshine is prevalent.

3.2.3. Slope Index Map

As slopes grow steeper, the rate of fire spread increases due to convective pre-ignition and ignition. Larger slopes had more importance given to them in the mathematical model. Slopes ranging from 50 to 60 degrees and above have a weight of 5 as the upper limit. The gradient of this hose goes up, so there is less water to lose and it is more effective at convective pre-heating (Figure 6).

3.2.4. Road Buffer Index Map

Due to the road density, area index values were assigned to the closer areas. It is possible to come close to these hotspots, which can be a fire hazard. The traffic intersections were set up at a distance of 500 m intervals (Figure 7). Approximately 200 miles of white markers were required to delineate the study area in order to cover the full extent of road lines at 1000 m (Appendix A).

3.2.5. Elevation Index Map

There is less risk of fire at higher altitudes due to the effects of climate. When this research was carried out, the possibility of fire was calculated. It was then discovered that in this region, with a maximum elevation of 4700 m and lower elevations averaging 750 m, the chances of fire were greater (Figure 8). Hence, the maximum weight was used as the base for the findings. The last weight for elevations lesser than 1000 m is set to five pounds (Appendix A).

3.2.6. Built Up Buffer Index Map

Accidental and intentional forest fires were the origin of the wild land fires. A value assignment buffer of 200 m with five rings (Figure 9) has been developed with the index value mention in Appendix A according to their risk of forest fire. The buffer rings surrounding the homes were the most critical in containing a fire; thus, they were given greater weights near the citizens and lesser on the periphery.

3.2.7. Land Use and Land Cover

Initially, the image on Sentinel 2 [32] was resembled at a 10-m spatial resolution. It was attempted to correct the findings, and it was visually compared to the forest management chart on several occasions. The imagery is made up of six spectral bands from the approximate spectrum of different central wavelengths: band 2 (blue (0.490 μm)), band 3 (green (0.560 μm)), band 4 (red (0.665 μm)), band 8 (NIR (0.842 μm)), band 11 and 12 (SWIR (1.610 μm and 2.190 μm)) and land feature are made accordingly as shown in Figure 10. The risk analysis is shown in Table A1 (Appendix A).

3.2.8. Drainage Buffer Index Map

The drainage, slope, aspect and elevation are calculated from the Cartosat 1 data of 30-m resolution downloaded from the Bhuvan ISRO [33]. When assigning index values to the buffer class, areas that were further away from the drainage points had higher indexes. It is possible to come close to these hotspots that can be a fire hazard. Drainage was invented at a unit distance of 200 m (Figure 11). Drainage and grout pipes of various sizes have been employed, from 100 m up to 1000 m in order to cover the area (Appendix A).

3.2.9. Geomorphology

The area in the Sirmaur district presents a dazzling geographic mosaic of high mountain ranges, hills, and valleys with elevations ranging from 300 m to 3000 m and the source of data from the Bhukosh site of Survey India in 1:50,000 scale (Geological Survey of India [22]). The mountain peaks of denudational hill remain snow-covered all year round. For the Siwalik range of hills, which takes up the southwestern part of the district, the words “low denuded” would be appropriate. Table 5 shows the index value according to the risk of forest fire of the geomorphological features. The Piedmont zone of the area covered with 0.13% is very high and with 88.31% is high (Figure 12).

3.2.10. Geology

The Sirmaur district is located in the Shiwalik and lesser Himalayan ranges and the source of data from the Bhukosh site of Survey India is presented on a 1:50,000 scale. The rocks in the region are all over 600 million years old and were deposited in a succession of sandstone, shale, limestone, and schist, with both limestone and schist forming at the same time during that time period. These geological formations are noted for their limestone deposits from the Krol epoch. The southern region is occupied largely by the Oligocene Sea deposit. Himalayan districts span into the southern area of the rivers’ main flow to the Indus Brahmaputra. A total of 26.42% of the area is considered to be in the low-risk zone and 0.62% of the area is in the very high-risk zone mentioned in Appendix A and Figure 13.

4. Results and Discussion

Appendix A shows the index value according to the risk of forest fire of the geomorphological features (Figure 12). The geomorphology shows that the Piedmont zone of area covering 0.13% is very high and 88.31% is high. Geology features of the Sirmaur district show 26.42% of the area in the low-risk As–sb deposits zone and 0.62% in the very high-risk granites area zone. The most important variable is the forest type, in which pine has the greatest risk of fires as it has an abundance of flammable materials as 18.38% of the forested land in this area is made up of grass and litter. It is highly vulnerable to wildfires. Pine forests, because of the high oil content of turbine resins, are more liable to be damaged by fires than other types of forests. After examination of the nature of the forest, the density of the wood, the physical characteristics of the forest, the animals, and the duration of fires to develop a forest fire risk, the findings are found in Figure 12 and Appendix A, with Table 4 showing the area in percent under the no- to low-risk zone and the total is classified as very high-risk, 14.28% as moderate, and 10.44% as low-risk. The elevation index map shows that 28.40% of the area comes under an elevation ranging from 0 to 700 m and 13.39% of the area is under 700 to 1000 m. Figure 9 and Appendix A show a comparison of the risks and return areas in various elevations, respectively. The slope variables show that 17.02% and 23.60% of all falls are from 0 to 10° and 10 to 20° angles, rather than just 10 to 30° (moderate risk), 30 to 40° (strong risk) and >40° (very high risk). There is a steepness to the gradient that is crucial to the fire’s spread, visible in the table and Figure 5 and Figure 10. The aspect index shows the direction of elevation in which the area exhibits about 14.00% in the south, 14.16% in the south-west, 11.97% facing the south-east, 11.53% towards the east, 23.88% under the north-west and west, and 24.47% under the north and the north-east. The risk index area chart data is shown in Figure 9. The closest-to-the-road buffers were given a 200 m limit. Figure 8 depicts the road network and buffer location for the purpose of analysis. Table 4 shows data according to their risk of a forest fire. The buffer rings surrounding the streams had no risk for containing a fire; thus, they were given lower weights near the drainage and higher on the periphery. Figure 7 depicts the results of the drainage and buffer planning calculations for forest fire risk. A value assignment buffer of 200 m with five rings has been developed with the index value shown in Table 4 according to their risk of a forest fire. The buffer rings surrounding the homes were the most critical in containing a fire; thus, they were given greater weights near the citizens and lesser on the periphery. Figure 6 depicts the results of the built-up land and buffer planning calculations. Figure 4 shows the land use and land cover classification risk zone; Figure 12 shows the area covered in the risk zone. A total of 51.86% of the area comes under forest cover, which is under the risk zone, 2.81% of the area is covered under no-risk zone, and 22.52% of the area comes under the low-risk zone.
The final result, according to the flow chart (Figure 3) methodology, shows that 12.13% of the area has a very high-risk of forest fire and 12.53% of the area has a high risk. A total of 53.92% of the area is covered under the low- to no-risk zone (Table 5). Red zones present in the map (Figure 14) are very high-risk-prone areas for forest fires. Hence, the Sirmaur district should be considered accordingly under this forest fire risk management plan by the government, local bodies, and NGOs, respectively. These classes shown in Table 5 could be validated using the data on wildfires for each category of forest fire risk index. Validation points are shown on Figure 15. A false alarm was detected in the north; it was one of 43 (or 2.33%) of the locations in the group, with an accuracy of 97.67%. MODIS data is downloaded from the Earthdata Firms site, from their archival data. By downloading MODIS data of 2019, the brightness value was converted into a shaped file. About 96% of the data shows accurate findings in the calculation of the forest fire risk zone. Table 5 shows that most of the points are shown in the moderately high- to very high-risk zone. Figure 15 shows the validation points on the forest fire risk zone and Figure 16 shows the brightness value range of the MODIS data. The three MODIS test site validation exercises completed in 2007 all had a final accuracy of 93.06% in the identification of forest fires. Additionally, the 2007–2008–2009 final organizational assessment in Thailand revealed a 95.64% level of confidence. In the India study of the Uttarakhand state through MODIS, the brightness value gives a 97% accuracy of wildfire zone assessment.

5. Conclusions

From the aforementioned result and discussion, we can conclude that geospatial technology can be used to quantify biomass, loss of forest, and loss of biodiversity due to forest fires in the future. Using high-resolution ground-satellite data, the forest fire risk-assessment will assist fire fighters serving in the location of administrative bodies such as local municipalities. We will measure the carrying capacity of herby habitat and use geospatial modeling to measure where the migratory animals reside in order to see whether there are any issues with these corridors being adequate. Such a model will work in the initial stages and will subsequently have to be devised to apply, for example, biomass density or degradation ratios, to the total human activity and biophysical parameters. Another kind of model could be created, allowing for both the total biomass loss and forest area change to be calculated including human-caused fire danger that exists in the region. Due to this, we have carried out research on the effect of climate changes on the current population, and socioeconomic location on people’s vulnerability. Census data will be used to produce a socio-economic analysis of the area, as well as an analysis of fire hazards through the computation of population density, agricultural laborers’ age demographics, and the number of children in the zero- to six-year-old age group. Every year in India, global warming has caused a rise in forest fires. There is an expected correlation between increased carbon emissions and boosts in burnt area and the number of fires. It can be shown that 12.13% of the area is under a very high risk of forest fires and 12.53% of the area is under a high risk. A total of 53.92% of the area is considered to be in the low- to no-risk zone. Validation of the final result by the MODIS brightness value shows the calculation of the forest fire risk zone and a 96% accuracy rate, which can be further utilized in these zones by not being affected from forest fire.
It has been shown that policies and activities aimed at empowering community participation are safer and produce significantly reduced risk than those previous wildfire preparations and executions. Research increasingly suggests that scientific research and planning also help communities in the West to deal with the danger of wildfires. In addition to the wildfire problems increasing, as more citizens migrate to riskier locations, it is important to address these inadequacies between the best practice and the importance of integrating wildfires in community decision-making. Thus, this research supports the notion that community development and advanced preparation strategies will help to deal with wildfire incidents. An adapted geospatial model was employed in the development of the forest fire hazard in the Sirmaur district of the Himalayan provinces of India. The greatest risk is seen in areas where the fire threat is moderate to high, and we need to be expeditious to prevent it from spiraling of control.

Author Contributions

Conceptualization, J.S.T. and N.K.; methodology, J.S.T.; software, S.K.; validation, N.K., B.Đ. and J.S.T.; formal analysis, S.K. and S.K.S.; investigation, J.S.T.; resources, J.S.T.; data curation, J.S.T.; writing—original draft preparation, J.S.T.; writing—review and editing, N.K.; visualization, S.K.; supervision, J.S.T.; funding acquisition, N.K. and B.Đ. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available at: European Space Agency web page Available online: http://www.esa.int/ (accessed on 24 May 2021). CARTOSAT-1 Available online: https://www.isro.gov.in/Spacecraft/cartosat-1 (accessed on 24 May 2021). Sentinel-2 Available online: https://sentinels.copernicus.eu/web/sentinel/missions/sentinel-2 (accessed on 24 May 2021). Geological map of India available online: https://www.gsi.gov.in/webcenter/portal/OCBIS/pageMAPS/pageMapsSeries?_adf.ctrl-state=eldqgtyzy_9&_afrLoop=22186242333769660#! (accessed on 24 May 2021).

Acknowledgments

The authors would like to express their gratitude to the reviewers for greatly improving the quality of the manuscript. The authors would like to thank Shruti and Suraj for improving the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Fire hazard classes.
Table A1. Fire hazard classes.
VariableClassesIntensity of ImportanceFire Hazards ClassesArea in km2Area in %
Land Use/Land coverAgriculture5High605.8521.42
Built-up Land2Low Risk12.110.43
Forest7Very High1466.9451.86
Grass Land/Grazing Land4Moderately High39.331.39
Waste Land3Low Risk624.6822.09
Water Bodies1No Risk79.492.81
GeologyAs–Sb Deposits2Low Risk372.96113.17
BlainiManjir formation1No Risk123.614.37
Granitoids7Very High17.60210.62
Jatog Groups4Moderately High257.0359.08
Shimla and Jaunsar groups3Moderately High375.18913.25
Siwalik Groups6High146151.60
Tal Kunzaml Thango5High120.4194.25
Unconsolidated Glacial1No Risk103.6713.66
GeomorphologyAlluvial Plain2Low Risk187.0526.63
Denudational Hills3Low Risk67.04372.38
Flood Plain1No Risk72.02142.55
Piedment Zone7Very High3.737770.13
Structural Hills5Moderately High2492.3988.31
Forest TypePlantation/TOF4Moderately High2.710.10
Lower or Siwalik chir pine forest9Very High518.3918.38
Non-Forest1No Risk1386.6649.16
Moist deodar forest2Low Risk178.896.34
Oak scrub3Low Risk103.483.67
Northern dry mixed deciduous forest5High203.467.21
Water1No Risk11.060.39
Moist temperate deciduous forest2Low Risk11.890.42
West Himalayan sub-alpine forest3Low Risk0.390.01
West Himalayan high-level dry blue pine forest7Very High6.730.24
Bhabar-dun sal forest4Moderately High397.1214.08
Drainage Buffer<200 m1No Risk483.24317.09
400 m2Low Risk444.22915.71
600 m3Moderately High400.58914.16
800 m4Moderately High354.95112.55
1000 m5High311.51311.01
>1000 m7Very High833.78329.48
Road Buffer<200 m7Very High907.49932.09
400 m5High582.69320.60
600 m3Moderate405.018147.61
800 m3Moderately High274.3929.70
1000 m1Low Risk185.1526.55
>1000 m1No Risk473.34416.74
Built-Up Buffer<200 m7Very High70.37442.49
400 m5High78.26072.77
600 m4Moderately High96.11693.40
800 m3Moderately High107.6463.81
1000 m2Low Risk113.9954.03
>1000 m1No Risk2362.0583.51
Slope<102Low Risk506.4617.02
10.1–203Moderately High702.18623.60
20.1–304Moderately High960.04432.27
30.1–405High646.02121.71
40>7Very High160.6285.40
AspectsNorth (0–22.5); North (337.5–360); Northeast (22.5–67.5)1No Risk727.49424.47
East (67.5–112.5)3Moderately High342.67911.53
Southeast (112.5–157.5)4Moderately High355.77511.97
South (157.5–202.5)7Very High416.38914.00
Southwest (202.5–247.5)5High420.90214.16
West (247.5–292.5); Northwest (292.5–337.5)2Low Risk710.09223.88
Elevation0–7006Very High844.95528.40
700.1–10004High398.44813.39
1000.1–15003Moderately High927.63431.18
1500.1–20002Low Risk570.53519.18
2000.1–35361No Risk233.7687.86

References

  1. Bahuguna, V.K. India: Forest Fire Prevention and Control Strategies (IFFN No. 20). Available online: https://gfmc.online/iffn/country/in/in_1.html#:~:text=Forest Fire prevention is looked after by the,provided financial assistance for fire prevention and control (accessed on 27 May 2021).
  2. United Nations. Sendai Framework for Disaster Risk Reduction 2015–2030; General Assembly, United Nations: Geneve, Switzerland, 2015. [Google Scholar]
  3. Dimopoulou, M.; Giannikos, I. Towards an integrated framework for forest fire control. Eur. J. Oper. Res. 2004, 152, 476–486. [Google Scholar] [CrossRef]
  4. Semwal, R.L. Forest Fires in India: Lessons from Case Studies; New Delhi World Wildlife Fund for Nature-India: New Delhi, India, 2003. [Google Scholar]
  5. Srivastava, P.; Garg, A. Forest Fires in India: Regional and Temporal Analyses. J. Trop. For. Sci. 2013, 25, 228–239. [Google Scholar]
  6. Jain, A.; Ravan, S.A.; Singh, R.K.; Das, K.K.; Roy, P.S. Forest fire risk modelling using remote sensing and geographic information system. Curr. Sci. 1996, 70, 928–933. [Google Scholar]
  7. Jaiswal, R.K.; Mukherjee, S.; Raju, K.D.; Saxena, R. Forest fire risk zone mapping from satellite imagery and GIS. Int. J. Appl. Earth Obs. Geoinf. 2002, 4, 1–10. [Google Scholar] [CrossRef]
  8. Alonso-Betanzos, A.; Fontenla-Romero, O.; Guijarro-Berdiñas, B.; Hernández-Pereira, E.; Inmaculada Paz Andrade, M.; Jiménez, E.; Luis Legido Soto, J.; Carballas, T. An intelligent system for forest fire risk prediction and fire fighting management in Galicia. Expert Syst. Appl. 2003, 25, 545–554. [Google Scholar] [CrossRef]
  9. Iliadis, L.S. A decision support system applying an integrated fuzzy model for long-term forest fire risk estimation. Environ. Model. Softw. 2005, 20, 613–621. [Google Scholar] [CrossRef]
  10. Bonazountas, M.; Kallidromitou, D.; Kassomenos, P.A.; Passas, N. Forest Fire Risk Analysis. Hum. Ecol. Risk Assess. Int. J. 2005, 11, 617–626. [Google Scholar] [CrossRef]
  11. Adab, H.; Devi, K.; Karim, K. Modeling forest fire risk in the northeast of Iran using remote sensing and GIS techniques. Nat. Hazards 2013, 65, 1723–1743. [Google Scholar] [CrossRef]
  12. Malik, T.; Rabbani, G.; Farooq, M. Forest Fire Risk Zonation Using Remote Sensing and GIS Technology in Kansrao Forest Range of Rajaji National Park, Uttarakhand, India. Int. J. Adv. Remote Sens. GIS 2013, 2, 86–95. [Google Scholar]
  13. Jung, J.; Kim, C.; Jayakumar, S.; Kim, S.; Han, S.; Kim, D.H.; Heo, J. Forest fire risk mapping of Kolli Hills, India, considering subjectivity and inconsistency issues. Nat. Hazards 2013, 65, 2129–2146. [Google Scholar] [CrossRef]
  14. Sowmya, S.V.; Somashekar, R.K. Application of remote sensing and geographical information system in mapping forest fire risk zone at Bhadra wildlife sanctuary, India. J. Environ. Biol. 2010, 31, 969–974. [Google Scholar]
  15. Ajin, R.; Loghin, A.-M.; Vinod, P.; Jacob, M. Forest Fire Risk Zone Mapping Using RS and GIS Techniques: A Study in Achankovil Forest Division, Kerala, India. J. Earth Environ. Health Sci. 2016, 2, 109–115. [Google Scholar] [CrossRef]
  16. Kanga, S.; Tripathi, G.; Singh, S.K. Forest Fire Hazards Vulnerability and Risk Assessment in Bhajji Forest Range of Himachal Pradesh (India): A Geospatial Approach. J. Remote Sens. GIS 2017, 8, 25–40. [Google Scholar]
  17. Government of Himachal Pradesh District Sirmaur-Demography. Available online: https://hpsirmaur.nic.in/demography/ (accessed on 24 May 2021).
  18. European Space Agency Web Page. Available online: http://www.esa.int/ (accessed on 27 May 2021).
  19. Government of India of CARTOSAT-1. Available online: https://www.isro.gov.in/Spacecraft/cartosat-1 (accessed on 24 May 2021).
  20. Rogeau, M.-P.; Armstrong, G.W. Quantifying the effect of elevation and aspect on fire return intervals in the Canadian Rocky Mountains. For. Ecol. Manag. 2017, 384, 248–261. [Google Scholar] [CrossRef]
  21. Agency, E.S. Sentinel-2. Available online: https://sentinels.copernicus.eu/web/sentinel/missions/sentinel-2 (accessed on 24 May 2021).
  22. Geological Survey of India Geological Survey of Geological Map of India. Available online: https://www.gsi.gov.in/webcenter/portal/OCBIS/pageMAPS/pageMapsSeries?_adf.ctrl-state=eldqgtyzy_9&_afrLoop=22186242333769660#! (accessed on 27 May 2021).
  23. Shakesby, R.A.; Doerr, S.H. Wildfire as a hydrological and geomorphological agent. Earth Sci. Rev. 2006, 74, 269–307. [Google Scholar] [CrossRef]
  24. Meyer, G.A.; Pierce, J.L. Climatic controls on fire-induced sediment pulses in Yellowstone National Park and central Idaho: A long-term perspective. For. Ecol. Manag. 2003, 178, 89–104. [Google Scholar] [CrossRef]
  25. Armenteras, D.; Barreto, J.S.; Tabor, K.; Molowny-Horas, R.; Retana, J. Changing patterns of fire occurrence in proximity to forest edges, roads and rivers between NW Amazonian countries. Biogeosciences 2017, 14, 2755–2765. [Google Scholar] [CrossRef] [Green Version]
  26. Tian, X.; Zhao, F.; Shu, L.; Wang, M. Distribution characteristics and the influence factors of forest fires in China. For. Ecol. Manag. 2013, 310, 460–467. [Google Scholar] [CrossRef]
  27. Fazlollahtabar, H.; Eslami, H.; Salmani, H. Designing a Fuzzy Expert System to Evaluate Alternatives in Fuzzy Analytic Hierarchy Process. J. Softw. Eng. Appl. 2010, 3, 409–418. [Google Scholar] [CrossRef] [Green Version]
  28. Saaty, T.L. How to make a decision: The analytic hierarchy process. Eur. J. Oper. Res. 1990, 48, 9–26. [Google Scholar] [CrossRef]
  29. Saaty, T.L. Analytic Heirarchy Process. In Wiley Stats Ref: Statistics Reference Online; Wiley Online Library, John Wiley & Sons: Hoboken, NJ, USA, 2014; ISBN 9781118445112. [Google Scholar]
  30. Saaty, T.L. Rank from comparisons and from ratings in the analytic hierarchy/network processes. Eur. J. Oper. Res. 2006, 168, 557–570. [Google Scholar] [CrossRef]
  31. Marinoni, O. Implementation of the analytical hierarchy process with VBA in ArcGIS. Comput. Geosci. 2004, 30, 637–646. [Google Scholar] [CrossRef]
  32. Copernicus Open Access Hub. Available online: https://scihub.copernicus.eu/ (accessed on 27 May 2021).
  33. National Remote Sensing Centre Bhuvan-3D. Available online: https://bhuvan.nrsc.gov.in/home/index.php (accessed on 27 May 2021).
Figure 1. Study area.
Figure 1. Study area.
Ijgi 10 00447 g001
Figure 2. Rainfall (in mm) in the Sirmaur district.
Figure 2. Rainfall (in mm) in the Sirmaur district.
Ijgi 10 00447 g002
Figure 3. Flow chart of forest fire area calculation.
Figure 3. Flow chart of forest fire area calculation.
Ijgi 10 00447 g003
Figure 4. Forest type map.
Figure 4. Forest type map.
Ijgi 10 00447 g004
Figure 5. Aspect buffer map.
Figure 5. Aspect buffer map.
Ijgi 10 00447 g005
Figure 6. Slope map.
Figure 6. Slope map.
Ijgi 10 00447 g006
Figure 7. Road buffer map.
Figure 7. Road buffer map.
Ijgi 10 00447 g007
Figure 8. Elevation map.
Figure 8. Elevation map.
Ijgi 10 00447 g008
Figure 9. Built up buffer risk index.
Figure 9. Built up buffer risk index.
Ijgi 10 00447 g009
Figure 10. Land use/land cover map.
Figure 10. Land use/land cover map.
Ijgi 10 00447 g010
Figure 11. Drainage buffer risk index.
Figure 11. Drainage buffer risk index.
Ijgi 10 00447 g011
Figure 12. Geomorphology map.
Figure 12. Geomorphology map.
Ijgi 10 00447 g012
Figure 13. Geology map.
Figure 13. Geology map.
Ijgi 10 00447 g013
Figure 14. Forest fire risk map.
Figure 14. Forest fire risk map.
Ijgi 10 00447 g014
Figure 15. Forest fire validation hotspots in forest fire risk zone map.
Figure 15. Forest fire validation hotspots in forest fire risk zone map.
Ijgi 10 00447 g015
Figure 16. Modis point with its brightness values.
Figure 16. Modis point with its brightness values.
Ijgi 10 00447 g016
Table 1. The ratio scale and definition of AHP [28].
Table 1. The ratio scale and definition of AHP [28].
ImportanceDefinitionExplanation
1Equally ImportantEqually vital to the target
3Moderately ImportantCompared to the overall profit or damage
5Strong ImportanceA strong preference for one factor over another
7Very Strong ImportanceThe one thing that has gained preeminence, considered to be above all the others and vastly superior in the real world, is the theory world of practice.
9Extreme ImportanceIf it is strongly proven with evidence and facts, then one element is favored in comparison to the other
2, 4, 6, 8Inter values
Table 2. Random Consistency Index (RI) [28].
Table 2. Random Consistency Index (RI) [28].
n123456789101112131415
RI000.580.91.121.241.321.411.451.491.511.531.561.571.59
Table 3. Pairwise comparisons.
Table 3. Pairwise comparisons.
Item Number12345678910
Item NumberItem DescriptionForest TypeAspectSlopeRoad BufferElevationBuilt up LandLand Use/Land CoverDrainage BufferGeomorphologyGeology
1Forest Type1.005.003.005.004.003.004.005.003.001.00
2Aspect0.201.002.001.000.504.001.005.001.002.00
3Slope0.330.501.001.001.003.001.006.003.001.00
4Road Buffer0.201.001.001.001.002.001.002.004.002.00
5Elevation0.252.001.001.001.001.001.002.002.002.00
6Built up land0.330.250.330.551.001.000.504.000.332.00
7Land use/land cover0.251.001.001.001.002.001.002.001.001.00
8Drainage Buffer0.200.200.170.500.500.250.501.003.000.33
9Geomorphology0.331.000.330.250.503.001.003.001.001.00
10Geology1.000.501.000.500.500.501.003.001.001.00
Sum4.1012.4510.8311.7511.0019.7512.0030.3319.3313.33
Table 4. Weighted values assigned to each class.
Table 4. Weighted values assigned to each class.
Item NumberVariableForest TypeAspectSlopeRoad BufferElevationBuilt up LandLand Use/Land CoverDrainage BufferGeomorphologyGeologyWeight
1Forest Type0.240.400.280.430.360.150.330.160.160.0825.9%
2Aspect0.050.080.180.090.050.200.080.160.050.1511.0%
3Slope0.080.040.090.090.090.150.080.200.160.0810.5%
4Road Buffer0.050.080.090.090.090.100.080.070.210.1510.0%
5Elevation0.060.160.090.090.090.050.080.070.100.159.4%
6Built up land0.080.020.030.040.090.050.040.130.020.156.6%
7Land use/land cover0.060.080.090.090.090.100.080.070.050.087.9%
8Drainage Buffer0.050.020.020.040.050.010.040.030.160.034.4%
9Geomorphology0.080.080.030.020.050.150.080.010.050.086.3%
10Geology0.240.040.090.040.050.030.080.100.050.088.0%
Table 5. Forest fire risk area distribution with validation points.
Table 5. Forest fire risk area distribution with validation points.
Forest Fire RiskArea (in km2)Area (in %)Index ValueValidation PointsBrightness Range Value
Very high339.15212.1355343.8–327.5
High350.44212.5347327.4–323.9
Moderately high599.35221.43311323.4–306.1
Low risk924.91033.0723303.1–303.6
No risk583.09520.8511301–303
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tomar, J.S.; Kranjčić, N.; Đurin, B.; Kanga, S.; Singh, S.K. Forest Fire Hazards Vulnerability and Risk Assessment in Sirmaur District Forest of Himachal Pradesh (India): A Geospatial Approach. ISPRS Int. J. Geo-Inf. 2021, 10, 447. https://doi.org/10.3390/ijgi10070447

AMA Style

Tomar JS, Kranjčić N, Đurin B, Kanga S, Singh SK. Forest Fire Hazards Vulnerability and Risk Assessment in Sirmaur District Forest of Himachal Pradesh (India): A Geospatial Approach. ISPRS International Journal of Geo-Information. 2021; 10(7):447. https://doi.org/10.3390/ijgi10070447

Chicago/Turabian Style

Tomar, Jagpal Singh, Nikola Kranjčić, Bojan Đurin, Shruti Kanga, and Suraj Kumar Singh. 2021. "Forest Fire Hazards Vulnerability and Risk Assessment in Sirmaur District Forest of Himachal Pradesh (India): A Geospatial Approach" ISPRS International Journal of Geo-Information 10, no. 7: 447. https://doi.org/10.3390/ijgi10070447

APA Style

Tomar, J. S., Kranjčić, N., Đurin, B., Kanga, S., & Singh, S. K. (2021). Forest Fire Hazards Vulnerability and Risk Assessment in Sirmaur District Forest of Himachal Pradesh (India): A Geospatial Approach. ISPRS International Journal of Geo-Information, 10(7), 447. https://doi.org/10.3390/ijgi10070447

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop