Saharan Hot and Dry Sirocco Winds Drive Extreme Fire Events in Mediterranean Tunisia (North Africa)

: With hot and dry summers, the Mediterranean basin is a ﬀ ected by recurrent ﬁres. While drought is the major driver of the seasonal and inter-annual ﬁre distribution in its northern and mildest climate conditions, some extreme ﬁre events are also linked to extreme winds or heat waves. The southern part of the Mediterranean basin is located at the driest range of the Mediterranean bioclimate and is inﬂuenced by Saharan atmospheric circulations, leading to extreme hot and dry episodes, called Sirocco, and potentially acting as a major contributor to ﬁre hazard. The recently created ﬁre database for Tunisia was used to investigate the ± 10-day pre- and post-ﬁre timeframe of daily weather conditions associated with ﬁre events over the 1985–2006 period. Positive anomalies in minimum and maximum temperatures, negative anomalies in air relative humidity, and a preferential south-eastern wind during ﬁre events were identiﬁed, which were characteristic of Sirocco winds. + 7 ◦ C anomalies in air temperature and − 30% in relative air humidity were the critical thresholds for the most extreme ﬁre conditions. In addition, meteorological anomalies started two days before ﬁre events and lasted for three days after for large ﬁres > 400 ha, which suggests that the duration of the Sirocco event is linked with ﬁre duration and ﬁnal ﬁre size. Lastly, the yearly number of intense Sirocco events better explained the inter-annual variability of burned area over the 1950–2006 period than summer drought based on Standardized Precipitation Evaporation Index (SPEI) indices.


Introduction
With hot and dry summers, the Mediterranean basin is affected by recurrent fires shaping landscapes [1,2] and species functional traits [3]. Because of a significant climate gradient from north to south and following the fire varying constraint hypothesis [4], fire hazard in this region is driven by both fuel availability at the southern part of the drought gradient, and the length of the drought period defined by seasonal fuel moisture content at the northern part [5]. This latitudinal gradient follows a constraint gradient, switching from moisture limitation to fuel limitation [6]. In non-fuel-limited ecosystems, recent studies illustrate the importance of summer drought on the inter-annual variation of burned area across the basin [7,8] and its value in fire hazard predictions [9]. However, the strength of the yearly drought/burned area relationship decreases going southward [10], which suggests the potential for spatial heterogeneity in the relative contribution of other keystone climate variables over the Mediterranean basin, as observed across continents within similar biomes [11] or between ecoregions [12]. Recent studies illustrate the significant influence of air temperature on fire activity in Spain [13], where summer drought is recurrent, so that inter-annual variability of heat wave frequency is more important than inter-annual variability of drought length. Most Mediterranean-scale projections of fire danger rely on the Canadian fire weather index (FWI) integrating drought, daily temperature, air humidity, and wind information [14]. More recently, fire events in the Mediterranean basin have been related to general circulation patterns and subsequent weather anomalies around these fire events [15]. Various fire weather types (named wind-, drought-, or heat-driven fire weathers) could then be identified to better characterize the heterogeneity of fire weather types in local fire regimes, and how the relative frequency of these weather types might change over long time-scales and from one year to the other [16]. These analyses were based on daily anomalies of keystone weather variables as temperature, relative air humidity, and wind speed to capture their temporal variability before and during fire events, and they could bring a more detailed understanding of the various weather conditions leading to large fires. We hypothesize that a particularly high frequency of heat waves would impact fire occurrence at the southern rim of the Mediterranean basin based on the identification of these various fire weather types potentially occurring in the Mediterranean basin, and the identified weaknesses on the applicability of a unique fire weather index. Indeed, the northward Sirocco wind is the most well-described climate phenomenon affecting air quality, health, and environment in this region, resulting from a depression centered in central Mediterranean, with northward winds bringing hot and dry air from Sahara [17]. Fires affect forested ecosystems that border arid areas worldwide during these extreme weather episodes and atmospheric circulation schemes blowing abnormally hot and dry winds for several consecutive days. The most documented phenomenon is the Santa Ana event in California with hot and dry winds blowing from the continental desert of Nevada toward the Californian coast and leading to the major fire events in this region [18]. The southern rim of the Mediterranean basin has been little investigated due to the lack of fire data, but it could highly contribute to a full understanding of fire hazard control along the whole Mediterranean gradient. Our study is based on the hypothesis that a process similar to Californian Santa Ana weather conditions, i.e., a shift in wind direction generating short-term weather anomalies, would be an indicator of Sirocco events [19]. The recently published quality-checked fire database in Tunisia for the 1985-2010 period [20] was used (i) to investigate the pre-and post-fire daily weather related to fire size occurrence over a +/−10-day timeframe and (ii) to test how southern winds from Sahara are associated with meteorological anomalies during fire events. The link between the intensity of weather anomalies and fire size was then tested in order to build a yearly burned area model based on the number of Sirocco events that was compared with a drought-based model.

Study Area
The study area covers the northern part of the Tunisian national territory where forest ecosystems occur in the humid Mediterranean and semi-arid bioclimates [20]. The region is mountainous, with elevation that ranges from 0 m on the coast to 1235 m at Mount Chaambi (Atlas Mountains range) ( Figure 1). The regional climate is Mediterranean, with precipitation mostly being distributed in autumn, winter, and spring, while summers are dry, with a strong north/south gradient getting drier upon reaching the Sahara desert in the south of the country. Tunisia covers a total area of 165,000 km 2 , with only 9416 km 2 of forests and shrublands. Quercus suber L. covers 70,000 ha on the northern part of the country where mean annual precipitation is 800 mm, varying between 700 mm and 1500 mm, and mean temperature is 16 • C [21]. Pinus halepensis Miller covers 361,000 ha in the central part of the country where mean annual precipitation is 400 mm [22]. Mediterranean shrublands cover 328,000 ha and they are composed of Arbutus unedo L., Erica arborea L. Cistus spp. L, Phyllirea latifolia L., and Pistacia lentiscus L. on the northern coastal part and Rosmarinus officinalis L., Cistus spp. L., and Lavendula spp. L. in the central part [23].

Fire Data
Daily fire data were collected from the Tunisian Forest services, starting in 1985 (Source Direction Générale des Forêts, Tunisie), and being checked and corrected with remote sensing [20]. For missing fires in the national database that were captured with Landsat images, burn date was determined by the MCD14ML global hotspot database from the MODIS sensor [24]. The exact burn date of these missing fires could not be identified due to the lack of MCD14ML data before 2000, and they were removed from the analysis. Each fire is defined by its ignition date, end date of the fire (so we could derive fire duration), location, burned area, and vegetation type. For this analysis, wildland fires (affecting forests and shrublands) were selected (so that agricultural burnings were removed) and then classified into small fires (>1 ha), medium (>100 h), and large fires (>400 ha), with 968, 94, and 23 fires, respectively, in each fire class.

Climate Data
Daily climate data from the Tunis meteorological station (latitude = 36.75° N, longitude = 10.25° E), being representative of the general circulation of the country, were provided by INM (Institut National de la Météorologie, Tunis) for the period 1950-2006. The following variables were used: maximum temperature Tmax (°C), minimum temperature Tmin (°C), air relative humidity HR (%), wind speed WS (m·s −1 ), and wind direction WDIR (in degrees of deviance from North). Climate reanalysis can be used as an alternative source of data [16], but our analysis covered the period back to 1950 with this long time series on temperature, and tested whether local weather data could be used for future early alert systems.
From these climate variables and daily precipitation, a drought index (DI) was computed from the drought code (DC) component of the fire weather index (FWI) [25]. This index has been developed for Canadian forests as an index of litter moisture, and further used globally and in the Mediterranean basin as a proxy for plant fuel moisture [26]. DI is 0 for the driest conditions and it reaches a maximum value of 1 when the soil water content is at field capacity.
As a global drought index reference, and for comparison with other studies, the Standardized Precipitation Evaporation Index (SPEI) was downloaded from the global SPEI database [27] in order to test for the relative importance of summer drought on the inter-annual burned area. SPEIs for the month of September (end of the dry season) and for the preceding timeframes of one to seven months

Fire Data
Daily fire data were collected from the Tunisian Forest services, starting in 1985 (Source Direction Générale des Forêts, Tunisie), and being checked and corrected with remote sensing [20]. For missing fires in the national database that were captured with Landsat images, burn date was determined by the MCD14ML global hotspot database from the MODIS sensor [24]. The exact burn date of these missing fires could not be identified due to the lack of MCD14ML data before 2000, and they were removed from the analysis. Each fire is defined by its ignition date, end date of the fire (so we could derive fire duration), location, burned area, and vegetation type. For this analysis, wildland fires (affecting forests and shrublands) were selected (so that agricultural burnings were removed) and then classified into small fires (>1 ha), medium (>100 h), and large fires (>400 ha), with 968, 94, and 23 fires, respectively, in each fire class.

Climate Data
Daily climate data from the Tunis meteorological station (latitude = 36.75 • N, longitude = 10.25 • E), being representative of the general circulation of the country, were provided by INM (Institut National de la Météorologie, Tunis) for the period 1950-2006. The following variables were used: maximum temperature T max ( • C), minimum temperature T min ( • C), air relative humidity HR (%), wind speed WS (m·s −1 ), and wind direction WDIR (in degrees of deviance from North). Climate reanalysis can be used as an alternative source of data [16], but our analysis covered the period back to 1950 with this long time series on temperature, and tested whether local weather data could be used for future early alert systems.
From these climate variables and daily precipitation, a drought index (DI) was computed from the drought code (DC) component of the fire weather index (FWI) [25]. This index has been developed for Canadian forests as an index of litter moisture, and further used globally and in the Mediterranean basin as a proxy for plant fuel moisture [26]. DI is 0 for the driest conditions and it reaches a maximum value of 1 when the soil water content is at field capacity.
As a global drought index reference, and for comparison with other studies, the Standardized Precipitation Evaporation Index (SPEI) was downloaded from the global SPEI database [27] in order to test for the relative importance of summer drought on the inter-annual burned area. SPEIs for the month of September (end of the dry season) and for the preceding timeframes of one to seven months were extracted for the 0.25 • grid cell of the Tunis station.

Statistical Analysis
For each fire event, daily weather conditions for all climate variables were extracted from 10 days before to 10 days after the fire. Daily weather conditions over the same timeframe were also extracted for each day over the fire season from July to September. The mean values and uncertainties for each climate variable were plotted on a temporal scale. The daily uncertainty was assessed through a bootstrapped subsampling with 1000 replicates due to a non-normal distribution of climate variables (Shapiro-Wilk W test, with a significance threshold fixed at 0.05, shapiro.test function in 'stats' Rcran package [28]). For each replicate, the mean value of climate variables on 10 events randomly sampled within the whole dataset was calculated. The upper and lower confidence levels were calculated with the 'quantile' function, respectively, at 0.9 and 0.1 significance threshold applied to the 1000 means obtained for each replicate.
From this analysis, the duration of weather anomaly was derived as the number of days when T max was above the upper confidence interval of the mean T max over the summer season for each fire size class. Plots of fire size (log10 transformed) as a function of the duration of the fire (as registered in the forest service database) and the duration of the weather anomaly were performed with the 'ggplot2' R cran package. A log10 model fitting was performed with a locally weighted least square regression 'loess' option in the 'geom_smooth' function with confidence interval being fixed at 0.95.
As a final step, the correlations between yearly burned area and summer drought or number of weather anomaly events (that were attributed to Sirocco meteorological phenomenon) during the fire season were evaluated. Summer drought was based on SPEI_september indices. The number of Sirocco events during the fire season was based on either a direct temperature anomaly as a difference with the mean summer T max value of 33 • C (Nsir_Tmax) or the abrupt temporal variation of T max between the day = i and day = i − 3 (Nsir_Tdiff). The effects of drought and Sirocco intensities were evaluated, respectively, from a one-month to a seven-month timeframe, and a + 1 • C to a + 7 • C T max anomaly. Correlations between yearly burned area and meteorological indices were performed with the Pearson test and the significance level was determined with a Monte Carlo randomization on 1000 replicates using the 'meancorr' function of the 'synchrony' R cran package [29]. The Monte Carlo randomizations were performed by independently shuffling the columns of matrix, thus destroying the potential temporal autocorrelation of time series. To capture only the qualitative synchrony of high and low burned area extremes (and not the burned area/weather relationship) with climate indices, the proportion of local minima/maxima common to both time series was also computed. The upper and lower peaks in each time series were identified when the derivative values changed signs. The significance level was estimated with the same Monte Carlo randomization as the correlation, by using the 'peaks' function of the 'synchrony' R cran package [29].
A linear model was finally applied to establish the relationship between annual burned area and annual climate indices of drought or the number Sirocco events (NSir_Tmax and NSir_Tdiff) using the 'dynlm' R cran function designed for time series regression. The modeled burned area was then simulated using the 'predict' function that was applied to the T max time series over the 1950-2006 period. Before analysis, the annual burned area and the number of Sirocco events were log10 transformed due to a non-normal distribution (p value < 0.05, Shapiro-Wilk test). The SPEI indices were normally distributed (p value > 0.05, Shapiro-Wilk test) and were kept non-transformed.
Wind roses for the whole summer season and for two fire size classes (fire size > 1 ha, and >400 ha) were built for the climate variables T max , wind speed WS, and air relative humidity HR, in order to illustrate the preferential directions of fire events and their associated weather conditions. Wind roses were built with the 'ggplot2' R cran Package. Figure 2 represents the mean and confidence interval for each climate variable during the 10 days before (−10) and after (+10) fire events and for the three fire size classes as compared to the whole fire season.

Time-Lagged Fire/Weather Analysis
For minimum temperature T min (Figure 2a), the mean value during the fire season (with or without fires) was 21 • C (±1 • C). For small fires (>1 ha), a positive non-significant anomaly of 22 • C(±1 • C) one day before the fire event was observed then reaching +23 • C (±1 • C) for fires > 100 ha and 24 • C (±1.5 • C) for fires > 400 ha. This anomaly was only observed for day (−1) for small fires, from day −1 to day +1 for medium fires > 100 ha, and from day −1 to day +3 for large fires > 400 ha.
For maximum temperature T max (Figure 2b), the mean value during the fire season was 33 • C (±1 • C). For small fires, the positive anomaly to 34 • C (±1 • C) during fire events was not significant, reaching significant differences of 37.5 • (±1 • C) and 38 • C (±2 • C) for medium and large fires, respectively. The anomaly was observed from day − 1 to day + 1 for small fires, from day − 2 to day + 3 for medium fires, and from day − 2 to day + 5 for large fires, respectively, leading to anomaly events prolonged over three, six, and eight days.
These positive temperature anomalies were associated with negative air humidity (HR) anomalies ( Figure 2c). The mean HR over the summer was 62% (±3%). For small fires, a non-significant negative anomaly down to 60% (±2.5%) was observed for day − 1 and during the fire event, significantly decreasing to 55% (±3%) and 52% (±2%) during medium and large fire events. A HR negative anomaly occurred from day − 2 until day + 1 for medium fires, and from day − 3 until day + 4 for large fires, respectively, lasting two, four, and eight days.
During these dry and hot events, wind speed was slightly higher than the mean summer value (Figure 2d). Wind speed over the summer season was 8.5 (±1) m·s −1 , not being significantly different during small and medium fires, but slightly significantly higher, up to 10 (±1) m·s −1 , from day − 1 to day + 3 for large fires.
All of these climate variable anomalies were associated with a significant change in wind direction (Figure 2e). The dominant incoming wind direction in the region was NW (245 • ± 10 • ). This wind direction progressively deviated to south according to fire size, reaching 235 • for small fires, 205 • (SW) for medium fires, and 175 • (SSE) for large fires. This wind direction anomaly started during the fire event and lasted until day + 4 for large fires, and until day + 3 for medium fires.
The drought conditions during fire events were longer lasting and did not show any specific temporal anomalies during or around the fire events (Figure 2f). The mean DI during the fire season was 0.4 (±0.05), decreasing to 0.35 for small and medium fires and to 0.3 for large fires, suggesting slightly dryer conditions around fire events than during the rest of the dry season.
(e) (f) Figure 2. Mean daily weather values from 10 days before to 10 days after fire events for fire size classes >1 ha (blue), >100 ha (yellow), >400 ha (red), and all days during the July-August fire season (black), for daily climate variables of (a) minimum temperature (Tmin in °C), (b) maximum temperature (Tmax in °C), (c) relative air humidity (HR in %) (d) wind speed (WS in m·s −1 ), (e) wind direction (WDIR in degree deviation from North), and (f) drought index (DI) varying from 0 (completely dry) to 1 (full moisture) soil water conditions. Colored shaded areas correspond to the 90% confidence interval obtained from 1000 bootstrapped subsamples.
For minimum temperature Tmin (Figure 2a), the mean value during the fire season (with or without fires) was 21°C (±1 °C). For small fires (>1 ha), a positive non-significant anomaly of 22 °C(±1 °C) one day before the fire event was observed then reaching +23 °C (±1 °C) for fires > 100 ha and 24 °C (±1.5 °C) for fires > 400 ha. This anomaly was only observed for day (−1) for small fires, from day −1 to day +1 for medium fires > 100 ha, and from day −1 to day +3 for large fires > 400 ha. Mean daily weather values from 10 days before to 10 days after fire events for fire size classes >1 ha (blue), >100 ha (yellow), >400 ha (red), and all days during the July-August fire season (black), for daily climate variables of (a) minimum temperature (T min in • C), (b) maximum temperature (T max in • C), (c) relative air humidity (HR in %) (d) wind speed (WS in m·s −1 ), (e) wind direction (WDIR in degree deviation from North), and (f) drought index (DI) varying from 0 (completely dry) to 1 (full moisture) soil water conditions. Colored shaded areas correspond to the 90% confidence interval obtained from 1000 bootstrapped subsamples.

Fire Weather Frequency and Fire Duration
Wind roses of maximum air temperature (T max ), wind speed (WS), and air relative humidity (HR) are presented in Figure 3 for all days during the summer season, and during small and large fires, to better understand the frequency distribution of fire weather during fire events. The dominant wind directions in the area were NW (285−315 • ), followed by NNE (15-45 • ) and ESE (105-135 • ). NW winds were usually cooler with 33% of days with T max < 20 • C, while ESE winds were warmer with 50% of days >30 • C with some extremes >40 • C. Regarding HR, all directions were quite similar, with 15% of days reaching HR > 80%, 40% of days 80% > HR > 70%, 40% of days 70% > HR > 60%, and 5% with HR < 60%. Wind speed was more related to wind direction. The fastest winds came from the NW direction with 50% of days reaching >10 m·s −1 , while 70% of days with wind direction NE to SE experienced wind speed <10 m·s −1 . We could then differentiate NW winds as cooler and faster as opposed to SE winds as hotter and slower in Tunisia during the fire season.   For small-fire days (>1 ha), the frequency of NW winds substantially decreased, so that NNE and ESE winds became dominant, with T max mostly > 30 • C and HR < 80%. Wind speed was faster in the ESE direction (70% of events > 8 m·s −1 ) than NNE (80% of events < 8 m·s −1 ). This result suggested that NW winds were less fire-prone as a consequence of cooler temperatures, even under faster wind speeds, and fires in Tunisia were mostly driven by warmer air temperature events.
Large fire events occurred almost only (75%) during E to SE wind directions. T max was always >30 • C with 25% of experiencing T max > 40 • C. Around 50% of days had HR < 60%, with the driest days (HR < 50%) occurring with ESE wind direction. Wind speed remained mild between 8 and 12 m·s −1 .
In turn, this analysis illustrated how a selective fraction of weather conditions were associated with fire events in Tunisia. It also confirmed the occurrence of ESE hot and dry winds in 75% of large fire events, complementing the mean fire/weather analysis previously performed, and characterizing the combination of weather anomalies occurring during Sirocco events. Figure 4 represents fire size (Y axis, log10 scale) as a function of fire duration in days. The two variables were significantly correlated (R 2 = 0.62), indicating that the final fire size depends on the duration of the fire. Particularly, fires smaller than 100 ha lasted for one or two days, while fires larger than 100 ha lasted from three to 12 days with a slight increase after that threshold. Plots of fire size as a function of Sirocco duration obtained from Figure 2b showed a similar relationship, with short Sirocco events (zero to three days) leading to small fires, and longer Sirocco events leading to fires >100 ha. Sirocco duration, fire duration, and final fire size were all highly interrelated, which suggested that yearly Sirocco event frequency and duration might drive the annual burned area.   . Relationship between fire size (Y-axis, Log10 scale in ha) and fire duration (X axis in days) over the 1985-2006 period in Tunisia (grey dots represent fire events and red dots the mean size for each duration with standard deviation intervals). The log10 regression fitted on mean values is also represented (Red line) with its 90% confidence interval (shaded red area). The duration of Sirocco events (in days) during fire size classes >1 ha, >100 ha and >400 ha is also represented (blue dots and 95% confidence interval bars). Figure 5 illustrates the burned area simulated by the number of Sirocco days based on temperature anomalies (Nsir_Tmax) for the 1950-2006 period (blue line) based on the previous results relating fire size and Sirocco duration. The yearly burned area over the 1950-2006 period was simulated with R = 0.44 (p value < 0.001) based on the number of Sirocco events defined by the highest maximum temperature anomaly (+7 • C). Correlation coefficients were lower, but still significant (p values < 0.01) for anomalies of +5 • C and +6 • C ( Table 1). The fraction of synchronized lower/upper extremes between burned area and Nsir_Tmax reached 0.45 for the highest temperature anomalies and decreased to 0.24 for the other significant correlations. However, the simulated burned area increased over time as a result of increasing temperatures over the region, leading to increasing model residuals in the recent years ( Figure 5c). As an alternative, the burned area was modeled by the number of Sirocco days based on the relative maximum temperature anomaly between the day of the fire event and three days before (Nsir_Tdiff). The burned area over the 1950-2006 period was simulated with a correlation coefficient R = 0.40 (p < 0.01) for a +7 • C anomaly, a slightly lower performance than the previous index, but better capturing the recent decrease in burned area in the context of increasing average air temperature. The model residuals of the observed burned area in recent years were lower than the previous model ( Figure 5c). The correlation coefficients decreased to 0.3 (p < 0.01) for temperatures anomalies ranging for +4 • C to +6 • C. The fraction of synchronized lower/upper extremes between burned area and Nsir_Tdiff reached 0.38 for the highest temperature anomaly, and decreased to 0.24 for other significant correlations. In turn, the yearly burned area in Tunisia was mostly correlated to the number of extreme Sirocco events. Sirocco events, defined as hot and dry air from the south, were correlated to absolute daily maximum air temperature, but with a temporal trend more particularly related to temperature variations within a short period of time. Besides the significant effect of Sirocco events, none of the SPEI indices covering one to seven months before September were significantly correlated to annual burned area ( Table 1). The fraction of synchronized low/upper extremes between burned area and SPEI was around 0.2, with the highest and significant (p < 0.05) value of 0.36 reached for SPEI_3 corresponding to the July-September period. Table 1. Correlation coefficients (and significance level 1 ) and fraction of low/upper synchrony between yearly burnt area anomaly and fire hazard indices modelled from the number of Sirocco events based on temperature anomalies three days before and during fire events (Nsir_Tdiff), temperature anomaly above the 33 • C threshold (Nsir_Tmax) and the drought index Standardized Precipitation Evaporation Index (SPEI) (September). Levels indicate the temperature anomalies (in • C) or the number of months used for SPEI indices.  Correlations between observed and simulated burned area anomalies for the temperature anomaly (red dots) and temperature threshold (blue dots) (c) Yearly model residuals with a temporal locally weighted least square regression and 0.9 confidence interval (shaded areas) calculated frim 'loess' option in the ggplot2 'geom_smooth' function. Table 1. Correlation coefficients (and significance level 1 ) and fraction of low/upper synchrony between yearly burnt area anomaly and fire hazard indices modelled from the number of Sirocco events based on temperature anomalies three days before and during fire events (Nsir_Tdiff), temperature anomaly above the 33 °C threshold (Nsir_Tmax) and the drought index Standardized Precipitation Evaporation Index (SPEI) (September). Levels indicate the temperature anomalies (in °C) or the number of months used for SPEI indices.   Correlations between observed and simulated burned area anomalies for the temperature anomaly (red dots) and temperature threshold (blue dots) (c) Yearly model residuals with a temporal locally weighted least square regression and 0.9 confidence interval (shaded areas) calculated frim 'loess' option in the ggplot2 'geom_smooth' function.

Heat Waves and Fire Events in the Mediterranean Basin
The results showed that anomalous dry and hot days significantly promoted large fire occurrence in Tunisia and these fires preferentially happened during the dry season under the driest conditions. Dry conditions were necessary for fires to happen, but not sufficient, as they had to be associated with combined dry and hot daily conditions. In Tunisia, fuel moisture conditions can be regionally variable in the forested wettest northern part of the country according to the climate gradient and topographical features [21,30,31]. In addition, shallow-rooted shrubland species and the herbaceous litter fuel at the surface layer comprise a large part of natural vegetation, and may be even more susceptible to frequent summer droughts [22,23,32]. Over the whole country, long-term drought reconstructions that are based on tree-rings revealed that recent decades appear as among the driest [33], and confirmed recurrent droughts in the studied period.
This study, which was the first one on the southern part of the Mediterranean basin, also confirmed previous studies performed in the European Mediterranean region, where annual burned area in the southernmost regions was least correlated to drought inter-annual variability [10]. In these regions, heat waves may be the most important climate-related feature for fire hazard. In Spain, for example, heatwaves have been the strongest drivers of fire occurrence in recent decades [34][35][36], as well as in Sardinia [37]. Extreme recent fire events in Portugal [38] and Greece [39] were also associated with heatwaves, as was southern France during the exceptional 2003 heat wave [40].
The plant ecophysiology underlying the important influence of heat waves on fire hazard could be leaf temperature and high evaporative demand desiccating leaves and low air humidity increasing heat transfer and decreasing energy used for water vaporization [41]. Recent studies illustrated a decoupling between photosynthesis and transpiration during heatwaves in order to maintain enough transpiration for thermal cooling, to keep leaf temperature below the limits for maintaining physiological functions at the expense of higher leaf desiccation [42]. Further studies should investigate plant functioning during heat waves under dry conditions in the Mediterranean to refine drought-driven fire weather indices.

Synoptic Conditions Associated to Sirocco Events and Climate Projections
In Tunisia and North Africa as a whole, hot and dry weather has been mostly associated with northward winds coming from Sahara [43,44]. This weather condition is called Sirocco and it is a recurrent regional atmospheric phenomenon that highly affects agriculture [45] and air quality with subsequent consequences on health and human activities [17,[46][47][48] and regional climate [49].
Meteorological characteristics of Sirocco events in Tunisia have been fully described as mostly occurring between May and August, during the dry season [45]. With an average of 13.3 Sirocco events per year between 1970 and 2004, they are subject to high inter-annual variability, from four events in 2004 to 33 events in 1994 in the city of Sfax (central eastern coast of Tunisia). Sirocco events are associated with slightly faster winds, low air humidity below 20%, and an increased maximum temperature of up to 7 • C, while nighttime minimum temperatures are similar to non-Sirocco days. Overall, these conditions lead to almost a doubling in potential evapotranspiration rate from 5-7 mm·day −1 to 10-11 mm·day −1 during Sirocco events. The duration of a Sirocco event can be short day-long events (60% of cases), but can also last from two to six days, particularly during the dry season. Our results are in full agreement with these observations, both regarding the frequency and intensity of the weather anomalies. The Sirocco climatological phenomenon has been described as resulting from three major atmospheric circulation types [45]: (i) when a depression occurs west of the Mediterranean basin and an anticyclone is found east (or either conditions separately) and leads to a northward wind corridor over Tunisia, (ii) when an anticyclonic situation occurs on the southeastern part of Tunisia and generates northward wind towards the depression in the northern part of the country, and (iii) when a depression circulates from west to east on the coastal area of North Africa. All of these situations bring dried and heated air passing over the Sahara desert toward the coastal northern part of Africa [43][44][45]. This Mediterranean coastal fringe of North Africa adapted to mild temperatures is then particularly highly susceptible to these unusual atmospheric conditions due to its geographic position at the front of these arid regions. The northern part of the Mediterranean basin is less affected by Sirocco winds bringing hot and dry conditions, because air humidity increases and temperature decreases when crossing the Mediterranean Sea. However, western Sirocco events crossing the Mediterranean sea through continental Spain keep their most extreme dry and hot conditions over southern Europe, as evidenced by the recent heatwaves in 2018 and 2019 [50]. General circulation analysis related to Sirocco events illustrates its link to the Northern Atlantic Oscillation (NAO) [51,52] as a synoptic index.
The synoptic conditions observed in North Africa during Sirocco events are actually very similar to the Santa Ana winds in California [18]. Continental dry and hot winds coming from inland desert blow over the coastal Mediterranean California and lead to extreme fire events [19,[53][54][55]. Similarly, southeastward winds coming from the central Australian desert highly influence southwestern Australia's Mediterranean climate, leading to similar dry and hot conditions, enhanced by Foehn desiccating effect when crossing the Australian Alps [56][57][58]. These conditions lead to the recent major fire events affecting the country. A similar influence of inland dry and hot winds on fire hazard was observed in the Mediterranean ecosystems of central and south central Chile [59], thus confirming the global propensity of Mediterranean ecosystems bordering hot and dry arid areas to be highly influenced by these dry heat waves.
In turn, most Mediterranean regions worldwide may be highly susceptible to arid weather conditions suddenly blowing over these ecosystems due to their situation at the interface between temperate and arid climates. The climate change-induced heat waves, demonstrated to be overall affected by global warming, may however be moderated or enhanced by local processes [60]. These local processes can lead to contrasted feedbacks, as a reduction in Santa Ana winds in California [61,62], or an increase, as in Australia [63]. Earth system models used in the fifth model intercomparison project (CMIP5) predict that heatwave amplitude, a proxy of heatwave intensity, will increase in the Mediterranean basin [13,[64][65][66][67], but the frequency of short term Sirocco events has been poorly investigated. Wind fields in climate models also remain an uncertain variable [68].

Future Relevant Indices for Fire Hazard Assessment and Modelling
Fire hazard in the Mediterranean basin and worldwide is mostly assessed through the Fire Weather Index, which combines drought, temperature, wind velocity, and air humidity, so that global and regional fire hazard trends are mostly based on this index [14,[69][70][71]. The drought component has been shown to properly describe the drought season through the fuel moisture status critical for allowing fire ignition and spread [26], and investigated from field measurements [72] and remote sensing [21]. This index has been discussed and modified for the Mediterranean area [73]. A similar approach combining an empirical drought code with weather variables is used in the fire modules of global dynamic vegetation models [74,75]. However, recent studies have also shown the potential regional [76] or temporal [38] variability in the relative contribution of each meteorological variable to properly model actual fire hazard due to changes in land cover, firefighting strategies, or anthropogenic factors. With no aerial intervention plans and only recently settled fire stations, Tunisia can be considered as much less equipped than European countries for fighting fires, so that all fire-prone climate conditions can actively contribute to fire hazard. European countries have increasingly managed to prevent large fires under non-windy conditions [40]. In turn, fire hazard models that are designed for Europe might fail on the southern shore of the Mediterranean basin. Here, our results showed that, in addition to recurrent drought with low inter-annual variability, heat waves during Sirocco events (that associate hot, dry, and windy conditions) were the most favorable conditions for large fires. The fire size was proportional to the duration of this event, as reported in California during Santa Ana events [53] and in the Mediterranean basin [15]. We concluded with a strong correlation between the timeframe of fire spread and the persistence of the fire-prone synoptic conditions.
The Hot-Dry-Windy (HDW) index has recently been developed and tested in the USA to better capture climate features that are associated with large fire events [77,78] as a combination of surface and high atmosphere synoptic and meso-alpha-scale conditions contribution to dangerous fire behavior, when compared to other indices. Local scale indices for early warning fire danger or fire hazard modelling could be developed for North Africa based on simple Sirocco events predictions. In our results, temperature alone can be misleading. Increasing temperature in the region did not lead to increasing fire activity in the last decade, but the extreme maximum temperature increases fairly captured Sirocco events and the subsequent burned area in Tunisia. Such events should also drive inter-annual variability in neighboring countries, such as Algeria [79] and Morocco [80], as well as Middle East countries [81,82] with similar arid/Mediterranean interface.

Conclusions
The regional effect of Saharan Sirocco events was shown to strongly influence the inter-annual variability of burned area in Tunisia. These atmospheric circulation events lead to a northward hot and dry wind coming from Sahara desert, modifying for three to 10 days the current weather, already hot and dry during the summer fire-prone season in the region. However, the increasing mean temperature in the last three decades did not lead to increasing burned area, so that fire hazard was more related to the number of yearly Sirocco events defined by short temporal scale anomalies in temperatures and air humidity than absolute temperature itself. In turn, Sirocco events in future climate projections should be better quantified with a standard index for their efficient characterization. Fire hazard modelling that is based on absolute temperature, such as the Fire Weather Index, should be evaluated in this region, as it might overestimate projected fire hazard and burned area. Further studies should also be devoted to plant functioning for desiccation processes during these particular events, and the subsequent impact on combustion processes. Funding: This research was funded by EU FP7 FUME project, contract grant N • 243888 and French observatory for Mediterranean environment OSU "OREME".