Assessment of the Characteristics of Recent Major Wildﬁres in the USA, Australia and Brazil in 2018–2019 Using Multi-Source Satellite Products

: This study analysed the characteristics of the recent (2018–2019) wildﬁres that occurred in the USA, Brazil, and Australia using Moderate Resolution Imaging Spectroradiometer (MODIS) active ﬁres (AF), ﬁre radiative power (FRP, MW) and burned area (BA) products. Meteorological and environmental parameters were also analysed. The study found various patterns in the spatial distribution of ﬁres, FRP and BA at the three sites, associated with various vegetation compositions, prevailing meteorological and environmental conditions and anthropogenic activities. We found signiﬁcant ﬁre clusters along the western and eastern coasts of the USA and Australia, respectively, while vastly distributed clusters were found in Brazil. Across all sites, signiﬁcant ﬁre intensity was recorded over forest cover (FC) and shrublands (SL), attributed to highly combustible tree crown fuel load characterised by leafy canopies and thin branches. In agreement, BA over FC was the highest in the USA and Australia, while Brazil was dominated by the burning of SL, characteristic of ﬁre-tolerant Cerrado. The relatively lower BA over FC in Brazil can be attributed to fuel availability and proximity to highly ﬂammable cover types such as cropland, SL and grasslands rather than fuel ﬂammability. Overall, this study contributes to a better understanding of wildﬁres in various regions and the underlying environmental and meteorological causal factors, towards better wildﬁre disaster management strategies and habitat-speciﬁc ﬁreﬁghting.


Introduction
Fires are essential for shaping most ecosystems around the globe; thus, they are crucial for conservation and ecosystems management. In most rangelands and forest ecosystems, fires are used to contain the extensive bush encroachment and invasive species [1], and for improving grasslands productivity [2]. Moreover, fire is an affordable clearance mechanism for agricultural expansion in most rural communities, especially in developing countries, where the adoption of machinery is limited due to the high cost and poor living standards [3,4]. In other communities, fires are a fundamental part of culture and heritage. While the above-mentioned fires are related to intentional or accidental anthropogenic activities, fires can also be ignited naturally through lighting and extreme and prolonged heat conditions. Regardless of the ignition source, wildfires can rapidly become widespread, uncontrollable, intense and disastrous due to the prevailing local meteorological and environmental territory, located at latitudes 34.095° S to 5.464° N and longitudes 74.431° W to 34.194° W. Brazil consists of six biomes, i.e., Amazon, Caatinga, Cerrado, Atlantic Forest, Pampa and Pantanal, with varying climatic patterns and fire regimes. Much of Brazil, i.e., 81.4% (including the Amazon), is characterised by a tropical climate [38]. The Amazon has a Humid Equatorial climate with annual precipitation of between 2050 mm and 2650 mm. Towards the south, the region is occupied by indigenous people, and dominated by increasing livestock and cultivation activities. Cerrado biome is characterised by a humid climate, with an average annual rainfall of between 800 to 2000 mm [39] during the wet season. The dry season in the Brazilian Cerrado is long, i.e., from May to September or October; thus, increases fuel flammability. Consequently, fires from natural (i.e., lightning) and anthropogenic (i.e., grazing management) sources are common during the dry season [40]. In the southernmost part of the country, the subtropical grasslands, i.e., known as Pampas, exists under a humid climate and borders the Atlantic Forest, which also has higher precipitation throughout the year, i.e., an annual average of 1800 mm. Consequently, the fires are very low in Pampas and Atlantic Forest biomes [41]. The 2019 fires resulted in a diplomatic crisis and controversy [33].
The third site constituted eastern Australia (Queensland, New South Wales, Victoria and Tasmania), located at latitudes 44.092° S to 9.872° S and longitudes 138.177° E to 153.864° E. The dry seasons are typically divided into early (i.e., April-July) and late (August-November) dry seasons in northern Australia [42]. Generally, eastern Australia is characterised by tropical Savannas and tussock grasslands (mainly in Queensland), and temperate forests in southeastern Australia. Both Savannas and grasslands burn extensively, on an annual basis, between March and December [43], as a result of declining rainfall from ~2000 mm to <400 mm and anthropogenic activities. Occasionally, i.e., around November-December, fires may be lit naturally by lightning due to the onset of monsoonal conditions [44]. The 2019 fires over temperate forest were unprecedented and reportedly the worst (i.e., largest) fires since the Black Friday fires in 1939 [35,45]. The fire season, in this region, starts from September to February.  The third site constituted eastern Australia (Queensland, New South Wales, Victoria and Tasmania), located at latitudes 44.092 • S to 9.872 • S and longitudes 138.177 • E to 153.864 • E. The dry seasons are typically divided into early (i.e., April-July) and late (August-November) dry seasons in northern Australia [42]. Generally, eastern Australia is characterised by tropical Savannas and tussock grasslands (mainly in Queensland), and temperate forests in southeastern Australia. Both Savannas and grasslands burn extensively, on an annual basis, between March and December [43], as a result of declining rainfall from~2000 mm to <400 mm and anthropogenic activities. Occasionally, i.e., Remote Sens. 2020, 12, 1803 4 of 18 around November-December, fires may be lit naturally by lightning due to the onset of monsoonal conditions [44]. The 2019 fires over temperate forest were unprecedented and reportedly the worst (i.e., largest) fires since the Black Friday fires in 1939 [35,45]. The fire season, in this region, starts from September to February.

Fire Parameters
Moderate Resolution Imaging Spectroradiometer (MODIS) Burned Area (MCD64A1) and Global Monthly Fire Location (MCD14ML) products at 500 m and 1 km [37,46], respectively, were used to characterise the wildfires at the three study sites, i.e., western USA, Brazil and eastern Australia. Specifically, burned area (BA), burn date (BD) and Fire Radiative Power (FRP) were used to determine the extent of the burned areas (km 2 ), the date of burning (DOY) and the fire intensity (i.e., FRP, MW), respectively. The MCD64A1 product provides monthly global BA with per-pixel quality information. This product is derived using the Collection 6 (C6) algorithm that uses a combination of Level 2G daily surface reflectance products (i.e., MOD09GHK and MYD09GHK) at 500 m spatial resolution, Level 3 active fire products (MOD14A1 and MYD14A1) at 1 km replicated to 500 m, and the most recent MODIS annual land cover product (i.e., MCD12Q1) at 500 m resolution. Generally, the C6 algorithm characterises burned and non-burned pixels based on persistent changes identified from a time series of burn-sensitive vegetation indices. Then, the spatial and temporal active fire information and burn-sensitive vegetation indices are used to generate dynamic thresholds that are applied to the composite data to determine if the pixel is burned or unburned. The BD is a single data layer consisting of the day-of-year (DOY) on which the burn occurred [37,46,47]. Collection 6 BA algorithm is an improvement of Collection 5.1 MCD64A1 and Collection 5.1 MCD45A1 mapping algorithms. According to Giglio et al. [47] and Humber et al. [48], MCD64A1.006 has a better detection rate for small burns and has reduced temporal uncertainty and extent of unmapped areas.
On the other hand, the MCD14ML product consists of geographic locations (i.e., centroids of 1-km pixels) of daily (i.e., 24-h) active fires and their FRP, based on an approach by Wooster et al. [28] and Wooster et al. [49], detected by MODIS sensors on-board Terra and Aqua satellites, at the time of the satellite overpass, based on the Collection 6 algorithm. The Collection 6 algorithm identifies one or more actively burning fires based on native brightness temperatures from MODIS bands centred at 4 µm, 11 µm and 12 µm and reflectance centred at 0.65 µm, 0.86 µm and 2.1 µm for daytime observations. On the other hand, FRP estimation is based on spectral radiances of a fire pixel and non-fire (i.e., background) pixels in the mid-infrared (MIR) band. Further details can be found in Wooster et al. [28].
Using 2500 images from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Giglio et al. [37] found the improved performance of the C6 AF algorithm. Specifically, the authors found reductions in omission errors over large fires, and commission error rates, associated with small forest clearings, in the Brazilian Amazon and other tropical ecosystems relative to previously reported errors by Schroeder et al. [50] and Cheng et al. [51]. The BA and AF products were accessed from Google Earth Engine (GEE) and LANCE Fire Information for Resource Management System (FIRMS, https://earthdata.nasa.gov/data/near-real-time-data/firms), respectively. AF points with < 10% detection confidence were discarded from the analysis.

Environmental Parameters
Vegetation forms the primary fuel load for wildfires. Thus, their frequency and intensity are influenced by vegetation availability, type and condition in a particular period [52]. In this study, the vegetation cover type was used to assess the extent of destruction caused by wildfires on the Remote Sens. 2020, 12, 1803 5 of 18 ecosystems. Specifically, the Vegetation cover types at 300 m spatial resolution were derived from the European Space Agency (ESA) Climate Change Initiative-Land Cover (CCI-LC) of 2018 [53]. The ESA CCI-LC was reclassified into five broad vegetation cover type classes (see Figure 1), namely, Cultivated land (CL), Forest cover (FC), Shrubland (SL), Grassland (GL) and other vegetation covers (OVC), while non-vegetation classes were discarded. The ESA CCI-LC is generated at a spatial resolution of 300 m, from global daily reflectance measurements acquired by MEdium Resolution Imaging Spectrometer (MERIS), Advanced Very High Resolution Radiometer (AVHRR), Satellite Pour l'Observation de la Terre Vegetation (SPOT-VGT), PROBA-Vegetation (PROBA-V) and Sentinel-3. The classification of land cover is based on the United Nations Land Cover Classification System (UN-LCCS) containing 37 classes, unsupervised and supervised algorithms and time-series data [53]. The use of multi-temporal, multi-sensor data allowed backwards and forward cross-checking to detect and verify land cover changes and eliminate false changes occurring due to inter-annual variability [54,55]. The ESA CCI-LC land cover product for 2018 was downloaded from http: //maps.elie.ucl.ac.be/CCI/viewer/download.php.
Moreover, the vegetation condition index (VCI) [56] from the Copernicus Global Land Service (CGLS, https://land.copernicus.eu/) was used to analyse the changes in vegetation conditions (%) for a period coinciding with the period of the wildfires. VCI is a measure of the current vegetation condition compared to historical maximum and minimum normalised difference vegetation index (NDVI). As Kogan [57] note, the VCI better portrays the impacts of weather dynamics on relative vigour of the vegetation. Consequently, VCI is one of the widely used satellite-based indices for detecting the onset, intensity, duration, and impacts of drought [56,[58][59][60]. CGLS VCI is provided at 1-km spatial resolution and based on a synthesis of 10-day NDVI: where VCI i is the VCI for period i, NDVI i is the observed 10-day NDVI values over period i and NDVI max and NDVI min are the long-term maximum and minimum NDVI, respectively.

Meteorological Parameters
The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) is a NASA reanalysis product using a new version of the Goddard Earth Observing System Data Assimilation System version 5 (GEOS-5). The Global Modelling and Assimilation Office (GMAO) uses its GEOS-5 atmospheric data assimilation system (ADAS) to synthesise the different observations collected over the satellite era into a dataset that is consistent over time as it uses a fixed assimilation system. MERRA-2 is created with version 5.2.0 of the GEOS-5 ADAS with a 0.5 • latitude × 0.625 • longitude × 72 layers model configuration [61]. The bottom 32 layers are terrain following, while the remaining model layers from 164 hPa to 0.01 hPa are constant pressure surfaces [62]. MERRA-2 assimilates bias-corrected AOD from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Very High-Resolution Radiometer (AVHRR) instruments [63,64]. More details on MERRA-2 can be found in Gelaro et al. [61]. In this work, MERRA-2 is used to produce precipitation (Prec.) and wind speed (Ws) data, averaged over the study sites.
Relative humidity (RH) and surface temperature (Ts) were measured using the Atmospheric Infrared Sounder (AIRS). AIRS Infrared Radiances (IR) have a spatial resolution of 13.5 km at nadir, while AIRS Visible and Near-Infrared Radiances have a spatial resolution of 2.3 × 1.8 km (across-track, along-track). With 2378 spectral channels, AIRS has a spectral resolution more than 100 times greater than previous infrared sounders and provides more accurate information on the vertical profiles of atmospheric temperature and moisture. MERRA-2 and AIRS data were accessed from the Geospatial Interactive Online Visualization and Analysis Infrastructure (Giovanni, http://giovanni.gsfc.nasa.gov) and averaged hourly (i.e., Prec. and Ws) and daily (i.e., Ts and RH) for further analysis.

Spatial Patterns of Fires
Anselin's local Moran's I statistic [65] is widely used for identifying statistically significant local spatial patterns (i.e., association, dispersion or randomness) and types of spatial clusters. A positive (i.e., >0) Moran's I statistic indicates clustered fires or a positive correlation in the spatial distribution of fires. In contrast, negative values (i.e., <0) indicate dispersed fires or a negative correlation in the spatial distribution of fires, and the values closer to zero indicate randomly distributed fires. In this study, spatial points of active fires (AF) and associated FRP, i.e., representing centroids of MODIS 1-km pixels, were used to identify spatial patterns and types of clusters in western USA (i.e., August-September 2018), Brazil (i.e., August-September 2019), and eastern Australia (i.e., November-December 2019). The z-score and p-values for each AF (i.e., the point representing centroid for 1-km pixel) were used to test the null hypothesis (H 0 ) of no spatial association between the AF. Consequently, two types of statistically significant spatial clusters, i.e., high-high (HH) and low-low (LL), and two statistically significant spatial outliers, i.e., high-low (HL) and low-high (LH) were identified. The HH (LL) clusters indicate several AF with similar high (low) intensity (FRP, MW) values, while the outliers, HL (LH) indicate a high (low) AF surrounded by several low (high) intensity values. Moreover, five fire categories proposed by Ichoku et al. [36] based on FRP (to facilitate fire rating by strength), i.e., category 1 (<100 MW), category 2 (100 to <500 MW), category 3 (500 to <1000 MW), category 4 (1000 to <1500 MW) and category 5 (≥1500 MW), were adopted in this study, to characterise the strength of fires (i.e., fire intensity) at the three study sites.

Statistical Analysis
A non-parametric test, i.e., Kruskal-Wallis ANOVA, was used to determine whether the fire intensity differed according to different vegetation types. Specifically, we tested the null hypothesis (H 0 ) that there is no significant difference between vegetation type-specific fire intensity, i.e., FRP per vegetation type. Moreover, Pairwise-Wilcox Test was used to determine whether class wise-comparisons of vegetation type-specific fire intensities (FRP) were different.

Analysis of Active Fires and Their Intensity over Various Vegetation Types
The results of the spatial distribution of Active Fires (AF) and Fire Radiative Power (FRP, MW) for the three study sites (Figure 2) show the vast distribution of high-density fires across all study sites. These dense fires indicate the areas (i.e., hotspots) where these fires were intense and became disastrous. For example, the fires with the highest FRP were found around 40 • and 48 • N and 120 • W (i.e., Washington and California) in the western USA, 10 • S and 54 • W in Brazil and 30 • S and 150 • E (i.e., New South Wales) in eastern Australia. Statistically, Moran's I statistic show the presence of statistically significant clusters and outliers, at all the three sites, at 95% confidence level (see Table 1). The results show that the most intense fire clusters (HH), with an average FRP of 762.02 MW (i.e., category 3) were recorded in western USA, followed by eastern Australia with an average FRP of 494.37 MW (i.e., category 2) and Brazil with an average FRP of 356.77 MW (i.e., category 2). In contrary, the less intense fire clusters (LL) with an average of 7.45 MW (i.e., category 1) were recorded in Brazil, followed by Australia with an average of 18.13 MW (i.e., category 1) and the USA with an average of 59.27 MW (i.e., category 1). Remote Sens. 2020, 12, x FOR PEER REVIEW 7 of 19  It is also evident from the summary statistics (Table A1) that forest cover (FC) had the highest fire intensity (FRP, MW), followed by shrublands (SL) across all sites, then cultivated land (CL) in Brazil and Australia and GL in the USA. The higher fire intensity in FC and SL can be attributed to highly combustible tree crown fuel load characterised by leafy-closed canopies and thin branches. Such fires can rapidly propagate, given favourable conditions, to nearby crowns especially in areas with non-fragmented vegetation types such as the Forests. This is consistent with Vadrevu et al. [66], who found that higher fire intensity for closed broadleaved evergreen and semi-deciduous forests in India than other vegetation cover types based on empirical cumulative distribution functions.  It is also evident from the summary statistics (Table A1) that forest cover (FC) had the highest fire intensity (FRP, MW), followed by shrublands (SL) across all sites, then cultivated land (CL) in Brazil and Australia and GL in the USA. The higher fire intensity in FC and SL can be attributed to highly combustible tree crown fuel load characterised by leafy-closed canopies and thin branches. Such fires can rapidly propagate, given favourable conditions, to nearby crowns especially in areas with non-fragmented vegetation types such as the Forests. This is consistent with Vadrevu et al. [66], who found that higher fire intensity for closed broadleaved evergreen and semi-deciduous forests in India than other vegetation cover types based on empirical cumulative distribution functions. Depending on their severity (not studied here), fires may result in higher tree mortality and permanent land-use conversion, especially in Brazil, where food production needs (i.e., agricultural expansion) surpass conservation [41]. Similarly, the number of AF over FC at the three sites were the highest, followed by SL and CL for the western USA and Brazil and GL in eastern Australia (Table A1). In contrast, grassland (GL) and other vegetation covers (OVC) showed the lowest fire intensity (FRP, MW) across all sites. This observation is expected since herbaceous vegetation types such as GL have relatively lower above-ground-biomass (AGB) than tree cover types such as FC and SL. Overall, the results indicate that the number of AF and fire intensity characteristics at the three sites varied by vegetation cover type. In fact, the null hypothesis that there is no significant difference in vegetation type-specific fire intensity was rejected at 95% confidence level (p-value < 2.2 × 10 −16 ) at all sites.
The class-wise comparisons of vegetation type-specific fire intensity using Pairwise Wilcox test ( Table 2) reveal even more exciting results, showing a significant difference in fire intensities between CL vs. GL, GL vs. FC, FC vs. OVC and FC vs. SL across all sites. In contrast, the fire intensities between CL vs. FC and GL vs. SL were significantly different in Brazil and eastern Australia only (see Table 2). These differences were mainly driven by different landscape compositions, such as different vegetation types and conditions, temperature and various sizes of the fires at the three sites. For example, the areas with a high number of AF in Brazil and eastern Australia are characterised by an extensive mosaic of CL, GL, SL and FC, while the predominant AF in western USA was over FC and SL. Therefore, fire intensity over GL was not significantly different from SL and OVC in western USA and OVC only in eastern Australia due to the low abundance of these vegetation types. The observed differences were due to relatively lower herbaceous biomass in GL and CL, which had corresponding lower intensity than FC and SL, which are characterised by thick, woody stems and branches and high density of leaves. When drier, these woody vegetation characteristics result in higher flammability and combustibility; thus, the propagation of wildfires to the broader spatial extent subject to fuel availability and favourable meteorological conditions. Although Brazilian fires were found to be intense, i.e., category 2, in this study, a recent related study [33] found that it is within long-term normal, while Kelley et al. [67] showed that meteorological conditions had an insignificant role, i.e., <7%, in their propagation. Previous studies [68,69], based on MCD14ML.005, reported possible omission errors in the detection of active fires as a result of clouds and smoke plumes. It is therefore acknowledged that these errors are inevitable; however, expected to be reduced due to improvements in MCD14ML.006 product [37] used in this study. Other sources of error, such as the omission of smaller fires, are probable; however, they were not critical in this study since the focus was on major wildfires characterised by large fires that are detectable at MODIS 500 m spatial resolution. Contrarily, errors in the spatial location of AF and detection efficiency of sensors may be problematic for fire authorities seeking to locate and extinguish actively burning fires and identify the burning vegetation types for habitat-specific fire fighting, to protect endangered fauna and flora. For example, the exact locations of AF are unknown, since each AF point represents the centroids of 1-km MODIS pixels; thus, several fires may be burning within 1-km pixel as the size of the fire relative to pixel size is always small [50]. Moreover, Hyer and Reid [70] showed that using MODIS, an accuracy of identifying the correct vegetation type of a single AF is approximately 88%. This error may be exacerbated by classification errors in land cover maps. Although algorithmic advances for characterising AF at Landsat resolution were achieved [71], the current temporal configurations (i.e., repeat cycle of ≥ 5 days) of medium to high resolution sensors limits their operational application for fire management.

Spatio-Temporal Variations in the Burned Area over Various Vegetation Types
The results (Figure 3) show various peaks in the burned area for specific dates. In western USA, three prominent peaks were observed in August 2018, with the highest BA per day of~500 km 2 . In contrast, the BA, in Brazil, increased consistently in August 2019 and peaked to the highest (i.e.,~4300 km 2 ) in September 2019. Although extensive BA are common in Brazil [68,72], this finding is consistent to Lizundia-Loiola et al. [33], who found that the total BA from these fires was more than twice the previous year (i.e., 2018); however, it was similar to an average BA over 17 years. Therefore, the 2019 Brazilian fires BA accumulation and distribution were not extraordinary. On the other hand, the highest peaks, of up to~2700 km 2 in BA per day, were recorded in November and December 2019 in eastern Australia. Comparatively, the peaks in BA per day in the western USA were the lowest among the three sites, followed by eastern Australia, while the highest peaks were recorded in Brazil. The observed higher daily BA can be attributed to its vast area with high continuous fuel load relative to other sites. Moreover, August and September are at the end of the dry season in most of Brazil, providing favourable conditions, i.e., drier and expansive fuel load, for propagation wildfires. Since these fires are mostly anthropogenic, i.e., mainly related to common clearing activities and pasture management [67], little or no effort is taken to combat their spread until they encroach on protected areas. On the other hand, western USA fires were confined to smaller areas, attributable weaker winds (Ws) (i.e., ranging between 4 to 6 m s −1 , see Figure 3 and Figure A1) that allowed prolonged burning of various vegetation types. Similar to Brazil, the fires in the western USA occurred mainly towards the end of the season, when preceding higher temperature have dried-out the vegetation. On the other hand, the fluctuating daily BA in eastern Australia, i.e., several peaks (maxima) and troughs (minima), can be attributed to varying fuel (i.e., vegetation) type and conditions. For example, the observed higher temperature (i.e., 30 • C to 40 • C) increased evapotranspiration, resulting in drier twigs and leaves especially in eucalyptus forest, thus supporting propagation of fires and higher BA [45]. Conversely, patches of rainforests impede the spread due to higher moisture content [45]. Across all sites, Ws were weak (i.e., 4 to 8 m s −1 ); thus they cannot be considered an important driver to the propagation of these fires (during the period under study).
Consequently, the fire scars (see Figure 4) were also variable per site. In agreement with the results in Figure 3, the fire scars were vastly distributed in Brazil and eastern Australia. In contrast, the fire scars were the lowest and widely scattered in western USA. In Brazil, the fire scars correspond to areas characterised as Cerrado (i.e., Brazilian Savanna) and fire-sensitive ecosystems such as the Amazon and Atlantic rainforests [52]. This is expected since these biomes occupy the majority of Brazilian territory, i.e., more than 70%. This finding is consistent with Moreira de Araújo et al. [68], who found the largest proportion of burned areas within Brazilian Cerrado and Amazon Rainforest, i.e., 73% and 14%, respectively. While Cerrado vegetation is more adapted to regular fires, these wildfires may result in considerable damage and loss of the Amazon and Atlantic rainforests which are not tolerant to fires [52].
On the other hand, the fluctuating daily BA in eastern Australia, i.e., several peaks (maxima) and troughs (minima), can be attributed to varying fuel (i.e., vegetation) type and conditions. For example, the observed higher temperature (i.e., 30°C to 40°C) increased evapotranspiration, resulting in drier twigs and leaves especially in eucalyptus forest, thus supporting propagation of fires and higher BA [45]. Conversely, patches of rainforests impede the spread due to higher moisture content [45]. Across all sites, Ws were weak (i.e., 4 to 8 m s -1 ); thus they cannot be considered an important driver to the propagation of these fires (during the period under study). Consequently, the fire scars (see Figure 4) were also variable per site. In agreement with the results in Figure 3, the fire scars were vastly distributed in Brazil and eastern Australia. In contrast, the fire scars were the lowest and widely scattered in western USA. In Brazil, the fire scars correspond to areas characterised as Cerrado (i.e., Brazilian Savanna) and fire-sensitive ecosystems such as the Amazon and Atlantic rainforests [52]. This is expected since these biomes occupy the majority of Brazilian territory, i.e., more than 70%. This finding is consistent with Moreira de Araújo et al. [68], who found the largest proportion of burned areas within Brazilian Cerrado and Amazon Rainforest, i.e., 73% and 14%, respectively. While Cerrado vegetation is more adapted to regular fires, these wildfires may result in considerable damage and loss of the Amazon and Atlantic rainforests which are not tolerant to fires [52]. As can be observed in Figure 5, SL, characteristic of Brazilian Cerrado vegetation, burned extensively with a BA of ~28,000 km 2 , followed by FC with ~14,000 km 2 , which also had the highest number of AF, then CL with ~8000 km 2 of BA. It is anticipated that the fires, commonly used for grassland management, hunting activities and 'slash and burn' agricultural practices in Brazil [52,73], spread into adjacent SL and FC, where they became difficult to control and widespread. However, the higher number AF in FC suggests relatively smaller and isolated fires that could not spread to more extensive areas due to constantly wet conditions in the Amazon rainforest. This is consistent As can be observed in Figure 5, SL, characteristic of Brazilian Cerrado vegetation, burned extensively with a BA of~28,000 km 2 , followed by FC with~14,000 km 2 , which also had the highest number of AF, then CL with~8000 km 2 of BA. It is anticipated that the fires, commonly used for grassland management, hunting activities and 'slash and burn' agricultural practices in Brazil [52,73], spread into adjacent SL and FC, where they became difficult to control and widespread. However, the higher number AF in FC suggests relatively smaller and isolated fires that could not spread to more extensive areas due to constantly wet conditions in the Amazon rainforest. This is consistent with Lizundia-Loiola et al. [33], who found newly burned pixels in the Amazon tropical forest in the 2001-2019 period, linked to new policies supporting agricultural expansion. Therefore, relatively lower BA over FC in Brazil can be attributed to fuel availability rather than flammability, as the forests have a deeper root system to withstand dry conditions. Nevertheless, the amount of FC loss to wildfires, in Brazil, was marked and similar to that in eastern Australia, i.e.,~16,000 km 2 . This is alarming considering that the entire Brazilian territory was analysed, relative to only a portion, i.e., eastern Australia. Similarly, FC was the most extensively burned vegetation type in the western USA with 1800 km 2 . However, it was markedly lower (i.e., six to seven times lower) than the burned FC in Brazil and Australia, respectively. In western USA and eastern Australia, the BA per vegetation type mostly corresponded to the high number of AF and higher fire intensity. It is expected that the antecedent higher Ts and Ws and low RH resulted in a rapid loss in canopy moisture and dry-up of the herbaceous plant components such as leaves, branches and grasses [74] (see Figure A2). At the time of the fires, higher Ts (i.e., 30 • C to 40 • C) coupled with lower Prec. (i.e., <0.4 mm hr −1 ), and availability of continuous fuel (Figures 1 and 3) influenced the establishment of new fires, propagation and intensity of the fires in FC, SL and GL, as evidenced by the higher number and clustering of AF and fire intensity (Figure 2, Table 1), and extensive losses in vegetation cover ( Figure 5). In fact, at the three sites, moderate to extreme dry (or drought) conditions were evident ( Figure A2), causing drier leaves. This is consistent with Nolan et al. [45], who showed that the widespread fires in eastern Australia were driven by dry fuel conditions, caused by prolonged drought, i.e., before and during the fire. In another study, Turco et al. [75] show that fire occurrence is significantly related to drought conditions. Similarly, Moreira de Araújo et al. [68] found that the higher concentration of fires in 2007 and 2010 over Brazilian Cerrado and Amazon Rainforest was related to low rainfall resulting from La Niña climatic pattern. Contrarily, Kelley et al. [67] showed that the 2019 Brazilian fires had a low meteorological influence, i.e., <7%, and concluded that land management was the most probable driver of the increase in the burned area. Overall, the interaction of various factors, i.e., meteorological, vegetation availability and conditions and human activities influenced the ignition, flammability, combustion, spatial patterns and intensity of wildfires at the three sites.
Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 19 lower BA over FC in Brazil can be attributed to fuel availability rather than flammability, as the forests have a deeper root system to withstand dry conditions. Nevertheless, the amount of FC loss to wildfires, in Brazil, was marked and similar to that in eastern Australia, i.e., ~16,000 km 2 . This is alarming considering that the entire Brazilian territory was analysed, relative to only a portion, i.e., eastern Australia. Similarly, FC was the most extensively burned vegetation type in the western USA with ~1800 km 2 . However, it was markedly lower (i.e., six to seven times lower) than the burned FC in Brazil and Australia, respectively. In western USA and eastern Australia, the BA per vegetation type mostly corresponded to the high number of AF and higher fire intensity. It is expected that the antecedent higher Ts and Ws and low RH resulted in a rapid loss in canopy moisture and dry-up of the herbaceous plant components such as leaves, branches and grasses [74] (see Figure A2). At the time of the fires, higher Ts (i.e., 30°C to 40°C) coupled with lower Prec. (i.e., <0.4 mm hr -1 ), and availability of continuous fuel (Figures 1 and 3) influenced the establishment of new fires, propagation and intensity of the fires in FC, SL and GL, as evidenced by the higher number and clustering of AF and fire intensity (Figure 2, Table 1), and extensive losses in vegetation cover ( Figure  5). In fact, at the three sites, moderate to extreme dry (or drought) conditions were evident ( Figure  A2), causing drier leaves. This is consistent with Nolan et al. [45], who showed that the widespread fires in eastern Australia were driven by dry fuel conditions, caused by prolonged drought, i.e., before and during the fire. In another study, Turco et al. [75] show that fire occurrence is significantly related to drought conditions. Similarly, Moreira de Araújo et al. [68] found that the higher concentration of fires in 2007 and 2010 over Brazilian Cerrado and Amazon Rainforest was related to low rainfall resulting from La Niña climatic pattern. Contrarily, Kelley et al. [67] showed that the 2019 Brazilian fires had a low meteorological influence, i.e., < 7%, and concluded that land management was the most probable driver of the increase in the burned area. Overall, the interaction of various factors, i.e., meteorological, vegetation availability and conditions and human activities influenced the ignition, flammability, combustion, spatial patterns and intensity of wildfires at the three sites. Although not studied here, the influence of topography on wildfires is considerable. For example, topographic variables such as elevation influence the temperature, while aspect influences the availability of moisture and heat, and slope influences the vegetation types, intensity and rate of spread of fires [76]. Moreover, Estes et al. [76] showed that lower slope positions have a higher probability of experiencing lower flame lengths and lower rates of spread in a backing downhill fire. Future studies should integrate more variables including meteorological, vegetation conditions and topographic variables. Moreover, there is a possibility of exclusion of smaller fires MODIS (MCD64A1, 500 m) product [19,77,78]; therefore, the BA and associated burned vegetation types, in this study, may have been underestimated considering vast climatic regions studied. In fact, a recent study, i.e., Boschetti et al. [79] found the lowest errors in Boreal Forest, Tropical and Temperate Savanna, and the highest errors in the Tropical Forest, Temperate Forest and Mediterranean biomes. Although not studied here, the influence of topography on wildfires is considerable. For example, topographic variables such as elevation influence the temperature, while aspect influences the availability of moisture and heat, and slope influences the vegetation types, intensity and rate of spread of fires [76]. Moreover, Estes et al. [76] showed that lower slope positions have a higher probability of experiencing lower flame lengths and lower rates of spread in a backing downhill fire.
Future studies should integrate more variables including meteorological, vegetation conditions and topographic variables. Moreover, there is a possibility of exclusion of smaller fires MODIS (MCD64A1, 500 m) product [19,77,78]; therefore, the BA and associated burned vegetation types, in this study, may have been underestimated considering vast climatic regions studied. In fact, a recent study, i.e., Boschetti et al. [79] found the lowest errors in Boreal Forest, Tropical and Temperate Savanna, and the highest errors in the Tropical Forest, Temperate Forest and Mediterranean biomes. The highest errors in these biomes were attributed to small and temporally non-persistent burns. Therefore, research into operational methods for characterising BA from higher resolution data such as from Landsat 8 or Sentinel-2 should be prioritised to enable representative BA estimate to support fire management and recovery efforts. Moreover, higher resolution and detailed land cover products from sensors such as Landsat-8 OLI Sentinel-2 MSI should be considered in future to characterise burned vegetation types. In the future, we will study the longer periods to characterise meteorological, fuel and environmental parameters before, during and after the wildfires as well as multifactor correlation analysis incorporating topographic (i.e., slope, aspect and elevation). Due to vast fire regimes at the sites chosen for this study, conclusions about local drivers of these wildfires cannot be made; thus, further exploration in future studies is needed.

Conclusions
This study found that AFs were densely clustered alongside the coast in the western USA and eastern Australia, while they were widely spread across the Brazilian landmass. The western USA had relatively lower BA and short-term burning. In contrast, Brazilian and Australian AF were denser, the BA per day was relatively higher and the fire scars were widespread, attributable to various vegetation distributions, types and drought conditions as well as increased human activities. This is consistent with Nolan et al. [45], and Lizundia-Loiola et al. [33]. Moreover, the highest FRP, at all sites, emanated from areas characterised mainly by FC, and SL in the western USA and eastern Australia and FC, SL, GL and CL in Brazil. In fact, the study found that the FRP over FC was significantly different (p < 2.2 × 10 −16 ) from other vegetation cover types in all study sites. This finding is consistent with Vadrevu et al. [66]. Moreover, this finding implies a more significant loss of carbon stocks and forest ecosystems, which may take a more extended period to recover or result in permanent deforestation and loss of habits especially in fire-sensitive biomes such as Amazon and Atlantic Forests. This necessitates habitat-specific wildfire disaster management and efficient firefighting resource allocation, to safeguard the rare and endangered fauna and flora and ecosystem services. Consequently, the BA per vegetation cover type was the highest in Brazil, owing to its vast area, widespread fires and high number of AF, driven by human interference and influenced by changing policies that favours agricultural expansion. Overall, this study contributes to a better understanding of wildfire behaviour and underlying environmental and meteorological causal factors in various regions, towards better habitat-specific wildfire disaster management plans, for better disaster preparedness, mitigation and recovery.