Next Article in Journal
Evacuation Simulation in a Cultural Asset Fire: Impact of Fire Emergency Evacuation Facilities for People with Disabilities on Evacuation Time
Next Article in Special Issue
Effects of Post-Fire Management on a Mediterranean Small Mammal Community
Previous Article in Journal
Combustion and Emission Analysis of Spent Mushroom Compost and Forestry Woodchip for Management and Energy Production
Previous Article in Special Issue
Remote Sensing-Based Monitoring of Postfire Recovery of Persistent Shrubs: The Case of Juniperus communis in Sierra Nevada (Spain)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Model to Estimate the Risk of Emission of Greenhouse Gases from Forest Fires

by
Victoria Lerma-Arce
1,*,
Celia Yagüe-Hurtado
1,
Helena Van den Berg
1,
Miguel García-Folgado
1,
Jose-Vicente Oliver-Villanueva
1,
Yacine Benhalima
2,3,
Inês Marques-Duarte
2,
Vanda Acácio
2,
Francisco C. Rego
2,
Eduardo López-Senespleda
4,
María Menéndez-Miguélez
4,
Ricardo Ruiz-Peinado
4,
Thomas Petillon
5,
Stéphanie Jalabert
5,
Ester Carbó-Valverde
6,
Eugenia Gimeno-García
6,
Rebeca Aleix-Amurrio
7 and
Edgar Lorenzo-Sáez
1
1
Institute of Information and Communication Technologies (ITACA), Universitat Politècnica de València (UPV), Camino de Vera s/n, 46022 Valencia, Spain
2
InBio, Centro de Ecologia Aplicada “Professor Baeta Neves”, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisbon, Portugal
3
Associated Laboratory TERRA, LEAF—Linking Landscape, Environment, Agriculture and Food—Research Center, Instituto Superior de Agronomia, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisbon, Portugal
4
Instituto de Ciencias Forestales (ICIFOR-INIA), CSIC, 28040 Madrid, Spain
5
Bordeaux Science Agro, UMR CNRS 5805 EPOC, 1 cours du Général De Gaulle, 33175 Gradignan, France
6
Department of Environmental Quality and Soils, Centro de Investigaciones sobre Desertificación—CIDE (CSIC-Universitat de Valencia-GV), 46113 Valencia, Spain
7
Asociación de Municipios Forestales de la Comunitat Valenciana (AMUFOR), 46810 Valencia, Spain
*
Author to whom correspondence should be addressed.
Submission received: 10 November 2022 / Revised: 2 December 2022 / Accepted: 17 December 2022 / Published: 29 December 2022
(This article belongs to the Special Issue Mediterranean Fires)

Abstract

:
While the Mediterranean basin is foreseen to be highly affected by climate change (CC) and severe forest fires are expected to be more frequent, international efforts to fight against CC do not consider forest fires’ greenhouse gas (GHG) emissions risk and the possibility of its mitigation. This is partly due to a lack of a methodology for GHG risk spatial assessment and consideration of the high value of carbon stocks in forest ecosystems and their intrinsic risk. To revert this, an innovative GHG emission risk model has been developed and implemented in a pilot forest area. This model considers geospatial variables to build up emission vulnerability based on potential fire severity and resistance of a landscape, value at risk and the hazard of a fire occurrence. The results classify low, moderate and high emission risks in the analysed areas. This identification of hotspots allows the prioritisation of fire prevention measures in a region to maximise the reduction of GHG emissions in the case of a fire event. This constitutes the first step in a holistic and consistent CC mitigation that not only considers anthropic GHG sources but also possible GHG emissions by forest fires that can be actively prevented, managed and reduced.

1. Introduction

Wildfires are a dominant disturbance factor in almost all forest vegetation zones throughout the world and represent important emissions of gaseous and particulate compounds into the atmosphere with serious physical, biological and environmental impacts [1].
The Mediterranean region is one of the “hot-spots” of climate change [2]. In this area, wildfires have a significant impact, burning approximately 0.5 million every year in the five southern EU member states (Portugal, Spain, Greece, Italy and France) [3].
The Mediterranean region is characterised by hot and dry summers, high diversity of plant species and unusual geographical/topographical variability related to the presence of a jagged coastline and many mountain ranges, often rather steep. Since 1960, wildfire occurrence has increased because of changes in land use, which resulted in extensive land abandonment, increases in the fuel load and continuity in the landscape [4]. These factors are considered to be the main drivers of increased wildfire occurrence and the resulting risk of soil degradation [5,6,7,8,9,10,11,12].
Moreover, the current climate has a large influence on wildfires, especially in extreme conditions [13,14]. An increase in the duration and severity of heat waves [15,16] has caused a dramatic increase in fire incidence in Mediterranean regions during recent decades, and fires are expected to become more prevalent in the future due to climate change [17].
When weather conditions are particularly severe (e.g., heat waves or very hot and dry summers combined with strong winds) forest fires may reach catastrophic proportions, so-called megafire events, producing not only a loss of biodiversity, soil erosion and desertification [9], large economic losses and threats to human lives but also noteworthy GHG emissions in the atmosphere [18].
The greenhouse gases and particulate matters directly influence climate [19,20] due to the release of various types of gases and aerosols into the atmosphere because of the combustion of biomass during forest fires, which have a significant impact on atmospheric chemistry, biogeochemical cycles and climate [21]. Not only do many of these gases contribute to climate change and the process of the greenhouse effect, but they can also trigger many other repercussions. This disturbance also implies a change in the dynamics of the gross primary production and ecosystem respiration. It also leads to net ecosystem production losses continuing for several years after the event, which is followed by a sustained multidecadal period of net ecosystem carbon uptake [22].
It is known that tremendous efforts are being performed at international political, regulatory and technological levels to reduce carbon emissions, with binding and non-binding compromises and high objectives of emission reductions on a worldwide scale [23,24,25] and at the European level [26,27,28], but there is no common policy on forest fire management and its prevention.
Recently, in 2021, the European Commission adopted a New Strategy for the European Union in favour of forests for 2035 [29] in which it is considered necessary to reduce the risks for forests in the context of uncertainty due to climate change through appropriate adaptation measures and forest management practices that strengthen resilience. For this, not only must technical knowledge be developed but also incentives and regulatory support must be offered.
There are many studies on ecological risk assessments [30] related to fire hazards, assessing the risks to ecosystem services from forest fires [31] or fire smoke risk on populations [32]. There is also very extensive literature on fire risk assessment [33,34,35,36,37] and fire danger in the Mediterranean basin using coarse-scale fire danger indices, such as the Canadian Forest Fire Weather Index (FWI) [38,39], the Wildland Fire Assessment System (WFAS) [40], the Forest Fire Danger Index (FFDI, also called Mark 5 the Burning Index, BI), the KBDI [41], etc. These indices rely solely on weather information or statistical models of climate and vegetation [42,43,44,45,46,47,48].
The fine-scale patterns in fire behaviour are important from a risk standpoint and can be assessed using a number of innovative approaches developed in recent decades for estimating and mapping wildfire risk and exposure [37,42,49,50,51,52]. These approaches use quantitative wildfire risk assessment [51] and employ a number of relatively new, large, fire modelling systems. Nevertheless, there are no studies assessing greenhouse gas emission risk from a preventive perspective with the aim of managing it to reduce emission risk impacts on the climate from forest fires. Therefore, the general objective of this study is to set the basis for a qualitative spatial assessment of forest fire emissions risk.

2. Materials and Methods

2.1. Basis of the Model

IPCC [53] defines the concept of risk as “the potential for adverse consequences for human or ecological systems, recognising the diversity of values and objectives associated with such systems”. According to [54], risk can be defined as a combination of hazard and damage; therefore, the Emission Risk assessment model can be defined as (Equation (1)):
Emission Risk = Hazard × Damage
where hazard (H) is the potential forest fire risk and shows the probability of a fire event [55] and depends on two variables: Statistical risk and site danger (composed of fire recurrence, cause and exposure). In this model, a hazard is taken as the probability of emissions due to forest fires (recurrence), expressed as the fire frequency as it indicates the risk of ignition. The cause is not relevant in terms of the risk of emission, and the exposure to fire is included as part of the other variables.
Damage (D) is the consequence in terms of GHG emissions from the stored carbon in the analysed ecosystem as a result of combustion produced by a forest fire [54]. The damage degree depends on the value at risk, in this case, total carbon from vegetation and soils stored in the ecosystem, and its vulnerability to emitting GHG emissions.
Damage = Value at risk × Emission vulnerability
Value at risk (VAR) is defined as the potential loss due to an adverse event [54]. In the REMAS context, the VAR variable includes carbon fixed in both vegetation and soils. It is the carbon stock of the ecosystem and represents the total carbon stored, which is available to be released into the atmosphere as CO2 emissions as a consequence of a forest fire [18].
Emission Vulnerability (EV), as defined above, is the degree of loss or damage that carbon stored in an ecosystem may cause in a wildfire event, materialised in GHG emissions. According to [56], this process depends on two main factors: The resistance of the ecosystem to fire ignition and propagation through the stand (potential fire resistance) and the degree of loss based on the intensity of the fire (potential fire severity).

2.2. Factors of Interest

The effects of carbon emissions are influenced by fire intensity and severity [57], which depend on an interaction between different factors that can vary a great deal depending on the place and the time of the year. These variables are considered in many studies [33,34,36,42,58] and are summarised as topography, weather conditions, forest fuel types and vegetation.
On the other hand, the impact of forest fires can be controlled and reduced with silvicultural treatments. Studies have demonstrated that fuel reduction treatments, among others, can reduce carbon losses in the case of a wildfire [59,60,61,62].

2.2.1. Fire Frequency

This is the indicator of the probability of a fire occurring based on the quantification of forest fires registered in a period of time and a determined spatial unit [63]. According to [18,58], the higher the fire frequency, the higher the emission risk envisaged. Its calculation is performed according to [64]. Its value assessment in terms of emission risk can be found in Table A1.

2.2.2. Ecosystem Carbon Stock

This variable is composed of the different carbon sinks that can be found and from which data are available in forest ecosystems. According to [65], the more carbon is stored, the more emission risk is envisaged. Its value assessment in terms of emission risk can be found in Table A2.
  • Vegetation carbon stock
It considers the total carbon stored in vegetation as tonnes (Mg) of carbon per hectare. Databases are obtained from the quantification of carbon stored in the main forest tree species representative of the study areas. The estimation of these vegetation carbon stocks can be obtained following several methodologies [66,67].
  • Soil organic carbon stock
With soil organic carbon, we define the carbon stocks, not the soil carbon content. Soil organic carbon is the carbon contained in soil organic matter, but stocks refer to a volume. Therefore, according to [68], it considers the mass of carbon in a sample of known bulk density. Soil organic carbon stocks are generally expressed in tonnes (Mg) per hectare for a nominated depth (from 0–30 cm depth) and commonly restricted to the fraction [68] n < 2 mm in size [68]. Available soil carbon stock maps were used to estimate it.

2.2.3. Preventive Silviculture (PS)

Preventive silviculture manages forests by enhancing their capacity to protect themselves from fires by creating discontinuities, avoiding very extensive, monospecific surface areas and creating a patchwork of different inflammability levels that disturb the fire [69]. This can be achieved through fuel treatments for biomass reduction, which are paramount to wildfire abatement [70]. In the model, two main factors have been considered:
  • Prescribed burns (PB)
This can be defined as the controlled use of fire to reduce vegetation under specific conditions that allow for setting the intensity of the fire and the amount of vegetable fuel to be eliminated according to a proposed objective [71]. The existence of prescribed burns is foreseen to lead to lower emission risk [17,71].
  • Vegetation treatments (VT)
The objective of vegetation treatment is to improve the conditions of life, growth and health of the forest stand. According to [72], it helps to reduce fire spread because it reduces vertical and horizontal biomass continuity. These practices include treatments applied to tree and shrub vegetation but also to the soil [72]. The main fire preventive treatments are auxiliary strips, tree clearing and thinning and shrub clearing.
Both PB and VT are assessed in terms of emission risk according to Table A3. The existence of vegetation treatments is envisaged to lead to lower emission risk [60,73,74].

2.2.4. Infrastructure (I)

This consists of artificial facilities for forest fire prevention. This infrastructure can be of help in the case of forest fires with dimensions within the extinction capacity of forest services [75]. In this case, firebreaks and fuelbreaks, road infrastructure and the water supply network accounted for the model with the following approach:
  • Firebreaks and fuelbreaks network (FN)
Firebreaks are defined as strips or elongated spaces with a width of 20 to 30 m in which all types of vegetation are removed, exposing mineral soil [72], while fuelbreaks have a biomass decrease associated with a fuel model change. Both nets are constructed on artificial lines, such as paths, boundaries of mountains or forest planning units; on natural lines such as maximum line slopes coinciding with the separation of rain gaps; and on summit lines, although in this case, they should be located not on the hills but rather in areas set back above the beginning of the slopes, where the wind speed is relatively low [69,72]. They have a double objective of acting as a barrier for the propagation of the fire and as a physical support for firefighters and extinction media to perform their work.
  • Road network (RN)
Different types of forest roads (road infrastructure) are distributed throughout the terrain. These paths have a double impact in relation to forest fires. On one hand, they act as a physical barrier to the spread of the fire. On the other hand, they act as an operative point where firefighters and extinction media can perform their work [76].
The firebreaks and forest road networks are merged in order to be included in the model since they often designate the same physical element. The existence of firebreaks and forest roads is envisaged to lead to lower emission risk [76,77].
  • Water supply network (WSN)
There are different types of water sources distributed throughout the terrain in order to cover the maximum surface to provide water in the case of fire [78]. This network is of help for fire extinction bodies in the case of a fire event, and therefore its presence contributes to fire suppression with a consequent emission reduction.
In order to calculate the area of influence of the water supply points, the buffer area established was 2500 m. The existence of water supply points is foreseen to lead to lower emission risk [78].
FN, RN and WSN were assessed in terms of the emission risk according to Table A3.

2.2.5. Spatial Heterogeneity (SH)

This is the spatial variability of different ecosystems in a determined landscape. This mosaic acts as a natural barrier opposing fire propagation [79]. It depends on two main factors: The spatial distribution of land uses in a determined unit (Landscape) and the continuity derived from the different ecosystem structures present in that landscape (Ecosystem). Landscape and Ecosystem factors have been calculated using the Fire selectivity Index (SI) [80], which is an indicator of land cover types that are positively selected by fire.
  • Landscape (L)
This is defined as the number of different values of the fire Selectivity Index (SI) in a determined spatial unit, determined by the area occupied by each of them [81]. It is measured through the “Shannon-Wiener Index” (H’), which is an indicator of landscape heterogeneity [81]. It is formulated as (Equation (3)):
H = i = 1 S ( P i × P i )
where:
S = number of different SI values.
Pi = occupied proportion of each SI with respect to the total surface.
The existence of a mosaic with different SI values (higher values of H’) is envisaged to lead to lower emission risk [79]. The landscape-assigned values in terms of the emission risk are classified according to Table A4.
  • Ecosystem (E)
This indicates the horizontal continuity of the different land uses, especially for those that are positively selected by fire, as a factor that facilitates the propagation of a fire. It is measured through the “Aggregation Index” (CONTAG), which expresses the degree of aggregation of the elements of a landscape [82]. It is formulated as (Equation (4)):
CONTAG = [1 + ∑∑[[(Pi)[gik/∑gik]] × [ln(Pi)[gik/∑gik]]]/2 × ln(m)]
where
m: The number of different Sis.
Pi: The proportion of the total surface occupied by a determined SI.
gik: The number of adjacencies between spot types (SI) i and k.
The bigger and more aggregated the SI patches are (lower values of CONTAG), the more horizontal continuity and emission risk are envisaged [72]. Table A5 shows the classification of emission risk according to the resulting ecosystem factor.

2.2.6. Fuel Model (FM)

The classification of different types of forest structures explains the behaviour of fire during a forest fire based on the length of the flame and the speed of propagation [83]. It was measured using the “Vegetation Index”, which expresses the risk derived from the fuel model types (see Table A6). These risk values were obtained following the methodology in Appendix B and are dependent on the site conditions. The higher the vegetation risk, the higher the emission risk [84], so a higher value is assigned as a consequence.

2.2.7. Climatology (C)

Fire behaviour depends, to a large extent, on atmospheric weather elements [84]. However, given that the present model is intended to be used in fire prevention plans, which are drawn up for a large period of time and have a predictive nature, it is not possible to consider variables as changing as the weather [84]. For this reason, climate variables will be included in their place, which are much more stable for long periods of time.
It is assessed based on several indexes proposed by [85] in the Worldwide Bioclimatic Classification System. The driest and warmest bioclimatic regions are envisaged to be the most dangerous and the least cold and humid subtypes [84]. To do so, three indexes are calculated:
  • Thermicity index (Itc)
Thermicity compensated index (for latitudes > 23°) [85] (Equation (5)):
If   Am < 8 ,   Itc = It 80 + 10   Am If   18 > Am > 8 ,   Itc = It If   Am > 18 ,   Itc = It + C 1 + C 2 + C 3 + C 4
where Am is the annual mean amplitude (Tmax-Tmin: Difference between the average temperature of the hottest and coldest months of the year, °C)
C1 = −(Am − 18) (0 < C1 < 15)
C2 = 1 − (Am − 21) (0 < C2 < 105)
C3 = 2 − (Am − 28) (0 < C3 < 450)
C4 = 3 − (Am − 46) (0 < C4 < 570)
The higher the thermicity index, the higher the emission risk is foreseen to be; consequently, Table A7 and Table A8 show the emission risk value for this factor.
  • Continentality index (Ic)
This expresses the difference between the average temperature of the hottest and coldest months of the year, measured in °C. According to [85] (Equation (6)):
Ic = −Tmax − Tmin
The higher the continentality (more extreme temperatures), the higher the emission risk envisaged as a consequence, as shown in Table A9.
  • Ombrothermic index (Io)
This relates to precipitation and temperature according to [85] (Equation (7)):
Io = Pp/Tp
where:
Pp: Positive precipitation (mm) (∑ of months with ti > 0 °C).
Tp: Positive temperature (°C) (∑ of months with ti > 0 °C).
Table A10 shows that the higher the ombrothermic index is (humid climates), the less emission risk is envisaged.

2.2.8. Topography

This describes the relief or morphology of the terrain. There are several factors to consider regarding wildfires, such as the slope, altitude, orientation and relief. Nevertheless, only the slope and orientation were assessed since they have a decisive influence on the behaviour of the fire [84].
  • Slope (S)
The slope is the topographic factor with the greatest influence on the speed of fire propagation [84]. Fire moves faster uphill. In fact, for every 10-degree increase in the slope, a fire will double in speed [86,87]. This is because the slope provides a similar effect to the wind, effectively laying the flames down into the slope and pre-heating the vegetation, allowing it to ignite more quickly [88].
Table A11 shows that the higher the slope (it accelerates the propagation of fire), the more emission risk is envisaged [84].
  • Orientation (O)
This is the geographical direction that the slope faces. The orientation of a slope influences a fire’s behaviour in several ways. Northern and western orientations receive more direct heat from the sun, drying both the soil and vegetation more than on southern or eastern slopes. Fuels are therefore usually drier and less dense on northern and western slopes than fuels on slopes with a different orientation [88]. In contrast, in the Northern hemisphere, Southern and Eastern orientations have drier vegetation and soils due to the more direct heat from the sun.
The Campbell Prediction System’s Flammability Card [89] illustrates the potential fuel temperature variations in the five primary orientations: N, S, E, W and flats. Southern orientations are those with higher fuel temperatures and flammability, followed by the west, east and northern ones. The higher the fuel temperatures and flammability, the more emission risk is envisaged [89]. Table A12 shows the values for the orientation factor.

2.3. Model Build-Up

Figure 1 presents a general scheme of the emission risk model. In the following paragraphs, the model is described and the component variables’ calculation is explained in detail.

2.3.1. Hazard (H)

The values of fire frequency calculated according to [64] were classified into five ranks for use as part of the Emission Risk formula adopting the values shown in Table A1.

2.3.2. Damage (D)

To calculate the damage factor, it is necessary to previously estimate the value at risk and the emission vulnerability. The final damage is the percentage of carbon that has been released by the fire event, and it depends on the ecosystem’s emission vulnerability.
  • Value at risk (VAR)
Carbon stock factor values have been classified into five ranks according to the percentile methodology (quintiles) to determine the classification intervals in the region.
Value at risk is formulated as Equation (8):
VAR = Vegetation Carbon Stock + Soil Organic Carbon Stock
  • Emission Vulnerability (EV)
EV is calculated according to Equation (9)
EV = Potential Fire Severity ÷ Potential Fire Resistance
-
Potential Fire Severity (FS): The term fire severity has traditionally been understood as the degree to which a site has been altered or disturbed by fire [90]. Nevertheless, most of the definitions do not consider operational metrics to assess it. According to [90], most empirical studies that have attempted to measure fire severity have had a common basis that focuses on the loss or decomposition of organic matter, both aboveground and belowground. Aboveground metrics are generally indicators of biomass loss [90,91]. On the other hand, soil characteristics include the loss of the litter and duff layers and ash characteristics, all of which reflect, to varying degrees, the level of organic matter consumed [92]. Additionally, the weather conditions and topography are the other two factors that contribute to the fire behaviour triangle that will be assessed [93]. Nevertheless, this model is intended to be used as a prediction of fire severity in future wildfire events, so other biotic and abiotic variables not related to measurements after a fire event will be considered [84].
Consequently, potential fire severity is calculated as follows (Equation (10)):
FS = a × Fuel Model + b × Climatology + c × Slope + d × Orientation
-
Potential Fire Resistance (FR): This represents the opposition to the propagation or spread of the fire through the forest biomass [56]. It depends on the vertical and horizontal continuity of the forest biomass, which is an indicator of its combustibility [72]. At the same time, combustibility depends on several factors such as the existence of preventive silviculture and infrastructures, heterogeneity of the landscape and extinction capability of the firefighting agents and media, but also meteorological conditions, moisture content of the biomass, terrain slope and orientation, etc. Nevertheless, only the first three variables were considered in the model, integrated into a formula as follows (Equation (11)):
FR = e × PS + f × PI + g × SH
In both equations (FS and FR), a–g coefficients are weighting parameters determined by a panel of 26 experts on silviculture, fire prevention and firefighting from academia, public and private institutions from Portugal, Spain and France, through the Analytic Hierarchy Process (AHP) [94,95].

2.3.3. Emission Risk Assessment

Applying Equation (1) with values from Hazard and Damage factors, the GHG Emission Risk of a territory ranges from 1 to 3, with 1 being a low risk of GHG emissions, 2 being moderate and 3 being a high emission risk from a forest fire. Equation (12) represents the complete GHG emission Risk (GER) model, and Equations (13) and (14) represent the development of the factors into their composing variables.
Equation (15) represents an adaptation of Equation (13) in a region when cartography on Preventive silviculture factors (PS) could be not obtained.
GER = H × D × EV
GER = H × (VAR) × (0.43FM + 0.36C + 0.15S + 0.06O)/(0.27PS + 0.11PI + 0.61SH)
GER = FF × (VCS + SOC) × (0.43FM + 0.36(Ict × Ic × Io) + 0.15S + 0.06O)/(0.27 ((PB + VT)/2) + 0.11(((FN + RN) +
WSN)/2) + 0.61 (L × E))
GER = FF × (VCS + SOC) × (0.43FM + 0.36 (Ict × Ic × Io) + 0.15S + 0.06O)/(0.17 (((FN + RN) + WSN)/2) + 0.83 (L × E))
After applying the final equations (Equations (14) or (15)), the GHG Emission Risk of a territory ranged from 1 to 5, with 1 representing the least risk of GHG emissions and 5 representing the riskiest areas if a fire event occurs. These values were reclassified into three groups to differentiate the zones with low risk, moderate risk and high risk. This classification used the mean ( x ¯ ) and standard deviation (SD) according to Table A12.
Note that this indicator is only valid for comparing intraregional risk values but not for interregional comparisons. The values from one region are not comparable to another, as the data sources and, consequently, the input data might be different. This GER indicator is useful and helpful in terms of identifying and locating the areas within a territory with a higher GHG emission risk in order to establish them as priority action areas and support decision making at all levels.

2.4. Chelva Forest District as Case Study

The Chelva Forest District is located in the northwestern part of the province of Valencia. It covers 183,497 ha, of which 134,121 ha is forest land. From the total area, 108,608 ha are under public management, 93.8% of which are Public Utility Forests, and the rest is forests under consortium or agreement.
Aleppo pine (Pinus halepensis) is the dominant tree species occupying 38.5% of the forest district according to [96]. It was greatly favoured by the reforestation of this area from the 1950s to the 1970s and its high capacity to grow in abandoned fields. Shrubland occupies 36% of the surface. The Chelva forest district falls within the climatic domains of mesomediterranean shrubland (Bupleuro rigidi-Querceto rotundifoliae sigmetum) and its degradation stages with the presence of kermes oaks (Rhamno lycioidi-Querceto cocciferae sigmetum) [97]. Agricultural land is mainly dominated by almond trees, olive groves, cereals, walnut trees and truffle oaks.
Massive, hard limestone and dolomite outcrop on 67% of the surface, while 21% of the surface outcrops limestone and marl in alternation or limestone and clay in alternation. Soils developed from hard rocks (limestone and dolomite) tend to be very shallow (an effective depth of less than 30 cm), stony and often with few carbonates or completely decarbonated.
Orographically, the region is a very rugged territory (average altitude of 500 metres) with a minimum of 400 m in the extreme south and 1400 m in the northwest, the central axis of which is the course of the Turia river, forming wide valleys. Mountainous reliefs dominate with moderate slopes (15–30%) in 43% of the area, followed by gentle slopes (<15%) in 34% of the area and 18% steep slopes (>30%).
2012 was one of the worst years in terms of forest fires for the Valencian Community. A total of 57,500 forest hectares were affected by 502 fires, 8 of which exceeded 100 ha. The Andilla fire was the second biggest fire in all Valencian communities. It started on 26 June 2012 and was under control seven days later, during which it burned 94% of forest land (approximately 19,688 ha) and 6% of agricultural land (1256).
In addition, other big fires occurred in a part of the same forest district, namely, the Chelva fire (started 1st of June 2012 and burned 670 ha) and the Chulilla fire (23 September 2012, which burned 7096 ha). Figure 2 shows its location within the Chelva forest district.
All the variables were studied and calculated to obtain hazard, emission vulnerability, damage and, finally, GHG emission risk values for Pinus halepensis forest lands, as the most representative forest ecosystem. Data sources can be consulted in the Supplementary Material.
The developed model was applied in the Chelva forest district before the big forest fires occurred in 2012 and in 2020 to assess GHG emission risks. To interpret the resulting values, several scale-ups of the areas affected by big forest fires in 2012 (Andilla fire, Chelva fire and Chulilla fire) were analysed and discussed, and the 2012 emission risk level was compared with the 2020 risk level in the same areas.

3. Results and Discussion

Table A13 shows the average risk values resulting from the application of Equations (13) or (14) to the Valencia study area. Figure 3 and Figure 4 show the resulting GHG emission risk in the Chelva forest district in 2012 and in 2020, respectively.
The results are subscribed to Pinus halepensis forest ecosystems. In an overview, it can be observed that GHG emission risk values in 2012 were lower in the central and southern parts of the study area, with higher values distributed mainly in the northeast and northwest. While lower values increased to the moderate emission risk level in these areas, in the Northeast, the emission risk has dropped from high to moderate due to the Andilla fire impacts and the Chelva fire impacts in the Southwest.
The carbon stock loss and consequent fuel model change after the Andilla fire affected the potential fire severity and spatial heterogeneity variables. Moreover, this affected the landscape’s potential resistance by increasing fuel discontinuity. This shows that a fire event has a high impact on future emission risk levels when there is no substantial modification of the rest of the influencing factors in the affected areas.

3.1. GHG Emission Risk Cartography

Figure 3a shows the result of the implementation of the GER on Pinus halepensis forest ecosystems of the Chelva forest district in Valencia (Spain) in 2012 before the big fire events that took place in June, July and August of 2012 (Andilla, Chelva and Chulilla fires, respectively).
Overlapping pictures 2 and 3a, it can be observed that the upper-right area where the Andilla fire later took place had a marked GHG emission risk, as well as Chelva lake surroundings where the Chelva fire took place (middle-left area). The case of the Chulilla fire area is not so easily observable due to the restricted Pinus halepensis forest ecosystems areas in that part (bottom right), but it can be appreciated that medium and high GHG emission risk also prevailed. This indicates that, despite the fire ignition point, the consequent fire spread or real emissions cannot be predicted by the GHG emission risk level. It shows the riskiest areas in terms of GHG emissions in the case of a fire event and how emission risk changes with the modification of model variables’ values over time, as can be observed in Figure 3b.
Figure 4 shows other areas that were used to interpret and validate the model results in detail with the scaling up of these areas and their GHG emission risks in 2012 and 2020.

3.1.1. Scale Up of GHG Emission Risk in Andilla Fire Burned Area (2012)

Figure 5a presents more in detail on the Andilla forest fire area previous to the fire in 2012 and its emission risk (Figure 5b), while Figure 5c shows the more recent state of risk (2020).
It can be observed that the areas with higher biomass density (Figure 6a) and therefore carbon stock together with steeper and higher continuous landscape areas showed high GHG emission risk values, while roads and firebreaks reduced this risk.
Later, in the 2020s, as can be observed in Figure 6c, the relative GHG emission risk decreased in the burned areas while it was maintained or increased in those that remained unburned. This may be due to the decrease in carbon stock, a change in fuel model types after the fire and the lower continuity of fuels in the burned area. In the case of the unburned areas, carbon stocks have increased during these ten years and have reached more emissive stages in the fuel models.

3.1.2. Scape Up of GHG Emission Risk in the Burned Area by Chelva Fire (2012)

Figure 6a shows the orthophoto image for 2010 before the Chelva fire (2012) and Figure 6b presents the GHG emission risk assigned for that area at that time, before the fire event.
The Northeast area with a higher emission risk was burned by the Chelva fire. Again, it can be observed that areas with an apparent higher amount of biomass had higher emission risks and preventive infrastructure such as the firebreak was decreasing the risk of the area.
In 2020, GHG emissions risk is shown in Figure 6c. It can be observed that the burned areas have decreased their emission risks.

3.1.3. Scale Up of GHG Emission Risk in Chulilla Fire Burned Area (2012)

Figure 7a shows the orthophoto image for 2010 before the Chulilla fire (2012) and Figure 7b shows the emission risk assigned for that area before the fire event.
Using Figure 7b,c to compare the GER results, it can be observed that the Chulilla-fire-burned area reduced the carbon stock. Moreover, the decrease in GHG emission risk was due to a change in spatial heterogeneity because of the higher diversity of soil land use than in 2020.

3.1.4. Scale Up of GHG Emission Risk in La Sazadilla

Figure 8a shows the orthophoto image for 2010 in the surroundings of La Sazadilla and Figure 8b presents the GHG emission risk assigned for that area. In the centre of Figure 8a, a high accumulation of carbon can be appreciated in forest ecosystems. Figure 8b shows that this area did not exhibit a high risk as it is diminished by the influence of water points for firefighting that create a risk reduction buffer around them (appreciable as a crown in Figure 8b in the bottom-right corner and the centre of the figure). Moreover, an extended low-risk area (green area) of young forest stands can be observed with relatively low carbon stocks. Forest heterogeneity and discontinuity can be also observed in the centre bottom area.
Nevertheless, Figure 8c, which is the most recent GHG emission risk map, shows a decrease in the reduction effect of the water point due to the increase in carbon stock, with the rest of the variables’ values being equal after ten years. Higher-risk fuel model types influenced this general increase in GHG emission risk in this area.

3.1.5. Scale Up of GHG Emission Risk in Aras de los Olmos

Again, Figure 9a shows the orthophoto image for 2010 in the surroundings of Aras de los Olmos and Figure 9b presents the emission risk assigned for that area at that time.
A general increase in risk can be observed between 2012 (Figure 9b) and 2020 (Figure 9c). This is due to the increase in carbon stock, fuel model change and spatial heterogeneity reduction due to the increase in the abandonment of agricultural terraces.

3.2. Framework, Restrictions and Future Research Envisaged

This model has been developed with the objective of providing an assessment of priority areas for strategic long-term planning to support decision making for preventing high emissions in the atmosphere from forest fires.
The GHG emission risk of a landscape or region is the relative value of a determined region that cannot be compared with other regions in terms of the amount of emission as it indicates an intrinsic regional risk level according to regional conditions and data availability.
Taking this into consideration, there are many variables affecting the risk of emissions that could not be included in the developed model as there are no regional cartographies of these variables or they need to be elaborated on; for instance, updated fuel models or preventive silviculture (not only treated areas but also treatment intensity). On the other hand, the quantification of carbon stocks highly depends on the data available, and National Forest Inventories (NFI) are too time spaced (more than 16 years have already passed since the last NFI in the Valencian Community) to rely on them. Therefore, new approaches such as remote sensing quantification of lidar need to be used to update data on carbon dynamics.
The potential fire severity is a dynamic variable that highly affects the GHG release by a fire and is highly dependent on the meteorological window (temperature, humidity and wind). It has not been contemplated in the scope of the present development and would need to be included for more tactical planning or firefighting use of the model.
On the other hand, this model has been restricted to the assessment of forest ecosystems without including a very important part of the landscape affecting fire behaviour, which is agricultural lands. Therefore, this model should be extended to all territories in a region in order to obtain more consistent results.
Deeper knowledge of the contribution of ecosystemic, landscape, climatological and fire behaviour variables in terms of emissions is needed to improve the results obtained by the present research. Climatic projections should be integrated into the model for envisaged future assessments and meteorology for operative risk level forecasting.

4. Conclusions

A novel methodology has been developed as a basis for the qualitative spatial assessment of forest fire GHG emissions risk. It provides three levels of GHG emissions risk (low, moderate and high) adapted to intrinsic regional conditions. Moderate- and high-emission-risk areas can be considered hotspots to focus forest fire preventive measures to highly reduce emissions in the case of a forest fire event.
This model has been implemented in the Chelva forest district in the Valencia Region (Spain) as a pilot case for 2012 and 2020. The results show high correspondence between a high emission risk level and the areas burned by big fires in 2012.
It provides a spatial assessment of hotspots of GHG emissions risk to implement an action protocol for reducing GHG emissions risk in the framework of the REMAS project for decision-making.
A detailed sensitivity analysis at the regional level allowed us to observe and quantify the effects of forest and landscape management on variables such as carbon stock, fuel models, preventive silviculture, preventive infrastructure and landscape spatial heterogeneity that affect potential fire severity and potential fire resistance and therefore are able to reduce the emission risk level in a region. This is key information for decision making and can be provided by the present model.
The inclusion and consideration of GHG forest fire emission risk in strategic and tactical regional risk and forest planning is a milestone aligned with the latest international objectives and commitments to emission reductions that will not only prevent high releases of GHG emissions but also increase resilience and reduce fire severity linked to megafire events. Moreover, this model contributes to measuring the objectives of the European Commission for 2035 [29] of reducing the risks for forests in a context of uncertainty due to climate change.
Being aware of the vulnerability of carbon stocks in forest ecosystems and the enormous natural capital of rural communities and the value at risk that represent the case of a forest fire event is the first step in holistic and consistent climate change mitigation action that not only considers anthropic GHG sources, but also possible GHG emissions caused by forest fires that can be actively prevented, managed and reduced.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fire6010008/s1.

Author Contributions

Conceptualization, V.L.-A., J.-V.O.-V., F.C.R. and E.L.-S. (Edgar Lorenzo-Sáez); Methodology, V.L.-A., C.Y.-H., Y.B., I.M.-D., F.C.R., S.J. and E.G.-G.; Software, H.V.d.B., M.G.-F., Y.B., E.L.-S. (Eduardo López-Senespleda) and T.P.; Validation, V.L.-A., H.V.d.B., M.G.-F., Y.B., E.L.-S. (Eduardo López-Senespleda), T.P., S.J., E.C.-V. and E.L.-S. (Edgar Lorenzo-Sáez); Formal analysis, V.L.-A., C.Y.-H., H.V.d.B., Y.B., E.L.-S. (Eduardo López-Senespleda), T.P. and E.C.-V.; Investigation, E.C.-V. and E.G.-G.; Data curation, H.V.d.B., E.L.-S. (Eduardo López-Senespleda), T.P., S.J. and E.C.-V.; Writing—original draft, V.L.-A., C.Y.-H. and H.V.d.B.; Writing—review & editing, V.L.-A., Y.B., V.A., M.M.-M., R.R.-P., E.G.-G., R.A.-A. and E.L.-S. (Edgar Lorenzo-Sáez); Funding acquisition, V.L.-A., J.-V.O.-V. and R.A.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financed by the REMAS project (co-financed by the Interreg Sudoe Programme through the European Regional Development Fund (ERDF), grant number “SOE3/P4/E0954”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors want to thank all the experts that voluntarily participated in model variable weighting and shared their wide knowledge and experience with the authors, namely: R Delgado Artés (UPV), JL Soriano, F Navarro, JR García, A. Cervera Montero and JL Soriano Sancho, forest fire analysts (VAERSA, Generalitat Valenciana GV), David Bordes Ortolà (Agència Valenciana de Seguretat i Resposta a les Emergències, GV) J Caamaño (Pau Costa Foundation), N López and M Aguilar (Junta de Comunidades de Castilla-La Mancha), J Madrigal (ICIFOR-INIA, CSIC), Bruno Moreira (CSIC, UV, GV), A López and E Hernández (Ministerio para la Transición Ecológica y el Reto Demográfico) and Cédric Barlet (DFCI Aquitaine). The authors want to especially thank Raul Quílez, former head of the Forest Firefighters from Valencian Province and Tecnosylva researcher and Miguel Ángel Botella and all forest fire analysts (VAERSA, GV) who contributed with their deep and extensive knowledge applied to the conceptualisation and improvement of the model. Finally, the authors want to thank Diego Marin and Mario Romero, director and service head, respectively, of the Valencian Fire Prevention Service from the Conselleria de Consellería de Agricultura, Desarrollo Rural, Emergencia Climática y Transición Ecológica for their support and guidance throughout the whole REMAS project development and Jorge Victoria who provided great support and guidance in the conceptualisation of the model.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Fire frequency assigned values [84].
Table A1. Fire frequency assigned values [84].
Fire Frequency (Fi)RatingValue
<0.2Very low0.8
0.2–0.49Low0.9
0.5–1.99Moderate1
2–3.99High1.1
4–5.99Severe1.2
>6Extreme1.3
Table A2. Vegetation and soil carbon stock assigned values for the Value at Risk rating.
Table A2. Vegetation and soil carbon stock assigned values for the Value at Risk rating.
VAR (Mg C/ha)Value
0–201
21–402
41–603
60–804
>805
Table A3. Preventive silviculture and infrastructure factors’ assigned values.
Table A3. Preventive silviculture and infrastructure factors’ assigned values.
PSValue
Non-existent0
Existent1
Table A4. Landscape assigned values.
Table A4. Landscape assigned values.
Shannon-Wiener Index (H’)Rating (Heterogeneity)Value
0–1Very low5/5
1–2Low4/5
2–3Moderate3/5
3–4High2/5
4–5Very high1/5
Table A5. Ecosystem assigned values.
Table A5. Ecosystem assigned values.
Shannon-Wiener Index (H’)Rating (Continuity)Value
0–20Very high1/5
20–40High2/5
40–60Moderate3/5
60–80Low4/5
80–100Very low5/5
Table A6. Vegetation index assigned values.
Table A6. Vegetation index assigned values.
FM [83]ValueFM [98]Value
TL1-TL41/101, 2, 81/5
TU1, TL5, TL7, SH12/10
GR1, GS1, SH2, SH33/105, 62/5
GR2, GR3, GS2, SH44/10
GS3, SH6, SH8, GR45/103, 7, 93/5
GR5, GR6, TU2, TL66/10
SB1, SH7, TL8, GS47/1010, 114/5
TU3, TU4, TL9, SB28/10
SH5, SB3, SB4, TU59/104, 12, 135/5
SH9, GR7-GR910/10
Table A7. Thermicity Index for Mediterraneas cases assigned values [85].
Table A7. Thermicity Index for Mediterraneas cases assigned values [85].
Thermicity index (Itc)ThermotypeValue
Itc > 120Itc < 120, Tp 1
-1–190Upper Cryoromediterranean1/12
-191–450Lower Cryoromediterranean2/12
-451–675Upper Oromediterranean3/12
-676–900Lower Oromediterranean4/12
120–150901–1200Upper Supramediterranean5/12
150–2201201–1500Lower Supramediterranean6/12
220–2851501–1825Upper Mesomediterranean7/12
285–3501826–2150Lower Mesomediterranean8/12
350–4002151–2300Upper Termomediterranean9/12
400–4502301–2450Lower Termomediterranean10/12
450–5152451–2650Upper Inframediterranean11/12
515–580>2650Lower Inframediterranean12/12
1 At any latitude, when the thermic index (Itc) is less than 120, or when the continental index (Ic) is equal to or greater than 21, the value of the annual positive temperature (Tp) is used to calculate the thermotype which represents the sum in tenths of degrees centigrade of the monthly mean temperatures (Ti) of the months with an average temperature above 0 °C (Tp = sum of Ti> = 0°).
Table A8. Thermicity Index for Temperate cases assigned values [85].
Table A8. Thermicity Index for Temperate cases assigned values [85].
Thermicity Index (Itc)ThermotypeValue
Itc > 120Itc < 120
-1–190Upper Cryorotemperate1/11
-191–380Lower Cryorotemperate2/11
-381–590Upper Orotemperate3/11
-591–800Lower Orotemperate4/11
-801–1100Upper Supratemperate5/11
120–1901101–1400Lower Supratemperate6/11
190–2401401–1700Upper Mesotemperate7/11
240–2901701–2000Lower Mesotemperate8/11
290–3502001–2175Upper Termotemperate9/11
350–4102176–2350Lower Termotemperate10/11
>410>2351Infra Template11/11
Table A9. Continentality Index assigned values [85].
Table A9. Continentality Index assigned values [85].
Continentality Index (Ic)Continentality TypesValue
0–10OceanicHyperoceanic1/6
10–15Euoceanic2/6
15–21Semioceanic3/6
21–27ContinentalSemicontinental4/6
27–46Eucontinental5/6
46–65Hypercontinental6/6
Table A10. Ombrothermic Index assigned values [85].
Table A10. Ombrothermic Index assigned values [85].
Ombrothermic Index (Io)RatingValue
<0.1UltrahyperaridUltrahyperarid16/16
0.1–0.2Hyper aridLower hyper arid15/16
0.2–0.3Upper hyper arid14/16
0.3–0.6AridLower arid13/16
0.6–1.0Upper arid12/16
1.0–1.5Semi-aridLower semi-arid11/16
1.5–2.0Upper semi-arid10/16
2.0–2.8DryLower dry9/16
2.8–3.6Upper dry8/16
3.6–4.8SubhumidLower subhumid7/16
4.8–6.0Upper subhumid6/16
6.0–9.0DampLower damp5/16
9.0–12.0Upper damp4/16
12.0–18.0HyperhumidLower hyperhumid3/16
18.0–24.0Upper hyperhumid2/16
>24.0UltrahyperhumidUltrahyperhumid1/16
Table A11. Slope assigned values.
Table A11. Slope assigned values.
Slope RangeValue
0.0–15.3%1/4
15.3–26.4%2/4
26.4–39.7%3/4
>35.7%4/4
Table A12. Orientation assigned values.
Table A12. Orientation assigned values.
Orientation RangeRatingValue
Flat areas-1/4
315–45°N1/4
45–135°E2/4
225–315°W3/4
135–225°S4/4
Table A13. GHG Emission Risk assigned values.
Table A13. GHG Emission Risk assigned values.
RangeRatingValue
0 − ( x _ S D ) Low risk1
( x _ S D ) ( x _ + 1 S D ) Moderate risk2
( x _ + S D ) − 5High risk3

Appendix B

In this appendix the methodology for rating fuel model risk is described. This methodology has been provided by fire analysts from VAERSA (public company dependent on the Conselleria d’Agricultura, Desenvolupament Rural, Emergència Climàtica i Transició Ecològica from Valencian Government, Generalitat Valenciana).
The methodology consists of comparing the different fuel models’ fire behaviour characteristics in order to assess the inherent danger of each of them. “BehavePlus (5.0.5)” software has been used for simulating fire severity of each fuel model category.
Concrete framework conditions have been established according to most common large forest fires standard conditions in Valencian Region for the simulations: Fuel Moisture 1 h = 3; 10 h = 5; 100 h = 7; Live Herbaceous = 70; Live Woody = 30.
These conditions have been applied to different fuel models, concretely [83,98], classifications.
For the simulation, slope steepness and wind speed variables will vary, establishing 4 scenarios:
  • plain surface and 30 km/h of wind speed (1);
  • plain surface and 50 km/h of wind speed (2);
  • 30% of terrain slope and 30 km/h of wind speed (3);
  • 50% of terrain slope and 50 km/h of wind speed (4).
Additionally, simulations have been run for both, sheltered and unsheltered forests. When simulating sheltered forest types, overstory conditions were: 70% of canopy cover and 12 m of canopy height 1 m of canopy base height (typical conditions for Pinus halepensis species-Valencia’s regional case), and 0.2 kg/m3 of canopy bulk density for Pinus halepensis according to [99,100]. With all these conditions and scenarios, surface forest fires have been simulated.
To assess the danger of a fuel model and to compare it among them, the most interesting parameter is the fireline intensity. According to [101], fireline intensity is defined as the heat energy release per unit time from a one-foot (one-meter) wide section of the fuel bed extending from the front to the rear of the flaming zone. This variable is a function of rate of spread and heat per unit area and is directly related to flame length. Consequently, it represents the potential biomass quantity that is susceptible to be combusted.
In case that the critical surface flame length value is overcome in a simulation, it is supposed that the fire will make a transition to the crowns. If that situation occurs, the heat generated by the crown combustion is added to the heat per unit area to calculate the total fireline intensity. If not, only heat per unit area has been considered.
Finally, the fireline intensity values for both Rothermel and Scott and Burgan have been divided into five and ten percentiles, respectively. According to these percentiles, the risk or danger has been ranged from 1 to 5, being all the fuel models included in percentile 20 or percentile 10 classified as lowest risk, for [83,98] classifications, respectively.

References

  1. Sutherland, E.R.; Make, B.J.; Vedal, S.; Zhang, L.; Dutton, S.J.; Murphy, J.R.; Silkoff, P.E. Wildfire Smoke and Respiratory Symptoms in Patients with Chronic Obstructive Pulmonary Disease. J. Allergy Clin. Immunol. 2005, 115, 420–422. [Google Scholar] [CrossRef] [PubMed]
  2. Giorgi, F. Climate Change Hot–Spots. Geophys Res. Lett. 2006, 33, 101029. [Google Scholar] [CrossRef]
  3. San-Miguel-Ayanz, J.; Moreno, J.M.; Camia, A. Analysis of Large Fires in European Mediterranean Landscapes: Lessons Learned and Perspectives. For. Ecol. Manag 2013, 294, 11–22. [Google Scholar] [CrossRef]
  4. Moreira, F.; Viedma, O.; Arianoutsou, M.; Curt, T.; Koutsias, N.; Rigolot, E.; Barbati, A.; Corona, P.; Vaz, P.; Xanthopoulos, G.; et al. Landscape Wildfire Interactions in Southern Europe: Implications for Landscape Management. J. Environ. Manag. 2011, 92, 2389–2402. [Google Scholar] [CrossRef]
  5. Caon, L.; Vallejo, V.R.; Ritsema, C.J.; Geissen, V. Effects of Wildfire on Soil Nutrients in Mediterranean Ecosystems. Earth Sci. Rev. 2014, 139, 47–58. [Google Scholar] [CrossRef]
  6. Giovannini, G.; Vallejo, R.; Lucchesi, S.; Bautista, S.; Ciompi, S.; Llovet, J. Effects of Land Use and Eventual Fire on Soil Erodibility in Dry Mediterranean Conditions. For. Ecol. Manag. 2001, 147, 15–23. [Google Scholar] [CrossRef]
  7. Pausas, J.G. Changes in Fire and Climate in the Eastern Iberian Peninsula (Mediterranean Basin). Clim. Chang. 2004, 63, 337–350. [Google Scholar] [CrossRef]
  8. Campo, J.; Andreu, V.; Gimeno-García, E.; González, O.; Rubio, J.L. Occurrence of Soil Erosion after Repeated Experimental Fires in a Mediterranean Environment. Geomorphology 2006, 82, 376–387. [Google Scholar] [CrossRef]
  9. Pausas, J.G.; Llovet, J.; Rodrigo, A.; Vallejo, R. Are Wildfires a Disaster in the Mediterranean Basin? A Review. Int. J. Wildland Fire 2008, 17, 713. [Google Scholar] [CrossRef]
  10. Raison, R.J.; Khanna, P.K.; Jacobsen, K.L.S.; Romanya, J.; Serrasolses, I. Effects of Fire on Forest Nutrient Cycles. In Fire Effects on Soils and Restoration Strategies; CRC Press: Boca Raton, FL, USA, 2009; pp. 225–256. [Google Scholar]
  11. Shakesby, R.A. Post-Wildfire Soil Erosion in the Mediterranean: Review and Future Research Directions. Earth Sci. Rev. 2011, 105, 71–100. [Google Scholar] [CrossRef]
  12. Castro, A.C.M.; Carvalho, J.P.; Meixedo, J.P. A Qualitative Description of Soil Parameters Variation Due to a Prescribed Fire in Portuguese Northwestern Forests Using Fuzzy Boolean Nets—The Case Study of Cabreira Mountain. Geoderma 2012, 191, 89–96. [Google Scholar] [CrossRef]
  13. Crockett, J.; Westerling, A. Greater Temperature and Precipitation Extremes Intensify Western, U.S. Droughts, Wildfire Severity, and Sierra Nevada Tree Mortality. J. Clim. 2018, 31, 341–354. [Google Scholar] [CrossRef]
  14. Holden, Z.A.; Swanson, A.; Luce, C.H.; Jolly, W.M.; Maneta, M.; Oyler, J.W.; Warren, D.A.; Parsons, R.; Affleck, D. Decreasing Fire Season Precipitation Increased Recent Western US Forest Wildfire Activity. Proc. Natl. Acad. Sci. USA 2018, 115, E8349–E8357. [Google Scholar] [CrossRef]
  15. Hoinka, K.; Carvalho, A.; Miranda, A. Regional-Scale Weather Pattern and Wildland Fires in Central Portugal. Int. J. Wildland Fire 2009, 18, 36–49. [Google Scholar] [CrossRef]
  16. Piñol, J.; Terradas, J.; Lloret, F. Climate Warming, Wildfire Hazard, and Wildfire Occurrence in Coastal Eastern Spain. Clim. Chang. 1998, 38, 345–357. [Google Scholar] [CrossRef]
  17. Vilén, T.; Fernandes, P.M. Forest Fires in Mediterranean Countries: CO2 Emissions and Mitigation Possibilities Through Prescribed Burning. Environ. Manag. 2011, 48, 558–567. [Google Scholar] [CrossRef] [PubMed]
  18. Chiriacò, M.V.; Perugini, L.; Cimini, D.; D’Amato, E.; Valentini, R.; Bovio, G.; Corona, P.; Barbati, A. Comparison of Approaches for Reporting Forest Fire-Related Biomass Loss and Greenhouse Gas Emissions in Southern Europe. Int. J. Wildland Fire 2013, 22, 730. [Google Scholar] [CrossRef]
  19. Menon, S.; Hansen, J.; Nazarenko, L.; Luo, Y. Climate Effects of Black Carbon Aerosols in China and India. Science 2002, 297, 2250–2253. [Google Scholar] [CrossRef]
  20. Ramanathan, V.; Carmichael, G. Global and Regional Climate Changes Due to Black Carbon. Nat. Geosci. 2008, 1, 221–227. [Google Scholar] [CrossRef]
  21. Alves, C.A.; Gonçalves, C.; Pio, C.A.; Mirante, F.; Caseiro, A.; Tarelho, L.; Freitas, M.C.; Viegas, D.X. Smoke Emissions from Biomass Burning in a Mediterranean Shrubland. Atmos. Environ. 2010, 44, 3024–3033. [Google Scholar] [CrossRef]
  22. Amiro, B.D.; Barr, A.G.; Barr, J.G.; Black, T.A.; Bracho, R.; Brown, M.; Chen, J.; Clark, K.L.; Davis, K.J.; Desai, A.R.; et al. Ecosystem Carbon Dioxide Fluxes after Disturbance in Forests of North America. J. Geophys. Res. 2010, 115, G00K02. [Google Scholar] [CrossRef]
  23. UNFCCC. Sharm el-Sheikh Implementation Plan. In Proceedings of the 27th Conference of the Parties of the UNFCCC (COP 27), Sharm El-Sheikh, Egypt, 6–18 November 2022. [Google Scholar]
  24. UNFCCC. Paris Agreeement. In Proceedings of the 21st Paris Climate Change Conference COP21, Paris, France, 30 November–12 December 2018. [Google Scholar]
  25. United Nation. Sustainable Development Knowledge Platform. Transforming Our World: The 2030 Agenda for Sustainable Development; United Nation: New York, NY, USA, 2015. [Google Scholar]
  26. European Commission, S.-G. COM/2019/640 Final the European Green Deal. Eur. Comm. 2019, 53, 24. [Google Scholar]
  27. European Parliament. Regulation (Eu) 2018/841 of the European Parliament and of the Council of 30 May 2018 on the Inclusion of Greenhouse Gas Emissions and Removals from Land Use, Land Use Change and Forestry in the 2030 Climate and Energy Framework, and Amending Regulation (EU) No 525/2013 and Decision No 529/2013/EU; European Parliament: Strasbourg, France, 2018. [Google Scholar]
  28. European Parliament; European Council. Decision No 529/2013/EU of 21 May 2013 on Accounting Rules on Greenhouse Gas Emissions and Removals Resulting from Activities Relating to Land Use, Land-Use Change and Forestry and on Information Concerning Actions Relating to Those Activities; European Parliament: Strasbourg, France, 2013. [Google Scholar]
  29. European Commission. Communication from the Commission to the European Parliament, the Council, the European Economic and Social Committee and the Committee of the Regions New EU Forest Strategy for 2030; COM/2021/572 Final; European Commission: Brussels, Belgium, 2021. [Google Scholar]
  30. USEPA. Guidelines For. Ecological Risk Assessment. Report No. EPA/630/R-95/002F; USEPA: Washington, DC, USA, 1998. [Google Scholar]
  31. Lecina-Diaz, J.; Martinez Vilalta, J.; Alvarez, A.; Vayreda, J.; Retana, J. Assessing the Risk of Losing Forest Ecosystem Services Due to Wildfires. Ecosystems 2021, 24, 1–15. [Google Scholar] [CrossRef]
  32. Dokas, I.; Statheropoulos, M.; Karma, S. Integration of Field Chemical Data in Initial Risk Assessment of Forest Fire Smoke. Sci. Total Environ. 2007, 376, 72–85. [Google Scholar] [CrossRef]
  33. 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. Geoinf. 2021, 10, 447. [Google Scholar] [CrossRef]
  34. Sivrikaya, F.; Küçük, Ö. Modeling Forest Fire Risk Based on GIS-Based Analytical Hierarchy Process and Statistical Analysis in Mediterranean Region. Ecol. Inform. 2022, 68, 101537. [Google Scholar] [CrossRef]
  35. You, W.; Lin, L.; Wu, L.; Ji, Z.; Yu, J.; Zhu, J.; Fan, Y.; He, D. Geographical Information System-Based Forest Fire Risk Assessment Integrating National Forest Inventory Data and Analysis of Its Spatiotemporal Variability. Ecol. Indic. 2017, 77, 176–184. [Google Scholar] [CrossRef]
  36. Chuvieco, E.; Aguado, I.; Jurdao, S.; Pettinari, M.L.; Yebra, M.; Salas, J.; Hantson, S.; Riva, J.; Ibarra, P.; Rodrigues, M.; et al. Integrating Geospatial Information into Fire Risk Assessment. Int. J. Wildland Fire 2014, 23, 606–619. [Google Scholar] [CrossRef]
  37. Chuvieco, E.; Allgöwer, B.; Salas, J. Integration of Physical and Human Factors in Fire Danger Assessment. In Wildland Fire Danger Estimation and Mapping: The role of Remote Sensing Data; World Scientific Publishing Company: Singapore, 2003; pp. 197–218. [Google Scholar] [CrossRef]
  38. Van Wagner, C.E. Development and Structure of the Canadian Forest Fire Weather Index System; Forestry Technical Report 35; Canadian Forestry Service: Ottawa, ON, Canada, 1987.
  39. Stocks, B.J.; Lynham, T.; Lawson, B.; Alexander, M.; Wagner, C.; McAlpine, R.; Dube, D. Canadian Forest Fire Danger Rating System: An Overview. For. Chron. 1989, 65, 258–265. [Google Scholar] [CrossRef]
  40. Burgan, R.E.; Bradshaw, L.S. WFAS Requires a Variety of Weather Information. Fire Manag. Notes 1997, 57, 18–21. [Google Scholar]
  41. Keetch, J.J.; Byram, G.M. A Drought Index for Forest Fire Control; Department of Agriculture, Forest Service, Southeastern Forest Experiment Station: Asheville, NC, USA, 1968; 35p.
  42. Chuvieco, E.; Aguado, I.; Yebra, M.; Nieto, H.; Salas, J.; Martín, M.P.; Vilar, L.; Martínez, J.; Martín, S.; Ibarra, P.; et al. Development of a Framework for Fire Risk Assessment Using Remote Sensing and Geographic Information System Technologies. Ecol. Modell 2010, 221, 46–58. [Google Scholar] [CrossRef]
  43. Turco, M.; Llasat, M.-C.; von Hardenberg, J.; Provenzale, A. Climate Change Impacts on Wildfires in a Mediterranean Environment. Clim. Chang. 2014, 125, 369–380. [Google Scholar] [CrossRef]
  44. Karali, A.; Hatzaki, M.; Giannakopoulos, C.; Roussos, A.; Xanthopoulos, G.; Tenentes, V. Sensitivity and Evaluation of Current Fire Risk and Future Projections Due to Climate Change: The Case Study of Greece. Nat. Hazards Earth Syst. Sci. 2014, 14, 143–153. [Google Scholar] [CrossRef]
  45. Lung, T.; Lavalle, C.; Hiederer, R.; Dosio, A.; Bouwer, L.M. A Multi-Hazard Regional Level Impact Assessment for Europe Combining Indicators of Climatic and Non-Climatic Change. Glob. Environ. Chang. 2013, 23, 522–536. [Google Scholar] [CrossRef]
  46. Carvalho, A.C.; Carvalho, A.; Martins, H.; Marques, C.; Rocha, A.; Borrego, C.; Viegas, D.X.; Miranda, A.I. Fire Weather Risk Assessment under Climate Change Using a Dynamical Downscaling Approach. Environ. Model. Softw. 2011, 26, 1123–1133. [Google Scholar] [CrossRef]
  47. Pellizzaro, G.; Ventura, A.; Arca, B.; Arca, A.; Duce, P.; Bacciu, V.; Spano, D. Estimating Effects of Future Climate on Duration of Fire Danger Season in Sardinia. In Proceedings of the VI International Forest Fire Research Conference, Coimbra, Portugal, 15–18 November 2010. [Google Scholar]
  48. Moriondo, M.; Good, P.; Durao, R.; Bindi, M.; Giannakopoulos, C.; Corte-Real, J. Potential Impact of Climate Change on Fire Risk in the Mediterranean Area. Clim. Res. 2006, 31, 85–95. [Google Scholar] [CrossRef]
  49. Martínez, J.; Vega-Garcia, C.; Chuvieco, E. Human-Caused Wildfire Risk Rating for Prevention Planning in Spain. J. Environ. Manag. 2009, 90, 1241–1252. [Google Scholar] [CrossRef] [PubMed]
  50. Botequim, B.; Garcia-Gonzalo, J.; Marques, S.; Ricardo, A.; Borges, J.; Tomé, M.; Oliveira, M. Developing Wildfire Risk Probability Models for Eucalyptus Globulus Stands in Portugal. IForest 2013, 6, 217–227. [Google Scholar] [CrossRef]
  51. Finney, M.A. The Challenge of Quantitative Risk Analysis for Wildland Fire. For. Ecol. Manag. 2005, 211, 97–108. [Google Scholar] [CrossRef]
  52. Verde, J.C.; Zêzere, J.L. Assessment and Validation of Wildfire Susceptibility and Hazard in Portugal. Nat. Hazards Earth Syst. Sci. 2010, 10, 485–497. [Google Scholar] [CrossRef]
  53. Reisinger, A.; Howden, M.; Vera, C.; Garschagen, M.; Hurlbert, M.; Kreibiehl, S.; Mach, K.J.; Mintenbeck, K.; O’Neill, B.; Pathak, M.; et al. The Concept of Risk in the IPCC Sixth Assessment Report: A Summary of Cross-Working Group Discussions; WMO: Geneva, Swizerland, 2020. [Google Scholar]
  54. Rego, F.C.; Colaço, M.D.C.A. Wildfire Risk Analysis. In Encyclopedia of Environmetrics; Wiley: Hoboken, NJ, USA, 2012. [Google Scholar] [CrossRef]
  55. BOE. REAL DECRETO 893/2013, de 15 de Noviembre, Por El Que Se Aprueba La Directriz Básica de Planificación de Protección Civil de Emergencia Por Incendios Forestales; BOE-A-2013-12823. 2013. Available online: https://www.boe.es/buscar/doc.php?id=BOE-A-2013-12823 (accessed on 7 December 2013).
  56. López-Poma, R. Vegetation Traits Modulate Resilience to Fire in Mediterranean Woodlands; University of Alicante: Alicante, Spain, 2014. [Google Scholar]
  57. Dale, V.H.; Joyce, L.A.; McNulty, S.; Neilson, R.P.; Ayres, M.P.; Flannigan, M.D.; Hanson, P.J.; Irland, L.C.; Lugo, A.E.; Peterson, C.J.; et al. Climate Change and Forest Disturbances: Climate Change Can Affect Forests by Altering the Frequency, Intensity, Duration, and Timing of Fire, Drought, Introduced Species, Insect and Pathogen Outbreaks, Hurricanes, Windstorms, Ice Storms, or Landslides. Bioscience 2001, 51, 723–734. [Google Scholar] [CrossRef]
  58. Garcia-Hurtado, E.; Pey, J.; Baeza, M.J.; Carrara, A.; Llovet, J.; Querol, X.; Alastuey, A.; Vallejo, V.R. Carbon Emissions in Mediterranean Shrubland Wildfires: An Experimental Approach. Atmos Environ. 2013, 69, 86–93. [Google Scholar] [CrossRef]
  59. Finkral, A.J.; Evans, A.M. The Effects of a Thinning Treatment on Carbon Stocks in a Northern Arizona Ponderosa Pine Forest. For. Ecol. Manag. 2008, 255, 2743–2750. [Google Scholar] [CrossRef]
  60. Hurteau, M.; North, M. Fuel Treatment Effects on Tree-Based Forest Carbon Storage and Emissions under Modeled Wildfire Scenarios. Front. Ecol. Environ. 2009, 7, 409–414. [Google Scholar] [CrossRef]
  61. North, M.; Hurteau, M.; Innes, J. Fire Suppression and Fuels Treatment Effects on Mixed-Conifer Carbon Stocks and Emissions. Ecol. Appl. 2009, 19, 1385–1396. [Google Scholar] [CrossRef] [PubMed]
  62. Stephens, S.L.; Moghaddas, J.J.; Hartsough, B.R.; Moghaddas, E.E.Y.; Clinton, N.E. Fuel Treatment Effects on Stand-Level Carbon Pools, Treatment-Related Emissions, and Fire Risk in a Sierra Nevada Mixed-Conifer ForestPublication No. 143 of the National Fire and Fire Surrogate Project. Can. J. For. Res. 2009, 39, 1538–1547. [Google Scholar] [CrossRef]
  63. Valenciana, G. Plan Especial Frente al Riesgo de Incendios Forestales de La Comunidad Valenciana. 2017, Volume 159. Available online: https://www.112cv.gva.es/documents/163565706/163566493/PE_Incendios.pdf/d615af2c-8655-4e39-9b0c-2ca4c251c1cf (accessed on 1 February 2017).
  64. VAERSA. Metodología Del Riesgo Estadístico Del Plan de Prevención de Incendios Forestales de La Demarcación de Chelva; VAERSA: València, Spain, 2017. [Google Scholar]
  65. le Noë, J.; Erb, K.-H.; Matej, S.; Magerl, A.; Bhan, M.; Gingrich, S. Altered Growth Conditions More than Reforestation Counteracted Forest Biomass Carbon Emissions 1990–2020. Nat. Commun 2021, 12, 6075. [Google Scholar] [CrossRef]
  66. Vinué-Visús, D.; Ruiz-Peinado, R.; Fuente, D.; Oliver-Villanueva, J.-V.; Coll-Aliaga, E.; Lerma-Arce, V. Biomass Assessment and Carbon Sequestration in Post-Fire Shrublands by Means of Sentinel-2 and Gaussian Processes. Forests 2022, 13, 771. [Google Scholar] [CrossRef]
  67. Lerma-Arce, V.; van den Berg, H.; Oliver-Villanueva, J.V.; Coll-Aliaga, E.P. Cartografía Territorial Del Stock de Carbono En La Comunitat Valenciana. 2020. Available online: https://politicaterritorial.gva.es/documents/20551069/174233262/Cartograf%C3%ADa+Territorial+del+Stock+de+Carbono+en+la+Comunitat+Valenciana.pdf/7e2501f8-2737-426e-80c5-617f33e98f36?t=1627301770832 (accessed on 19 January 2020).
  68. FAO. Measuring and Modelling Soil Carbon Stocks and Stock Changes in Livestock Production Systems: Guidelines for Assessment (Version 1). Livestock Environmental Assessment and Performance (LEAP) Partnership; FAO: Rome, Italy, 2019. [Google Scholar]
  69. Vélez, R. Selvicultura Preventiva de Incendios Forestales; ICONA. 1990. Available online: https://www.researchgate.net/publication/337331889_Selvicultura_preventiva_de_incendios (accessed on 20 January 1990).
  70. Omi, P.N.; Martinson, E.J. Effects of Fuels Treatment on Wildfire Severity. Final Report to the Joint. Fire Science Program Governing Board; Colorado State University: Fort Collins, CO, USA, 2002. [Google Scholar]
  71. Narayan, C.; Fernandes, P.M.; van Brusselen, J.; Schuck, A. Potential for CO2 Emissions Mitigation in Europe through Prescribed Burning in the Context of the Kyoto Protocol. For. Ecol. Manag. 2007, 251, 164–173. [Google Scholar] [CrossRef]
  72. Serrada, R. Apuntes de Selvicultura; Escuela Universitaria de Ingeniería Técnica Forestal; Universidad Politécnica de Madrid. Fundación Conde del Valle Salazar: Madrid, Spain, 2011. [Google Scholar]
  73. Graham, R.T.; Harvey, A.E.; Jain, T.B.; Tonn, J.R. The Effects of Thinning and Similar Stand Treatments on Fire Behavior in Western Forests.; U.S. Department of Agriculture, Forest Service, Pacific Northwest Research Station: Portland, OR, USA, 1999. [CrossRef]
  74. Crecente-Campo, F.; Pommerening, A.; Rodríguez-Soalleiro, R. Impacts of Thinning on Structure, Growth and Risk of Crown Fire in a Pinus Sylvestris L. Plantation in Northern Spain. For. Ecol. Manag. 2009, 257, 1945–1954. [Google Scholar] [CrossRef]
  75. Generalitat Valenciana. Infraestructuras de Prevención. Available online: http://agroambient.gva.es/es/web/prevencion-de-incendios/infraestructuras-de-prevencion (accessed on 1 November 2022).
  76. VAERSA. Manual de Ingeniería. Infraestructuras de Prevención de Incendios Forestales. Norma Técnica de Viales Forestales. Red Eléctrica de España, Ed.; Dirección General de Prevención, Extinción de Incendios y Emergencias Generalitat Valenciana. 2015. Available online: https://agroambient.gva.es/documents/162905929/163206728/NT+Viales+forestales+%289%2C7Mb%29/5996f665-3fd2-46ed-a8e6-99f5e953576f (accessed on 13 February 2015).
  77. VAERSA. Manual de Ingeniería. Infraestructuras de Prevención de Incendios Forestales. Norma Técnica de Áreas Cortafuegos; Red Eléctrica de España, Ed.; Dirección General de Prevención, Extinción de Incendios y Emergencias. Generalitat Valenciana, 2015. Available online: https://agroambient.gva.es/documents/162905929/163206728/NT+Puntos+de+agua+%285%2C6Mb%29/91ee562d-08c1-440b-807e-7ab83ffce0f3 (accessed on 13 February 2015).
  78. VAERSA. Manual de Ingeniería. Infraestructuras de Prevención de Incendios Forestales. Norma Técnica de Puntos de Agua. Red Eléctrica de España, Ed.; Dirección General de Prevención, Extinción de Incendios y Emergencias. Generalitat Valenciana, 2015. Available online: https://agroambient.gva.es/documents/162905929/163206728/NT+Puntos+de+agua+%285%2C6Mb%29/91ee562d-08c1-440b-807e-7ab83ffce0f3 (accessed on 13 February 2015).
  79. Alcasena-Urdíroz, F.J.; Vega-García, C.; Ager, A.A.; Salis, M.; Nauslar, N.J.; Mendizabal, F.J.; Castell, R. Metodología de Evaluación Del Riesgo de Incendios Forestales y Priorización de Tratamientos Multifuncionales En Paisajes Mediterráneos. Cuad. De Investig. Geográfica 2019, 45, 571–600. [Google Scholar] [CrossRef]
  80. Rego, F.C.; Bunting, S.C.; Strand, E.K.; Godinho-Ferreira, P. Applied Landscape Ecology; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  81. European Landscape Character Assessment Initiative; European Landscape Character Areas. Typologies, Cartography and Indicators for the Assessment of Sustainable Landscapes; Wascher, D.M., Ed.; Landscape Europe: Wageningen, The Netherlands, 2005. [Google Scholar]
  82. Vallés Planells, M.C. Definición de Variables Medioambientales Para La Clasificación Jerárquica de Unidades de Paisaje. Aplicación a Casos En La Comunidad Valenciana. Doctoral Dissertation, Universitat Politècnica de València, Valencia, Spain, 2009. [Google Scholar] [CrossRef]
  83. Scott, J.H.; Burgan, R.E. Standard Fire Behavior Fuel Models: A Comprehensive Set for Use with Rothermels Surface Fire Spread Model; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Baltimore, MD, USA, 2005. [CrossRef]
  84. DOGV. ORDEN 30/2017, de 20 de Noviembre, de La Consellera de Agricultura, Medio Ambiente, Cambio Climático y Desarrollo Rural, Por La Que Se Unifican y Aprueban Las Normas Técnicas Para La Redacción de Planes Locales de Prevención de Incendios Forestales (PLPIF). 2017, pp. 44256–44370. Available online: https://dogv.gva.es/datos/2017/11/30/pdf/2017_10697.pdf (accessed on 3 November 2017).
  85. Rivas-Martínez, S. Clasificación Bioclimática de La Tierra. Global Bioclimatics; Centro de Investigaciones Fitosociológicas, Ed.; Universidad Complutense de Madrid: Madrid, Spain, 2004. [Google Scholar]
  86. McArthur, A.G. Fire Behaviour in Eucalypt Forests; Forestry and Timber Bureau: Australia, 1967. [Google Scholar]
  87. Wagner, C.E. van. Effect of Slope on Fires Spreading Downhill. Can. J. For. Res. 1988, 18, 820–822. [Google Scholar] [CrossRef]
  88. Government of South Australia. Fire Behaviour. Department for Environment and Water. Government of South Australia. Available online: https://www.environment.sa.gov.au/topics/fire-management/fire-science/fire-behaviou (accessed on 21 February 2021).
  89. Campbell, D. The Campbell Prediction System. A Wildland Fire Prediction and Communication System. Chapter 3: Solar Preheating. 1991. Available online: http://cps.emxsys.com/ (accessed on 6 September 1991).
  90. Keeley, J.E. Fire Intensity, Fire Severity and Burn Severity: A Brief Review and Suggested Usage. Int. J. Wildland Fire 2009, 18, 116. [Google Scholar] [CrossRef]
  91. Wagner, C.E. van. Height of Crown Scorch in Forest Fires. Can. J. For. Res. 1973, 3, 373–378. [Google Scholar] [CrossRef]
  92. Wells, C.G.; Campbell, R.E.; DeBano, L.F.; Lewis, C.E.; Fredriksen, R.L.; Franklin, E.C.; Froelich, R.C.; Dunn, P.H. Effects of Fire on Soil: A State of Knowledge Review; U.S. Department of Agriculture, Forest Service: Washington, DC, USA, 1979.
  93. National Park Service. Wildland Fire Behaviour Series. Wildland Fire-Learning in Depth. US. Department of the Interior. Available online: https://www.nps.gov/articles/wildland-fire-behavior.htm#:~:text=Topography%20can%20have%20an%20influence,or%20than%20on%20flat%20terrain.&text=These%20topographical%20features%20can%20help,wide%20gap%20of%20open%20space (accessed on 21 February 2021).
  94. Saaty, R.W. The Analytic Hierarchy Process—What It Is and How It Is Used. Math. Model. 1987, 9, 161–176. [Google Scholar] [CrossRef]
  95. Aznar, J.; Guijarro, F. Nuevos Métodos de Valoración. Modelos Multicriterio, 2nd ed.; Editorial Universitat Politècnica de València: Valencia, Spain, 2012. [Google Scholar]
  96. Ministerio de Agricultura, A. y M. A. Tercer Inventario Forestal Español. 2006. Available online: https://www.miteco.gob.es/es/biodiversidad/servicios/banco-datos-naturaleza/informacion-disponible/ifn3.aspx (accessed on 24 November 2006).
  97. Rivas Martínez, S. Memoria Del Mapa de Series de Vegetación de España; Serie Técnica; ICONA; Ministerio Agricultura, Pesca y Alimentación: Madrid, Spain, 1987.
  98. Rothermel, R.C. A Mathematical Model for Predicting Fire Spread in Wildland Fuels; Intermountain Forest and Range Experiment Station: Ogden, UT, USA, 1972; 44p, Pamphlet. Research Paper INT-116. [Google Scholar]
  99. Andersen, H.; McGaughey, R.; Reutebuch, S. Estimating forest canopy fuel parameters using LIDAR data. Remote Sens. Environ. 2005, 94, 441–449. [Google Scholar] [CrossRef]
  100. Mitsopoulos, I.; Dimitrakopoulos, A. Canopy fuel characteristics and potential crown fire behavior in Aleppo pine (Pinus halepensis Mill.) forests. Ann. For. Sci. 2007, 64, 287–299. [Google Scholar] [CrossRef]
  101. Alexander, M.; Cruz, M. Fireline Intensity. In Encyclopedia of Wildfires and Wildland-Urban Interface (WUI) Fires; Springer: Berlin/Heidelberg, Germany, 2019; pp. 1–8. [Google Scholar] [CrossRef]
Figure 1. General scheme of the emission risk model.
Figure 1. General scheme of the emission risk model.
Fire 06 00008 g001
Figure 2. Location of the largest fires in Chelva forest district in 2012.
Figure 2. Location of the largest fires in Chelva forest district in 2012.
Fire 06 00008 g002
Figure 3. (a) GHG emission risk in Chelva forest district in 2012, previous to the fires. (b) GHG emission risk in Chelva’s forest district in 2020.
Figure 3. (a) GHG emission risk in Chelva forest district in 2012, previous to the fires. (b) GHG emission risk in Chelva’s forest district in 2020.
Fire 06 00008 g003
Figure 4. GHG emission risk scale-up areas.
Figure 4. GHG emission risk scale-up areas.
Fire 06 00008 g004
Figure 5. (a) Scale up of Andilla forest fire (2010 orthophoto). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020 (2020 orthophoto).
Figure 5. (a) Scale up of Andilla forest fire (2010 orthophoto). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020 (2020 orthophoto).
Fire 06 00008 g005
Figure 6. (a) Scaling up of Chelva forest fire (2010 orthophoto). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Figure 6. (a) Scaling up of Chelva forest fire (2010 orthophoto). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Fire 06 00008 g006aFire 06 00008 g006b
Figure 7. (a) Scale-up of Chulilla forest fire (2010 orthophoto). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Figure 7. (a) Scale-up of Chulilla forest fire (2010 orthophoto). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Fire 06 00008 g007
Figure 8. (a) Scale-up of La Sazadilla (orthophoto 2010). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Figure 8. (a) Scale-up of La Sazadilla (orthophoto 2010). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Fire 06 00008 g008aFire 06 00008 g008b
Figure 9. (a) Scale-up of Aras de los Olmos (orthophoto 2010). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Figure 9. (a) Scale-up of Aras de los Olmos (orthophoto 2010). (b) GHG emission risk in 2012. (c) GHG emission risk in 2020.
Fire 06 00008 g009aFire 06 00008 g009b
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lerma-Arce, V.; Yagüe-Hurtado, C.; Van den Berg, H.; García-Folgado, M.; Oliver-Villanueva, J.-V.; Benhalima, Y.; Marques-Duarte, I.; Acácio, V.; Rego, F.C.; López-Senespleda, E.; et al. Development of a Model to Estimate the Risk of Emission of Greenhouse Gases from Forest Fires. Fire 2023, 6, 8. https://doi.org/10.3390/fire6010008

AMA Style

Lerma-Arce V, Yagüe-Hurtado C, Van den Berg H, García-Folgado M, Oliver-Villanueva J-V, Benhalima Y, Marques-Duarte I, Acácio V, Rego FC, López-Senespleda E, et al. Development of a Model to Estimate the Risk of Emission of Greenhouse Gases from Forest Fires. Fire. 2023; 6(1):8. https://doi.org/10.3390/fire6010008

Chicago/Turabian Style

Lerma-Arce, Victoria, Celia Yagüe-Hurtado, Helena Van den Berg, Miguel García-Folgado, Jose-Vicente Oliver-Villanueva, Yacine Benhalima, Inês Marques-Duarte, Vanda Acácio, Francisco C. Rego, Eduardo López-Senespleda, and et al. 2023. "Development of a Model to Estimate the Risk of Emission of Greenhouse Gases from Forest Fires" Fire 6, no. 1: 8. https://doi.org/10.3390/fire6010008

APA Style

Lerma-Arce, V., Yagüe-Hurtado, C., Van den Berg, H., García-Folgado, M., Oliver-Villanueva, J. -V., Benhalima, Y., Marques-Duarte, I., Acácio, V., Rego, F. C., López-Senespleda, E., Menéndez-Miguélez, M., Ruiz-Peinado, R., Petillon, T., Jalabert, S., Carbó-Valverde, E., Gimeno-García, E., Aleix-Amurrio, R., & Lorenzo-Sáez, E. (2023). Development of a Model to Estimate the Risk of Emission of Greenhouse Gases from Forest Fires. Fire, 6(1), 8. https://doi.org/10.3390/fire6010008

Article Metrics

Back to TopTop