Next Article in Journal
An Empirical Modelling and Simulation Framework for Fire Events Initiated by Vegetation and Electricity Network Interactions
Previous Article in Journal
Kinetic Study of Pyrolysis of Coniferous Bark Wood and Modified Fir Bark Wood
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Fire
Volume 6
Issue 2
10.3390/fire6020060
fire-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Recent Trends in Fire Regimes and Associated Territorial Features in a Fire-Prone Mediterranean Region

by
Francisco Moreira
1,2,3,
Miguel Leal
4,5,*,
Rafaello Bergonse
5,
Maria João Canadas
4,
Ana Novais
4,
Sandra Oliveira
5,
Paulo Flores Ribeiro
4,
José Luís Zêzere
5 and
José Lima Santos
4
1
CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Instituto Superior de Agronomia, Universidade de Lisboa. Tapada da Ajuda, 1349-017 Lisbon, Portugal
2
CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Campus de Vairão, Universidade do Porto, 4485-661 Vairão, Portugal
3
BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO, Campus de Vairão, 4485-661 Vairão, Portugal
4
Forest Research Centre and Associate Laboratory TERRA, School of Agriculture, Universidade de Lisboa, Tapada da Ajuda, 1349-017 Lisbon, Portugal
5
Centre of Geographical Studies and Associate Laboratory TERRA, Institute of Geography and Spatial Planning, Universidade de Lisboa, Edifício IGOT, Rua Branca Edmée Marques, 1600-276 Lisbon, Portugal
*
Author to whom correspondence should be addressed.
Fire 2023, 6(2), 60; https://doi.org/10.3390/fire6020060
Submission received: 14 December 2022 / Revised: 3 February 2023 / Accepted: 6 February 2023 / Published: 8 February 2023
Browse Figures
Versions Notes

Abstract

:
Fire regimes in Mediterranean countries have been shifting in recent decades, including changes in wildfire size and frequency. We sought to describe changes in fire regimes across two periods (1975–1995 and 1996–2018) in a fire-prone region of central Portugal, explore the relationships between these regimes and territorial features, and check whether these associations persisted across periods. Two independent indicators of fire regimes were determined at parish level: fire incidence and burn concentration. Most parishes presented higher values of both indicators in the second period. Higher values of fire incidence were associated with lower population densities, lower proportions of farmland areas and higher proportions of natural vegetation. Higher levels of burn concentration were associated with smaller areas of farmland and natural vegetation. These associations differed across periods, reflecting contrasting climatic and socio-economic contexts. Keeping 40% of a parish territory covered by farmland was effective to buffer the increased wildfire risks associated with different management and climate contexts. The effectiveness of higher population densities in keeping fire incidence low decreased in the last decades. The results can improve the knowledge on the temporal evolution of fire regimes and their conditioning factors, providing contributions for spatial planning and forest/wildfire management policies.

1. Introduction

Fire regimes can be understood as the average fire conditions (frequency, density, intensity, severity, seasonality, etc.) in a particular area and over a long period [1,2,3]. Past and present fire regimes are crucial for understanding fire spatial and temporal behavior and for estimating future changes [1,4,5,6]. This is particularly important in Mediterranean Europe, one of the most fire-prone regions worldwide due to its climate and vegetation characteristics [7,8].
Fire regimes are affected by land use/land cover (LULC), human activities, climate, and topographical context (e.g., [6,7,9,10,11,12,13,14,15,16,17,18]). Understanding and modeling these factors, in particular how they interact with each other to influence fire regimes, remains a key challenge [6,19,20]. LULC has a crucial role on driving fire regimes, by regulating moisture content and the availability, continuity and flammability of fuel loads [21,22]. In fact, several studies demonstrated that fires are selective regarding land cover (e.g., [22,23,24]), and in the case of the Mediterranean region, scrublands, forest plantations (mainly maritime pine (Pinus pinaster) and eucalyptus (Eucalyptus globulus)), grasslands and woodlands are the most fire-prone, while built-up surfaces, farmland areas and agroforestry systems are less likely to burn [22,23,24,25,26]. The continuity/contiguity of fuel load is also strongly dependent on LULC patterns. A diverse and heterogeneous landscape is capable of modifying the fuel characteristics and conditions in space, hampering wildfire propagation and reducing the fire extent and its capacity to cause damage [13,24,27,28]. On the contrary, a more homogeneous landscape is less fire-resistant and promotes a higher wildfire hazard [22].
Population is another major driver of fire regimes. Fire ignitions are often associated with human presence, agricultural activity, and livestock grazing [7,9,16,20,22,29,30,31,32,33,34]. Intentional or negligent/accidental ignitions frequently result from using fire for agricultural and grazing purposes, i.e., harvest waste removal and brushwood/abandoned land clearance [20,30]. On the other hand, the abandonment of rural areas and consequent population loss has caused the decline of traditional land uses (farmland, silviculture, and livestock), enabling the growth of unmanaged forests and scrublands, promoting fuel accumulation and resulting in more frequent and larger wildfires [7,12,30,35,36,37,38,39].
Finally, the temporal trends of wildfires and how they spread over landscapes are greatly influenced by annual variations in temperature, rainfall and relative humidity, with wildfire occurrence and severity being favored by long-term droughts and heat waves [5,9,10,11,15,40,41,42,43,44].
The effect of LULC and human population on fire regimes is therefore likely to vary across periods in the magnitude of decades, due to temporal changes in climatic factors and socioecological context [6,45]. For example, some authors state that the fire regimes in the western Mediterranean have already shifted from fuel-limited to drought-driven, which led to larger burned areas [7,12,46,47]. However, the temporal changes in fire regime drivers and their relationship with fire incidence have been considered in few studies [1,10,11,37,48], and are not yet fully understood.
Here, we examine how spatial patterns of LULC and population density are associated with local fire regimes in a wildfire-prone region of central Portugal, and how different periods influence these relationships. Portugal stands out among the European countries, recording the highest wildfire density and percentage of burned areas in the last decades [16,29,47,49,50,51,52]. We expected that the effects of LULC changes, rural abandonment, and lack of land management in this region affect the severity of the fire regime across time. To address this, we compared LULC and population density associations to fire regimes across two periods with similar durations (1975–1995 and 1996–2018), ensuring a representative number of years in each period. Environmental and socioeconomic changes are progressive along time, so we expected that this simple division would contrast the average values and trends of features in these two periods. Specifically, we aimed: (1) to characterize and describe changes in fire regimes and their associated territorial features (LULC and population) across periods; and (2) to model how LULC and population were associated with local fire regimes, and to what extent the impacts of these drivers changed between the two periods. The obtained results provide insights on the temporal evolution of fire regimes and their conditioning factors, with implications for spatial planning and forest/wildfire management policies.

2. Materials and Methods

2.1. Study Area

The study area (Figure 1) is located in central Portugal, a highly wildfire-prone region, and covers 15,242 km2 and 568 civil parishes, our unit of analysis that corresponds to the smallest administrative units in the country. The average area of these parishes is 27 km2 (range = 2–172 km2). The climate is Mediterranean (Köppen classes Csa and Csb), with dry and hot/warm summers. Rainfall is mainly concentrated in the autumn and winter months (around 70% of the annual values), although the spring months represent on average 20–25% of the annual values. Elevation ranges between 1993 m a.m.s.l. in the central section of the region and 23 m a.m.s.l. in the south. In 2015, more than 3/4 of the landscape was occupied by forests, in which maritime pine (Pinus pinaster), scrublands and eucalyptus (Eucalyptus globulus) stand out, representing 28%, 20% and 14% of the study area, respectively. Built-up areas are mostly scattered and interspersed with wildland, corresponding only to 3% of the study area. In 2011, ca. 730,000 people lived in the region (48.1 inhabitants/km2), but the population has been declining and aging over recent decades, along with progressive rural abandonment. During the 1975–2018 period, the total burned area reached 2,113,528 ha (1.39 times the study area’s surface).

2.2. Fire Regime Characterization

Fire regimes were characterized from a set of seven variables expressing how fire affected each of the parishes across time and space during each of the two considered periods: 1975–1995 and 1996–2018 (Table 1). These were estimated based on maps of yearly burned area, obtained in vector format from the National Forests Service (Portuguese Institute for the Conservation of Nature and Forests—ICNF) (Lisbon, Portugal). ArcMap 10.7.1 (ESRI Inc.) (Redlands, CA, USA) software was used for all spatial analysis operations.
The total burned percentage describes the accumulated burned area during each period, expressed as percentage of the parish area. It was used as an indicator of the extensiveness of wildfires. The proportions of years in which a parish burned 25%, 50% and 75% of its area were likewise used as complementary indicators of fire extent.
Mean and maximum times burned were used as indicators of the frequency of burning. They were obtained by dividing each parish into 25 m pixels and calculating the mean and the maximum number of times that the pixels within each parish burned during each period.
Finally, we used a variable to describe the temporal distribution of wildfire damage within the fire regime. This was the Gini concentration index (Gini), applied to the yearly burned areas of each parish over each period. This index quantifies (between 0 and 100) the temporal concentration of the total burned area for each parish, allowing the differentiation of parishes where most of the burned area is concentrated in a small number of years (higher values) from those where the burned area is more regularly distributed over time (lower values). It was originally developed in the field of economics [53] and it has been used to measure the concentration of wildfire damage over wildfire events [54] or as the basis for concentration indices with the same purpose [55,56,57].

2.3. Territorial Features

Population density in each parish for each period was calculated from the average of two national population censuses (1981 and 1991 for the first period; 2001 and 2011 for the second period), produced by Statistics Portugal (Table 2).
LULC was estimated considering three main categories (Table 2): farmland (including annual crops, permanent crops, and pastures), tree plantations (including pine and eucalypt plantations), and natural vegetation (including scrublands and oak forests). For each parish, the proportional cover by these three categories was estimated from detailed LULC maps provided by the Portuguese Directorate-General for Territory in 1995 (used to characterize the first period) and 2015 (used to characterize the second period). The date of the first LULC map (1995) is the criterion for the division between periods in 1995/1996. The LULC maps were created through the interpretation of digital orthorectified aerial photos, with a minimum cartographic unit of 1 hectare.
The Castelo Branco rain gauge, located in the southeast of central Portugal, was used to determine the average annual values of temperature and rainfall for the study area. This rain gauge belongs to the Portuguese Institute for Sea and Atmosphere (IPMA) and is the only one whose temperature and rainfall data cover the entire studied period.

2.4. Data Analyses

All analyses were performed using R software [58]. We used principal component analysis (PCA) [59] to summarize the main trends in fire regimes. Each parish was characterized in both periods so that resulting parish coordinates in the PCA could be comparable across time. Differences across periods were visualized using kernel density plots to show the probability density function of coordinate distribution using the ggplot2 package [60]. The significance of the observed differences was tested using linear mixed models with parish ID as a random effect. The package lme4 [61] was used for fitting these models.
The PCA derived two main axes representing two indicators of fire regime. The relationship between these indicators (as response variable) and LULC variables, population density (log-transformed) and period were analyzed using linear mixed models using parish ID as a random effect. Before running the models, collinearity between explanatory variables was assessed using correlation and variance inflation factors (VIF) using the usdm package [62]. The proportion of natural vegetation and proportion of plantations were highly correlated (r = −0.82); therefore we discarded the latter, after which no variable showed collinearity problems, with all VIF smaller than 2 [63]. The two-way interaction between period and each of the other variables was included in an initial full model that was simplified by discarding nonsignificant variables (and interactions), using the step function, until all remaining variables were significant. Predicted values were visualized using the package jtools [64], controlling for the effects of other variables in the model (by keeping them at their mean values or base level for categorical variables). A few parishes (n = 15) did not burn in any of the two periods and were discarded from the analyses.

3. Results

3.1. Fire Regimes

Yearly burned area in the region between 1975 and 2018 averaged 48,000 ha (41,739 ha in the first period and 53,783 ha in the second period), ranging from a minimum of 1866 ha (in 2018) to 383,000 ha (in 2017), one of the worst fire years in the country due to the occurrence of two extreme events (June and October). The three worst fire years occurred in the second period (Figure 2). The sum of the three years with the largest burned area (2003, 2005 and 2017) reached 677,263 ha, which means that almost a third of the total burned area in 44 years was concentrated in just three years.
The PCA yielded two principal components with eigenvalues > 1, which accounted for 81.3% of the total variance (Table 3). PC1 (hereafter designated “fire incidence”) represented a gradient of increasing total burned area and fire frequency, with increasing parish PC scores representing increased total accumulated burned area, higher average number of years a pixel was burned, and a higher proportion of years when a parish was burned in at least 25% or 50% of its total area. PC2 (hereafter designated “burn concentration”) expressed mainly the Gini index, i.e., represented a gradient of increased temporal concentration of the total burned area (regardless of its extent), as opposed to a more homogeneous distribution of burned area throughout the period. It therefore expresses the likelihood that a parish will be significantly affected in specific years, rather than being progressively burned at regular intervals. Associated with this variable, parishes with higher PC2 coordinates also had a higher proportion of years where they burned at least 75% of the total area.
A graphical representation of the values of both PCs for both periods (Figure 3) shows that the range of PC2 values was wider for parishes with higher PC1 coordinates, regardless of the period. This demonstrates that a higher variability in the temporal concentration of total burned area is associated with parishes with higher fire incidence. At the same time, the figure shows that the second period is marked by increased variability in the values of both PCs. It is worth highlighting the upper right quadrant of the Figure 3, which reveals the most severe situation, where the parishes with higher scores of fire incidence and burn concentration are represented. It should also be noted that the second period dominates this quadrant.
This contrast between both periods can also be seen when considering the distributions of the PC score values in each period (Figure 4) with coordinates in both PC1 and PC2 less concentrated and shifted towards the right in the second period, showing an increasing trend in both fire incidence and burn concentration. In the case of PC2, the second period is also marked by the increased number of parishes with relatively high burn concentration (peak at PC2 scores higher than 2).
The spatial distribution of fire regime indicators is shown in Figure 5. The increase in fire incidence across periods was widespread in the whole area, although particularly concentrated in the center and southern sectors. The increase in burn concentration was mostly found in the western and southern sectors of the study area, with a reduced expression in its eastern limit.

3.2. Territorial Features Associated with Study Periods

Basic statistics on landscape features across periods are shown in Table 2 and Figure 6. Population density was quite variable within the region, but with a significant declining trend (ca. 10 inhabitants/km2) when comparing the two periods (Figure 6a). The main land cover (on average, ca. 40% of parish area) was pine and eucalypt plantations, which suffered a small (although significant) 3% decline across time (Figure 6b). Both farmland and native vegetation (oak forests and scrublands) covered ca. 25% of the parishes, with significant but inverse trends across time: a 4% decrease in the former (Figure 6c) and a similar increase in the latter (Figure 6d).

3.3. Influence of Period and Territorial Features on Fire Regime Indicators

The intersection between the cumulative burned areas and LULC classes in both periods is shown in Table 4. Plantations and scrubland together represented ca. 80% of the total burned areas in both periods, while oak forests and farmland accounted for 11% and 8%, respectively. Burned areas increased in the second period for all LULC classes, but their proportions of total burned areas remained practically unchanged between periods. The second period also recorded higher values in the ratio between burned areas and the LULC class area, with emphasis on scrubland (157%; 37% more than in the first period) and plantations (80%; 26% more when compared to the first period). Oak forests recorded the smallest difference between periods (3% higher in the second period). These results confirm the higher propensity to burn of scrubland and plantations, and the lower propensity to burn of farmland and oak forests.
Fire incidence at parish level was related to four variables (Figure 7; Table 5): period, population density, and the proportion of farmland and natural vegetation.
Parishes with higher population density were associated with lower fire incidence, but only during the first period, as the trend was no longer visible in the second period. Except for parishes with a very low population density, for any given value of density, fire incidence was always higher in the second period. Increasing proportional cover by farmland was associated with decreased fire incidence, although the effect was stronger in the second period. Above a threshold of ca. 40% cover by farmland, there was no difference in fire incidence across periods, but under this threshold fire incidence was higher in the second period. Parishes with higher fire incidence were also associated with a higher proportional cover by natural vegetation, regardless of the period.
Burn temporal concentration at parish level was related to three variables (Figure 8; Table 6): proportional cover by farmland, proportional cover by natural vegetation, and period.
The strong interaction between period and proportional cover by farmland was reflected in the fact that the latter variable was strongly associated with a decline in burn concentration, but only in the second period, this effect being negligible in the first period (Figure 8a). Consequently, above a threshold of ca. 50–60% cover by farmland, the period variable did not have a strong effect on this indicator. At smaller values, however, burn concentration was much higher in the second period.
Higher cover by natural vegetation was associated with lower levels of burn concentration, i.e., a more equal distribution of the total burned area across years (Figure 8b). The second period usually registered higher burn concentration, and the difference between periods was higher for parishes with lower cover by natural forests and scrublands. Overall, for the same amount of cover by natural vegetation, burn concentration was higher in the second period.

4. Discussion

4.1. Changes in Fire Regime and Territorial Features

The PCA summarized the used fire regime variables for the 1975–2018 period into two major independent indicators that we designated fire incidence and burn concentration. Fire incidence indicated the overall amount of area burned over a certain period, separating parishes with high and low fire prevalence (expressed in total area burned and fire frequency). Burn concentration expressed how the total burned area accumulated over the period (regardless of being high or low) was distributed over time, i.e., whether it was concentrated in a few years versus dispersed over the whole period. This use of the Gini coefficient provided a relevant view of an important aspect of fire regime, as parishes with total burned area concentrated in fewer years tended to have a higher proportion of years in which their total surface was more than 75% affected by wildfires, which may be an indication of higher severity [55,65] and higher socioecological damage [11,66,67,68]. Parishes with higher wildfire hazard are expected to be the ones with both higher fire incidence and higher burn concentration, which in our case are more prevalent in the western and southern sectors of the region (Figure 5), regardless of the period. For the second period, this could be explained by the large wildfires of 2017, which affected mainly these sectors.
There was a trend for both fire incidence and burn frequency at parish level to increase across time. This is concordant with most studies in the Mediterranean region describing changes in fire regimes across time, and is likely linked to differences in climate and territorial features [2,6,7,12,13,14,15,46]. For climate, a series of studies report trends of rising temperature and decreasing rainfall in the Iberian Peninsula and Mediterranean basin [7,69,70]. This is also registered for central Portugal [71,72]. In the Castelo Branco rain gauge, both the average annual temperature and the average annual rainfall increased in the second period: from 15.8 °C to 16.2 °C and from 724 mm to 789 mm, respectively. However, 2017 is simultaneously the year with the highest average annual temperature (17.5 °C) and the lowest average annual rainfall (418 mm) during the 1975–2018 period in this rain gauge, and that year also showed the largest burned areas (383,000 ha) in the study area. Thus, more detailed data would be necessary to identify significant changes in climate between the two studied periods and to draw definitive conclusions about its possible effects on fire regimes.
The two periods also differed in some key territorial features, whose effects on fire regimes cannot be clearly separated from climate effects. Population density declined, a trend also occurring at the country level [36,47,73]. Farmland areas also decreased, a likely consequence of the loss of population involved in the primary sector and the low economic profitability of marginal agricultural areas, leading to land abandonment [22,30,36,37,38,39,51]. The area covered by pine and eucalypt plantations, the dominant land cover in the region, also declined in the second period, likely due to the impact of wildfires on these plantations and their replacement by secondary succession vegetation, mainly scrubland or regeneration plantations with scrubland-like physiognomy. This could explain the increase in the surface occupied by oak and scrubland found for the second period (Table 4).

4.2. Changes in Territorial Features—Fire Regime Associations across Period

Three territorial features were correlated with the two fire regime indicators: population density, proportion of parish area covered by farmland, and proportion of parish area covered by natural vegetation. However, the way these were associated with fire regimes differed between periods.
Parishes with higher population density usually had lower fire incidence, but this association was only visible in the first period. Higher population density has traditionally been considered an important indicator of the intensity of landscape management (including in existing farmland or forests) that may lead to decreased fuel accumulation and fire hazard [9,13,16,46,74,75]. Additionally, it is also an indicator of the likelihood of early detection and more effective suppression of wildfires, in spite of the higher human-caused ignition probability [10,26,36,42,48]. However, this association was no longer visible in the second period. At first sight, this could suggest that the impact of human management on fire incidence in the territory has decreased. However, it is more likely that the lack of association could be due to changes in the population itself, such as aging (impairing its ability to intervene in the territory) or a reduced level of engagement with farmland and forest management, as suggested by prior studies [47]. In fact, the decrease in importance of farmland and forest for the employment and income of the rural population between the first and second periods has been largely documented [76,77,78], impacting that engagement for the same population density. Therefore, and except for parishes with a very low population density, fire incidence was always higher in the second period. Studies [10,79] stated that human factors are losing importance in driving fire regimes when compared to climate/weather factors and fuel moisture. Adverse weather and moisture conditions also favor the occurrence of extreme fire events [9,80,81,82,83], whose size and magnitude exceed the human capacity to suppress them [84,85]. In this sense, and according to the obtained results, population density at a parish level revealed: (1) a nonsignificant relationship with burn concentration (PC2) in both periods, confirming that this variable does not have a relevant association with expected fire severity; and (2) a negligible association with fire incidence (PC1) during the second period, when the fire regime was characterized by a few major events that contributed to an increase in overall burned area. However, the year 2017, which was marked by the occurrence of two extreme fire events (in June and October), but also the events of 2003 and 2005, certainly prevented a stronger relationship between population density and fire incidence in the second period.
The increased proportion of farmland in the parish area was strongly and negatively associated with both fire incidence and burn concentration. With regard to fire incidence, this result is in accordance with the low fire selectivity for farmland already observed by several authors [24,86,87,88,89]. Regarding burn concentration, the effect of farmland can be interpreted in two complementary ways. On the one hand, farmland areas are likely to have an increased human presence in comparison to forest or scrubland, which will promote early detection and response to ignitions, promoting smaller fires over time and leading to the temporal dispersion of area burned. On the other hand, it is likely that the increased percentage of farmland has caused a greater fragmentation in the second period for non-farmland LULC classes within a parish, such as scrubland and forest, constraining fire propagation and the occurrence of large fires that would lead to higher burn concentration. These associations were stronger in the second period, mainly for burn concentration. The difference across periods was also more obvious for parishes with less farmland: the differences between periods became negligible above a threshold of ca. 40–50% of farmland. As previously, the interaction effect between period and farmland proportion may be explained by a mixture of climate differences (e.g., causing higher flammability to the same farmland land use) and management differences (e.g., increased flammability or human presence associated with transitions between crop types, or to transitions between more intensively managed annual crops and less intensively managed grasslands). Future studies should address these hypotheses, as well as a more detailed analysis of the impacts of landscape configuration (rather than simply composition) on fire regimes.
Regardless of the period, higher fire incidence was associated with parishes with higher proportions of natural vegetation (and therefore less pine and eucalypt plantations). Possible contributors to this pattern include: (a) an effect of wildfires on land cover—high fire incidence will promote fewer forests and a higher transition from forests to scrublands [22]. As the land cover information relates to a relatively late time frame in each period (the last year of the first period and 3 years before the end of the second period), it might express the land cover effects of fire regimes; (b) the elevated fire-proneness and high rate of postfire regeneration of scrubland [25,86,90,91], which will promote both extensive and frequent fires, even though this LULC class was aggregated in this study with the less fire-prone oak forests. The interpretation of the negative relationship between natural vegetation and burn concentration is not straightforward, as this LULC class includes both oak forest and scrubland. However, considering that areas with lower cover by this LULC have higher cover by tree plantations (these LULC classes are highly negatively correlated at a parish level), this pattern may show that areas with more plantations will have a higher burn concentration. Again, for the same cover by plantations, burn concentration will be higher in the second period (and the difference increases the higher the cover by plantations).

5. Conclusions

We described the temporal change in fire regime in a fire-prone region and registered a trend for increased fire incidence and burn concentration. Three main findings seem relevant for the implementation of wildfire management policies in this highly fire-prone region.
Firstly, the key role of farmland—in a context of farmland abandonment, keeping farmland in the territory is a big challenge, but it is the key tool to decrease fire incidence and burn concentration. Our results show that a minimum of 40–50% of a parish territory covered by farmland was highly effective to buffer the increased risks associated with different management and climate contexts of the considered periods.
Secondly, the challenge of keeping population active in the primary sector—the effectiveness of a higher population density on keeping fire incidence low—was much decreased in the latter period. We hypothesized that this was due to persisting population being less engaged in farmland and forest management. If this assumption is correct, policies to keep or increase population in these marginal territories may not by themselves solve the problem. Attracting people to these areas will not necessarily be reflected in a decreased fire hazard if these “new rurals” will only stay in larger urban centers with access to services (internet, schools, health care) and not devoted to activities in the primary sector. Therefore, fostering the socioeconomic viability of farmland in this region will require coordinated policies across different sectors [92].
Finally, the negative impact of large areas with tree plantations, which cause an increase in burn concentration, make a higher level of socioecological damage more likely.

Author Contributions

Conceptualization, F.M.; data curation, F.M., M.L., R.B. and P.F.R.; formal analysis, F.M., M.L., R.B. and P.F.R.; investigation, F.M. and M.L.; methodology, F.M.; project administration, F.M. and J.L.S., software, F.M.; supervision, J.L.Z. and J.L.S.; validation, R.B., M.J.C., A.N., S.O., P.F.R., J.L.Z. and J.L.S.; visualization, F.M. and M.L.; writing—original draft preparation, F.M., M.L.; writing—review and editing, F.M., M.L., R.B., M.J.C., A.N., S.O., P.F.R., J.L.Z. and J.L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by national funds through FCT—Portuguese Foundation for Science and Technology, I.P., under the framework of the project People&Fire: Reducing Risk, Living with Risk (PCIF/AGT/0136/2017). This research was also supported by the Forest Research Centre (UID/AGR/00239/2020) and Centre of Geographical Studies (UIDB/00295/2020 and UIDP/00295/2020). Francisco Moreira was supported by contract IF/01053/2015. Sandra Oliveira was funded through FCT, I.P., under the program Stimulus of Scientific Employment—Individual Support within the contract 2020.03873.CEECIND.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. These data can be found here: URL https://www.ine.pt and https://www.dgterritorio.gov.pt (accessed on 18 March 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carlucci, M.; Zambon, I.; Colantoni, A.; Salvati, L. Socioeconomic Development, Demographic Dynamics and Forest Fires in Italy, 1961–2017: A Time-Series Analysis. Sustainability 2019, 11, 1305. [Google Scholar] [CrossRef]
  2. Chuvieco, E.; Giglio, L.; Justice, C. Global characterization of fire activity: Toward defining fire regimes from Earth observation data. Glob. Change Biol. 2008, 14, 1488–1502. [Google Scholar] [CrossRef]
  3. Krebs, P.; Pezzatti, G.B.; Mazzoleni, S.; Talbot, L.M.; Conedera, M. Fire regime: History and definition of a key concept in disturbance ecology. Theory Biosci. 2010, 129, 53–69. [Google Scholar] [CrossRef] [PubMed]
  4. Morgan, P.; Hardy, C.C.; Swetnam, T.W.; Rollins, M.G.; Long, D.G. Mapping fire regimes across time and space: Understanding coarse and fine-scale fire patterns. Int. J. Wildland Fire 2001, 10, 329–342. [Google Scholar] [CrossRef]
  5. Eskandari, S.; Miesel, J.R.; Pourghasemi, H.R. The temporal and spatial relationships between climatic parameters and fire occurrence in northeastern Iran. Ecol. Indic. 2020, 118, 106720. [Google Scholar] [CrossRef]
  6. Jiménez-Ruano, A.; Mimbrero, M.R.; de la Riva Fernández, J. Understanding wildfires in mainland Spain. A comprehensive analysis of fire regime features in a climate-human context. Appl. Geogr. 2017, 89, 100–111. [Google Scholar] [CrossRef]
  7. Chergui, B.; Fahd, S.; Santos, X.; Pausas, J.G. Socioeconomic Factors Drive Fire-Regime Variability in the Mediterranean Basin. Ecosystems 2017, 21, 619–628. [Google Scholar] [CrossRef]
  8. Amraoui, M.; Pereira, M.G.; DaCamara, C.C.; Calado, T.J. Atmospheric conditions associated with extreme fire activity in the Western Mediterranean region. Sci. Total. Environ. 2015, 524–525, 32–39. [Google Scholar] [CrossRef]
  9. Pinto, G.; Rousseu, F.; Niklasson, M.; Drobyshev, I. Effects of human-related and biotic landscape features on the occurrence and size of modern forest fires in Sweden. Agric. For. Meteorol. 2020, 291, 108084. [Google Scholar] [CrossRef]
  10. Rodrigues, M.; Jiménez-Ruano, A.; Peña-Angulo, D.; de la Riva, J. A comprehensive spatial-temporal analysis of driving factors of human-caused wildfires in Spain using Geographically Weighted Logistic Regression. J. Environ. Manag. 2018, 225, 177–192. [Google Scholar] [CrossRef] [Green Version]
  11. Calheiros, T.; Nunes, J.P.; Pereira, M. Recent evolution of spatial and temporal patterns of burnt areas and fire weather risk in the Iberian Peninsula. Agric. For. Meteorol. 2020, 287, 107923. [Google Scholar] [CrossRef]
  12. Fernandes, P.M.; Loureiro, C.; Guiomar, N.; Pezzatti, G.B.; Manso, F.T.; Lopes, L. The dynamics and drivers of fuel and fire in the Portuguese public forest. J. Environ. Manag. 2014, 146, 373–382. [Google Scholar] [CrossRef] [PubMed]
  13. Mancini, L.D.; Corona, P.; Salvati, L. Ranking the importance of Wildfires’ human drivers through a multi-model regression approach. Environ. Impact Assess. Rev. 2018, 72, 177–186. [Google Scholar] [CrossRef]
  14. Molina, C.M.; Martín, O.K.; Martín, L.G. Regional fire scenarios in Spain: Linking landscape dynamics and fire regime for wildfire risk management. J. Environ. Manag. 2018, 233, 427–439. [Google Scholar] [CrossRef] [PubMed]
  15. Mueller, S.E.; Thode, A.E.; Margolis, E.Q.; Yocom, L.L.; Young, J.D.; Iniguez, J.M. Climate relationships with increasing wildfire in the southwestern US from 1984 to 2015. For. Ecol. Manag. 2020, 460, 117861. [Google Scholar] [CrossRef]
  16. Nunes, A.; Lourenço, L.; Castro-Meira, A. Exploring spatial patterns and drivers of forest fires in Portugal (1980–2014). Sci. Total Environ. 2016, 573, 1190–1202. [Google Scholar] [CrossRef] [PubMed]
  17. Oliveira, S.; Gonçalves, A.; Zêzere, J.L. Reassessing wildfire susceptibility and hazard for mainland Portugal. Sci. Total. Environ. 2020, 762, 143121. [Google Scholar] [CrossRef]
  18. Pezzatti, G.B.; Zumbrunnen, T.; Bürgi, M.; Ambrosetti, P.; Conedera, M. Fire regime shifts as a consequence of fire policy and socio-economic development: An analysis based on the change point approach. For. Policy Econ. 2013, 29, 7–18. [Google Scholar] [CrossRef]
  19. Elia, M.; D’Este, M.; Ascoli, D.; Giannico, V.; Spano, G.; Ganga, A.; Colangelo, G.; Lafortezza, R.; Sanesi, G. Estimating the probability of wildfire occurrence in Mediterranean landscapes using Artificial Neural Networks. Environ. Impact Assess. Rev. 2020, 85, 106474. [Google Scholar] [CrossRef]
  20. Ganteaume, A.; Camia, A.; Jappiot, M.; San-Miguel-Ayanz, J.; Long-Fournel, M.; Lampin, C. A Review of the Main Driving Factors of Forest Fire Ignition Over Europe. Environ. Manag. 2012, 51, 651–662. [Google Scholar] [CrossRef] [Green Version]
  21. Saura-Mas, S.; Paula, S.; Pausas, J.G.; Lloret, F. Fuel loading and flammability in the Mediterranean Basin woody species with different post-fire regenerative strategies. Int. J. Wildland Fire 2010, 19, 783–794. [Google Scholar] [CrossRef]
  22. 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] [PubMed]
  23. Bajocco, S.; Ricotta, C. Evidence of selective burning in Sardinia (Italy): Which land-cover classes do wildfires prefer? Landsc. Ecol. 2007, 23, 241–248. [Google Scholar] [CrossRef]
  24. Oliveira, S.; Moreira, F.; Boca, R.; San-Miguel-Ayanz, J.; Pereira, J.M.C. Assessment of fire selectivity in relation to land cover and topography: A comparison between Southern European countries. Int. J. Wildland Fire 2014, 23, 620–630. [Google Scholar] [CrossRef]
  25. Carmo, M.; Moreira, F.; Casimiro, P.; Vaz, P. Land use and topography influences on wildfire occurrence in northern Portugal. Landsc. Urban Plan. 2011, 100, 169–176. [Google Scholar] [CrossRef]
  26. Calviño-Cancela, M.; Chas-Amil, M.L.; García-Martínez, E.D.; Touza, J. Wildfire risk associated with different vegetation types within and outside wildland-urban interfaces. For. Ecol. Manag. 2016, 372, 1–9. [Google Scholar] [CrossRef]
  27. Curt, T.; Borgniet, L.; Bouillon, C. Wildfire frequency varies with the size and shape of fuel types in southeastern France: Implications for environmental management. J. Environ. Manag. 2013, 117, 150–161. [Google Scholar] [CrossRef]
  28. Lloret, F.; Calvo, E.; Pons, X.; Díaz-Delgado, R. Wildfires and landscape patterns in the Eastern Iberian Peninsula. Landsc. Ecol. 2002, 17, 745–759. [Google Scholar] [CrossRef]
  29. Catry, F.; Rego, F.C.; Moreira, F.; Bacao, F. Characterizing and modelling the spatial patterns of wildfire ignitions in Portugal: Fire initiation and resulting burned area. Trans. Ecol. Environ. 2008, 119, 213–221. [Google Scholar] [CrossRef]
  30. Vilar, L.; Gómez, I.; Martínez-Vega, J.; Echavarría, P.; Riaño, D.; Martín, M.P. Multitemporal Modelling of Socio-Economic Wildfire Drivers in Central Spain between the 1980s and the 2000s: Comparing Generalized Linear Models to Machine Learning Algorithms. PLoS ONE 2016, 11, e0161344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Price, O.; Bradstock, R. Countervailing effects of urbanization and vegetation extent on fire frequency on the Wildland Urban Interface: Disentangling fuel and ignition effects. Landsc. Urban Plan. 2014, 130, 81–88. [Google Scholar] [CrossRef]
  32. Syphard, A.D.; Radeloff, V.C.; Keeley, J.E.; Hawbaker, T.J.; Clayton, M.K.; Stewart, S.I.; Hammer, R.B. Human influence on California fire regimes. Ecol. Appl. 2007, 17, 1388–1402. [Google Scholar] [CrossRef]
  33. Alcasena, F.J.; Ager, A.A.; Bailey, J.D.; Pineda, N.; Vega-García, C. Towards a comprehensive wildfire management strategy for Mediterranean areas: Framework development and implementation in Catalonia, Spain. J. Environ. Manag. 2018, 231, 303–320. [Google Scholar] [CrossRef]
  34. Costafreda-Aumedes, S.; Comas, C.; Vega-Garcia, C. Human-caused fire occurrence modelling in perspective: A review. Int. J. Wildland Fire 2017, 26, 983–998. [Google Scholar] [CrossRef]
  35. Benali, A.; Sá, A.; Pinho, J.; Fernandes, P.; Pereira, J. Understanding the Impact of Different Landscape-Level Fuel Management Strategies on Wildfire Hazard in Central Portugal. Forests 2021, 12, 522. [Google Scholar] [CrossRef]
  36. Nunes, A.N. Regional variability and driving forces behind forest fires in Portugal an overview of the last three decades (1980–2009). Appl. Geogr. 2012, 34, 576–586. [Google Scholar] [CrossRef]
  37. Rodrigues, M.; Jiménez-Ruano, A.; de la Riva, J. Fire regime dynamics in mainland Spain. Part 1: Drivers of change. Sci. Total. Environ. 2019, 721, 135841. [Google Scholar] [CrossRef]
  38. Torres, F.; Fernandes, P.; Pinto, R.; Botelho, H.; Monzon, A. Regional livestock grazing, human demography and fire incidence in the Portuguese landscape. For. Syst. 2014, 23, 15. [Google Scholar] [CrossRef]
  39. Canadas, M.J.; Leal, M.; Soares, F.; Novais, A.; Ribeiro, P.F.; Schmidt, L.; Delicado, A.; Moreira, F.; Bergonse, R.; Oliveira, S.; et al. Wildfire mitigation and adaptation: Two locally independent actions supported by different policy domains. Land Use Policy 2023, 124, 106444. [Google Scholar] [CrossRef]
  40. Parente, J.; Pereira, M.; Amraoui, M.; Fischer, E. Heat waves in Portugal: Current regime, changes in future climate and impacts on extreme wildfires. Sci. Total. Environ. 2018, 631–632, 534–549. [Google Scholar] [CrossRef]
  41. Sarricolea, P.; Serrano-Notivoli, R.; Fuentealba, M.; Hernández-Mora, M.; de la Barrera, F.; Smith, P.; Meseguer-Ruiz, Ó. Recent wildfires in Central Chile: Detecting links between burned areas and population exposure in the wildland urban interface. Sci. Total. Environ. 2019, 706, 135894. [Google Scholar] [CrossRef] [PubMed]
  42. Sousa, P.M.; Trigo, R.M.; Pereira, M.G.; Bedia, J.; Gutiérrez, J.M. Different approaches to model future burnt area in the Iberian Peninsula. Agric. For. Meteorol. 2015, 202, 11–25. [Google Scholar] [CrossRef]
  43. Da Camara, C.C.; Pereira, M.G.; Calado, T.J.; Calheiros, T. Impacts of Climate Change on the Fire Regime in Portugal. In Advances in Forest Fire Research; Viegas, D.X., Ed.; Imprensa da Universidade de Coimbra: Coimbra, Portugal, 2014; pp. 1743–1748. [Google Scholar] [CrossRef]
  44. Bergonse, R.; Oliveira, S.; Gonçalves, A.; Nunes, S.; DaCamara, C.; Zêzere, J.L. Predicting burnt areas during the summer season in Portugal by combining wildfire susceptibility and spring meteorological conditions. Geomat. Nat. Hazards Risk 2021, 12, 1039–1057. [Google Scholar] [CrossRef]
  45. Silva, J.M.N.; Moreno, M.V.; Le Page, Y.; Oom, D.; Bistinas, I.; Pereira, J.M.C. Spatiotemporal trends of area burnt in the Iberian Peninsula, 1975–2013. Reg. Environ. Change 2018, 19, 515–527. [Google Scholar] [CrossRef]
  46. Pausas, J.G.; Fernández-Muñoz, S. Fire regime changes in the Western Mediterranean Basin: From fuel-limited to drought-driven fire regime. Clim. Change 2011, 110, 215–226. [Google Scholar] [CrossRef]
  47. Oliveira, S.; Zêzere, J.L.; Queirós, M.; Pereira, J.M. Assessing the social context of wildfire-affected areas. The case of mainland Portugal. Appl. Geogr. 2017, 88, 104–117. [Google Scholar] [CrossRef]
  48. Jiménez-Ruano, A.; Fernández, J.D.L.R.; Rodrigues, M. Fire regime dynamics in mainland Spain. Part 2: A near-future prospective of fire activity. Sci. Total. Environ. 2019, 705, 135842. [Google Scholar] [CrossRef]
  49. San-Miguel-Ayanz, J.; Durrant, T.; Boca, R.; Maianti, P.; Libertá, G.; Artés-Vivancos, T.; Oom, D.; Branco, A.; de Rigo, D.; Ferrari, D.; et al. Forest Fires in Europe, Middle East and North Africa 2019; Joint Research Centre, European Commission: Luxembourg, 2020. [Google Scholar] [CrossRef]
  50. Rodrigues, M.; Costafreda-Aumedes, S.; Comas, C.; Vega-García, C. Spatial stratification of wildfire drivers towards enhanced definition of large-fire regime zoning and fire seasons. Sci. Total. Environ. 2019, 689, 634–644. [Google Scholar] [CrossRef]
  51. Tonini, M.; Parente, J.; Pereira, M.G. Global assessment of rural–urban interface in Portugal related to land cover changes. Nat. Hazards Earth Syst. Sci. 2018, 18, 1647–1664. [Google Scholar] [CrossRef]
  52. Parente, J.; Pereira, M.G.; Tonini, M. Space-time clustering analysis of wildfires: The influence of dataset characteristics, fire prevention policy decisions, weather and climate. Sci. Total. Environ. 2016, 559, 151–165. [Google Scholar] [CrossRef]
  53. Gini, C. Measurement of Inequality of Incomes. Econ. J. 1921, 31, 124. [Google Scholar] [CrossRef]
  54. Miralles Ortega, R. Anàlisi Climàtica i Ambiental dels Incendis Forestals de Catalunya (1968–2008); Universitat de Barcelona: Barcelona, Spain, 2016. [Google Scholar]
  55. Royé, D.; Tedim, F.; Martin-Vide, J.; Salis, M.; Vendrell, J.; Lovreglio, R.; Bouillon, C.; Leone, V. Wildfire burnt area patterns and trends in Western Mediterranean Europe via the application of a concentration index. Land Degrad. Dev. 2019, 31, 311–324. [Google Scholar] [CrossRef]
  56. Loepfe, L.; Martinez-Vilalta, J.; Oliveres, J.; Piñol, J.; Lloret, F. Feedbacks between fuel reduction and landscape homogenisation determine fire regimes in three Mediterranean areas. For. Ecol. Manag. 2010, 259, 2366–2374. [Google Scholar] [CrossRef]
  57. Barreal, J.; Jannes, G. Spatial and Temporal Wildfire Decomposition as a Tool for Assessment and Planning of an Efficient Forest Policy in Galicia (Spain). Forests 2020, 11, 811. [Google Scholar] [CrossRef]
  58. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 21 April 2022).
  59. Wigtil, G.; Hammer, R.B.; Kline, J.D.; Mockrin, M.H.; Stewart, S.I.; Roper, D.; Radeloff, V.C. Places where wildfire potential and social vulnerability coincide in the coterminous United States. Int. J. Wildland Fire 2016, 25, 896–908. [Google Scholar] [CrossRef]
  60. Wickham, H. ggplot2: Elegant Graphics for Data Analysis, 2nd ed.; Springer: New York, NY, USA, 2016. [Google Scholar] [CrossRef]
  61. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 2015, 67, 48. [Google Scholar] [CrossRef]
  62. Naimi, B.; Hamm, N.A.S.; Groen, T.A.; Skidmore, A.K.; Toxopeus, A.G. Where is positional uncertainty a problem for species distribution modelling? Ecography 2013, 37, 191–203. [Google Scholar] [CrossRef]
  63. Zuur, A.F.; Ieno, E.N.; Elphick, C.S. A protocol for data exploration to avoid common statistical problems. Methods Ecol. Evol. 2010, 1, 3–14. [Google Scholar] [CrossRef]
  64. Long, J.A. jtools: Analysis and Presentation of Social Scientific Data. 2020. Available online: https://cran.r-project.org/package=jtools (accessed on 2 May 2022).
  65. Tedim, F.; Leone, V.; Amraoui, M.; Bouillon, C.; Coughlan, M.R.; Delogu, G.M.; Fernandes, P.M.; Ferreira, C.; McCaffrey, S.; McGee, T.K.; et al. Defining Extreme Wildfire Events: Difficulties, Challenges, and Impacts. Fire 2018, 1, 9. [Google Scholar] [CrossRef]
  66. Liu, N.; Lei, J.; Gao, W.; Chen, H.; Xie, X. Combustion dynamics of large-scale wildfires. Proc. Combust. Inst. 2021, 38, 157–198. [Google Scholar] [CrossRef]
  67. García-Llamas, P.; Suárez-Seoane, S.; Taboada, A.; Fernández-Manso, A.; Quintano, C.; Fernández-García, V.; Fernández-Guisuraga, J.M.; Marcos, E.; Calvo, L. Environmental drivers of fire severity in extreme fire events that affect Mediterranean pine forest ecosystems. For. Ecol. Manag. 2018, 433, 24–32. [Google Scholar] [CrossRef]
  68. Mourão, P.R.; Martinho, V.D. The choices of the fire—Debating socioeconomic determinants of the fires observed at Portuguese municipalities. For. Policy Econ. 2014, 43, 29–40. [Google Scholar] [CrossRef]
  69. Vicente-Serrano, S.M.; Lopez-Moreno, I.; Beguería, S.; Lorenzo-Lacruz, J.; Sanchez-Lorenzo, A.; García-Ruiz, J.M.; Azorin-Molina, C.; Morán-Tejeda, E.; Revuelto, J.; Trigo, R.; et al. Evidence of increasing drought severity caused by temperature rise in southern Europe. Environ. Res. Lett. 2014, 9, 044001. [Google Scholar] [CrossRef]
  70. Pausas, J.G. Changes in Fire and Climate in the Eastern Iberian Peninsula (Mediterranean Basin). Clim. Change 2004, 63, 337–350. [Google Scholar] [CrossRef]
  71. Portela, M.; Espinosa, L.; Zelenakova, M. Long-Term Rainfall Trends and Their Variability in Mainland Portugal in the Last 106 Years. Climate 2020, 8, 146. [Google Scholar] [CrossRef]
  72. Santo, F.E.; de Lima, M.I.P.; Ramos, A.M.; Trigo, R.M. Trends in seasonal surface air temperature in mainland Portugal, since 1941. Int. J. Clim. 2013, 34, 1814–1837. [Google Scholar] [CrossRef]
  73. Meira Castro, A.C.; Nunes, A.; Sousa, A.; Lourenço, L. Mapping the Causes of Forest Fires in Portugal by Clustering Analysis. Geosciences 2020, 10, 53. [Google Scholar] [CrossRef] [Green Version]
  74. Curth, M.D.T.; Biscayart, C.; Ghermandi, L.; Pfister, G. Wildland–Urban Interface Fires and Socioeconomic Conditions: A Case Study of a Northwestern Patagonia City. Environ. Manag. 2012, 49, 876–891. [Google Scholar] [CrossRef]
  75. Canadas, M.J.; Novais, A. Bringing local socioeconomic context to the analysis of forest owners’ management. Land Use Policy 2014, 41, 397–407. [Google Scholar] [CrossRef]
  76. Baptista, F.O. O Espaço Rural: Declínio da Agricultura; Celta: Lisbon, Portugal, 2010. [Google Scholar]
  77. Rolo, J.C.; Cordovil, F. Rural, Agriculturas e Políticas; ANIMAR: Lisboa, Portugal, 2014. [Google Scholar]
  78. Radich, M.C.; Baptista, F.O. Terra e Agricultura: Século XX; 100LUZ: Lisboa, Portugal, 2021. [Google Scholar]
  79. Rodrigues, M.; Jiménez, A.; de la Riva, J. Analysis of recent spatial–temporal evolution of human driving factors of wildfires in Spain. Nat. Hazards 2016, 84, 2049–2070. [Google Scholar] [CrossRef]
  80. Barbero, R.; Curt, T.; Ganteaume, A.; Maillé, E.; Jappiot, M.; Bellet, A. Simulating the effects of weather and climate on large wildfires in France. Nat. Hazards Earth Syst. Sci. 2019, 19, 441–454. [Google Scholar] [CrossRef]
  81. Turco, M.; von Hardenberg, J.; AghaKouchak, A.; Llasat, M.C.; Provenzale, A.; Trigo, R.M. On the key role of droughts in the dynamics of summer fires in Mediterranean Europe. Sci. Rep. 2017, 7, 81. [Google Scholar] [CrossRef] [PubMed]
  82. Parente, J.; Amraoui, M.; Menezes, I.; Pereira, M. Drought in Portugal: Current regime, comparison of indices and impacts on extreme wildfires. Sci. Total. Environ. 2019, 685, 150–173. [Google Scholar] [CrossRef] [PubMed]
  83. Ruffault, J.; Curt, T.; Martin-StPaul, N.K.; Moron, V.; Trigo, R.M. Extreme wildfire events are linked to global-change-type droughts in the northern Mediterranean. Nat. Hazards Earth Syst. Sci. 2018, 18, 847–856. [Google Scholar] [CrossRef]
  84. Bergons, R.; Oliveira, S.; Gonçalves, A.; Nunes, S.; da Câmara, C.; Zêzere, J.L. A combined structural and seasonal approach to assess wildfire susceptibility and hazard in summertime. Nat. Hazards 2021, 106, 2545–2547. [Google Scholar] [CrossRef]
  85. Duane, A.; Brotons, L. Synoptic weather conditions and changing fire regimes in a Mediterranean environment. Agric. For. Meteorol. 2018, 253–254, 190–202. [Google Scholar] [CrossRef]
  86. Moreira, F.; Vaz, P.; Catry, F.; Silva, J. Regional variations in wildfire susceptibility of land-cover types in Portugal: Implications for landscape management to minimize fire hazard. Int. J. Wildland Fire 2009, 18, 563–574. [Google Scholar] [CrossRef]
  87. Bistinas, I.; Harrison, S.P.; Prentice, I.C.; Pereira, J.M.C. Causal relationships versus emergent patterns in the global controls of fire frequency. Biogeosciences 2014, 11, 5087–5101. [Google Scholar] [CrossRef]
  88. Andela, N.; Morton, D.C.; Giglio, L.; Chen, Y.; van der Werf, G.R.; Kasibhatla, P.S.; DeFries, R.S.; Collatz, G.J.; Hantson, S.; Kloster, S.; et al. A human-driven decline in global burned area. Science 2017, 356, 1356–1362. [Google Scholar] [CrossRef]
  89. Calviño-Cancela, M.; Chas-Amil, M.L.; García-Martínez, E.D.; Touza, J. Interacting effects of topography, vegetation, human activities and wildland-urban interfaces on wildfire ignition risk. For. Ecol. Manag. 2017, 397, 10–17. [Google Scholar] [CrossRef]
  90. Nunes, M.C.; Vasconcelos, M.J.; Pereira, J.M.; Dasgupta, N.; Alldredge, R.J.; Rego, F.C. Land Cover Type and Fire in Portugal: Do Fires Burn Land Cover Selectively? Landsc. Ecol. 2005, 20, 661–673. [Google Scholar] [CrossRef]
  91. Marques, S.; Borges, J.G.; Garcia-Gonzalo, J.; Moreira, F.; Carreiras, J.; Oliveira, M.; Cantarinha, A.; Botequim, B.; Pereira, J.M.C. Characterization of wildfires in Portugal. Eur. J. For. Res. 2011, 130, 775–784. [Google Scholar] [CrossRef]
  92. Lomba, A.; Moreira, F.; Klimek, S.; Jongman, R.H.; Sullivan, C.; Moran, J.; Poux, X.; Honrado, J.P.; Pinto-Correia, T.; Plieninger, T.; et al. Back to the future: Rethinking socioecological systems underlying high nature value farmlands. Front. Ecol. Environ. 2019, 18, 36–42. [Google Scholar] [CrossRef] [Green Version]
Fire 06 00060 g001 550
Figure 1. Study area in central Portugal, civil parish boundaries and land use/land cover (LULC) of 2015 in the region. LULC source: Portuguese Directorate-General for Territory.
Figure 1. Study area in central Portugal, civil parish boundaries and land use/land cover (LULC) of 2015 in the region. LULC source: Portuguese Directorate-General for Territory.
Fire 06 00060 g001
Fire 06 00060 g002 550
Figure 2. Variation in the total annual burned area (in thousand hectares) across the two studied periods (first = 1975–1995; second = 1996–2018).
Figure 2. Variation in the total annual burned area (in thousand hectares) across the two studied periods (first = 1975–1995; second = 1996–2018).
Fire 06 00060 g002
Fire 06 00060 g003 550
Figure 3. Plot of the values of the two extracted principal components (PC) by parish for the two studied periods (first = 1975–1995; second = 1996–2018).
Figure 3. Plot of the values of the two extracted principal components (PC) by parish for the two studied periods (first = 1975–1995; second = 1996–2018).
Fire 06 00060 g003
Fire 06 00060 g004 550
Figure 4. Kernel density estimates of PC1 (a) and PC2 (b) scores in the two periods (first = 1975–1995; second = 1996–2018). Vertical lines represent median values for each period.
Figure 4. Kernel density estimates of PC1 (a) and PC2 (b) scores in the two periods (first = 1975–1995; second = 1996–2018). Vertical lines represent median values for each period.
Fire 06 00060 g004
Fire 06 00060 g005 550
Figure 5. Spatial patterns of parish coordinates in the PCA axes, across periods. Above, PC1 (fire incidence), below, PC2 (burn concentration). Color codes represent the quartiles of the PC scores values.
Figure 5. Spatial patterns of parish coordinates in the PCA axes, across periods. Above, PC1 (fire incidence), below, PC2 (burn concentration). Color codes represent the quartiles of the PC scores values.
Fire 06 00060 g005
Fire 06 00060 g006 550
Figure 6. Kernel density estimates of territorial features in the two periods (first = 1975–1995; second = 1996–2018). (a) Population density (log-transformed); (b) proportion of farmland (%); (c) proportion of plantations (%); and (d) proportion of natural vegetation (%). Vertical lines represent median values for each period.
Figure 6. Kernel density estimates of territorial features in the two periods (first = 1975–1995; second = 1996–2018). (a) Population density (log-transformed); (b) proportion of farmland (%); (c) proportion of plantations (%); and (d) proportion of natural vegetation (%). Vertical lines represent median values for each period.
Fire 06 00060 g006
Fire 06 00060 g007 550
Figure 7. Predicted value (and 95% confidence intervals) of significant variables associated with PC1 parish coordinates, expressing fire incidence. For each variable, all the others are set to their mean (or base level) values. (a) Population density (log-transformed); (b) proportion of farmland (%); and (c) proportion of natural vegetation (%). See model detail in Table 5.
Figure 7. Predicted value (and 95% confidence intervals) of significant variables associated with PC1 parish coordinates, expressing fire incidence. For each variable, all the others are set to their mean (or base level) values. (a) Population density (log-transformed); (b) proportion of farmland (%); and (c) proportion of natural vegetation (%). See model detail in Table 5.
Fire 06 00060 g007
Fire 06 00060 g008 550
Figure 8. Predicted value (and 95% confidence intervals) of significant variables associated with PC2 parish coordinates, expressing burn concentration over time. For each variable, all the others are set to their mean (or base level) values. (a) Proportion of farmland (%) and (b) proportion of natural vegetation (%). See model detail in Table 6.
Figure 8. Predicted value (and 95% confidence intervals) of significant variables associated with PC2 parish coordinates, expressing burn concentration over time. For each variable, all the others are set to their mean (or base level) values. (a) Proportion of farmland (%) and (b) proportion of natural vegetation (%). See model detail in Table 6.
Fire 06 00060 g008
Table
Table 1. Median, mean, and range (min–max) for the fire regime variables estimated for each parish, in each of the two periods.
Table 1. Median, mean, and range (min–max) for the fire regime variables estimated for each parish, in each of the two periods.
VariableDescriptionFirst Period
(1975–1995)		Second Period
(1996–2018)
Total burned percentageTotal accumulated burned area expressed as % of the parish area.49.3/58.1
(0.1–259.3)		73.8/80.1
(0.3–329.5)
GiniGini index of the yearly distribution of the total burned area81.5/81.1
(58.8–95.2)		87.7/85.4
(53.9–95.7)
Mean times burnedAverage number of years a pixel inside the parish was burned0.49/0.58
(0.0–2.6)		0.74/0.80
(0.0–3.3)
Max times burnedMaximum number of years a pixel inside the parish was burned3/3.5
(1.0–8.0)		3/3.6
(1.0–10.0)
Prop years 25%Proportion of years that a parish burned over 25% of its area0/0.02
(0.0–0.24)		0.04/0.04
(0.0–0.22)
Prop years 50%Proportion of years that a parish burned over 50% of its area0/0.0
(0.0–0.14)		0/0.02
(0.0–0.13)
Prop years 75%Proportion of years that a parish burned over 75% of its area0/0.0
(0.0–0.05)		0/0.01
(0.0–0.09)
Table
Table 2. Median and range for territorial variables estimated for each parish in each of the two periods.
Table 2. Median and range for territorial variables estimated for each parish in each of the two periods.
Variable NameFirst Period
(1975–1995)		Second Period
(1996–2018)
Population density (inhabitants/km2)45.1 (7.4–810.6)35.3 (4.2–741.1)
Proportion of area covered by farmland (%)26.4 (1.8–79.3)22.4 (0.7–75.3)
Proportion of area covered by plantations (%)41.5 (0–92.2)39.1 (0–88.8)
Proportion of area covered by natural vegetation (%)24.8 (0.6–79.8)28.9 (0.9–91.3)
Table
Table 3. Variable loadings for the two extracted principal components (PC) and explained variance.
Table 3. Variable loadings for the two extracted principal components (PC) and explained variance.
VariablePC1PC2
Total burned percentage0.98−0.05
Gini−0.290.88
Mean times burned0.98−0.06
Max times burned0.64−0.67
Prop years 25%0.830.17
Prop years 50%0.700.56
Prop years 75%0.460.59
Explained variance (%)53.927.4
Table
Table 4. Burned areas per LULC class in each of the two periods.
Table 4. Burned areas per LULC class in each of the two periods.
LULC ClassesFirst Period (1975–1995)Second Period (1996–2018)
LULC AreaBurned AreasLULC AreaBurned Areas
ha%ha% of
Total		% of
LULC Area		ha%ha% of
Total		% of
LULC Area
Farmland368,57024.270,9408.119.2308,42020.298,3708.031.9
Plantations675,59044.3363,19041.453.8650,84042.7520,48042.180.0
Natural vegetationOak forests148,6509.892,41010.562.2202,06013.3130,78010.664.7
Scrubland281,86018.5338,95038.7120.3298,25019.6468,68037.9157.1
Total (study area)1,524,190100876,51010057.51,524,1901001,237,02010081.2
Table
Table 5. Significance of variables and interactions included in final mixed model to predict PC1. Model fit information: AIC = 4306.89, BIC = 4352.06. Pseudo-R2 (fixed effects) = 0.27. Pseudo-R2 (total) = 0.49. See model predicted values in Figure 7.
Table 5. Significance of variables and interactions included in final mixed model to predict PC1. Model fit information: AIC = 4306.89, BIC = 4352.06. Pseudo-R2 (fixed effects) = 0.27. Pseudo-R2 (total) = 0.49. See model predicted values in Figure 7.
VariableF-ValueSignificance
prop_naturalveg174.2<0.0001
log_pop_dens:period7.000.0083
period:prop_farmland7.480.0064
Table
Table 6. Significance of variables and interactions included in final mixed model to predict PC2. Model fit information: AIC = 3528.18, BIC = 3568.33. Pseudo-R2 (fixed effects) = 0.28. Pseudo-R2 (total) = 0.52. See model predicted values in Figure 8.
Table 6. Significance of variables and interactions included in final mixed model to predict PC2. Model fit information: AIC = 3528.18, BIC = 3568.33. Pseudo-R2 (fixed effects) = 0.28. Pseudo-R2 (total) = 0.52. See model predicted values in Figure 8.
VariableF-ValueSignificance
period:prop_naturalveg13.10.0003
period:prop_farmland21.6<0.0001
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

Moreira, F.; Leal, M.; Bergonse, R.; Canadas, M.J.; Novais, A.; Oliveira, S.; Ribeiro, P.F.; Zêzere, J.L.; Santos, J.L. Recent Trends in Fire Regimes and Associated Territorial Features in a Fire-Prone Mediterranean Region. Fire 2023, 6, 60. https://doi.org/10.3390/fire6020060

AMA Style

Moreira F, Leal M, Bergonse R, Canadas MJ, Novais A, Oliveira S, Ribeiro PF, Zêzere JL, Santos JL. Recent Trends in Fire Regimes and Associated Territorial Features in a Fire-Prone Mediterranean Region. Fire. 2023; 6(2):60. https://doi.org/10.3390/fire6020060

Chicago/Turabian Style

Moreira, Francisco, Miguel Leal, Rafaello Bergonse, Maria João Canadas, Ana Novais, Sandra Oliveira, Paulo Flores Ribeiro, José Luís Zêzere, and José Lima Santos. 2023. "Recent Trends in Fire Regimes and Associated Territorial Features in a Fire-Prone Mediterranean Region" Fire 6, no. 2: 60. https://doi.org/10.3390/fire6020060

APA Style

Moreira, F., Leal, M., Bergonse, R., Canadas, M. J., Novais, A., Oliveira, S., Ribeiro, P. F., Zêzere, J. L., & Santos, J. L. (2023). Recent Trends in Fire Regimes and Associated Territorial Features in a Fire-Prone Mediterranean Region. Fire, 6(2), 60. https://doi.org/10.3390/fire6020060

Article Metrics

Back to TopTop