Next Article in Journal
Combining LiDAR, SAR, and DEM Data for Estimating Understory Terrain Using Machine Learning-Based Methods
Previous Article in Journal
Three Ophiostomatalean Fungi Associated with Bark Beetles from Pinus thunbergii Infested by Bursaphelenchus xylophilus in Laoshan Mountain (Shandong, China)
Previous Article in Special Issue
Research Trends in Wildland Fire Prediction Amidst Climate Change: A Comprehensive Bibliometric Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forest Fire Severity and Koala Habitat Recovery Assessment Using Pre- and Post-Burn Multitemporal Sentinel-2 Msi Data

by
Derek Campbell Johnson
*,
Sanjeev Kumar Srivastava
and
Alison Shapcott
School of Science Technology and Engineering, Centre for Bioinnovation, University of the Sunshine Coast, Maroochydore, QLD 4558, Australia
*
Author to whom correspondence should be addressed.
Forests 2024, 15(11), 1991; https://doi.org/10.3390/f15111991
Submission received: 19 September 2024 / Revised: 20 October 2024 / Accepted: 31 October 2024 / Published: 11 November 2024

Abstract

:
Habitat loss due to wildfire is an increasing problem internationally for threatened animal species, particularly tree-dependent and arboreal animals. The koala (Phascolartos cinereus) is endangered in most of its range, and large areas of forest were burnt by widespread wildfires in Australia in 2019/2020, mostly areas dominated by eucalypts, which provide koala habitats. We studied the impact of fire and three subsequent years of recovery on a property in South-East Queensland, Australia. A classified Differenced Normalised Burn Ratio (dNBR) calculated from pre- and post-burn Sentinel-2 scenes encompassing the local study area was used to assess regional impact of fire on koala-habitat forest types. The geometrically structured composite burn index (GeoCBI), a field-based assessment, was used to classify fire severity impact. To detect lower levels of forest recovery, a manual classification of the multitemporal dNBR was used, enabling the direct comparison of images between recovery years. In our regional study area, the most suitable koala habitat occupied only about 2%, and about 10% of that was burnt by wildfire. From the five koala habitat forest types studied, one upland type was burnt more severely and extensively than the others but recovered vigorously after the first year, reaching the same extent of recovery as the other forest types. The two alluvial forest types showed a negligible fire impact, likely due to their sheltered locations. In the second year, all the impacted forest types studied showed further, almost equal, recovery. In the third year of recovery, there was almost no detectable change and therefore no more notable vegetative growth. Our field data revealed that the dNBR can probably only measure the general vegetation present and not tree recovery via epicormic shooting and coppicing. Eucalypt foliage growth is a critical resource for the koala, so field verification seems necessary unless more-accurate remote sensing methods such as hyperspectral imagery can be implemented.

Graphical Abstract

1. Introduction

Habitat loss due to anthropogenic impacts and climate-induced effects including wildfires is becoming an increasingly serious problem internationally for threatened animal species [1,2]. Tree-dependent and arboreal animals are at risk of displacement. Examples include the African Forest Elephant (Loxodonta cyclotis) [3]; the Woodland Caribou in Canada [2], the Mount Graham Red Squirrel in the USA (Tamiasciurus hudsonicus grahamensis) [4] if fires increase in scale, the Giant Noctule in Europe (Nyctalus lasiopterus) [5], and in Australia, the Koala (Phascolartos cinereus) [6]. Koalas are a threatened species mainly because of habitat loss from clearing and fire, but also because of drought, predators and disease [3]. Although koalas are arboreal, they are not as mobile as some other arboreal animals when they need to escape from a wildfire. Birds are able to fly to safety, and some Australian arboreal mammals can glide between trees, such as the Greater Glider (Petauroides volans), while some mammals can move between tree canopies, such as some monkeys [7]. The koala, however, usually has to climb down from a tree and move across the ground to move to other trees [8]. Dense shrub layers and undergrowth can make traversal difficult and expose koalas to predators, particularly dogs. This probably makes them more vulnerable to fires than many other arboreal species.
Large areas of Australian forests were burnt by widespread wildfires in Australia in 2019/2020 [9,10,11]. These fires were preceded by drought [12], which likely increased the flammability of these forests [13]. Most of these forests and woodlands are dominated by eucalyptus trees, which provide habitats for the koala (Phascolarctos cinereus); hence, these fires potentially impacted large areas of their habitats, a phenomenon that occurs across much of Eastern Australia [14]. Shabani et al. [15] found that nearly 40% of the total koala habitat in Australia has a high or very-high fire susceptibility rating and predicted that this figure will increase to nearly 45% by 2070 due to climate change. Eucalypt forests can burn at frequent intervals (for example, annually) depending on climatic conditions and anthropogenic activities [16]. If these fires increase in intensity and frequency, then the koala will be at increased risk.
The recovery responses of eucalypt communities have the potential to significantly affect the suitability of habitats for koala populations [17]. Many eucalypts can recover by means of epicormic shooting on the stems and/or coppicing from the base of the trunk, which are both parts that are edible for koalas. If the fire intensity is too severe, then this will result in top kill [16], but the root system can survive and subsequently be coppiced. Although the most widespread severe fires of 2019/2020 in Australia occurred in the states of New South Wales (NSW) and Western Australia, severe fires were also widespread in Victoria and Queensland [13]. The koala habitat in South-East Queensland (SEQ) was therefore also impacted by the fire in 2019/20 [18], and the recovery of that habitat was found to vary with fire severity [18]. Koalas are known to utilise a variety of different eucalypt-dominated forest types in Queensland, Australia; these fine-scale differences are mapped as Regional Ecosystems (REs; Queensland Department of Environment and Science) [19]. Each RE is defined by a descriptive code consisting of three numbers, which denote biogeographical region, broad geomorphological unit (land zone) [20], and unique vegetation type, respectively (e.g., 12.8.16; Table 1).The post-fire recovery of eucalypt forests has been found to vary with forest type, and this is influenced by the varying post-fire responses of the constituent tree species [18,20,21,22]. The efficacy of these responses will affect habitat availability for koalas. Koala habitat mapping [14] has been used to identify the primary koala habitat REs (Supplementary Table S1).
The mapping of koala habitats while also incorporating fire resilience and vulnerability is becoming increasingly necessary due to climate change and the increasing frequency and severity of wildfires [17], which may destroy koala habitats temporarily or permanently [23,24]. Knowledge of where these koala habitat areas with fire vulnerability and fire resilience are located will assist in the planning of future habitat management.
The extent and severity of fires can be mapped using remote sensing data acquired from Sentinel-2 satellites [25]. To investigate the extent of and damage to koala habitats, the spectral bands of satellite imagery can be manipulated to identify burnt areas using indices such as the Normalised Burn Ratio and calibrated with field data such as the Composite Burn Index [26]. Sentinel-2 MSI satellite data [27] are useful for analysing localised fire impacts due to their relatively high spatial (10–20 m) and temporal (5 days) resolutions [28]. By extrapolating findings from a local study area (LSA), regional assessments can be made. Regional vegetation fire studies using Sentinel-2 data have been successful worldwide [29,30,31]. A 230 Ha koala habitat study area [18] in SEQ was extensively burnt by a severe wildfire in mid-November 2019. Regional Ecosystem (RE) mapping (Queensland Department of Environment and Science) [32] provides information on described forest types and other vegetation types in Queensland [19] and can therefore be used in conjunction with Sentinel-2 data to expand a local fire study to a regional scale. Each Sentinel-2 scene covers about 12,000 km2, so one scene covering a great expanse of SEQ can be classified according to the results from a local fire study, and this can in turn be applied to each of the relevant koala habitat REs.
This study identifies the distribution and extent of koala habitats in a large part of SEQ and quantifies how much of each koala habitat forest type was burnt, how badly each was burnt, and how well each one recovered. This study was conducted to identify an effective way of mapping wildfires’ impacts on koala habitats and subsequent recovery at a regional scale. This study asks the following question: How do koala habitat forest types respond to different fire severities at a regional scale? The following specific questions were addressed: (a) How much of each koala habitat forest type was burnt on a regional scale, and how severely were they burnt? (b) Where are the fire-free refugia for koalas in this study area? (c) Do some koala habitat forest types take longer to recover than others in this region?

2. Methods

2.1. Local Study Area and Forest Types

The local on-ground survey area (LSA) [18] was a 230 Ha property about 10 km south-east of Crows Nest, Queensland, Australia, located within SEQ (Figure 1) and the SEQ bioregion. The LSA has several native eucalypt forest types, comprising remnant vegetation on steep hills, and several small creeks [18]. A survey of the LSA found that the types of forests, which were all eucalypt forests (Supplementary Figure S1), were generally consistent with RE classifications [19]. Koala habitat mapping is available for the Australian distribution of the koala (Koala Habitat Atlas; KHA) [14] and based on the constituent tree species of the vegetation-mapping units produced by state government authorities. It is classified according to the known tree preferences of koalas. In Queensland, the KHA is based on RE mapping [32] at 1:25,000 scale for SEQ. It includes a variety of vegetation types classified by floristics and structure, including forests, woodlands, wetlands, rainforest, etc. Three REs in SEQ have been identified as ‘Primary’ koala habitats by the Australian Koala Foundation [14]. Two ‘Primary’ Koala habitats are found in the LSA; they are mnemonically coded as GBS (mountain blue gum and stringybark on granite; RE 12.12.23) and BFA (blue gum on alluvial flats) (Table 1). There are also two REs present in the LSA classified as ‘Secondary A’ [12], also constituting valuable koala habitats. These are IPR (ironbark and pink bloodwood on granite ridges; RE 12.12.12) and IBM (ironbark and blue gum on microgranite; RE 12.8.16) (Table 1). The third ‘Secondary A’ habitat is very limited in extent and devoid of koala food trees (as defined by AKF) [12], and this RE was excluded from this study. The sixth RE present in the LSA is a wetland and was also excluded.
By adding the third bioregional primary koala habitat RE to this study, which was not observed within the LSA, the complete coverage of all primary koala habitats within a large part of the SEQ bioregion was achieved. This habitat was BCF (blue gum on alluvial flats closer to the coast and not within the LSA: RE 12.3.11) (Table 1). The REs in the LSA were assessed and mapped from field survey sites (Figure 1), and then the RE mapping [32] was used to investigate regional impacts for these REs. Only BCF, the RE of interest outside the LSA, lacked field data.
The LSA was assessed using 88 plotless sample sites (Figure 1) in a stratified random arrangement to include all terrain types and all fire severity and recovery types [18]. The LSA was first surveyed in early 2020 to identify the REs within it, validate and revise existing RE mapping via QDES [32], and assess fire impact severity. Revised RE mapping within the LSA was performed at a finer scale of approximately 1:5000 so that sample sites could be correctly located and classified [18].
The fire in the LSA occurred in November 2019 (late spring), and the impact survey was conducted in February 2020 (soon after in late summer). The plotless sites were revisited annually to monitor recovery in late 2020 and then in 2021 and 2022. Fire recovery monitoring for our study was limited to three years, set as the minimum anticipated time to observe recovery but also to allow the acquisition of adequate data to yield timely research findings. It is acknowledged that further monitoring could be of additional benefit, but the results will indicate if benefits decrease beyond this time. The influence of terrain (slope and aspect) was assessed to determine the consistent impact among similar terrain classes [18]. The criteria for fire severity ratings (FSRs) [18] were assigned to a range from 0 to 5 (Table 2, Figure 2) based on post-fire observations in the LSA [18] and similar fire studies in Australia conducted by other researchers [33,34].
A more detailed post-fire assessment was made using the GeoCBI (geometrically structured composite burn index) [25] from detailed field photographs. The GeoCBI is a modification/improvement of the Composite Burn Index [26], which is an on-ground assessment of the effect of fire on vegetation and soil and measures burn severity on each of the vegetation strata on a site. The GeoCBI adds a fraction of cover (FCOV) for each stratum and the leaf area index (LAI) for the intermediate and tall-tree strata, and a rating is applied to each site. In our study, it was not possible to record changes in the leaf area index (LAI), as pre-burn data were not available. The use of other remote-sensing-based indices instead of GeoCBI was also considered, including the Burn Area Index (BAI) or Relativised Burn Ratio (RBR). The BAI is dependent on persistence of charcoal, whose presence is more typical in Mediterranean areas [35]. The RBR was tested in the Western United States [36], but unlike our study, the plant communities tested are not eucalypt forests, so this metric was not preferred, whereas the GeoCBI has been tested successfully in Australia with eucalypts [37]. The GeoCBI uses a straightforward set of field-collected criteria, making it suitable for calibration of remotely sensed data, but the leaf area index (LAI) needs to be measured before and after a fire, so pre-fire measurement may not have been foreseen in some instances. Goodness-of-fit of the GeoCBI with remotely sensed data can be tested using regression. Extending the use of the GeoCBI to a regional scale, based on locally collected data, may have limitations due to regional environmental variation, but this also applies to other methods such as those using the BAI and RBR. In our study, a comparison of the FSR and the GeoCBI with the dNBR using regression [38] confirmed that both methods correlated with the dNBR based on 88 plotless sample sites in the LSA. For FSR, R = 0659, and for GeoCBI, R = 0.716. GeoCBI values were accordingly categorised into six classes (0–5) similar to FSR (Table 1). Within the LSA, the GeoCBI values ranged from 0.0002 to 2.8667, with a mean value of 1.777.

2.2. Regional Context

The regional study area, including the LSA, is located in the southern part of SEQ, which forms the SEQ bioregion [20]. There are currently 172 Regional Ecosystems (REs; plant communities) recognised and mapped in this bioregion with mapping [32] and descriptions [19]. The RE mapping covers the whole of Queensland and is publicly available in digital format at 1:100,000 scale statewide and at 1:25,000 scale for the SEQ bioregion. At the national scale, RE classes are generalised into Broad Vegetation Groups at 1:1,000,000 scale [39,40,41], but this scale was too coarse for our study.
The five REs of interest (Table 1) were overlaid on five Sentinel-2 scenes covering a significant portion of the southern half of the SEQ bioregion. Each scene, or tile, covers an area of slightly greater than 100 km × 100 km (12,000 km2; scenes overlap). The coverage excluded the state of NSW to the south, for which a different vegetation-mapping system is used; the Pacific Ocean to the east; and the Brigalow Belt bioregion to the west (markedly different vegetation types) [19,32]. The regional area of interest was reduced to a single Sentinel-2 scene, which included the LSA, and all five of the REs of interest (Figure 1). This area constitutes almost 20% of the SEQ bioregion, which covers approximately 66,000 km2 [20]. The LSA is located towards the north-west corner of the Sentinel-2 scene, and adjacent scenes were not used because the Brigalow Belt bioregion is located about 20 km west, and there were a limited number of REs of interest to the north. The LSA was far away enough from the Brigalow Belt bioregion to eliminate edge effects, and the five REs of interest were well distributed across the selected Sentinel-2 scene. A map of the five Sentinel-2 scene extents, and the REs of interest in relation to these extents, is shown in Supplementary Figure S2. Within the regional area covered by the Sentinel-2 scene, RE mapping was clipped into each RE of interest by using ArcGIS v10.8 [42]. RE polygons include both single (pure, 100%) RE types and mosaics, which are mixed with other RE types as a proportion (e.g., RE type 1, 70%; RE type 2, 30%). The satellite imagery (to be processed) was considered reliable for 100% (pure) polygons but less reliable as the proportion of RE types mixed. Excluding all mosaic polygons would have significantly under-represented each RE of interest, so a minimum of 70% dominance was used [41,43,44]. This allowed up to 30% of confounding RE types to occur within mosaics, which could reduce the reliability of the dNBR but generally increase the inclusion of the RE of interest.

2.3. Data Sets and Sources

Datasets used in this study (other than field sample sites) were sourced from public web servers (Table 3). The RE mapping [32] was obtained from a geodatabase polygon file (ArcGIS v10.8) [42] with several attributes, and those used were RE and mosaic percentage of each RE (up to five REs per polygon). Field data included the GeoCBI, Fire Severity Ratings (FSRs), and the RE mapping revision (forest types), and all were stored as shapefiles (ArcGIS v10.8) [42]. Sentinel-2 MSI data [27] (Table 4) was selected for analysis due to their relatively high temporal and spatial resolutions compared to other earth observation products. These datasets were analysis-ready data from Digital Earth Australia (DEA) [45]. Optical surface reflectance data were standardised using pre-processing models already processed by DEA to correct inconsistencies in upper-atmospheric reflectance values, including aerosol optical thickness and ozone, atmospheric correction coefficients, MODTRAN, and FMASK [45]. Cloud-free images were selected. Inconsistencies in terrain were already corrected by DEA using nadir-corrected bidirectional-reflectance-distribution-function-adjusted reflectance (NBAR) with additional terrain illumination correction (NBART) [45]. The quality assurance processes are detailed by FrontierSI [46]. These corrections were applied uniformly across all temporal images for a single scene. Queensland Fire Scar Mapping [47] is a geodatabase [42] of data collected based on methods developed by Goodwin and Collett [48] and Hardtke et al. [49] (Supplementary Table S2). Fire history mapping [50] was in shapefile format [42], based on reports of fire timing (start to finish) and extent, and used to confirm which fires within the regional study area were burning in between the pre- and post-fire Sentinel-2 dates in November 2019. Regional koala records for the last 10 years were obtained from the Atlas of Living Australia [51].

2.4. Remote Sensing Data Analysis

The NIR and SWIR2 spectral bands from Sentinel-2 (Table 5) were used for the calculation of the Normalised Burn Ratio (NBR), the differenced NBR (dNBR), and the Normalised Difference Vegetation Index (NDVI). The imagery used was a series of scenes from 2017 through to late 2022 (used to provide years before fire), then 10 days before fire in mid-November 2019, then 20 days post-fire, and then one image per year after 2019. Sentinel-2 datasets free of clouds (with <10% cloud cover, and no clouds over the LSA) were selected to match dates of field sampling as closely as possible (Table 4). A cloud-free image taken at a date as close to before the fire as possible was selected, and the timing of images in the years preceding the fire was intended to be as close as possible to an exact year, with some variation due to cloud-free-image availability. Following the fire, a cloud-free image taken as close to immediately after the fire as possible was selected, and then recovery-year images as close to the field-work sampling dates as possible were selected (these images were cloud-free).
Whilst vegetation recovery can be assessed with the NDVI and the Differenced NDVI (dNDVI), the dNBR was preferred for this study [52,53,54] to maintain consistency when comparing years of recovery to the initial fire impact, which was measured with dNBR [55]. We tested our field measurements of fire impact and recovery [18] against the Sentinel-2 data (dNBR and NDVI) in order to extrapolate our findings from the LSA (230 Ha) to a regional scale (Table 1). In our study, the NDVI was preferred over the dNDVI as it was used as a condition assessment tool, as an independent snapshot in time, and not a measurement of change. We used the dNBR to measure change.
The NDVI for each year of Sentinel-2 imagery was calculated using ArcGIS v10.8 [42]. The NDVI is a normalised ratio of near-infrared (NIR) to red:
NDVI = (NIR − Red)/(NIR + Red)
Sentinel-2 band for NIR = band 8a, and for Red = band 4 (Table 5). Narrow NIR from band 8a (20 m resolution) was preferred to the broader spectrum infrared from band 8 (10 m resolution). A narrower spectral resolution for NIR was considered more important for this analysis than the spatial resolution, the latter of which was adequate for a regional-scale study. The NDVI was only used for general comparison of vegetative growth for each year, including years prior to the fire as a reference.
The differenced NBR (dNBR) calculation for each year of Sentinel-2 imagery was performed using ArcGIS v10.8 [42]. The NBR [56] is a normalised ratio of near infrared (NIR) to shortwave infrared (SWIR):
NBR = (NIR − SWIR)/(NIR + SWIR)
The Sentinel-2 bands for NIR = band 8a, and for SWIR2 = band 12 (Table 5). Narrow NIR from band 8a (20 m resolution) was preferred to the broader-spectrum infrared from band 8 (10 m resolution). A narrower spectral resolution for NIR was considered more important for this analysis than the spatial resolution, the latter of which was adequate for a regional-scale study. This study was conducted at 20 m resolution so band 8a did not decrease overall resolution in this study. The differenced NBR (dNBR) is the difference between two NBRs calculated at different times, before and after a fire:
dNBR = pre-fire NBR − post-fire NBR
Resolution of the NBR was retained at 20 m, as this was sufficiently precise for regional analysis, so no resampling to 10 m resolution was necessary. This also significantly improved computation speed for repetitive analyses of entire Sentinel-2 scenes and eliminated system crashes. A resolution of 20 m is approximately equivalent to the nominal 1:25,000 scale of the base mapping of REs [32] in SEQ, where line-work precision is no finer than 25 m. The NDVI was also processed at 20 m resolution due to the regional scale of the analysis. Differentiation of small vegetation patches below 20 m resolution was not considered desirable at regional scale due to the potentially fine-scale stochastic nature of burn patterns within each patch of vegetation [57], including influence of vehicle tracks, fence lines, rocky outcrops, etc.
The dNBR was classified for fire impact and fire recovery using ArcGIS v10.8 [42]: For fire impact, the GeoCBI was used based on the 88 plotless sites in the local study area (Supplementary Figures S3 and S4). However, the GeoCBI was inadequately sensitive to differentiate very-low levels of fire severity (or relatively healthy vegetation) throughout the region. A different method of classification was therefore required to fully assess regional recovery (GeoCBI was only used for the dNBR of the year of impact, 2019). For the years of fire recovery (2020, 2021, 2022) and also impact (2019) for comparison, a system of manual breaks of the dNBR was used (Table 6). The dNBR, GeoCBI, and FSR data were explored to determine their distribution, and various categorisation methods were investigated. Finally, to draw a comparison of dNBR with GeoCBI and FSR on a temporal scale (corresponding to the year of the 2019 fire and then three subsequent years of recovery), we categorised the dNBR into 22 incremental classes (manual breaks; ArcGIS v10.8) [42] (Table 6). This was defined as a Comparable Incremental Scale Reference (CISR). Other methods of classification included in ArcGIS v10.8 [42], which were not suitable, included natural breaks (Jenks), equal interval, quantile, geometric interval, and standard deviation (ArcGIS v10.8, Supplementary Figure S5) [42]. They were unsuitable for three reasons: Firstly, the classes created did not have a common break point between years, based on an dNBR of zero, so years could not be directly compared. Scales of dNBR would also have been different, so high dNBR percentiles from one year might not be the same as high dNBR percentiles in another year. Secondly, all of those automatic classifications were unsupervised, so none of the classes could be assigned meaningful fire severities or recovery statuses (whereas GeoCBI can). Therefore, multi-temporal comparison of categories was performed based on CISR, as opposed to the use of GeoCBI for the fire impact alone. The GeoCBI intervals were different to those of the CISR, but its six FSR classes (0–5) could be aligned with those increments (Supplementary Table S3). This therefore provided some impact description to the burnt section of the CISR. Thirdly, although some of the other methods (natural breaks, quantile and standard deviation) could be used to divide the distribution of dNBR values into meaningful classes for the year of the fire (2019), those break values could not be meaningfully transferred to later years of recovery (2020, 2021, 2022) because the distribution of dNBR values in most cases moved outside the bounds of the classes set for 2019 (note the graphical stepwise explanation in Supplementary Figure S6). Other methods were also considered, such as k-means clustering and machine learning-based approaches, but the CISR method produced the most easily interpretable and intuitive graphical results and was able to be processed within ArcGIS v10.8 [42], along with all other data in our study.
False positives from the dNBR were not a problem in this study because the RE mapping of koala habitat forest types was unlikely to contain significant wet areas such as wetlands or water bodies. The main source of false positives was agricultural areas due to changes in cropping between pre-fire and post-fire Sentinel-2 images, which generated a non-zero dNBR value. For example, a harvested crop paddock would generate a positive value, and this was observed abundantly in the Lockyer Valley, an extensive agricultural area between Toowoomba and Brisbane (Supplementary Figure S7). An agricultural mask [58] coincided with most, or possibly all, of these false positives. Nevertheless, the mapped REs of interest do not coincide with agricultural or cultivated areas.

2.5. GeoCBI Classification

The GeoCBI [25] was assigned to each of the 88 plotless sites within the LSA and then classified to enable comparison of different burn severities for each forest type. The existing classification tool from the field proforma [25] is a Burn Severity Scale ranging from 0 to 3, but it has seven sub-classes within it. In order to make this scale compatible with previous research in this LSA conducted by Johnson and Shapcott [18] (for which a Fire Severity Rating was used; Table 2), two of the intermediate sub-classes were combined to form a 0–5 scale with generalised descriptors: 0 = none, 1 = low, 2 = moderate low, 3 = moderate high, 4 = high, and 5 = severe. This GeoCBI could then be applied regionally to the dNBR.

2.6. Data Validation

A digital elevation model [59] was used to test for terrain as a confounding factor for this study by comparing dNBR values from the Sentinel-2 data with elevation, slope, and aspect values (terrain) using ArcGIS v10.8 [42] and linear regression with SPSS [38]. Weather conditions in the weeks leading up to the fire in mid-November 2019 and for the remainder of that month were consistently dry, with minimal rainfall [10,18]. Regional Ecosystem mapping (QDES 2023a) was ground-truthed in key regional locations for each RE that occurred within the LSA. In some areas that could not be visited, Google Maps Street View [60] was used. Within the LSA, GeoCBI classes from the 88 plotless sites correlated with the dNBR, with a Pearson value of 0.716 (c.f. [37,61]). The GeoCBI-classified dNBR was therefore accepted as a reliable indicator of wildfire severity for this study [61]. Results of burnt area mapping and total burnt areas were compared with Queensland Fire Scar Mapping [47].

2.7. Analysis

2.7.1. How Much of Each Koala Habitat Forest Type Was Burnt on a Regional Scale, and How Severely Was Each Type Burnt?

Fire impact was assessed by comparing the areas of fire severity classes (GeoCBI) within each RE. Total area calculations were obtained from raster analysis (ArcGIS v10.8) [42]. Total areas of each fire severity for each RE were also converted to percentages. These results provided information on how many koala habitats were affected by fire and how severe the impact was.
We compared total areas of each fire severity (GeoCBI classes 0–5) within each RE using the Kruskal–Wallis test in SPSS [38]. We also compared the highest fire severity (GeoCBI classes 4 and 5 combined) between REs. No statistical comparison was needed for subsequent years (of recovery) using the 2019 GeoCBI because it was not able to differentiate between areas of low fire severity at the regional scale.

2.7.2. Where Are the Fire-Free Refugia for Koalas in This Study Area?

Koala habitat areas suitable for refuge from fire were identified using the following criteria: (1) the definition of koala habitat forest type as specified in our study (i.e., REs were GBS, IPR, IBM, BFA, and BCF); (2) fire severity class 0 = unburnt; (3) proximity to burnt areas of moderate, high, and severe ratings—classes 2, 3, 4, and 5; (4) proximity to recorded koala sightings. Only koala habitat forest type areas with a koala sighting recorded less than 10 years ago were included to increase the likelihood of a koala being present in these mapped forest types.
Buffers around burnt areas ranged from being directly contiguous with the burnt area to a maximum distance equivalent to the average area of a koala’s maximum home range. This was initially set at a buffer of 700 m based on mean home-range area of 39.5 Ha [62] in the event of a koala relocating to the refuge permanently or for the long term. However, home ranges were not necessarily round; they could be elongated [63] to suit resources. The distance (length) was therefore increased to 4 km (corresponding to a rectangle with a width of 100 m), which is within known travel distances of koalas [62]. This area could represent a gully, for example. All koala sighting records were buffered by the average radius of this home range area and therefore set to 350 m. The distances selected were intended to be practical minima and could be extended for various reasons, such as connectivity, or locally reduced for reasons such as land-use tenure or drought. Polygons of each criterion were buffered and then converted to raster in ArcGIS v10.8 [42]. All unburnt (fire severity class 0) areas that were overlapped by burnt area buffer were selected as refuges, provided that at least some koala sighting buffer area overlapped with the burnt area buffer.

2.7.3. Do Some Koala Habitat Forest Types Take Longer to Recover than Others in This Region?

Fire impacts based on the classified GeoCBI from the LSA were initially intended to also be used as the reference for fire recovery across the regional study area, but we found that it was not an effective regional metric when based on data from the LSA, which was relatively severely burnt. The proportional amounts of fire impact at different levels of severity in the year of the fire (2019) in the region were easily distinguished, but the subsequent years of recovery (2020, 2021, 2022) showed that almost all of the region had suffered no fire impact or had completely recovered or was in a state of healthy growth (Table 6). This finding was similar to Queensland Fire Scar Mapping [47] and indicated that a different scaling of fire impact and recovery (CISR, based on the dNBR) needed to be developed for regional study of fire recovery to account for where fires were minimal or absent. We calculated the total areas of each fire severity level within each forest type (RE) and the total percentages of each forest type using ArcGIS v10.8 [42]. The classified dNBR data were analysed with the Bivariate Correlation package in SPSS v28.0.1.1 [38]. We used Spearman’s rank correlation to test between slower and faster recovery.
As an overall assessment of original pre-fire condition for all vegetation and the impact of fire, followed by overall recovery, the NDVI was calculated for two years prior to the fire in 2019, then for 2019, and then for three years after the fire, across 88 plotless sites in the LSA [18]. This gave us a baseline from which to determine if vegetation had apparently returned to its original condition after the fire.

2.7.4. Reliability Assessment of the dNBR and NDVI for Forest Fire Recovery

Reliability of the dNBR for forest fire recovery was tested post hoc because the dNBR is a well-recognised method of assessing fire impact [54,55]. The dNBR was preferred over the NDVI for this, although both methods are suitable for monitoring recovery [64]. There are various ways of improving the reliability of the NBR (and hence the dNBR), such as reducing false positives from water bodies and clouds [65], but distinguishing between different types of vegetative growth (e.g., epicormic shooting and coppicing) may be more difficult [66]. In the case regarding vegetation recovering from fire, the dNBR can quantify overall vegetative recovery [52,53]. There are other indices for burned-area mapping, but their evaluation requires more detailed field data. For example, the Burned Area Index for Sentinel-2 (BAIS2) [67] requires leaf area index data, which we could not record in our study. Our study focusses on koala habitats and therefore requires the assessment of specific types of tree regrowth responses, namely, epicormic shooting, coppicing, and apparent tree mortality. These responses were recorded at 88 plotless sites within the LSA over three years [18] as a percentage of trees and then compared using the dNBR and the NDVI to see if these responses could be identified and quantified using regression in SPSS [38].

3. Results

3.1. How Much of Each Koala Habitat Forest Type Was Burnt on a Regional Scale, and How Severe Was the Burning?

The total burnt area of the five koala habitat REs studied in the regional area, including fire severity classes, 1–5, was 2400 Ha (Table 7), which is 10% of the total area of these REs (wherein those REs occupy >70% of the forest). Some additional area would have been burnt in REs comprising less than 70%, but this is likely to have been relatively minimal due to the sparsity of such regions. The five koala habitat REs only occupied 2% of the area of this region. The remainder of the study area was composed of other REs of probably less habitat value for koalas, non-remnant (cleared) areas, and water bodies.
The areas of each koala habitat forest type in the regional study area ranged from approximately 3000 Ha to 7000 Ha (Table 7), with GBS corresponding to 5744 Ha (mid-range), but it had approximately 10 times more area that was severely burnt in 2019 (class 5; 566 Ha; Table 8; Figure 3) than the other forest types combined (54.7 Ha). It also had almost four times more area badly burnt in 2019 (class 4; 332 Ha) than the other forest types combined (89 Ha). The fire impacts on koala habitat forest types in 2019 in the regional study area were greatest for the forest type (RE) GBS, followed by the REs IBM and IPR (Table 7; Table 1; Figure 3). There was 2910 Ha of IPR within the regional study area, of which 161 Ha was burnt, and there was 5848 Ha of IBM within the regional study area, of which 374 Ha was burnt. In contrast, for the alluvial areas, there was 6759 Ha of BFA within the regional study area, of which only 11 Ha was burnt, and there was 3199 Ha of BCF within the regional study area, of which only 1 Ha was burnt. For GBS, at least 15% (898 Ha) experienced high-to-severe fire impacts (Table 8), with approximately 30% of the areas of GBS burnt to some degree (1855 Ha). All other REs studied (IPR, IBM, BFA, and BCF) were burnt less (only 7% burnt combined, totalling 550 Ha) and suffered a minimal impact from high-to-severe fire intensities (3%, 140 Ha). The overall impacts on the two alluvial REs BFA and BCF were negligible. In subsequent years (2020, 2021, 2022), almost nothing was classified as burnt or unrecovered based on the GeoCBI (Table 8, Figure 3).

3.2. Where Are the Fire-Free Refugia for Koalas in This Study Area?

The fire-free refugia for koalas in the regional study area were found to be widespread, with the largest areas being situated around the LSA near Crows Nest (Regional: Figure 4; Inset: Figure 5) and on the Main Range, a mountain range between Allora and Boonah. There were several smaller areas throughout the region based on smaller fires identified.

3.3. General Recovery Trend

The NDVI measurements indicated that vegetation cover returned to previous extents during the recovery phase after the fires (Figure 6). A significant drop in NDVI values in 2019 was confirmed (p < 0.001, Kruskal–Wallis H = 330.277).
The mean dNBR values for each forest type studied were all positive (>0) in the year of the fire (2019) and then in all cases became increasingly negative in progressive years of recovery (2020, 2021, 2022; Figure 7). The results for all the forest types are similar, except for GBS, which experienced the greatest fire impact (dNBR > 0.2), and BCF, which showed the least apparent recovery from fire, with dNBR values for each recovery year of about −0.3, whereas the other forest types showed dNBR values of about −0.4.

3.4. Do Some Koala Habitat Forest Types Take Longer to Recover than Others in This Region?

The distribution of areas across the 0.1 dNBR increments of fire severity (Figure 8) shows modes (highest scores) for each RE, and in 2019, they mostly have a slightly positive dNBR value (zero is neutral; Table 6). Although GBS has the lowest mode (see Supplementary Tables S4–S7 of results that support the graphical figures), it is bimodal and has a second peak at class 15, which is a higher fire severity.
In contrast to the year in which the fire had an impact, i.e., 2019, for the first year of recovery (2020) and both subsequent years, all REs were found to have a unimodal distribution of growth response (Figure 9, Supplementary Figures S8 and S9). All REs showed a slightly negative peak response (zero is neutral; Table 6), indicating most areas were not burnt.
A comparison was also made of the yearly changes in the total areas for each 0.1 increment in dNBR (CISR) for each forest type (RE) across the regional study area (Figure 10). For all the REs studied, the values and peak values of dNBR were highest in the year of the fire (2019) and became progressively lower in the subsequent years of recovery. The forest type GBS was the only RE studied that had a bimodal distribution of high fire severity, with a higher fire severity impact overall. All forest-type REs then experienced a reduction in their dNBR values in the first year of recovery (2020; e.g., IPR, Figure 10) and then a further reduction in dNBR in the second and third years (2021 and 2022). For all forest types studied, there was no further reduction in dNBR values from 2021 to 2022, indicating that minimal change occurred from 2021 to 2022 (Supplementary Figures S10–S12). This correlation of dNBR distributions between 2021 and 2022 was confirmed to be highly significant (p < 0.001), with Spearman values of rs = 0.953 for GBS, rs = 0.948 for IPR, rs = 0.969 for IBM, rs = 0.967 for BFA, and rs = 0.976 for BCF. All other rs values were less than 0.9, and, in most cases, they were less than 0.8, indicating that the dNBR values were similar between other years but not as closely correlated as between 2021 and 2022.
Queensland Fire Scar Mapping [47] (Supplementary Figure S13) did not differentiate between forest types (REs) or fire severities, but within our study area, it had similar overall indications of recovery (Figure 10), with approximately 4% of the regional study area (clipped to the same Sentinel-2 scene) being burnt in 2019 and less than 1% being burnt in the following years, i.e., 2020 and 2021. These were new fires and not an indication of recovery from previous fires.

3.5. Reliability Assessment of the dNBR and NDVI for Forest Fire Recovery

Tree recovery responses, namely, epicormic shooting, coppicing, and apparent tree mortality, did not correlate with the dNBR values for each of the three years of recovery (Figure 11). The dNBR values were highest in 2020, the year immediately after the fire, and then reduced in subsequent years. However, high values of field-observed tree responses (in this case, epicormic shooting) did not correlate with low dNBR values. A linear regression found no correlation between the two variables, A similar recovery response was recorded for coppicing and for trees with no signs of life, with similar non-significant regression results (Supplementary Figures S14 and S15).

4. Discussion

4.1. How Much of Each Koala Habitat Forest Type Was Burnt on a Regional Scale, and How Severe Was the Burning?

Fires impact different types of forests worldwide; for example, they affect the habitats of the woodland caribou [2] in Canada. Fire also impacts several arboreal mammals in the Brazilian Cerrado [68], which is not as fire-resilient as other vegetation types in that region, such as savannah and woodlands. In eastern Australia, different eucalypt forest types are likely to provide varying habitat quality and resources for the koala [63], so their responses to wildfire in terms of vulnerability and resilience could be important for the koala. We found that the five preferred koala habitat forest types in our study comprised only 2% of a 12,000 km2 area, which illustrates the paucity of preferred koala habitats [14] in the region. A total of 10% of all koala habitats that we studied were burnt. An even more significant proportion of GBS was burnt (32%), and 10% was severely burnt. In a study on koala habitats on the NSW North Coast, Australia, modelling by Phillips et al. [17] revealed that koalas were five times more likely to survive in areas with unburnt or partially burnt forest canopies compared to areas with fully burnt canopies. Of the five forest types we studied, GBS was far more significantly affected by fire and contains Eucalyptus tereticornis subsp. basaltica (mountain blue gum), regarded as a primary food tree for koalas [14,69]. Therefore, the fire event in 2019 was likely to have driven relatively large numbers of koalas from this habitat type to adjacent refuges, if available. At least 15% (898 Ha) of GBS experienced a high to severe fire impact, while only a combined 3% (140 Ha) of the other koala habitat REs studied (IPR, IBM, BFA, and BCF) were impacted by high-to-severe burns. Reducing the severity of large, catastrophic fires improves the survival chances of koala populations [70,71]. Other Australian arboreal mammal species are also at risk of displacement from their preferred habitats, such as the greater glider (Petauroides volans), which prefers mountain ash (Eucalyptus regnans) forests to other eucalypt type forests [72]. These mountain ash forests are at an increased risk of wildfire [73], and hence the habitat of the greater glider may also be at risk.
The presence of stringybark eucalypts, as opposed to smooth-barked eucalypts, presents a likely flammability hazard for many forest types [74,75], so a forest without stringybark species may provide a safer habitat for koalas and many other fauna species under fire-season conditions. Gill and Ashton [74] found that the relatively coarse, stringy bark of Eucalyptus obliqua caused greater heat penetration to the cambium than two other eucalypt species that had subfibrous bark and decorticating bark, respectively. Collins [75] found that eucalypts with fibrous bark were more prone to topkill. The vulnerability of fibrous-barked eucalypts to fire is uncertain, though, as Nolan et al. [21] found, the eucalypts most resistant to topkill have thick bark but low density. Further, in a South African study, Higgins et al. [76] did not find a relationship between bark thickness and topkill, but the species studied were not eucalypts. In our study, we found that the forest type IBM, which does not have stringybark as one of its major component species, was affected by fire to a much lesser extent, with less than 1% being severely burnt, even though it occupies similar terrain to GBS. There are two possible reasons that IBM burnt less: the trees in this forest type likely have lower bark flammability, and this forest type grows on finer-grained substrates like microgranite rather than granite (Queensland Department of Natural Resources, Mines and Energy) [77], so the corresponding soils are probably less well-drained, more fertile, and promote a greener and denser ground layer, which would be less flammable [78,79]. This forest type therefore serves as a less hazardous habitat than GBS, and this appears to be what we observed. It also occupies about twice as much area across the region than GBS does (33,000 Ha vs. 15,000 Ha), so this is a benefit for the koala.
In a study on an area of eucalypt forest in Tasmania, Ndalila et al. [80] found that fire severity was strongly influenced by terrain and vegetation type, which are likely to affect local arboreal fauna in the same way koalas are affected on the Australian mainland. In our study, the ridge-top forest type IPR was affected by fire to a much lesser extent than RE GBS, with only about 0.3% of it having been severely burnt, but it is subject to the inherent hazard of more limited routes of escape for koalas since fires generally tend to run uphill towards ridges [80], and access to the refuge of waterways and cooler gullies is more difficult. The optimal koala habitat therefore needs to have more of the advantageous qualities of REs to reduce fire hazards and increase the chance of escape to refugia. The upland forest types in our study (GBS, IPR, and IBM) have the inherent disadvantage for koala habitats of being located on hills or ridges and therefore vulnerable to ascending flames, as well as having flammable stringybark or a drier ground layer. Therefore, an alluvial landscape is likely to provide a less hazardous habitat for the koala. By far, the least fire-affected forest types in our study were the two on alluvium, namely, BFA and BCF, where much less than 1% of the area was burnt in each. These two forest types have all of the advantages that the others in our study lack and occupy relatively large areas (about 40,000 Ha each) across the bioregion. However, the inland alluvial BFA is under pressure from agricultural and pastoral clearing, and the coastal and sub-coastal equivalent of it, BCF, is more under threat from residential development [19]. Nevertheless, these low-lying areas of forest on floodplains and in valleys and gullies are likely to be of higher habitat value to the koala after a fire. In a study conducted in South-Eastern Australia, which is also within the distribution of the koala, Chia et al. [81] found a higher abundance of arboreal mammals in forested gullies rather than on slopes.

4.2. Where Are the Fire-Free Refugia for Koalas in Our Study Area?

Refuges allow immediate survival for animals escaping from fires, providing habitats for individuals and populations post-fire, and enable long-term population recovery [82]. Koalas have been documented to recolonise nearby burnt areas within months [62]. However, severe fires can limit the amount of available refuge to more sheltered areas [83]. We confirmed this with our observations of GBS, which occurs on drier upland areas. Of the five forest types we studied, GBS occupies most of these upland areas, and about 67% of this RE remained unburnt. Although this constitutes a relatively large proportion of the overall RE area, the other REs studied were almost completely unburnt. Because of the overall negligible impacts of fire on the two lowland alluvial forest types BFA and BCF in our study, these present themselves as more suitable areas of post-fire refuge. Other areas of refuge, as identified by our mapping (Figure 4 and Figure 5), include some of the other REs, but the alluvial REs appear to be the most suitable. However, these alluvial areas are not always near the burnt upland areas, so parts of these unburnt upland areas may function as the most immediate refuge [84].
In a study on Leadbeater’s Possum (Gymnobelideus leadbeateri) in South-Eastern Australia, Durkin et al. [85] found that population recovery depended mainly on in situ survival and not external recolonisation. A study in Victoria conducted by Berry et al. [86] found that refuges of unburnt peninsulas extending into burnt areas helped some species to persist within extensively burnt landscapes (including the mountain brushtail possum, Trichosurus cunninghami) but not others (such as the greater glider, Petauroides volans). This suggests that the home ranges of individual species are important when responding to a wildfire, and the koala is fortunate to have a home range [62] that is relatively large compared to the ranges of other arboreal mammals such as the greater glider [86]. Koalas’ home ranges are variable and have been observed to account for as much as approximately 14 Ha [63] of remnant forest near Kenilworth, Queensland, Australia, or as much as approximately 45 Ha [87] of an urban area in Lismore, NSW, Australia. Koala home ranges can be even larger, as found by Matthews et al. [62] in a remnant forest in Port Stephens, NSW, Australia, with up to approximately 70 Ha for males and up to approximately 34 Ha for females. This variation is probably due to overall habitat quality and other environmental influences. Most 4 km buffers around the burn areas in our study had koala records [51], so it was assumed that most of the refuge areas within those buffer zones were useful as koala refuges, as they are within the average koala’s home range [62]. We found that there were no cases where a koala would have had to travel more than 4 km to escape from a burnt area (Figure 4 and Figure 5), and in most cases where koalas had been recorded [51], the distance was much less. This is fortunate for the koala, but larger, more extensive wildfires may create unfeasible travel distances if climatic conditions are unfavourable.
The mapping in our study is intended to assist with the prioritisation of areas of refuge such as that established in the work carried out by Wilson et al. [88]. Landscape characteristics also need to be taken into account [89,90,91], such as cleared-edge effects and available water. Heterogeneity of burn patterns within the site may also be a factor in refuge selection [92], so in some cases in our study, adjacent refugia were possibly not as important as in situ refugia. Burnt areas can also be used for shelter, as has been observed by Matthews et al. [93]. Tree utilisation can be influenced by factors other than food value, such as security [17] and height [63]. Conservation area planning accounting for these considerations would be best supported with tools like Marxan software [94,95] and can take current koala populations into account when prioritising refuge areas.

4.3. Do Some Koala Habitat Forest Types Take Longer to Recover than Others in This Region?

Different forest types can be expected to respond differently to fire, both in terms of vulnerability and resilience [18,22]. This difference in response can result in different recovery times for each particular forest type [96], and this is of significance for fauna habitats, both resident and transitory. In our study, the forest type (RE) GBS experienced the most severe fire impact in 2019 and provided a relatively poor habitat for koalas in the first year of recovery. In that time, koalas needed to find areas of refuge in other REs until the GBS recovered. However, the comparatively rapid recovery of the GBS indicates that it probably became a suitable source of food again relatively rapidly and was most likely in a state of new foliar growth from epicormic shooting and coppicing [18]. This rapid recovery was also observed in studies by Matthews et al. [62,93] conducted in Port Stephens, NSW, Australia, and by Law et al. [24] in North-East NSW, Australia. In our study, the trees themselves were most likely unsuitable for roosting due to a lack of shade (Johnson, pers. obs.), so the use of these habitats was probably transitory, with concurrent use of adjacent habitats for roosting. This was indicated by the widespread presence of koala scat without simultaneous sightings of koalas (Johnson, pers. obs.). Studies by others on forest recovery found that koalas returned to burnt habitats within several months. Matthews et al. [93] found that within three months after a wildfire, koalas were seen among the epicormic growth. Law et al. [24] noted that koalas were temporarily extirpated wherever high fire severity dominated the landscape, but some localised recovery was evident after one year. This is consistent with our results.
Eucalypt forests can recover from fires via epicormic shooting, coppicing, and seedling production, depending on the eucalypt species [18,22] and fire severity [97]. All of our dNBR analyses found that whilst all the forest types studied had recovered noticeably in the first and second years, the amount of vegetative recovery detected did not change in the third year of recovery, which suggested that foliage production had slowed or matured. As with the results from the NDVI, the dNBR can probably only measure general vegetation conditions, including understorey growth, ground layer, and possibly weeds, in addition to tree canopies, and do not necessarily indicate tree recovery via epicormic shooting and coppicing [18]. This under-representation of the dominant tree stratum was also observed by Cansler and McKenzie [98]. The two metrics were also unable to reliably detect the proportion of epicormic shooting, coppicing, and apparently dead trees killed by the fire. This means that the NDVI and dNBR should not be solely relied upon for an indication of koala habitat recovery, as it may not be eucalypt foliage growth that is being detected, and this is the critical resource for the koala. Field verification is necessary to confirm what the increases in NDVI and dNBR actually represent, unless more-definitive measures can be used to differentiate the various plant functional traits, such as hyperspectral analysis [99]. The use of remote sensing with resolution higher than that of Sentinel-2 may allow the detection of these particular types of regrowth via patterns and/or wavelengths. There is not much more in the way of field data that could be collected to improve the current relationship between the Sentinel-2 data and field-derived indices such as the GeoCBI. Drone imagery or LiDAR may also be appropriate, but the regional scale of this study indicates that hyperspectral imagery is probably the most practical and affordable. An improved ability to discern the extent and amount of recovery via epicormic shooting and coppicing would enable habitat recovery mapping to identify more fire-resilient areas of koala habitats.

5. Conclusions

The use of GeoCBI as a dNBR classifier effectively separated fire severity impacts and clearly identified the forest type GBS as the koala habitat most severely affected by fire. The study of forest recovery required manual classification (CISR) that was sensitive enough to detect foliar regrowth. Our study found that the most suitable koala habitat in SEQ, Australia, occupies only about 2% of the region, and of that limited area, about 10% of it was burnt by wildfire in late 2019. This represents a very limited safe habitat for the koala, and refuges in wetter alluvial areas are important for its ongoing survival (e.g., the forest type REs BFA and BCF). Unfortunately, these alluvial areas are under pressure from agriculture, grazing, and clearing for development. However, the conservation of the other forest types studied is equally important. Forest types on hills were found to be more susceptible to fire, but they are good refuges from the extensive land clearing occurring on the alluvial areas. Each of these three forest types on hills has its own advantage: GBS was most severely and extensively burnt but recovered quickly; IPR on ridges would be hazardous for the koala in terms of escape from fire but accounted for the most area of all the REs studied within the bioregion; and IBM also comprises quite a large area that did not burn very extensively. We found that the dNBR can probably only measure general vegetation present and not specifically tree recovery from epicormic shooting and coppicing. Eucalypt foliage growth is the critical resource for the koala, so field verification seems necessary unless more-accurate remote sensing methods such as hyperspectral imagery can be employed.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/f15111991/s1.

Author Contributions

D.C.J.: conceptualization, methodology, investigation, writing—original draft, formal analysis, and visualization. S.K.S.: methodology, writing—review and editing, and supervision. A.S.: methodology, writing—review and editing, supervision, and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

Funding was provided by the Australian Koala Foundation (Save The Koala Fund ABN 76 400 543 983). Additional funding was provided by the University of the Sunshine Coast. Both are part of a University of the Sunshine Coast Scholarship Agreement dated 24 June 2021.

Data Availability Statement

The data that support this study cannot be publicly shared at this time as further, different analyses from these data are intended. Some data elements may be shared upon reasonable request to the corresponding author if appropriate.

Acknowledgments

We thank Deborah Tabart OAM (AKF Chief Executive Officer) and the AKF Board for their dedication and support. We thank Dave Mitchell from the AKF for useful discussions. We thank Shaun Wilson, Rani Amaral Lustosa, and Kellie Newport for field assistance. We thank Andrew Crook for his generous donation to AKF to assist with funding this research.

Conflicts of Interest

The authors declare no conflicts of interest. While the AKF provided funding for Derek Johnson’s stipend and project costs, the AKF had no role in study design, data analysis, decision to publish, or preparation of the manuscript.

References

  1. Mpakairi, K.S.; Ndaimani, H.; Tagwireyi, P.; Zvidzai, M.; Madiri, T.H. Futuristic climate change scenario predicts a shrinking habitat for the African elephant (Loxodonta africana): Evidence from Hwange National Park, Zimbabwe. Eur. J. Wildl. Res. 2020, 66, 1. [Google Scholar] [CrossRef]
  2. Nagy-Reis, M.; Dickie, M.; Calvert, A.M.; Hebblewhite, M.; Hervieux, D.; Seip, D.R.; Gilbert, S.L.; Venter, O.; Demars, C.; Boutin, S.; et al. Habitat loss accelerates for the endangered woodland caribou in western Canada. Conserv. Sci. Pract. 2021, 3, e437. [Google Scholar] [CrossRef]
  3. IUCN, 2024. Red List. International Union for Conservation of Nature; Switzerland. Version 2021-1. Threats, Habitat and Ecology. Available online: http://www.iucnredlist.org (accessed on 7 November 2024).
  4. Leonard, K.M.; Koprowski, J.L. Effects of Fire on Endangered Mount Graham Red Squirrels (Tamiasciurus hudsonicus grahamensis): Responses of Individuals with Known Fates. Southwest. Nat. 2010, 55, 217–224. [Google Scholar] [CrossRef]
  5. Popa-Lisseanu, A.G.; Bontadina, F.; Ibáñez, C. Giant noctule bats face conflicting constraints between roosting and foraging in a fragmented and heterogeneous landscape. J. Zool. 2009, 278, 126–133. [Google Scholar] [CrossRef]
  6. Legge, S.; Woinarski, J.; Garnett, S.; Nimmo, D.; Scheele, B.; Lintermans, M.; Mitchell, N.; Whiterod, N.; Ferris, J. Rapid Analysis of Impacts of the 2019–20 Fires on Animal Species, and Prioritisation of Species for Management Response. Report Prepared for the Wildlife and Threatened Species Bushfire Recovery Expert Panel. 14 March 2020. Department of Agriculture, Water and the Environment. Australian Government. 2020. Available online: https://researchoutput.csu.edu.au/ws/portalfiles/portal/389233067/378471994_published_report.pdf (accessed on 7 November 2024).
  7. McGraw, W.S. Comparative locomotion and habitat use of six monkeys in the Tai Forest, Ivory Coast. Am. J. Phys. Anthropol. 1998, 105, 493–510. [Google Scholar] [CrossRef]
  8. Jackson, S. Koala: Origins of an Icon; Allen and Unwin: Crows Nest, NSW, Australia, 2007; ISBN 978-1-74237-323-2. Available online: https://books.google.com.au/books?hl=en&lr=&id=1lzE8IxPYy4C&oi=fnd&pg=PR3&dq=jackson+koala+origins+icon&ots=nFSmgUVort&sig=2RcptUrMzCRhfiV46q_zUPX1vs8&redir_esc=y#v=onepage&q=isbn&f=false (accessed on 7 November 2024).
  9. Boer, M.M.; de Dios, V.R.; Bradstock, R.A. Unprecedented burn area of Australian mega forest fires. Nat. Clim. Chang. 2020, 10, 171–172. [Google Scholar] [CrossRef]
  10. Ehsani, M.R.; Arevalo, J.; Risanto, C.B.; Javadian, M.; Devine, C.J.; Arabzadeh, A.; Venegas-quiñones, H.L.; Dell’oro, A.P.; Behrangi, A. 2019–2020 Australia fire and its relationship to hydroclimatological and vegetation variabilities. Water 2020, 12, 3067. [Google Scholar] [CrossRef]
  11. Collins, L.; Bradstock, R.A.; Clarke, H.; Clarke, M.F.; Nolan, R.H.; Penman, T.D. The 2019/2020 mega-fires exposed Australian ecosystems to an unprecedented extent of high-severity fire. Environ. Res. Lett. 2021, 16, 044029. [Google Scholar] [CrossRef]
  12. Bureau of Meteorology. Daily Rainfall: Perseverance Dam. Australian Government. 2020. Available online: https://reg.bom.gov.au/qld/index.shtml (accessed on 7 November 2024).
  13. Ward, M.; Tulloch, A.I.T.; Radford, J.Q.; Williams, B.A.; Reside, A.E.; MacDonald, S.L.; Mayfield, H.J.; Maron, M.; Possingham, H.P.; Vine, S.J.; et al. Impact of 2019–2020 mega-fires on Australian fauna habitat. Nat. Ecol. Evol. 2020, 4, 1321–1326. [Google Scholar] [CrossRef]
  14. Australian Koala Foundation. Koala Habitat Atlas. 2024. Available online: https://www.savethekoala.com/our-work/koala-habitat-atlas (accessed on 7 November 2024).
  15. Shabani, F.; Shafapourtehrany, M.; Ahmadi, M.; Kalantar, B.; Özener, H.; Clancy, K.; Esmaeili, A.; Da Silva, R.S.; Beaumont, L.J.; Llewelyn, J.; et al. Habitat in flames: How climate change will affect fire risk across koala forests. Environ. Technol. Innov. 2023, 32, 103331. [Google Scholar] [CrossRef]
  16. Williams, R.J.; Gill, A.M.; Bradstock, R.A. (Eds) Flammable Australia: Fire Regimes, Biodiversity 1173 and Ecosystems in a Changing World; CSIRO Publishing: Clayton, Australia, 2012. [Google Scholar]
  17. Phillips, S.; Wallis, K.; Lane, A. Quantifying the impacts of bushfire on populations of wild koalas (Phascolarctos cinereus): Insights from the 2019/20 fire season. Ecol. Manag. Restor. 2021, 22, 80–88. [Google Scholar] [CrossRef]
  18. Johnson, D.C.; Shapcott, A. Koala habitat forest recovery varies with fire severity. For. Ecol. Manag. 2024, 556, 121704. [Google Scholar] [CrossRef]
  19. Queensland Department of Environment and Science. Regional Ecosystem Description Database Maintained by Queensland Herbarium, Brisbane. Version 13. 2023. Available online: https://apps.des.qld.gov.au/regional-ecosystems// (accessed on 7 November 2024).
  20. Sattler, P.S.; Williams, R.; The Conservation Status of Queensland’s Bioregional Ecosystems. Queensland Environmental Protection Agency, Queensland Government. 1999. Available online: https://catalogue.nla.gov.au/catalog/2956630 (accessed on 7 November 2024).
  21. Nolan, R.H.; Rahmani, S.; Samson, S.A.; Simpson-Southward, H.M.; Boer, M.M.; Bradstock, R.A. Bark attributes determine variation in fire resistance in resprouting tree species. For. Ecol. Manag. 2020, 474, 118385. [Google Scholar] [CrossRef]
  22. Prior, L.D.; Foyster, S.M.; Furlaud, J.M.; Williamson, G.J.; Bowman, D.M.J.S. Using permanent forest plots to evaluate the resilience to fire of Tasmania’s tall wet eucalypt forests. For. Ecol. Manag. 2022, 505, 119922. [Google Scholar] [CrossRef]
  23. Fairman, T.A.; Nitschke, C.R.; Bennett, L.T. Too much, too soon? A review of the effects of increasing wildfire frequency on tree mortality and regeneration in temperate eucalypt forests. Int. J. Wildland Fire 2016, 25, 831–848. [Google Scholar] [CrossRef]
  24. Law, B.S.; Gonsalves, L.; Burgar, J.; Brassil, T.; Kerr, I.; O’Loughlin, C. Fire severity and its local extent are key to assessing impacts of Australian mega-fires on koala (Phascolarctos cinereus) density. Glob. Ecol. Biogeogr. 2022, 31, 714–726. [Google Scholar] [CrossRef]
  25. De Santis, A.; Chuvieco, E. GeoCBI: A modified version of the Composite Burn Index for the initial assessment of the short-term burn severity from remotely sensed data. Remote Sens. Environ. 2009, 113, 554–562. [Google Scholar] [CrossRef]
  26. Key, C.H.; Benson, N.C. Landscape Assessment (LA). In FIREMON: Fire Effects Monitoring and Inventory System; Lutes, D.C., Keane, R.E., Caratti, J.F., Key, C.H., Benson, N.C., Sutherland, S., Gangi, L.J., Eds.; General Technical Reports RMRS-GTR-164-CD; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2006; Available online: https://www.fs.usda.gov/treesearch/pubs/24066 (accessed on 7 November 2024).
  27. European Space Agency. Sentinel Online: MultiSpectral Instrument (MSI) Overview. 2023. Available online: https://sentiwiki.copernicus.eu/web/sentinel-2 (accessed on 7 November 2024).
  28. Gibson, R.; Danaher, T.; Hehir, W.; Collins, L. A remote sensing approach to mapping fire severity in south-eastern Australia using sentinel 2 and random forest. Remote Sens. Environ. 2020, 240, 0034–4257. [Google Scholar] [CrossRef]
  29. Han, A.; Qing, S.; Bao, Y.; Na, L.; Bao, Y.; Liu, X.; Zhang, J.; Wang, C. Short-Term Effects of Fire Severity on Vegetation Based on Sentinel-2 Satellite Data. Sustainability 2021, 13, 432. [Google Scholar] [CrossRef]
  30. Seydi, S.T.; Akhoondzadeh, M.; Amani, M.; Mahdavi, S. Wildfire Damage Assessment over Australia Using Sentinel-2 Imagery and MODIS Land Cover Product within the Google Earth Engine Cloud Platform. Remote Sens. 2021, 13, 220. [Google Scholar] [CrossRef]
  31. White, L.A.; Gibson, R.K. Comparing Fire Extent and Severity Mapping between Sentinel 2 and Landsat 8 Satellite Sensors. Remote Sens. 2022, 14, 1661. [Google Scholar] [CrossRef]
  32. Queensland Department of Environment and Science. Vegetation Management Act Regional Ecosystem and Remnant Map Queensland Government digital data. Version 13. 2023. Available online: https://www.qld.gov.au/environment/land/management/vegetation/maps (accessed on 7 November 2024).
  33. Chafer, C.J.; Noonan, M.; Macnaught, E. The post-fire measurement of fire severity and intensity in the Christmas 2001 Sydney wildfires. Int. J. Wildland Fire 2004, 13, 227–240. [Google Scholar] [CrossRef]
  34. Edwards, A.; Russell-Smith, J.; Meyer, M.; Burrows, N. Northern Fire Mapping Developing Robust Fire Extent and Severity Mapping Products for the Tropical Savannas. Presentation—Charles Darwin University, with Bushfire CRC and Bushfires NT. 2013. Available online: https://www.bushfirecrc.com/sites/default/files/managed/aedwards_fire_severity_jun13.pdf (accessed on 7 November 2024).
  35. Chuvieco, E.; Martín, M.P.; Palacios, A. Assessment of different spectral indices in the red-near-infrared spectral domain for burned land discrimination. Int. J. Remote Sens. 2002, 23, 5103–5110. [Google Scholar] [CrossRef]
  36. Parks, S.A.; Dillon, G.K.; Miller, C. A New Metric for Quantifying Burn Severity: The Relativized Burn Ratio. Remote Sens. 2014, 6, 1827–1844. [Google Scholar] [CrossRef]
  37. Parker, B.M.; Lewis, T.; Srivastava, S.K. Estimation and evaluation of multi-decadal fire severity patterns using Landsat sensors. Remote Sens. Environ. 2015, 170, 340–349. [Google Scholar] [CrossRef]
  38. IBM Corp. IBM SPSS Statistics for Windows, Version 28.0; IBM Corp: Armonk, NY, USA, 2021. [Google Scholar]
  39. Neldner, V.J.; Wilson, B.A.; Dillewaard, H.A.; Ryan, T.S.; Butler, D.W.; McDonald, W.J.F.; Richter, D.; Addicott, E.P.; Appelman, C.N. Methodology for Survey and Mapping of Regional Ecosystems and Vegetation Communities in Queensland Version 7.0. Updated December 2023b. Queensland Herbarium, Queensland Department of Environment and Science, Brisbane. 2023. Available online: https://www.qld.gov.au/environment/plants-animals/plants/herbarium/mapping-ecosystems (accessed on 7 November 2024).
  40. Neldner, V.J.; Niehus, R.E.; Wilson, B.A.; McDonald, W.J.F.; Ford, A.J.; Accad, A. The Vegetation of Queensland. Descriptions of Broad Vegetation Groups. Version 6.0. Queensland Herbarium and Biodiversity Science, Department of Environment and Science. 2023. Available online: https://www.qld.gov.au/environment/plants-animals/plants/ecosystems/broad-vegetation (accessed on 7 November 2024).
  41. NVIS Technical Working Group. Australian Vegetation Attribute Manual: National Vegetation Information System, Version 7.0; Bolton, M.P., deLacey, C., Bossard, K.B., Eds.; Department of the Environment and Energy: Canberra, Austria, 2017; Available online: https://www.dcceew.gov.au/environment/land/publications/australian-vegetation-attribute-manual-version-7 (accessed on 7 November 2024).
  42. Esri. ArcGIS, Version 10.8; Environmental Systems Research Institute, Inc.: Redlands, CA, USA, 2020. [Google Scholar]
  43. Avolio, M.L.; Forrestel, E.J.; Chang, C.C.; La Pierre, K.J.; Burghardt, K.T.; Smith, M.D. Demystifying dominant species. New Phytol. 2019, 223, 1106–1126. [Google Scholar] [CrossRef]
  44. Draper, F.C.; Costa, F.R.; Arellano, G.; Phillips, O.L.; Duque, A.; Macía, M.J.; Ter Steege, H.; Asner, G.P.; Berenguer, E.; Schietti, J.; et al. Amazon tree dominance across forest strata. Nat. Ecol. Evol. 2021, 5, 757–767. [Google Scholar] [CrossRef]
  45. Geoscience Australia. Digital Earth Australia. 2023. Available online: https://knowledge.dea.ga.gov.au/guides/reference/analysis_ready_data_corrections/ (accessed on 7 November 2024).
  46. rontierSI, 2021. AusCalVal Report: Establishing Australia as a Global Leader in Delivering Quality Assured Satellite Earth Observation Data. Available online: https://frontiersi.com.au/wp-content/uploads/2021/05/FrontierSI_AusCalVal_27052021_Final.pdf.pdf (accessed on 7 November 2024).
  47. Queensland Department of Environment and Science, 2023c. Fire Scar Mapping: Sentinel-2 Annual Fire Scars Queensland. Available online: https://www.qld.gov.au/environment/land/management/mapping/statewide-monitoring/firescar (accessed on 7 November 2024).
  48. Goodwin, N.R.; Collett, L.J. Development of an automated method for mapping fire history captured in Landsat TM and ETM+ time series across Queensland, Australia. Remote Sens. Environ. 2014, 148, 206–221. [Google Scholar] [CrossRef]
  49. Hardtke, L.A.; Blanco, P.D.; del Valle, H.F.; Metternicht, G.I.; Sione, W.F. Semi-automated mapping of burned areas in semi-arid ecosystems using MODIS time-series imagery. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 25–35. [Google Scholar] [CrossRef]
  50. Queensland Parks and Wildlife Service. Fire history—Queensland Parks and Wildlife Service. In History of planned burns and wildfires on and around the Protected Areas and Forests of Queensland; 2021. Available online: https://www.data.qld.gov.au/dataset/fire-history-queensland-parks-and-wildlife-service (accessed on 7 November 2024).
  51. Atlas of Living Australia. Supported by the Australian Government through the National Collaborative Research Infrastructure Strategy (NCRIS) and is hosted by the Commonwealth Scientific and Industrial Research Organisation (CSIRO). 2023. Available online: https://ala.org.au/ (accessed on 7 November 2024).
  52. Qarallah, B.; Othman, Y.A.; Al-Ajlouni, M.; Alheyari, H.A.; Qoqazeh, B.A.A. Assessment of Small-Extent Forest Fires in Semi-Arid Environment in Jordan Using Sentinel-2 and Landsat Sensors Data. Forests 2023, 14, 41. [Google Scholar] [CrossRef]
  53. Shvetsov, E.G.; Kukavskaya, E.A.; Buryak, L.V.; Barrett, K. Assessment of post-fire vegetation recovery in Southern Siberia using remote sensing observations. Environ. Res. Lett. 2019, 14, 055001. [Google Scholar] [CrossRef]
  54. Cocke, A.; Fulé, P.; Crouse, J. Comparison of burn severity assessment using Differenced Normalized Burn Ratio and ground data. Int. J. Wildland Fire 2005, 14, 189–198. [Google Scholar] [CrossRef]
  55. Escuin, S.; Navarro, R.; Fernández, P. Fire severity assessment by using NBR (normalized burn ratio) and NDVI (normalized difference vegetation index) derived from LANDSAT TM/ETM images. Int. J. Remote Sens. 2008, 29, 1053–1073. [Google Scholar] [CrossRef]
  56. Roy, D.R.; Boschetti, L.; Trigg, S.N. Remote sensing of fire severity: Assesing the performance of the normalized Burn ratio. IEEE Geosci. Remote Sens. Lett. 2006, 3, 112–116. [Google Scholar] [CrossRef]
  57. Wan, H.Y.; Cushman, S.A.; Ganey, J.L. The effect of scale in quantifying fire impacts on species habitats. Fire Ecol. 2020, 16, 9. [Google Scholar] [CrossRef]
  58. Queensland Department of Environment and Science. Queensland Agricultural Land Classes—Land Class A and B with Urban Mask. 2023. Available online: https://www.qld.gov.au/environment/land/management/soil (accessed on 7 November 2024).
  59. Queensland Department of Resources. Digital Elevation Model. 2023. Available online: https://www.arcgis.com/index.html (accessed on 7 November 2024).
  60. Google. Google Maps: Street View. 2023. Available online: https://www.google.com.au/maps (accessed on 8 November 2024).
  61. Thorley, J.; Srivastava, S.K.; Shapcott, A. What type of rainforest burnt in the South East Queensland’s 2019/20 bushfires and how might this impact biodiversity. Austral Ecol. 2023, 48, 616–642. [Google Scholar] [CrossRef]
  62. Matthews, A.; Lunney, D.; Gresser, S.; Maitz, W. Movement patterns of koalas in remnant forest after fire. Aust. Mammal. 2016, 38, 91–104. [Google Scholar] [CrossRef]
  63. Mitchell, D.L.; Soto-Berelov, M.; Jones, S.D. Remote sensing shows south-east Queensland koalas (Phascolarctos cinereus) prefer areas of higher tree canopy height within their home ranges. Wildl. Res. 2023, 50, 939–953. [Google Scholar] [CrossRef]
  64. Navarro, G.; Caballero, I.; Silva, G.; Parra, P.-C.; Vázquez, Á.; Caldeira, R. Evaluation of forest fire on Madeira Island using Sentinel-2A MSI imagery. Int. J. Appl. Earth Obs. Geoinf. 2017, 58, 97–106. [Google Scholar] [CrossRef]
  65. Alcaras, E.; Costantino, D.; Guastaferro, F.; Parente, C.; Pepe, M. Normalized Burn Ratio Plus (NBR+): A New Index for Sentinel-2 Imagery. Remote Sens. 2022, 14, 1727. [Google Scholar] [CrossRef]
  66. Gibson, R.K.; Hislop, S. Signs of resilience in resprouting Eucalyptus forests, but areas of concern: 1 year of post-fire recovery from Australia’s Black Summer of 2019–2020. Int. J. Wildland Fire 2022, 31, 545–557. [Google Scholar] [CrossRef]
  67. Filipponi, F. BAIS2: Burned Area Index for Sentinel-2. Proceedings 2018, 2, 364. [Google Scholar] [CrossRef]
  68. Mendonca, A.F.; Armond, T.; Camargo AC, L.; Camargo, N.F.; Ribeiro, J.F.; Zangrandi, P.L.; Vieira, E.M. Effects of an extensive fire on arboreal small mammal populations in a neotropical savanna woodland. J. Mammal. 2015, 96, 368–379. [Google Scholar] [CrossRef]
  69. Callaghan, J.; McAlpine, C.; Mitchell, D.; Thompson, J.; Bowen, M.; Rhodes, J.; De Jong, C.; Domalewski, R.; Scott, A. Ranking and mapping koala habitat quality for conservation planning on the basis of indirect evidence of tree-species use: A case study of Noosa Shire, south-eastern Queensland. Wildl. Res. 2011, 38, 89–102. [Google Scholar] [CrossRef]
  70. Lunney, D.; Gresser, S.E.; O’Neill, L.; Matthews, A.; Rhodes, J. The impact of fire and dogs on Koalas at Port Stephens, New South Wales, using population viability analysis. Pac. Conserv. Biol. 2007, 13, 189–201. [Google Scholar] [CrossRef]
  71. Orlando, C.G.; Montague-Drake, R.; Turbill, J.; Crowther, M.S. Megafires and koala occurrence: A comparative analysis of field data and satellite imagery. Aust. Mammal. 2024, 46, AM23054. [Google Scholar] [CrossRef]
  72. Lindenmayer, D.B.; McBurney, L.; Blanchard, W.; Marsh, K.; Bowd, E.; Watchorn, D.; Taylor, C.; Youngentob, K. Elevation, disturbance, and forest type drive the occurrence of a specialist arboreal folivore. PLoS ONE 2022, 17, e0265963. [Google Scholar] [CrossRef]
  73. Burns, E.L.; Lindenmayer, D.B.; Stein, J.; Blanchard, W.; McBurney, L.; Blair, D.; Banks, S.C. Ecosystem assessment of mountain ash forest in the Central Highlands of Victoria, south-eastern Australia. Austral Ecol. 2015, 40, 386–399. [Google Scholar] [CrossRef]
  74. Gill, A.; Ashton, D. The role of bark type in relative tolerance to fire of three central Victorian Eucalypts. Aust. J. Bot. 1968, 16, 491–498. [Google Scholar] [CrossRef]
  75. Collins, L. Eucalypt forests dominated by epicormic resprouters are resilient to repeated canopy fires. J. Ecol. 2020, 108, 310–324. [Google Scholar] [CrossRef]
  76. Higgins, S.I.; Bond, W.J.; Combrink, H.; Craine, J.M.; February, E.C.; Govender, N.; Lannas, K.; Moncreiff, G.; Trollope WS, W. Which traits determine shifts in the abundance of tree species in a fire-prone savanna? J. Ecol. 2012, 100, 1400–1410. [Google Scholar] [CrossRef]
  77. Queensland Department of Natural Resources, Mines and Energy. Esk Surface Geology Sheet 9343 Compilation Map The Geology on This Map Depicts Data from the Current MinesOnline Geodatabase Release Compiled on 16 July 2018. 2018. Available online: https://geoscience.data.qld.gov.au/data/report/cr039219 (accessed on 8 November 2024).
  78. Gibson, R.K.; Bradstock, R.A.; Penman, T.; Keith, D.A.; Driscoll, D.A. Climatic, vegetation and edaphic influences on the probability of fire across mediterranean woodlands of south-eastern Australia. J. Biogeogr. 2015, 42, 1750–1760. [Google Scholar] [CrossRef]
  79. McColl-Gausden, S.C.; Bennett, L.T.; Duff, T.J.; Cawson, J.G.; Penman, T.D. Climatic and edaphic gradients predict variation in wildland fuel hazard in south-eastern Australia. Ecography 2020, 43, 443–455. [Google Scholar] [CrossRef]
  80. Ndalila, M.N.; Williamson, G.J.; Bowman, D.M.J.S. Geographic Patterns of Fire Severity Following an Extreme Eucalyptus Forest Fire in Southern Australia: 2013 Forcett-Dunalley Fire. Fire 2018, 1, 40. [Google Scholar] [CrossRef]
  81. Chia, E.K.; Bassett, M.; Nimmo, D.G.; Leonard SW, J.; Ritchie, E.G.; Clarke, M.F.; Bennett, A.F. Fire severity and fire-induced landscape heterogeneity affect arboreal mammals in fire-prone forests. Ecosphere 2015, 6, art190. [Google Scholar] [CrossRef]
  82. Robinson, N.M.; Leonard, S.W.J.; Ritchie, E.G.; Bassett, M.; Chia, E.K.; Buckingham, S.; Gibb, H.; Bennett, A.F.; Clarke, M.F. REVIEW: Refuges for fauna in fire-prone landscapes: Their ecological function and importance. J. Appl. Ecol. 2013, 50, 1321–1329. [Google Scholar] [CrossRef]
  83. Berry, L.E.; Driscoll, D.A.; Stein, J.A.; Blanchard, W.; Banks, S.C.; Bradstock, R.A.; Lindenmayer, D.B. Identifying the location of fire refuges in wet forest ecosystems. Ecol. Appl. 2015, 25, 2337–2348. [Google Scholar] [CrossRef]
  84. Shaw, R.E.; James, A.I.; Tuft, K.; Legge, S.; Cary, G.J.; Peakall, R.; Banks, S.C. Unburnt habitat patches are critical for survival and in situ population recovery in a small mammal after fire. J. Appl. Ecol. 2021, 58, 1325–1335. [Google Scholar] [CrossRef]
  85. Durkin, L.K.; Moloney, P.D.; Cripps, J.K.; Nelson, J.L.; Macak, P.V.; Scroggie, M.P.; Collins, L.; Emerson, L.D.; Molloy, J.; Lumsden, L.F. Unburnt refugia support post-fire population recovery of a threatened arboreal marsupial, Leadbeater’s possum. For. Ecol. Manag. 2024, 551, 121487. [Google Scholar] [CrossRef]
  86. Berry, L.E.; Driscoll, D.A.; Banks, S.C.; Lindenmayer, D.B. The use of topographic fire refuges by the greater glider (Petauroides volans) and the mountain brushtail possum (Trichosurus cunninghami) following a landscape-scale fire. Aust. Mammal. 2015, 37, 39–45. [Google Scholar] [CrossRef]
  87. Goldingay, R.L.; Dobner, B. Home range areas of koalas in an urban area of north-east New South Wales. Aust. Mammal. 2014, 36, 74–80. [Google Scholar] [CrossRef]
  88. Wilson, K.; Pressey, R.L.; Newton, A.; Burgman, M.; Possingham, H.; Weston, C. Measuring and Incorporating Vulnerability into Conservation Planning. Environ. Manag. 2005, 35, 527–543. [Google Scholar] [CrossRef] [PubMed]
  89. Rhodes, J.R.; McAlpine, C.; Peterson, A.; Callaghan, J.G.; Lunney, D.; Possingham, H.P.; Mitchell, D.L.; Curran, T. Linking landscape ecology to planning for koala conservation. Aust. Plan. 2008, 45, 24–25. [Google Scholar] [CrossRef]
  90. Reside, A.E.; Briscoe, N.J.; Dickman, C.R.; Greenville, A.C.; Hradsky, B.A.; Kark, S.; Kearney, M.R.; Kutt, A.S.; Nimmo, D.G.; Pavey, C.R.; et al. Persistence through tough times: Fixed and shifting refuges in threatened species conservation. Biodivers. Conserv. 2019, 28, 1303–1330. [Google Scholar] [CrossRef]
  91. Queensland Department of Environment and Science. Koala Habitat Maps for South East Queensland. 2023. Available online: https://environment.des.qld.gov.au/wildlife/animals/living-with/koalas/mapping/koalamaps (accessed on 7 November 2024).
  92. Chia, E.K.; Bassett, M.; Leonard SW, J.; Holland, G.J.; Ritchie, E.G.; Clarke, M.F.; Bennett, A.F. Effects of the fire regime on mammal occurrence after wildfire: Site effects vs. landscape context in fire-prone forests. For. Ecol. Manag. 2016, 363, 130–139. [Google Scholar] [CrossRef]
  93. Matthews, A.; Lunney, D.; Gresser, S.; Maitz, W. Tree use by koalas (Phascolarctos cinereus) after fire in remnant coastal forest. Wildl. Res. 2007, 34, 84–93. [Google Scholar] [CrossRef]
  94. Ball, I.R.; Possingham, H.P.; Watts, M. Marxan and relatives: Software for spatial conservation prioritisation. In Spatial Conservation Prioritisation: Quantitative Methods and Computational Tools; Moilanen, A., Wilson, K.A., Possingham, H., Eds.; Oxford University Press: Oxford, UK, 2009; Volume 14, pp. 185–196. Available online: https://www.oup.com.au/books/others/9780199547760-spatial-conservation-prioritization (accessed on 7 November 2024).
  95. Marxan Software. 2024. Available online: https://marxansolutions.org/software/ (accessed on 7 November 2024).
  96. Bartels, S.F.; Chen HY, H.; Wulder, M.A.; White, J.C. Trends in post-disturbance recovery rates of Canada’s forests following wildfire and harvest. For. Ecol. Manag. 2016, 361, 194–207. [Google Scholar] [CrossRef]
  97. Bennett, L.T.; Bruce, M.J.; MacHunter, J.; Kohout, M.; Tanase, M.A.; Aponte, C. Mortality and recruitment of fire-tolerant eucalypts as influenced by wildfire severity and recent prescribed fire. For. Ecol. Manag. 2016, 380, 107–117. [Google Scholar] [CrossRef]
  98. Cansler, C.A.; McKenzie, D. How robust are burn severity indices when applied in a new region? Evaluation of alternate field-based and remote-sensing methods. Remote Sens. 2012, 4, 456–483. [Google Scholar] [CrossRef]
  99. Kattenborn, T.; Fassnacht, F.E.; Schmidtlein, S. Differentiating plant functional types using reflectance: Which traits make the difference? Remote Sens. Ecol. Conserv. 2019, 5, 5–19. [Google Scholar] [CrossRef]
Figure 1. Map showing the regional study extent and the location of the local study area (LSA) therein (centre panel and right-hand legend and inset). The forest types comprising 70% or greater of the vegetation in an area are as follows: GBS = grey gum, mountain blue gum, and stringybark; IPR = ironbark on ridges; IBM = ironbark and mountain blue gum on microgranite; BFA = blue gum flats on alluvium inland; BCF = blue gum flats on alluvium closer to the coast. The upper left-hand panel shows the LSA with the locations of the ‘Plotless sites’ used to identify forest types (Regional Ecosystems) and collect the GeoCBI samples. Raw GeoCBI continuous values range from 0 to 3 and were derived from GeoCBI values from the sites, but they were reclassified in this study into classes 0–5. The lower left-hand panel shows the nearest towns and road network.
Figure 1. Map showing the regional study extent and the location of the local study area (LSA) therein (centre panel and right-hand legend and inset). The forest types comprising 70% or greater of the vegetation in an area are as follows: GBS = grey gum, mountain blue gum, and stringybark; IPR = ironbark on ridges; IBM = ironbark and mountain blue gum on microgranite; BFA = blue gum flats on alluvium inland; BCF = blue gum flats on alluvium closer to the coast. The upper left-hand panel shows the LSA with the locations of the ‘Plotless sites’ used to identify forest types (Regional Ecosystems) and collect the GeoCBI samples. Raw GeoCBI continuous values range from 0 to 3 and were derived from GeoCBI values from the sites, but they were reclassified in this study into classes 0–5. The lower left-hand panel shows the nearest towns and road network.
Forests 15 01991 g001
Figure 2. Typical burn severities of eucalypt forest three months after the fire in the local study area. Ratings range from 0 to 5. (0) None. (1) Low (<2 m scorch)—note recovery of the ground layer. (2) Moderate low (2–5 m scorch)—note some epicormic shooting. (3) Moderate high (just trunks scorched >5 m OR a minority of crowns scorched). (4) High (most crowns scorched)—note recovery via epicormic shooting. Also note recovery of ground layer, which may look like canopy recovery from remote sensing. (5) Severe (crown foliage gone, ground bare)—note topkill and recovery via coppicing.
Figure 2. Typical burn severities of eucalypt forest three months after the fire in the local study area. Ratings range from 0 to 5. (0) None. (1) Low (<2 m scorch)—note recovery of the ground layer. (2) Moderate low (2–5 m scorch)—note some epicormic shooting. (3) Moderate high (just trunks scorched >5 m OR a minority of crowns scorched). (4) High (most crowns scorched)—note recovery via epicormic shooting. Also note recovery of ground layer, which may look like canopy recovery from remote sensing. (5) Severe (crown foliage gone, ground bare)—note topkill and recovery via coppicing.
Forests 15 01991 g002
Figure 3. Regional fire severity percentages of total area for each forest type tested (GBS, IPR, IBM, BFA, BCF) for 2019, the year of the fire. The regional study area defined by the Sentinel-2 scene is 12,056.04 km2. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11). GeoCBI burn classes: 0 = none, 1 = low, 2 = moderate low, 3 = moderate high, 4 = high, and 5 = severe.
Figure 3. Regional fire severity percentages of total area for each forest type tested (GBS, IPR, IBM, BFA, BCF) for 2019, the year of the fire. The regional study area defined by the Sentinel-2 scene is 12,056.04 km2. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11). GeoCBI burn classes: 0 = none, 1 = low, 2 = moderate low, 3 = moderate high, 4 = high, and 5 = severe.
Forests 15 01991 g003
Figure 4. Koala habitat areas suitable for refuge from fire. View covers the entire Sentinel-2 scene. Criteria are (1) koala habitat forest type; (2) fire severity class 0—unburnt; (3) proximity to burnt areas with moderate, high, and severe ratings—classes 2, 3, 4, and 5; and (4) proximity to recorded koala sightings.
Figure 4. Koala habitat areas suitable for refuge from fire. View covers the entire Sentinel-2 scene. Criteria are (1) koala habitat forest type; (2) fire severity class 0—unburnt; (3) proximity to burnt areas with moderate, high, and severe ratings—classes 2, 3, 4, and 5; and (4) proximity to recorded koala sightings.
Forests 15 01991 g004
Figure 5. Koala habitat areas suitable for refuge from fire. View is local study area and its surrounds. Criteria are (1) koala habitat forest type; (2) fire severity class 0—unburnt; (3) proximity to burnt areas with moderate, high, and severe ratings—classes 2, 3, 4, and 5; and (4) proximity to recorded koala sightings.
Figure 5. Koala habitat areas suitable for refuge from fire. View is local study area and its surrounds. Criteria are (1) koala habitat forest type; (2) fire severity class 0—unburnt; (3) proximity to burnt areas with moderate, high, and severe ratings—classes 2, 3, 4, and 5; and (4) proximity to recorded koala sightings.
Forests 15 01991 g005
Figure 6. Normalised Difference Vegetation Index (NDVI) measured for two years prior to the fire in 2019, and three years after the fire, across 88 plotless sites in the local study area.
Figure 6. Normalised Difference Vegetation Index (NDVI) measured for two years prior to the fire in 2019, and three years after the fire, across 88 plotless sites in the local study area.
Forests 15 01991 g006
Figure 7. Mean dNBR values for each forest type in each year. A negative dNBR indicates no burning, a dNBR of 0 is neutral (with no change between two years), and a dNBR >= 1 in this study indicates severe burning. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11).
Figure 7. Mean dNBR values for each forest type in each year. A negative dNBR indicates no burning, a dNBR of 0 is neutral (with no change between two years), and a dNBR >= 1 in this study indicates severe burning. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11).
Forests 15 01991 g007
Figure 8. Burn severity classes for koala habitat forest types (Regional Ecosystems) at the regional scale across one whole Sentinel-2 scene (covering an area of 100 km × 100 km) in 2019, immediately after the fire. Negative dNBR indicates unburnt, a dNBR of 0 means neutral (with no change between two years), and dNBR ≥ 1 in this study indicates severely burnt. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11).
Figure 8. Burn severity classes for koala habitat forest types (Regional Ecosystems) at the regional scale across one whole Sentinel-2 scene (covering an area of 100 km × 100 km) in 2019, immediately after the fire. Negative dNBR indicates unburnt, a dNBR of 0 means neutral (with no change between two years), and dNBR ≥ 1 in this study indicates severely burnt. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11).
Forests 15 01991 g008
Figure 9. Burn severity classes for koala habitat forest types (Regional Ecosystems) at the regional scale across one whole Sentinel-2 scene (an area of 100 km × 100 km) in 2020, one year after the fire. Negative dNBR indicates unburnt, an dNBR of 0 means neutral (with no change between two years), and dNBR ≥ 1 in this study indicates severely burnt. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11).
Figure 9. Burn severity classes for koala habitat forest types (Regional Ecosystems) at the regional scale across one whole Sentinel-2 scene (an area of 100 km × 100 km) in 2020, one year after the fire. Negative dNBR indicates unburnt, an dNBR of 0 means neutral (with no change between two years), and dNBR ≥ 1 in this study indicates severely burnt. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11).
Forests 15 01991 g009
Figure 10. Burn severity and recovery trend for forest type GBS (RE 12.12.23) and forest type IPR (RE 12.12.12) in the study region. Colour coding: red indicates year of fire, and other colours are years of recovery. Negative dNBR indicates unburnt, an dNBR of 0 means neutral (with no change between two years), and dNBR ≥ 1 in this study indicates severely burnt.
Figure 10. Burn severity and recovery trend for forest type GBS (RE 12.12.23) and forest type IPR (RE 12.12.12) in the study region. Colour coding: red indicates year of fire, and other colours are years of recovery. Negative dNBR indicates unburnt, an dNBR of 0 means neutral (with no change between two years), and dNBR ≥ 1 in this study indicates severely burnt.
Forests 15 01991 g010
Figure 11. Tree recovery responses from epicormic shooting versus dNBR value for each of 88 plotless sites within the local study area over three years for all forest types.
Figure 11. Tree recovery responses from epicormic shooting versus dNBR value for each of 88 plotless sites within the local study area over three years for all forest types.
Forests 15 01991 g011
Table 1. Regional Ecosystems observed within the study area that are koala habitats and were analysed in this study. * BVG = Broad Vegetation Group at 1:1,000,000 map scale [19]. Upland REs are listed first, followed by alluvial REs. Sample numbers are from the total of 88 plotless sample sites.
Table 1. Regional Ecosystems observed within the study area that are koala habitats and were analysed in this study. * BVG = Broad Vegetation Group at 1:1,000,000 map scale [19]. Upland REs are listed first, followed by alluvial REs. Sample numbers are from the total of 88 plotless sample sites.
Regional Ecosystem CodeForest Type Analysis CodeShort Description (QDES 2023b) and CommentsCommon Names of DominantsLandscapeBVG1 m *Number of Samples
12.12.12IPREucalyptus tereticornis, Corymbia intermedia, E. crebra +/− Lophostemon suaveolens—woodlands on Mesozoic to Proterozoic igneous rocks.ironbark and pink bloodwoodgranite ridges9 g9
12.12.23GBSEucalyptus tereticornis subsp. tereticornis or E. tereticornis subsp. basaltica +/− E. eugenioides—woodlands on crests, upper slopes, and elevated valleys and plains on Mesozoic to Proterozoic igneous rocks.stringybark and mountain blue gumgranite hills9 g66
12.8.16IBMEucalyptus crebra +/− E. melliodora, E. tereticornis—woodlands on Cainozoic igneous rocks.blue gum and ironbarkbasalt hills11 a8
12.3.3BFAEucalyptus tereticornis—woodlands and open forests on alluvial plains. Grow away from the coast.blue gumalluvial16 c5
12.3.11BCFThe only RE outside of the local study area that is also classified as a primary koala habitat (AKF 2023) and therefore included in this analysis. Eucalyptus tereticornis +/− Eucalyptus siderophloia, Corymbia intermedia—open forests on alluvial plains usually near the coast.blue gum with or without ironbark and pink bloodwoodalluvial16 c0
Table 2. Criteria for fire severity ratings (FSRs) [18].
Table 2. Criteria for fire severity ratings (FSRs) [18].
ClassFire Impact Severity Criterion
0none
1low (<2 m were scorched)
2moderately low (2–5 m were scorched)
3moderately high (just trunks were scorched at >5 m OR a minority of crowns were scorched)
4high (most crowns were scorched)
5severe (crown foliage is gone; ground is bare)
Table 3. Datasets used in this study.
Table 3. Datasets used in this study.
DatasetMap ScaleSourceDescription
Regional Ecosystem mapping1:25,000Queensland Herbarium, QLD Dept. Environment and Science 2023 Vegetation mapping of Queensland
Sentinel-220 m pixelsGeoscience Australia 2023Analysis-ready data from Digital Earth Australia
GeoCBI sample sites1:5000Johnson and Shapcott, 202488 plotless field sites
Fire severity rating (FSR) sample sites1:5000Johnson and Shapcott, 202488 plotless field sites
Fire recovery sample sites1:5000Johnson and Shapcott, 202488 plotless field sites
Fire history mapping1:100,000 nominal depends on sourceQueensland Parks and Wildlife Service 2021Mapping of reported fires in Queensland 1930 to 2023
Fire scar mapping20 m pixelsQueensland Dept. Environment and Science 2023Mapping of fires in Queensland based on dNBR 2019 to 2020
Table 4. Dates the field data were collected and Sentinel-2 scenes that were analysed.
Table 4. Dates the field data were collected and Sentinel-2 scenes that were analysed.
EventMonth and YearBest Image Date (No Clouds or Smoke)Satellite
Pre-fire (2 years)September 20177 September 20172B
Pre-fire (1 year) December 201821 December 20182B
Pre-fire (10 days)November 20196 November 20192B
Fire (not used)November 201916 November 2019 (mid-burn not analysed)2B
Fire (just after)December 20196 December 20192B
Recovery Year 1November 2020 20 November 20202B
Recovery Year 2May 202114 May 20212A
Recovery Year 3February 202218 February 20222A
Table 5. Sentinel-2 band wavelengths. Bands are grouped here in order of resolution (pixel or cell size, square.) Blue, green, and red are in the visible spectrum. NIR is near infrared. Veg red edge is for vegetation. Narrow NIR is a narrower bandwidth of NIR. SWIR is shortwave infrared and divided into two bands. SWIR—cirrus detects cirrus cloud reflection. Sentinel-2A (S2A) and Sentinel-2B (S2B) orbit on opposite sides of the Earth used to halve the data collection interval, decreasing it from 10 days to 5 days.
Table 5. Sentinel-2 band wavelengths. Bands are grouped here in order of resolution (pixel or cell size, square.) Blue, green, and red are in the visible spectrum. NIR is near infrared. Veg red edge is for vegetation. Narrow NIR is a narrower bandwidth of NIR. SWIR is shortwave infrared and divided into two bands. SWIR—cirrus detects cirrus cloud reflection. Sentinel-2A (S2A) and Sentinel-2B (S2B) orbit on opposite sides of the Earth used to halve the data collection interval, decreasing it from 10 days to 5 days.
Band NumberResolution (m)NameS2A Central Wavelength (nm)S2A Bandwidth (nm)S2B Central Wavelength (nm)S2B Bandwidth (nm)
210blue492.466492.166
310green559.836559.036
410red664.631664.931
810NIR832.8106832.9106
520veg red edge704.115703.816
620veg red edge740.515739.115
720veg red edge782.820779.720
8a20narrow NIR864.721864.022
1120SWIR1613.7911610.494
1220SWIR2202.41752185.7185
160coastal aerosol442.721442.221
960water vapour945.120943.221
1060SWIR—cirrus1373.5311376.930
Table 6. Classification scale for dNBR for use with Comparable Incremental Scale Reference (CISR). Class 1 is a conglomerated class including several unburnt classes. Classes 11 and 12 are central classes indicative of neutral conditions (not burnt, not recovering). Class 22 is a conglomerated class including several severely burnt classes. The dNBR value immediately post-burn will show areas with no change (unburnt or non-vegetated) and high values for burnt areas. The dNBR for subsequent years provides an estimate of vegetation recovery when compared with immediate post-burn dNBR.
Table 6. Classification scale for dNBR for use with Comparable Incremental Scale Reference (CISR). Class 1 is a conglomerated class including several unburnt classes. Classes 11 and 12 are central classes indicative of neutral conditions (not burnt, not recovering). Class 22 is a conglomerated class including several severely burnt classes. The dNBR value immediately post-burn will show areas with no change (unburnt or non-vegetated) and high values for burnt areas. The dNBR for subsequent years provides an estimate of vegetation recovery when compared with immediate post-burn dNBR.
ClassdNBR Central ValueBurn SeverityClassdNBR Central ValueBurn Severity
1≤−1Vigorous growth120.1Neutral or minimal change
2−0.9 130.2Slightly burnt
3−0.8 140.3
4−0.7 150.4
5−0.6 160.5
6−0.5 170.6
7−0.4 180.7
8−0.3 190.8
9−0.2 200.9
10−0.1Slight recovery211
110Neutral or minimal change22>1Severely burnt
Table 7. Summary of areas burnt within each forest type (Regional Ecosystem (RE)) within the regional study area (with an area of approximately 12,000 km2) in 2019. Burnt areas are classes 1–5 and exclude class 0, unburnt. Forest types: GBS = grey gum and mountain blue gum and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11). The REs studied were restricted to those mapped areas [32] of greater than 70% of a particular RE due to the limitations of the dNBR in detecting small proportions of the RE (when mixed with other REs).
Table 7. Summary of areas burnt within each forest type (Regional Ecosystem (RE)) within the regional study area (with an area of approximately 12,000 km2) in 2019. Burnt areas are classes 1–5 and exclude class 0, unburnt. Forest types: GBS = grey gum and mountain blue gum and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11). The REs studied were restricted to those mapped areas [32] of greater than 70% of a particular RE due to the limitations of the dNBR in detecting small proportions of the RE (when mixed with other REs).
Forest TypeGBSIPRIBMBFABCFTotal
RE12.12.2312.12.1212.8.1612.3.312.3.11
Total area in bioregion Ha ^15,00053,00033,00038,00040,000179,000
Total study area of RE > 70% *5744291058486769319924,470
Proportion within bioregion %38.05.517.717.88.013.7
Burnt study area Ha *18551613741112402
Burnt study area % *32.35.56.40.10.039.8
^ Total area of all of the REs in the whole bioregion (QDES 2023b), including smaller percentages of mapped patches of remnant vegetation (from 100% to as little as 5% of the patches). * Based on REs constituting 70% or greater of a mapped patch of remnant vegetation within the regional study area covered by one Sentinel-2 scene (approximately 12,000 km2).
Table 8. Regional fire severity impact with classes based on the GeoCBI from the local study area for the year of the fire, 2019, and those same classes applied to recovery in the following years. The regional study area defined by the Sentinel-2 scene is 1,2056.04 km2. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11). GeoCBI burn classes: 0 = none, 1 = low, 2 = moderate low, 3 = moderate high, 4 = high, and 5 = severe. Total areas for each RE were subject to minor variation in later years. Values are rounded.
Table 8. Regional fire severity impact with classes based on the GeoCBI from the local study area for the year of the fire, 2019, and those same classes applied to recovery in the following years. The regional study area defined by the Sentinel-2 scene is 1,2056.04 km2. Forest types: GBS = grey gum, mountain blue gum, and stringybark (RE 12.12.23); IPR = ironbark on ridges (RE 12.12.12); IBM = ironbark and mountain blue gum on microgranite (RE 12.8.16); BFA = blue gum flats on alluvium inland (RE 12.3.3); BCF = blue gum flats on alluvium closer to the coast (RE 12.3.11). GeoCBI burn classes: 0 = none, 1 = low, 2 = moderate low, 3 = moderate high, 4 = high, and 5 = severe. Total areas for each RE were subject to minor variation in later years. Values are rounded.
Total Ha All GeoCBI Classes5744291058486769319924,470
2019
GeoCBI classGBS HaIPR HaIBM HaBFA HaBCF HaTotal Ha
03889274954736759319822,069
13776612021566
2321398610448
3259236210345
4332226420421
5566104140621
Total burnt area in Ha (classes 1–5)18551613741112402
2020
GeoCBI classGBS HaIPR HaIBM HaBFA HaBCF HaTotal Ha
05744290958506769317824,451
1000099
2000055
3000055
4000033
5000000
Total burnt area in Ha (classes 1–5)00002121
2021
GeoCBI classGBS HaIPR HaIBM HaBFA HaBCF HaTotal Ha
05744291058506769319824,472
1000000
2000000
3000000
4000000
5000000
Total burnt area in Ha (classes 1–5)000000
2022
GeoCBI classGBS HaIPR HaIBM HaBFA HaBCF HaTotal Ha
05651290858496761319324,362
128104134
224101127
319101121
419102122
5301218
Total burnt area in Ha (classes 1–5)9321115113
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Johnson, D.C.; Srivastava, S.K.; Shapcott, A. Forest Fire Severity and Koala Habitat Recovery Assessment Using Pre- and Post-Burn Multitemporal Sentinel-2 Msi Data. Forests 2024, 15, 1991. https://doi.org/10.3390/f15111991

AMA Style

Johnson DC, Srivastava SK, Shapcott A. Forest Fire Severity and Koala Habitat Recovery Assessment Using Pre- and Post-Burn Multitemporal Sentinel-2 Msi Data. Forests. 2024; 15(11):1991. https://doi.org/10.3390/f15111991

Chicago/Turabian Style

Johnson, Derek Campbell, Sanjeev Kumar Srivastava, and Alison Shapcott. 2024. "Forest Fire Severity and Koala Habitat Recovery Assessment Using Pre- and Post-Burn Multitemporal Sentinel-2 Msi Data" Forests 15, no. 11: 1991. https://doi.org/10.3390/f15111991

APA Style

Johnson, D. C., Srivastava, S. K., & Shapcott, A. (2024). Forest Fire Severity and Koala Habitat Recovery Assessment Using Pre- and Post-Burn Multitemporal Sentinel-2 Msi Data. Forests, 15(11), 1991. https://doi.org/10.3390/f15111991

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop