Wildfires cause substantial losses to socio-economic and natural assets, especially in Mediterranean climate regions. Despite human activity being the main cause of wildfires in Mediterranean European countries, lightning-ignited wildfires should also be considered a major disruptive agent as they can trigger large fires. In addition, recent studies on the potential climate change effects on wildfires pointed out that lightning-ignited wildfires may gain relevance in Mediterranean areas in the years to come. The present study analyses the synoptical weather patterns favouring lightning-ignited wildfires in Catalonia (NE Iberian Peninsula). Being able to identify areas with an elevated lightning-ignition survival at daily timescales would be of great assistance to wildfire management agencies, i.e. locating ignitions and potential holdover fires, preparing for days with multiple ignitions or routing detection flight paths. It is worth noticing that one of the reasons that lightning-caused wildfires are difficult to manage is that they can survive for several days after the ignition, emerging days later once surface vegetation becomes dry enough to support sustained combustion. For this reason, in a first step, a reliable lightning–wildfire association is needed to properly identify the date and time of the fire starter for each wildfire. Afterwards, the circulation types on the days of ignition are analysed. The study relies on a dataset of 870 lightning-ignited wildfires, gathered by the Forest Protection Agency of the autonomous government of Catalonia between 2005 and 2020. Lightning data were provided by the Lightning Location System operated by the Meteorological Service of Catalonia. Results show that lightning-ignited wildfires in Catalonia are related to a typical synoptic weather pattern dominated by a short-wave trough at 500 hPa, with three distinct associations: an Iberian thermal low (51 % of the fires), a northern flow (24 %) and prefrontal convection (13 %).
Wildfires cause substantial losses to socio-economic and natural assets, especially in Mediterranean climate regions. Causality records show that the cause of wildfires in Mediterranean countries is mostly related to human activity (Ganteaume et al., 2013). Even if lightning is a minor cause of forest fires, lightning-ignited wildfires (hereafter LIWs) must be considered a major disruptive agent in Mediterranean climate regions as they can trigger large fires (San-Miguel-Ayanz et al., 2013). Wildfires generally occur under extreme meteorological conditions (Pereira et al., 2005; Oliveras et al., 2009). Recurrent dry spells lead to extended periods of low fuel moisture, allowing for more opportunities for fire ignitions to end in severe fires (Jain et al., 2018; Westerling et al., 2006). Under these conditions, thunderstorms can cause multiple, relatively simultaneous LIWs (Keeley and Syphard, 2021). In addition, since lightning ignitions can occur anywhere, ignitions in complex topography can turn into complicated LIW episodes that may overwhelm fire brigades (Costafreda-Aumedes et al., 2016).
To eventually become a LIW, ignitions must overcome a complex process, involving vegetation type, fuel moisture and weather conditions conducive to both fire survival and growth (Latham and Williams, 2001; Anderson, 2002). Forest fuel characteristics, especially moisture content of the upper organic layers, play a key role in the survival of the lightning-caused ignition (Flannigan and Wotton, 1991; Duff et al., 2017; Pineda et al., 2022). Fuel moisture content is increased by groundwater availability, atmospheric humidity and especially precipitation (Flannigan and Wotton, 1991; Evett et al., 2008; Morin et al., 2015). Therefore, weather characteristics such as rainfall, relative humidity and temperature on the ignition spot will have a direct impact on the probability of survival (Hall, 2007; Benson et al., 2009; Amraoui et al., 2013).
Although thunderstorms produce thousands of potential fire starters every year (lightning with long continuous currents), only few lightning-caused ignitions survive to eventually reach a flaming stage (Latham and Williams, 2001; Podur et al., 2003). The relative timing between the lightning strike and precipitation is critical to whether lightning ignitions will survive and turn into flaming combustion. Past research has dealt with the relationship between lightning-ignited wildfires and precipitation (Hall, 2007, 2008; Dowdy and Mills, 2012; Pineda and Rigo, 2017; Soler et al., 2021). Most lightning-caused ignitions are rapidly extinguished by thunderstorms' concurrent heavy rainfall; only a small fraction of lightning strikes outside the main rain shaft survive. Rorig and Ferguson (1999) defined these strikes as “dry lightning”, as they occur with less than 2 mm of precipitation. Dry lightning also tends to occur under high cloud base thunderstorms (Nauslar, 2013). Still, LIWs can occur with precipitation higher than 2 mm. Pineda and Rigo (2017) revealed that, in Catalonia, 25 % of the lightning strokes related to wildfire ignitions had no associated precipitation at all, 40 % had less than 2 mm of precipitation and 90 % had less than 10 mm. When the fuel moisture content is high but less than the moisture content of extinction, the lightning-caused ignition can smoulder as a “holdover fire”, surviving for hours to several days, until the surface vegetation becomes dry enough to support sustained combustion (Flannigan and Wotton, 1991). In the region of study, Pineda and Rigo (2017) found that latent fires above 24 h are not common (12 %), and only 5 % of the LIWs have a holdover period above 3 d.
In addition to the significance of dry lightning in starting fires, classifications of the synoptic weather patterns (hereafter SWPs) prevailing at large scales offer the boundary conditions for integrating weather-related factors into our understanding of fire regime attributes at regional scales (Duane and Brotons, 2018). Indeed, atmospheric factors are key issues, both in the processes of ignition and wildfire propagation (Pyne et al., 1996). A detailed classification of the synoptic situations on the days with wildfires due to natural causes will be a powerful tool to improve fire risk forecast and to establish more effective alarm systems (Rorig et al., 2007; Garcia-Ortega et al., 2011; Resco de Dios et al., 2022). SWP classifications based on the use of multivariate statistics make it possible to obtain objective and reliable results. Relations between large-scale circulation patterns and wildland fire severity have been studied in Canada (Skinner et al., 2002), California (Rorig and Ferguson, 1999; van Wagtendonk and Cayan, 2008), the Alpine region of central Europe (Wastl et al., 2013), the Iberian Peninsula (Millán et al., 1998; Rasilla et al., 2010; Garcia-Ortega et al., 2011), Portugal (Pereira et al., 2005) and Greece (Kassomenos, 2010), among others.
This study was undertaken to further explore the relationship between SWPs and LIWs. Although the datasets and methods of this study partially overlap with those in prior studies, their scientific objectives are different. The same methodology was previously applied to identify SWPs favouring hail episodes (Aran et al., 2011) and strong wind episodes in Catalonia (Peña et al., 2011). Similarly, it was used to analyse SWPs leading to heatwaves in the city of Barcelona (Peña et al., 2015). On the other hand, LIWs in the region have been already studied in Pineda et al. (2014), Pineda and Rigo (2017) and Soler et al. (2021). The present study extends the science of these earlier works, in this case focusing on the SWPs conducive to LIW episodes.
The study area was Catalonia, a region of
The spatial distribution of lightning in Catalonia is governed by the
distance to the Mediterranean and local orographic features like the
Pyrenees, leading to flow deviations and, consequently, convergence zones
(Campins et al., 1995; Pascual and Callado, 2002; Pineda et al., 2011).
Lightning activity in Catalonia peaks during the summer months (JJA). The
thunderstorm season gently starts around late April as surface heating
becomes the main source of instability and convection. Lightning activity
starts moving to the seashore from mid-September, where it becomes dominant
in autumn. This change is related to the evolution of the land and sea average
temperature, as the sea surface is warmer compared with land from
mid-September to April. Regarding the daily cycle of lightning activity in
Catalonia, the displayed pattern is clearly related to the solar heating
cycle, with an increase
This study relies on 16 years (2005–2020) of lightning data gathered by the Lightning Location System (LLS) of the Servei Meteorològic de Catalunya (SMC) (for details, see Pineda and Montanyà, 2009). Cloud-to-ground (CG) lightning data were used to define a thunderstorm event as a period of 6 h with more than 100 CG strikes recorded in the study area. In doing so, a total of 2351 significant events were identified.
Convective indices like the CAPE (convective available potential
energy) index and the convective condensation level (CCL) were retrieved from the
Barcelona sounding station (N 41
Synoptic weather patterns (10) linked to lightning activity in Catalonia. Each SWP is characterised by the average sea level pressure (left column), the geopotential at 500 hPa (middle column) and the averaged lightning distribution (right column).
Information on the cause of ignition and the date, time and coordinates of
the point of ignition of the LIWs that occurred in Catalonia was retrieved from
the wildfire database (WFDB), managed by the Forest Protection Agency of the
autonomous government of Catalonia (Servei Prevenció Incendis Forestals, SPIF). According to the SPIF-WFDB
(GENCAT, 2021), in the period of study (2005–2020), there were 561 wildfires a
year in Catalonia, which burned 2732 ha of forest area each year on average.
Large fires (
NCEP/NCAR Reanalysis data (Kalnay et al., 1996) at 2.5
A reliable lightning–wildfire association is first needed to properly describe the natural fire regime. As in Soler et al. (2021) and Pineda et al. (2022), each LIW recorded in the SPIF-WDB was crossed with the SMC-LLS lightning database to designate the most probable candidate (MPC) as the fire starter. The MPC is designated through a “proximity index” (Larjavaara et al., 2005). The proximity index (A) is a combined probability that a lightning with a temporal distance and spatial distance to ignition fire could have caused the wildfire. See Pineda et al. (2022) for details. The MPC is the lightning stroke with the highest score in A. In total, we matched 870 LIWs using A (93 %). The selection of a fire starter allows a date and time of ignition to be assigned to each LIW.
Lightning-related SWPs were obtained using an objective methodology developed by Aran et al. (2011). Firstly, a principal component analysis (PCA) was used to reduce data dimensions (Richman, 1986). The PCA is not only designed as a data reduction technique, but also as a method which ensures that only the fundamental variation modes of the data are considered for the clustering process. The PCA in S mode was applied separately to VT, with relative humidity at 700 and at 850 hPa. In this mode, variables are grid points and days are observations. The scree test was used to determine the number of components involved (Cattell, 1966), and orthogonal varimax rotation was applied to minimise the dependence between the common principal components (Richman, 1986).
Secondly, cluster analysis (CA) was used to determine the main SWPs
associated with lightning activity in Catalonia. CA was applied to the score
matrix formed by scores for each variable retained in the PCA analysis. The
non-hierarchical
Discriminant analysis (DA) was applied to the score matrix retained in the PCA, as independent variables, and a categorical dependent variable, i.e. the cluster label from CA. The DA was used to (i) obtain the discriminant functions from a stepwise selection criterion, Wilks' lambda criterion (e.g. Diab et al., 1991); (ii) validate the SWPs to determine whether CA is an effective tool for predicting the category membership; (iii) use the discriminant functions to reclassify any borderline cases if necessary (Sioutas and Flocas, 2003; Michailidou et al., 2009); and (iv) also to classify and/or predict future events in Catalonia.
The multivariate analysis produced 10 SWPs. To analyse the synoptic features, composites (averages) of all the events included in each SWP were constructed using mean sea level pressure (SLP) and geopotential at 500 hPa (Z500) data grids (Fig. 2). The main characteristics of the resulting SWPs are detailed below, sorted in descending order according to the percentage of events belonging to each type.
SWPs conducive to thunderstorms (TSs) and to lightning-ignited wildfires (LIWs). Columns show event distribution per SWT; LIW events per SWP; 6 h periods with LIWs; 6 h periods with one single LIW, multiple LIWs or five or more LIWs; and average number of LIWs per 6 h period.
The 15 August 2017 12:00 UTC sounding from the Barcelona
sounding station. This episode, corresponding to the SWT-TL category,
produced six LIWs between 12:30 and 17:30 UTC. CAPE: 1837 J km
Considering the scale of dynamics and convection type, the 10 SWPs can be grouped into three broad categories. Thermal dynamics categories (SWT-TL and SWT-NF) contribute 50 % of the total SWPs, mesoscale dynamics (SWT-PC, SWT-IL and SWT-IB) 32 % and synoptic dynamics (IrL, SWT-WML, COL, DB-EF and SL-NF) the remaining 16 %.
All in all, results show that SWPs conducive to thunderstorm activity in Catalonia are mainly related to atmospheric convection caused by a trough oriented from north to south at the 500 hPa level (up to 90 % of the total, namely all categories with SWT). This general configuration leads to thunderstorm development in the region of study, often coupled with severe weather, as stated by other authors (e.g. Aran et al., 2007; Mateo et al., 2008). Additionally, regional differences are modulated by the predominant surface flow.
Once the SWPs supporting thunderstorm activity in Catalonia have been described, we now focus on the SWPs that favour lightning-caused wildfires (summarised in Table 1). Most wildfires belong to SWT-TL (51.4 %). The other SWP with a relevant number of LIWs is SWT-NF (23.7 %) and SWT-PC (13.3 %). Together, these three categories account for 88.4 % of the LIWs. Compared to the distribution for regular thunderstorm events, the frequency of the category SWT-TL is almost doubled in the case of LIWs. Moreover, in 67 % of the SWT-TL episodes, a LIW took place. At the other end, categories like COL, SWT-IB and SL-NF are irrelevant to LIW production.
Convective indices like the CAPE and the CCL derived from the Barcelona radiosonde were used to analyse seven representative episodes, showing multiple LIWs, corresponding to the three SWPs that account for the majority of the LIWs. These three SWPs, which occurred in July and August, are characterised by the arrival of a short-wave trough at 500 hPa. Under these conditions, we have seen that thunderstorms that trigger LIWs during the warmer hours of the day generally develop under relative high CAPE values. Still, a high CAPE is not exclusive of the LIW episodes but to thunderstorm development. Contrarily, the CCL index is more representative, since in five of the seven cases analysed it is above 1000 m and even around 2000 or higher in three cases (Fig. 3 shows an example).
Such conditions would also be conducive to the type of storm described in Pérez-Invernón et al. (2021), who characterised fire-producing thunderstorms by clouds with a high base and by weak updrafts between 300 and 450 hPa pressure levels. Isolated afternoon storms are more likely to generate a fire, compared to storms associated with frontal systems. In fact, small thunderstorm cells produce scattered, low-intensity precipitation, decreasing the likelihood that precipitation would either extinguish the fire or provide enough moisture to inhibit ignition (Larjavaara et al., 2005; Hall, 2008; Soler et al., 2021).
Table 1 shows that in about half of the 6 h episodes, there were multiple LIWs,
sometimes more than five. The three main categories linked to the LIWs had
a similar average of
Lightning-ignited wildfires (LIWs) are not easy to analyse, since the ignition can remain latent from hours to days until the wildfire eventually flares up and is finally reported. Such dating discrepancies pose problems for researchers while establishing the environmental factors involved in LIW processes. The selection of a fire starter, by crossing the wildfire database with a lightning database, allows a date and time of ignition to be assigned to each LIW. Once the origin of the LIW is established, the atmospheric conditions leading to the lightning ignition can be examined. This way, synoptic weather patterns favouring thunderstorm occurrence were analysed, with a particular focus on the thunderstorm episodes that triggered LIWs. The present analysis has shown that LIWs in Catalonia are related to a typical synoptic weather pattern dominated by a short-wave trough at 500 hPa, with three distinct associations: an Iberian thermal low (SWT-TL) (51 % of the LIW), a northern flow (SWT-NF) (24 %) and prefrontal convection (SWT-PC) (13 %).
Further research should investigate how these particular SWPs create a favourable environment for the lightning-caused ignitions to survive and ultimately develop into an active wildfire. Ideally, instability conditions allowing for convection and finally lightning should be coupled with dry air mass at lower levels. Such conditions would also be conducive to high-based clouds, increasing the likelihood that precipitation will evaporate before reaching the ground. The scarcity of moisture in the environment strongly influences the fuel to combust. Being able to identify areas with an elevated lightning-ignition survival at daily timescales would be of great assistance to the forest protection tactical decision-making process, i.e. locating ignitions and potential holdover fires, preparing for days with multiple ignitions or routing detection flight paths.
The meteorological data that support this study were obtained from Servei Meteorològic de Catalunya. Data from the wildfire database were obtained from Servei de Prevenció d'Incendis Forestals by permission. Data shall be requested from the corresponding institution.
All authors contributed to the writing of this paper. NP initiated the study and wrote the first draft. JPC and NPZ prepared the data, designed the methodology, and carried out the classification and the analysis output. MA and XS developed the meteorological interpretation of the results. All authors discussed the results and reviewed the final version of the paper.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This article is part of the special issue “21st EMS Annual Meeting – virtual: European Conference for Applied Meteorology and Climatology 2021”.
We acknowledge the Servei Prevenció Incendis Forestals (Forest Protection Agency) of the Generalitat de Catalunya (autonomous government of Catalonia) for providing access to the wildfire database.
This paper was edited by Gert-Jan Steeneveld and reviewed by Marc Castellnou and one anonymous referee.