Abstract
Wildland fire can cause significant damage but is also a natural process that is key to the healthy functioning of many ecosystems worldwide. Primary fuels for a wildland fire are the dead foliage and small branches which accumulate as litter on the ground. A cone calorimeter was used to measure the various aspects of these fuels. A single sample of preignition gases from the live leaves of seven plant species were vacuum collected on quality filters and within super-chilled solvent mixtures. GC-TOFMS (1D) and GCxGC-TOFMS (2D) were used to characterize the “white” smoke emissions. The vegetation chemicals were separated into 4 categories: hydrocarbons (CH), oxygenated organics (CHO), unknown peaks (UNK), and organic non-metals (ONM). The multivariate paired Hotelling T2 test determined that the composition of the white smoke as described by the relative number of peaks in the four chemical groups differed significantly between 1D and 2D (Prob > F3,4 = 0.00004). In contrast, the relative peak area percentages in the four chemical groups did not differ between 1D and 2D (Prob > F3,4 = 0.1258). The Molecular Chemical Maps (MCMs) were used to identify chemical trends between the known and unknown chemicals in live oak and longleaf pine. Application of the 2D technique may provide more detailed information necessary to improve the numerical modeling of wildland fire behavior and emissions production.
Download PDF
Full Article
Utilizing Two-Dimensional Gas Chromatography Time of Flight Mass Spectrometry (GCxGC ToFMS) to Characterize Volatile Products from Pyrolysis of Living Vegetation Foliage
Roderquita K. Moore,a,* Mark A. Dietenberger,a Doreen H. Mann,a Patricia K. Lebow,a and David R. Weise b
Wildland fire can cause significant damage but is also a natural process that is key to the healthy functioning of many ecosystems worldwide. Primary fuels for a wildland fire are the dead foliage and small branches which accumulate as litter on the ground. A cone calorimeter was used to measure the various aspects of these fuels. A single sample of preignition gases from the live leaves of seven plant species were vacuum collected on quality filters and within super-chilled solvent mixtures. GC-TOFMS (1D) and GCxGC-TOFMS (2D) were used to characterize the “white” smoke emissions. The vegetation chemicals were separated into 4 categories: hydrocarbons (CH), oxygenated organics (CHO), unknown peaks (UNK), and organic non-metals (ONM). The multivariate paired Hotelling T2 test determined that the composition of the white smoke as described by the relative number of peaks in the four chemical groups differed significantly between 1D and 2D (Prob > F3,4 = 0.00004). In contrast, the relative peak area percentages in the four chemical groups did not differ between 1D and 2D (Prob > F3,4 = 0.1258). The Molecular Chemical Maps (MCMs) were used to identify chemical trends between the known and unknown chemicals in live oak and longleaf pine. Application of the 2D technique may provide more detailed information necessary to improve the numerical modeling of wildland fire behavior and emissions production.
DOI: 10.15376/biores.17.1.862-889
Keywords: GCxGC TOFMS; Molecular chemical mapping; Live fuel properties; Wildfire smoke
Contact information: a: USDA Forest Service, Forest Product Laboratory, One Gifford Pinchot Drive, Madison WI 53726; b: USDA Forest Service, PSW Research Station, Fire and Fuels Program, Riverside, CA 92507; *Corresponding author: roderquita.k.moore@usda.gov
INTRODUCTION
In the United States for the last 10 years there have been an average of 7.5 million acres burned by 71,950 fires, causing loss of life and significant damage to structures (Butry et al. 2001; Kukhar et al. 2020). The smoke generated by these fires can usually be seen from outer space, is a major human health concern, and has significant impact on the atmosphere. Similarly, in many locations, fires have released energy at rates higher than the current vegetation is adapted to, potentially causing long-term effects including soil loss, conversion to other vegetation types, loss of plant and animal habitat, and decreased water quality. Prescribed burning can be an effective tool to reduce hazardous fuel loads, which decreases the risk of larger fires while providing ecological benefits; however, it produces smoke, which is regulated in the United States by the Clean Air Act and its amendments (Martin and Dell 1978; Green 1981; Waldrop and Goodrick 2012; National Wildfire Coordinating Group 2018).
The burning of vegetation is a complex phenomenon (Bytnerowicz et al. 2009) composed of two main processes – pyrolysis (thermal decomposition) and combustion (oxidation), which are still not well-described even after a century of study. Pyrolysis of the solid material (vegetation) produces gases and a solid material called char. The gases can react (combust) with oxygen, producing a flame. The char can also combust. Lignin, cellulose, and hemicellulose are the major components of vegetation. In live plants, pyrolysis (Shafizadeh 1968; Yokelson et al. 1996; Ward 2001; Dietenberger et al. 2020) includes several stages: (1) evaporation of water, (2) distillation of terpenes and light volatile organic carbons (VOCs), (3) pyrolysis of aboveground plant parts (foliage, stem, bark) (Matt et al. 2020) with associated chemical reactions that produce volatiles such as VOCs, trace gases, tar, and particles, and (4) charring of these same plant parts. Pyrolysis of belowground parts (roots) may also occur under certain ground fire settings (e.g. Rein 2013). The pyrolysis process of these components is complex and temperature dependent (Neves et al. 2011). Furthermore, there are two main zones of combustion, the first within the flame sheets of the diffusion flame that partially or entirely capture the volatiles and react them with entrained air, and the second within the organic char layer where the penetrating air oxidizes the char interior surfaces. Therefore, emissions contain a large component of uncombusted gases and aerosols (grouped as volatiles) as sourcing from pyrolysis that circumvent the combustion zones or partake in inefficient combustion processes (Yokelson et al. 1997; Leonellia et al. 2017; Barboni et al. 2020).
The VOC formation mechanism is temperature dependent (Sekimoto et al. 2018; Yokelson et al. 1996, 1997; Zhao et al. 2017). The depolymerization of lignin occurs in a low temperature range of 300 to 500 °C. The degraded lignin monomers consist of substituted methoxy compounds such as guaiacols, isoeugenol, and syringols. The cellulose and hemicellulose temperature ranges for depolymerization are 300 to 400 °C. There are several reactions occurring simultaneously that emit VOCs, for example hydroxyacetone, acetaldehyde, and acetic acid, during fragmentation from the depolymerization of cellulose and/or hemicellulose. During pyrolysis of biopolymers, VOC emissions are equally distributed during low and high temperatures. The concentration of chemicals may vary as temperature increases. The release of oxygenated compounds increases carbon content in the solid fuel, which contributes to the strong bonds present in char. Aromatization of nonsubstituted and polycyclic aromatic hydrocarbons (PAHs) increases during high temperatures. The functional groups represented in the high-temperature profile are alkenes, alkynes, dienes, PAHs, nitriles, amides, nitrites, nitrates, and thiols/sulfides. The low-temperature functional groups are alcohols, ketones, acids, methoxy groups, esters, and amines. Description and measurement of the chemistry of pyrolysis and combustion is challenging because it is dynamic and requires multiple instruments to measure the many compounds. While numerous instruments have been used effectively, spatial, temporal and chemical resolution of measurements continues to be a challenge. Application of instruments capable of high-resolution measurements under tightly controlled combustion experiments applied to the chaotic conditions associated with the diffusion flames of a wildland fire is also a challenge, but this may yield information not previously known. Over the past decade, several instrument intercomparisons have been performed to identify and compare various components of smoke (Christian et al. 2004; von Bobrutzki et al. 2010; Warneke et al. 2011; Li et al. 2019). More characterization and quantitative comparison are needed to better understand VOCs and primary organic aerosols.
Smoke emissions from fires can be grouped into (A) air-borne particles, known as brown carbon (organic carbon or aerosols) and black carbon (soot) (Olson et al. 2015; Leonellia et al. 2017), (B) high molecular weight organic compounds, known as heavy tar (condensable on filters) and light tar (condensable in cooled solvent) (Leonellia et al. 2017; Knorr et al. 2019), (C) fuel gases (particularly methane and other VOCs) (Yokelson et al. 1996), and (D) combustion gaseous products such as CO2, CO, and H2O. Indeed emissions A-D have suggested common measurement methods, such as non-dispersive infrared (NDIR) used for combustion gases, Fourier transform infrared (FTIR) used for fuel gases, gas chromatography mass spectrometry (GC/MS) used for light tar (or bio-oil), and laser light extinctions or impactors used for brown and black carbon. However, a significant amount of emission mass is still unaccounted for by these particular measurement methods, particularly for smoldering, which may be due to undetected or unidentified aerosols and tar as sourced from the primary pyrolysis volatile products of condensed matter, such as that of forest litter, live foliage, and wood.
There are several existing smoke emission repositories containing information from the field and laboratory of wildland fires. While emissions measurement has occurred for the past 50 years, the repository lacks diversity among vegetative fuel sources; more data exists to support flaming combustion, i.e., wildfire and, prescribed burning (Prichard et al. 2020). This lack of diversity may hinder prediction for the emissions of pollutants. More studies are needed for all pollutants in the area of wildfires, residual smoldering and organic soils, and specific vegetation types, including grasslands, northern forests and broadleaf deciduous forests. Certain trace gases and aerosols such as non-methane organic compounds (NMOCs) were underrepresented. The limited detection of trace aerosols or NMOCs are because of detection sensitivity and the resolving power of the instrumentation (Moore et al. 2015a,b, 2017a,b,c).
Yokelson et al. (2013) investigated the quantitative measurements on NMOCs emissions. During this investigation, it was discovered that approximately 72% of NMOCs emitted by some fuel types were unidentified using the current technology. The emission gap of NMOC limits the possibility of discovering the impact of biomass burning on the atmosphere, modeling the impact of the local to global chemistry of both natural and anthropogenic fires, and developing preventative solutions to deactivate chemicals that contributes to the formation of ozone (O3) or organic aerosols. Most of the unidentified species have molecular mass greater than 90 m/z and thus may be semi-volatile or of intermediate volatility. State of the art instrumentation was used such as: open-path Fourier transform infrared (OPFTIR) spectrometer, proton-transfer-reaction mass spectrometry (PTR-MS), proton-transfer ion-trap mass spectrometry (PIT-MS), negative-ion proton-transfer chemical-ionization mass spectrometry (NI-PT-CIMS), and gas chromatography with MS detection (GC-MS). However, efforts to better characterize the chemistry and oxidation products of the species with high resolution mass spectrometry is still needed for complete mass balance characterization of organic trace gases and investigation of the effects of oxidation, temperature, and humidity on full mass scans.
Several analytical instruments have been used to characterize vegetation smoke emissions. GC/MS is often used; however, GCxGC is expected to provide an enhanced separation (typically one to two orders of magnitude) compared with one-dimensional GC. GCxGC uses a two-tandem column set that has two different stationary phases to separate complex mixtures. The chemical phases (column packing) of these columns are nonpolar for the 1st dimension and semipolar for the 2nd dimension. The samples are injected into the first column and travel to a thermal modulator with cryogenic jets and are then injected into the second-dimension column (Liu and Phillips 1991). The separation power of this technique has the potential to classify known and unknown compounds through chromatographic retention patterns.
When coupled with time of flight mass spectrometry (TOFMS), this technique will allow constituents of unresolved complex mixtures to take advantage of continuous full range nonskewed mass spectral information and fast acquisition rates. GCxGC has been successfully applied to a number of complex mixtures such as petroleum (Mao et al. 2009; Kohl et al. 2010), tobacco (Lu et al. 2003), flavors (Bordiga et al. 2009), fragrances (Shellie et al. 2004), essential oils (Dimandja et al. 2000), and metabolic (Welthagen et al. 2005) profiling. These applications have shown that GCxGC is able to separate mixtures that contain up to 10,000 compounds (Dimandja 2004). Panic and Gorecki (2006) characterized trace toxic compounds in the environment. The characterization and identification of complex mixtures increase the discovery of new chemicals for medicinal purposes and may provide solutions to grand challenges that are facing our environment.
Although GCxGC has been applied to the pyrolysis of a dried, ground sample of a grass (Miscanthus giganteum) (Groenewold et al. 2017) and to biomass burning smoke in a limited fashion (Hatch et al. 2015; Jen et al. 2019), to our knowledge the GCxGC TOFMS has not been applied to intact pyrolyzing of live wildland fuels. The sample collection, sample preparation and instrument parameters play a significant role in the characterization of the smoke samples. The present study used a two-phase continuous collection process to assist with lowering temperature of the smoke volatiles in a super chilled solution. The super-chilled solution was used to stop degradation or any high volatiles from escaping. In literature, smoke samples went through desorption, and trapping during sample preparation which requires additional heating before injection.
The aim of this study was to utilize GCxGC TOFMS methodology to characterize primary volatiles generated directly from heating via radiation of a thin layer of foliage in a cone calorimeter (pyrolysis gases). Adapting GCxGC-TOFMS offered an opportunity to characterize pyrolytic emissions from wildland fuels in greater detail. This information, coupled with measurement from related studies (Safdari et al. 2018, 2019, 2020; Amini et al. 2019; Scharko et al. 2019a,b; Phillips et al. 2020) can provide (1) a description of primary pyrolysis volatiles helpful for numerical pyrolysis modeling in the presence of air; and (2) a description of more organic compounds that are unfettered by destruction or changes in the close proximity to hot flames. This paper will identify light tar volatiles and characterize these volatiles through molecular chemical mapping. Analysis of heavy tars and aerosols collected on the filters will be presented elsewhere.
EXPERIMENTAL
Materials
Nursery-grown seedlings of seven plant species found in wildland fuel beds in the southeastern United States were shipped to the USDA Forest Service Forest Products Laboratory where they were maintained in a greenhouse until being tested (Table 1). Dichloromethane (80%) acetone (20%) solution (labeled as DCM), and tetracosane 99% were purchased from Sigma-Aldrich (St. Louis, MO, USA), isopropanol (IPA) was purchased from Thermo Fisher Scientific (Waltham, Massachusetts, USA), and dried ice was purchased locally to keep the liquids super-chilled during DCM collection of the white smoke.
Table 1. Leaf Type Names, and Moisture Content (MC) of Nursery Plants for the Southern United States
Plant Sample and Preparation
The same plants were used in this project as in Matt et al. (2020). High heating values as well as gross heat of combustion calculated from elemental composition of the plants were presented in Matt et al. (2020). The living leaf from a plant contains major components such as water, lipids, glucose, fructose, protein, pectin, hemicellulose, starch, cellulose, phenolic compounds, lignin, silicates, and minerals (Raven et al. 1981; Matt et al. 2020). This complexity poses additional challenges to identification of pyrolysis products compared with that of wood. As a live leaf is heated, it begins to emit water vapor and volatiles from the leaf components. Contrary to the assumption made with dead fine fuels that all moisture is volatilized before other volatile products are emitted, live foliage may contain a significant amount of moisture (Table 1) when pyrolysis begins (Fletcher et al. 2007; Prince and Fletcher 2014; Yashwanth et al. 2016; Matt et al, 2020). The leaf typically has structural layers of waxy cuticle, upper epidermis, mesophyll, vein, lower epidermis, and lower cuticle, (Biology Libre texts 2020) that ultimately affect the thermal and mass transports of the heated leaf (Raven et al. 1981). The complexity of a live leaf coupled with the large diversity present in the plant kingdom may require a corresponding complex analysis of plant pyrolytic properties. Mature leaves from a plant were snipped with scissors and placed in a pre-weighed petri dish. After about 3 grams of foliage had been harvested (some plants had multiple leaves that added up to 3 grams), the petri dish was covered, weighed, and delivered promptly to the cone calorimeter. Mature leaves from a plant were snipped with scissors and placed in a pre-weighed petri dish. After about 3 grams of foliage had been harvested (some plants have multiple leaves that add up to 3 grams), the petri dish was covered, weighed, and delivered promptly to the cone calorimeter. In the particular case of longleaf pine needles, the typical lengths (203 to 457 mm) were trimmed to the diameter of the cone calorimeter sample holder (80 mm), and approximately two or three layers were required to adequately cover the area. Foliage from the other species was processed similarly except for trimming.
Enhanced Cone Calorimetry Method for White Smoke Collection
Several authors have reported the emission of white smoke prior to ignition, which results from pyrolysis and distillation (Yokelson et al. 1996; Leonellia et al. 2017; Pallozzi et al. 2018). A well-developed method for collecting cigarette smoke for analysis with GCxGC-TOFMS (Knorr et al. 2019) was adapted to capture this white smoke. Key features of the collection technique were a vacuum pump to draw the white smoke through an inert GC quality filter to capture the heavy tar and then through two in-series impinging bottles with chilled DCM solution to capture the light tar and condensable combustible gases. Since the cone calorimeter was operated in the unpiloted mode, the leaves were completely consumed by smoldering down to the nil ash mass that started when the pyrolysis commenced as the leaves were drying down to nil ash mass.
An enhanced cone calorimetry test was devised to examine leaf drying and volatile properties closely. The vegetation cone holder and smoke sampler (Fig. 1) was modified by adding six 36-gauge Type K thermocouples (T/C) to the cone calorimeter data acquisition board. Three thermocouples were laid on the bottom steel mesh that is on the top of the sample holder. The leaves were placed singly in a thin bunched layer on the steel mesh. A leaf was placed over each thermocouple, and the circular exposed area was filled in by the leaf. The remaining three thermocouples were placed on top of the leaves and the top steel mesh was placed on top and was pressed/hooked to the bottom mesh to keep all the leaves flat so that they all received the same irradiance (35 kWm-2) from the cone’s heater. The sample holder was placed in such a fashion to minimize the impact of the six T/Cs on the mass loss measurement. The cone calorimeter was in a waiting mode to begin the test at the selected irradiance.
Fig. 1. Enhanced cone calorimeter setup used to capture preignition “white” smoke from live foliage samples
Preliminary testing showed that an irradiance of 35 kWm-2 produced temperatures within the leaf of 400 to 600°C, which have been found to produce maximum yield of tars from the leaf (Amini et al. 2019); these lower temperatures avoid the high-temperature secondary pyrolysis process. The initial radiant heating of the leaves showed a plateau near 100°C on both the exposed and shaded side of the leaves, which indicates both moisture desorption and a high thermal conductivity due to the high moisture content of the live leaves. This differs from other temperature profiles reported for single leaves at higher heating rates (Fletcher et al. 2007). After drying, the exposed surfaces (via T/C measurements) show a fairly rapid rise up to around 500°C, thereby initiating both pyrolysis and char oxidation (i.e., smoldering) on the exposed surface at around 200°C. Meanwhile, the shaded T/Cs had a delayed temperature rise, due to the low thermal conductivity of the dried leaf components. Because of this temperature differential across the leaf thickness, the shaded side could still be evaporating the absorbed water, while simultaneously, the exposed surface could also be emitting water vapor that is primarily from the combustion products of organic char oxidation (e.g., Shotorban et al. 2018).
However, this overlap of water vapor sources appears to be present for a relatively short time interval, such that white smoke was primarily the products of both pyrolysis and organic char oxidation. Furthermore, as the white smoke transitions to the barely visible black smoke at a later degradation time, the leaf residue is undergoing primarily char oxidation, as evidenced by the continuing mass loss and heat release rate and leaf temperatures remaining high. All these processes were measured in the cone calorimeter as mass loss rate and heat release rate as functions of time. During the moisture desorption phase, the cone calorimeter measured leaf mass loss while there were nil measured CO2 and CO production, nor O2 consumption, nor red laser smoke extinction recorded, but the H2O sensor signal was recorded. Therefore, the additional advantage of the 35kWm-2 heat flux exposure was the ability to separate leaf degradation events into the phases of (1) moisture desorption, (2) pyrolysis and char oxidation (white smoke), and (3) char glowing. In the case of longleaf pine live needles, white smoke emitted was thus coming from about 67% of the leaf dry mass in these tests. This was calculated as initial mass minus moisture mass, and then the dried residue mass was separated into pyrolysis (white smoke source) and char masses. Samples were collected during the appearance of white smoke Because numerous leaves were used per test, abundant white smoke was produced, resulting in an averaged quantity for each test. A total of seven measurements (one per species) was collected.
GC and GCxGC Sample Preparation
After each calorimetry test, the tubing and impinger bottles were rinsed with IPA. The soiled filters were sealed in separate glass jars and stored in the refrigerator. The smoke samples were filtered, evaporated using a rotary evaporator, and transferred to preweighed vials, sealed, and stored in the freezer. The injections of the smoke samples were removed from the freezer, thawed, spiked with tetracosane (internal standard), mixed, and injected. The tetracosane was used for calibration and alignment analyses. There were several injections done for development of the GC and GCxGC method. However, only two injection data were used. A rotary evaporator was used to dry the samples. The mass of the concentrated samples ranged from 1.5 to 15.4 mg. Samples were spiked with tetracosane and dissolved with dichloromethane. National Institute of Standards and Technology (NIST) and Wiley libraries were used to identify the chemicals detected in the samples. Only chemical masses greater than 90 were considered.
Instrumentation Operations
The instrument parameters for the analyses of smoke samples are shown in Table 2. The injection split and the modulation time varied for different plant species. Samples were injected into the same system with different parameters. The acquisition rates for the GC TOFMS and GCxGC TOFMS were 10 and 150 spectra/s, respectively. The split ratio varied among samples; LO, SP, and DP required a lower split ratio. The maximum modulation time was 5 seconds for each smoke sample.
Table 2. GC-TOFMS and GCxGC-TOFMS (Pegasus BT 4D) Instrument Parameters
GC and GCxGC Data Preparation
GC/MS is one of the instrument combinations commonly used to characterize VOCs (ASTM 2018; ISO 2011). GC-TOFMS (1D) and GCxGC-TOFMS (2D) were used to characterize VOCs (DCM solvent collection) from the pyrolyzed vegetation. Characterizing some of these chemicals can be difficult because of the detectability, high volatility, and/or lack of volatility of high mass chemicals found in light tars.
To compare the methods, a diverse group of seven gas samples (1 per plant species) were collected. Due to no replication, it is not possible to identify any effect due to plant species or any interaction between species and analysis method. From each of the seven samples, four extractions were made. Two were analyzed by GC and two were analyzed by GCxGC. NIST and Wiley libraries were used to determine the similarity matches using special algorithms that find the similarities between the libraries and the samples’ mass spectra. A higher similarity number indicates a better match between library and sample spectra. The number of peaks detected and the peak area percentage were summarized into four groups of compounds – hydrocarbons (CH), oxygenated organics (CHO), organic non-metals (ONM), and unknown (UNK). For each sample, method, and compound group the repeat extractions were averaged for further comparisons; Figs. 3 and 4 show the extraction measurement error.
Relationships between the 1D and 2D methods for the count data initially were evaluated based on paired comparisons assuming either Poisson or negative binomial distributions, while percent area was evaluated assuming beta distributions (Table 3). A generalized linear mixed model with a fixed effect for GC method and a random effect for sample was fit using PROC GLIMMIX from SAS/STAT® V14.1 (SAS Institute Inc., Cary, NC, 2015). However, to fully understand whether the relationships between chemical group identified by the two different methods (i.e., 2D vs. 1D) shift, the data, which are compositional (i.e., constrained), were further analyzed using compositional data analysis (CoDA) techniques (Aitchison 1986). This adjustment for counts to relative counts also partially removes effects due to the samples having different total counts; thus, models tend to indicate the random effect for sample is not significant. It should be noted that geometric means (not accounting for sample effects) resulted in similar values as given in Table 3, but the values in Table 3 are from models that include random sample effects.
The data averaged over method and compound group were analyzed as paired samples to test for differences in the number of peaks and the peak area percentages between the two methods. Because the data were compositional, they were transformed into log-ratios using an ilr-transformation, a standard approach for compositional data (Egozcue et al. 2003). The particular ilr-transformation used, called a balance, formed scientifically meaningful log-ratios. Specifically, the log ratios of the geometric means of 1) CH and ONM vs. CHO and UNK, 2) CH vs. ONM, and 3) CHO vs. UNK were calculated (Egozcue and Pawlowsky-Glahn 2005). Once transformed, the data are assumed to have a normal distribution and then can be analyzed using common statistical techniques, as has been done recently with smoke emissions (Weise et al. 2020). To compare the 1D method to the 2D method, multivariate paired Hotelling’s T2 statistic was used followed by confidence intervals on the differences (Hotelling 1931, 1951).
Observations are further made on the typical functional groups for each species. Although the species were not replicated for this study, the groupings from GCxGC are important to understand their reliability in the context of other research in the literature. Details and discussion are provided that further explain GCxGC plots and visually trending data relative to the wildfire smoke’s typical functional groupings for the two species Live Oak and Longleaf Pine. Molecular chemical maps (MCM) are developed and explained to aid groupings. Understanding how these MCMs are developed will aid in understanding the statistical summaries of the data when more samples can be evaluated.
RESULTS AND DISCUSSION
One-Dimensional versus Two-Dimensional Gas Chromatography
Table 3 shows a wide range in the number of peaks identified by the two methods. Due in part to the small sample size and diversity of samples, the sizes of the 95 percent confidence intervals were quite large. The multivariate paired Hotelling T2 test indicated that the composition of the white smoke as described by the relative number of peaks in the four chemical groups differed significantly between the GC methods (Prob > F3,4 = 0.00004). In contrast, the relative peak area percentages in the four chemical groups did not differ between the GC methods (Prob > F3,4 = 0.1258).
The difference plots with 95% Bonferroni confidence intervals (Fig. 2) can be used to determine where the estimated effect of GC method on the particular balances occurred. Recall only the number of peaks was significantly affected by the GC analysis method across the balances. In Fig. 2a, the 95% confidence intervals for the log-ratio of the CH compounds to ONM and the log-ratio of the sum of CH and ONM to the sum of CHO and UNK both narrowly include zero. It appears that the differences based on these two balances were correlated, which could be impacting the overall vs. individual tests. The average difference greater than 0 indicates that relatively more peaks were observed in the CH compared to ONM in the 1D method as opposed to the 2D method.
Table 3. Summary Statistics for two GC Methods Used to Determine Chemical Composition of “White Smoke” Samples Collected from Live Fuel Samples in a Cone Calorimeter. Estimated Means, Standard Errors and 95% Confidence Intervals are Given
The log-ratio of CHO to UNK is closer to 0, suggesting that their proportions were not impacted by GC method. While the 95% Bonferroni intervals for peak area differences (Fig. 2b) indicate a difference between GC method for the log-ratio of the sums, the overall Hotelling T2 indicated that the collection of log-ratios were not significant at the 0.05 level.
Fig. 2. Influence of GC analysis method on selected balances of chemical groupings. Values plotted are the average difference of log-ratios between 1D versus 2D for the indicated balance with 95% Bonferroni interval for that balance: (a) relative number of peaks, (b) area of peaks. Average difference less than 0 indicates either relatively less in the numerator or relatively more in the denominator of the balance.
It is important to remember that these data are relative and not absolute values. In summary, these limited data tend to suggest that the 1D GC analysis method yielded a relatively larger number of peaks in the CH category compared to the ONM category, or equivalently that the 2D GC analysis method yielded a relatively larger number of peaks in the ONM category compared to the CH category and that the log-ratio of sums of the chemical groups may be impacted as well.
Figures 3 and 4 show the raw data and suggest that there was a significant difference between 1D and 2D GC methods for the number of chemicals. The increase in diversity of chemicals will assist with the NMOC aerosols chemical gap, as discussed in the literature (Yokelson et al. 2013). The maximizing of tar production was achieved in the cone calorimeter under irradiance of 35 kWm-2 at a final temperature range of 400 to 600°C along with the full smoldering in air to nearly zero ash. These conditions provided a severe smoke profile containing primary pyrolysis products. These organic tars were separated into various categories and characterized in the 2D analysis.
The observance of nil ash weight after the tests was confirmed by the high number of ONM compounds. The observance of white smoke penetrating through the filters into the DCM solution seemed sufficiently rapid for avoiding the condensation reactions that form creosote in alternate smoke sampling schemes.
Fig. 3. 1D and 2D analysis of the number of peaks separated in DCM solution collected from seven vegetation pyrolysis white smoke tests. Chemicals are grouped by hydrocarbons (CH), oxygenated organics (CHO), unknowns (UNK), and organic non-metals (ONM) for the seven samples.
Fig. 4. The number of chemicals detected in DCM solution utilizing GC-TOFMS compared with GCxGC-TOFMS
Although the white smoke collected does deviate from other smoke sampling schemes, the overall purpose was to increase number of pyrolysis products likely to be produced in the field and to aid the goals of the project for numerical pyrolysis modelling.
Categorizing GCxGC Oxygenated Organics by Functional Groups
Figures 5 and 6 are refined separations of the raw data from CH and CHO categories for each sample. The endings of the chemical names were used to group items together to classify the chemical categories. Functional group (FG) classifications were -one (ketones), -ol (alcohols), -ane (alkanes), -ene (alkenes), -yne (alkynes), -ose (sugars), -al (aldehyde), and PAH, esters, acids, ide, furans, and miscellaneous (misc). Figure 5 and 6 show the number of chemicals and the peak area (PA%) of the FG classifications for the samples.
Fig. 5. The number of chemicals identified and their functional groups types for each vegetation sample from the pyrolysis analyses
LO, LLP, and WO were the vegetation with the greatest number of chemicals from CH and CHO, as shown in Fig. 5. The classes of chemicals identified with the most diversity were ketones, alcohols, and acids. The number of chemicals identified as PAHs were not numerous. This is to be expected because aromatization occurs at higher temperatures (Zhao et al. 2017; Sekimoto et al. 2018). Methoxy and phenolic hydroxyl were the major chemical groups found in lignin and (poly)phenolic compounds. The phenolic chemicals were grouped with alcohols, which could be a strong component of alcohols being one of the most abundant chemicals identified. The alkanes and alkenes classification included aromatic compounds. This is consistent with the literature (Amini et al. 2019). Tar from the pyrolysis of similar vegetations was composed mostly of ketones, alcohols, and acids. The number of such chemicals appear to be significant.
Fig. 6. Functional groups PA% in each vegetation sample
Figure 6 shows the PA% comparison for the seven plant samples. Within each vegetation, the PA% showed that there were major differences between the alcohols, ketones, and acids. Acids were the most abundant in six of the vegetations. SP, DP, and YP had the highest acid PA%, which was 61, 60, and 52%, respectively. LO, SB, and WO had the most abundant ketones at 24, 24, and 25 PA%, respectively. The 44% of acids reported in the literature was for carboxylic acids (Wang et al. 2017). The alkenes in LLP stood out among the vegetation at 24%. The contributing chemicals to this increase were aromatic compounds, and as reported later in the paper, there were several benzene compounds identified in this group. During pyrolysis, benzene rings were a major contributor to the formation of PAHs at higher temperatures. Another interesting observation was that there were three vegetation samples that had no alkynes detected. These were SP, LLP, and SB. This is especially noteworthy because alkynes are known to be a part of the chemical makeup of LLP extractives. During pyrolysis, many chemicals were degraded to form other chemicals. Alkynes are triple-bonded carbons; this bond structure is very unstable during pyrolysis. Benzene and alkynes were a part of a mechanism to form PAHs (Norinaga et al. 2014). The absence of alkynes and the abundance of alkenes (benzene) could contribute to the formation of PAHs. The 2D analysis provided a qualitative perspective, which can assist with characterizing the chemicals in the vegetation samples. The chromatograms produced from the data create a visual image of the locations of the unidentified and identified chemicals.
GCxGC Plots and Trending Data
GCxGC chromatograms are a three-dimensional view of the samples. Figure 7 contains different chromatogram styles for viewing the injected samples. The top is the surface plot, and the bottom is the contour plot with the 1D chromatogram to the left. The surface plot shows the 1st and 2nd retention times and the intensity of the peaks. Red represents the most intense peak for the different chemicals. The surface plot also gives a better view of coeluting chemicals. The ripples in the blue surface are trace chemicals that have been detected. From the contour view, any coeluting chemicals and any chemicals that would be hiding because of concentration can be seen. Yellow and red represent the intensity. Red is the most concentrated, and yellow is the least. The black dots are trace chemicals that have been either identified or detected. The 1D chromatogram shows high concentrated and coeluting peaks, making it difficult to identify (Dimandja et al. 2000; Dimandja 2004). These characteristics are what sets the 2D analysis apart from the 1D analysis.