Biomass-burning-derived particles from a wide variety of fuels - Part 1: Properties of primary particles

. Relationships between various optical, physical, and chemical properties of biomass-combustion-derived particles are characterized for particles produced in the laboratory from a wide range of fuels and burn conditions. The modiﬁed combustion efﬁciency (MCE), commonly used to parameterize biomass particle emissions and properties, is shown to generally have weak predictive capabilities, especially for more efﬁcient combustion conditions. There is, however, a strong relationship between many intensive optical properties (e.g., single-scatter albedo, Ångström absorption exponent, mass absorption efﬁciency) and the organic aerosol-to-black carbon ([OA] / [BC]) mass ratio over a wider range than previously considered (0.3 to 10 5 ). The properties of brown carbon (BrC, i.e., light-absorbing organic carbon) also vary with [OA] / [BC]. Coating-induced enhancements (i.e., “lensing”

Abstract. Relationships between various optical, physical, and chemical properties of biomass-combustion-derived particles are characterized for particles produced in the laboratory from a wide range of fuels and burn conditions. The modified combustion efficiency (MCE), commonly used to parameterize biomass particle emissions and properties, is shown to generally have weak predictive capabilities, especially for more efficient combustion conditions. There is, however, a strong relationship between many intensive optical properties (e.g., single-scatter albedo, Ångström absorption exponent, mass absorption efficiency) and the organic aerosol-to-black carbon ([OA] / [BC]) mass ratio over a wider range than previously considered (0.3 to 10 5 ). The properties of brown carbon (BrC, i.e., light-absorbing organic carbon) also vary with [OA] / [BC]. Coating-induced enhancements (i.e., "lensing" effects) contribute only a minor amount to BC absorption for all of the burns despite some burns producing particles having large ensemble-average coating-to-core mass ratios. The BC-OA mixing state varies strongly with [

Introduction
While it is understood that both open and controlled biomass combustion are major sources of particles to the atmosphere (Andreae and Merlet, 2001), questions remain regarding the properties of the emitted particles, their relationship with combustion conditions and fuel type, and their atmospheric evolution. Particles emitted from biomass combustion impact the global radiation budget and contribute to poor air quality in impacted regions. The emitted primary particles are primarily composed of organic aerosol (OA) and black carbon (BC), in varying amounts, with trace inorganic species (Reid et al., 2005;McMeeking et al., 2009;Levin et al., 2010). Particle intensive properties are often compared against the modified combustion efficiency (MCE ∼ [CO 2 ]/( [CO]+ [CO 2 ])), which provides a measure of the combustion efficiency of a burn. For example, various particle properties show some relationship with MCE, but often these relationships are weak, especially for more efficient combustion (higher MCE, corresponding typically to flaming conditions) (McMeeking et al., 2009Liu et al., 2013). Understanding the diversity in the chemical, physical, and optical properties of the emitted particles is important for establish-Published by Copernicus Publications on behalf of the European Geosciences Union. ing the fire-or region-specific emissions and subsequent impacts.
The emitted OA from biomass combustion is somewhat light absorbing (Kirchstetter et al., 2004). Absorbing OA is commonly referred to as brown carbon (BrC), with properties that appear to depend on the fuel and combustion conditions (Saleh et al., 2014;Laskin et al., 2018), which affect particle organic composition (Jen et al., 2019). However, the properties of primary BrC absorption and, especially, understanding of the relationships between BrC absorption and other particle properties and burn conditions are only beginning to be unraveled. Additionally, it is established from theory and laboratory experiments that nonabsorbing coatings on black carbon and other strongly absorbing particles can enhance the absorption (commonly referred to as the "lensing" effect but more accurately termed here the coating-induced enhancement) (Fuller et al., 1999;Lack et al., 2009;Shiraiwa et al., 2010;Cappa et al., 2012). Yet, the extent to which coating-induced enhancements impact absorption by ambient particles or for mixed-component particles from complex sources, such as biomass burning, remains contentious (Cappa et al., 2012;Healy et al., 2015;Liu et al., 2015Liu et al., , 2017Peng et al., 2016).
Here, we expand on current understanding of the relationships between various primary particle properties and burn conditions by analyzing measurements of primary biomass burning particles produced from combustion of a variety of fuel types, many of particular relevance to the western United States. We demonstrate that various optical properties exhibit a strong relationship with the [OA] / [BC] mass ratio, much stronger than their relationship with the MCE. We use the measurements to quantify the individual contributions of BC and BrC, as well as from internal mixing of BC to the observed light absorption, and we examine the variability in the properties of BrC specifically. We uniquely characterize the mixing state of BC and OA, as well as how mixing states vary between individual burns and depend on the mean properties of the emitted particles. We characterize the variability of OA-specific properties, including OA volatility, bulk chemical composition (characterized by the O : C and H : C atomic ratio, as well as the presence of key marker ions), and, uniquely, the relative abundance of organic nitrate species. We also examine the variability in the emitted particle size distribution, both for the total particles and for the BC particles specifically. Some of our analysis serves to support and extend previously determined relationships by considering a wider range of conditions, while other aspects are unique to this study. These observations provide a foundation for understanding and interpretation of experiments on the influence of photochemical aging on biomass particle properties, discussed in a related paper .

Methods
All experiments were conducted during the Fire Influence on Regional to Global Environments Experiment (FIREX) lab study, which took place at the Missoula Fire Sciences Laboratory in Missoula, MT, USA, during November, 2016. Numerous types of biomass were combusted in a large chamber (12 m × 12 m × 19 m) and the smoke was sampled to provide information on the physical, chemical, and optical properties of the resulting smoke (i.e., particulate and gas emissions). The general fuels types combusted included (exclusively or in combination) duff, dung, excelsior, straw, litter, untreated lumber, rotten debris, woody debris, shrub, herbaceous, and canopy biomass. A complete list of fuels and types is provided in Table S1 in the Supplement, with further details available on the U.S. National Oceanic and Atmospheric Administration (NOAA) data archive (https://esrl.noaa.gov/csd/ projects/firex/, last access: 24 January 2020). All data used in this publication are also available on the NOAA archive, with the processed data summarized in the complementary data repository (Cappa et al., 2019a).
Both "room" and "stack" burns were conducted, although here we include results only from stack burns. During stack burns, the smoke was mixed with background room air and funneled up a large cylindrical stack (2 m diameter ×15 m height) where it was sampled into a high-flow transfer line at ca. 0.27 m 3 s −1 . This flow rate corresponded to sampling approximately 10 % of the stack flow. Smoke was transferred to an adjacent room via the high-flow transfer line (residence time ca. 2 s) where it was subsampled through a PM 2.5 cyclone and injected into a 0.25 m 3 Teflon photochemical reaction chamber (the mini chamber). Details on the construction and operation of the mini chamber can be found in Lim et al. (2019). Here, we focus exclusively on the properties of particles sampled prior to initiation of photochemical oxidation; results of the photochemical oxidation experiments are discussed in a series of papers Lim et al., 2019). In brief, prior to each burn, the chamber was flushed with clean air with a relative humidity (RH) of approximately 40 %. To fill the chamber, smoke was subsampled from the high-flow inlet and injected across the entire burn (typically lasting for 10-20 min) or until the chamber concentration reached a maximum. A suite of instruments sampled from the mini chamber at a flow rate of approximately 6 L min −1 . This flow rate varied from burn to burn due to the exact suite of instrument sampling. Clean makeup air was being injected simultaneously from a zero-air generator to equal the air being sampled out of the chamber. The sampled smoke was diluted by a factor of ca. 7 relative to the air in the high-flow inlet. Subsequent dilution after filling was characterized by the decay of acetonitrile (ACN). Properties of the primary particles are averaged over the 5-10 min period after filling but before the initiation of photochemistry.
Particle-phase instrumentation sampled alternatingly every 2 min through a thermodenuded or ambient sample line.
The thermodenuder was operated at 300 • C with a residence time of approximately 5 s and volatilized semivolatile components, including those that are internally mixed with BC. The ambient line was lined with a charcoal cloth that removed excess gases (such as VOCs, NO x , and O 3 ) that could interfere with particle-phase measurements. Comparison of thermodenuded versus ambient particles allowed for the investigation of coating amount and volatility. The gasphase composition in the mini chamber was similar to that sampled directly from the fire (Koss et al., 2018;. Particle-phase instrumentation included a multiwavelength cavity ring-down photoacoustic absorption spectrometer (CRD-PAS) and a photoacoustic absorption spectrometer (PASS-3) for characterization of light absorption and extinction coefficients at 405, 532, and 781 nm; a highresolution aerosol mass spectrometer (HR-ToF-AMS) for characterization of nonrefractory submicron particulate matter (NR-PM 1 ) components (i.e., OA, NO 3 , SO 4 , NH 4 , Cl, K); a soot photometer AMS (SP-AMS) in laser-only mode for characterization of refractory BC and the NR components that are internally mixed with BC; a single-particle soot photometer (SP2) for characterization of refractory BC mass concentrations and size distributions; and a scanning electrical mobility spectrometer (SEMS) for measurement of particle mobility size distributions. Further details regarding instrument operation and calibration are provided in the Supplement and in Lim et al. (2019).

Bulk optical property relationships
Due to the wide variety of biomass fuels and types used during FIREX, there was a substantial diversity in the properties of primary particles produced. Previous studies have shown both the single-scatter albedo (SSA) and wavelength dependence of absorption (the absorption Ångström exponent, AAE) depend on the modified combustion efficiency (MCE) (Liu et al., 2013;McMeeking et al., 2014;Pokhrel et al., 2017). The MCE is defined here as .
The SSA is defined as where b ext is the wavelength-specific extinction coefficient and b abs is the wavelength-specific absorption coefficient. The AAE is defined as where λ 1 and λ 2 indicate two different wavelengths, here 405 and 532 nm. The MCE characterizes the overall combustion efficiency, with values closer to unity indicating more complete combustion. In general, higher MCEs correspond to more flaming combustion conditions while smaller MCEs correspond to more smoldering conditions. We find a similar relationship between SSA 405 nm , AAE, and [OA] / [BC] with MCE to previous studies ( Fig. 1) (McMeeking et al., 2009Liu et al., 2013;Pokhrel et al., 2017). Specifically, the SSA 405 nm is relatively constant and near unity for MCE <∼ 0.9, but above this value it exhibits a rapid decline, albeit with a substantial amount of scatter (Fig. 1a). The AAE is also relatively constant when MCE < 0.9, with very large values (AAE ∼ 8). There is a rapid, scattered decrease in the AAE as MCE increases further (Fig. 1b). The relationship between [OA] / [BC] and MCE is similar, with values generally decreasing as MCE increases but with a large amount of scatter (Fig. 1d). There is also a general relationship between the mass absorption coefficient referenced to BC (MAC BC ) at 405 nm and the MCE, but with similar scatter to the other properties (Fig. 1c). The MAC BC is defined as The MAC BC,405 nm includes contributions from absorption by BC and BrC, as well as from coating-induced enhancement of BC absorption. These results, along with the literature, indicate that MCE can provide guidance as to the general magnitude of these particle properties but that the MCE is ultimately a fairly imprecise metric, especially for the SSA 405 nm . However, we find a very strong relationship between the ∼ 33 (the largest value reported), which is consistent with our findings.
These observations demonstrate that the optical properties of the primary particles depend on the relative amount of OA versus BC. This is as expected because OA is generally more scattering, compared to BC, and light-absorbing OA  (Table S2). (also known as BrC) typically exhibits a much stronger wavelength dependence than BC. Based on these relationships, we divide the individual burns into different classes (Table 1). We have chosen to classify particles based on the observed SSA 405 nm values; use of [OA] / [BC] for classification yields largely similar results, given the strong relationship between the two. The dividing lines between classes are selected to yield six classes that span the entire range of SSA 405 nm values, from 0.23 (Class 1) to 0.97 (Class 6), with approximately equal numbers of individual burns in each class (ca. 8-10). Partitioning the observations into different particle classes facilitates the interpretation of the photochemical evolution of the particles, to be discussed in future work. In addition, we find that use of the class average properties versus MCE generally provides more representative fits to the observations (visually apparent in Fig. 1 and supported by the reduced χ 2 for the fits).

OA composition and volatility
Variability in the bulk composition of the OA is characterized by the O : C and H : C atomic ratios and the fractional abundance (f x ) of two marker ions, m/z = 44 and m/z = 60. The f 44 is complementary to O : C and larger values generally indicate a greater degree of oxygenation and the presence of carboxylic acids. The f 60 is often taken as a marker ion for biomass burning, in particular a signature of levoglucosan and similar molecules (Schneider et al., 2006;Alfarra et al., 2007). The high-resolution ion C 2 H 4 O + 2 contributes to and exhibits a similar behavior to f 60 ; the slope for f C 2 H 4 O + 2 against f 60 is 0.98. While it is known that properties such as f 60 vary in different biomass burning samples (Schneider et al., 2006) or between near-source intercepts of different ambient plumes (Garofalo et al., 2019), the specific dependence on burn conditions or overall particle composition (e.g., [OA] / [BC]) has not been systematically explored to our knowledge.
The average f 60 = 0.022 ± 0.01 (1σ ). The f 60 values vary nonmonotonically with [OA] / [BC], exhibiting a slight increase from Class 1 to Class 3 and then a decrease from Class 4 to Class 6 ( Fig. 2a). This indicates that, while f 60 is overall a useful marker ion for biomass burning, it cannot be used to distinguish between different burn conditions. The f 44 generally decreases with [OA] / [BC] ( Fig. 2b; r 2 = 0.33). However, the average f 44 values for particle classes 2-5 differ negligibly, suggesting that f 44 might be useful in discriminating between extreme cases (e.g., Class 1 versus Class 6) but that it is of limited general use in distinguishing between burn conditions and fuel types. The O : C atomic ratio (average = 0.37 ± 0.09) exhibits similar behavior -expected as = 0.47). This indicates that a larger proportion of ON species and functionalities are produced when particles are, on average, more BCrich. This does not reflect differences in fuel nitrogen content as there is no relationship between fuel N and f ON−OA (Fig. S2b). Therefore, it seems that the relationship between f ON−OA and [OA] / [BC] is related more so to the burn conditions than the fuel N content, although as with many other properties the relationship with [OA] / [BC] is clearer than with the MCE (Fig. S2c).
The OA volatility is characterized as the ratio between the OA concentration after thermodenuding to that without thermodenuding (the mass fraction remaining, MFR OA ). The MFR OA decreases as [OA] / [BC] increases (Fig. 2f), indicating that the OA at lower [OA] / [BC] values is less volatile than the OA at higher values. This observation provides support for the proposal by Saleh et al. (2014) that less volatile, more absorbing species are preferentially formed under conditions where BC formation is favored, which is discussed further in Sect. 3.4.2. The relationship between MFR OA and [OA] / [BC] is reasonably described by an exponential function.

BC mixing state
As discussed above, the relative amounts of OA and BC vary greatly between fuel types and combustion conditions. However, the distribution of BC and OA between particles, and how this varies between very different burn conditions, has not been previously explored in detail to our knowledge. The bulk average fraction of OA that is internally mixed with BC versus OA that is externally mixed from BC is determined using the HR-ToF-AMS and SP-AMS measurements. The HR-ToF-AMS quantifies OA independent of mixing state, whereas the SP-AMS (as operated here) quantifies only the OA that is internally mixed with BC. The fraction of OA that is internally mixed with BC (f OA,int ) is where the subscript int indicates the OA that is internally mixed with BC and the subscript tot indicates the total OA. , also with a sigmoidal relationship (Fig. 3b). Thus, while a smaller fraction of the total OA is internally mixed with BC for larger total [OA] / [BC], the amount of OA that coats BC increases. Most likely, this behavior reflects the fact that BC and OA are generated with different efficiencies in different parts of the combusting biomass. BC is more efficiently generated from flaming combustion, while OA is more efficiently generated from smoldering combustion. These observations demonstrate that the extent to which atmospheric models can assume that all OA is internally mixed with or externally mixed from BC at the point of emission will depend on the combustion conditions.

Coating-induced absorption enhancement
Non-or weakly absorbing coatings on black carbon particles can theoretically increase the absorption by BC (Fuller et al., 1999;, an effect which has been confirmed by laboratory experiments (Lack et al., 2009;Shiraiwa et al., 2010;Cappa et al., 2012). The extent to which coat- ings on BC actually enhance absorption by BC in the atmosphere remains unclear. Some studies indicate minor coatinginduced enhancements while others indicate substantial enhancements (Cappa et al., 2012(Cappa et al., , 2019bHealy et al., 2015;Liu et al., 2015Liu et al., , 2017Peng et al., 2016;Zhang et al., 2016). Understanding the nature of the coating-induced enhancement is important for quantifying the radiative impacts of BC (Jacobson, 2001;Bond et al., 2013). Further, these coatinginduced absorption enhancements (E abs.coat ) complicate the determination of brown carbon (BrC) absorption, and the two must be separated. Here, we examine the extent to which coatings on BC for primary biomass burning particles enhance the BC absorption. Theoretically, the magnitude of E abs,coat for an individual particle depends primarily on the coating thickness and secondarily on the size of the BC core Fuller et al., 1999). Thus, the extent to which coatings enhance BC absorption for a given situation can be assessed through the relationship between the observed MAC BC and the coating-to-core mass ratio (R coat-rBC = [NR-PM] int / [BC], where int indicates that the coating material is internally mixed with BC). The expectation is that the MAC BC increases with R coat-BC . However, absorption by BrC can also lead to an apparent increase in the normalized absorption with R BC if the BrC abundance correlates with the total coating amount. Because BrC absorbs more strongly at shorter wavelengths, the wavelength dependence of the MAC BC -to-R BC relationship can be used to further separate the influence of coating versus BrC absorption. The MAC BC exhibits a wavelength-dependent relationship with R coat-rBC for fresh biomass particles (405, 532, and 781 nm) (Fig. 4a-c). The MAC BC increases notably with R coat-rBC at 405 nm and to a lesser extent at 532 nm. At 781 nm the MAC BC is essentially independent of R coat-rBC up to R coat-rBC values as large as 10 but does exhibit some increase at R coat-rBC > 10. However, this is most likely a result of absorption by OA at 781 nm and not indicative of an increase in the coating-induced enhancement, discussed further below. The wavelength dependence provides clear evidence of BrC absorption at shorter wavelengths.
That the MAC BC at 781 nm is nearly independent of R coat-rBC up to such large R coat-rBC values indicates that there is only a minor coating-induced enhancement for the primary biomass particles, the magnitude of which is discussed below. Our observations are consistent with , who also investigated the relationship between the MAC BC and R coat-rBC for primary biomass particles from multiple fuel types. Most likely, this lack of a substantial coating-induced enhancement results from a noneven distribution of non-BC mass across the population of BC particles (Fierce et al., 2016;Liu et al., 2017) and from the morphology of BC-containing particles not conforming to an idealized core-shell structure (Adachi et al., 2010). The influence of photochemical aging on the coating-induced enhancement will be examined in future work.
The relationship between MAC BC and the coating amount (R coat-rBC ) can be contrasted with the relationship between MAC BC and the total [OA] / [BC]. At all three wavelengths the MAC BC 's exhibit strong, sigmoidal relationships with [OA] / [BC] (Fig. 4d-f). That MAC BC,781 nm exhibits such a clear relationship with [OA] / [BC] suggests that even the small, apparent coating-induced enhancement, implied above from the very weak with R coat-rBC , is largely driven by absorption by BrC rather than from the impact of coating on BC. Pokhrel et al. (2017) (Fig. 3b), extrapolation against [OA] / [BC] to zero effectively removes both contributions from BrC and any coating-induced enhancement. The derived MAC BC,pure values are 11.8 m 2 g 1 at 405 nm, 8.8 m 2 g −1 at 532 nm, and 5.5 m 2 g −1 at 781 nm, with estimated fit-based uncertainties of ∼ 10 %. The absolute uncertainties on the MAC BC,pure are primarily dependent on the uncertainty in the b abs and [rBC] measurements and are ∼ 35 %. The derived MAC BC,pure values are very similar to those recently reported by Forestieri et al. (2018) for fresh BC particles: MAC BC,pure = 11.9 m 2 g −1 at 405 nm and 8.8 m 2 g −1 at 532 nm, with an extrapolated value at 781 nm of 5.7 m 2 g −1 . The value at 532 nm is somewhat higher than that suggested by  (7.75 m 2 g −1 at 532 nm). Our derived MAC BC,pure values yield an AAE = 1.17, determined from a fit to the three wavelengths. An AAE close to unity indicates absorption is dominated by BC, as expected.
Values for the absorption enhancement at 781 nm are calculated as the ratio between the observed MAC BC in Fig. 4 and the derived MAC BC,pure . The derived E abs values range from 0.95 to 31. Values greater than 2 occur only for the particles having particularly large [OA] / [BC], > 400. As E abs values much greater than 2 at 781 nm are unlikely to result from mixing-induced enhancements (Chakrabarty and Heinson, 2018), this again suggests that the OA is somewhat absorbing at this wavelength. For the burns where [OA] / [BC] < 400, the median E abs = 1.14 and the arithmetic mean E abs = 1.19 ± 0.14 (1σ ). Given that some of this enhancement may result from BrC absorption at 781, these values can be considered upper limits on E abs,coat , and the small magnitude is consistent with our conclusion above that, while likely greater than 0, the mixing-induced enhancement is generally minor. It is possible that the E abs,coat values when [OA] / [BC] > 400 are substantially larger. However, given the general lack of a dependence of the MAC BC,781 nm for R BC-coat < 10, this seems unlikely.

Primary brown carbon absorption
The absorption due to brown carbon is determined by the difference as b abs,BrC = b abs,obs − MAC BC,pure · [BC] · E abs,coat , where b abs,BrC is the absorption due to BrC specifically. Importantly, the use of study-specific MAC BC,pure values serves to reduce systematic biases in the b abs,BrC , compared to direct use of literature MAC BC,pure values. Assuming E abs,coat = 1 provides an upper limit on the BrC absorption, which we note is likely most appropriate for the particles sampled here, as discussed in the previous section. Therefore, we use the upper-limit values throughout the analysis that follows, unless otherwise stated. However, a lower limit for BrC absorption can be determined at 405 and 532 nm assuming that all of the enhancement at 781 nm results from coatings and not from BrC. The resulting E abs,obs (= MAC BC,obs /MAC BC,pure ) at 781 nm averages 1.19 for R BC-coat < 10. Using E abs,coat = 1.19 in Eq. (7) yields a lower limit for the BrC absorption at the two shorter wavelengths, which is appropriate since E abs,coat generally has only a small wavelength dependence. A fit to the coatingcorrected (lower-limit) versus upper-limit b abs,BrC yields a slope of 0.97 at 405 nm and 0.88 at 532 nm (Fig. S3). The smaller difference at 405 nm results from the fractional contribution of BrC to the total absorption being larger at this wavelength. Brown-carbon-specific mass absorption coefficients (MAC BrC ) are determined as the ratio between b abs,BrC and the total OA concentration: The MAC BrC values from Eq. (7) Fig. 5a). For example, the average MAC 405 nm = 2.3 ± 1 m 2 g −1 for Class 1 and 0.35 ± 0.09 m 2 g −1 for Class 6. Although the uncertainties on the derived MAC BrC increase substantially as [OA] / [BC] decreases -because BrC absorption contributes to a smaller extent at longer wavelengths -the observations nonetheless indicate that the BrC absorptivity depends on the combustion conditions. The relationship at 405 nm is well-described by a sigmoidal function in log-log space, with limiting values of 0.35 m 2 g −1 at large [ Saleh et al. (2018), who argue that there is a continuum of BrC properties that depends on the combustion conditions, as demonstrated in that study for low-temperature benzene and toluene combustion. The range of the MAC BrC values observed here, including that there is notable absorption at 781 nm, encompass many previous measurements, summarized in Table S3. This likely reflects the wide diversity of fuel types and burn conditions considered here, as exemplified by the very large range of [OA] / [BC].
Estimated values of the imaginary component of the refractive index (RI) for BrC (k BrC ) are determined from Mie theory via optical closure , assuming a real part of the refractive index of 1.5 and a particle diameter of 150 nm, a typical value for these experiments. Imaginary RI values are of use in atmospheric models for calculation of BrC absorption. There is a linear relationship between MAC BrC and k BrC (Fig. S4a). Thus, the k BrC exhibits a similar correlation with [OA] / [BC] as does the MAC BrC (Fig. 5a).
The wavelength dependence of absorption, i.e., the AAE 405−532 , also varies with [OA] / [BC], in this case with a positive relationship between the two (Fig. 5b). The relationship is reasonably described by a sigmoidal function. This implies that, while the MAC BrC  The MAC BrC values also correlate with the organonitrogen fraction of OA, the latter of which, as noted above, also correlates with the [OA] / [BC] (Fig. 6a). This observation suggests that organic nitrate and nitro functionalities may be at least somewhat responsible for the increase in absorption. Laskin et al. (2018) performed offline molecular level analyses of primary OA collected during FIREX. They found that nitroaromatics and N-containing polycyclic aromatic hydrocarbons (PAHs) contribute notably to the total light absorption by BrC, although there are many non-N-containing species that also contribute to BrC absorption. The variability between particle classes is consistent with the results of Lin et al. (2016), which show that the abundance of N-containing chromophores varies between particles produced from different biomass fuels. Additionally, Mohr et al. (2013) observed a relationship between the concentration of nitrated phenols and short-wavelength absorption by BrC, although it is possible that for their measurements these species were produced from chemical processing, as opposed to being directly emitted. Altogether, our results provide support for the idea that organonitrogen functionalities are an important contributor to BrC absorption. However, it is very likely that other functional groups also contribute to the total absorption.
The MAC BrC,405 nm exhibits an inverse correlation with the f 60 /f 44 ratio of the OA, although there is substantial scatter in the f 60 /f 44 ratio for a given particle class (Fig. 6b). (The f 44 and f 60 have no discernable relationship.) The observed MAC BrC,405 nm relationship with f 60 /f 44 is opposite that reported by Lack et al. (2013) for ambient measurements of particles of a biomass burning plume, as they found a reasonable positive correlation. This difference in behavior results from our sampling primary particles directly -thereby focusing on the inherent variability in the properties of the emitted particles -while Lack et al. (2013) sampled ambient particles. For ambient sampling, the observed relationship will be sensitive to mixing of biomass burning particles with background or aged biomass particles, which are known to have a smaller f 60 . Thus, the relationship observed by Lack et al. (2013) can best be viewed as a mixing line between the fresh primary particles (having large MAC BrC,405 nm and large f 60 /f 44 ) and background or aged biomass particles (having small MAC BrC,405 nm and small f 60 /f 44 ), rather than providing information on the inherent variability in the absorptivity of the fresh particles.

Size distributions
Total particle mobility size distributions and BC-only size distributions were measured (Fig. 7). Primary particle size distributions are important parameters specified in regional and global models. The number-weighted and volumeweighted size distribution are generally described by either one or two lognormal modes for individual burns; a twomode fit provides a more robust solution across all modes. The mass-weighted BC size distributions are similarly described by one or two lognormal modes. A fit to the average number-weighted distribution across all particle classes yields geometric median diameters (d p,N ) and widths (σ g ) of 60.3 nm and 1.76, respectively, for the smaller mode and 153 nm and 1.64 for the larger mode (Fig. 8). The amplitude of the smaller mode is 4.6 times the larger mode. A singlemode fit yields d p,N = 68 nm and σ g = 1.93, although the fit is poorer. Mann et al. (2014) report d p,N values used by a variety of global models for biofuels. The models tend to use either 80 nm or 150 nm, although a few use other values (30, 60, 100 nm). Those using 80 nm typically use σ g = 1.80 while those using 150 nm typically use σ g = 1.59, although there are exceptions. Our observations indicate that use of Figure 7. Class-specific total particle number-weighted (top) and volume-weighted (middle) mobility size distributions, and the BC-only mass-weighted (bottom) size distribution. Individual burns are shown in gray and class averages are shown as colors. Bimodal lognormal fits are thin black lines. Note that the number-weighted and volume-weighted distributions are graphed versus mobility diameter and the BC mass-weighted distribution against the BC volume equivalent diameter. a bimodal distribution within models would be more representative, but that a single mode can perform acceptably. We find that the volume-weighted distribution calculated from a single-mode fit to the number-weighted distribution is similar to the observed volume-weighted distribution (Fig. 8). Thus, the use of a single mode to represent biomass burning size distributions thus appears acceptable, so long as the appropriate parameters are used. In this context, the widths of the distribution used by the various global models appear somewhat too small. However, we note that the microphysics occurring in the fresh smoke sampled here, which will gov-ern the size distributions, may differ from that in atmospheric plumes.
The average BC-specific mass-weighted size distribution mode is at 148 nm (Fig. 8). A bimodal fit yields values for the mass median diameter (d p,M ) and σ g of 137.2 nm and 1.62, respectively, for the smaller mode and 197.1 nm and 1.24 for the larger mode, with most of the mass contained in the smaller mode. May et al. (2014) report d p,M from laboratory biomass combustion ranging from 140 to 190 nm, averaging 170 nm. Their average is somewhat larger than ours, likely reflecting differences in the exact fuels sampled. The mode Figure 8. Average total particle number-weighted (red) and volume-weighted (blue) size distributions and the BC-specific mass-weighted size distributions. Black dashed lines are bimodal lognormal fits. The dashed blue line is the total particle volumeweighted distribution calculated from a single-mode fit to the number-weighted distribution. Note that the number-weighted and volume-weighted distributions are graphed versus mobility diameter and the BC mass-weighted distribution against the BC volume equivalent diameter. diameter for the BC-specific distribution is especially smaller than observed for biomass burning particles from some ambient observations, which tend to give values closer to 200 nm (Schwarz et al., 2008;Kondo et al., 2011;Sahu et al., 2012;May et al., 2014;Cappa et al., 2019b). This difference between lab and field observations was also noted by May et al. (2014). We speculate that the influence of coagulation may be suppressed in our experiments relative to what occurs in the atmosphere due to slower overall dilution, leading to smaller BC size distributions. To the extent this is the reason for the difference, the total particle distributions would also be biased towards particles that are too small compared to the atmosphere. However, there is no relationship between d p,N and the total particle number concentration for our experiments. Formation of secondary aerosol in the near-field of a sampled ambient plume could also contribute to this difference.
There is substantial variability between individual burns within a given particle class in terms of the shape of the size distributions (Fig. 7). This variability is most evident for Class 1, 2, and 5 but present for all classes. Nonetheless, the number-weighted mean diameter (d p,N,mean ) appears to decrease somewhat with MCE ( Fig. 9; r 2 = 0.38). However, the relationship is largely driven by the Class 6 particles, which generally have lower MCE values, having larger d p,N,mean values. A lack of any particularly clear relationship is consistent with Hosseini et al. (2010), who observed the d p,N,mean to exhibit a complex relationship with combustion conditions. The d p,N,mean varies nonmonotonically  (Fig. 9). This is despite the notable burn-to-burn variability. It is important to note that the mobility-based size is particle shape dependent; very-BC-rich particles are more likely to have nonspherical shapes and thus have larger mobility diameters. This could explain the minimum in d p,N around Class 3 particles, for which [OA] / [BC] = 10.
Some of the variability within a class appears related to the presence of different fuel types within a class. Numberweighted and BC-specific mass-weighted size distributions by fuel type are shown in Fig. 10. For the number-weighted distributions, leaf litter and rotten logs exhibit the greatest variability between different burns, although we note that multiple burns were not performed for all fuels. The shapes of the leaf litter, peat and "other" fuel types, differ most notably from the other fuel types, with the presence of more than one mode more apparent. (The "other" category here includes nontraditional biofuels, specifically building materials and excelsior.) For the BC-specific size distributions, the litter, canopy, and duff exhibited the greatest intrafuel variability. For most fuels, the BC-specific distribution peaks around 150 nm, as noted above. However, for a subset of burns (eight of them) the BC-specific distribution peaks around 100 nm Figure 10. Normalized total particle number-weighted (top) and the BC-only mass-weighted (bottom) size distributions shown by fuel type (see legend). Individual burns are gray and averages for a fuel type in color. For some fuels there is only one size distribution. Bimodal lognormal fits are the black lines. The "other" category includes nontraditional biofuels, specifically building materials and excelsior. (Fig. 10). These small BC-mode distributions occur for the OA-rich particle classes 4, 5 and 6 ( Fig. 7), although there is no clear pattern to their occurrence.

Conclusions and implications
Measurements of primary particles produced from combustion of a variety of biomass fuel types indicate the optical, physical, and chemical properties of the emitted particles exhibit wide variability. We show that variability in many optical properties (e.g., single-scatter albedo, wavelength dependence of absorption, mass absorptivity of black and brown carbon) is directly linked to the [OA] / [BC] ratio of the emitted particles; the relationships with [OA] / [BC] are much stronger than with the commonly used modified combustion efficiency, and mathematical relationships between the various properties are determined. However, the absorption enhancement due to coating of BC (the so-called lensing effect) is shown to be minor and essentially independent of the amount of coating up to large coating-to-BC mass ratios. The brown carbon mass absorptivity correlates with the organonitrogen fraction of OA, suggesting that organonitrogen species contribute to BrC absorption. Many bulk chemical properties (i.e., O : C, H : C, and the relative concentra-tions of key marker ions such as f 60 ) exhibit limited dependence on the burn conditions and the [OA] / [BC] ratio. However, both the OA volatility and organonitrogen fraction of OA decrease with [OA] / [BC]. The fraction of OA that is internally mixed with BC was shown to decrease strongly with the [OA] / [BC] ratio, from nearly all OA being internally mixed with BC when the particles are overall BC-rich to only a few percent of OA being mixed with BC when OA dominates. Yet, the relative amount of OA coating the BC increases with [OA] / [BC]; that is, when more of the OA is externally mixed from BC those particles that do contain BC nonetheless have thicker OA coatings. The observed total particle size distributions are reasonably well described by a single lognormal mode but are better fit using a bimodal distribution. The BC-specific size distributions are similarly best fit using a bimodal distribution, although a single mode provides a reasonable representation. The dependence of the geometric median mobility diameter on the burn conditions or particle state (i.e., the [OA] / [BC]) is complicated by the mobility diameter being sensitive to variations in particle shape, which depend on the [OA] / [BC] ratio. Overall, these results expand on previous observations of primary biomass burning particle properties, considering a wider range of [OA] / [BC] and associated properties. Further, they provide a foundation for understanding the postemission evolution of biomass burning smoke due to photochemical oxidation as discussed in Lim et al. (2019).
Data availability. All data are available from the NOAA FIREX-AQ data repository (https://esrl.noaa.gov/csd/projects/firex/firelab/, last access: 24 January 2020) (National Oceanic and Atmospheric Administration Chemical Sciences Division, 2020). This includes a summary of the fuel types used for each burn and the measurement time series for each burn. The primary particle averages used in this work are additionally collected in the UC DASH data repository (Cappa et al., 2019a).
Author contributions. CDC and JHK designed the experiments. CDC, CYL, and DHH carried out the measurements and data processing. CDC, CDM, and CYL analyzed data. CDC and CDM wrote the paper, with contributions from all coauthors.
Competing interests. The authors declare that they have no conflict of interest.
Acknowledgements. The entire FIREX team, especially Bob Yokelson and Jim Roberts and the staff of the Missoula Fire Sciences Laboratory, is acknowledged for their assistance. Putting together the community inlet was a community effort -thank you to all who contributed. Shuka Schwarz and Gavin McMeeking are also thanked for their assistance with the SP2, and Tim Onasch is thanked for his assistance with the SP-AMS.
Financial support. This research has been supported by the National Oceanic and Atmospheric Administration's Climate Program Office (grant nos. NA16OAR4310111 and NA16OAR4310112). CYL was additionally supported by the National Science Foundation Graduate Research Fellowship Program.
Review statement. This paper was edited by Manabu Shiraiwa and reviewed by two anonymous referees.