40-Year Fire History Reconstruction from Landsat Data in a Mediterranean Area of North Africa Following International Standards

1. Introduction

Fire is an intrinsic process of the Earth system [1], which is projected to increase in burned areas (BAs) owing to climate change and escalating anthropogenic effects [2]. Among all forest disturbances, fire is the major forest destructive agent in the Mediterranean Basin [3], where heat wave-related fires may increase in the future [4]. When looking across the African side of the Mediterranean Basin, Algeria is the main fire hotspot [5]. Throughout history, this country has witnessed an unprecedented series of large extreme fires [6], with a record-breaking heatwave in the summer of 2023, which affected ecological and socio-economic assets [7,8], including human casualties. These fires may seriously degrade forest habitats in this country, a large part of which may not be restored [9]. However, accurate and spatially explicit BA data about this region, as in several middle-eastern and north-African countries, are scarce or even lacking [10].

Remote sensing has become the most efficient tool for addressing all fire management aspects, including the generation of BA products [11]. Unlike ground-based fire datasets that are often biased or incomplete, or exhibit inconsistencies [12,13,14], satellite-derived BA products provide spatially and temporally consistent and reliable information about fires on regional and global scales [15]. In practice, several pixel-level BA products have been employed in a wide range of research works, including global fire trends [16], characterization of fire regimes [17,18], climate impact on fire patterns [19], fire emission modelling [20] and fire model benchmarking [21], and to derive global databases of single fire events [22].

Early global BA products were based on the coarse resolution data from the SPOT-VGT, ERS2-ATSR2, ENVISAT-AATSR, NOAA-AVHRR, PROBA-V and MODIS sensors [15]. In the last few years, major efforts have been made to develop comprehensive global and regional BA products, mainly according to two major programs: the ESA Fire Disturbance Climate Change Initiative (FireCCI) and the NASA MODIS Land Science Team. The current global BA products from the ESA FireCCI project include: FireCCI51 (2001–2020; 250 m), which derives from the MODIS surface reflectance imagery coupled with thermal anomaly observations [23]; FireCCILT11 (1982–2018; 5 km) from AVHRR LTDR [24], including new developments and sensors as in FireCCIS310 (2019; 300 m) from Sentinel-3 SYN coupled with VIIRS active fire hotspots [25]. Perspectives with newly delivered medium-resolution sensors have been tested regionally, such as FireCCISFD11 and FireCCISFD20 (2016 and 2019, respectively; 20 m) from the Sentinel-2 MSI coupled with the active fire data for sub-Saharan Africa [26,27]. On the other hand, the MCD64A1 collection 6.1 (2000–present; 500 m) is NASA’s current standard global BA product, which derives from MODIS daily surface reflectance imagery combined with MODIS active fire data [28].

Other available coarse-resolution products from different agencies include the Copernicus Climate Change Service Burned Area product, version 1.1 (C3SBA11) (2017–2022; 300 m) from Sentinel-3 OLCI data [29] and the European Forest Fire Information System (EFFIS) BA product (2000–present; 250 m and 20 m) from the MODIS and Sentinel-2 imagery [30]. Additional efforts have been made to provide finer-resolution BA products, a major end-user request [31]. Landsat-based BA mapping includes the novel 30-m resolution Global Annual Burned Area Maps (GABAM 1985–2019; 30 m), which derived from the Landsat dense time-series data by means of a global automated BA mapping approach [32] in the Google Earth Engine (GEE) [33]. In the same context, albeit with limited spatial coverage, other products include the Monitoring Trends in Burn Severity (MTBS 1984–2022; 30 m) across the whole of the U.S. [34], and the Landsat Collections 1 and 2 BA products for CONUS (1984–2022; 30 m) [35].

Although freely accessible for the scientific community with widespread use on different spatial and temporal scales, the above-mentioned BA products exhibit certain limitations. These are mainly caused by the inherent coarse spatial resolution that results in very high omission rates of small burned patches [26,27,36,37], particularly in the Mediterranean Basin where smaller fires happen [10,38], poor temporal fire reporting accuracy to prevent a fire seasonality analysis [32], and limited spatial and temporal coverage [34,35], which restrict their usage in other areas of the globe.

In Algeria, the available ground-based fire dataset provides invaluable information that can hardly be obtained by satellite-based systems. This includes, among others, the exact date and time of ignition/intervention/extinction, burned vegetation type and species, origin, and cause of fires. Nevertheless, this dataset is acknowledged to be incomplete, lacks fire perimeters and displays discrepancies in fire extent terms. This is attributed mainly to the visual estimation of fire-affected areas, which is often conservative especially in inaccessible areas. Furthermore, a standardized BA estimation protocol across local forest services in the country is lacking.

Considering these limitations in both national statistics and the performance of existing BA datasets, the development of a reliable and long-term BA product for such an insufficiently investigated part of the Mediterranean Basin would strengthen future research into forest management plans and for understanding Mediterranean fire hazards in this southern part of the Mediterranean Basin [5,13]. This would allow accurate in-depth analyses of the fire regime over lengthy time periods, and to learn the factors that underlie fire occurrence and propagation in this region with a Mediterranean climate, but with substantial socio-economic and political differences compared to the more studied European side of the Mediterranean Basin.

In recent years, several BA-mapping approaches have been developed for different study regions using medium-resolution data [39,40,41], including the BA Mapping Tools (BAMTs) [42]. By leveraging the powerful capabilities of the GEE’s cloud computing platform [33], the BAMTs constitute not only a significant stride as innovative, time-efficient, and resource-conserving tools for accurate multi-year BA mapping [10,43], but also the creation of independent reference data for validation exercises [44,45,46].

Based on these premises, we aim to exploit these efficient tools for systematically reconstructing the fire history in NE Algeria. Specifically, we aim to: (1) generate a BA product from the Landsat Collection-2 Surface Reflectance (LC2SR) product covering the 1984–2023 period; (2) assess its spatio-temporal consistency; (3) provide pieces of evidence for a significantly revised BA estimate compared to existing BA products (GABAM, FireCCI51, C3SBA11, MCD64A1, and EFFIS) and a ground-based fire dataset. This work constitutes the mandatory initial step for creating a spatially explicit BA database following international standards for the whole of Algeria to further open a major research field for fire hazard, impacts and vulnerability assessments that lead to firefighting and fire management policies [47].

2. Materials and Methods

2.1. Study Area

The study area is located in NE Algeria, on the southern rim of the Mediterranean Basin, and spans longitudes from 3.71° E to 8.68° E and latitudes from 36.21° N to 37.08° N. The total study area covers 17,036 km². It includes the coastal “Wilayas” (administrative provinces) of Annaba, Béjaïa, El-Tarf, Jijel, Skikda, and Tizi Ouzou in a narrow strip that stretches 440 km in length and up to 96 km in depth, with altitudes ranging from m.s.l. to 2 291 m (Figure 1a and b). Landscapes are typically Mediterranean dominated by forests, shrubs, and grasslands, where the mean annual rainfall varies from 725 to 924 mm (Figure 1c and Table 1). The main forest types are composed of oaks, such as cork-oak (Quercus suber L.), Algerian oak (Quercus canariensis Willd.), African oak (Quercus afares Pomel.) and holm oak (Quercus ilex L.), and of conifers, such as maritime pine (Pinus pinaster Aiton.), Aleppo pine (Pinus halepensis Miller.), and Atlas cedar (Cedrus atlantica (Endl.) Carrière.). Our study area is acknowledged as being one of the regional plant biodiversity hotspots in the Mediterranean Basin [48]. It harbors the densest, richest, and most pristine forest ecosystems, which are reflected in four National Parks (Djurdjura, Gouraya, Taza and El-Kala) that form part of the UNESCO World Network of Biosphere Reserves (WNBR) [49]. However, this region experiences the highest incidence of fires in the country [50,51] owing to high fuel availability, extreme climate conditions in summer (prolonged drought and Saharan heat waves), and increasing anthropogenic activities [5].

2.2. Input Data

2.2.1. Landsat and Sentinel-2 Imagery

The NASA/USGS Landsat program has launched eight Earth Observation (EO) satellites to date. These satellites have been continuously collecting data on the Earth’s land surface for over 50 years, the world’s longest EO data series (1972–present-day) [54]. Different multispectral and thermal sensors have been developed, including Landsat 1–5 MSS (1972–1999), Landsat-4 TM (1982–1993), Landsat-5 TM (1984–2012), Landsat-7 ETM⁺ (1999–2021), and the currently operational Landsat-8 OLI/TIRS (2013–present-day) and Landsat-9 OLI-2/TIRS-2 (2021–present-day) [55]. Landsat instruments from TM to OLI-2/TIRS-2 eras have been generating scenes over 185 km swath width using the Worldwide Reference System-2 (WRS-2) with seven, eight and 11 spectral bands at 15 m (panchromatic), 30 m (visible, NIR, SWIR, coastal/aerosol, and cirrus), and 60–120 m (thermal bands) spatial resolutions every 16 days, reduced to 8 days when two satellites are combined (currently, Landsat 8 and 9). The ESA/Copernicus Sentinel-2 Mission comprises a constellation of two identical EO satellites, Sentinel-2A and 2B, which have been simultaneously orbiting the Earth since July 2015 and July 2017, respectively. The Sentinel-2 MSI sensor captures scenes of 110 × 110 km² using the Military Grid Reference System (MGRS) with 13 spectral bands at 10 m (visible and NIR), 20 m (red edge and SWIR) and 60 m (atmospheric bands) every 10 days, or every 5 days with both satellites combined [56]. In this work, we used the LC2SR images from the last five Landsat satellites (1984–2023) and Sentinel-2 Level-1C Top of Atmosphere (TOA) reflectance images (2017 and 2021), both available as image collections in the Earth Engine Data Catalogue (EEDC) (https://developers.google.com/earth-engine/datasets, last accessed December 2023).

2.2.2. Global Burned Area Products

To evaluate our generated BA product, we made a comparison to different existing BA products. For this purpose, we used five readily available BA products with different spatial resolutions:

• The GABAM is the first and only medium-resolution BA product at 30-m resolution with a global scale and long-term BA data to date. A novel automatic pipeline [32] was applied to generate the annual global BA maps from the Landsat time-series on the GEE platform [33] from 1985 to 2019. Nevertheless, 1986, 1988, 1990, 1991, 1993, 1994, 1997 and 1999 are unavailable. Yearly GABAM composites were downloaded as 10° × 10° tiles in GeoTIFF format at ftp://124.16.184.141/GABAM (accessed February 2023).
• The FireCCI51 provides monthly global BA maps at the 250-m spatial resolution based on a hybrid algorithm coupling daily surface reflectance imagery and the active fire data from MODIS for 2001–2020 [23]. The FireCCI51 pixel product in GeoTIFF format was downloaded at https://doi.org/10.5285/58f00d8814064b79a0c49662ad3af537 (accessed April 2023).
• The C3SBA11 delivers monthly global BA maps at 300-m spatial resolution since 2017 onwards, using an adaptation of the FireCCI51 BA algorithm to Sentinel-3 OLCI data [25]. The pixel BA product was obtained from the C3S Climate Data Store (CDS) in NetCDF files at https://doi.org/10.24381/cds.f333cf85 (accessed April 2023).
• The MCD64A1 collection 6.1 BA mapping approach combines daily surface reflectance images with the active fire data from MODIS to produce monthly global BA maps at a spatial resolution of 500 m from 2000 onwards [28]. The product is provided by the USGS Land Processes Distributed Active Archive Centre (LP-DAAC) and was obtained in GeoTIFF format from the Application for Extracting and Exploring Analysis Ready Samples (AppEEARS) service (https://lpdaacsvc.cr.usgs.gov/appeears/, accessed April 2023).
• The EFFIS BA product provides daily updates of fire contours with information on the initial and final dates of burn detection, BAs, administrative units, and vegetation type for 49 countries in Europe, the Middle East, and North Africa (MENA) from 2000 to the present-day. The EFFIS Rapid Damage Assessment (RDA) module performs BA delineation by processing daily imagery from the MODIS instrument at a spatial resolution of 250 m [30]. Since 2018, EFFIS fire contours have been generated using the 20-m resolution Sentinel-2 imagery. For Algeria, EFFIS fire contours are available only for 2004 and 2005, and from 2009 onwards. The EFFIS BA product was provided in the ESRI Shapefile format from the EFFIS of the European Commission Joint Research Centre (https://effis.jrc.ec.europa.eu, accessed April 2023).

2.2.3. Active Fire Products

Active fire observations from the MODIS and VIIRS sensors were leveraged to examine the temporal fire reporting accuracy of our BA product by comparing detected burn dates to active fire dates. In line with this, the standard MCD14ML collection 6.1 product from MODIS onboard Terra and Aqua provides the systematic detection of active fires worldwide at 1-km resolution from November 2000 to the present-day [57]. In contrast, the 375 m standard VIIRS VNP14IMGML collection 1 product offers better accuracy in detecting active fires from January 2012 onwards [58]. These datasets were provided as feature point data in the ESRI Shapefile format by the Fire Information for Resource Management System (FIRMS) (https://firms.modaps.eosdis.nasa.gov/download/, accessed January 2024).

2.2.4. Ground-Based Fire Dataset

Ground-based annual BA estimates over the 1985–2023 period for the six wilayas under study were provided by the Directorate General of Forests (DGF), which operates under the authority of the Algerian Ministry of Agriculture and Rural Development (MADR), Algeria. The DGF oversees fire prevention and control, and develops national strategies and action plans for fire management across the country. Fire statistics are collected from 40 local forest services of the different wilayas that are regularly affected by fires [50]. However, detailed single-fire datasets have been recently compiled only for the 2012–2023 period, where detailed information is provided: smallest administrative division, forest name or locality, latitudes/longitudes, date and time of ignition/intervention/extinction, burned vegetation type, and more. Here we note that the minimum recorded fire size is not standardized across years and wilayas, a common feature that has been previously reported in Spain [59], France [60] or Tunisia [13], for which a standard threshold of 0.1 ha is usually assumed.

2.2.5. Land Cover Map

The ESA WorldCover project produced 10 m resolution global land cover maps for 2020 and 2021 from Sentinel-1 and Sentinel-2 data with 11 land cover classes conforming to the UN-FAO's Land Cover Classification System (LCCS) [52]. We obtained the GeoTIFF land cover map of 2021 for our study area from ESA at https://worldcover2021.esa.int/downloader (accessed May 2023).

2.3. Burned Area Generation

In 2014, Bastarrika et al. [61] developed a Landsat interactive BA mapping procedure, named the BA Mapping Software (BAMS). This software allowed to supervise the BA mapping between single pre- and postfire images by adopting a two-phase supervised strategy [62]. This approach has been widely employed to create fire reference perimeters to validate global and regional BA products [23,26,27,63,64,65], and to operationally reconstruct historical fires in Greece (1984–1991 and 1999–2011) [66], Argentina (1999–2011) [67] and Tunisia (1984–2010) [13]. However, it heavily depended on the system’s computational capacity and had several limitations, including the inability to work with higher resolution Sentinel-2 data and maintenance difficulties due to changes in the Landsat data format and metadata, which led to frequent malfunctioning and uncontrolled commissions over agricultural and cloud-shadowed areas [42]. To address these limitations, Roteta et al. [42] proposed the BA Mapping Tools (BAMTs). The current version (BAMTs v1.7) presents a set of six tools implemented in the GEE [33] that perform BA mapping (BA Cartography tool), the selection of validation areas (VA tool) and image dates (VA Dates tool), reference data creation (RP tool), validation (Assessment tool) and image data visualization (Image Viewer). The JavaScript API scripts for these tools with a user guide are publicly available at https://github.com/ekhiroteta/BAMT (accessed December 2023).

To generate our BA product, yearly pre- and postfire consecutive multitemporal composites were derived using 10,406 surface reflectance images from the LC2SR product from 1984 to 2023 with the BA Cartography tool to cover the whole study area. This tool is a semi-automatic BA algorithm that performs the supervised classification of Landsat or Sentinel-2 data based on a Random Forest (RF) machine-learning classifier [68]. For each year, the prefire period was defined from January 1^st to April 30^th, and the postfire period from May 1^st to December 31^st, following Roteta et al. [42]. This was done to employ the maximum number of cloud-free scenes and to, thus, reduce BA omissions, and to capture all the fires with detection delays caused by the Landsat revisit time (8–16 days).

In addition to six reflectance bands over the visible near-infrared (VNIR) and short-wavelength infrared (SWIR), for each Landsat image during the pre- and postfire periods, the BA Cartography tool computes three spectral indices: the normalized difference vegetation index (NDVI) [69]; the normalized burn ratio (NBR) [70]; the normalized burn ratio 2 (NBR2) [71]. The date with the lowest NBR value (date of the strongest burn signal) was used to generate image mosaics. Then a difference image was generated by subtracting the post- and prefire image composites. Subsequently, the RF classifier was iteratively trained using burned and unburned polygons, which were defined incrementally over both temporal composites until the desired BA delineation accuracy was achieved. On average, 38,894 burned and 45,723 unburned training samples (pixels) were manually and visually delineated every year. In the BA cartography tool, a shapefile vector layer was used as a mask to limit the processing extent to our study area. The RF model is widely applied for BA extractions from medium-resolution data [32,46]. Next, a two-phase supervised strategy to balance errors of commission and omission [61,62] was employed on the RF probability Images. In the first phase, seed pixels with strong burned signals were extracted using a threshold based on the average probabilities of the burned training samples. In the second phase, a region-growing process was applied up to a 50% threshold outwardly from the seed pixels to capture the entire burned patch using Rook’s case contiguity, where only the burned patches containing at least one seed pixel were labelled as burned.

Finally, each annual BA map was exported as four 2° × 2° tiles at the 30-m resolution in the ESRI Shapefile format, and a thorough fire-by-fire interactive visual quality check was performed in a GIS environment to ensure the highest quality of the final burned patches. This allowed a few anomalies to be eliminated, such as very few continuous strip lines (one to two pixels wide), often starting from the edge of detected burned patches. These errors are caused mainly by the Landsat-7 ETM⁺ SLC failure. Each burned polygon has a single attribute: indicating either the approximate burn detection date or null values for unobserved areas (i.e., clouds and cloud shadows).

At this step, it is important to note that the supervised BA contouring process under the BA Cartography tool may not be accomplished because of GEE’s scaling errors [42]. This issue was specifically encountered with years holding large amounts of BA. To avoid this shortcoming, we opted for a balanced and conservative approach when training the burned/unburned samples by avoiding very large polygons and complex geometries. The procedural steps for the generation and validation of the BA product (Phase I and II), and the intercomparison analysis (Phase III) are given in Figure 2.

2.4. Spatio-Temporal Validation

2.4.1. Spatial Validation

Validation plays a pivotal role in ensuring the integrity and reliability of satellite-derived BA products. The Committee on Earth Observation Satellites (CEOS) Working Group on Calibration and Validation (WGCV) defines validation as “the process of assessing, by independent means, the quality of the data products derived from those system outputs” [72]. The spatial validation of our BA product was performed using two validation years: 2017 and 2021. These years were selected because they hold the largest annual BA during Sentinel-2 era (2015–present). Sentinel-2 data are used in Section 2.4.3 below to create independent reference data. To assess the temporal fire reporting accuracy of the BA product, we compared the burn detection dates from the BA Cartography tool to the active fire dates from MODIS and VIIRS.

Five validation areas were sampled for 2017 and 2021 by the VA tool by following a stratified random sampling methodology [73,74]. At each sampling unit, this tool incorporates two key stratification criteria: the predominant Olson biomes [75] and fire activity according to FireCCI51 or MCD64A1. To ensure comprehensive and robust validation, long sampling units (LUs) lasting 100 days, and consisting of 110 × 110 km² tiles with a minimum temporal frequency of 30 days between consecutive Sentinel-2 images and cloud cover below 10%, were employed. This approach ensures a long data series with frequently available cloud-free images to create reliable reference data from the continuous and frequent ground observations throughout validation years. As the generation of high-quality reference data in validation areas of 110 × 110 km² (~ 12 000 km²) is time-intensive, the VA tool was designed to create smaller sampling units at the center of each Sentinel-2 tile [23,27,46,64,76]. However, this subsampling resulted in sampling units not being representative across the study area. To solve this problem, we manually generated a tessellated grid with 168 cells of 20 × 20 km² (400 km²) over the originally sampled validation areas [77]. Thereafter, we selected the cells with the greatest fire activity based on the BA estimates from FireCCI51 in the corresponding validation years. In all, 10 cells were selected as new validation sites, with two cells assigned to each sampled Sentinel-2 tile, and all of them lying in the Mediterranean forest biome (Figure 3).

To properly date the validation period, the CEOS Land Product Validation Subgroup (CEOS-LPVS) recommends generating reference data from multitemporal pairs of images rather than from single pairs of pre- and postfire images [42,76]. Therefore, 227 Sentinel-2 consecutive images with cloud cover below 10% were identified for each validation site using the VA Dates tool. Validation periods were about 5–7 months for 2017 and 4–8 months for 2021 (Table A1).

Satellite-derived BA product validation requires independent reference data at a higher spatial resolution than the product to be validated [78]. Accordingly, the RP tool was used to extract the 10-m resolution reference perimeters (RPs) from the selected pairs of consecutive Sentinel-2 images at each validation site. The RP tool performs the same BA mapping exercise as the BA Cartography tool, previously used in Section 2.3, albeit with some modifications. Instead of temporal composites, single-image pairs from the pre- and postfire periods were employed. Additionally, the SLC image or QA60 and B1 bands were used to mask clouds and cloud shadows. Likewise, the RF classifier was iteratively trained by defining burned and unburned training pixels until the optimal delineation of burned patches was achieved.

It is important to note that the “cloud cover” information used by the VA Date tool to select the cloud-free Sentinel-2 images at each validation site is that of the entire Sentinel-2 scene (110 x 110 km²) stored in the image metadata rather than being specific to the manually sampled validation sites (20 x 20 km²). Consequently, very few consecutive dates were excluded in the reference data creation step owing to high cloud cover that obstructed ground observations. The resulting Sentinel-2 reference data (hereinafter referred to as “S2RD”) were exported in the ESRI Shapefile format. A single attribute is given to the S2RD vector layers to indicate burned and unobserved polygons. Finally, a visual quality control of the short units of the RPs was performed before merging them into long units at each validation site [46,76,79] (Figure A1 and Figure A2).

At each validation site, the 2017 and 2021 BA maps were compared to their counterparts from S2RD. The commission error (CE), omission error (OE), overall accuracy (OA), Dice coefficient (DC) [80], and Relative bias (RelB) were computed employing the assessment tool based on the confusion matrix [81]. These accuracy metrics are widely used to assess the accuracy of satellite-derived BA products [45,79,82]. We considered additional metrics to measure the total area correctly detected as burned (SurfBA) or unburned (SurfUB). Similarly, SurfCE and SurfOE represent the total committed and omitted burned surfaces, respectively.

2.4.2. Temporal Validation

The burn detection dates in the BAMTs-derived product were determined based on the most frequently occurring date (the mode) for all the pixels in each detected burned patch in the Landsat post-fire composite [42]. Consequently, these dates do not imperatively correspond to the effective burn dates. To assess our BA product’s temporal reporting accuracy, the temporal delay in days was measured between the burn detection dates and the active fire dates from the MODIS and VIIRS hotspots [26,27,46]. Nevertheless, this analysis was performed by bearing in mind all the burned patches in the entire study area for a temporal validation period spanning from 2001 to 2023 rather than at the validation sites because taking specific years as the number of burned patches to be compared was insufficient. To accomplish this, 1 × 1 km² windows were created around the yearly active fire hotspots from MCD14ML from 2001 to 2011, which represents the approximate MODIS pixels. Similarly, windows of 375 × 375 m² were defined around the VNP14IMGML hotspots from 2012 to 2023. Next, we extracted the most frequent dates from the active fire windows that intersected the individual burned patches in the corresponding years then we compared them to the burn detection dates (Figure A3).

This analysis was performed on 33.68% of the burned patches (n = 22,236), which represented the intersection rate with the active fire windows from MODIS and VIIRS. Burned patches were found in 49.14% of the active fire windows (n = 9,030 from MCD14ML and n = 34,745 from VNP14IMGML) between 2001 and 2023. The remaining burned patches (66.32%) were not found in the active fire windows, which can be explained by active fire sensors’ native coarse spatial resolution or satellite overpass. In addition, environmental conditions like thick smoke plumes and dense vegetation cover significantly reduce fire detection. Ideally, the MODIS sensor routinely detects active fires of ~ 0.1 ha and, under ideal observational conditions, fires up to ~ 50 m² in size can be captured, but those of lesser extents or with low emitted fire radiative power (FRP) may not be detected [83]. Conversely, the VIIRS sensor has a minimum theoretical detectable fire size of ~ 20 m² at 1 200 K fires burning in the daytime and one half at 900 K instantaneous fires at night time [58].

2.5. Intercomparison Analysis

In this section, we made an intercomparison between our validated BA product (hereinafter referred to as NEALGEBA; NE Algeria BA) and five existing BA products with different spatial resolutions and input sensor data in terms of omissions and commissions. We also conducted validation and performed temporal trends using the ground-based BA estimates from the DGF.

2.5.1. Spatial Accuracy

To highlight the differences in terms of burned patch delineation between NEALGEBA and existing BA products, we compared the accuracy metrics (CE, OE, DC, RelB) obtained for all the BA products. To this end, we used the same Sentinel-2 reference data at the same previously selected validation sites to validate the NEALGEBA map for 2017 to assess GABAM, C3SBA11, FireCCI51, MCD64A1 and EFFIS because this is the common year between all these BA products [45,46]. The 2017 BA maps from FireCCI51 and MCD64A1 were already available as image collections in the GEE’s data catalogue, and were directly evaluated against S2RD in the assessment tool. The GABAM, C3SBA11 and EFFIS maps were firstly adjusted to the required format and uploaded to the assessment tool. It is important to note that GABAM does not provide the burn date, as required by the assessment tool, but indicates only the burn year. To overcome this situation, we assigned approximate burn detection dates to the burned patches on the 2017 GABAM map by overlaying the BA map from NEALGEBA. The burned patches with the largest overlap were assigned the same burn detection date, whereas those with no overlap were not dated.

2.5.2. Cross-Correlation with Ground-Based Fire Dataset

To assess NEALGEBA performance in BA estimation terms, we performed a cross-correlation analysis among the total annual BA estimates from our product, existing BA products and the ground-based fire dataset from the DGF over the 13 available years during the 2004–2019 period using the “PerformanceAnalytics” R-package [84]. However, the C3SBA11 product was not considered since only three years were available (2017, 2018, and 2019) among the overlapping years between the BA products. Here and in the subsequent step (Section 2.5.3), exclusively, the land cover class 10 (tree cover), class 20 (shrubland), and class 30 (grassland) from the ESA WorldCover map 2021 [52] were assumed to represent the spatial extent of natural vegetation landscapes and were used to mask out detected burned patches outside these classes from NEALGEBA and existing BA products. This was done following the definition of wildland fires in Algeria [50] based on which the fraction of cropland fires is not taken into consideration when reporting BA estimates by the local forest services of the DGF.

2.5.3. Temporal Burned Area Trends

We analyzed the temporal trends of the annual BA estimates from filtered NEALGEBA and the DGF datasets for the 1985–2023 period (the DGF dataset lacks data for 1984) using the Mann-Kendall trend test [85,86] and Sen’s slope estimator [87]. These two non-parametric methods are widely used to test the significance of trends in fire regimes [88].

3. Results

3.1. Analysis of the Generated BA Product

3.1.1. Spatial and Temporal Patterns

The BA mapping phase resulted in the generation of spatially explicit annual BA maps covering all fires occurring in the natural landscapes and rural areas of the six wilayas in the NE coast of Algeria at the 30-m resolution for the 1984–2023 period (Figure 4a). In the last four decades, the total mapped BA was 1.16 M ha with an average of 29,130.5 ha year⁻¹ affecting forests (65.2%), shrublands (9.30%), grassland (22.2%), croplands (2.79%) and other types (< 0.6%). This BA is double the official BA estimate from the DGF (0.65 M ha). The total annual BA showed wide temporal variability in the last 40 years, with 68.8% of the total BA occurring in 13 peak years above the mean. The most extreme years were 1994, 2000, 2017 and 2021, when more than 15.0%, 6.6%, 5.3% and 8.4% of the total BA were respectively observed. These years accounted for up to 35.3% of the total BA in the entire time series. Conversely, the lowest annual BAs were found in 1984, 1996 and 2018, which represent less than 1% of the total BA over the entire study period (Figure 4b). The study area showed a high fire recurrence rate with more than 52.7% of the total BA affected by two fire events or more in the last four decades (Figure 4c). Indeed, the total mapped BA (1.16 M ha) was twice the extent of the land affected by at least one fire (0.58 M ha). Only 47.3% of this land was burned once, while 25.9% was burned twice and 25.8% was affected by three fires or more (Figure 4c).

Besides temporal variability, the annual BA showed wide spatial variability between the different wilayas. The fraction of burned burnable land represented 101.2% (2.53% year⁻¹) in El-Tarf, 94.5% (2.36% year⁻¹) in Béjaïa, 81% (2.025% year⁻¹) in Tizi Ouzou, 78.6% (1.96% year⁻¹) in Annaba, 67.3% (1.68% year⁻¹) in Skikda and 61.8% (1.54% year⁻¹) in Jijel. When looking at the fire event spatial pattern independently of their surface, the kernel density estimation (KDE) analysis showed fire hotspots in the Tizi Ouzou and Béjaïa wilayas, which are high ignition rates in these wilayas (Figure 5a). However, when fires were weighted by their sizes, the WKDE showed more fire hotspots in the Skikda, Annaba, and El-Tarf wilayas (Figure 5b).

3.1.2. Fire-Size Distribution

The minimum detectable fire size equaled ~ 0.1 ha, which corresponded to the areal unit covered by a Landsat pixel on the ground. The maximum fire size was 26,300 ha, which occurred in the Tizi Ouzou wilaya in 2021. Fire size distribution using all the fire patches (n = 98,259) followed the self-organized criticality hypothesis [89] up to a fire size threshold of 1.5 ha, below which some small fires could be missing (Figure 6).

3.1.3. Fire Seasonality

Fire seasonality, based on fire number and BA, could be generated on 22,115 fire patches, which could be dated out of the 66,018 fire patches identified over the 2001–2023 period when the VIIRS/MODIS hotspots were available. NE Algeria displays typical Mediterranean fire activity, and July and August are the most burnable months, which prolong to September and October (Figure 7).

3.2. Spatio-Temporal Validation

3.2.1. Spatial Validation

The results of the accuracy assessment of the NEALGEBA maps of 2017 and 2021 at the 10 validation sites using S2RD are presented in Table 2. The corresponding validation images are found in Figure B1 and Figure B2. For the validation year 2017, the OA was 98.22%, with a CE of 7.96%, an OE of 8.19% and a DC of 91.92%. Notably, validation site 31SEA-E3 had the highest OE with 30.16%, which resulted in an underestimation of the total BA of 23.19% and the lowest DC (79%) of all the validation sites. The highest commissions (CE = 9.65%) occurred at 32SMF-Y3. The overall RelB indicated a slight underestimation of the total BA of 0.24%. The total area correctly detected as burned was 195.51 km², while 1,700 km² of unburned surface was correctly detected at all the validation sites (covering 2,000 km²). Additionally, the product misclassified 16.92 km² of the unburned surface as burned and omitted 17.43 km² of the burned surfaces. For 2021, the product achieved similar results to those for 2017 in terms of OA (98.15%) and commissions (7.92%), albeit with a twofold decrease in omissions (4.76%) that resulted in a higher DC (93.63%). Of all the validation sites, the highest commissions and omissions were for 32SKF-P3 with values of 25.38% and 14.88%, respectively, and an RelB value of 14.08%. The DC was the lowest at all the validation sites with 79.52%. The total area correctly detected as burned was 270.70 km² and 1,689.84 km² was correctly detected as unburned. There were instances where the product incorrectly classified 25.53 km² of unburned surface as burned and omitted 13.52 km² of the burned surface. Overall, cloud coverage was the primary source of BA omissions, whereas most commissions were caused by the BA algorithm over land covers spectrally akin to the burned surfaces, such as harvested croplands, labored lands, water bodies and shadowed areas.

3.2.2. Temporal Validation

The NEALGEBA product showed significant fire detection delay compared to the active fire hotspots from MODIS and VIIRS. On average, burned pixels were detected 42.55 days after hotspot detection: 21.07% on the first 15 days, 46.47% in the first month and 75.67% in a 2-month window (Figure 8). This disparity was attributed mainly to Landsat’s low temporal revisit frequency (8–16 days), which could be worsened by the unavailability of cloud-free images for BA extraction purposes.

3.3. Intercomparison Analysis

3.3.1. Spatial Accuracy

The evaluation of the NEALGEBA product versus existing BA products showed our product’s overall outperformance (Figure 9). Indeed, NEALGEBA had the lowest CE (7.96%), the highest DC (91.92%) and the best total BA estimation compared to S2RD (RelB = -0.24%). Despite being generated using the same input data (30 m Landsat images) and having a lower OE (3.89%), GABAM obtained a higher CE (13.33%) than NEALGEBA (7.96%), which resulted in an overestimated BA (RelB = 10.89%) compared to our product (RelB = -0.24%). As expected, coarser resolution products based on MODIS and Sentinel-3 OLCI: FireCCI51, MCD64A1, EFFIS and C3SBA11, produced the highest CEs and OEs, the lowest DCs with a significant lack of estimation of the total BA compared to S2RD. Specifically, FireCCI51, C3SBA11 and MCD64A1 had CEs ≥ 37.43%, which resulted in a marked overestimation of the total BA that reached 58.39%, 43.10 %, and 36.03%, respectively. However, these three products had relatively lower omissions ≤ 14.89%, compared to the EFFIS product (≥ 33.26%). The last one had the highest OE of all the assessed BA products, which led to the highest BA underestimation rate (RelB = -12.69%). Detailed results of the accuracy assessment of the GABAM, FireCCI51, C3SBA11, MCD64A1 and EFFIS BA maps of 2017, along with their corresponding validation images, appear in Tables B1 to B5, and in Figures B3 to B7, respectively.

The visual analysis clearly shows that all the detected large burned patches were consistent in all the BA products, although most relatively small fires were captured only by medium spatial resolution products (NEALGEBA and GABAM) (Figure 10a). GABAM demonstrated good spatial matching with S2RD, but obtained a high CE (13.33%) at all the validation sites (Table B1). This resulted in relatively extended burned patch limits compared to NEALGEBA. In other words, GABAM tended to overestimate the BA and included additional areas that may not have been burned. NEALGEBA showed higher OEs at all the validation sites, which resulted in slightly reduced burned patches compared to S2RD (Table 2). This indicates a tendency to exclude some burned patches, which leads to reduced BA compared to GABAM. These products tend to slightly underestimate or overestimate burned areas, respectively, compared to the higher resolution reference data from Sentinel-2. Furthermore, both omitted a few small burned patches and misclassified a few unburned pixels as burned. It is important to consider the trade-off between the CEs and OEs when evaluating the performance of BA maps.

Conversely, the larger pixel size in FireCCI51 (250 m), C3SBA11 (300 m), and MCD64A1 (500 m) led to size of burned patches being significantly overestimated. These products tended to overlook small fires and were unable to accurately delineate larger ones. They also displayed difficulty in discerning neighboring fire patches by perceiving them as single fire events, and they omitted the unburned islands in the main fire patches, which meant that the total BA was perpetually overestimated. EFFIS displayed a smoothing effect that oversimplified fire perimeters, with loss of details inside burned patches (unburned islands), which led to significant CEs. Particularly in this product, no fires were captured at validation site 31SEA-E3, an important part of the BA extent was omitted at validation site 32SKF-V2, and relatively small and spatially fragmented fires, at a distance from large burned patches, were completely omitted at the remaining validation sites. This resulted in significant OEs of the BA. As Figure 10b depicts, the total detected amount of BA at the five validation sites in S2RD of 2017 was 212.9 km². This was slightly reduced in NEALGEBA with 204.0 km² and slightly higher in GABAM with 231.4 km². Due to high commissions in FireCCI51, C3SBA11 and MCD64A1, the amount of BA was overestimated with 346.1 km², 316.1 km², and 296.9 km², respectively, whereas EFFIS significantly detected less BA with only 188.8 km² due to high omission values.

3.3.2. Cross-Correlation with Ground-Based Fire Dataset

The evaluation of the different BA products compared to the DGF dataset was done using the log10-transformed total annual BA estimates over the 13 overlapping years that comprised the 2004–2019 period between the different products. In general, the distinct BA products showed a strong correlation with the DGF data over the entire study area (Figure 11a). NEALGEBA, FireCCI51, and MCD64A1 obtained the highest correlations (r = 0.96), while GABAM yielded a slightly lower value (r = 0.94) even though it had the same spatial resolution as NEALGEBA (30 m). EFFIS had the lowest correlation with the DGF dataset. On the scale of the different wilayas, NEALGEBA demonstrated a better correspondence with the DGF estimates than existing BA products (Figure 11b).

3.3.3. Temporal Trends of the Burned Area

Generally, the temporal evolution of total annual BA estimates from NEALGEBA and DGF displayed a consistent increasing trend over the 1985–2023 period (1984 is unavailable in the DGF dataset), yet differ in the annual rate, with NEALGEBA showing significant trends. For all the wilayas combined, NEALGEBA showed a significant increase rate of 532.4 ha year⁻¹ at 10% (p-value = 0.053) in the total annual BA while DGF showed lower increase rate (197.1 ha year⁻¹; Figure 12a). At the wilaya scale, the annual BA from the two products showed significant positive trends at 5% (p-value ≤ 0.05) for the two products in Béjaïa and only significant for NEALGEBA in Jijel (p-value = 0.05) and in El-Tarf at 10% (p-value = 0.056). In general, NEALGEBA showed higher increased BA rates (182.2, 105.6 and 67.9 ha year⁻¹ in Béjaïa, Jijel, Skikda and El-Tarf, respectively; Figure 12b).

4. Discussion

In this analysis, we reconstructed and validated 40 years (1984–2023) of historical fire events at fine spatial resolution in typical Mediterranean ecosystems of NE Algeria. The newly generated NEALGEBA product represents the first and most extensive time series of BA in this part of the Mediterranean Basin, which faces a substantial fire occurrence threat. The BA product generation (phase I) proved the high potential and reliability of the BA Cartography tool in generating spatially consistent annual BA maps based on a Random Forest supervised classification and a two-phased strategy [42,62]. Despite being labor-intensive and heavily relying on the visual interpretation of pre- and postfire temporal composites, this semi-automatic procedure enabled high control over CEs and OEs and, thus, improved the BA product’s quality. The analysts’ expertise is more involved in selecting representative burned (seeds) and unburned samples with an iterative analysis of BA delineation [77], a considerable asset that is not provided with fully automated methods [32]. Additionally, the visual quality control and manual refinement of the generated fire perimeters allows to reduce potential anomalies such as those caused by the sensor. The RP, VA, VA dates, and the assessment tools, greatly facilitated the spatial validation exercise of satellite-derived BA products compared to previous studies [64,90], and all in accordance with the BA assessment standardized protocol endorsed by the CEOS. These tools ensured the creation of high-quality reference data (RP tool) from consecutive 10-m cloud-free Sentinel-2 images (VA Date) located at the validation sites preselected by stratified random sampling (VA tool). The assessment tool allowed a full-automatic comparison of the BA maps to the Sentinel-2 reference data, and reported accuracy metrics (CE, OE, DC and RelB) at each validation site.

The accuracy assessment of the 2017 and 2021 NEALGEBA maps showed remarkable results with CEs of 7.96% and 7.92%, OEs of 8.19% and 4.76% and DCs of 98.22% and 98.15%, respectively. These metrics are consistent with better performance than those obtained in the original case study in south-eastern Australia, in which a BA product for the 2019 / 2020 fire season was generated and validated using the same input data with a CE of 11.80%, an OE of 8.90% and a DC of 89.60% [42]. However, the larger pixel size (30 m) in NEALGEBA led to a subtle alteration in the extent of the burned patches, which meant that their boundaries slightly extended outwardly compared to the reference perimeters from the 10-m Sentinel-2 independent reference data. We also observed that almost all the spatially isolated small burned patches were misclassified as burned, and mainly in 2021. This was due to the algorithm’s limitation to discriminate small spectrally confusing surfaces with a similar spectral response to the burned surfaces. For the 2021 fire season, most of the BAs were in mountainous areas, which made it quite challenging to select representative and sufficient burned seed pixels to capture the entire burned patches, and thus, to reduce omissions. We attempted to avoid burned pixels in shadowed areas to reduce commissions on the classification map. We found that the algorithm failed to ensure the continuity of some large burned patches, and omissions occurred mostly on the edges of the main burned patches and on unburned islands with very few small isolated burned patches that were completely omitted. Overall, the obtained accuracy metrics indicated the NEALGEBA product’s spatial consistency. However, the temporal validation using the active fire hotspots from MODIS and VIIRS underlined its limitation for accurately reporting fire events over time. This was explained by the long revisit time of the Landsat satellites (8–16 days), atmospheric conditions (i.e., clouds) and temporal compositing, which could significantly delay BA detection. This is not uncommon in medium-resolution products [42,46] compared to MODIS-derived products that incorporate active fire information [23,28], and underscores the need for further development to improve temporal uncertainty. Employing data from satellite sensors with a higher observation frequency, such as Sentinel-2 and the active fire hotspots from VIIRS, would reduce temporal reporting delays [26].

Coarse resolution BA products were found to significantly overestimate the total BA on a finer spatial scale due to a larger pixel size (≥ 250 m), unlike the continental scale on which the total BA was overly underestimated when compared to more accurate data from Sentinel-2 MSI sensor [26,27,36,37]. In addition, their limited temporal coverage (2001–present) prevents long-term fire studies compared to the BA products generated from the Landsat data archive dating back to 1984. The validation of the 2017 EFFIS BA map showed that the latter presented the highest omissions of all the assessed BA products, which resulted in a considerable underestimation of the total BA, plus a smoothing effect on fire perimeters that roughly delineated the burned patches. These inconsistencies have been previously reported in [40,91,92], and are attributed mainly to the 250-m coarse-resolution input data from MODIS used to generate the EFFIS product. On the other hand, GABAM is, to date, the only available global high-resolution BA product to provide BA mapping at finer spatial resolution to reliably detect smaller burned patches [32]. However, this product was generated in yearly composites by providing only the burn year without indicating the approximate burn detection date, which prevents a fire seasonal analysis. Moreover, while GABAM has the longest time span amongst global BA products (1985–2019), some years remain unavailable. In addition to the reported commissions over agricultural lands, significant systematic errors were observed when we examined the GABAM annual maps versus the corresponding Landsat post-fire image composites in a Long SWIR/NIR/Red color composition and NEALGEBA. The former represents BA commissions over water bodies, unburned forest areas, clouds, flares in oil/gas facilities and Landsat strip errors (Figure B8 a to c). Additionally, significant errors occurred in 2002 (Figure B8 d) and were perhaps caused by the significant alteration of the primary functioning mode of Landsat-5 TM’s scan mirror, known as the scan angle monitor (SAM), which led to internal synchronization problems. This failure caused diagonal patches of anomalous observations with very high reflectance values in the long shortwave infrared (SWIR2) towards the scene footprint edges, and led to false fire detections. The SAM system was then switched to an alternative one called the bumper mode [93], which overcame this anomaly. The semi-automatic BA extraction procedure, which uses the BA Cartography tool along with a thorough visual inspection of the mapped burned patches, allowed these anomalies to be mitigated and consistent results to be obtained on the NEALGEBA annual maps. Overall, these limitations highlight the challenges and complexities involved in using existing BA products to accurately characterize local and regional fire regimes. The newly generated BA product herein presented serves as a surrogate to existing BA products by offering improved spatio-temporal resolutions allowing for a thorough assessment of fire impacts on forest ecosystems and, in turn, assists in designing strategies and adapted action plans to mitigate their severity in NE Algeria. Additionally, BAMTs can be easily leveraged by forests and protected areas managers in Algeria to operationally extract burned perimeters and to integrate complementary field data, especially the day and time of ignition, which reduces temporal uncertainties.

Our first evaluation of the newly generated NEALGEBA could properly address the fire seasonal distribution that spans from July to October and peaks in August. This result is in accordance with the seasonal drought and fire weather index seasonality characterizing our region [5]. Similar fire season length is reported in neighboring Tunisia [13] and in Portugal [94], and a slightly longer than the usual July-September fire season reported in Morocco [95], Italy [96,97], Greece [98], Bulgaria [99], and the Iberian Peninsula [100]. The length of the fire season in our region may be shaped by the early and late fire-prone weather conditions favored by climate change effects and socio-economic factors [5,9]. Regarding affected vegetation, we obtained a fraction of burnable areas affected by fires reaching 2.95% year⁻¹ in fair agreement with Portugal (3.31%), making our BA estimates on the highest range of variability observed in the Mediterranean Basin. Only 0.19% was observed from Tunisia [13], Lebanon (0.58%), France (0.53%), Greece (0.57%), Spain (0.84%) and Italy (1.14%), as reviewed in [10]. We also detected an increasing trend in BA (532.4 ha year⁻¹) over the region. This trend is different from that in the northern part of the Mediterranean, where a general decreasing fire trend has been observed [101]. More precisely, an abrupt decrease has been observed in France in 1990 with increasing firefighting expenditures [60], and an increase in the 1980s and 1990s then a decrease since 2000 in Spain owing to land abandonment [59,102], and a decreasing trend was observed in Greece [14]. In the middle east and north Africa, no trends have been particularly observed since the 1980s [10,13], consistent with our study. However, the recent collapse of political regimes that have led to an abrupt increase in the BA in Tunisia [103], that we did not detect in Algeria not affected by this political collapse. Throughout history, Algeria has although been marked by significant political events that have led to heavy burning period, such as the Algerian Civil War (the Black Decade) in the 1990s [6], which was reflected by the highest peak in BA observed in 1994. Algeria has remained quite stable since 2011, when the Arab Spring started in the southern Mediterranean Basin and did not, thus, affect the BA during this recent period. Here we note the exceptional heat wave that hit the region in 2023 with record-breaking temperatures in April [104] and in July 2023, and with some casualties during fire events, which were widely reported in the media, but did not lead to the most extreme fire year compared to 1994 and 2021. Hence, NEALGEBA appears as a keystone database to provide accurate information on burned area and its temporal trends, and as a reference database to allow for the objective contextualization of fire years, to be further used in fire weather analyses, fire impacts assessments, fire model benchmarking or euro-Mediterranean initiatives of fire-related issues.

In its current version, NEALGEBA covers all types of fires that have occurred across all landscapes in NE Algeria from 1984 through 2023. However, this product, as the case with all satellite-derived BA products, exhibits specific limitations that necessitate reporting for future improvements. First, the burn dates indicated in our product do not match the effective fire dates, when the fires were actively burning. In fact, the BA Cartography tool computes the modal date from all pixels within each detected burned patch in the yearly Landsat post-fire composite. This results in significant fire detection delay, which impacts the product’s temporal fire reporting accuracy, especially in large fires lasting several days. We highlight here that using data at a higher temporal resolution from Sentinel-2 MSI could significantly reduce this disparity [26,27,43]. Second, the product’s commission errors were primarily observed over agricultural lands (harvested or ploughed croplands) which exhibit similar spectral features as to that of burned areas, characterized by abrupt changes in the reflectance data particularly in the near and shortwave infrared bands [27,32,35,61,62,105]. Here, we should also emphasize the high uncertainties associated with the detection of cropland fires [106]. Third, very large fire events enduring several days may not be effectively captured as single burned patches due to the low revisit frequency of Landsat satellites and availability of cloud-free image. Some large fire patches may not be spatially contiguous owing to low burn signal over shadowed areas, sparse vegetation, or discontinuity in vegetation cover. Additionally, spotting fires can result in spatially isolated burned islands from the main fire patches. Fourth, Landsat sensor anomalies were one of the main challenges especially the Landsat-7 ETM⁺ SLC failure, which affected most of its time coverage. Post-processing procedures have been applied to mitigate commissions caused by these anomalies. However, some omissions or late fire detections (strips within an actual fire patch) should be acknowledged. Regarding future work, it is envisaged to complement the NEALGEBA product with detailed information contained in the ground-based fire inventories from local forest services of the DGF, especially for extreme fire events. For instance, burn detection dates can be corrected, thus reducing the product’s temporal uncertainties. Possible attributes include: forest name or locality, date and exact time of ignition/intervention/extinction, burned vegetation type and species, land ownership, cause of ignition, perpetrator of the fire, fire reporter, weather conditions, participating bodies in fire suppression, damage assessment, and investigation. Ongoing efforts involve expanding NEALGEBA to a country-level BA product with continuous mapping of fire affected areas for the upcoming years using higher resolution imagery data from Sentinel-2 MSI. This aims to provide accurate characterization of the spatio-temporal patterns of fires across a larger geographical scale.

5. Conclusions

This paper presents a new spatially explicit fire database (NEALGEBA) in Mediterranean-type ecosystems of NE Algeria. This database is a reconstruction of 40 years (1984–2023) of historical fire patterns from Landsat data. NEALGEBA product offers the Algerian forest managers and policymakers valuable insights into the spatio-temporal distribution and magnitude of fires in this region which is prevalently and substantially vulnerable to extreme and large fires. Its provision of accurate and long-term BA mapping surpasses the limitations of existing global and regional BA product and is, thus, adequate for disentangling fire regime and identifying the key driving factors of fire occurrence and propagation in this part of the Mediterranean Basin. Furthermore, this dataset significantly contributes to further national and international fire hazard and impact assessments, and acts as a reference for contextualizing future weather extremes, such as the 2023 exceptional heat wave, which we show not to have led to the most extreme fire year over the last 40 years.