Ali
Akherati‡
a,
Yicong
He
a,
Lauren A.
Garofalo
b,
Anna L.
Hodshire
b,
Delphine K.
Farmer
b,
Sonia M.
Kreidenweis
c,
Wade
Permar
d,
Lu
Hu
d,
Emily V.
Fischer
c,
Coty N.
Jen
e,
Allen H.
Goldstein
f,
Ezra J. T.
Levin
c,
Paul J.
DeMott
c,
Teresa L.
Campos
g,
Frank
Flocke
g,
John M.
Reeves
g,
Darin W.
Toohey
h,
Jeffrey R.
Pierce
*c and
Shantanu H.
Jathar
*a
aDepartment of Mechanical Engineering, Colorado State University, Fort Collins, CO, USA. E-mail: shantanu.jathar@colostate.edu
bDepartment of Chemistry, Colorado State University, Fort Collins, CO, USA
cDepartment of Atmospheric Science, Colorado State University, Fort Collins, CO, USA. E-mail: jeffrey.pierce@colostate.edu
dDepartment of Chemistry and Biochemistry, University of Montana, Missoula, MT, USA
eDepartment of Chemical Engineering, Carnegie Mellon University, Pittsburgh, PA, USA
fDepartment of of Environmental Science, Policy & Management, University of California Berkeley, Berkeley, CA, USA
gEarth Observing Laboratory, National Center for Atmospheric Research, Boulder, CO, USA
hDepartment of Atmospheric and Oceanic Science, University of Colorado Boulder, Boulder, CO, USA
First published on 15th June 2022
Wildfires are an important atmospheric source of primary organic aerosol (POA) and precursors for secondary organic aerosol (SOA) at regional and global scales. However, there are large uncertainties surrounding the emissions and physicochemical processes that control the transformation, evolution, and properties of POA and SOA in large wildfire plumes. We develop a plume version of a kinetic model to simulate the dilution, oxidation chemistry, thermodynamic properties, and microphysics of organic aerosol (OA) in wildfire smoke. The model is applied to study the in-plume OA in four large wildfire smoke plumes intercepted during an aircraft-based field campaign in summer 2018 in the western United States. Based on estimates of dilution and oxidant concentrations before the aircraft first intercepted the plumes, we simulate the OA evolution from very close to the fire to several hours downwind. Our model results and sensitivity simulations suggest that dilution-driven evaporation of POA and simultaneous photochemical production of SOA are likely to explain the observed evolution in OA mass with physical age. The model, however, substantially underestimates the change in the oxygen-to-carbon ratio of the OA compared to measurements. In addition, we show that the rapid chemical transformation within the first hour after emission is driven by higher-than-ambient OH concentrations (3 × 106–107 molecules per cm3) and the slower evolution over the next several hours is a result of lower-than-ambient OH concentrations (<106 molecules per cm3) and depleted SOA precursors. Model predictions indicate that the OA measured several hours downwind of the fire is still dominated by POA but with an SOA fraction that varies between 30% and 56% of the total OA. Semivolatile, heterocyclic, and oxygenated aromatic compounds, in that order, were found to contribute substantially (>90%) to SOA formation. Future work needs to focus on better understanding the dynamic evolution closer to the fire and resolving the rapid change in the oxidation state of OA with physical age.
Environmental significanceWildfires are an important atmospheric source of primary organic aerosol (POA) and precursors for secondary organic aerosol (SOA) at regional and global scales. However, there are large uncertainties surrounding the emissions and physicochemical processes that control the transformation, evolution, and properties of POA and SOA in large wildfire plumes. In this work, we develop a plume version of a kinetic model to simulate the dilution, oxidation chemistry, thermodynamic properties, and microphysics of organic aerosol (OA) in wildfire smoke. We find that dilution-driven evaporation of POA and simultaneous photochemical production of SOA are likely to explain the observed evolution of OA in wildfire plumes. Further, we show that there is rapid chemical transformation of wildfire smoke aerosol within the first hour after emission, driven by high concentrations of the hydroxyl radical. |
Wildfires directly emit particles that are dominated by primary organic aerosol (POA) with smaller amounts of black carbon (BC) and inorganic species (e.g., sulfate, nitrate, potassium).4,19,20 Wildfire POA has been shown to be semivolatile and evaporates with dilution.21–23 Furthermore, wildfire POA and the vapors in equilibrium with the POA (or semivolatile organic compounds, SVOC) are also known to be reactive.23–25 Wildfires emit OA precursors that include intermediate-volatility and volatile organic compounds (IVOC and VOC)26–28 and these oxidize in the atmosphere to form secondary organic aerosol (SOA).29 In addition to emissions of reduced hydrocarbons such as alkanes, aromatic, and biogenic VOCs,4 wildfires also emit oxygenated IVOCs and VOCs.23,27,28,30–34 These include oxygenated aromatic and heterocyclic organic compounds that have recently been shown to be important precursors for SOA formation.35–39 Although a lot has been understood about the composition, oxidation chemistry, and properties of OA and OA precursors from biomass burning emissions, especially in laboratory experiments, a detailed understanding of the physicochemical evolution of OA and OA precursors in real wildfire plumes is less well understood.
Recently, Hodshire et al.40 undertook a comprehensive review of wildfire OA data gathered in four laboratory campaigns and thirteen field campaigns performed over the past two decades. An analysis of these data suggested that while photooxidation of biomass burning emissions resulted in an enhancement in OA mass in laboratory experiments (mean = 1.25, median = 1.44, and interquartile range or IQR = 1.1–1.54),37–39,41–43 a similar enhancement was missing for the OA tracked in real wildfire plumes (mean = 1.1, median = 1.0, and IQR = 0.77–1.0).20,44–60 In addition, the OA measured within wildfire plumes closest to the fire (<1 h) was more oxidized (i.e., a higher oxygen-to-carbon (O:C) ratio for the OA) (mean = 0.36, median = 0.38, and IQR = 0.18–0.54) than the initial OA measured in laboratory experiments (mean = 0.23, median = 0.23, and IQR = 0.17–0.29), and the OA from wildfires exhibited a stronger O:C enhancement with photochemical age (mean = 1.82, median = 1.85, and IQR = 1.5–2.0) than the OA from laboratory fires (mean = 1.6, median = 1.5, and IQR = 1.3–1.8). Similar to the trends in OA O:C, a different enhancement was also observed for the ratio of two mass fragments (m44/m60) measured by the aerosol mass spectrometer between the field (mean = 4.5, median = 3.8, and IQR = 2.3–6.7) and the laboratory (mean = 3.5, median = 3.0 and IQR = 1.8–3.6). m44 refers to the mass fragment that is associated with the low-volatility organic compounds found in SOA61 while m60 refers to the mass fragment arising from fragmentation of primary emissions of anhydrous sugars such as levoglucosan.54
Based on this analysis, Hodshire et al.40 put forth four hypotheses to explain this field versus laboratory difference. First, they argued that the dilution of the wildfire plume could result in evaporation of POA and this evaporation could be approximately balanced by SOA production to result in a small change in the OA mass with photochemical aging. Since the condensing SOA is likely to have a higher O:C ratio than the evaporating POA, this POA-SOA swap would not affect the OA mass but instead produce an increase in the O:C ratio with aging. This hypothesis has been supported by theoretical calculations.62,63 Second, the authors hypothesized that because wildfire plumes are typically only sampled at least 15 minutes after emission (often longer), the plume measurements could have missed the OA evolution that happened prior to the first measurement. Early sampling poses safety issues as well as the prospect of measuring poorly mixed plumes. If rapid chemistry does indeed occur early on, this may explain both the higher O:C ratio measured in the field and the reduced propensity to form SOA after the aircraft first intercepts a plume. Third, they postulated that the field and laboratory differences could arise from differences in the emissions and chemistry of OA and OA precursors, attributed to differences in the fuel mixtures, burn conditions, combustion efficiency, and environmental conditions (e.g., temperature, relative humidity) and regimes (e.g., photolytic rates, NOX). Finally, they noted several experimental artifacts linked to laboratory experiments, including losses of OA and OA precursors to transfer ducts38 and walls of the environmental chamber,64–66 that could drive differences between laboratory and field measurements. For the most part, these hypotheses remain untested and will serve as the basis for this work.
Several studies have tested the first hypothesis above, i.e. that dilution-driven evaporation of POA is roughly balanced by SOA production in wildfire plumes.33,67 In one of these studies based on wildfire plumes sampled in the western United States (US), Palm et al.33 quantified the evolution of OA mass as a function of dilution and photochemical age. They found that plumes that had diluted but not yet undergone photooxidation showed evidence for POA evaporation, while plumes that had both diluted and undergone photooxidation showed evidence for replacement of the lost POA mass with SOA. However, Palm et al.33 did not account for the photochemical evolution prior to the first sample (i.e. <30 minutes since emission), nor did they account for the oxidation of all potential precursors (i.e., SVOCs) leading to SOA formation. Based on closure calculations, they argued that nearly 90% of the SOA came from oxidation of evaporated POA vapors with marginal contributions from other VOCs. To the best of our knowledge, no bottom-up, detailed models have been used to simulate and evaluate the OA mass and composition evolution in wildfire plumes.
In this study, we simulated the physicochemical evolution of OA in a subset of wildfire plumes sampled during a recent airborne field campaign based in the western US. The OA evolution was simulated using a kinetic model that accounts for the dilution, oxidation chemistry, thermodynamic properties, and microphysics of OA. A novel contribution of this work is that we incorporate estimates of dilution and photochemical age before the first airborne plume transect, to simulate the OA evolution from very close to the fire to several hours downwind. Model results of mass and composition were compared to plume measurements to study the contribution of precursors and processes to OA evolution.
Here, we focused on the evolution of the smoke emitted from the following four wildfires: Taylor Creek Fire (southwest Oregon), Sharps Fire (southern Idaho), Bear Trap Fire (eastern Utah), and Silver Creek Fire (northwest Colorado). The aircraft completed two distinct pseudo-Lagrangian sampling efforts for the Bear Trap Fire and each of these transect sets were analyzed separately. In total, the analysis focused on five transect sets. The location and dominant fuel(s) for each wildfire can be found in Table 1 of Lindaas et al.68 These wildfires were chosen because of their ideal sampling strategy, clear trends in the evolution of OA mass and O:C with time, and because all the required measurements were available to perform the proposed OA modeling.
Fire | OH concentration (molecules per cm3) | OH exposure (molecules h per cm3) | ||
---|---|---|---|---|
Before 1st transect | After 1st transect | At 1st transect | At final transect | |
Taylor Creek | 8.9 × 106 | 0.97 × 106 | 3.1 × 106 | 5.7 × 106 |
Sharps | 4.0 × 106 | 0.47 × 106 | 3.0 × 106 | 6.5 × 106 |
Bear Trap 1 | 3.0 × 106 | 0.57 × 106 | 2.7 × 106 | 4.6 × 106 |
Bear Trap 2 | 2.8 × 106 | 0.65 × 106 | 2.6 × 106 | 4.7 × 106 |
Silver Creek | 3.5 × 106 | 1.0 × 106 | 1.9 × 106 | 5.5 × 106 |
A suite of gas, particle, meteorological, and remote sensing instrumentation was deployed on the NSF/NCAR C-130 to comprehensively characterize the smoke emissions and its environment. Below, we briefly discuss the measurements and data products that relate to the scope of this work. Pressure and temperature were measured using base instrumentation available onboard the aircraft.68,69 CO was quantified using a Quantum-Cascade Tunable Infrared Laser Direct Absorption Spectrometer (QC-TILDAS, Aerodyne Research) and reported at 1 Hz.70 CO was assumed to be an inert tracer and used to model dilution with the background air and to develop convenient metrics to characterize the OA evolution with physical age. VOC mixing ratios were quantified using a Proton-Transfer-Reaction Time-of-Flight Mass Spectrometer (PTR-ToF-MS 4000, Ionicon Analytik) at 1 Hz.34 The PTR-ToF-MS quantified concentrations for 122 VOCs, including the most important SOA precursors that have been identified for biomass burning in previous work.37,38 A High-Resolution Aerosol Mass Spectrometer (HR-AMS; Aerodyne Research) measured the mass concentrations and composition of the sub-micron, non-refractory aerosol every 5 seconds.20 Measurements of aerosol mass concentrations and OA elemental ratios of H:C and O:C from the HR-AMS measurements were used here. A Single Particle Soot Photometer (SP2; Droplet Measurement Technologies) quantified BC mass concentrations every 10 seconds.20 Finally, aerosol size distributions were measured using a Passive Cavity Aerosol Spectrometer Probe (PCASP) at 1 Hz (ref. 71) and these were used to inform the shape of the initial aerosol size distribution (median and geometric standard deviation). We used the 10 seconds merge data for all the measurements mentioned above and conversions were performed when necessary to calculate concentrations at ambient conditions (instead of at standard temperature and pressure conditions). Links to the data are provided in the ‘Data Availability’ section.
We determined transect-average values for OA, non-OA (sulfate + nitrate + ammonium + chloride + BC), CO, and VOC concentrations as well as the OA O:C and the number size distributions. Following earlier published work with the WE-CAN data,20,33,68,69 we identified the plume-sampling periods by visually examining the CO time series for increases and decreases in mixing ratios. These averages were then corrected for background conditions by computing a universal background value for each measured species (e.g., CO, OA, benzene) using all data from outside the plume while the aircraft performed a transect set. This method assumed that the composition of the gas and aerosol species in the background air remained constant in time and space. We also calculated a normalized excess mixing ratio (NEMR) for OA and VOCs by ratioing its background-corrected concentration (in μg m−3 for OA and ppbv for VOCs) with the background-corrected CO (in ppbv). The background-corrected OA O:C was calculated as a ratio of the background-corrected molar concentrations for O and C; equations can also be found in Hodshire et al.72 For this calculation, we assumed the OA to be only composed of C, H, and O such that the molar concentrations of C, H, and O were determined from the available OA, O:C, and H:C data. While wildfire OA is likely to be composed of nitrogen-containing organic compounds, the low N:C values measured during WE-CAN (∼0.02)20 meant that accounting for nitrogen had a negligible impact on the reported OA mass concentrations and O:C and H:C values. Background corrections were performed separately for each wildfire plume transect set. We studied the sensitivity of the findings from this work to other ways in which the background values can be calculated and these are described later (Section 3.2).
Results for the first transect from the five transect sets analyzed in this work are shown in Fig. 1, as an example. The ratio of VOC NEMRs exhibited an inverse relationship with kOH. The likely explanation for this inverse trend was that the oxidation chemistry prior to this transect resulted in stronger depletion of VOCs with a higher kOH and weaker depletion of VOCs with a lower kOH. An alternate explanation is that the ratio of VOC NEMRs was biased lower for the higher kOH species compared to the lower kOH species on account of differences in fuel and burn conditions as well as the timing of the emissions, but we are not aware of why there might be systematic changes in the emissions with the species kOH. Assuming that chemistry prior to the first transect explains the trends in Fig. 1, the slope in the ratio of VOC NEMRs with kOH was fit to determine an effective OH exposure for the time period before this transect. We should note that we did not assume that the field and laboratory VOC NEMRs at the time of emission were the same, which would require setting the intercept for the fit to 0, which we did not do. Rather, we assumed that there were no systematic kOH-dependent differences in emissions between the field and laboratory cases. This exercise was repeated for all transects to calculate a time-varying OH exposure estimate for all transect sets, results for which are shown in Fig. 2. A line was fit through the OH exposure data to calculate an average OH concentration for the wildfire plume after the first transect. We assumed that the OH exposure changed linearly from zero, close to the fire, to the OH exposure estimated at the first transect, which resulted in a separate, time-invariant OH concentration estimate for the time period before the first transect.
Fig. 1 Logarithm of the ratios of the VOC NEMR in the first transect of the wildfire plume to the average value of the VOC NEMR measured in laboratory experiments plotted against kOH. Results are shown for the (a) Taylor Creek, (b) Sharps, (c) Bear Trap 1, (d) Bear Trap 2, and (e) Silver Creek Fire transect sets. Solid red lines represent the linear fit to the data and the red bands capture the standard error. The table lists statistics based on the linear fit including the OH concentration and exposure. The time-varying OH exposures for all five transect sets are shown in Fig. 2 and the OH concentrations are listed again in Table 1. |
The OH concentrations before and after the first transect are tabulated for the five transect sets in Table 1. For the five transect sets and four fires analyzed, the total mean OH exposure varied between 4.6 × 106 (Bear Trap 2) and 6.5 × 106 (Sharps) molecules h per cm3 to produce 3.1 to 4.4 hours of photochemical aging at an average OH concentration of 1.5 × 106 molecules per cm3. The OH concentrations varied between 2.8 × 106 (Bear Trap 2) and 8.9 × 106 (Taylor Creek) before the first transect and between 0.47 × 106 (Sharps) and 1.0 × 106 (Silver Creek) beyond the first transect. This analysis approach suggests that OH concentrations, on average, were a factor of six larger before the first transect compared to after the first transect suggesting that photochemical oxidation was as relevant to the wildfire plume in the near field (between the point of emission and the first transect) as it was for the far field (beyond the first transect) after accounting for differences in the physical age before and after the first transect; the first transect was estimated between 21 and 56 minutes of the physical age of the wildfire plume.
We assessed our estimates of OH concentrations in the wildfire plumes after the first transect and found them to be consistently lower than those estimated in previous work. For instance, our estimates for OH concentration in the Taylor Creek and Bear Trap wildfire plumes after the first transect (0.97 × 106 and ∼0.6 × 106 molecules per cm3, respectively) were a factor of eight and five smaller than previous OH estimates for the same wildfire plumes calculated from the observed decay for a handful of VOCs in the core of the plume (∼8 × 106 and ∼3 × 106 molecules per cm3, respectively).69 Similarly, our estimates of OH concentration after the first transect, across all five transect sets, were a factor of five to ten smaller than those estimated in Mexican and African biomass burning plumes44,48,53 where OH concentrations were calculated using techniques similar to those outlined in Juncosa Calahorrano et al.69 We tentatively argue that the historical estimates of OH concentrations in wildfire plumes are larger and different than those measured here because earlier work has only used a few species (<5) to calculate OH concentrations. We hypothesize that our approach, which uses a much larger number of species (106), provides a more robust estimate for mean OH concentrations in the plume. However, this discrepancy in OH concentrations needs to be investigated in detail in future work. Finally, our estimated OH concentrations were found to be lower than the average OH concentrations typically found in urban environments and regional backgrounds (∼1.5 × 106 molecules per cm3).73 In this specific urban and regional comparison, the low OH concentrations in the wildfire plume might be a result of the large OH sinks present therein as well as the low photochemical activity expected in optically dense plumes.74
We note that the OH estimates used in this work are uncertain. The OH estimates are likely to be unbiased if the field and laboratory emissions ratios have no systematic relationship with kOH. In case there is a relationship and that this relationship is different between the field and laboratory, our OH estimates could be biased lower or higher. Hence, there is no clear way in which we could argue that our OH estimates bound the range of expected values in a wildfire plume. Since the OH estimates were developed based on the observed decay of all species including reactive oxygenated VOCs, any chemical production of oxygenated species should bias our OH estimate lower. Hence, based on the inclusion of the reactive oxygenated VOCs alone, our OH estimates potentially present a lower bound estimate. Although there is limited evidence, Coggon et al.75 have found in laboratory experiments performed on biomass burning smoke that there was little to no production of reactive oxygenated VOCs, such as phenol and furan, relative to these species' primary emissions. When OH was calculated from hydrocarbon VOCs measured by the PTR-ToF-MS to eliminate the influence from including oxygenated VOCs, the inverse relationship was weakened and produced OH concentrations that were at least a factor of 2 lower than those listed in Table 1, both before and after the first transect. We also calculated OH concentrations beyond the first transect by using the VOC NEMRs at the first transect as the reference instead of using the VOC NEMRs from the laboratory. The concentrations so calculated and limited to the time period beyond the first transect were found to be only slightly higher to those listed in Table 1 (0.67 × 106–1.2 × 106 molecules per cm3) but still lower than those in historical studies mentioned earlier. The VOC NEMRs at the first transect, by definition, cannot be used to determine the OH concentrations prior to the first transect. Given the uncertainty, we studied the sensitivity of the findings from this work to the OH concentrations before and after the first transect and these are described later (Section 3.2).
The SOM tracks the chemical evolution of OA and its precursors using a two-dimensional, carbon (NC) and oxygen (NO) number grid. The properties of each model species (e.g. reactivity (kOH), volatility (c*)) are parameterized based on their NC and NO. The SOM has six adjustable parameters that govern the oxidation chemistry and thermodynamic properties of the model species: (i–iv) pf,1–pf,4, the yields of four functionalized products that add one, two, three, and four oxygen atoms to the carbon backbone, respectively; (v) mfrag, the parameter that characterizes the fragmentation probability (Pfrag); and (vi) ΔLVP, the decrease in the c* of the model species per addition of an oxygen atom. The particle-phase species in SOM are tracked in 36 TOMAS size sections that span dry diameters between 3 to 10000 nm. TOMAS simulates coagulation between size sections and the kinetic condensation/evaporation of mass between the particle and vapor phases for all SOM model species. The SOM-TOMAS model also accounts for formation of highly oxygenated organic molecules (HOMs) and formation/dissociation of oligomers.84 The SOM-TOMAS model was also updated recently with the diffusive-reactive framework described in Zaveri et al.87 to model the influence of phase state on the kinetic gas/particle partitioning of OA.84 Although SOA precursors found in biomass burning emissions (e.g., monoterpenes) are known to form HOMs88,89 and oligomers90,91 and certain biomass burning particles can be viscous, neither HOM nor oligomer formation was modeled and the OA was assumed to be liquid-like with a diffusion coefficient of 10−10 m2 s−1. These assumptions surrounding HOMs, oligomers, and phase state will need to be examined in future work.
For each wildfire plume, the SOM-TOMAS model was used to simulate the OA evolution from the time just after emission up to the last measured transect. While this should include both vertical plume rise and horizontal plume transport after reaching the equilibrium height, we do not explicitly model the change in the pressure and temperature of the air parcel during vertical plume rise. The background-corrected CO concentrations past the first transect were first fit to an exponential function and the CO values including those extrapolated to t = 0 were then used to determine dilution of the wildfire plume with background air starting at t = 0. Since we only used background-corrected values in this work, we assumed the background air to be free of any trace gases and particles. This assumption affects non-linear processes, such as equilibrium partitioning and coagulation; however, because background concentrations are generally much lower than in-plume concentrations, these are expected to have a minimal effect.63 While the pressure and temperature values changed modestly between wildfire plumes, they were found to be in a relatively narrow range within each individual wildfire plume (Fig. S2†). Hence, an average pressure and temperature value was used for the entire wildfire plume. Model predictions and measurements of concentrations and mixing ratios in this work were expressed at the plume pressure and temperature.
The VOC NEMRs at the first transect were extended to t = 0 by correcting for the OH exposure before the first transect (see Table 1 for transect-set-specific OH exposures). These field VOC NEMRs at t = 0 and the first transect were aggregated by precursor class and compared to the VOC NEMRs measured by Koss et al.28 during a recent laboratory campaign. Results of this comparison, combined over all five transect sets and four wildfires, are shown in Fig. 3. Results when separated by transect set were similar to those presented in Fig. 3 and, hence, are not shown or discussed. The VOC NEMRs for the five SOA precursor classes at the first transect were substantially lower than those measured in the laboratory. The median values at the first transect were 74%, 58%, 28%, 15%, and 75% lower for the biogenic, oxygenated aromatic, heterocyclic, aromatic hydrocarbons, and alkane classes, respectively, compared to the laboratory values. The largest differences were observed for the most reactive SOA precursors with respect to OH (e.g., oxygenated aromatics) and the smallest differences were observed for the least reactive SOA precursors with respect to OH (e.g., aromatic hydrocarbons), a feature that was leveraged to determine the OH concentrations and exposure in Section 2.2. Even after correcting for OH exposure prior to the first transect, the VOC NEMRs for three of the organic classes (i.e., biogenic, oxygenated aromatic, and alkane) were still modestly lower compared to those measured in the laboratory. The median values at t = 0 were 18%, 38%, and 61% lower for the biogenic, oxygenated aromatic, and alkane classes, respectively, compared to the laboratory values. The median values were 32% and 15% higher for the heterocyclic and aromatic hydrocarbon classes, respectively, compared to the laboratory values.
There are two important implications of the differences in the ratios of VOC NEMRs shown in Fig. 3 and described in the previous paragraph. First, field emissions of SOA precursors (those at t = 0) are lower than those measured in laboratory experiments, for at least a few of the important organic classes. It is expected then that the SOA precursor mixture would be relatively less potent in forming SOA in wildfire plumes than what has been observed in laboratory experiments.37–39,41–43 Second, a significant fraction of the SOA precursor emissions are depleted between t = 0 and the first transect (between 15% and 74%). This implies that not only is SOA being actively produced in the near field after emission (<1 hour) but also that there is significantly reduced potential for continued production of SOA beyond the first transect. These comparisons partly explain the field versus laboratory differences outlined by Hodshire et al.,40 for the OA evolution after the first transect.
We used the detailed speciation data for POA and SVOCs from Jen et al.25 to create a substitute volatility distribution, which was then used to study the sensitivity in model results to the volatility distribution assumed for POA and SVOCs. Jen et al.25 used two-dimensional gas chromatography with time-of-flight mass spectrometry (GC × GC-ToF-MS) on derivatized samples to measure the detailed chemical composition of POA and SVOC emissions from 29 laboratory-generated fires performed on 24 different fuels. They measured ∼150 unique species in the sample collected on quartz filters and these amounted to between 10 and 65% of the total POA/SVOC mass. The averaged and normalized volatility distribution and mass distribution in carbon-oxygen space for the POA and SVOC emissions is shown in Fig. S4(a) and (b),† respectively. The volatility distribution was constructed by binning the ∼150 species by c*, which was estimated for each species using EPISuite 4.11,98 a numerical model that estimates physical and chemical properties of pure substances.
Jen et al.25 found that the SVOC emissions were dominated by sugars, phenols, and other complex organic compounds (e.g., terpenoids, heterocyclics), organic compound classes also measured by the PTR-ToF-MS. In theory, there should be little if any overlap between the SOA precursors quantified by the PTR-ToF-MS and GC × GC-ToF-MS since the PTR-ToF-MS exclusively sampled gas-phase compounds with a c* higher than 105 μg m−328 and the GC × GC-ToF-MS measurement was performed on quartz filters that generally trap gas- and particle-phase compounds lower than a c* of 105 μg m−3.22 In this work, we assumed that the SVOCs and VOCs represent a mutually exclusive set of SOA precursors despite similarities in the compound classes within these two categories. SVOCs have been hypothesized to be important precursors of SOA formation from biomass burning emissions,20,33 but there are few laboratory datasets that can be called upon to inform the oxidation chemistry of SVOCs in our model. In addition to Jen et al.,25 several studies have measured semivolatile multi-ring aromatics in biomass burning emissions and they are expected to serve as SOA precursors.23,35,99 Hence, in the absence of any model SOA precursors for which laboratory data are available, we simulated the oxidation chemistry of SVOCs in the base configuration of the SOM-TOMAS model assuming that they have a similar potential to form SOA as multi-ring aromatics (i.e., naphthalene). However, given the uncertainty in this assumption, we performed sensitivity simulations where we modeled the oxidation chemistry of SVOCs as oxygenated aromatics (i.e., phenol), heterocyclics (i.e., 2-methylfuran + dimethylfuran), linear alkanes (i.e., n-dodecane), or biogenics (i.e., α-pinene) with the species in parentheses used as the surrogate to model SOA formation (Section 3.2).
(1) In addition to the base simulations (SV-POA + FullChem), we performed a systematic sequence of simulations where we tested assumptions about POA volatility and SVOC and VOC chemistry: (i) non-volatile POA and no chemistry (NV-POA + NoChem), (ii) non-volatile POA and VOC chemistry (NV-POA + Chem), (iii) semivolatile POA and no chemistry (SV-POA + NoChem), and (iv) semivolatile POA and SVOC chemistry (SV-POA + NoVOCChem). Results from these simulations are presented in Fig. 4 and 5.
(2) Simulations were performed to assess the sensitivity in model results to the treatment of POA and SVOCs: (i) we simulated the oxidation chemistry for SVOCs using the SOM parameters for oxygenated aromatics (i.e., phenol), heterocyclics (i.e., 2-methylfuran + dimethylfuran), linear alkanes (i.e., n-dodecane), or biogenics (i.e., α-pinene), with the species in parentheses used as the surrogate to model SOA formation, (ii) we used the POA volatility distribution determined from the work of Jen et al.25 in combination with the five surrogate species to model the oxidation chemistry of SVOCs, and (iii) we performed a thousand Monte-Carlo simulations where we randomly chose eight species in the SOM grid to represent the POA and SVOC mass and fit a mass distribution with these eight species that reproduced the POA volatility behavior observed by May et al.22 Results from these simulations are presented in Fig. 6, S9 and S10.†
(3) Simulations were performed to assess the sensitivity in model results to the assumed OH concentrations, ‘first’ transect, and background correction: (i) we assumed a low or high OH concentration over the entire wildfire plume with the low value based on the OH estimate after the first transect and the high value based on the OH estimate before the first transect (Table 1), (ii) we assumed a time-varying OH concentration that was determined by fitting a power function through the OH exposure data presented in Fig. 2, (iii) we assumed a constant OH concentration of 1.5 × 106 molecules per cm3 that is commonly used to reflect average OH concentrations in polluted environments, (iv) instead of using the transect closest to the fire as the first transect, we assumed the second or third closest transects to be the ‘first’ transect to perform the simulations, and (v) measurements and model inputs were calculated by performing background corrections by using the transect-specific background concentrations or assuming the concentrations upwind of the fire to be representative of the true background. Results from these simulations are presented in Fig. 7, S11 and S12.†
For the simulations with non-volatile POA and no oxidation chemistry (NV-POA + NoChem), the model reproduced observations of OA NEMR (Mean Bias Error (MBE) = -0.009, Mean Absolute Error (MAE) = 0.011; μg m−3 ppbv−1) but, by definition, produced no change in the OA O:C (MBE = −0.060, MAE = 0.062). With the oxidation chemistry turned on to produce SOA from VOCs (NV-POA + Chem), the model only produced a marginal increase in the OA NEMR compared to the NV-POA + NoChem simulation and, hence, reproduced observations of OA NEMR (MBE = −0.007, MAE = 0.011; μg m−3 ppbv−1). The NV-POA + Chem simulation produced a slight increase in OA O:C with physical age compared to the NV-POA + NoChem model but still significantly lower than the observed increase in OA O:C (MBE = −0.055, MAE = 0.057). The relatively small increase in OA NEMR and OA O:C can be explained by the marginal amounts of SOA formed from the VOC mixture beyond the first transect. As noted in Section 2.4, this VOC mixture was substantially depleted in important SOA precursors by the first transect. For the simulation with semivolatile POA with no oxidation chemistry (SV-POA + NoChem) and the simulation with semivolatile POA with oxidation chemistry for SVOCs alone (SV-POA + NoVOCChem), the model appeared to underestimate the OA NEMR (MBE = −0.035 and −0.028, MAE = 0.035 and 0.028, respectively; μg m−3 ppbv−1) with a substantial increase in OA O:C. The predicted decrease in OA NEMR compared to the NV-POA + NoChem model stemmed from the evaporation of POA with dilution that was only partly recovered through SOA formation in the SV-POA + NoVOCChem model. Results from the SV-POA + NoChem simulation suggested that POA evaporation alone could explain the increase in the OA O:C with physical age as the lower-volatility material left in the particle phase after evaporation had a relatively higher O:C than the semivolatile material that had evaporated (Fig. S3(b)†).
The base simulation that assumed a semivolatile POA and oxidation chemistry for both SVOCs and VOCs (SV-POA + FullChem) underestimated the OA NEMR by 15% compared to observations (MBE = −0.024, MAE = 0.024; μg m−3 ppbv−1) but produced a large increase in OA O:C consistent with the observations (MBE = −0.019, MAE = 0.032; μg m−3 ppbv−1). The base simulation predicted a lower increase in O:C compared to the SV-POA + NoChem simulation because the SOA being formed in the base simulation had a lower O:C than the remaining POA. Overall, two of the simulations appeared to be the most consistent with observations of OA NEMR but significantly underestimated the observations of OA O:C (NV-POA + NoChem and NV-POA + Chem), and two of the simulations came close to reproducing the increase in OA O:C with physical age but underestimated the observations of OA NEMR (SV-POA + NoChem and SV-POA + NoVOCChem). The base or SV-POA + FullChem simulation seemed to offer a balanced comparison with observations where both the OA NEMR and O:C were only slightly underestimated.
Results from simulations for the remaining fires and transect sets (Sharps, Bear Trap, and Silver Creek) are shown in Fig. 5(a–h), where the relative trends across the five different simulations were similar to those presented for Taylor Creek in Fig. 4. Therefore, the differences in these simulations are not described further. However, as discussed below, the absolute performance of the base simulation for these other transect sets was mixed. The base simulation could not reproduce the initial increase and later decrease in observed OA NEMR for the Sharps and Silver Creek Fires. Both of these fires were sampled by the aircraft much faster than the physical age (Fig. S1†), suggesting that the measurements may partly reflect changes in emissions rather than those from their physicochemical evolution. The base simulation produced a mildly increasing OA NEMR for the Sharps Fire (MBE = −0.004, MAE = 0.006; μg m−3 ppbv−1) and a sharply increasing OA NEMR for the Silver Creek Fire (MBE = −0.004, MAE = 0.012; μg m−3 ppbv−1). Unlike the comparison for the Taylor Creek Fire, the base simulation appeared to reproduce the relatively constant observations of OA NEMR for the Bear Trap Fire transects. The base simulations produced a different trend in the modeled OA NEMR with physical age across the five different transect sets because presumably, in each of these transect sets, there were differences in the dilution rate, environmental conditions, and absolute concentrations of the OA, OA precursors, and oxidants. The base simulation consistently underestimated the change in OA O:C with physical age, with the strongest comparison for the Bear Trap 1 Fire (MBE = −0.044, MAE = 0.044) and the weakest comparison for the Silver Creek Fire (MBE = −0.087, MAE = 0.087).
Taken together, we draw the following conclusions from the model-measurement comparisons presented in Fig. 4 and 5. First, a non-volatile and non-reactive treatment of POA, regardless of the inclusion of SOA produced from VOC oxidation, is unlikely to explain the combined observations of OA NEMR and O:C. Second, POA evaporation with dilution alone can potentially explain the change in the OA O:C with time but results in loss of OA mass (decrease in OA NEMR) that is not consistent with observations. Third, we argue that POA evaporation with dilution and SOA formation from both SVOCs and VOCs, as captured by the base model, best explains the trends in both the OA NEMR and OA O:C across the five transect sets and four wildfires: average MBE = −0.007 and MAE = 0.011 (μg m−3 ppbv−1) for OA NEMR and average MBE = −0.046 and MAE = 0.055 for OA O:C across all five transect sets. We note that the model treatment in the base simulations closely reflects our updated understanding of the POA and SOA system from wildfire emissions.
Model predictions of the normalized OA composition with physical age from the base simulations are shown in Fig. 4(c) for the Taylor Creek Fire and in Fig. 5(i–l) for all other transect sets while model predictions of the POA-SOA split with photochemical age for all transect sets are shown in Fig. S7.† There are several interesting features to note. The OA composition changed rapidly starting at t = 0 from the dilution-driven evaporation of directly emitted POA and the SOA produced from the oxidation of SVOCs and VOCs. By the first transect, the average OA across all five transect sets, 21 to 56 minutes after emission, was 73% POA and 27% SOA. As the base-simulation-predicted OA NEMR was found to universally increase during the period before the first transect for all five transect sets, the changes in OA composition were largely driven by SOA condensation rather than POA evaporation. The rapid SOA production was facilitated by the higher OH concentrations experienced before the first transect (2.8–8.9 × 106 molecules per cm3; see Table 1). Generally speaking, the observed and base-simulation-predicted OA NEMRs did not vary much past the first transect for any of the transect sets. Despite that fact, the base simulations predicted a modest change in the OA composition with physical age after the first transect, suggesting a roughly equal replacement of POA with SOA. Over all five transect sets, POA continued to dominate the total OA mass beyond the first transect (>45%) but there was continued production of SOA over this time period. By the last transect, the SOA contribution to the total OA varied between 30% for the Taylor Creek Fire and 56% for the Silver Creek Fire.
The base simulations predicted that the majority of the SOA was formed from the oxidation of SVOCs, heterocyclics, and oxygenated aromatics, in that order. The contribution of the different VOC classes to SOA formation was similar between the different transect sets, although oxygenated aromatics contributed much more to SOA formation in the Silver Creek Fire than in the other Fires. On average, these three precursor classes accounted for 45, 25, and 21% of total SOA and 18, 11, and 9% of the total OA by the last transect. Heterocyclics and oxygenated aromatics have already been implicated as important SOA precursors in laboratory experiments performed on biomass burning emissions37–39,100 and our results here confirm their relevance for wildfire plumes as well. Biogenic VOCs were found to be significantly less influential compared to the precursor classes just discussed where they accounted for less than 5% of total SOA and 2% of total OA by the last transect.
These simulations provide model-based evidence for dilution-driven evaporation of POA mass being replaced by SOA mass formed from the oxidation of SVOCs and VOCs in wildfire plumes during WE-CAN. This conclusion is consistent with the theoretical findings of Bian et al.62 and Hodshire et al.,63 who showed that POA evaporation can be approximately replaced with SOA condensation under certain conditions pertaining to the fire size, background concentrations, and atmospheric stability. The base simulations predicted a mean POA-SOA split of 59–41% by the last transect over all five transect sets. These model-predicted POA-SOA splits agreed well with the theoretical findings of Hodshire et al.,63 who predicted a POA-SOA split of ∼50–50% for 1 km2 fires and ∼75–25% for 100 km2 fires, and the findings of Palm et al.,33 who analytically determined a maximum POA-SOA split of 66–33% for OA measured over several wildfire plumes during WE-CAN. Note that Palm et al.33 estimated their maximum POA-SOA split assuming that the OA at the first transect was exclusively POA. In contrast to the precursor-resolved findings discussed above, Palm et al.33 proposed that SVOCs were responsible for nearly 90% of the SOA formed within the plume. However, their approach would overestimate the SVOC contribution to SOA because the underlying closure calculation subtracted the SOA estimated from non-SVOC precursors from the total SOA formed.
On using the substitute volatility distribution based on the speciation data of Jen et al.,25 the model produced a relatively larger spread in the OA NEMR predictions compared to the sensitivity result discussed above. The use of an oxygenated aromatic (i.e., phenol/guaiacol) and biogenic (i.e., α-pinene) surrogate seemed to agree the best and the use of a heterocyclic (i.e., 2-methylfuran + dimethylfuran), multi-ring aromatic (i.e., naphthalene), and linear alkane (i.e., n-dodecane) surrogate seemed to agree the least with the OA NEMR observations. All simulations produced a flat response in the OA O:C with physical age despite a gradual change in the POA-SOA split and dilution-driven evaporation of the semivolatile material. This flat response, which was inconsistent with the observed change in the OA O:C, was primarily from the lower-volatility material that was left in the particle phase being less oxidized than the higher volatility material that had evaporated (Fig. S4(c)†). The substitute volatility distribution resulted in a larger POA-SOA split (average of 83–17%) compared to the first set of sensitivity results (average of 75–25%) likely because a smaller fraction of the fresh POA mass was lost to evaporation from the use of a less volatile volatility distribution (Fig. S4(a) and (d)†). The mixed comparisons for OA NEMR and O:C suggest that model inputs determined from speciation data hold promise but might be limited because, in this specific case, the speciated compounds only represented a fraction (10–65%) of the total POA + SVOC mass.
Finally, we performed a thousand Monte-Carlo simulations where we randomly specified the mass distribution of the POA + SVOC mass in the SOM grid while ensuring that this mass distribution reproduced the POA volatility behavior observed by May et al.22 (Fig. S8†). The iterations produced a relatively narrower spread in the OA NEMR predictions compared to the previous sensitivity simulations but they all seemed to underestimate the observed trends. By the last transect, the predicted OA NEMR was 17 to 29% lower than the average observed OA NEMR. Compared to the OA NEMR, there was a much larger spread in the predicted OA O:C with iterations predicting a decrease in O:C with physical age at the one end (from 0.41 to 0.31) to reproducing the observed increase at the other end (0.38 to 0.54). It appears that the mass distribution of POA + SVOC in the SOM grid had a significant, non-linear influence on POA evaporation, SOA production, and subsequently on the OA O:C evolution with photochemical age. We note that the mass distribution of POA + SVOC in the SOM grid across the thousand simulations was always constrained to observations of POA volatility but had enough degrees of freedom to produce a substantial spread in model predictions. The iteration that resulted in the largest absolute OA NEMR at the end of the simulation and then the largest increase in OA O:C was chosen to represent the base simulation results presented in Fig. 4 and 5. In other words, the distribution of POA + SVOC mass in the SOM grid for the base simulation was chosen from amongst those used in the Monte-Carlo simulations that produced the most optimum comparison against measurements of OA NEMR and O:C for the Taylor Creek Fire. There is some indirect evidence for this POA + SVOC mass distribution in that the O:C dependence with OA mass loading seemed to agree qualitatively with a subset of laboratory and field observations of biomass burning OA (not shown).21,97
Results from simulations performed to assess the sensitivity to the OH estimates are presented in Fig. 7(a, d and g). The use of a power function fitted to the OH exposure data to determine OH concentrations (OH-power fit) produced results that were slightly higher compared to those from the base simulation. Similarly, if we assumed that the lower OH concentration after the first transect was also relevant to the time period before the first transect (OH-low; 0.97 × 106 molecules per cm3), the model predictions of OA NEMR and O:C were slightly higher than those from the base simulation. When using a constant OH concentration of 1.5 × 106 molecules per cm3 (OH-ambient) or the higher OH concentration from before the first transect (OH-high; 8.9 × 106 molecules per cm3) for the entire evolution, the model predicted a higher OA NEMR compared to predictions from the base, OH-power fit, and OH-low simulations. This was because the OH concentrations after the first transect in both of these instances were larger than those used in the base, OH-power fit, and OH-low simulations and these higher OH concentrations, which were ∼50% larger in OH-ambient and a factor of ∼10 larger in OH-high, promoted SOA formation.
For the OH-high simulation, the increase in OA NEMR was found to be relatively consistent with the evolution in the observations indicating that the OH concentrations may continue to be elevated even after the first transect. However, the OH-high simulations overestimated the OA NEMR compared to the observations for the other Fires (Fig. S11†), where the base, OH-power fit, and OH-low simulations produced results that were more in line with the observations. As the Taylor Creek Fire dataset is the most Lagrangian amongst all Fires, the OH sensitivity simulation results presented here provide some evidence that our OH concentration estimates after the first transect (Table 1) may be biased low and would need to be revised in future work to be consistent with the higher OH concentrations estimated in earlier work (see Section 2.2 for a longer discussion). Interestingly, a higher OA NEMR in the OH-high simulation did not change predictions for OA O:C presumably because the additional SOA formed had an O:C similar to the existing OA's O:C. Overall, the model predictions appeared to be somewhat sensitive to the OH concentration inputs that produced a significant spread in the OA NEMR and POA-SOA splits but not so much in the OA O:C.
Qualitatively, the modeled trends in the OA NEMR and OA O:C were not very different when we assigned different transect sets (i.e., ‘2nd Transect’, ‘3rd Transect’) to be the ‘1st transect’ where model predictions of OA NEMR and OA O:C were anchored to observations; results are presented in Fig. 7(b, e and h). There appeared to be some tradeoff in the model-measurement comparisons for OA NEMR and OA O:C with the transect chosen. For instance, anchoring the model predictions to the information at the third transect seemed to produce better agreement with observations of OA NEMR but underestimated observations of OA O:C. The opposite was found to be true when anchoring the model predictions to the information at the second transect. Regardless of the variability in the OA NEMR and O:C, the POA-SOA splits were nearly identical between the three simulations. The use of ‘Transect-Specific’ or ‘Upwind’ data to calculate values to perform the background corrections did not seem to have any significant effect on model predictions; results are presented in Fig. 7(c, f and i).
First, fresh emissions of POA were found to be similar in magnitude to the sum of SOA-forming SVOCs and VOCs. For reference, for mobile sources, SOA precursor emissions are easily an order of magnitude larger than those for POA.29,101 This means that direct emissions of POA are likely to be an important constituent of smoke aerosol downwind of the fire, even as some fraction of it is lost to evaporation and surface reactions (not modeled in this work) and SVOCs and VOCs oxidize to form SOA and add to OA mass.
Second, up to half of the SOA precursor mass was rapidly oxidized to form SOA early on, which primarily resulted in an increase in the OA NEMR and O:C prior to the first transect. This result implies that a substantial fraction of the total SOA is actually formed very close to the fire and that the OA measured on the first aircraft transect is likely to be enhanced compared to the fresh OA emissions measured in laboratory experiments, after accounting for differences in gas/particle partitioning at different OA mass loadings. Furthermore, the rapid evolution of the OA system close to the fire is bound to alter the chemical, microphysical, and optical properties of smoke aerosols early on and confound comparisons of aerosol measurements at the first transect with similar measurements made on fresh emissions in laboratory environments.
Third, for most of the modeled and measured transect sets, except for the model prediction for the Silver Creek Fire, the OA NEMR did not vary much between the first and last transects, but there was a gradual change in the modeled OA composition with POA evaporation and SOA formation from SVOCs and VOCs. This is consistent with both theoretical and analytical findings in Bian et al.,62 Hodshire et al.,63 Palm et al.,33 and Liang et al.67 The change in composition implies that while the OA mass may remain constant, its atmospheric properties will continue to evolve with physical age. Furthermore, as the SOA precursors were heavily depleted and the OA mass concentrations at the last transect (5–30 μg m−3) were only marginally larger than the background concentrations (1–20 μg m−3), we postulate that the OA NEMR and O:C are unlikely to change dramatically with additional photochemical aging, any different than the changes experienced by background aerosol.
Fourth, by the last transect, SVOCs and VOCs contributed about equally to SOA formation in our wildfire plumes. The dominant SVOCs and VOCs contributing to SOA formation are very likely to be oxygenated organic compounds (e.g., sugars, heterocyclics, oxygenated aromatics), classes that are not explicitly included or are represented too coarsely in emissions inventories and chemical mechanisms, part of atmospheric models. Hence, the representation of SOA formation, as studied here, needs to be reflected for biomass burning sources in atmospheric models.
Model results were found to be moderately sensitive to the treatment for POA and SVOCs. Hence, continued work to fully speciate the POA and SVOC mass to inform the volatility properties of POA and to identify surrogate species to model the oxidation chemistry of SVOCs will likely lead to improvements in model predictions. A point of contention for SVOCs is that they have not been explicitly considered when studying the SOA formation from biomass burning emissions in laboratory experiments.37,38 Akherati et al.38 observed that lower-volatility SOA precursors, especially in the SVOC range, were susceptible to loss in transfer ducts used to direct smoke emissions into environmental chambers. If this is indeed true, this might be one reason why SVOCs remain highly relevant for wildfire plumes but may not have been for laboratory experiments. Regardless, the chemical composition and oxidation chemistry of SVOCs relevant to SOA formation needs to be studied in the future.
We acknowledged a significant discrepancy in OH concentrations in the wildfire plume based on techniques used in this work and OH concentrations estimated in earlier work. In addition, model predictions were found to be somewhat sensitive to the OH concentrations assumed in the wildfire plume. Hence, ongoing and future work needs to focus on developing and applying analytical and modeling techniques to better estimate and evaluate OH concentrations in wildfire plumes. For instance, recently, Peng et al.74 calculated HOX (OH + HO2) production rates in wildfire plumes sampled during WE-CAN from the photolysis of nitrous acid (HONO), O3, and other smaller aldehydes (e.g., formaldehyde) and ozonolysis of alkenes. These HOX production rates could be used to inform OH concentrations. Similarly, OH concentrations could be constrained by applying explicit gas-phase chemical mechanisms to reproduce the time-dependent evolution of VOCs and their oxidation products in wildfire plumes.
In addition to the uncertainties alluded to in this work, there are several additional aspects to consider while modeling the OA evolution in wildfire plumes. First, the model parameterizations in this work (e.g., POA volatility, SOA parameters) were based on simpler model systems studied in laboratory environments and these parameters may not accurately represent the processes in real wildfire plumes. Most obviously, differences in the fuel complex, burn conditions, combustion efficiency, and environmental conditions (e.g., temperature, relative humidity) and regimes (e.g., photolytic rates, NOX) could produce differences in the emissions, chemistry, and properties of OA and OA precursors between the laboratory and the field.34,68 A more specific example is that the parameters that we used to model the SOA formation from SVOCs, heterocyclics, and oxygenated aromatics (precursor classes that contributed the most to SOA production) came from environmental chamber experiments performed on a handful of surrogate species (phenol, guaiacol, and syringol, 2-methylfuran + dimethylfuran, and naphthalene, respectively) under relatively dry (relative humidity < 20%) and high NOX conditions (200–800 ppbv).86,94,102 While this extrapolation is typical for how laboratory data are translated into parameters for use in atmospheric models, these laboratory versus field differences need to be considered when evaluating model predictions against measurements.
Second, the physicochemical evolution modeled prior to the first transect remains extremely uncertain as there are no observations to evaluate those model predictions. Aircraft campaigns in the future should aim to characterize the near-field evolution in the hour after emission by performing transects closer to the fire when conditions allow. Moreover, campaigns should also accommodate repeated sampling of the near-field to assess changes in emissions over the same timescales used to perform the transect set. Any emissions changes would then need to be considered in interpreting the plume evolution inferred from the transect dataset.
Third, the current version of the SOM-TOMAS model does not simulate the photolysis or aqueous chemistry of OA or OA precursors. Photolysis has been shown to be an important loss pathway for SOA formed from monoterpenes.103,104 Oxygenated aromatics that include phenols, methoxyphenols, and phenolic carbonyls, after uptake into aerosol water, can participate in aqueous reactions to form low-volatility and light absorbing SOA.26,105 Both of these chemical processes are likely occurring in wildfire plumes and, hence, need to be included in future modeling efforts.
Fourth, the model initialization and evaluation in this work only relied on a subset of measurements made during WE-CAN. Future work could certainly benefit from leveraging an extended set of measurements gathered during WE-CAN and similar field campaigns focused on studying biomass burning emissions (e.g., BBOP,106 LASIC,107 FIREX-AQ (https://csl.noaa.gov/projects/firex-aq/)). For example, model predictions could be compared against measurements of the evolving composition (e.g., oligomers), size distribution, and thermodynamic (e.g., volatility), optical (e.g., scattering, extinction), and climate (e.g., cloud condensation nuclei) properties.
Fifth, Peng et al.74 and Hodshire et al.72 were recently able to study the distinct evolution of trace species resolved over the width of the wildfire plume. Both found evidence for increased photochemical activity near the edges and wings of the plume since these regions diluted much faster and were less optically dense compared to the core of the plume. Modeling in the future could use the information inherent in gradients within the transect to constrain the OA evolution under varying dilution and photochemical conditions.
And finally, the modeling in this work focused on simulating the plume evolution in a subset of large, daytime fires in the western US. In the future, the model will need to be applied to study a diversity of fires in terms of size, fuels, and geography (e.g., agricultural fires in the southeast US) to assess the broader applicability of the findings presented in this work.
Footnotes |
† Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d1ea00082a |
‡ Now at the Department of Civil and Environmental Engineering, University of California Davis, Davis, CA, USA. |
This journal is © The Royal Society of Chemistry 2022 |