Fire as a removal mechanism of pyrogenic carbon from the environment: effects of fire and pyrogenic carbon characteristics

ORIGINAL RESEARCH published: 10 October 2018 doi: 10.3389/feart.2018.00127 Frontiers in Earth Science | www.frontiersin.org 1October 2018 | Volume 6 | Article 127 Edited by: Satish Myneni, Princeton University, United States Reviewed by: Ronny Lauerwald, Free University of Brussels, Belgium Dimitris Poursanidis, Foundation for Research and Technology Hellas, Greece *Correspondence: Stefan H. Doerr [email protected] Specialty section: This article was submitted to Biogeoscience, a section of the journal Frontiers in Earth Science Received: 15 December 2017 Accepted: 10 August 2018 Published: 10 October 2018 Citation: Doerr SH, Santín C, Merino A, Belcher CM and Baxter G (2018) Fire as a Removal Mechanism of Pyrogenic Carbon From the Environment: Effects of Fire and Pyrogenic Carbon Characteristics. Front. Earth Sci. 6:127. doi: 10.3389/feart.2018.00127 Fire as a Removal Mechanism of Pyrogenic Carbon From the Environment: Effects of Fire and Pyrogenic Carbon Characteristics Stefan H. Doerr1*, Cristina Santín1,2, Agustín Merino3, Claire M. Belcher4and Greg Baxter5 1Department of Geography, College of Science, Swansea University, Swansea, United Kingdom, 2Department of Biosciences, College of Science, Swansea University, Swansea, United Kingdom, 3Department of Soil Science and Agricultural Chemistry, University of Santiago de Compostela, Lugo, Spain, 4Wildfire Lab, Hatherly Laboratories, University of Exeter, Exeter, United Kingdom, 5Wildfire Operations Research, FPInnovations, Hinton, AB, Canada Pyrogenic carbon (PyC, charcoal) is produced during vegetation fires at a rate of ∼116–385 Tg C yr−1globally. It represents one of the most degradation-resistant organic carbon pools, but its long-term fate and the processes leading to its degradation remain subject of debate. A frequently highlighted potential loss mechanism of PyC is its consumption in subsequent fires. However, only three studies to date have tested this hypothesis with reported losses of <8–37%, with the effects of PyC chemical characteristics and fire conditions on PyC loss in wildfires remaining unexplored. To address this, we placed materials with different degrees of thermal and chemical recalcitrance (A: wildfire charcoal, B: slash-pile charcoal, C: pine wood and D: cedar wood) on the ground surface just prior to a high-intensity and a low-intensity boreal forest wildfire. Mass losses were highly variable and dependent on fire- and sample characteristics. Mass losses across both fires (as % of dry weight) were for A: 66.5 ±25.2, B: 41.7 ±27.2, C: 78.2 ±14.9, and D: 83.8 ±18.9. Mass loss correlated significantly with maximum temperature (Tmax) recorded on sample surfaces (r=0.65, p=0.01), but only weakly (r=0.33) with time >300◦C. Mass losses also showed a significant negative correlation (r= −0.38, p=0.05) with thermal recalcitrance (T50) determined using Differential Scanning Calorimetry (DSC) and Tmax with charcoal reflectance (Ro) determined after the fires (r=0.46, p=0.05). Losses in the high-intensity fire were significantly higher (p=0.05) than in the low-intensity fire, but the latter had a higher rate of conversion of fuel to PyC. Our results demonstrate that exposure to fire can indeed be a significant removal mechanism for PyC that remains exposed on the ground after a previous fire. The losses found, however, are likely to represent an extreme upper range as most PyC produced in a fire would not remain exposed on the ground surface by the time the next fire occurs. Our data also demonstrate, for real wildfire conditions, the (i) contrasting resistance of different PyC types to combustion and (ii) contrasting net PyC losses between different fire intensities. The DSC and reflectance (Ro) results support the usefulness of these analyses in reflecting thermal degradation resistance and temperature exposure under actual wildfire conditions. Keywords: black carbon, boreal forest, carbon balance, charcoal, charcoal reflectance, wildfire, management burn, thermal analysis Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon INTRODUCTION Pyrogenic carbon (PyC, also known as charcoal, pyrogenic organic matter or black carbon) is a ubiquitous organic residue resulting from the incomplete combustion of organic fuels during vegetation fires (both natural and anthropogenic fires). It comprises a continuum of pyrogenic organic compounds ranging from partly charred material to soot (Santín et al., 2016a). Its chemical properties result in some of the PyC materials being highly resistant to biological breakdown, making PyC one of Earth’s most degradation-resistant organic carbon pools. Estimates of its mean residence time range from decades or centuries (e.g., Steinbeiss et al., 2009; Bird et al., 2015) to millennia (e.g., Thevenon et al., 2010; de Lafontaine et al., 2011), and, importantly, consistently one or two orders of magnitude longer than those of their unburnt precursors (Bruun and EL- Zehery, 2012; Santos et al., 2012; Naisse et al., 2014). The most recent estimates of PyC global production range from between 63 and 140 Tg C yr−1(Bird et al., 2015) and 116–385 Tg C yr−1(Santín et al., 2016a), which is equivalent to ∼0.1–0.6% of the annual terrestrial net primary productivity (Huston and Wolverton, 2009). Production of PyC and its subsequent fate, therefore, have important implications for the global carbon cycle (Landry and Matthews, 2017). However, its long-term fate in the environment and the specific biotic and abiotic processes leading to its degradation remain the subject of much debate (Bird et al., 2015; Santín et al., 2016a). A potentially major mechanism highlighted in the literature for the abiotic terrestrial breakdown of PyC has been its in situ combustion in subsequent fires (e.g., Ohlson and Tryterud, 2000; Czimczik et al., 2005; Preston and Schmidt, 2006; Czimczik and Masiello, 2007; Kane et al., 2010; Foereid et al., 2011). PyC produced in a fire and not transported off site by wind or water, or protected through its transfer to sufficient depths below litter or soil surface layers, can act as fuel and be subject to combustion in subsequent fires. This abiotic loss mechanism for PyC could be particularly important in environments with short fire-return intervals, low bioturbation and topography, which limit, respectively, its vertical movement and off-site transport. The efficacy of this proposed loss mechanism, however, has remained largely untested and was based, until recently, on little more than speculation. Indeed, to our knowledge, this “fire loss” hypothesis has, to date, been tested empirically in only the following three studies. In a high-intensity experimental boreal forest fire, Santín et al. (2013) found median losses of charcoal pieces (2–4 cm size) placed in the forest floor layer (∼2 cm depth) to be <15%. Saiz et al. (2014) reported mean losses of <8% for charcoal pieces (0.5–1.5 cm size) placed on the soil surface in a prescribed fire in an open savannah woodland. Using coniferous fuels typical of the western US, Tinkham et al. (2016) found that repeated burning of litter beds in the laboratory reduced the PyC mass (>6 mm size) formed in the first fire, by 37% in the subsequent fire. Given that vegetation fires as well as PyC materials vary widely in their characteristics, the current study aimed to explore the role of fire and PyC characteristics as potentially important variables in determining PyC removal rates during fires. We hypothesized that PyC with different degrees of carbonization exposed to fires of contrasting intensities would result in different PyC losses during burning. Thus, two different types of PyC (wildfire and slash-pile charcoal from jack pine wood) and, for comparison, uncharred wood (jack pine and western red cedar), were exposed to two contrasting types of wildfires: a high-intensity crown fire and a low-intensity surface fire. Temperature exposure of the samples during the fires was determined using thermocouples attached to sample surfaces, and the chemical and physical characteristics of the exposed and, where relevant, reference (unburnt) samples determined using thermogravimetric-differential scanning calorimetry, elemental- and charcoal surface reflectance analyses. MATERIALS AND METHODS Study Area and Experimental Setup Field experiments were carried out in the Great Slave Lake region of the Canadian Northwest Territories in summer 2015. This boreal region has a dry, sub-humid continental climate with an annual precipitation of ∼300 mm, a wildfire season usually between May and September, and a fire recurrence interval in the region of 150 to 200 years. The forest is dominated by black spruce (Picea mariana) and jack pine (Pinus banksiana) with often thick organic soil layers. The topography is largely flat with an elevation between 100 and 200 m.a.s.l. (Alexander et al., 2004; Stocks et al., 2004). Two types of PyC and two types of uncharred wood materials were used in the fire exposure experiments (Figures 1A,B): (A) Wildfire charcoal generated during an experimental wildfire (June, 2012) in a mature jack pine stand in this region and collected in June 2013 (61◦34’55” N; 117◦11’55” W, see Santín et al. (2015), for site and fire characteristics). (B) Charcoal from a jack pine slash-pile (a pile of slash generated in clearing/harvesting) burn in the winter of 2011/12 and collected in June 2012 at a site adjacent to where wildfire charcoal was sampled. Due to its production in a longburning slash pile fire, this type of charcoal was speculated in previous work (Santín et al., 2013) to exhibit a higher degree of charring and therefore greater resistance to combustion when exposed to a subsequent fire than wildfire charcoal. (C) Unburned jack pine down wood (small dead branches) collected in 2012 from an unburned area in the same study region. (D) Unburned blocks cut from western red cedar (Thuja plicata) wood were included for comparison as a low-density wood used in a complementary study focusing specifically on charcoal reflectance (Belcher et al., 2018). A–C pieces were 2–3 cm in diameter and 2–4 cm in length. D blocks were 3 ×3×3 cm. All were free of bark. Sets were made with one piece of each type (A+B+C+D), and the four pieces were linked together by thin wire and with a metal label added to facilitate relocation after fire (Figure 1B). A total of 30 sets were placed onto the forest floor <24 h prior to exposure to two types Frontiers in Earth Science | www.frontiersin.org 2October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon FIGURE 1 | (A) A sample set placed on the forest floor prior to the high-intensity fire. The samples are connected by a thin wire to facilitate relocation. The thicker wires are thermocouples. (N.B. The pine wood sample appears dark, but was free of bark. The wildfire charcoal is not visible and was located just below the image frame.) (B) Examples of charred samples after the high-intensity fire. (C) The high intensity fire. (D) Sampling after the low-intensity fire where much of the standing vegetation remained unaffected. The figure is published with the consent of the depicted individual (the lead author) in (D). of boreal wildfire conditions: a high-intensity crown fire and a low-intensity surface fire (i.e., 15 sets for each fire). Individual pieces were situated at least 5 cm apart (Figure 1A). K-type thermocouples (≤1 mm diameter) were attached directly to the surface of the samples and connected to data loggers (Lascar, Easylog) buried in the adjacent soil to record the temperature at the sample surface-flame interface during the fire (at 10 s intervals). Given the limited number of loggers available, we chose to attach thermocouples to only one type of charcoal (A: wildfire charcoal) on the basis that any differences in temperature/duration data between the two types of charcoal might be expected to be less than between the charcoal and the contrasting types of wood. All B, C, and D samples had thermocouples attached. Immediately after the fires samples were identified (Figure 1B), thermocouples and loggers removed and samples placed in sealed containers for subsequent laboratory analysis. Wildfire Burn Experiments (1) The high-intensity crown fire was an experimental fire conducted on 15/6/2015 as part of the Canadian Boreal Community FireSmart Project at Ft. Providence, (61◦34’55” N; 117◦10’13” W) and aimed at simulating wildfire conditions. The experimental plot (∼2,000 m2) was a mature stand of black spruce and jack pine, originating from a stand-replacing wildfire in 1931, with a tree density (live and dead, >1.5 m high) of approximately 5,500 stems ha−1and an average canopy tree height of 11 m (a typical mature stand height in this boreal region). The understorey was very sparse, the forest floor composed of litter, mosses and lichens, with some down wood present. Fifteen sets of A-D samples (Figures 1A,B) were placed randomly at least 1 m apart from each other within an area of 15 ×15 m. The fire was fast moving with a total burn time of <5 mins and flame lengths >5 m above canopy height (Figure 1C). The overall fireline intensity was estimated as ∼8,000 kW m−1, which is within the typical range of high-intensity boreal crown fires in Canada (de Groot et al., 2009). The fire resulted in the homogenous and complete consumption of the canopy, thin branches and the understory, with no unburnt or low severity patches remaining within the plot. (2) The low-intensity fire was situated ∼170 km east of the high-intensity fire site (60◦49′38′′ N; 114◦24′28′′ W). It was a slow-moving, surface wildfire (Fire 28, Hay River) affecting a mixed black spruce and jack pine stand of comparable stand characteristics to the high-intensity fire, but with a peaty soil layer of up 50 cm depth. This lightning-caused wildfire burnt across the site on 29/6/2017 during calm atmospheric conditions, advancing at <5 m h−1. This allowed sufficient time to place samples and thermocouples with loggers ahead of one of the flanks of the fire. Fifteen sets of A-D samples were placed 2 m apart along a transect parallel to the fire front. This surface fire led to the patchy smoldering combustion of the forest floor to a depth of up to 30 cm (Figure 1D), charring of trees to ∼1 m height and the Frontiers in Earth Science | www.frontiersin.org 3October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon “torching” of a few individual trees, with some areas remaining unburned. The fireline intensity was estimated to not exceed 500 kW m−1based on flame-length/intensity relationships for jack pine forest (Cruz and Alexander, 2010). Due to the patchy nature of this burn, not all samples were exposed to fire and, but in some cases, deep burning of the organic soil layer also led to the destruction of loggers and loss of data. Further details about the fires are given in Table 1. Laboratory Analyses The air-dry weight of all samples (N=120) was determined before and after exposure to fire, and mass loss calculated as: Mass loss (%)=100∗(1 −postfire weight prefire weight ) Of the total 120 samples, 90 had thermocouples attached, and of those, 19 loggers were lost or did not record properly during the fires. A subset of 40 samples (5 of each type A-D for each fire) spanning a representative range of mass losses and temperature/duration records, and pre-fire reference samples (n =1–5 for each type A–D) were selected for further analyses to allow the potential relationships between mass loss, sample characteristics, and temperature exposure during the fires to be examined. Differential Scanning Calorimetry (DSC) These analyses were carried out using a Mettler Toledo instrument to determine the thermolability (Merino et al., 2014) of the different samples before and after the exposure to the fires. Approximately 4 mg of ground sample was placed in aluminum pans under a flow of dry air (under O2flux; flow rate, 50 mL s−1) at a scanning rate of 10◦C min−1. The temperature ranged between 50 and 600◦C. Samples of indium (mp: 156.6◦C) were used to calibrate the calorimeter. Samples were analyzed in duplicate. The heat of combustion (Q) was determined by integrating the DSC thermographs with respect to time over the exothermic region (150<T<600◦C). The region T<150◦C was not considered as it is dominated by endothermic reactions associated with water loss. The areas under the thermographs were divided into three temperature regions representing different levels of resistance of the organic material to thermal oxidation (Merino et al., 2014, 2015): between 150 and 375◦C (Q1), between 375 and 475◦C (Q2) and between 475 and 600◦C (Q3). The temperature at which 50% of the total energy is released under the given conditions (T50) was also determined (Rovira et al., 2008). Elemental Composition Total carbon (C), nitrogen (N) and hydrogen (H) contents (%) were determined using a LECO elemental analyzer. Values reported were corrected by mineral ash content (as determined by DSC): e.g., C ash-corrected =(C∗100)/Total OM loss. The percentage of the C remaining in samples after fire (%Crem) compared to their initial C content was also calculated as follows using the ash-corrected C% values: %Crem =100 ∗postfire C concentration g/g∗postfire mass (g) prefire C concentration g/g∗prefire mass(g) For the cedar and pine wood pieces (C & D sample types) that were fully charred following fire, this figure is equivalent to the conversion rate of fire-affected C to PyC (%PyC). This rate quantifies the fraction of fire-affected C that is retained as PyC in charred material rather than emitted to the atmosphere (Santín et al., 2015). TABLE 1 | Experimental conditions and resulting PyC losses in previous- and the current studies (n.d., no data; T, ambient air temperature, and RH, relative humidity immediately prior to ignition). Study Environment PyC source and size of pieces Surface fuel load (kg m−2) PyC size range and location Atmospheric conditions (T, RH, wind speed) & fire intensity Mean PYC losses (% mass) Saiz et al., 2014 Savanna Acacia aneura, ∼700◦C slash burn charcoal 0.8–1.0 ∼0.5–1.5 cm soil surface 26–36◦C; 17–45%; ∼5 m s−1;n.d. <8 Santín et al., 2013 Boreal forest Pinus banksiana, slash-pile charcoal 9.4 2–4 cm inside litter 28◦C; 22%; 3 m s−1, ∼8,000 kW m−1 23 (median <15) Tinkham et al., 2016 Laboratory Pseudotsuga menziesii, Pinus monticola, P. contorta wood chips, experimental lab burn* 2.1 <6 mm−7 cm inside litter 17◦C; 36%; <1 m s−1; n.d. 36* This study (high-intensity fire) Boreal forest Pinus banksiana, slash-pile charcoal) >7 1–3 cm litter surface 28◦C; ∼32%; ∼2.5 m s−1, 6,000–8,000 kW m−1 64 Pinus banksiana, wildfire charcoal 84 This study (low-intensity fire) Boreal forest Pinus banksiana, slash-pile >10 1–3 cm litter surface ∼28◦C; ∼25%; <1 m s−1, <500 kW m−1 17 Pinus banksiana, wildfire charcoal 50 *PyC was quantified only in the fraction >6 mm using a method which only considers part of the PyC spectrum. Frontiers in Earth Science | www.frontiersin.org 4October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon Charcoal Reflectance Charcoals have been observed to have varying abilities to reflect light when studied under oil using a reflectance microscope, which has been shown to evolve throughout pyrolysis (e.g., Belcher and Hudspith, 2016). Therefore, we anticipated that charcoal reflectance might vary between the types of fires that created or modified the PyC examined in this study. In order to assess this, the PyC samples were sliced into sections using a band saw and embedded in polyester resin, then ground and polished using a Buehler MetaServ 250 grinder–polisher (Buehler, Neckar, Germany). The samples were placed in the resin such that the top surface exposed to the fire was the polished face following the method described in Belcher and Hudspith (2016). The polished samples were analyzed under oil (RI 1.514) at 238◦C, using a Zeiss Axio-Scope A1 optical microscope, with a TIDASMSP 200 microspectrometer (SMCS Ltd, Baldock, UK). Thirty manual photometric reflectance measurements were taken at cell-wall junctions across the polished surface of all 35 sufficiently sized post-fire samples and the mean charcoal reflectance (Ro%) calculated for each sample. Statistical Analysis Statistical analyses were performed with the software IBM SPSS Statistics 20. The level of significance used for all tests was 5% (i.e., α=0.05). Pearson product-moment correlation coefficients (r) were calculated to test linear correlations between pairs of variables. Variables examined were: mass loss (%), TG-DSC parameters (T50, Q1, Q3), C, N, H contents (%), initial C remaining (%Ci), and Ro (%), as well as maximum temperature recorded on samples (Tmax), time >300◦C (s), integral of duration >300◦C (◦Cs). The threshold of 300◦C was chosen because is widely considered the temperature above which charring begins in woody fuels (Drysdale, 2011), the main chemical transformations of biomass-derived materials during heating start around 200–300◦C (Keiluweit et al., 2010), and it having been used in previous studies focusing on PyC (e.g., Santín et al., 2016b, 2017). To analyze the differences in mass loss (%) between sample types, one-factor ANOVA was performed using sample type as the independent factor. Before the ANOVAs, the equality of variances for the studied groups was tested by the Levene’s homoscedasticity test and the post-hoc Duncan’s multiple range test (groups had equal variances) were performed to identify classes with significantly different means. RESULTS Sample Properties Before the Fires The two unburned wood samples used in this study (types C and D) showed the lowest % C, the highest contents of thermolabile biomass and the lowest thermal recalcitrance (Q1 >33%; T50 <420◦C), with the cedar wood (D) being more thermolabile than the pine wood (Table 2,Figure 2). In comparison, both charcoal types (A and B) showed much higher C content and thermal recalcitrance, and lower thermolability (Q1 <21%; T50 >447◦C). As expected, the slash-pine charcoal (B) presented, by far, the highest C% and thermal recalcitrance of all the samples before fire (Table 2). TABLE 2 | Arithmetic mean values and standard deviation (italicized) of elemental composition and outcomes of thermal analysis (Differential Scanning Calorimetry; DSC) before and after fire, and of mass loss and reflectance (Ro) after fire, of the four different sample types of samples for the high-intensity (H) and low-intensity (L) fires. Pre-fire Post-fire C (%) N (%) H (%) T50 (◦C) Q1 (%) Q3 (%) C (%) N (%) H (%) T50 (◦C) Q1 (%) Q3 (%) Ro (%) Mass loss (%) Wildfire charcoal 72.5 0.23 2.9 447 20.5 43.8 H 80.2 0.26 2.31 486 9.5 50.6 1.73 77.9 0.9 0.03 0.2 4.1 0.1 3.6 4.3 0.03 1.04 14.5 1.0 12.3 0.74 14.1 L 78.8 0.43 2.43 468 14.1 43.9 0.92 53.7 6.4 0.30 0.51 19.0 6.6 14.0 0.36 31.0 Slash-pile charcoal 80.2 0.16 2.6 510 6.4 71.4 H 86.8 0.11 1.29 523 6.9 76.4 2.80 52.9 n.a n.a n.a n.a n.a n.a 2.0 0.17 0.24 22.5 1.8 9.1 0.54 19.4 L 85.8 0.17 1.36 519 5.9 76.3 2.37 24.8 2.2 0.14 0.24 14.7 1.4 5.8 0.98 22.7 Pine wood 46.1 0.07 5.4 417 33.8 25.8 H 76.3 0.29 2.74 470 13.8 46.9 1.47 89.0 n.a n.a n.a n.a n.a n.a 2.2 0.03 0.36 7.1 2.5 4.8 0.26 6.6 L 73.2 0.38 2.83 459 18.4 39.6 0.97 65.9 2.3 0.06 0.41 10.4 4.2 6.9 0.27 9.0 Cedar wood 50.9 0.23 n.d. 395 40.5 4.6 H 75.6 0.32 2.65 458 19.1 35.2 1.46 96.2 0.2 0.04 1.3 0.5 1.6 3.6 0.03 0.20 3.1 3.6 3.7 0.34 3.8 L 72.9 0.34 2.75 453 21.1 32.5 1.33 70.4 1.2 0.09 0.28 10.5 5.2 7.1 0.16 19.3 T50 index: loss of 50% of energy during DSC; Q1 and Q3: percentage contributions of combustion heats obtained by DSC analysis. Note that mean mass loss values given here represent a subset of samples analyzed for specific characteristics given in this table. Mean mass losses for all samples are given in Table 4. Frontiers in Earth Science | www.frontiersin.org 5October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon FIGURE 2 | Comparison of the differential scanning calorimetry curves (arithmetic mean values of 5 samples each) of the four samples types before and after the fire. TABLE 3 | Pearson product-moment correlation coefficients (r) between selected chemical and reflectance characteristics of the studied samples, and key fire parameters. MaxT IntSecs >300◦C Mass loss% Ro% C% H% T50 MaxT 1 0.228 0.644** 0.462* 0.351 −0.486** 0.374 IntSecs>300◦C 1 0.483** 0.548** 0.168 −0.322 0.216 Mass loss% 1 −0.148 −0.418** 0.190 −0.379* Ro% 1 0.768** −0.785** 0.722** C% 1 −0.900** 0.908** H% 1 −0.856** T50 1 Only the most relevant correlations are shown here. For correlations among other variables analyzed see Table S1. Number of samples: 28–47 (see Table S1) *Significant at the 0.05 level. **Significant at the 0.01 level. Mass Losses and Their Correlations With Temperature/Duration Records Considering data from both fires together (N=101), wildfire charcoal (type A), exhibited a significantly higher mean mass loss (66.5 ±25.2%) compared to slash-pile charcoal (Type B, 41.7 ±27.2%). Pine wood (type C) showed higher mean mass losses (78.2 ±14.9%) compared to both charcoal types A&B, and cedar wood (type D) the greatest losses (83.8 ±18.9%). However, given the high variability in the data, the differences between mean mass losses between the groups of wildfire charcoal, pine wood and cedar wood were not significant. For the combination of all the samples with temperature data associated (N=71) and of the three fire parameters analyzed, Tmax showed the strongest positive correlation (r=0.64) with mass loss (Table 3,Figure 3). Correlation between mass loss and the integral >300◦C was also significant, but weaker (r=0.48). No significant correlation was found between mass loss and time >300◦C (r=0.33) (Supplementary Table 1). When comparing the two fires (Figure 4), the mean maximum temperature (Tmax) recorded on all samples in the high-intensity fire was 855◦C (range 661–1005◦C; N=39), whereas in the low-intensity fire it was only 434◦C (range 40– 792◦C; N=32). Indeed nearly all Tmax values recorded on samples in the high-intensity fire exceed those of low-intensity fire (Figure 4). The high-intensity fire also exhibited longer mean residence times (>300◦C) than the low-intensity fire, Frontiers in Earth Science | www.frontiersin.org 6October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon FIGURE 3 | Mass loss versus maximum temperature recorded for all cedar wood, pine wood and wildfire charcoal samples (N=71; slash-pile charcoal samples did not have temperature data associated). Green circles, high-intensity fire; blue circles, low-intensity fire. although these differences between fires were less pronounced (Figure 4). Mean sample exposure durations >300◦C were 245 s (range 70–570 s) and 142 s (0–620 s) for the high- and low-intensity fire respectively, and the mean integrals >300◦C were 53,656◦C s (range 20,689–122,798◦C s) and 18,589◦C s (range 0–99,327◦C s) for these respective fires. All samples exposed to the high-intensity fire, and most samples exposed to the low-intensity fire experienced mass losses, with these mean losses being significantly higher in the wood pieces (80.7 ±16.8%; n=49) than in charcoal pieces (54.7 ± 28.5%; n=52) and in the high-intensity than in the lowintensity fire (80.4 ±18.6%; n=47 vs. 53.7 ±27.4%; n =54). These significant differences in losses between fire types agree with their contrasting temperature records (Table 4, Figure 4). Regarding the differences in mass losses between the two fires for the different sample types, mass losses were consistently smaller in the low-intensity than in the highintensity fire and in both cases following the order D>C>A>B (Figure 5). The variability of these losses was higher in the lowintensity fire in all cases except type D, the slash-pile charcoal (Figure 5). Thermal-, Elemental-, and Reflectance Characteristics Results of the elemental-, DSC- and reflectance analyses for each fire are given in Table 2, and DSC curves for each sample type in Figure 2. As highlighted before, the DSC analyses show that, prior to the fires, the slash-pile charcoal was the most thermally recalcitrant (i.e., least thermolabile) material used in this experiment (mean T50 =510◦C) followed by wildfire charcoal (447◦C), pine wood (417◦C), and cedar wood (395◦C). This order is reflected also in the increasing values for Q1 and decreasing values for Q3 (Table 2). After sample exposure to fire, the same order is retained, but with substantial increases in thermal recalcitrance for the wood samples following charring, and only modest increases for the charcoal samples (Table 2,Figure 2). These increases in recalcitrance were greater in the highthan the low-intensity fire, which matches the greater thermal exposure of the samples in the former (Table 4). The C concentration also increased in all sample types following fire exposure, but this change was much more pronounced for the wood samples than charcoal samples (Table 2). Given that exposure to fire does not only lead to mass- and therefore C losses through combustion, but also to relative enrichment in C in the remaining PyC through pyrolysis, we also calculated the % of PyC remaining in relation to the total C content of the samples prior to exposure (Table 4). Mean values of PyC% remaining were, respectively for high- and low-intensity fires: 26.2 vs. 48.3 (wildfire charcoal), 44.8 vs. 80.8 (slash charcoal), 15.7 vs. 54.2 (pine wood) and 8.2 vs. 30.8 (cedar wood). As highlighted in the methodology, for the fully charred wood samples C&D, this parameter represents the fuel C to PyC conversion rate used in Santín et al. (2015). It is interesting to note that the low-intensity fire resulted in a higher mean PyC production ratio (54.2 and 30.8%) than the highintensity fire (15.7 and 8.2%) for the pine and cedar wood samples respectively. Regarding N and H concentrations, N increased and H decreased after fire for most sample types (Table 2). Charcoal reflectance (Ro) following the fires was highest for slash-pile charcoal with values for the other sample types being around half this value (Table 2). Values following the high-intensity fire overall exceed those following the low-intensity fire with respectively mean Ro values of 2.80 vs. 2.37 (slash charcoal), 1.73 vs. 0.92 (wildfire charcoal), 1.47 vs. 0.97 (pine wood), and 1.46 vs. 1.33 (cedar wood). Correlations Between Sample Characteristics, Charcoal Reflectance and Fire Properties When considering all samples exposed to the fires together, significant correlations were found between some sample chemical and reflectance characteristics, and with fire properties (Table 3). For example, T50 showed positive correlations with C% (r=0.91), Ro (r=0.72) and a negative correlation with mass loss (−0.38) and with H% (r= −0.86). C% correlated with Ro (r=0.72) and mass loss (−0.42). Maximum temperature at the fuel-flame interface (Tmax) showed not only a positive correlation with mass loss (r=0.64) as presented in more detail above, but also with Ro (r=0.46) and H% (r= −0.49). The duration integral of exposure >300◦C was similarly well correlated to Ro (r=0.55), and weaker (but still significantly) correlated to mass loss than Tmax (r=0.48). All other relationships between fire characteristics, including exposure duration >300◦C, showed weaker and non-significant relationships with sample chemical and reflectance characteristics (Table S1). Frontiers in Earth Science | www.frontiersin.org 7October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon FIGURE 4 | Maximum temperatures (◦C) (left) and residence times of T >300◦C (right) recorded at the sample-flame interface during the low- and high-intensity fires for all sample types combined (N=71). Circles indicate outliers. TABLE 4 | Mean values (and range) of temperature characteristics, mass losses and % PyC remaining, separated by sample and fire types. Tmax (◦C) Time >300◦C (s) Integral >300◦C (◦Cs) Mass loss (%) PyC remaining (%) HIGH-INTENSITY FIRE Wildfire charcoal (n=12; 6) 882 283 61197 82 26.2 (735–1005) (70–410) (20689–101113) (54–100) (6.0–53.6) Slash-pile charcoal (n=15; 5) n.d. n.d. n.d. 63 44.8 (33–100) (19.4–73.7) Pine wood (n=13; 5) 833 249 53,132 86 15.7 (661–1003) (80–570) (21406–122798) (57–100) (5.1–31.7) Cedar wood (n=14; 3) 851 208 47681 98 8.2 (747–1002) (90–460) (22932–97238) (90–100) (4.7–14.2) LOW-INTENSITY FIRE Wildfire charcoal (n=10; 5) 390 16 13685 56 48.3 (50–603) (0–590) (0–43290) (2–82) (17.2–97.1) Slash-pile charcoal (n=15; 5) n.d. n.d. n.d. 17 80.8 (12–64) (39.5–100) Pine wood (n=11; 5) 432 90 14286 41 54.2 (85–715) (0–300) (0–51011) (5–80) (31.5–65.8) Cedar wood (n=11; 4) 475 172 27349 55 30.8 (40–792) (0–620) (0–99327) (5–99) (16.3–44.1) % PyC remaining is the fraction of C remaining compared to the total C content of the samples prior to fire exposure. (The post-fire wood samples selected for chemical analysis were fully charred, therefore all the C in them is considered PyC). Sample size (n) applies to all parameters except for PyC remaining where the second italicized figure applies. n.d., no data. DISCUSSION Mass Losses Compared to Previous Studies The only two previous studies examining PyC recombustion under field conditions reported mean mass losses of <8% for charcoal pieces of ∼0.5–1.5 cm (Saiz et al., 2014) and of 23% (median value <15%) for charcoal pieces of 2–4 cm (Santín et al., 2013). In the current study, mass losses were substantially higher except for the slash-pile charcoal exposed to the low intensity fire, which had a mean loss of 17.0% (median 11.8%) (Table 4). In the high-intensity fire, slash-pile charcoal showed mean losses of 63.6% (median 58.6). Wildfire charcoal exhibited losses of 84.0% (median 81.9%) and 49.7% (median 58.6%) in the high- and low-intensity fires, respectively (Table 4). Thus, all losses recorded here across the two types of PyC and fires (i) exceed those previously reported from in-situ field studies and, except for the slash-pile charcoal in the low intensity fire, (ii) exceed also the mean losses of 36% for PyC fragments >6 mm reported in the laboratory study by Tinkham et al. (2016). The data presented in Frontiers in Earth Science | www.frontiersin.org 8October 2018 | Volume 6 | Article 127 Doerr et al. Fire as a Removal Mechanism of Pyrogenic Carbon FIGURE 5 | Mass losses for the different type of samples for the (A) high - and (B) low-intensity fire. Circles and stars indicate outliers. Total N=101. Table 1 allows comparison of the current study with fire and PyC characteristics in these previous studies. The comparatively low losses in the study by Saiz et al. (2014) (<8%; Table 1) could be a reflection of a relatively high thermal resistance of their PyC combined with a low capacity of the fire to lead to PyC combustion. The source material, Acacia aneura, has roughly twice the density (1 g cm−3) of pine wood and was produced in a slash burn, which typically have much longer flame residence times than wildfire and therefore lead to a greater thermal degradation resistance of the charcoal produced than in a wildfire (Massman et al., 2008). The fire, although considered intense in the context of savanna fires, had the lowest fuel load and highest wind speed of all experimental fires listed in Table 1 and hence is likely to have led to a burn of comparatively low intensity and flame residence time. The experimental crown fire examined by Santín et al. (2013) was carried out in the same fuel complex and study region as the current study and used the same source of slash-pile charcoal. Fuel load, atmospheric conditions and fireline intensity (∼8,000 kW m−2) were also broadly comparable to the high intensity fire of current study (Table 1). The main difference between these two studies is that, in the current study, PyC pieces were placed on the forest floor surface, whereas in Santín et al. (2013) they were placed at ∼2 cm depth in the forest floor layer. In the 2013 study, this could have resulted in the initial protection from thermal exposure of the charcoal pieces until the litter layer itself was consumed to this depth. This may have contributed to the substantially lower combustion of the slash-charcoal of 23% reported by Santín et al. (2013) compared to the 64% in the high-intensity fire of the current study (Table 1). The mass losses of both slash-pile charcoal in the highintensity fire and of wildfire charcoal in both fires of the current study are substantially greater than the 36% found in the laboratory burns of Tinkham et al. (2016) (Table 1). Their experiment, however, is not directly comparable to the current study not only due to it being carried out under laboratory conditions using masticated fuels. They also considered only PyC pieces >6 mm and within those, only the fraction that is resistant to chemical degradation based on the CTO375 method (Hatten and Zabowski, 2009). In addition, the sample placement on the litter surface, the substantially higher fuel loads and warmer and drier atmospheric conditions in the current study are all more conducive to a greater combustion completeness than those in the laboratory burn by Tinkham et al. (2016). Mass Losses and Their Correlations With Temperature/Duration Records Our research design allowed us to link PyC losses directly to temperature recorded on the fuel-flame interface of the respective PyC pieces. That mass losses show a reasonably strong correlation with Tmax (r=0.65; Table 3) is perhaps not surprising. The thermocouple tip (≤1 mm diameter) reflects the temperature at the fuel-flame interface, which is influenced by the energy received from the surrounding combustion combined with that released from the burning PyC piece itself. Tmax could therefore be expected to be at least somewhat related to the mass loss of a PyC piece, with higher temperatures reflecting greater losses. What may be surprising is that mass losses are only weakly correlated with exposure duration >300◦C (r =0.33, not significant; Table S1). Three Hundred degree is widely considered the threshold above which charring begins in woody fuels (Drysdale, 2011) and hence the duration above this threshold could be expected to correlate more closely with mass loss than Tmax. Given that the correlation between mass loss and the integral >300◦C was stronger (r=0.48) than for time >300◦C (and significant), suggests that the faster rate of combustion associated with the process that generates higher temperature readings in thermocouples is more important in the overall mass loss process observed here than the total time during which mass loss occurs. This notion is supported by a study examining organic matter in a soil block affected by experimental laboratory burns, where the time >300◦C also showed no correlation with C loss (Merino et al., 2018). Frontiers in Earth Science | www.frontiersin.org 9October 2018 | Volume 6 | Article 127