Figures
Abstract
The influence of plant traits on forest fire behaviour has evolutionary, ecological and management implications, but is poorly understood and frequently discounted. We use a process model to quantify that influence and provide validation in a diverse range of eucalypt forests burnt under varying conditions. Measured height of consumption was compared to heights predicted using a surface fuel fire behaviour model, then key aspects of our model were sequentially added to this with and without species-specific information. Our fully specified model had a mean absolute error 3.8 times smaller than the otherwise identical surface fuel model (p < 0.01), and correctly predicted the height of larger (≥1 m) flames 12 times more often (p < 0.001). We conclude that the primary endogenous drivers of fire severity are the species of plants present rather than the surface fuel load, and demonstrate the accuracy and versatility of the model for quantifying this.
Citation: Zylstra P, Bradstock RA, Bedward M, Penman TD, Doherty MD, Weber RO, et al. (2016) Biophysical Mechanistic Modelling Quantifies the Effects of Plant Traits on Fire Severity: Species, Not Surface Fuel Loads, Determine Flame Dimensions in Eucalypt Forests. PLoS ONE 11(8): e0160715. https://doi.org/10.1371/journal.pone.0160715
Editor: Gil Bohrer, The Ohio State University, UNITED STATES
Received: April 1, 2016; Accepted: July 22, 2016; Published: August 16, 2016
Copyright: © 2016 Zylstra et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All data not contained in the manuscript and supporting materials are available from https://github.com/pzylstra/ffm_cpp/tree/Zylstra2016. A link to this is provided in the methods and in S1 Table.
Funding: All work was funded by the NSW Office of Environment and Heritage (http://www.environment.nsw.gov.au) as part of an ongoing contract with the University of Wollongong Centre for Environmental Risk Management of Bushfires. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Plants may affect forest fire behaviour by influencing the quantity and flammability of surface fuel, the three dimensional structure of the forest, and through the flammability of their live parts. The effect of variation in dominant species on surface fine litter load is well established [1] and evidence of species effects on flammability is accumulating [2–5]. The importance of forest structure to fire behaviour and severity, and the role of species in determining this, has however received less attention. Consequently, although there is broad acceptance of a link between this and stand flammability, only a small number of studies have demonstrated this link for a limited range of structural traits, e.g. [6–8]. Crown fire modelling has led the way in this regard [9,10].
The role of leaf traits in influencing fire behaviour via flammability remains uncertain. While strong evidence exists for the influence of chemistry, moisture and morphology on leaf flammability [11–13] and laboratory studies have found correlations with some observed fire behaviour characteristics [14–16], the outcomes of variation in these leaf attributes has not been mechanistically defined for individual plants [17], and remains largely qualitative and incomplete for forest stands [18]. In some cases, conclusions drawn from observations appear to conflict. Belcher et al for example characterised sites with broader leaves as less flammable due to a positive correlation between leaf width and time to ignition [19], while others noted an opposite correlation with site flammability, due to the influence of longer or larger leaves in aerating litter layers [4,5].
Given these uncertainties, it has been argued that fire behaviour models can disregard the influence of certain flammable properties of plant parts. Alexander and Cruz [16] for instance proposed that variations in ignitability produced by different leaf moisture values will be overwhelmed by the larger heat fluxes involved in fires under field conditions, and therefore have little influence on rates of spread. However empirical evidence to support this prediction may be lacking [20]. The mechanisms underpinning this prediction also require further elucidation [21]. Schwilk [22] has argued that scaling from plant traits to ecosystem effects is a fundamental goal of functional ecology, but mechanisms for extrapolating flammable traits of plants to large scales, as encompassed in fire behaviour science, are lacking. Here we address this problem through the presentation and validation of a model that incorporates mechanisms for scaling from flammability of individual leaves to forest stand-level fire behaviour attributes such as flame height.
Our aim was to show that stand-level characteristics of fire behaviour could be successfully predicted on the basis of knowledge of the flammable properties of plant species and the interaction of these properties that arise from considering the three dimensional arrangement of whole plants and their relevant, flammable features, such as leaves and stems. Specifically, we predicted flame height using a modelling approach that accounted for the arrangement of plants and their flammable parts plus the influence of the morphology of these parts on flammability. We then tested predictions derived from the model against field estimates of char height or the height to which leaves are consumed by fire following a wildfire, in a range of forest sites across a mountainous region of south eastern Australia.
The study
Plant traits can potentially affect fire characteristics at the whole plant and stand scale via several fundamental processes. First, the properties of individual leaves such as size, shape, moisture content and chemical composition [11–13,23] determine their ignitability, combustibility and sustainability [24]. Second, the three-dimensional arrangement and spacing of leaves and fine stems either as litter or as live foliage, will affect the likelihood of fire propagating or spreading through the plant via various mechanisms of drying and heat transfer (radiation, conduction, convection—(radiation, conduction, convection, [9,25]). Third, the three dimensional arrangement and spacing of leaves will affect the nature of environmental influences on ignition and spread of fire among fuel elements via effects on insolation and wind [26–28]. Many fire behaviour models only partially represent these processes and arrangements, if at all. The Rothermel model [29] for example includes some determinants of leaf flammability but models heat transfer using only radiation and without explicitly considering the nature of gaps between individual fuel elements. Australian fire behaviour models based on surface ‘fuel load’ or the weight of fine dead material on the ground surface e.g. [30,31] do not explicitly consider any of these influences, whereas the more recent “Project Vesta” model utilises a small subset of vegetation structure and flammability properties to subjectively rank the flammability of some strata into scores [6].
In this study we tested the hypothesis that flammability properties of leaves will significantly affect characteristics of fire at the scale of a forest stand via the mechanisms listed above. Specifically, we examined how variations in flame propagation and resultant height of flames across multiple strata of plants within forest stands emerge from variations in characteristics of component species and their inherent effects on flammability and the spatial arrangement of leaves. We explored this using a comparative modelling approach to examine the sensitivity of predicted variations in the propagation and dimensions of flames to different representations of the structure and composition of eucalypt forest stands. We compared predictions from the models against ground-measured and remotely-sensed estimates of fire severity, for which we used height of combustion in vegetation as an index of flame dimensions [32]. Char height [33], or the height to which leaves were blackened, has been used as an estimate of flame height in some studies, but post-fire measurements were required immediately after the fire [34]; an option that was not available to us. Scorch height differs from char height in that air temperatures above a fire may be sufficient to kill the cells in leaves without producing char. This is frequently considered in severity measurements [32], but relating it back to flame height modelling requires the use of a scorch height model of unknown reliability.
Our measurements were taken within the perimeter of a major wildland fire which burned across diverse terrain, forest and weather conditions. These comparisons allowed us to test whether incorporation of leaf characteristics of individual species and their effects on ignitability, combustibility, sustainability, heat transfer and wind profiles, improved predictions of flame dimensions.
Materials and Methods
We conducted our study in two stages to quantify the influence of plant species composition and structural manifestation on forest fire behaviour. First, we proposed and tested a site-based biophysical, mechanistic model of fire behaviour as influenced by the species-driven properties of forest structure and leaf traits, and in the second we analysed the model outputs to identify the conditions under which leaf traits were important determinants of fire behaviour, and which groups of traits had the most influence
Model description and validation
The Forest Flammability Model [35,36], (FFM) was used to explore effects of vegetation structure and leaf flammability on flame propagation and dimensions. Full code, instructions and data required to replicate this study are available at https://github.com/pzylstra/ffm_cpp/tree/Zylstra2016, and further explanation is provided in S1 Table.
The FFM predicts fire behaviour from leaf traits and plant structure in the following way:
- When a leaf ignites it becomes a heat donor, producing a convective plume that decreases in temperature with distance from the source along a vector defined by the flame angle. The heat output of the burning leaf (combustibility) determines the dimensions of the heat plume and the nature of the temperature gradient within it (Fig 1). The duration of burning (i.e. the sustainability of flaming in that leaf) defines the period for which a given temperature is maintained.
- Leaves ignite if the period of flame duration produced by the donor leaf exceeds the time to ignition for the receiver. This is a function of both the temperature at that point in the plume and the ignitability of the leaf (Fig 1). Leaf ignitability in this study is represented by the Ignitability Coefficient, where IC = leaf moisture (% Oven Dry Weight) * thickness (mm) / number of sides on the leaf. As IC increases, ignitability decreases.
- This process is repeated at escalating scales from leaf to branch, plant and plant stratum.
The flame from the donor produces a convective plume following a direction described by the flame angle (broken arrow), where the temperature of the air in the plume decreases with distance from the donor (solid curve) in a pattern determined by the flame produced from that leaf. The time of heating required for ignition of the receiver increases as the temperature decreases, at a rate determined by the ignitability of the leaf. The plume temperature model is taken from [37], and the time to ignition modelled from [35], where ignitability is a function of plume temperature and the Ignitability Coefficient (IC = leaf moisture (% Oven Dry Weight) * thickness (mm) / number of sides on the leaf).
Combustibility and ignitability therefore interact in the model to determine the depth of ignition (distance from the flame at which leaves can be ignited), which together with the density of foliage gives the number of leaves ignited for each one-second time step. The angle at which the burning plume intersects a plant crown or stratum defines a plume pathway or potential depth of foliage that can be ignited. This angle is adjusted at each time step as wind speed and flame dimensions affect the flame angle, so that the plume pathway evolves over time. The total number of leaves burning in a time step is the sum of leaves burning in the previous step plus those newly ignited, minus the leaves that have extinguished as determined by their sustainability properties. The resulting flame length for that step is a function of the number of leaves burning, the combustibility of those leaves and physical processes of heat transfer and air entrainment between them.
In vegetation with multiple plant strata, the severity of a fire expressed as the height of combustion [32] will be determined by the capacity of flames from each stratum to ignite strata above them (Fig 2A). This will be a function of the heat produced by the donor strata, the heat required for ignition by the receiver stratum and the spacing between donor and receiver, as convective heat dissipates over distance. The structure and composition of the forest influences all these processes via the effects of the arrangement and spacing, dimensions, moisture and chemistry of leaves (Fig 2B and 2C).
Wind speed is shown above and below the canopy by the solid arrows, with thicker and longer arrows showing greater wind speed. The trajectory of convective heat transfer is shown by the broken line. Four scenarios b to e are shown relative to the left scene a. The convective plume produced by the donor plant in a. intersects the receiver; however it is insufficient to ignite it. This is changed when in b the donor flammability is increased to give a larger flame that ignites the receiver, and in c when the donor flame is the same as in a but the receiver ignitability is greater. The flammability of the plants in scenarios d and e is the same as in a, but the wider tree spacing has reduced the overstorey sheltering so that the wind speed is greater at the level of the flame. This directs the plume through its neighbouring plants in d so that the flame depth is increased and the resulting larger flame ignites the receiver stratum. In scenario e, the plume passes over neighbouring shrubs or elevated stratum [7] so that they are not ignited and the flame dimensions remain unchanged from a The more acute angle of the plume, however, increases the distance to the receiver stratum so that the heat dissipates and that stratum is even less likely to ignite than in a.
Plants in strata above a flame but not burning may influence the behaviour of the fire via effects on wind speed in the burning strata. Such effects will be a function of the amount of shelter provided by the overhead foliage, which in the FFM is modelled from the Leaf Area Index (LAI) calculated from leaf size, density on the branch, ramification and spacing of branches in the non-burning strata above the flames, along with the physical size and spacing of plant crowns. These effects on wind may, in turn, affect the angle of the heat plume and thereby the length of the plume pathway, the length of flame produced, and the angle and distance to receiver leaves (Fig 2D and 2E). Thus leaf traits may affect flame dimensions and propagation via the endogenous mechanisms of donor flammability and receiver ignitability, and the exogenous (i.e. environmental effects) mechanism of overstorey sheltering. The model represents both these exogenous and endogenous mechanisms and their effects on flame characteristics and propagation.
The FFM is initiated using a model of surface fire spread that produces a pilot flame to potentially ignite above-ground fuel strata. Species-specific leaf traits are known to affect surface fire behaviour (see introductory discussion), however at this point no model exists that incorporates these effects. In the absence of such a model, this study uses the Burrows model [38] for this purpose, as its development in a controlled laboratory setting was considered to more accurately reproduce the behaviour of flames burning in eucalypt leaf litter alone, as opposed to other leaf litter fire behaviour models developed in field situations [30,31]. When incorporated into the FFM, this model was constrained to the approximate domain of the data from which it was constructed (Figure 7 in [39]), and flame heights were modelled using an angle derived from flame length and wind speed [40] to remain consistent with the rest of the FFM.
Model validation.
We validated the model by investigating whether variations in structure and species composition of forest stands affected stand-level flame dimensions as predicted. We compared outputs from three modelling treatments with field observations of the height of combustion (height above the ground at which leaves and stems were consumed by fire) as a proxy for flame height in eucalypt forests of south-eastern Australia. The comparative modelling approach enabled us to systematically contrast predictions of stand-level flame dimensions produced by models that treated the representation of fuel with successively increasing levels of complexity. These ranged from surface litter fuel only to full consideration of structural arrangement and species-level leaf dimensions (Tables 1 and 2). This approach provided a test of the null hypothesis that explicit representation of the flammable characteristics of leaves of individual species within a spatial framework does not explicitly improve the capacity to model flame dimensions.
We tested the null hypothesis by comparing the three sets of predictions against estimates of height of combustion across a range of sites burnt by a major wildland fire in south-eastern Australia (Fig 3). Mountainous terrain and diverse forest types, along with variation in ambient weather conditions at the time of burning, resulted in marked variation in fire behaviour among sites. This provided a wide range of fire severity for the evaluation of predictions representing the influences of vegetation structure and leaf traits affecting fire properties at the stand level.
Three scenarios were modelled to address the hypothesis that vegetation structure and leaf traits are important to fire behaviour; these were:
- Surface fuel load only (F)–no structural or leaf traits effects of standing plants considered;
- Surface fuel load and stand structure (FS)–as for (F), but considering structural effects of above ground plants but not effects of leaf traits; and
- Surface fuel load, stand structure and leaf traits (FSL)–as for (FS), but also considering leaf traits of component species.
Study sites.
The derivation of various model inputs and estimation of fire severity or height of combustion were made following a large wildfire in January 2003 in the Brindabella Range near Canberra (Fig 1). A subset of 58 sites was sampled from a larger set of 130 that were originally established prior to the fire for a vegetation survey [42–45]. Sites were selected that had been burned at a known time, and for which adequate data existed to provide fuel inputs. The examined sites were situated in eight different montane to sub-alpine classes of forest dominated by Eucalyptus species of varying height (8 to 35 m, Table 3). The date and time (morning / afternoon) which the sites were burnt was estimated from maps of fire spread based on on-ground observations and remote sensing (multispectral line-scans from aircraft) (Australian Capital Territory Emergency Services unpublished data). All sites were re-surveyed and photographed after the fire in 2003 [42].
Weather conditions at the time of burning in each site were estimated using gridded weather surfaces for parameters used in the various models (Table 4). Slope and wind speed were calculated relative to the direction of fire spread at each site, so that fire flanks burning at right angles to the wind were given a wind value of zero, and downslope spread was given a negative slope value. Fuel and structural parameters required by the models were calculated from the initial site surveys and photographs (Tables 1 and 2). Leaf characteristics of 53 dominant species used in the FFM were taken from published sources [35,46–49], measurements of herbarium specimens and measurements taken in the field during 2013 at the study sites and in the Australian National Botanic Gardens at Canberra (S5 Table).
Field estimation of flame heights.
In the absence of definitive evidence, we estimated a range for the flame heights at each site rather than assigning a single value. The minimum was set by the height of combustion, which was measured at the study sites from multiple sources: i.e. remotely sensed difference Normalised Burn Ratio (dNBR,e.g. [32]) and estimates of fine stem consumption from the post-fire survey and photographs in 2003 [42]. In the field survey, a stratum was classed as burnt if both leaves and fine stems were absent at the time of post-fire surveys, and the minimum flame height was derived from either this or the presence of char on smooth barked trees as discussed by Alexander and Cruz [52]. Flame height was estimated as being within the bounds defined by the height of combustion (loss of leaves and fine twigs) and the lowest of either the surveyed scorch height or the base of the next unconsumed stratum. Estimates of fire severity derived from dNBR using remote sensing (i.e. LANDSAT imagery at 25 m resolution) were obtained from a major survey after the fires [53], and where field records and photographs were unclear, the dNBR was used to identify whether elevated and crown strata were burnt or unburnt using the classes in the Barrett survey. “Low” severity indicated 50% burn in the elevated layer and all higher classes indicated complete burn, and “high’ indicated 50% burn of the canopy, so that classes below this were unburnt and those above were completely burnt.
Flame height prediction.
Predictions of flame height for each site were produced using site-specific exogenous conditions (Table 4) and addressing the three treatments already described. Treatment F was modelled using [38] alone, FS was modelled using the FFM, but leaf traits were averaged into the two groups monocotyledon and dicotyledon and internal structure of plants was not characterised. FSL was also modelled using the FFM, with incorporation of species-specific estimates of leaf traits and using internal structure as summarised in Table 2 and listed in S5 Table.
All treatments used wind speed adjusted for overstorey sheltering by vegetation, and this was summarised for each site using the modelled wind reduction factors (WRF: above-canopy wind speed / wind speed at 1.5 m above ground level) as approximations of the full vertical wind-speed profiles used by the FFM for each site.
The predictions produced from each approach were then compared with the site-based upper and lower limits of flame height to determine whether the inclusion of leaf traits and their arrangement within various strata improved the accuracy of predictions.
The absolute error (AE) of predictions was estimated to identify treatment bias, and the mean absolute error (MAE) to provide a natural and unambiguous measure of model performance [54]. To better identify the nature of treatment errors, the proportion of correct predictions (PCP) was calculated for all sites, and PCP1 for a subset of sites where measured flames were at least 1m in height. The MAE, PCP and PCP1 were then compared between successive inclusions of structure and leaf traits using paired t-tests to determine whether any apparent improvements were statistically significant for those statistics. This approach was chosen in preference to a more standard nested-model comparison, as in addition to identifying the better treatment, it quantified the accuracy of each with statistics that we considered meaningful to model application.
Sensitivity of flame height to environmental conditions and specific leaf traits
In the second stage, we examined how the inclusion of leaf traits affected the relative importance of the three key mechanisms represented in the FFM (i.e. overstorey sheltering, donor flammability, receiver ignitability) in predicting estimated flame height. This was done by comparing the way that overall flame height changed between the FS and FSL treatments in response to parameters representing the three key processes, along with a range of other influences representing characteristics of the study sites, weather at the time of burning and other measures of vegetation/fuel structure (Table 5).
To achieve this, we performed Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis [55] to identify which of the three mechanisms most influenced differences in predicted flame height from the mean and species models. A set of candidate predictors was derived to represent each of the three mechanisms (Table 6); further predictors were added to represent site and environmental influences, giving 12 candidate predictors. All predictors were continuous variables.
LASSO regression was used to simultaneously fit models and select those subsets of predictors having most influence, with the coefficients for other predictors being forced to zero. The response variable was predicted flame height under FS minus that under FSL. Additive regression models were fitted using the “glmnet” package [56] for R [57], with the single tuning parameter chosen to either minimize prediction error as estimated by random k-fold cross-validation, or to find the most parsimonious regression model with a prediction error within one standard error of the minimum value. We used the cv.glmnet function with k = 10.
Trial regressions showed that there was substantial variation in the predictors chosen when the fitted model was selected using minimum prediction error (“min” rule) in repeated runs of the cross-validation procedure. In contrast, only one predictor (Delta donor Fh) was present in models selected with the one-standard-error rule (“1se” rule). This indicated a single dominant variable with a number of secondary variables of much less influence, whose inclusion depended on the random allocation of sites to training and test folds during cross-validation. Since our purpose was interpretation rather than prediction, we did not wish to simply discard these secondary variables. To investigate them further we undertook an ensemble regression analysis in which 1000 LASSO models were selected using the “min” rule, and the frequency of inclusion of each predictor in the fitted additive models was recorded along with statistics (mean, upper 97.5% and lower 2.5% quantiles) describing the characteristic regression coefficients for each.
Results
Field estimations of flame height
Observed fire severity ranged from surface (i.e. consumption of litter fuel only) at two sites to crown fire at four sites (Fig 4, S6 Table). Estimated flame heights ranged from 10 cm to more than 20 m with an average of 4.0 m (σ = 5.9 m). The range in estimated flame height varied between sites, with large levels of uncertainty evident when the gap between consumed vegetation strata and either the scorch height or the unconsumed strata above was large. The majority of sites (62%) were burnt by fires spreading either up or downslope with the wind behind them (usually head fires), 19% burnt downslope and against the wind (tail fires) and the remainder burnt either against the wind or across it, regardless of slope (flank fires).
Using surface fuel only (F), all predictions were for low flame heights (R2 = 0.11). Including structure with that fuel (FS) enabled the prediction of large flames, but the accuracy was low (R2 = 0.24). The inclusion of leaf traits (FSL) significantly improved the accuracy of predictions (R2 = 0.80), producing a MAE 3.80 times smaller than F, (p < 0.01, paired t-test), and 4.67 times smaller than FS (p < 0.01). The line of exact agreement is shown as solid, and the trend of the data is shown as a broken line, with R2 reported under the treatment name.
Flame height prediction
The successive inclusion of information on structure and other plant attributes improved the accuracy of predictions of flame dimensions relative to field estimates (Fig 4, Table 7). Flames predicted on the basis of surface fuel load only (F) were consistently small resulting in a large proportion (i.e. 31%) less than 1 m in height and a ME of -1.47 (Fig 5). As a result of the frequency of small observed flames, the MAE for these predictions was relatively low (1.53 m) despite the fact that it failed to predict any of the large flames that were observed.
As per [58], the box plots indicate data range, quartiles, and median, and the shaded area shows the density trace.
The inclusion of vegetation structural attributes without leaf trait information (FS) resulted in the prediction of large flames not previously predicted on the basis of F. However, in several sites where small flames were estimated, very large flames were predicted (Fig 4). As a result the MAE was higher for FS than for F (1.23 times the error, N.S., Fig 5). The PCP for FS was, however, 1.94 times higher than for F (p < 0.001, Fig 6), and PCP1 was 9.50 times higher (p < 0.001), which was indicative of the inability of F to predict these flames.
PCP calculations were based on all 58 sites; PCP1 calculations were based on the 40 sites with flame heights ≥1m. Error bars show one standard error above and below the mean.
The accuracy of predictions was considerably improved with the inclusion of leaf traits (FSL), removing most of the false predictions of large flames (Fig 4). This provided an improvement of the MAE (0.40 m), which was 3.80 times more accurate than F (p < 0.01) and 4.67 times more accurate than FS (p < 0.01, Fig 5). PCP and PCP1 were respectively 2.25 and 12.00 times higher for FSL than for F (p < 0.001 for both), and 1.16, 1.26 times larger than for FS (p < 0.1 for both, Fig 6).
Sensitivity of predicted flame height to environmental conditions and specific leaf traits
The most frequent (>70%) predictors included in the LASSO regression analyses were donor flammability (delta donor flame height), wind velocity, overstorey shelter (delta sum LAI) and dead fuel moisture content (DFMC); with delta canopy LAI and slope affecting the regression at moderate (>30%) frequency (Table 8). Receiver ignitability affected the predictions to only a minor degree (< 10%). Donor flammability (delta donor flame height) and overstorey shelter (delta LAI) were the most influential of these predictors, based on the size of the mean regression coefficients (Table 8). In comparison, effects of other influential predictors were weak, based on the size of their regression coefficients (Table 8)
Predictor groups are DF–donor flammability, RI–receiver ignitability, OS–overstorey shelter and General–exogenous, structural and surface parameters. Results are shown from 1000 regressions.
Donor flammability, DFMC and slope had the effect of reducing flame height predicted using species leaf traits (FSL) relative to predictions based on structure only (FS) (i.e. positive regression coefficients, Table 8). By contrast, wind velocity and overstorey shelter (e.g. delta sum LAI) tended to increase predicted flame height when leaf traits were included (FSL) relative to predictions based on structure only (FS) (i.e. negative regression coefficients, Table 8).
Discussion
Predictions of flame height in eucalypt forests were improved nearly four-fold when effects of plant community structure and species-level leaf traits were considered using the Forest Flammability Model, compared with use of surface-fuel load variables alone. The study therefore showed that floristic composition directly affects fire behaviour at the stand level through vegetation structure and leaf traits. As a result, changes to floristic composition have the potential to alter fire behaviour in a manner that cannot be accounted for by considering only conventional fuel parameters such as surface fuel load or understorey cover.
Variations in surface-fuel variables had almost no explanatory power for larger (≥1 m) flames but produced satisfactory predictions when flames were small. This suggests that surface fuel variables were important drivers of low intensity fires in these forests, but ignition of whole plants was the main cause of larger flames. Flame propagation across gaps between plants to produce larger flames was either limited or enabled by plant size and spacing and the flammability of leaves. Species-specific leaf traits are known to affect surface fire behaviour but models for representing these effects on the mass and packing of surface litter and subsequent fire behaviour are lacking. In addition, leaves of the Eucalyptus spp. that dominated the forests sampled in this study exhibited relatively small morphological variations (authors’ unpublished data). Thus we omitted consideration of effects of leaf variations on the flammability of surface litter in the current study Exploration of such effects will be an important priority for future research.
Our analyses indicated that donor flame height was the most influential mechanism on overall predicted flame height (Fig 7). For example, predicted flame heights produced by FSL were substantially increased, relative to FS, from inclusion of specific leaf traits in sites 1, 56, 45, 6, 91 (Fig 7A). In these sites major increases in donor flame height in the near surface, elevated and midstorey layers were associated with increases in flame height of varying magnitude among a broad range of species (e.g. Acacia, Bossiaea, Bursaria, Carex, Eucalyptus, Leptospermum spp.) resulting from inclusion of specific leaf traits (Fig 7A). By contrast, in sites 94, 48, 23, 115, 56, predicted flame heights from FSL were substantially reduced relative to those predicted from FS (Fig 7B). Large reductions in donor flame height in many of the same species noted above, plus other species (Coprosma, Kunzea, Lomatia, Phebalium), resulting from the inclusion of leaf traits, were associated with the overall reductions in predicted flame height for these sites (Fig 7B). While these results illustrated the strong sensitivity of predicted flame height to leaf traits via changes to donor flammability, such effects were not evident for the same species and others in the majority of sites (S3 Fig). Expression of species-level effects was therefore conditional on other influences, including some of the other structural and environmental effects indicated in the ensemble regression analyses (Table 8).
Showing the five sites where the difference in prediction (FSL—FS) was most positive (a), and the five sites where it was most negative (b). Grey bars show flame heights per species from FSL and horizontal lines show flame heights from mean species in FS. Sites are ordered by delta flame height and the site number is given at the top of each plot.
Effects of overstorey sheltering (i.e. various measures of delta LAI, Table 8) on predicted flame height produced from FSL relative to FS were consistently negative. This indicated that LAI had a stronger dampening effect on overall predicted values when leaf traits were included (FSL) compared with FS. These trends possibly reflected the effects of inclusion of a much wider range of leaf sizes on wind penetration within and beneath plant canopies. Nonetheless the influence of delta LAI effects was small relative to that of donor flame height (Table 8).
As expected from previous analyses of fire severity variations in eucalypt forests [59–61], various aspects of weather influenced the sensitivity of predicted flame height. Wind velocity and DFMC were frequently influential because of the capacity to propagate flames across strata in the former case, and the availability of litter to burn in the latter case. Wind velocity had a positive influence on predicted flame height for FSL, so that flame height predictions were boosted for FSL compared with FS under stronger winds (Table 8). By contrast, increasing DFMC and slope tended to reduce flame height predicted by FSL relative to FS. In the models, DFMC not only affected litter but also suspended dead fuel in surface and elevated strata. DFMC in these above-ground layers may have been more influential on flame propagation when combined with leaf traits such as size and thickness. Effects of slope may have functioned similarly. In both cases, the magnitude of effects was weak (Table 8).
The net effect of the key mechanisms frequently included in the regression models was that a model not using leaf traits (i.e. FS) was more likely to over-predict flame height. This was exacerbated on steep slopes, when dead fuel was moister and when wind speeds were slower. In contrast, a model without leaf traits (FS) tended to under-predict flame height when LAI and consequent wind reduction were overestimated. Leaf traits therefore become increasingly important for predicting large flames under dry, windy conditions, and to a lesser extent, on shallow slopes.
Given the complexity of these mechanisms it was unsurprising that the sites exhibited a wide range of variation in comparative predictions of flame height (i.e. flame height FS–flame height FSL), including under-prediction and over-prediction (Fig 7) and relative insensitivity to the component species mix and abundances. Overall, variations in flame height predicted and validated in this study therefore emerged from complex interactions among species and traits. For example, fine leaves may be characterised as more flammable than thick leaves due to their shorter time to ignition [11,35,62]. The effect of this was apparent at site 45, where the receiver ignitability of fine leaves in the shrub Bossiaea foliosa was greater than the mean values used in FS, so that they ignited and thus propagated flames into the midstorey using FSL (Fig 7). Thinner leaves however also burn for less time [35], so that at site 66, the fine leaves of the shrub Cassinia aculeata reduced flame duration compared with mean values used in FS. This result contributed to a reduction in donor flame height (see S1 Table, S2 Fig and S4 Table for an explanation of the mechanism), preventing the upward propagation of flame into the higher strata using FSL (Fig 7). Thus the expression of basic patterns of flammability associated with particular leaf traits will be contingent on complex interactions among species and their three dimensional arrangement in any community. Higher order consequences of fundamental properties of flammability of particular plant traits may therefore not be readily predictable without detailed knowledge of the make-up of communities. The FFM provides a mechanistic understanding of such complexity and the resultant fire behaviour.
Ecological implications
Our study found that for these sites, plant traits were more important for predicting flame height than was surface fuel load. This has important implications for fire ecology and management as it verifies some elements of a potential feedback relationship between fire behaviour and fire ecology that have been identified in conceptual models e.g. [63,64] but less commonly demonstrated and quantified in field studies. Invasive species such as grasses can alter fire regimes post-invasion [65] and differing plant communities can produce in situ micro-climates that are likely to alter the probability of fire e.g. [66]; however the ability of relatively subtle variations in species composition to differentially influence fire behaviour is less well explored.
While it is well known that particular fire regimes can select for certain species and thus alter the structure and composition of plant communities [67–69]; our study shows that variations in, for example, key shrub and herb species that may ensue from fire regime variations have the potential to influence subsequent fire behaviour via changes in plant traits and stand structure. Such effects on one aspect of fire behaviour—flame height, have the potential to cause flow-on effects on structure and composition via differential effects on survival and recovery of individuals of differing species. For example, the difference between a fire confined to surface or near-surface fuel strata versus a fire in the mid-storey or tree crowns may determine survival or mode of resprouting in various species of shrubs and trees in these forests [70,71]. Flame height will also have other important ecological influences in eucalypt forests, such as the nature of tree injuries and hollow formation, as well as the survival of arboreal mammals e.g. [72].
The results of the study indicate the need to more widely examine implications of forest structure and composition for fire behaviour. For example, the need to explore a wide range of plant trait variation in global vegetation models has been advocated [73,74]. Associated fire behaviour models may need to be structured in a way that more fully allows expression of the influences of such traits on flammability and fire behaviour.
Model limitations and future development
An important limitation of the FFM is the time to ignition sub-model [35], which does not include the possible influences of leaf phosphate content or tissue density [11,13,62]. The FFM framework could incorporate such information once suitable relationships are developed and adequate data for individual species become widely available. The FFM utilises many inputs to describe the structure and leaf traits of component plant species, and this has the potential to constrain its use in operational applications. Nonetheless, the predictive accuracy of the current study was achieved using only some structural and leaf trait data measured in the field. Estimates for missing parameters were derived from photogrammetry, herbarium records and published databases e.g. [47]. Along with these approximations, satellite-based multi and hyper-spectral scanning offer possibilities for measuring parameters such as leaf moisture [63,64,75], particularly when supplemented with easily accessible information on species characteristics [76]. LiDAR has the potential to be used to measure structural parameters e.g. [77,78] and some leaf traits e.g.[79,80]. Further work should quantify model sensitivity to the accuracy of these inputs.
It is likely that a considerable proportion of the error in flame height predictions may be due the use of photogrammetric measurement of fuels, standard live moisture values, and the predictions of weather parameters from weather stations not located on the fire ground. Considerations of such sources of uncertainty and error will be a priority for future studies.
While we have demonstrated that the FFM more accurately predicted variations in flame length as a function of variations in the composition and structure of some eucalypt forests, further work is required to explore the accuracy of the FFM in other communities, and its capacity to represent other aspects of fire behaviour such as rate of spread. Such insights will help to determine the types of applications where the use of a complex model such the FFM may be most viable, given the demands of parametrisation of fuel compared with other fire behaviour models.
Conclusions
The inclusion of species level traits in a modelling framework that represented flammability of plant parts and the resultant capacity of flames to propagate through their three dimensional arrangement (i.e. as represented in the FFM) considerably enhanced the ability to predict flame height in eucalypt forest communities. Conventional approaches to modelling fire behaviour based on the mass of surface litter and simple measures of above-ground fuel strata may therefore be unable to predict aspects of fire behaviour that arise from variations in forest composition. The FFM mechanistically explains the effects of plant traits on fire behaviour, thereby providing a basis for exploring the way in which variations in community composition can affect fire behaviour in the short-term and, perhaps, longer term fire regimes. The structure of the model also allows for the integration of new theoretical work on flammability and fire physics, providing an adaptive framework that can provide both direction for new research and application to evolutionary, ecological and management questions.
Supporting Information
S2 Table. Sub-models used in the FFM for this study.
https://doi.org/10.1371/journal.pone.0160715.s005
(PDF)
S4 Table. Plant and flame traits, and resultant heating as shown in S2 Fig.
https://doi.org/10.1371/journal.pone.0160715.s007
(PDF)
S5 Table. Species leaf traits used in the FFM for this study.
https://doi.org/10.1371/journal.pone.0160715.s008
(PDF)
S6 Table. Observed and predicted flame heights.
https://doi.org/10.1371/journal.pone.0160715.s009
(PDF)
Acknowledgments
The FFM was developed for a PhD by P.Z. as part of the NSW Office of Environment and Heritage, supervised by R.W., M.G., G.C., R.B. and Roger Good with coding and error-checking by Christopher Thomas and M.B. Paul Carmen of the Australian National Botanical Gardens and Donwen Zylstra assisted with field work, and Rick McRae of the ACT Emergency Services Agency provided historical data on weather and fire spread.
This paper is dedicated to the legacy of the late Mr Roger Good, who was an early and powerful advocate of the need to manage fire with consideration to plant flammability and ecology. Roger championed the development of the FFM in its initial stages and largely made it possible within the NSW National Parks and Wildlife Service, then followed this with ongoing support and insight. His love and voice for the mountains and wild places was an inspiration to many, and he will be missed as a powerful voice for conservation, a mentor and a friend.
Author Contributions
- Conceived and designed the experiments: PZ RB RW AMG GC.
- Performed the experiments: PZ.
- Analyzed the data: PZ RB MB TP.
- Contributed reagents/materials/analysis tools: PZ MD.
- Wrote the paper: PZ RB MB TP MD RW AMG GC.
- Designed the software used in analysis: PZ MB.
References
- 1. Olsen TS. Energy storage and the balance of producers and decomposers in ecological systems. Ecology. 1963;44: 322–331.
- 2. Cornwell WK, Elvira A, van Kempen L, van Logtestijn RSP, Aptroot A, Cornelissen JHC. Flammability across the gymnosperm phylogeny: the importance of litter particle size. New Phytol. 2015;206: 672–681. pmid:25675853
- 3. de Magalhães RMQ, Schwilk DW. Leaf traits and litter flammability: evidence for non-additive mixture effects in a temperate forest. J Ecol. 2012;100: 1153–1163.
- 4. Schwilk DW, Caprio AC. Scaling from leaf traits to fire behaviour: community composition predicts fire severity in a temperate forest. J Ecol. 2011;99: 970–980.
- 5. Scarff FR, Westoby M. Leaf litter flammability in some semi-arid Australian woodlands. Funct Ecol. 2006;20: 745–752.
- 6. Cheney NP, Gould JS, McCaw WL, Anderson WR. Predicting fire behaviour in dry eucalypt forest in southern Australia. For Ecol Manage. Elsevier B.V.; 2012;280: 120–131.
- 7. McCaw WL, Gould JS, Cheney NP, Ellis PFM, Anderson WR. Changes in behaviour of fire in dry eucalypt forest as fuel increases with age. For Ecol Manage. Elsevier B.V.; 2012;271: 170–181.
- 8. Taylor C, McCarthy MA, Lindenmayer DB. Nonlinear effects of stand age on fire severity. Conserv Lett. 2014;7: 355–370.
- 9. Cruz MG, Alexander ME, Wakimoto RH. Development and testing of models for predicting crown fire rate of spread in conifer forest stands. Can J For Res. 2005;35: 1626–1639.
- 10. Van Wagner CE. Conditions for the start and spread of crown fire. Can J For Res. NRC Research Press Ottawa, Canada; 1977;7: 23–34.
- 11. Grootemaat S, Wright IJ, van Bodegom PM, Cornelissen JHC, Cornwell WK. Burn or rot: leaf traits explain why flammability and decomposability are decoupled across species. Funct Ecol. 2015; n/a–n/a.
- 12. Madrigal J, Marino E, Guijarro M, Hernando C, Díez C. Evaluation of the flammability of gorse (Ulex europaeus L.) managed by prescribed burning. Ann For Sci. 2011;69: 387–397.
- 13. Scarff FR, Gray BF, Westoby M. Exploring phosphate effects on leaf flammability using a physical chemistry model. Int J Wildl Fire. 2012;21: 1042–1051.
- 14. Plucinski MP, Anderson WR, Bradstock RA, Gill AM. The initiation of fire spread in shrubland fuels recreated in the laboratory. Int J Wildl Fire. 2010;19: 512–520.
- 15. Simpson KJ, Ripley BS, Christin P-A, Belcher CM, Lehmann CER, Thomas GH, et al. Determinants of flammability in savanna grass species. J Ecol. 2015;104: 138–148. pmid:26877549
- 16. Alexander ME, Cruz MG. Assessing the effect of foliar moisture on the spread rate of crown fires. Int J Wildl Fire. 2013;22: 415–427.
- 17. Pérez-Harguindeguy N, Díaz S, Garnier E, Lavorel S, Poorter H, Jaureguiberry P, et al. New handbook for standardised measurement of plant functional traits worldwide. Aust J Bot. 2013;61: 167–234.
- 18. Gill AM, Zylstra P. Flammability of Australian forests. Aust For. 2005;68: 87–93.
- 19. Belcher CM, Mander L, Rein G, Jervis FX, Haworth M, Hesselbo SP, et al. Increased fire activity at the Triassic/Jurassic boundary in Greenland due to climate-driven floral change. Nat Geosci. 2010;3: 426–429.
- 20.
Zylstra P. Live fuel moisture and wildland fire behaviour. In: Viegas DX, editor. Advances in forest fire research. 1st ed. Coimbra, Portugal: Imprensa da Universidade de Coimbra; 2014. pp. 326–335. http://dx.doi.org/10.14195/978-989-26-0884-6
- 21. Fernandes PM, Cruz MG. Plant flammability experiments offer limited insight into vegetation–fire dynamics interactions. New Phytol. 2012;194: 606–609. pmid:22288940
- 22. Schwilk DW. Dimensions of plant flammability. New Phytol. 2015;206: 486–488. pmid:25800615
- 23. Zylstra P. Physical properties determining flammability in sclerophyllous leaves and flame propagation within shrubs and trees in the Australian Alps. For Ecol Manage. 2006;234: S81.
- 24. Anderson HE. Forest fuel ignitability. Fire Technol. 1970;6: 312–322.
- 25. Zylstra P. The historical influence of fire on the flammability of subalpine Snowgum forest and woodland. Vic Nat. 2013;130: 232–239.
- 26. Cassiani M, Katul GG, Albertson JD. The effects of canopy Leaf Area Index on airflow across forest edges: large-eddy simulation and analytical results. Boundary-Layer Meteorol. 2007;126: 433–460.
- 27. Cionco RM. A wind-profile index for canopy flow. Boundary-Layer Meteorol. 1972;3: 255–263. Available: http://link.springer.com.ezproxy.uow.edu.au/article/10.1007/BF02033923
- 28. Matthews S. Dead fuel moisture research: 1991–2012. Int J Wildl Fire. 2014;23: 78–92.
- 29.
Rothermel RC. A mathematical model for predicting fire spread. Research Paper INT-115. Ogden, Utah; 1972.
- 30. Burrows ND. Fire behaviour in Jarrah forest fuels: 2. Field experiments. CALMScience. 1999;3: 57–84.
- 31. Noble IR, Bary GA V, Gill AM. McArthur’s fire-danger meters expressed as equations. Austral Ecol. 1980;5: 201–203.
- 32. Keeley JE. Fire intensity, fire severity and burn severity: a brief review and suggested usage. Int J Wildl Fire. 2009;18: 116–126. Available: http://www.publish.csiro.au/?paper=WF07049
- 33.
Gill AM, Moore PHR. Some ecological research perspectives on the disastrous Sydney fires of January 1994. Proceedings of the 2nd International Conference on Forest Fire Research Volume I. Coimbra; 1994. pp. 63–72. Available: https://publications.csiro.au/rpr/pub?list=BRO&pid=procite:4551cd49-2ed5-44d4-9888-f196a2be8dd5
- 34. Williams RJ, Gill AM, Moore PHR. Seasonal changes in fire behaviour in a tropical savanna in northern Australia. Int J Wildl Fire. 1998;8: 227–239.
- 35.
Zylstra P. Forest Flammability: Modelling and Managing a Complex System [Internet]. University of NSW, Australian Defence Force Academy. 2011. https://doi.org/10.13140/2.1.3722.0166
- 36. Zylstra P. Modelling leaf, plant and stand flammability for ecological and operational decision making. Geophys Res Abstr. 2014;16. Available: http://meetingorganizer.copernicus.org/EGU2014/EGU2014-1688.pdf
- 37. Weber RO, Gill AM, Lyons PRA, Moore PHR, Bradstock RA, Mercer GN. Modelling wildland fire temperatures. CALMScience Suppl. 1995;4: 23–26.
- 38. Burrows ND. Fire behaviour in Jarrah forest fuels: 1. Laboratory experiments. CALMScience. 1999;3: 31–56.
- 39.
Gould JS, McCaw WL, Cheney NP, Ellis PF, Matthews S. Field guide. Fuel assessment and fire behaviour prediction in dry eucalypt forest. Interim. Canberra ACT, Perth WA: Ensis-CSIRO, Department of Environment and Conservation; 2007.
- 40. Van Wagner CE. Height of crown scorch in forest fires. Can J For Res. 1973;3: 373–378.
- 41. Philpot CW. Influence of mineral content on the pyrolysis of plant materials. For Sci. 1970;16: 461–471.
- 42.
Doherty MD, Wright GT. Preliminary assessment of the impact of the 2003 wild fires on the flora and vegetation of Brindabella National Park, Bimberi Nature Reserve and Burrinjuck Nature Reserve. A report prepared for the Reserve Conservation Unit, New South Wales National Parks. Canberra, ACT; 2004.
- 43.
Doherty MD, Wright GT. Continuing assessment of the impact of the 2003 wildfires on the flora and vegetation of Brindabella National Park, Bimberi Nature Reserve and Burrinjuck Nature Reserve, NSW. Second year post-fire sampling report. A report prepared for the Reserve Conserv. Canberra, ACT; 2005.
- 44.
Doherty MD. Vegetation survey and mapping of Brindabella National Park and adjacent vacant Crown Lands. A report prepared for NSW National Parks and Wildlife Service, Queanbeyan District. Canberra, ACT; 1997.
- 45.
Doherty MD. Vegetation survey and mapping of Burrinjuck Nature Reserve. A report prepared for NSW National Parks and Wildlife Service, Queanbeyan District. Canberra, ACT; 1998.
- 46.
Gill AM, Moore PHR. Ignitability of leaves of Australian plants. Contract report to the Australian flora foundation. Canberra, ACT; 1996.
- 47. Kattge J, Díaz S, Lavorel S, Prentice IC, Leadley P, Bönisch G, et al. TRY—a global database of plant traits. Glob Chang Biol. 2011;17: 2905–2935.
- 48. Onoda Y, Westoby M, Adler PB, Choong AMF, Clissold FJ, Cornelissen JHC, et al. Global patterns of leaf mechanical properties. Ecol Lett. 2011;14: 301–312. pmid:21265976
- 49. Dickinson K, Kirkpatrick JB. The flammability and energy content of some important plant species and fuel components in the forests of southeastern Tasmania. J Biogeogr. 1985;12: 121–134. Available: http://www.jstor.org/discover/10.2307/2844836?uid=2&uid=4&sid=21105117046943
- 50.
Butler BW, Forthofer JM, Finney MA, Bradshaw LS, Stratton R. High resolution wind direction and speed information for support of fire operations. In: Aguirre-Bravo C, Pellicane PJ, Burns DP, Draggan S, editors. Monitoring science and Technology symposium: Unifying knowledge for sustainability in the Western Hemisphere RMRS-P-42CD. Fort Collins, CO; 2004. pp. 595–602.
- 51. Sharples JJ, McRae RHD, Weber RO. Wind characteristics over complex terrain with implications for bushfire risk management. Environ Model Softw. Elsevier Ltd; 2010;25: 1099–1120.
- 52. Alexander ME, Cruz MG. Interdependencies between flame length and fireline intensity in predicting crown fire initiation and crown scorch height. Int J Wildl Fire. 2012;21: 95–113.
- 53.
Barrett TW. Modelling burn severity for the 2003 NSW/ACT unplanned fires using Landsat imagery. Bushfire Conference 2006 Life in a fire-prone environment: translating science into practice. Brisbane, Qld; 2006.
- 54. Willmott CJ, Matsuura K. Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance. Clim Res. 2005;30: 79–82.
- 55. Tibshirani R. Regression shrinkage and selection via the Lasso. J R Stat Soc Ser B. 1996;58: 267–288.
- 56. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33: 1–22. Available: http://www.jstatsoft.org/v33/i01/ pmid:20808728
- 57.
R Core Team. R: A language and environment for statistical computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2016. Available: http://www.r-project.org/
- 58. Hintze JL, Nelson RD. Violin Plots: A box plot-density trace synergism. Am Stat. 1998;52: 181–184.
- 59. Bradstock RA, Hammill KA, Collins L, Price OF. Effects of weather, fuel and terrain on fire severity in topographically diverse landscapes of south-eastern Australia. Landsc Ecol. 2010;25: 607–619.
- 60. Price OF, Bradstock RA. The efficacy of fuel treatment in mitigating property loss during wildfires: Insights from analysis of the severity of the catastrophic fires in 2009 in Victoria, Australia. J Environ Manage. 2012;113: 146–157. pmid:23025983
- 61. Collins L, Bradstock RA, Penman TD. Can precipitation influence landscape controls on wildfire severity? A case study within temperate eucalypt forests of south-eastern Australia. Int J Wildl Fire. 2013;
- 62.
Grootemaat S, Zylstra P, van Bodegom PM. Models for leaf ignitability based on leaf traits. Plant traits and their effect on fire and decomposition. North Ryde: Macquarie University; 2015. pp. 130–160. https://doi.org/10.1017/CBO9781107415324.004
- 63. Yebra M, Dennison PE, Chuvieco E, Riaño D, Zylstra P, Hunt ER, et al. A global review of remote sensing of live fuel moisture content for fire danger assessment: Moving towards operational products. Remote Sens Environ. 2013;136: 455–468.
- 64. Chuvieco E, Cocero D, Riaño D, Martin P, Martínez-Vega J, de la Riva J, et al. Combining NDVI and surface temperature for the estimation of live fuel moisture content in forest fire danger rating. Remote Sens Environ. 2004;92: 322–331.
- 65. Bowman DMJS, Perry GLW, Marston JB. Feedbacks and landscape-level vegetation dynamics. Trends Ecol Evol. Elsevier Ltd; 2015;30: 255–260. pmid:25837918
- 66. Little JK, Prior LD, Williamson GJ, Williams SE, Bowman DMJS. Fire weather risk differs across rain forest-savanna boundaries in the humid tropics of north-eastern Australia. Austral Ecol. 2012;37: 915–925.
- 67.
Keeley JE, Bond WJ, Bradstock RA, Pausas JGH, Rundel PW. Fire in Mediterranean ecosystems. Ecology, evolution and management. Cambridge, United Kingdom: Cambridge University Press; 2012.
- 68.
Scott AC, Bowman DMJS, Bond WJ, Pyne SJ, Alexander ME. Fire on earth: an introduction. 1st ed. Chichester: Wiley Blackwell; 2014.
- 69. Gill AM. Fire and the Australian flora. Aust For. 1975;38: 4–25.
- 70.
Gill AM, Catling PC. Fire regimes and biodiversity of forested landscapes of southern Australia. In: Bradstock RA, Williams JE, Gill AM, editors. Flammable Australia: fire regimes and biodiversity of a continent. Melbourne: CSIRO Publishing; 2012. pp. 351–369.
- 71.
Gill AM. Bushfires and biodiversity in southern Australian forests. In: Gill AM, Bradstock RA, Williams RJ, editors. Flammable Australia: fire regimes and biodiversity in a changing world. CSIRO Publishing; 2012. pp. 235–252.
- 72. Lindenmayer DB, Blanchard W, McBurney L, Blair D, Banks SC, Driscoll D, et al. Fire severity and landscape context effects on arboreal marsupials. Biol Conserv. Elsevier Ltd; 2013;167: 137–148.
- 73. van Bodegom PM, Douma JC, Verheijen LM. A fully traits-based approach to modeling global vegetation distribution. Proc Natl Acad Sci U S A. 2014;111: 13733–8. pmid:25225413
- 74. Higgins SI, Langan L, Scheiter S. Progress in DGVMs: A comment on “Impacts of trait variation through observed trait-climate relationships on performance of an Earth system model: A conceptual analysis” by Verheijen et al. (2013). Biogeosciences. 2014;11: 4357–4360.
- 75. Caccamo G, Chisholm LA, Bradstock RA, Puotinen MLA, Pippen BG. Monitoring live fuel moisture content of heathland, shrubland and sclerophyll forest in south-eastern Australia using MODIS data. Int J Wildl Fire. 2012;21: 257–269.
- 76. Yebra M, Chuvieco E. Linking ecological information and radiative transfer models to estimate fuel moisture content in the Mediterranean region of Spain: Solving the ill-posed inverse problem. Remote Sens Environ. 2009;113: 2403–2411.
- 77. Lovell JL, Jupp DLB, Culvenor DS, Coops NC. Using airborne and ground-based ranging lidar to measure canopy structure in Australian forests. Can J Remote Sens. 2003;29: 607–622.
- 78. Andersen H, McGaughey RJ, Reutebuch SE. Estimating forest canopy fuel parameters using LIDAR data. Remote Sens Environ. 2005;94: 441–449.
- 79. Jensen JLR, Humes KS, Vierling LA, Hudak AT. Discrete return lidar-based prediction of leaf area index in two conifer forests. Remote Sens Environ. 2008;112: 3947–3957.
- 80. Korhonen L, Korpela I, Heiskanen J, Maltamo M. Airborne discrete-return LIDAR data in the estimation of vertical canopy cover, angular canopy closure and leaf area index. Remote Sens Environ. Elsevier Inc.; 2011;115: 1065–1080.