New In-Flame Flammability Testing Method Applied to Monitor Seasonal Changes in Live Fuel

Improving the accuracy of fire behavior prediction requires better understanding of live fuel, the dominant component of tree crowns, which dictates the consumption and energy release of the crown fire flame-front. Live fuel flammability is not well represented by existing evaluation methods. High-flammability live fuel, e.g., in conifers, may maintain or increase the energy release of the advancing crown fire flame-front, while low-flammability live fuel, e.g., in boreal deciduous stands, may reduce or eventually suppress flame-front energy release. To better characterize these fuel–flame-front interactions, we propose a method for quantifying flammability as the fuel’s net effect on (contribution to) the frontal flame energy release, in which the frontal flame is simulated using a methane diffusion flame. The fuel’s energy release contribution to the methane flame was measured using oxygen consumption calorimetry as the difference in energy release between the methane flame interacting with live fuel and the methane flame alone. In-flame testing resulted in fuel ignition and consumption comparable to those in wildfires. The energy release contribution of live fuel was significantly lower than its energy content measured using standard methods, suggesting better sensitivity of the proposed metric to water contentand oxygen deficiency-associated energy release reductions within the combustion zone.


Introduction
The efficiency of wildland fire management in protecting values at risk and addressing emerging climate change-related environmental challenges depends on the ability to predict wildfire behavior, which is controlled by the fire environment [1]. Increasing the understanding of the fuel component of the fire environment, in particular, the flammability of live plant material available for combustion (i.e., live fuel) can improve the accuracy of fire behavior predictions [2,3]. As a dominant component of crown fuel consumption, live fuel drives crown fires, which are difficult to predict and control, and which make up the largest part of the area burned in the North American boreal forest [4]. Ever since considered the heat (energy) release rate to be the most important variable in determining fire hazard. The available energy release per unit area within the flaming front (energy release component) is an important output of the NFDRS [38]. The energy release from burning fuel into the unburned fuel ahead (horizontal propagating flux) has been linked with the rate of fire spread, the preheat ignition energy, and the fuel bulk density in the heat balance equation [39]. This relationship is used for calculating the rate of spread in the Rothermel surface fire spread model within the NFDRS [40], predicting the initiation of crowning [31], and modifying crown fire rate of spread [33] within the FBP. A fire-front's energy release rate, or fire intensity, directly affects firefighter safety [41], the probability of initial attack success [42], and the range of adequate strategies and tactics for wildfire control [32]. Fire intensity in Byram's [43] formula is calculated as a product of the dry mass of fuel consumed per unit area in the active flaming zone, the rate of fire spread, and heat of combustion (H) as a measure of fuel flammability.
Heat of combustion has been traditionally measured using oxygen bomb calorimetry as the gross ("high") heat of combustion (H gross , kJ g −1 ). By testing oven-dry plant material in a pure oxygen environment (e.g., [25]), H gross overestimates potential energy release. To evaluate a more realistic "lower" heat of combustion, or "heat yield", H gross requires adjustment for losses and reductions in energy release that occur within real flame-fronts and are related to the significant and variable water content of live plants [43]. For instance, the FMC of white spruce ranges from 75% to 480% [44,45] or even 500% (as measured in this study) of dry mass. The combustion of live fuel occurs simultaneously with the evaporation of water present in substantial quantities [46,47] because high internal leaf pressure during burning allows live plant tissue to retain water within the temperature range of 160 • C to 220 • C, which is substantially higher than the normal boiling point (100 • C) of water [48]. Additionally, the high heating rates of a typical fire often create temperature gradients within the fuel, with surfaces reaching ignition temperatures while water is still being evaporating from the much cooler internal regions [27]. High energy losses for fuel preheating as well as for evaporation of water of the reaction and water contained in the fuel [43] in turn result in a reduction in reaction temperature and energy release. When compared with rehydrated dead foliage of similar moisture content, live foliage reaches a lower temperature during preheating and drying within an incoming flame (175 • C versus 200 • C), exhibits a noticeably slower temperature increase, and takes longer to ignite (10 s versus 5 s) [48]. Further reductions in energy release are caused by the dilution of the gaseous products of pyrolysis and oxygen by water vapor [43,47,49], the oxygen deficiency due to increased oxygen consumption [47,50], and the flow dynamic alteration of interacting flames [50] resulting in an incomplete combustion and a substantial release of unburned hydrocarbons from high-intensity fires [43]. The FBP System does not take all these factors into account, and instead applies a heat of combustion of 18 kJ g −1 [51] estimated as H gross with a single 5-10% deduction for energy lost via the evaporation of only water of the reaction but not water contained in the fuel [52]. This can lead to a substantial error in fire intensity estimation.
The effective heat of combustion (H eff , kJ g −1 ) [18] measured using oxygen consumption calorimetry better accounts for reductions associated with oxygen deficiency and water content by testing fresh plant material in an open-air environment (instead of pure oxygen). It also produces an "effective" value that accounts for incomplete char oxidation, which is observed in real wildfires due to the short duration of flaming combustion during fire front passage. Within the 80-170% range of moisture content typical for live conifers, Babrauskas' method produced values of H eff at 7-12 kJ g −1 (fresh mass basis), or approximately 19 kJ g −1 dry mass basis (at an average of 100% shoot water content) for live fuel [18], which is close to that measured as H gross or assumed by the FBP model at 18 kJ g −1 . The method only partially represents real fires because it utilizes radiative-only heating, while fuels in a wildfire setting are exposed to both radiative and convective heat transfer. Convective heating has been considered [31,53,54] and has been shown to be the dominant energy transfer mechanism in many cases [55], especially in the mid-upper portion of the canopy [56].
Additionally, H eff , as measured by Babrauskas [18], is insensitive to additional reductions in energy release resulting from the interaction of the live fuel flame (reacting flammable gases emitted by the recently ignited live fuel element) with the incoming frontal flame (the combined flame comprised of reacting flammable gases emitted by the already burning fuel elements) (Figure 1, left) within the flames interaction zone (FIZ). This interaction of the flames results in the creation of an oxygen-deficient gaseous mixture caused by the increased oxygen consumption and flow dynamic alteration for the live fuel flame. Therefore, the energy release of the live fuel flame within the FIZ (Figure 1, left) is most likely lower than the energy release of the live fuel flame alone, tested traditionally as H eff using oxygen consumption calorimetry, in which oxygen deficiency is nearly absent because the fuel is surrounded by air due to the use of radiant only heating (Figure 1, right). Moreover, the energy release of the incoming frontal flame itself (Figure 1, left) in the FIZ may also be reduced-both by oxygen deficiency and by the high water content of the live fuel. These reductions are not accounted for by the existing techniques due to the complexity of the multiple factors involved such as fuel water content, dry matter content, the rates of heating, pyrolysis, and water evaporation, as well as the concentrations of oxygen, pyrolysates, and water vapor. Thus, the net change in the energy release of the incoming frontal flame resulting from its interaction with the live fuel element burning within the flame, or the fuel's net contribution to the frontal flame energy release, may be substantially smaller than H gross or H eff . The region where these flames interact-the flames interaction zone (FIZ)-includes energy release reductions that are unaccounted for by current methods. The vertical direction of flame propagation on the diagram, from bottom to top instead of forward-upward as in real crown fire flame-fronts, represents the experimental setup and apparatus. Right: combustion of a live fuel element in traditional tests out of a frontal flame where the live fuel flame is surrounded by atmospheric air as performed, for example, when measuring effective heat of combustion using standard oxygen consumption calorimetry test setup, e.g., [57].
Considering these issues, the main objective of this study was to introduce a new experimental methodology, developed by Melnik et al. [58] and Paskaluk et al. [59], which utilizes in-flame flammability testing (1) to better represent ignition heat transfer within wildfires and (2) to physically represent and account for the additional energy release reductions resulting from the interaction of flames within the FIZ discussed above. Instead of separately estimating the energy release reductions that exist in real flame-fronts and subtracting them from the gross heat of combustion to evaluate "heat yield" [43], the proposed methodology directly measures the heat yield added to the flame-front by the fuel element as the fuel element's energy release contribution to the incoming frontal flame.

Fuel Element's Energy Release Contribution to the Incoming Frontal Flame
The incoming frontal flame was simulated using a non-premixed methane diffusion flame. Although methane is one of the combustible gases released by wildland fuels, methane was used in the tests primarily due to being a readily available laboratory gas with a well-known and consistent composition and heat of combustion. The known flow rate of methane in the tests allows to calculate its energy release potential to verify HRR values measured with oxygen consumption calorimetry. Importantly, using a nonpremixed diffusion methane flame facilitates the creation of exposure conditions similar in temperature and oxygen concentration to those encountered in the diffusion rate-limited wildfire flame with a temperature typically around 1000 • C. The flame was produced by a 10 × 10 cm open burner that provided approximately 40 kW m −2 total heat flux measured at the bottom-center of the sample holder with radiative heating comprising about 35% of this value, which is higher than the 15-20% radiative heat fraction reported in the literature for smaller methane flames [60,61]. The fuel element was represented in the tests by a live fuel sample. Therefore, the fuel element's contribution to the energy release of the incoming frontal flame ∆e + was measured using oxygen consumption calorimetry [57] as the net difference in energy release between the methane flame interacting with the live fuel sample (Figure 2a) and the methane flame alone (Figure 2b), as in Equation (1) where: ∆e + is the fuel element's/sample's contribution to the energy release of the incoming flame (kJ), Q (incoming flame + fuel) is the total energy release of the methane flame interacting with the live fuel sample burning within it (kJ), and Q (incoming flame alone) is the total energy release from the methane flame alone (kJ). In a wildfire, live fuel interacts with the passing flame-front and contributes to its energy release, propagation, and behavior only during the time in which the flame-front is present and the fuel is exposed to it. This flame-front exposure time includes fuel preheating/ignition time and flame-front residence (flaming combustion) time. By analyzing existing literature and fire-front video recordings [56,62,63], the flame-front exposure time is estimated to be on average 61 s, including 29 s of fuel preheating/ignition time and 32 s of flame-front residence time ( Table 1). The relevance of this analysis for the particular fuel type used in the present study (conifers) was confirmed by our preliminary experimental results [59] showing that burning fuel contributes significantly to the heat release rate (HRR) of the incoming methane flame for only a short period (55-65 s). Therefore, to adequately represent the fuel element's contribution to the energy release of the passing flame-front, the duration of 60 s was chosen for the evaluation of the total energy release of the incoming methane flame with the fuel sample burning within the flame (Equation (1)) by integrating its HRR measured with oxygen consumption calorimetry [57]. This 60-s integration window ensured that fuel ignition, flaming combustion, and, partially, char oxidation phases were included in the evaluation. The total energy release of the methane flame alone passing an empty sample holder was determined by integrating its measured HRR over the same time period. Methane-only tests were conducted both at the beginning and at the end of each day of testing, and these two results were averaged to provide a reference HRR over the 60 s period. This measurement was compared to the theoretical value calculated from the measured methane flow rate to confirm the result. In a wildfire, live fuel interacts with the passing flame-front and contributes to its energy release, propagation, and behavior only during the time in which the flame-front is present and the fuel is exposed to it. This flame-front exposure time includes fuel preheating/ignition time and flame-front residence (flaming combustion) time. By analyzing existing literature and fire-front video recordings [56,62,63], the flame-front exposure time is estimated to be on average 61 s, including 29 s of fuel preheating/ignition time and 32 s of flame-front residence time ( Table 1). The relevance of this analysis for the particular fuel type used in the present study (conifers) was confirmed by our preliminary experimental results [59] showing that burning fuel contributes significantly to the heat release rate (HRR) of the incoming methane flame for only a short period (55-65 s). Therefore, to adequately represent the fuel element's contribution to the energy release of the passing flame-front, the duration of 60 s was chosen for the evaluation of the total energy release of the incoming methane flame with the fuel sample burning within the flame (Equation (1)) by integrating its HRR measured with oxygen consumption calorimetry [57]. This 60s integration window ensured that fuel ignition, flaming combustion, and, partially, char The proposed test set-up for measuring the energy release contribution of the sample to the incoming flame using an oxygen consumption cone calorimeter. The energy release of (a) the incoming methane flame interacting with live fuel was greater (even when only judged visually by the volume of flames), compared with (b) the energy release of the methane flame alone. This difference in energy release represents the fuel element's/sample's contribution to the energy release of the incoming methane flame ∆e + (kJ) in Equation (1). From bottom to top in (a): load cell, a methane burner, incoming methane flame, wire-mesh sample holder containing a live fuel sample burning within incoming methane flame, and outgoing flame (methane flame mixed with the flame of the burning live fuel sample). Vertical (upward) direction of flame propagation represented by the experimental setup of the apparatus is slightly different from that occurring in real crown fire flame-fronts, where it is forward-upward. HRR calculations were performed as specified for oxygen consumption cone calorimeter tests in [57], based on O 2 and CO 2 measurements using a Fire Testing Technology (East Grinstead, West Sussex, UK) bench scale mass loss calorimeter instrumented with a Model 300 analyzer from California Analytical Instruments (Orange, CA, USA) with IR CO and CO 2 detectors. The exhaust was sampled at 0.017 L s −1 through the paramagnetic oxygen analyzer. Data were collected at 10 Hz per channel with a USB-2416 data acquisition device (Measurement Computing Corporation, Norton, MA, USA) and a PC using DASYLab 11 data acquisition software (Hoskin Scientific, Burnaby, BC, Canada). The mass loss of the fuel samples during the tests was measured using a 600 g load cell with a custom-made lightweight methane burner mounted on top. Methane flow to the burner at 0.15 L s −1 during testing was controlled with a rotameter to provide a diffusion flame with a nominal heat release rate of 5.5 kW. During each test, after a wire-mesh sample holder containing a pre-weighed sample was placed on the methane burner, data acquisition was initiated, and the methane flow was started and ignited using a gas lighter. Since data acquisition continued for four minutes, the approximately 15 s delay between the start of data acquisition, ignition, and the gas analyzer sampling resulted in 225 s of test data. The measurement uncertainty for standard oxygen consumption calorimetry in an open system, where the products of combustion are diluted with ambient air, includes the uncertainty associated with oxygen concentration measurements (oxygen analyzer accuracy), the assumed heat of combustion (calorimetric coefficient E), measurements of the mass flow rate of the exhaust products, and the assumed combustion expansion factor, which depends on oxygen depletion [64]. The total range of uncertainties in the measured heat release rate could be as high as 20% primarily due to analyzer uncertainty at low oxygen depletion levels, which decreases with moderate oxygen depletion and increases again from 12% to 15% with growth in oxygen depletion, ambient air dilution, and higher contribution of the assumed expansion factor [64]. The uncertainty of CO and soot corrections is lower when the composition of the fuel is known [64] and is expected to be ±5% or lower for conventional organic fuels when the 13.1 kJ/g constant, widely accepted for oxygen-consumption calorimetry, is used [65]. Instead of focusing on measuring the heat release rate directly, this study rated the heat release rate of the combined flame (methane plus forest fuels) relative to a methane flame alone. Since most of the measurement uncertainties are the same for both measurements and are negated in the relative measurement, the uncertainty in this study can be expected to be even smaller than reported by Huggett [65]. More details on the apparatus, procedure, as well as data acquisition and analysis can be found in Paskaluk et al. [59].
The energy release contribution ∆e + in Equation (1) measures the difference in energy release for the incoming flame that resulted from the interaction with the fuel element and, similarly to the effective heat of combustion H eff in [18], considers an "effective" value that accounts for incomplete char oxidation. Therefore, the fuel element's/sample's energy release contribution ∆e + (kJ) in Equation (1) will be referred to as the fuel's differential effective heat of combustion (∂H eff , kJ g −1 ) when expressed either on a mass loss basis or a fresh mass basis. Both of these metrics were compared to determine which one better represented and predicted the flammability of live fuel. However, in the rest of the study, only fresh mass basis ∂H eff was considered and used to investigate the differences in flammability attributable to the age of the plant material, seasonal changes in live fuel flammability, and factors affecting these changes, as well as to evaluate the substantiality of the energy release reductions within the flames interaction zone.

Fuel Samples
Previous studies varied in terms of what was consumed during the flame-front passage in crown fires in coniferous forests-fresh foliage alone or with some fresh branchwood [31], or fresh foliage with fresh branchwood of varying diameters (0-3 mm [66], 0-10 mm [24], and 0-30 mm [8]). The plant material tested in this study consisted of fresh twigs/branchwood 0-9 mm thick with the attached foliage, collectively referred to as shoots. The flow of combustion gases around and through thin, sparsely distributed fuels, such as the live coniferous shoots consumed in a crown fire, significantly differs from the flow above surface fuels such as the layer of needles/foliage on the forest floor. The arrangement of needles removed from branches and laid flat in the sample holder is more representative of surface fuels than fuels consumed in a crown fire, potentially resulting in very different preheating, ignition, and burning behavior. Consequently, it is important to preserve the fuel structure in tests as much as possible, as was achieved in this study by testing the flammability and biophysical properties of the exact same live plant material as is consumed by crown fire-fresh shoots rather than foliage alone. To emphasize the focus on shoots, the terms shoot flammability and shoot water content (SWC) will be used in this study rather than the more traditional terms foliar flammability and foliar moisture content (FMC) or fuel flammability and fuel moisture content.
To adequately represent the spatial structure and flammability at a tree branch scale but to avoid variation in the results due to the irregular spatial distribution of shoots within the branch, fuel samples were standardized for fuel mass, spatial arrangement, and bulk density using the concept of "flat" fuel sample of defined bulk density introduced in this study. A plant canopy is a porous fuel where a fuel element of given mass burns within the average per fuel element combustion space of given volume, which determines fuel porosity and bulk density. The volume of the roughly 30 × 30 × 30 cm combustion/testing space was 0.027 m 3 , which included an approximately 10 × 10 × 30 cm combined burner and sample flame with some surrounding air space since the flame is non-premixed ( Figure 2). To standardize and represent in tests the typical canopy bulk density of full-density conifers at 0.2 kg m −3 [31,66,67], the fuel sample mass within the 0.027 m −3 space should be 0.0054 kg dry mass, or 0.011 kg fresh mass (at an average 100% shoot water content on a dry mass basis; see Section 2.5 below). These 11 g (mean value) samples were further used for flammability testing. To prepare a flat fuel sample of defined bulk density, approximately 9-13 g of shoots were arranged in a single layer (instead of many layers as on an actual tree branch) and placed into a sample basket, ensuring that the spatial arrangement of plant material resembled that in a real tree branch of white spruce ( Figure 3a) and was as consistent as possible from test to test. The design of the wire-mesh sample holder allowed for a constant distance (5 cm) from plant material to the base of the methane flame and its unrestricted flow through the sample (Figure 3b).

Field Sampling
Sampling was performed in a 50 to 70-year-old mixed stand of white spruce (Picea glauca (Moench) Voss) and trembling aspen (Populus tremuloides Michx.) located in the ecological reserve of the University of Alberta Botanic Garden, 15 km SW of Edmonton, Canada. Eighteen white spruce trees 15-20 m tall were selected across the site to represent a variety of local soil moisture conditions. Sampling occurred between 12:00 and 16:00 during 11 sampling days without precipitation or visible moisture on the surface of the plants from May to October 2014. Each sampling day, three to five trees out of the 18 identified were randomly selected and one tree branch from each tree within lower-one-third outer southfacing part of the crown was harvested using a pole pruner ( Figure A1). Tree branches were stored in resealable plastic bags in a refrigerator at 4 • C before flammability tests were performed. Full sampling protocols are described in [7].
Fire 2022, 4, x FOR PEER REVIEW of plant material resembled that in a real tree branch of white spruce ( Figure 3a) and was as consistent as possible from test to test. The design of the wire-mesh sample holder allowed for a constant distance (5 cm) from plant material to the base of the methane flame and its unrestricted flow through the sample (Figure 3b).

Field Sampling
Sampling was performed in a 50 to 70-year-old mixed stand of white spruce (Picea glauca (Moench) Voss) and trembling aspen (Populus tremuloides Michx.) located in the ecological reserve of the University of Alberta Botanic Garden, 15 km SW of Edmonton, Canada. Eighteen white spruce trees 15-20 m tall were selected across the site to represent a variety of local soil moisture conditions. Sampling occurred between 12:00 and 16:00 during 11 sampling days without precipitation or visible moisture on the surface of the plants from May to October 2014. Each sampling day, three to five trees out of the 18 identified were randomly selected and one tree branch from each tree within lower-onethird outer south-facing part of the crown was harvested using a pole pruner ( Figure A1). Tree branches were stored in resealable plastic bags in a refrigerator at 4 °C before flammability tests were performed. Full sampling protocols are described in [7].

Test Sequence
Each fresh tree branch was separated into shoots of different ages, and their respective mass proportions in the branch composition were measured. Four sample types were considered: new shoots (if present, N = 42), 1-year-old shoots (N = 48), 2+ year-old shoots (N = 48), and tree branch (made up of new, 1-year-old and 2+ year-old shoots according to their respective mass proportions in the composition of a given branch, N = 47). For each sample type, one fuel sample was prepared, and its differential effective heat of combustion (∂Heff) was tested (185 fuel samples in total). The remaining shoots of a given age were subsampled to determine water content, dry matter content, and fresh mass basis energy content. For the tree branch sample, these biophysical characteristics were estimated as a weighted average of new, 1 year, and 2+ year shoots from the same branch according to their proportions in the branch composition. Since three to five branches (one from each selected tree) were harvested on each sampling day, three to five individual measurements of ∂Heff and biophysical characteristics were performed for each of the four

Test Sequence
Each fresh tree branch was separated into shoots of different ages, and their respective mass proportions in the branch composition were measured. Four sample types were considered: new shoots (if present, N = 42), 1-year-old shoots (N = 48), 2+ year-old shoots (N = 48), and tree branch (made up of new, 1-year-old and 2+ year-old shoots according to their respective mass proportions in the composition of a given branch, N = 47). For each sample type, one fuel sample was prepared, and its differential effective heat of combustion (∂H eff ) was tested (185 fuel samples in total). The remaining shoots of a given age were subsampled to determine water content, dry matter content, and fresh mass basis energy content. For the tree branch sample, these biophysical characteristics were estimated as a weighted average of new, 1 year, and 2+ year shoots from the same branch according to their proportions in the branch composition. Since three to five branches (one from each selected tree) were harvested on each sampling day, three to five individual measurements of ∂H eff and biophysical characteristics were performed for each of the four sample types for any given sampling day. Daily average results were used for calculating the data points in the seasonal time series.

Biophysical Characteristics
Considering the fact that substantial seasonal variation in dry matter content can lead to the misrepresentation of water content when measured on a dry mass basis [2] (e.g., FMC), shoot water content SWC was calculated on a dry mass basis [68], fresh mass basis, and volume basis [7] (see Nomenclature). Dry matter content was calculated on a fresh mass basis. Gross heat of combustion H gross on a dry mass basis was measured using a model 1341 Plain Jacket Bomb Calorimeter (Parr Instrument Company, Moline, IL, USA) and the standard oxygen bomb calorimetry test method [69]. The H gross , when expressed on a fresh mass basis [19] is referred to in our study as fresh mass basis energy content (EC) because, assuming that water content is an inert diluent [18], this metric represents the theoretical maximum amount of energy that can be released by a unit of live fuel's fresh mass with combustion in pure oxygen after it has been oven dried.

Heat Transfer
In existing oxygen consumption calorimetry methods, a heat flux of 25-50 kW m −2 is within the range observed in wildland fires: 13-140 kW m −2 peak convective and 20-132 kW m −2 peak radiative heat fluxes for surface and mixed (surface/crown) fires and 32-42 kW m −2 peak convective and 120-300 kW m −2 peak radiative heat fluxes in crown fires, with noticeably lower time-averaged values [55] (Table 2). In wildfires, heat transfer is both radiative and convective [55], and the direction of convective heating coincides with the direction of flame propagation (tilted sideways-upward in crown fire flame-fronts). In contrast, in traditional oxygen consumption/cone calorimetry, heat transfer is practically radiative-only. Unlike in real fires, the flame in cone calorimetry tests propagates downward through the fuel sample because energy is emitted by a radiant source above the sample and is directly received only by the upward-facing outer portion of the sample, which ignites first. The opposite direction of the upward flow of hot combustion products from the already burning fuel results in only a slight contact with the unburned fuel in the lower portion of the sample and a negligible element of convective heat transfer. These test conditions result in the partial and variable consumption of fresh plant material due to inconsistent delayed ignition at 52-555 s versus 1-50 s in wildfire flame-fronts [56,63,70]. The prolonged ignition leads to variability in test results largely driven by water evaporation and pyrolysis rather than combustion and, unlike within real flame-fronts [46], substantially reduces fuel water content before ignition, therefore masking water content-related energy release reductions when the fuel finally ignites. In our study, these issues were resolved by using combined radiative and convective heating from the methane flame where the direction of convective flux coincided with the direction of flame propagation (upward: the sample was ignited from below/sides). This is similar to real conditions in crown fires where the direction of heat transfer and flame propagation also coincides (though it is sideways-upward rather than upward, as in the tests). Although the heat flux of 40 kW m −2 we used was comparable to that of existing methods, the changes listed above resulted in rapid and consistent ignition times of 10-30 s and near-complete consumption (on average 87.1%) of tested fresh 0-9 mm thick branchwood with the attached foliage, which closely represents the live fuel consumed within real flame-fronts [8,24,31,66]. Prince and Fletcher [48] achieved a similarly fast (~10 s) and consistent ignition of fresh live leaves by using a similar upward convective heating test setup. Table 2. Convective and radiative heat transfer in wildland fires for different fuel and fire types. It is important to note that convective heat flux is usually inferred from measurements of a total heat flux gauge, the geometry of which is not representative of wildland fuels, so these values must be considered with caution.

Energy Release Reductions
The test method presented in this paper quantifies flammability as the differential effective heat of combustion (∂H eff ), which is an "effective" value that accounts for reduced energy release with incomplete char oxidation during the flame-front passage. Due to the in-flame testing setup, ∂H eff directly accounts for the energy release reductions caused by fuel water content [49] and oxygen deficiency [50] with the interaction of flames in the flames interaction zone as discussed in the last paragraph of the Introduction. The described method produced a considerably lower and broader range of values for live fuel flammability compared with traditional methods, suggesting that the energy release reductions within the flames interaction zone are substantial. The mean ∂H eff for new shoots measured with our method was 0.23 kJ g −1 (Table 3), showing a 97% reduction in energy release compared to the more traditional fresh mass basis energy content (EC) measured here with a bomb calorimeter at 7.55 kJ g −1 . With a 65% reduction in energy release compared to the EC of 9.70 kJ g −1 , the combined mean ∂H eff for all ages of shoots measured with our method was 3.38 kJ g −1 (on a fresh mass basis), or approximately 6.8 kJ g −1 on a dry mass basis (at average 100% shoot water content). In contrast, the FBP model uses a constant of 18 kJg −1 [51] for the "lower" heat of combustion [52], which is almost three times higher and likely substantially over-predicts fire intensity and the resulting spotting distance in live fuels where convective energy is directly calculated [74], while also missing seasonal variation in live fuel conditions. Table 3. Seasonal variation in energy content and flammability. Minimum, maximum, mean, and standard deviations of fresh mass basis energy content (EC, kJ g −1 ) and flammability measured as differential effective heat of combustion (∂H eff , kJ g −1 ). When measured as effective heat of combustion using oxygen consumption calorimetry in the open air with radiant-only heating [18], live fuel flammability ranged from 7 kJ g −1 to 12 kJ g −1 , depending on water content within the 80-170% (dry mass basis) range typical for most live conifers. In contrast, the differential effective heat of combustion (∂H eff ), measured in our study using the same oxygen consumption calorimetry equipment, but with the added in-flame testing setup, showed values for all ages of shoots that were on average lower by 9 kJ g −1 and ranged from a positive 7.10 kJ g −1 to a negative −6.33 kJ g −1 , depending on water content. For new shoots, ∂H eff similarly varied from a positive 6.48 kJ g −1 to a negative value of −6.33 kJ g −1 . In some cases, the ∂H eff was negative for the whole tree branch (Table 3). Since the ∂H eff represents the energy release contribution of the fuel to the incoming flame, its negative values indicated a reduction in the energy release of the incoming methane flame resulting from the interaction with the live fuel sample of high water content and the associated substantial energy release reductions within the flames interaction zone. New shoots had substantial negative ∂H eff at the beginning of the season in Figure 4a and suppressed the energy release of the methane flame (Figure 5a), in contrast to the 1-year-old shoots (Figure 5b). Traditional measurements of energy content using oxygen bomb or radiant heating oxygen consumption calorimetry cannot be negative because they represent the fuel's potential energy release and are insensitive to energy release reductions within the flames interaction zone.   ). Therefore, the live fuel sample's contribution to the energy release of the incoming methane flame expressed as ∂Heff was negative. In the case of (b) highly flammable 1 year-old shoots, the ∂Heff was positive, where the combined volume (and hence energy release) of the incoming methane flame interacting with burning live fuel was larger compared to that of the incoming methane flame alone (indicated by white dashed line).

Flammability Definition and Numerical Fuel Classification
Traditionally, flammability is always a positive quantity because it is defined as the fuel's ability to burn as represented by the ease/time of ignition (ignitibility), as well as combustion rate (combustibility), duration (sustainability) [11], and completeness (consumability) [12]. As a contribution to this broad mostly time/mass-based set, we introduce an energy release-based criterion. Flammability in our study is defined as the ability of a fuel or material to sustain flame propagation, or a fuel element's energy release contribution to the incoming flame ∆e + expressed on a mass loss basis or fresh mass basis as the differential effective heat of combustion (∂Heff). Therefore, the observed variation from a positive value of 7.10 kJ g −1 to a negative −6.33 kJ g −1 in ∂Heff clearly indicates that the contribution of the burning live fuel element to the incoming flame energy release can vary from highpositive to low-negative. The sensitivity of the ∂Heff to these positive or negative effects allows for the development of a numerical classification of materials and substances. Rather than arbitrarily classifying them into fuels, non-fuels, and suppressants, their flammability can be directly measured using ∂Heff as the positive, neutral, or negative value of their contribution to the energy release of the incoming flame. This is especially important for evaluating suppressants and fire chemicals as well as fuel-to-suppressant transitioning materials such as live fuel. Live plant tissue substantially changes the proportions of "combustibles" (dry matter) and "suppressants" (water) in its composition during the  ). Therefore, the live fuel sample's contribution to the energy release of the incoming methane flame expressed as ∂H eff was negative. In the case of (b) highly flammable 1 year-old shoots, the ∂H eff was positive, where the combined volume (and hence energy release) of the incoming methane flame interacting with burning live fuel was larger compared to that of the incoming methane flame alone (indicated by white dashed line).

Flammability Definition and Numerical Fuel Classification
Traditionally, flammability is always a positive quantity because it is defined as the fuel's ability to burn as represented by the ease/time of ignition (ignitibility), as well as combustion rate (combustibility), duration (sustainability) [11], and completeness (consumability) [12]. As a contribution to this broad mostly time/mass-based set, we introduce an energy release-based criterion. Flammability in our study is defined as the ability of a fuel or material to sustain flame propagation, or a fuel element's energy release contribution to the incoming flame ∆e + expressed on a mass loss basis or fresh mass basis as the differential effective heat of combustion (∂H eff ). Therefore, the observed variation from a positive value of 7.10 kJ g −1 to a negative −6.33 kJ g −1 in ∂H eff clearly indicates that the contribution of the burning live fuel element to the incoming flame energy release can vary from high-positive to low-negative. The sensitivity of the ∂H eff to these positive or negative effects allows for the development of a numerical classification of materials and substances. Rather than arbitrarily classifying them into fuels, non-fuels, and suppressants, their flammability can be directly measured using ∂H eff as the positive, neutral, or negative value of their contribution to the energy release of the incoming flame. This is especially important for evaluating suppressants and fire chemicals as well as fuel-to-suppressant transitioning materials such as live fuel. Live plant tissue substantially changes the proportions of "combustibles" (dry matter) and "suppressants" (water) in its composition during the season depending on the phenophase and the level of physiological drought. During June, new shoots of white spruce showed the highest seasonal water content and the lowest fresh mass basis energy content ( Figure 6) resulting in negative values of ∂H eff (Figure 4a) and actually suppressing the energy release of the incoming methane flame (Figure 5a), in contrast to late summer, when the new shoots' flammability is similar to that of older growth (Figure 4a).

Energy Balance
Through in-flame testing, the fuel's energy release contribution expressed as ∆e + (per fuel element) and the differential effective heat of combustion ∂H eff (per unit of fuel element's fresh mass) better represent the processes and conditions within a flame-front including fuel ignition and the interaction of flames within the flames interaction zone. By measuring the fuel element's contribution to the energy release of the incoming flame, ∆e + directly quantifies the gain or reduction in energy release at a given fuel element, which may or may not be sufficient to compensate for the energy losses from that fuel element (∆e − ) into the environment and into the horizontal propagation flux for preheating the next fuel elements. Higher, similar, or lower values of ∆e + relative to /∆e − / indicate increases, no effect, or declines in the horizontal propagation flux for the preheating of the next fuel elements and, hence, the growth, steady propagation, or decline of the incoming flame (see Figure 7 for details). Therefore, ∆e + and ∂H eff more accurately represent the energy generation component of the energy balance of the incoming flame at a fuel element scale and can be used as a more adequate flammability input for flame propagation and fire behavior modelling based on energy balance rather than FMC or time-to-ignition. may or may not be sufficient to compensate for the energy losses from that fuel element (∆e − ) into the environment and into the horizontal propagation flux for preheating the next fuel elements. Higher, similar, or lower values of ∆e + relative to /∆e − / indicate increases, no effect, or declines in the horizontal propagation flux for the preheating of the next fuel elements and, hence, the growth, steady propagation, or decline of the incoming flame (see Figure 7 for details). Therefore, ∆e + and ∂Heff more accurately represent the energy generation component of the energy balance of the incoming flame at a fuel element scale and can be used as a more adequate flammability input for flame propagation and fire behavior modelling based on energy balance rather than FMC or time-to-ignition. The characteristics of the spatial structure of live fuel can alter the complex boundary layer flow of hot combustion gases around and through thin fuels such as fresh live shoots, thus affecting the heat transfer coefficients from combustion gases to the fuel element. This can shift the energy balance (∆E , shown in Figure 7, by affecting energy generation defined in our study as the fuel's energy release contribution ∆e + as well as energy losses from a fuel element ∆e − and the proportions of its two components-losses to the environment and energy used for preheating the next fuel elements. The positive or negative shift in energy balance will affect the propagation of flame from one fuel element to the next and the resulting fire behavior. In addition, live fuels "burst" and shoot jets of gases [48] and burning needles (observed in our study) during combustion due to high , which is the sum of energy generation (∆e + ) and energy losses (∆e − ). Flame propagates from fuel element F1 to fuel element F4; flame from each previous element represents incoming frontal flame. Vertical direction of flame propagation, from the bottom to the top (instead of tilted sideways-upward as in real crown fire flame-fronts) represents the experimental setup and the apparatus. Depending on the weather conditions, and the physical, chemical, and spatial properties of the particular fuel bed, the value of ∆e + may or may not be sufficient to compensate for the energy losses from a fuel element ∆e − to the environment and into the horizontal propagation flux for preheating the next fuel elements. The frontal flame propagates steadily (middle image, equilibrium state) when ∆e + = |∆e − | because the horizontal propagation flux for the preheating of the next fuel elements (which is the "useful" part of ∆e − ) is compensated by the sufficient part of energy generation ∆e + . The frontal flame declines (left image) if ∆e + < |∆e − | because lower values of ∆e + relative to |∆e − | indicate declines in the horizontal propagation flux for the preheating of the next fuel elements, which is now insufficiently compensated by ∆e + . The frontal flame grows (right image) when ∆e + > |∆e − | because higher values of ∆e + relative to |∆e − | indicate increases in the horizontal propagation flux for the preheating of the next fuel elements.
The characteristics of the spatial structure of live fuel can alter the complex boundary layer flow of hot combustion gases around and through thin fuels such as fresh live shoots, thus affecting the heat transfer coefficients from combustion gases to the fuel element. This can shift the energy balance (∆E, shown in Figure 7, by affecting energy generation defined in our study as the fuel's energy release contribution ∆e + as well as energy losses from a fuel element ∆e − and the proportions of its two components-losses to the environment and energy used for preheating the next fuel elements. The positive or negative shift in energy balance will affect the propagation of flame from one fuel element to the next and the resulting fire behavior. In addition, live fuels "burst" and shoot jets of gases [48] and burning needles (observed in our study) during combustion due to high internal leaf pressures [48], which potentially also changes the boundary layer flow and may or may not contribute to the ignition of the neighbouring fuel elements and flame propagation. Since no in-depth analysis of boundary layer fluid motion was undertaken, and consequently the effects of fuel properties on heat transfer coefficients are not known, live fuel flammability testing should be phenology-and species-specific with a special attention to preserving the spatial structure of the fuel.

Stand-Scale Flammability
Although energy release is directly related to the fuel mass loss [11] and, theoretically, the traditional mass loss basis approach should have an obvious advantage, the fresh mass basis approach introduced in this study was equally successful in predicting variation in flammability measured as ∂H eff (Table 4).Therefore, the species-specific ∂H eff for live fuel can be predicted at the forest stand scale using remote sensing-derived predictor variables such as shoot water content and others in Table 4 (see also Figure 8). With further research on the effects of heat transfer intensity, this will allow for operationally predicting the potential energy release of live fuel for the forest stand. It can be calculated as the fresh mass of live fuel in the forest stand available for high-intensity crown fire (typically fresh 0-9 mm thick branchwood with the attached foliage [8,24,31,66]) multiplied by its potential energy output-the fresh mass basis ∂H eff of the same live plant material determined using our method. The amount of live fuel available for crown fire can be measured using standard fuel inventory protocols. This approach, when applied for live and dead fuel, allows for the operational calculations of a maximum possible energy release under extreme fire-weather conditions or the potential net heat content (PNHC) of the forest stand. As a numerical measure of the potential forest stand flammability, the PNHC can be further used in the development of a new numerical stand characteristics-based fuel classification within a new generation of crown fire models. The PNHC, when reduced from potential to actual value depending on the severity of fire-weather conditions, represents the actual net heat content (ANHC) of the forest stand that can be further used as a numerical input of the actual forest stand flammability for energy release-based fire behavior modelling. Table 4. Adjusted R-squared values for the predictor variables in modelling flammability as differential effective heat of combustion on a fresh mass basis (∂H eff ) using traditional and proposed approaches. The proposed fresh mass basis approach introduced in this study showed same or better results in predicting flammability compared with the traditional mass loss basis approach.

Seasonal Variation and Drivers of Flammability
The seasonal trend of live fuel flammability for white spruce observed in 2014 differs substantially from that assumed by the FBP model ( Figure 9) and better matches the historical seasonality of extreme wildfire in Canada (see Figure A2 and data set in Table A1). According to the FBP, extreme crown fire behavior can be expected around 1 June, during the "spring dip", when the FMC is assumed to be the lowest [75] and the corresponding live fuel flammability represented by FMC-derived Crown Spread Factor is the highest [33,76]. However, most extreme wildfires in Canada since 1825 started either substantially earlier (early April to late May) or later (July-August and mid-fall). In this study, the first seasonal peak in live fuel flammability was observed in early May, three weeks earlier than was predicted by the FBP (Figure 9 Figure A2). In contrast, at this time of the season, the FBP predicts the lowest seasonal values of live fuel flammability represented by the Crown Spread Factor [33,76] for conifer stands common in Canada. Figure 9. Seasonal changes in live fuel flammability as measured in our study as differential effective heat of combustion (∂H eff ) for 2014 (solid line) and as assumed by the FBP model when expressed as an FMC-derived Crown Spread Factor (dashed line). The FPB model assumes only one seasonal maximum in live fuel flammability around 1 June. Flammability measured in this study indicates the first seasonal maximum three weeks earlier, in early May, which closely matches the historical seasonality of extreme crown fire behavior in Canada ( Figure A2 and Table A1). In agreement with the historical seasonality of extreme crown fire behavior, flammability measured in this study indicates the second and the third seasonal maximums around 1 July and 1 August when the FBP model assumes lowest values of the season.
Shoot age had a significant effect on live fuel flammability (ANOVA, p < 0.001, F = 60.081, n = 42). New shoots played an important role at the beginning of the season. Their flammability was on average lower (0.23 kJ g −1 against 4.75 kJ g −1 or 4.76 kJ g −1 for 1-year and 2+year old shoots respectively (Table 3)) and varied more widely compared with older growth (Figure 4). The timing and magnitude of the observed seasonal maximums in live fuel flammability for a tree branch (Figure 4a) were best explained by the opposite seasonal trend of the shoot water content volume basis (SWC vol , Figure 6a). The observed "early-August dip" in SWC vol and the simultaneous resulting spike in flammability were likely caused by a summer-fall drought [77]. The first seasonal maximum in flammability, observed in early May, was less accurately compared with SWC vol , indicated by the corresponding minimum in the traditional shoot water content, on a dry mass basis, (SWC dm ) only in the end of June, which is almost three weeks later (Figure 6b). This suggests that SWC dm and its analog FMC alone cannot fully represent the flammability of live fuel. Moreover, since substantial seasonal variation in dry matter content is a major issue in measuring two-variable water content on a dry mass basis [2], such as SWC dm and FMC, the use of single-variable shoot water content on a fresh mass basis (SWC fm ) or SWC vol may be advantageous.
As in previous studies [17,18,29,78], the flammability of live fuel was strongly inversely related to water content (SWC fm in Figure 8a and more traditional SWC dm in Figure 8b). The differential effective heat of combustion was negative for new shoots with SWC dm over 210%. Flammability was strongly directly related to dry matter content (Figure 8c). Confirming the findings of [19], the traditional gross heat of combustion (dry mass basis, H gross ) was unable to satisfactorily explain variation in live fuel flammability (adjusted R 2 = −0.005 in Table 4). In contrast to their results, a non-standard fresh mass basis energy content (Figure 8d), measured in this study as H gross on a fresh mass basis, was as successful in explaining the variation in flammability as water or dry matter content. Since fresh mass basis energy content is determined by both chemical composition and water content, this also supports the conclusions of [79] concerning the importance of these two variables in predicting live fuel flammability.

Limitations and Future Research
To improve the understanding of the effects of canopy spatial structure and fireweather conditions on wildfire behavior, the proposed method requires further exploration of the effects of the amount, arrangement, and bulk density of the tested plant material, the intensity and duration of the methane flame exposure, and the distance from the flame base to the sample. The oxygen consumption calorimetry method [57], which was used as a part of the experimental methodology for measuring differential effective heat of combustion in our study, is insensitive to direct energy losses with fuel preheating and water desorption and evaporation (latent heat). These losses need to be accounted for in further studies. The substantial differences between the seasonal pattern of live fuel flammability assumed by the FBP model and that measured in 2014 suggest the necessity of further investigations over multiple seasons. Different regions, species, and age-classes should also be represented. The water content for some samples taken in May and early June was likely underestimated due to prolonged storage; close-to-real-time testing will improve the representation of seasonal changes in water content and flammability. A greater ability to explain seasonal changes in flammability and the higher sensitivity to drought of the shoot water content volume-basis metric, as compared with more traditional shoot water content metrics (on a dry mass basis), suggests the necessity of further studies on quantifying the flammability of live fuel using a volumetric approach.

Conclusions
The present study was the first to use in-flame flammability testing for quantifying energy release; previously, in-flame testing was only used for quantifying time-to-ignition, e.g., [46], and for studying increased oxygen consumption and flow dynamic alteration within the flames interaction zone of burning fuel elements [50]. An in-flame test setup with upward convective heating similar to that in our study was also used by Borujerdi et al. and Prince and Fletcher [47,48] for testing live leaves; however, only combustion temperature was monitored, rather than energy release measured in the current study. Determining energy release in conditions similar to those within a flame-front, i.e., directly in the flame, allows for more realistic conditions of heat transfer, ignition, and combustion. The samples tested were representative of live fuel consumed by crown fire flame-front, and consisted of fresh branchwood 0-9 mm thick with attached foliage. Fast and consistent ignition and almost complete consumption of tested fuel reinforces the validity of the method. By using in-flame testing, the experimental methodology documented here directly accounts for the additional water content-and oxygen deficiency-associated energy release reductions caused by the interaction of the flames. The values of live fuel flammability measured in our study were almost three times lower compared with those currently used in the FBP System and on average 9 kJ g −1 lower than the values measured traditionally, suggesting an important effect of the energy release reductions within the flames interaction zone.
The observed seasonal trend of live fuel flammability for white spruce in 2014 substantially differs from that assumed by the FBP model and better matches the historical seasonality of extreme wildfire in Canada. At the tree branch-scale, changes in live fuel flammability were dictated by phenology-associated changes in the relative amount and flammability of new shoots during spring and by drought-induced changes in flammability of all ages of shoots throughout the season. Variation in live fuel flammability was equally well explained using water content, dry matter content, and fresh mass basis energy content (the latter is not typically used in wildfire applications). Similar models developed for main forest species should provide stand-specific input of live fuel flammability that can be directly linked with the existing FBP modules as a replacement of the fixed seasonal pattern of variation in FMC.
Using differential effective heat of combustion, flammability in this study was quantified as the fuel's net contribution to the energy release of the incoming flame, that showed both positive and negative values. Therefore, rather than arbitrarily classifying materials and substances into fuels, non-fuels, or suppressants, their flammability can be directly measured using the proposed method as a positive, neutral, or negative energy release contribution to the incoming flame. This is especially important for characterization of suppressants, fire chemicals, and fuel-to-suppressant transitioning materials such as live fuel.
With in-flame testing, our method more accurately measures the energy-generation component of energy balance at a fuel element scale because it better represents the processes and conditions within real flame-fronts and directly quantifies changes in the energy release of the incoming flame. For instance, low and negative values for live crown fuels in leafed-out deciduous and mixedwood boreal stands will indicate a reduction in fire intensity and the eventual suppression of the incoming crown fire flame-front. High values for live crown fuels in coniferous stands (e.g., during drought) will indicate growth in the intensity of the incoming crown fire flame-front, while low and negative values will suggest crown fire weakening.
Successful modelling of energy release on a fresh mass basis instead of a traditional mass loss basis, with further research, will allow for the operational prediction of the potential energy release of a whole forest stand as a measure of its flammability. This variable, determined for extreme fire-weather conditions, can be used in the development of a new numerical stand characteristics-based fuel classification and-when reduced for the actual fire-weather conditions-can be used in energy balance-based fire behavior modelling. Using a more adequate value to represent the flammability of live fuel and forest stand will contribute to improving the accuracy of fire behavior predictions and increasing the efficiency of forest and wildfire management in the face of increasingly complex environmental challenges arising from changes in climate and fire regimes.   Backcountry [95], 2017 Kenow [96], and 2017 Verdant Creek wildfires [97], as well as the 2010 British Columbia, 2016 Peace Region BC, and 2017 British Columbia [86] wildfires. Figure A2. Historical seasonality of extreme crown fire behavior in Canada for 1825-2017 compared with seasonal trend of live fuel flammability measured in 2014 and that assumed by the FBP model. The FPB model represents live fuel flammability expressed as a FMC-derived Crown Spread Factor (dashed line) at a rough temporal scale assuming only one seasonal maximum in live fuel flammability around 1 June. Live fuel flammability measured in our study as differential effective heat of combustion ∂H eff for 2014 (solid line) showed seasonal trend of higher temporal resolution with three seasonal maximums that match the maximums in historical seasonality of extreme crown fire behavior (gray histogram representing number of extreme fires started during a given 5-day period) according to the dataset in Table A1: early May (three weeks earlier compared with the assumptions of the FBP), early July, and early August. By contrast, the Crown Spread Factor by the FBP model assumes close to the lowest values for the season in early July and early August. Table A1. List of extreme wildfires in Canada for 1825-2017. Data from the list were used to build a histogram of historical seasonality of extreme crown fire behavior in Canada for 1825-2017 in Figure A2. When the information differed between two referenced sources, data from both were reported, separated by a slash (/).