require(pacman)
p_load(sf, here, terra, geodata, tidyverse, ggspatial, mapview, elevatr,
rasterVis, tidyterra, cptcity, RColorBrewer, spatialEco, report, lubridate, RStoolbox, flextable, leafpop)
Analysis of summer 2022 wildfires blast in Folgoso do Courel, Spain
1 Introduction
During the summer of 2022, the region of Galicia was devastated by wildfires, with over 40,000 hectares burned in just July and August. The fire outbreak began with an intense dry storm on the evening of July \(15^{th}\), which brought with it over 6,000 lightning strikes in less than 4 hours. The same week, temperatures in Ourense reached a record high of 44ºC, with similar heat waves reported throughout the province. These extreme weather conditions created the perfect storm for the spread of devastating wildfires.
The year 2022 was marked by a severe drought, which affected not only the summer months but also late winter and spring, a crucial time for replenishing the trees’ water supply. This shortage of water caused significant stress to the vegetation, which was further exacerbated by the ongoing rural abandonment. With no human intervention, dead biomass accumulated in large quantities, creating dangerous conditions for wildfires. The unusual weather conditions combined with the high amount of fuel available led to an unprecedented number of wildfires, which overwhelmed the firefighting resources and were difficult to contain. The existing firefighting infrastructure was not equipped to handle such a catastrophic situation.
In this study, I studied how the vegetation stress evolved over time. To relate it, a comparison was made of a relatively normal (average precipitation) year as it was 2021 with 954 registered wildfires and only 4,400 ha burned, against a dry year as 2022 with more than 40,000 ha burned in only 2 months (July and August).
The aim of this work was to study the spatio-temporal dynamics using vegetation indices derived from remote sensing images. The specific objectives were: (1) to study spatio-temporal dynamics of vegetation statues using the well-known Enhanced Vegetation Index (EVI); (2) to study the spatio-temporal dynamics of vegetation moisture through the Normalized Difference Water Index (NDWI); (3) to define the burned area and its severity.
2 Materials and Methods
2.1 Packages
2.2 Study area
The selected study area is the subregion of Folgoso de Courel (Galicia) which was severely affected by wildfires.
# Download spanish level 4 boundaries
<- gadm(country = 'ESP',
spain level = 4,
path = here('inputs/'),
resolution = 2)
# Filter study area
<- spain |>
aoi st_as_sf() |>
filter(NAME_4 == 'Folgoso do Courel')
In Fig. 1 we can see the location of the study area. Setting on the Esri.WorldImagery map, we can see that the region is very montanious, and mostly covered by forests.
Code
mapview(aoi,
viewer.suppress = mapviewGetOption("viewer.suppress"),
color = 'red',
lwd = 3,
alpha.regions = 0,
label = 'Folgoso do Courel',
legend = FALSE
)
We can illustrate this by plotting the slope of the study area. To do so, we can download the DEM with the next chnuk of code, and compute the slope with the terra::terrain()
function.
# Download the DEM
<- get_elev_raster(
dem locations = aoi,
z = 14,
clip = "bbox"
|>
) rast()
# Compute the slope
<- terrain(dem, v = 'slope') slope
Next, a graphical representation of the DEM and the slope is presented in Fig. 2 and Fig. 3, respectively.
Code
<- brewer.pal(n = 10, name = 'RdYlBu')
pal plot(dem, col = pal, main = 'DEM')
plot(slope, col = pal, main = 'Slope')


3 Resultados
Ardieron 40,000 hectáreas