Synoptic weather patterns conducive to lightning-ignited wildfires in Catalonia

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 %).


Introduction
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 . Wildfires generally occur under extreme meteoro-logical 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).
N. Pineda et al.: Synoptic weather patterns conducive to lightning-ignited wildfires 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 lightningignited wildfires and precipitation (Hall, 2007(Hall, , 2008Dowdy 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  and strong wind episodes in Catalonia . Similarly, it was used to analyse SWPs leading to heatwaves in the city of Barcelona . 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.

Area of study
The study area was Catalonia, a region of ∼ 32 000 km 2 located in the north-east Iberian Peninsula and in the northwest of the Mediterranean basin (Fig. 1a). Air flows through the region are modified by local effects, mainly due to two geographical features: the Pyrenees range (∼ 3000 m a.s.l.), creating a natural border to the north, and the Mediterranean Sea to the east (Fig. 1b). The spatial variability of rainfall is high; annual precipitation values range from 300 mm in the south to 1200 mm in the Pyrenees (Martín-Vide et al., 2008). Annual average temperatures also vary greatly due to the considerable altitude range, from an average of 17 • C on the southern coast to about 0 • C along the Pyrenees. In terms of thunderstorm activity, Catalonia and its surroundings are among the regions with moderate lightning activity within western Europe, with an average lightning flash density of ∼ two flashes km −2 yr −1 (Anderson and Klugmann, 2014;Poelman et al., 2016).
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 ∼ 12:00 UTC and a maximum at 15:00 UTC, followed by a slow decrease (Pineda et al., 2011;Aran et al., 2015). Contrarily, the daily cycle of the Balearic Sea is more homogeneous, with minimum activity at midday and maximum activity in the evening and first night hours. The warm waters of the Mediterranean at the end of the summer and beginning of autumn, coupled with favourable thermodynamic conditions for convection, ensure sustained moisture supply, favouring long-lasting thunderstorms over the Mediterranean Sea (Pineda et al., 2009;Kotroni and Lagouvardos, 2016).

Meteorological data
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 • 23 4.08 E 2 • 7 3.36 -WMO code 08190).

Reanalysis data
NCEP/NCAR Reanalysis data (Kalnay et al., 1996) at 2.5 • resolution were used to classify SWPs on a 6 h period basis. The gridded area covered was 30 to 70 • N and 30 • W to 30 • E. Variables chosen in this study were selected to discriminate, as far as possible, convective and non-convective events. Firstly, the vertical totals index, defined by VT = T850-T500, was chosen as it accounts for the static stability (Miller, 1972). Secondly, relative humidity at 700 and at 850 hPa were selected since they consider the water content, which is necessary for assessing thunderstorms associated with potential instability (Haklander and Van Delden, 2003). We used the normalised anomalies of variables, subtracting from each grid point value the temporal average of the grid points separately for each day and each 6 h and corrected daily anomalies at each grid point by the square root of the cosine of the latitude.

Lightning-ignition candidate selection
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.

Synoptic classification
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 K-means method (McQueen, 1967) was used as the clustering algorithm. To determine the number of clusters, the elbow method (Thorndike, 1953cited by Zhang et al., 2016 was applied to a hierarchical Ward cluster (Ward, 1963), which gives the number of initial groups for K means .
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.

Results
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.

Synoptic weather patterns
i. Short-wave trough-thermal low (SWT-TL). This is the most frequent SWP (28 % of the 2351 6 h events). It is characterised by the presence of an Iberian thermal low and a short-wave trough at 500 hPa. Thunderstorms under this synoptic pattern mostly affect the Pyrenees and the Pre-Pyrenees area, during summer months and between 12:00 and 18:00 UTC. Southerly flows provide moisture from the Mediterranean and the Pyrenees the orographic lift.
ii. SWT-northern flow (SWT-NF) (21.2 % of the events). This SWP is defined by an Atlantic short-wave trough associated with thermal waves at 500 hPa level. An anticyclonic ridge reaches Catalonia at the surface. It is frequent from June to September, and thunderstorms take place mostly between 12:00 and 18:00 UTC, the spatial distribution also showing a maximum in the Pyrenees.
iii. SWT-prefrontal convection (SWT-PC) (20.0 %). The SWP shows a short-wave trough at the 500 hPa level with relatively low pressures over Catalonia, especially on the south coast. Such a SWP takes place throughout the year, except for winter. It mostly affects the Pyrenees, the Pre-Pyrenees and the Catalan Pre-Coastal ranges. The SWT-PC is characterised by a bimodal dis-  iv. SWT-Iberian low (SWT-IL) (8.9 %). A trough is located north-west of the Iberian Peninsula at 500 hPa level. At the surface, low pressures affect the entire Iberian Peninsula. Moreover, a relative low in the Balearic Sea leads to easterly winds over the Catalan coast. Lightning activity is mostly located in the south and western parts of Catalonia. It takes place during spring and late summer in the mountain areas but also during autumn on the southern coastal areas.
v. Ireland low (IrL) (8.5 %). The mid-tropospheric trough is associated with a low centred in Ireland. Low pressure at the surface extends from Great Britain to the north-east Iberian Peninsula. This SWP favours lightning activity throughout Catalonia, with a maximum on the coast and in the neighbouring mountain ranges and a minimum in the north-western part of Pyrenees. This SWP is bimodal, with maximum activity during spring and autumn.
vi. SWT-western Mediterranean low (SWT-WML) (4.1 %). The trough at 500 hPa is located in central Europe until the north of Catalonia. Western circulation at mid-levels predominates over the region. At the surface, a low located in Italy reinforces the north-east flow over the north Catalan coast and Pyrenees. This type is predominant in spring.
vii. SWT-Balearic low (SWT-IB) (3.9 %). The configuration of this SWP is characterised by a surface low pressure centred at the Balearic Sea. Lightning activity concentrates along the coastline and Pre-Coastal range, with a maximum located on the central south coast. This SWP mainly is most common in September and early October.
viii. Cut-off low (COL) (3.4 %). In the average composite map of the synoptic types (Fig. 2), it appears that the trough at 500 hPa is disconnected from the main low and deepening in the south Iberian Peninsula. At the surface, a low is located in the Balearic Sea. This SWP is concentrated in October and November and mainly affects the Catalan coastline.
ix. Diffluent block-eastern flow (DB-EF) (1.6 %). At the mid-tropospheric level, the circulation has an omega block configuration over the British Isles. A cut-off low-pressure system is observed over the south-west of Europe. At the surface, in the study area, the synoptic circulation is defined by an eastern flux linked to the anticyclone centred in the Great Britain and the low-pressure centre in the north of the African continent. Under this SWP, lightning activity affects the Pre-Pyrenees, Pyrenees and, locally, the centre and south sectors along the coast. The DB-EF is a residual type that can occur throughout the year and has no preference of time of day.
x. Scandinavian low-northern flow (SL-NF) (0.4 %). This SWP is characterised at the surface by a strong north flow over western Europe due to the anticyclone centre in the Atlantic and two depressions, one of them located in the Scandinavian Peninsula and the other one at the south of the Alps range, which keeps the northerly flow at mid-latitudes. At the mid-tropospheric level, the configuration is defined by a trough in high latitudes. Lightning activity is likely to occur in September and October.
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.

Synoptic weather patterns related to lightning-ignited wildfires
Once the SWPs supporting thunderstorm activity in Catalonia have been described, we now focus on the SWPs that favour lightning-caused wildfires (summarised in 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, lowintensity 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 ∼ two LIWs per 6 h period. From an operational perspective, simultaneous fires impede the coordination of the fire extinction brigades. So, especially in these circumstances, prediction is particularly important. In addition, lightning ignitions can occur in remote locations, and therefore these complex LIW episodes may overwhelm fire brigades, resulting in longer response times and more difficult firefighting campaigns (Podur et al., 2003;Wotton and Martell, 2005;Morin et al., 2015).

Final remarks
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 lightningcaused 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 decisionmaking process, i.e. locating ignitions and potential holdover fires, preparing for days with multiple ignitions or routing detection flight paths. Data availability. 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.
Author contributions. 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.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement.
This article is part of the special issue "21st EMS Annual Meeting -virtual: European Conference for Applied Meteorology and Climatology 2021".