Spatiotemporal Analysis of Fire Foci and Environmental Degradation in the Biomes of Northeastern Brazil

Forest fires destroy productive land throughout the world. In Brazil, mainly the Northeast of Brazil (NEB) is strongly affected by forest fires and bush fires. Similarly, there is no adequate study of long-term data from ground and satellite-based estimation of fire foci in NEB. The objectives of this study are: (i) to evaluate the spatiotemporal estimation of fires in NEB biomes via environmental satellites during the long term over 1998–2018, and (ii) to characterize the environmental degradation in the NEB biomes via orbital products during 1998–2018, obtained from the Burn Database (BDQueimadas) for 1794 municipalities. The spatiotemporal variation is estimated statistically (descriptive, exploratory and multivariate statistics) from the Normalized Difference Sustainability 2022, 14, 6935. https://doi.org/10.3390/su14116935 https://www.mdpi.com/journal/sustainability Sustainability 2022, 14, 6935 2 of 19 Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and Standardized Precipitation Index (SPI) through the Climate Hazards Group InfraRed Precipitation Station (CHIRPS). Moreover, we identify 10 homogeneous groups of fire foci (G1–G10) with a total variance of 76.5%. The G1 group is the most extended group, along with the G2 group, the exception being the G3 group. Similarly, the G4–G10 groups have a high percentage of hotspots, with more values in the municipality of Grajaú, which belongs to the agricultural consortium. The gradient of fire foci from the coast to the interior of the NEB is directly associated with land use/land cover (LULC) changes, where the sparse vegetation category and areas without vegetation are mainly involved. The Caatinga and Cerrado biomes lose vegetation, unlike the Amazon and Atlantic Forest biomes. The fires detected in the Cerrado and Atlantic Forest biomes are the result of agricultural consortia. Additionally, the two periods 2003–2006 and 2013–2018 show periods of severe and prolonged drought due to the


Introduction
Forest fires are a global natural disaster with a greater impact on forest ecosystem carbon cycling, and they are due to global warming, the annual increase in the amount of combustible materials and difficulties in controlling fire sources [1][2][3]. The climate and weather conditions are the dominant drivers of forest fires. Wu et al. [4] presented that global climate change can affect the frequency, intensity and the length of the wildfire season. Similarly, Li et al. [5] evaluated the climate change impact on future wildfire risk on the basis of Global Climate Models (GCMs) and suggested that climate change will increase the wildfire risk across the United States, particularly in the South Central region. However, the authors reported considerable variations across the future climate scenarios projected by the GCMs under the different Representative Concentration Pathway (RCP) scenarios proposed by Intergovernmental Panel on Climate Change (IPCC) Reports. The meteorological anomalies of surface temperature, surface pressure, precipitation, cloud water content, surface horizontal winds and 500 hPa vertical velocity were the main factors related to the devastating wildfires in California that occurred in October 2017 [5].
The wildfire risk and forest fire spread are also influenced by several other factors, such as topography, changes in LULC, distance from the road, source of ignition and moisture content of the fuel [6][7][8][9][10][11][12][13]. Moreover, the fire spread along the direction of the slope, elevation and topographic characteristics of the surface area are also important. For example, fires spread more quickly up-slope than down-slope. The up-slope fires accelerate the combustion process by the heat rise in front of the fire and dry the fuel faster up in the mountain more effectively [10][11][12][13]. Gheshlaghi [14] developed a fire risk map from slope, aspect, altitude, land cover, NDVI, yearly rainfall, temperature, distance to settlements and distance to roads in a Geographic Information System (GIS) based on analytical techniques. Furthermore, Salavati et al. [13] assessed fire risk potential in the city of Sanandaj in Western Iran using a statistical approach. They considered the local wildfire risk, including altitude, slope percentage, slope direction, distance from the road, distance from the river, LULC, average annual rainfall and average annual temperature as main factors. The above findings can assist decision-makers in taking appropriate precautions to prevent forest and rangeland fires and/or to minimize the fire damages.
Several approaches have been used to understand the forest fire propagation and wildfire risk to improve the fire control [15][16][17]. Cellular Automata (CA) have been used for a long time to model the complicated mechanism of fire spread [18]. Such simulations are driven by the Rate of Fire Spread (ROS), which is difficult to assess, hindering its implementation in the forest fire management system so far. Li et al. [18] also proposed two modified ROS models in order to improve the risk management and effectively implement fire suppression policies. They used a normal Long Short-Term Memory (LSTM)-based speed model with a one-dimensional matrix consisting of wind speed and ROS as inputs.
As an important conclusion from the study, the authors indicated that the process of forest fire spread is very complex, where the ROS was influenced by many factors, and indicated the need for more studies on this issue.
Forest fires annually destroy millions of hectares of productive land worldwide [19,20]. Brazil is among the countries that are the most affected by forest fires and burned cases [21,22]. Forest fires differ from those occurring in Brazil due to the extent of their effects, changes in LULC and use for deforestation [23,24]. It is worth noting that the use of fire is a cultural practice in Brazil, being applied in the management of agricultural crops, in the creation and maintenance of pastures for livestock and for cleaning urban waste [21,25]. Today, the impact of wildfires on the Earth's climate is a global concern due to the global emissions of water vapor (H 2 O), carbon dioxide (CO 2 ), methane (CH 4 ), nitrogen dioxide (N 2 O) and aerosols that are changing the Earth's radiative balance [26]. The aerosols reduce the transmission of solar energy to the Earth's surface, while greenhouse gases (GHGs) trap solar radiation [26,27]. Forest fire is a primary process to influence the composition and structure of the vegetation of the affected site and further exert a direct impact on biogeochemical cycles and the landscape [28,29].
In Brazil, forest fires are a concern among society and a governance problem, mainly due to the emission of GHGs that contribute to global warming [23,30]. GHG emissions via fires contribute to increased air pollution and worsening air quality in several metropolitan regions (MR) of Brazil [31,32]. Fires also contribute to the reduction of water in forests, affect human health and cause the intensification of extreme rainfall via increased condensation nuclei [33,34]. Since 1997, Brazil has been monitoring fire foci on the basis of the calculation and forecasting of the fire risk using environmental satellites from the Burn Database (BDQueimadas) developed by the Weather Forecast and Climate Studies Center (CPTEC)/Institute National Space Research Institute (INPE) [35]. BDQueimadas helps to prevent and minimize the environmental impacts caused by fires, especially in Environmental Preservation Areas (EPA), Conservation Units (CU) and Legal Reserves (LR) in the country [22,23,25,36]. Facing a growing number of fires and burnings, Brazil created the National Center for the Prevention and Fighting of Forest Fires (PREVFOGO) [37].
Previously, the spatiotemporal analysis of the fire foci in the NEB was evaluated via multivariate analysis (cluster analysis-CA) using the Ward method [21]. NEB was the only homogeneous group of fire foci in Brazil formed during 1998-2011, and this group showed a strong influence on multi-scale meteorological systems and economic activities, including livestock, agriculture and agropastoral activities [21,22]. However, the dynamics of the environmental degradation were not addressed for the fire foci in the biomes existing in the NEB (Atlantic Forest, Caatinga, Cerrado and Amazon). Therefore, characterizing the variability of fire foci in the NEB from space and ground-based products is fundamental for the implementation of new public policies and environmental management. Therefore, the objectives of this study are to first evaluate the space-temporality of the fire foci in the biomes within the NEB via environmental satellites, and secondly to characterize the environmental degradation in the NEB biomes via satellite data. Section 2.2 of this paper contains detailed information on fire foci data obtained for the years 1998-2018 via BDQueimadas. Section 2.3 describes the detailed formulation of the multivariate analyses for the distribution and counting of fire foci divided into homogeneous groups. Section 2.4 shows the products obtained via remote sensing regarding precipitation and vegetation indices in order to evaluate the behavior of vegetation and associate it with fire foci, in search of confirmation of the relationship between LULC and fire foci. Similarly, Section 3 includes the results and information related to the distribution of groups of fire foci in a temporal, annual, monthly and spatial manner. Furthermore, the spatiotemporal distribution is estimated from the NDVI, Enhanced Vegetation Index (EVI), Standardized Precipitation Index (SPI) and precipitation. Section 4 presents a precise discussion on all the results and Section 5 refers to the conclusions of the research in this study.

Study Area
The NEB region has an area of approximately 1,558,196 km 2 and a population of 53,081,510 inhabitants [38]. The NEB consists of nine states, namely: Maranhão (MA), Piauí (PI), Ceará (CE), Rio Grande do Norte (RN), Paraíba (PB), Pernambuco (PE), Alagoas (AL), Sergipe (SE) and Bahia (BA) (Figure 1). The NEB has high inter-annual and intra-annual rainfall variability, with the thermal influence of the Pacific and tropical Atlantic oceans, which control and modulate years of drought and floods in the region [39][40][41].
fire foci. Similarly, Section 3 includes the results and information related to the distribution of groups of fire foci in a temporal, annual, monthly and spatial manner. Furthermore, the spatiotemporal distribution is estimated from the NDVI, Enhanced Vegetation Index (EVI), Standardized Precipitation Index (SPI) and precipitation. Section 4 presents a precise discussion on all the results and Section 5 refers to the conclusions of the research in this study.

Study Area
The NEB region has an area of approximately 1,558,196 km 2 and a population of 53,081,510 inhabitants [38]. The NEB consists of nine states, namely: Maranhão (MA), Piauí (PI), Ceará (CE), Rio Grande do Norte (RN), Paraíba (PB), Pernambuco (PE), Alagoas (AL), Sergipe (SE) and Bahia (BA) (Figure 1). The NEB has high inter-annual and intraannual rainfall variability, with the thermal influence of the Pacific and tropical Atlantic oceans, which control and modulate years of drought and floods in the region [39][40][41].

Fire Foci Data
Fire foci data were obtained from CPTEC/INPE via BDQueimadas [35]. Currently, CPTEC/INPE uses several environmental polar orbit and geostationary satellites in the composition of its observation network for the monitoring of fire foci in South America. The data are available in detail at the BDQueimadas website through the link http://pirandira.cptec.inpe.br/queimadas/, accessed on 20 January 2019. The long-term time series data of fire foci in 1794 municipalities in the NEB during 1998-2018 are analyzed in this study. The data were extracted by an algorithm in annual and monthly format.

Fire Foci Data
Fire foci data were obtained from CPTEC/INPE via BDQueimadas [35]. Currently, CPTEC/INPE uses several environmental polar orbit and geostationary satellites in the composition of its observation network for the monitoring of fire foci in South America. The data are available in detail at the BDQueimadas website through the link http://pirandira. cptec.inpe.br/queimadas/, accessed on 20 January 2019. The long-term time series data of fire foci in 1794 municipalities in the NEB during 1998-2018 are analyzed in this study. The data were extracted by an algorithm in annual and monthly format.

Descriptive and Multivariate Analysis
For the characterization of homogeneous regions based on the records of fire foci of municipalities and states, CA was applied using the Euclidean distance measure [21,23], as shown in Equation (1): In the above equation, X aj is the ajth characteristic of the ath observation; X bj is the jth characteristic of the bth observation; N ab is the number of municipalities. Ward's method [21,23] was used in the study, according to Equation (2): In the above equation, X ijk is the value of the jth variable for the ith of the N k observations of the K groups and X ·jk is the average value. The dot represents the average value of the sub-index.
After obtaining the constituted homogeneous groups, an exploratory analysis of the time series data was carried out, as well as the distribution of data on a monthly and annual scale for each group. This procedure allowed the detection of outliers, based on the expected average variation ranges of the data, and the identification of the temporal variability existing in the data. In descriptive statistics analysis part, the mean, the standard deviation, the frequency, the total and the maximum monthly values of fire foci were computed for each group, using the R programming environment (version 3.4-1).

Vegetation Indices
The images acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor were used for the temporal detection of the heat sockets on the biomes within in the NEB. The MODIS produces 16-day interval images at several spatial resolutions of the Vegetation Index (VI) products, values of the MOD13A1 and MYD13A1 products on the TERRA and AQUA satellites Didan [42]. Moreover, it is obtained for the years 2003-2006 and 2013-2018. The MOD13A1 and MYD13A1 products provide VI values at 500 m spatial resolution per pixel.
The two main vegetation layers, EVI and NDVI, are calculated from Equation (3) and Equation (4), respectively, to assess the attenuation of the effects of the soil. Furthermore, the effects of the Soil Adjusted Vegetation Index (SAVI) and atmosphere (adjusted by the Atmospherically Resistant Vegetation Index (ARVI)) are adjusted by the vegetation mapping [42].
where ρ IVP , ρ V and ρ A are the reflectances in the spectral range of the near-infrared, red and blue, respectively. g (gain factor) is 2.5, c 1 and c 2 are the correction coefficients for red (6) and blue (7.5) atmospheric effects, respectively, and L is the correction factor for soil interference (1).

Precipitation
We used the monthly precipitation Climate Hazards Group InfraRed Precipitation with Stations (CHIRPS) product, the Climate Hazards Group (CHG) from the University of California at Santa Barbara (UCSB) (http://chg.geog.ucsb.edu/data/chirps/index.html, accessed on 20 January 2019) [43,44]. CHIRPS data are nearly real-time, updated, cover almost the entire globe (Earth only, 50 • S-50 • N) and incorporate a 0.05 • resolution satellite (with 5.6 km) at different time scales (daily, monthly, decadal and pentadal). CHIRPS data are based on precipitation estimates from cold infrared (CCD) duration observations calibrated using the Multi-Satellite Precipitation Analysis by Tropical Rain Measurement Mission (TRMM/MSPA) with in situ data for the station from a variety of sources, including national and regional meteorological services [43,44]. In addition, the spatial resolution of CHIRPS is superior compared to other satellite-based global precipitation data sets.

Standardized Precipitation Index (SPI)
In order to identify the occurrence of drought and wet events, we used the SPI based on the CHIRPS precipitation product. The SPI formulation is based on the density and gamma probability function (Equation (5)) that is calculated every month, where α is the shape parameter (α > 0), β is the scale parameter (β > 0), which is determined using the maximum likelihood method, and x is the amount of rainfall, which can vary according to α and β [45].
where Γ(α) is the gamma function.
For the calculation of the SPI, time series of 12-month accumulations were generated and fitted to a two-parameter gamma distribution function. Finally, an equiprobability transformation from the cumulative density functions to the standard normal distribution with mean 0 and variance 1 was performed to obtain the SPI. Therefore, the SPI values are expressed in standard deviations, with positive or negative SPI values indicating greater or less than median precipitation, respectively. Three categories-dry, wet and normal-can be defined based on the SPI values, according to the methodology used by [45].

Temporal
Cluster analysis highlighted the formation of 10 homogeneous sub-groups of fire foci (G1-G10) in the NEB, with a variance of 76.5% and a total number of 6,618,858 fire foci records in the studied period. High variability was observed between the groups in the absolute or relative number of fire foci, mean, standard deviation and maximum monthly values (Table 1). In these groups, G1 has the largest territorial extension, consisting of 1575 municipalities, and covers all biomes existing in the NEB, with a total of 2,130,317 foci (32% of the total events). Moreover, G2 is the second largest group, consisting of 144 municipalities, with a total of 1,459,133 foci (22% of the total number), and mainly covers Cerrado and Caatinga, with less coverage in the Amazon and the Atlantic forest. However, G3 is the smallest territorial extended group, formed of 37 municipalities, with a total of 815,221 foci (12.3% of the total of events), with coverage only in Cerrado. Groups 4 to 10 (G4-G10) consist of 38 municipalities, with a total of 2,214,187 fire foci, and hold the highest percentage (33.7% of the total of events), with coverage in Cerrado and the Amazon. In addition, the municipality of Grajaú (MA) presented the highest number of records (139,807 foci) and the maximum monthly number (18,812 foci) in relation to the other municipalities evaluated in the study. It is worth mentioning that the municipality of Grajaú is seriously affected by burnings in MA, most likely due its livestock and agriculture activities.

Annual
The annual analysis of fire foci is displayed in two distinct periods of the total study data, where period 1 (P1) denotes the years 2003-2006 and period 2 (P2) denotes 2013-2018 ( Figure 2). Both periods were also characterized in the largest (G1 and G2) and smallest (G8 and G10) groups. According to the annual analysis, groups G3-G6 were all part of P2, whereas groups G7 and G9 did not have a defined period and displayed the highest number of outliers. Despite G1 being the largest group in terms of territorial extension, it had a lower number of outliers (1999,2006  The analysis of the monthly fire foci in the NEB using a boxplot exhibited four distinct periods (P3-P6). These four periods refer to three, four, five and six consecutive months of the highest records, respectively ( Figure 3). All periods mentioned deal with the spatiotemporal variability of precipitation in the NEB. P4 corresponded to groups G3, The analysis of the monthly fire foci in the NEB using a boxplot exhibited four distinct periods (P3-P6). These four periods refer to three, four, five and six consecutive months of the highest records, respectively ( Figure 3). All periods mentioned deal with the spatiotemporal variability of precipitation in the NEB. P4 corresponded to groups G3, G6, G7 and G8, and P5 was related to groups G1, G4, G9 and G10, while the exceptions were G2 and G5, with associations with P3 and P6, respectively. In G1, the months showing high variability were August to December (P5), including an outlier in September. Comparable to the annual analysis, G1 had a lower number of outliers, mainly in January, February and April. Lastly, the monthly maximum values identified were greater than 30 foci per month (Figure 3a). In G2, the variable months were from July to December (P6), with outliers in almost each month, with the exceptions of January, February, August, September and October. Regarding the interquartile range, G2 was larger than G1, showing maximum values of more than 300 foci per month (Figure 3b). G3 presented a significant number of outliers, with October being the only exception, similar to the frequency of outliers in the annual scale. G3 had maximum records above 900 foci per month, even though it was lower than the consecutive months (P4) (Figure 3c). G4-G10 were basically formed by P4 and P5. These groups had an excessive number of outliers, excluding September and October. The maximum values varied between 3000 and 5000 foci per month (Figure 3d-i). One exception was the municipality of Grajaú (MA), with a maximum of 7500 foci per month (Figure 3j). Furthermore, the monthly outliers corresponded to monthly droughts or prolonged droughts in the NEB.

Spatial
The spatial distribution of the fire foci showed high variability along the NEB area, mainly on the coast, close to the Borborema Plateau (Figure 1), and within the existing biomes. It is worth noting the formation of a gradient of fire foci from the coast to the interior of the NEB, which is possibly associated with the changes in land use and occupation, and the variability of evapotranspiration (ET) (Figure 4). Group G1 (4.7 ± 6.5 foci) was concentrated close to the coastal region, followed by Agreste and Sertão in the NEB, mainly in the Amazon, Caatinga and Atlantic Forest biomes, with the exception of the Cerrado biome, with the lowest number of records (Figure 4a). Group G2 (40.6 ± 71.7 foci) was partly concentrated in the coastal region, over the states of MA and PI, and in the interior of the NEB, mainly in the Cerrado biome, with the exceptions of the Caatinga and Atlantic Forest biomes, with the smallest records ( Figure 4b). Group G3 (94.7 ± 212.0 foci) was concentrated inside the NEB, including the states of BA, MA and PI, and fully inserted in the Cerrado biome, with the exceptions of the Amazon and the Caatinga biomes ( Figure 4c). Finally, groups G4-G10 were similar to G2, with the concentration of fire foci in part of the coastal region of the states of MA and PI, followed by the northwest (NW) of BA, while the Amazon and Cerrado biomes were exclusive to these groups (Figure 4d).

NDVI
Thematic maps of NDVI and EVI were generated based on the results obtained from the annual analysis (P1 and P2). The spatial variability of environmental degradation from the dense and abundant vegetation categories via NDVI was highlighted in the Amazon and Atlantic Forest biomes in P1, without significant variations in the (NW) portion, only in the Atlantic Forest biome ( Figure 5). On the contrary, P2 exhibited a significant reduction in both biomes. The central region of the NEB showed the predominance of the barren land/built-up/rock categories at the beginning of P1 and along P2, while the sparse vegetation category increased significantly in P2, a period of severe drought in the NEB. Additionally, the category of water bodies remained unchanged regardless of the period evaluated. Comparatively, P1 and P2 showed a significant change in the NEB landscape, in the percentage of categories and in the coverage area.
of outliers in the annual scale. G3 had maximum records above 900 foci per month, even though it was lower than the consecutive months (P4) (Figure 3c). G4-G10 were basically formed by P4 and P5. These groups had an excessive number of outliers, excluding September and October. The maximum values varied between 3000 and 5000 foci per month (Figure 3d-i). One exception was the municipality of Grajaú (MA), with a maximum of 7500 foci per month (Figure 3j). Furthermore, the monthly outliers corresponded to monthly droughts or prolonged droughts in the NEB.

Spatial
The spatial distribution of the fire foci showed high variability along the NEB area, mainly on the coast, close to the Borborema Plateau (Figure 1), and within the existing biomes. It is worth noting the formation of a gradient of fire foci from the coast to the interior of the NEB, which is possibly associated with the changes in land use and occupation, and the variability of evapotranspiration (ET) (Figure 4). Group G1 (4.7 ± 6.5 foci) was concentrated close to the coastal region, followed by Agreste and Sertão in the NEB, was partly concentrated in the coastal region, over the states of MA and PI, and in the interior of the NEB, mainly in the Cerrado biome, with the exceptions of the Caatinga and Atlantic Forest biomes, with the smallest records (Figure 4b). Group G3 (94.7 ± 212.0 foci) was concentrated inside the NEB, including the states of BA, MA and PI, and fully inserted in the Cerrado biome, with the exceptions of the Amazon and the Caatinga biomes ( Figure  4c). Finally, groups G4-G10 were similar to G2, with the concentration of fire foci in part of the coastal region of the states of MA and PI, followed by the northwest (NW) of BA, while the Amazon and Cerrado biomes were exclusive to these groups (Figure 4d).

NDVI
Thematic maps of NDVI and EVI were generated based on the results obtained from the annual analysis (P1 and P2). The spatial variability of environmental degradation from the dense and abundant vegetation categories via NDVI was highlighted in the Amazon and Atlantic Forest biomes in P1, without significant variations in the (NW) portion, only in the Atlantic Forest biome ( Figure 5). On the contrary, P2 exhibited a significant reduction in both biomes. The central region of the NEB showed the predominance of the barren land/built-up/rock categories at the beginning of P1 and along P2, while the sparse vegetation category increased significantly in P2, a period of severe drought in the NEB. Additionally, the category of water bodies remained unchanged regardless of the period In general, the greatest changes were observed in the Caatinga biome, being directly associated with the rainfall pattern, which alters the entire landscape. The vegetation of Caatinga has high sensitivity to changes in weather and climate; with this, any variation in the pluviometric regime causes a direct impact on the vegetation. In addition, human interference contributes to intense degradation, as shown in Figure 5.   In general, the greatest changes were observed in the Caatinga biome, being directly associated with the rainfall pattern, which alters the entire landscape. The vegetation of Caatinga has high sensitivity to changes in weather and climate; with this, any variation in the pluviometric regime causes a direct impact on the vegetation. In addition, human interference contributes to intense degradation, as shown in Figure 5.

EVI
Based on the EVI records, it was possible to identify changes in the vegetation of the biomes within the NEB, particularly in areas with an increase or decrease in vegetation (without vegetation). In the thematic maps, EVI values ranging from 0.3 to 0.5 indicate the presence of areas with vegetation, while areas of decreased or sparse vegetation and exposed soil have EVI values ranging from 0.0 to 0.3 ( Figure 6). Significant accumulation of water, also known as water bodies, where the vegetation is typically poor, characterized by low EVI values (between 0.0 and 0.10), was observed throughout the NEB semiarid region. Low EVI values were observed mainly in the NE/SW areas of the NEB, which corresponds to the Caatinga and Cerrado biomes, respectively, with the intensified loss of vegetation over the years. However, the presence of areas with high vegetation density (>0.40) was observed in the Amazon and Atlantic Forest biomes.

Precipitation (CHIRPS)
The CHIRPS annual rainfall archive showed high interannual variability in the NEB, highlighting the coastal region, with some sectors south, east and NW of the NEB, followed by the border regions of the Amazon and Cerrado biomes, with the highest total annual rainfall (between 1641 and 2200 mm). In the central and northern regions of the NEB, corresponding to the Caatinga and Atlantic Forest biomes, the total annual rainfall varied between 240 and 800 mm (Figure 7). Several of the studied years, located mainly in the central and northern regions of the NEB, were dry-for example, 2003 and the period 2013-2017. It is known that the interannual rainfall variability in the NEB is highly conditioned by the ENOS climate variability mode, being intensified by the Pacific Decadal Oscillation (PDO) and by the influence of the Interhemispheric Gradient of the Atlan-

Precipitation (CHIRPS)
The CHIRPS annual rainfall archive showed high interannual variability in the NEB, highlighting the coastal region, with some sectors south, east and NW of the NEB, followed by the border regions of the Amazon and Cerrado biomes, with the highest total annual rainfall (between 1641 and 2200 mm). In the central and northern regions of the NEB, corresponding to the Caatinga and Atlantic Forest biomes, the total annual rainfall varied between 240 and 800 mm (Figure 7). Several of the studied years, located mainly in the central and northern regions of the NEB, were dry-for example, 2003 and the period 2013-2017. It is known that the interannual rainfall variability in the NEB is highly conditioned by the ENOS climate variability mode, being intensified by the Pacific Decadal Oscillation (PDO) and by the influence of the Interhemispheric Gradient of the Atlantic Sea Surface Temperature (GITSMA).

Discussion
According to the CA, the behavior of the fire foci could be categorized into 10 homogeneous groups (Table 1). It is worth mentioning that the CA technique has been applied in several studies of fire foci, with satisfactory results in Brazil [21,23] and worldwide [16,46]. The group G1 is located on the coast and in the interior of the NEB, and consists of the Caatinga, the Atlantic Forest and the Amazon biomes, emphasizing the lower occurrence in Cerrado. Group G1 also showed the highest population density, the lowest Human Development Indexes (HDI) [47,48] and the strongest LULC [49,50], being identified with the influence of socioeconomic data and LULC in other studies on fires around the world [2,[4][5][6]10,11], since it is a region where large and rain-fed agriculture is concentrated [51,52]. The group G2 is located on a part of the north and inland NEB coast, corresponding to the Caatinga, the Cerrado and the Amazon biomes. Both G1 and G2 make up most of the NEB municipalities and are strongly influenced by the variability of ET, since the groups formed of ET are similar to groups of fire foci [53,54]. Rocha Júnior et al. [54] evaluated the influence of potential evapotranspiration (PET) on Northeastern Cerrado, and identified that the highest values occurred in the north of MA, with values ranging between 1250 and 2000 mm per year, while, in the Cerrado portion located in Bahia, the values oscillate between 600 and 1500 mm/year. The group G3 is fully inserted within Northeastern Cerrado, while the remaining groups, G4-G10, are fully inserted within the Amazon and Northeastern Cerrado biomes. The rainfall dynamics in groups G4-G10 are strongly conditioned by the biomes [55,56] and by the intense LULC in the border regions [57], mainly by the MATOPIBA agricultural consortium (acronym for the states: MAranhão, TOcantins, PIauí and BAhia) [58,59].

Discussion
According to the CA, the behavior of the fire foci could be categorized into 10 homogeneous groups (Table 1). It is worth mentioning that the CA technique has been applied in several studies of fire foci, with satisfactory results in Brazil [21,23] and worldwide [16,46]. The group G1 is located on the coast and in the interior of the NEB, and consists of the Caatinga, the Atlantic Forest and the Amazon biomes, emphasizing the lower occurrence in Cerrado. Group G1 also showed the highest population density, the lowest Human Development Indexes (HDI) [47,48] and the strongest LULC [49,50], being identified with the influence of socioeconomic data and LULC in other studies on fires around the world [2,[4][5][6]10,11], since it is a region where large and rain-fed agriculture is concentrated [51,52]. The group G2 is located on a part of the north and inland NEB coast, corresponding to the Caatinga, the Cerrado and the Amazon biomes. Both G1 and G2 make up most of the NEB municipalities and are strongly influenced by the variability of ET, since the groups formed of ET are similar to groups of fire foci [53,54]. Rocha Júnior et al. [54] evaluated the influence of potential evapotranspiration (PET) on Northeastern Cerrado, and identified that the highest values occurred in the north of MA, with values ranging between 1250 and 2000 mm per year, while, in the Cerrado portion located in Bahia, the values oscillate between 600 and 1500 mm/year. The group G3 is fully inserted within Northeastern Cerrado, while the remaining groups, G4-G10, are fully inserted within the Amazon and Northeastern Cerrado biomes. The rainfall dynamics in groups G4-G10 are strongly conditioned by the biomes [55,56] and by the intense LULC in the border regions [57], mainly by the MATOPIBA agricultural consortium (acronym for the states: MAranhão, TOcantins, PIauí and BAhia) [58,59].
The annual fire foci boxplot for the NEB identified two periods, P1 (2003)(2004)(2005)(2006) and P2 (2013-2018) (Figure 2), both classified within the largest (G1 and G2) and the smallest groups (G8 and G10) of fire foci. It is worth mentioning that both periods are influenced by the phases of ENSO and the variability of GITSMA, which influence rain dynamics in the NEB [39,45]. During the El Niño years, severe droughts occur, as discussed by Marengo et al. [60], while, during the La Niña years, large rainfall occurs [61,62], being duly characterized by the CHIRPS product (rainfall and SPI) (Figures 7 and 8). Due to the lack of consistent rainfall data in the NEB, the CHIRPS data were validated for the region, with satisfactory results [62]. It is known that the distribution of fire foci has a direct relationship with the inter-annual variability of rainfall in Brazil, and mainly in the NEB area [21,23,30], and this relationship is striking in the existing biomes in the NEB [22,56]. The monthly boxplot exhibited four distinct periods (P3, P4, P5 and P6-3, 4, 5 and 6 months, respectively) of the fire foci in the NEB (Figure 3). The abovementioned periods are influenced by the spatiotemporal variability of precipitation in the NEB [39][40][41]45]. Both periods defined in the NEB show the variability of the fire foci, where the largest (smallest) records of fire foci correspond to the dry season (transition and rainy seasons) [21,37,38], and it is therefore directly associated with the rainfall intra-annual variability and the multiscale meteorological systems that operate in the NEB [21][22][23]42]. In both the annual and the monthly scales, outliers that corresponded to severe or prolonged droughts in the NEB were identified [38,45,46].
The environmental degradation was identified by the NDVI and EVI vegetation indexes, widely used in the scientific literature [23,24,26,29,30]. The dense vegetation and the abundance categories, as identified via the NDVI, occurred in the Amazon and Atlantic Forest biomes in P1; however, in P2, a significant reduction in both biomes was recorded. Other studies [30,49,50] identified that drought events across the country were more severe and widespread between 2011 and 2017 based on the Integrated Drought Index (IDI), most likely owing to the increase in extreme drought events in P2 (El Niño and neutral years) [60][61][62][63][64][65], and the intensification of agricultural production of commodities (corn and soybeans) in the MATOPIBA (Amazon) consortia [58,59] and the SEALBA (acronym for the states: SErgipe, ALagoas and BAhia) in the region of the coastal NEB [30,57]. In the central region of the NEB, a high proportion of the barren land, built-up and rock categories were registered at the beginning of P1 and along P2; however, the sparse vegetation increased significantly in P2. Notably, the central region, which includes the Caatinga biome, is characterized by a highly variable rainfall regime [53,55,56], and areas under a desertification process [47,50]. It is known that during extreme droughts, food security in the NEB semiarid region is under serious risk due to agricultural collapse, as highlighted recently by [40,64,65]. There are different statistical analyses for estimating fire foci throughout the world during long-term and short-term periods [62,[66][67][68].
The limitations of the study are as follows: (i) it does not have a historical series of hotspots over 30 years and (ii) we do not have an ROS for the study region. The main recommendations for other studies are (i) to relate the hotspots to socioeconomic and environmental indicators, as well as to compare them with other drought indices existing in the literature, such as the PDSI and Reconnaissance Drought Index (RDI) among others; (ii) to characterize hotspot events with meteorological systems that can intensify fires in biomes-for example, Upper Tropospheric Cyclonic Vortices (UTCV) and the South Atlantic subtropical high (SASH).

Conclusions
Cluster analysis of the fire foci of biomes within in the NEB revealed 10 homogeneous subgroups. The different groups of fire foci presented both space and time variability, depending on the inter-annual and intra-seasonal rain regimes and alterations associated with land use and occupation. G1 covers all biomes present in the NEB, unlike G2, which mainly covers the Cerrado and the Caatinga biomes. The G3 group covered only Cerrado, and G4-G10 cover both the Cerrado and Amazon biomes. It is worth mentioning the fire foci gradient formation starting from the coast towards the inner areas of the NEB. Such a gradient is directly linked to changes in land use and occupation, which were identified correspondingly by the vegetation indexes accounting for the sparse vegetation category (NDVI) and areas deprived of vegetation (EVI). This study made it possible to identify the loss of the vegetation over the years in the Caatinga and Cerrado biomes, unlike the significant increase in vegetation in the Amazon and Atlantic Forest biomes. The fire foci detected in the Cerrado biome are related to the expansion of agricultural consortia (MATOPIBA), while the fire foci detected in the Atlantic Forest biome and coastal areas are due to the SEALBA consortium.