Chemical composition, structures, and light absorption of N-containing aromatic compounds emitted from burning wood and charcoal in household cookstoves

N-containing aromatic compounds (NACs) are an important group of light-absorbing molecules in the atmosphere. They are often observed in combustion emissions, but their chemical formulas and structural characteristics remain uncertain. In this study, red oakwood and charcoal fuels were burned in cookstoves using the standard water-boiling test (WBT) procedure. Submicron aerosol particles in the cookstove emissions were collected using quartz (Qf) and polytetrafluoroethylene (PTFE) filter membranes positioned in parallel. A backup quartz filter (Qb) was also installed downstream of the PTFE filter to evaluate the effect of sampling artifacts on NAC measurements. Liquid chromatography– mass spectroscopy (LC–MS) techniques identified 17 NAC chemical formulas in the cookstove emissions. The average concentrations of total NACs in Qb samples (0.37± 0.31– 1.79±0.77 μgm−3) were greater than 50 % of those observed in the Qf samples (0.51±0.43–3.91±2.06 μgm−3), and the Qb-to-Qf mass ratios of individual NACs had a range of 0.02–2.71, indicating that the identified NACs might have substantial fractions remaining in the gas phase. In comparison to other sources, cookstove emissions from red oak or charcoal fuels did not exhibit unique NAC structural features but had distinct NAC composition. However, before identifying NAC sources by combining their structural and compositional information, the gas-particle partitioning behaviors of NACs should be further investigated. The average contributions of total NACs to the light absorption of organic matter at λ= 365 nm (1.10 %–2.57 %) in Qf and Qb samples (10.7 %–21.0 %) are up to 10 times larger than their mass contributions (Qf: 0.31 %–1.01 %; Qb: 1.08 %– 3.31 %), so the identified NACs are mostly strong light absorbers. To explain more sample extract absorption, future research is needed to understand the chemical and optical properties of high-molecular-weight (e.g., molecular weight, MW> 500 Da) entities in particulate matter. Published by Copernicus Publications on behalf of the European Geosciences Union. 14078 M. Xie et al.: N-containing aromatic compounds emitted from burning wood and charcoal

Household solid-fuel combustion is a leading human health risk, especially for women and children, who tend to spend more time indoors than men (Anenberg et al., 2013).Estimates show that exposures to PM 2.5 from domestic solidfuel combustion caused 3.9 million premature deaths and a loss of healthy life years of ∼ 4.8 % (Smith et al., 2014).In addition, the emissions of carbonaceous aerosols from cookstoves can affect the Earth's radiative balance by absorbing and scattering incoming solar radiation (Lacey and Henze, 2015;Aunan et al., 2009).BC is the most efficient light absorber in the atmosphere, while the total aerosol absorption, including that from OC, is still highly uncertain (Yang et al.,2009;Park et al., 2010;Feng et al., 2013;Wang et al., 2014;Tuccella et al., 2020).Multiple field and laboratory studies have demonstrated that OC in both primary PM emissions (e.g., biomass and fossil fuel combustions) and secondary organic aerosol (SOA) feature a range of absorptivity in the near-ultraviolet (UV) and short-visible-wavelength regions (Nakayama et al., 2010;Forrister et al., 2015;Lin et al., 2015;De Haan et al., 2017;Xie et al., 2017aXie et al., , b, 2018)).The light-absorbing OC fraction is often referred to as "brown carbon" (BrC).Unlike open BB (e.g., forest, grassland, and cropland fires) -one of the most important primary sources for organic aerosols (Bond et al., 2004) -the light absorption of BrC from household cookstove emissions is rarely investigated.Sun et al. (2017) found that the BrC absorption from residential coal burning accounted for 26.5 % of the total aerosol absorption at 350-850 nm.BrC from wood combustion in cookstoves has a greater mass-specific absorption than that from open BB over the wavelength range of 300-550 nm (Xie et al., 2018).These results suggest that cookstove emissions may also be an important BrC source, which needs to be accounted for separately from open BB.
Organic molecular markers (OMMs) are commonly used in receptor-based source apportionment of carbonaceous aerosols (Jaeckels et al., 2007;Shrivastava et al., 2007;Xie et al., 2012).Polycyclic aromatic hydrocarbons (PAHs) and their derivatives are a group of OMMs with light absorption properties dependent on ring number or the degree of conjugation (Samburova et al., 2016).As discussed in Xie et al. (2019), PAHs are generated from a multitude of combustion processes (e.g., BB, fossil fuel combustion; Chen et al., 2005;Riddle et al., 2007;Samburova et al., 2016), and their ubiquitous nature makes them less than ideal OMMs for BrC source attribution.Because of the specific toxicological concern raised by PAHs -they are mutagenic and carcinogenic (International Agency for Research on Cancer, 2010) -source emission factors (EFs), ambient levels, and potential health effects of PAHs are investigated exhaustively (Ravindra et al., 2008;Kim et al., 2013).Similar to PAHs, Ncontaining aromatic compounds (NACs) are a group of BrC chromophores commonly detected in ambient PM and source emissions.Zhang et al. (2013) and Teich et al. (2017) calculated the absorption of individual NACs in aqueous extracts of ambient PM, the total of which explained ∼ 3 % of the bulk extract absorption at 365-370 nm.With the same approach, Xie et al. (2017aXie et al. ( , 2019) ) found that the absorbance due to NACs in BB or secondary OC was 3-10 times higher than their mass contributions.Lin et al. (2016Lin et al. ( , 2017) ) estimated an absorbance contribution of 50 %-80 % from NACs in BB OC directly from their high-performance liquid chromatographer (HPLC)-photodiode array (PDA) signals, which are subject to considerable uncertainty due to the coelution of other BrC chromophores (e.g., PAHs and their derivatives).These results indicate that NACs are strong BrC chromophores, but the estimation of their contributions to BrC absorption depends largely on how well they are chemically characterized.Nitrophenols, methyl nitrophenols, nitrocatechols, and methyl nitrocatechols (including isomers) are typical atmospheric NACs (Claeys et al., 2012;Desyaterik et al., 2013;Zhang et al., 2013).These NACs can be generated from BB (Lin et al., 2016(Lin et al., , 2017;;Xie et al., 2019), fossil fuel combustion (Lu et al., 2019), and the reactions of aromatic volatile organic compounds (VOCs) with reactive nitrogen species (e.g., NO x ; Xie et al., 2017a) and are not unique to specific sources (e.g., BB).By using an HPLC interfaced to a diode array detector (DAD) and quadrupole (Q) time-of-flight mass spectrometer (ToF-MS), Xie et al. (2019) found that BB NACs contain methoxy and cyanate groups.Nitronaphthol, nitrobenzenetriol, and methyl nitrobenzenetriol are characteristic NACs for NO x -based chamber reactions of naphthalene, benzene, and m-cresol, respectively (Xie et al., 2017a).Yet, few studies have investigated the composition of NACs from household cookstove emissions (Fleming et al., 2018;Lu et al., 2019).
The present study aims to characterize NACs in PM 2.5 from burning red oak and charcoal in a variety of cookstoves and calculate their contributions to bulk OC absorption.The absorption of OC in solvent extracts of cookstove emissions was measured in our previous work (Xie et al., 2018).Presently, NACs are identified and quantified using an HPLC-DAD-Q-ToF-MS system described earlier.In addition, the NACs adsorbed on a backup quartz filter downstream of a polytetrafluoroethylene (PTFE) membrane filter are analyzed to evaluate the potential for sampling artifacts of PM 2.5 NACs on the bare quartz filter in parallel.This work unveils BrC composition at a molecular level and increases the understanding of BrC chromophores and their sources.It also shows that further identification of large molecules (e.g., > 500 Da) may better explain BrC absorption in the particle phase.

Cookstove emission sampling
The cookstove emission test facility, fuel-cookstove combinations, water-boiling test (WBT) protocol, and PM 2.5 emission sampling were described previously in Jetter and Kariher (2009) and Jetter et al. (2012).Briefly, the cookstove emission tests were performed at the United States Environmental Protection Agency (US EPA) cookstove test facility in Research Triangle Park, NC, USA.Red oakwood and lump charcoal were burned in fuel-specific cookstoves under controlled conditions.Emission tests for each fuelcookstove combination were performed in triplicate.The WBT protocol (version 4; Global Alliance for Clean Cookstoves, 2014) is designed to measure cookstove power, energy efficiency, and fuel use and contains cold-start (CS) high-power, hot-start (HS) high-power, and simmer (SIM) low-power phases.Both CS and HS phases are defined by the duration between the ignition and the water boils.The CS phase starts with the cookstove, pot, and water at ambient temperature; the HS immediately follows the CS with the cookstove hot but the pot and water at ambient temperature; and the SIM phase is defined by a 30 min time period with the cookstove hot and water temperature maintained at 3 • C below the boiling point.Low-moisture (∼ 10 %) oak and charcoal fuels were burned with five specifically designed cookstove types (Tables S1 and S2 in the Supplement); highmoisture (∼ 30 %) oak fuels were burned in one cookstove (Jiko Poa, BURN Manufacturing, Kenya).A brief description of each fuel-specific cookstove is given in the Supplement (Sect.S1).Gaseous pollutant (e.g., CO; methane, CH 4 ) emissions were monitored continuously, and PM 2.5 filter samples were collected during each test phase of the WBT protocol.The modified combustion efficiency (MCE), defined as CO 2 /(CO 2 + CO) on a molar basis, was calculated and discussed in Xie et al. (2018).A quartz-fiber filter (Q f ) and a PTFE membrane filter positioned in parallel collected PM 2.5 isokinetically at a flow rate of 16.7 L min −1 .The adsorption artifact of Q f was evaluated using a quartzfiber backup filter (Q b ) installed downstream of the PTFE filter during PM 2.5 sampling.

Chemical analysis
The OC and elemental-carbon (EC) emissions and UVvis light absorption properties (BrC) of methanol-extracted cookstove particles were reported in Xie et al. (2018).Details for determinations of OC-EC concentrations and BrC absorption are provided in the Supplement (Sect.S2).Except the three-stone fire, EFs of OC and EC in the SIM phase were substantially lower than those in high-power phases (CS and HS), so the BrC absorption from red oak and charcoal burning was primarily measured for CS-and HS-phase samples in Xie et al. (2018).The SIM-phase samples were analyzed only for red oak burning in a three-stone fire.This test had comparable OC emissions between CS-and SIM-phase combustions, and CS and HS phases of the three-stone fire were typically similar and could not be separated (Xie et al., 2018).In the current work, the same emission samples were selected for the analysis of NACs, and the three SIM-phase samples from the three-stone fire were treated as HS-phase samples of other cookstove tests.Tables S1 and S2 summarized the measurement results of Q f and Q b , respectively, for each fuelcookstove combination, including concentrations of carbon contents and light-absorbing properties of sample extracts.As the light absorption of BB BrC is expected to depend largely on burn conditions (Saleh et al., 2014;Pokhrel et al., 2016), the MCE and EC/OC ratio, two indicators of burn conditions, are also given in Table S1.
The Q f and Q b sample extraction and subsequent analysis for NACs were conducted as described in Xie et al. (2019).In brief, an aliquot of each filter sample was prespiked with 250 ng nitrophenol-d4 (internal standard) and extracted ultrasonically twice for 15 min in 3-5 mL of methanol.After filtration (30 mm diameter × 0.2 µm pore size; PTFE filter, National Scientific Co. Ltd, TN, USA), the extract volume was reduced to ∼ 500 µL with rotary evaporation prior to HPLC-DAD-MS (Q-ToF) analysis.The NACs targeted in this work were chromatographed using an Agilent 1200 Series HPLC equipped with a Zorbax Eclipse Plus C18 column (2.1 mm × 100 mm, 1.8 µm particle size; Agilent Technologies, CA, USA).The gradient separation was performed using water (eluent A) and methanol (eluent B) containing 0.2 % acetic acid (v/v) with a total flow rate of 0.2 mL min −1 .The eluent B fraction was held at 25 % for 3 min; increased to 100 % over the next 7 min, where it was held for 22 min; and then returned to 25 % over 5 min.An Agilent 6520 Q-ToF MS equipped with a multimode ion source operating in electrospray ionization (ESI) negative (−) mode was used to determine the chemical formula, molecular weight (MW), and quantity of each target compound.All sample extracts were analyzed in full-scan mode over 40-1000 Da.A mass accuracy of ±10 ppm was selected for compound identification and quantification.Samples with individual NACs exhibiting the highest MS signal intensities in full-scan mode were re-examined in targeted MS-MS mode using a collisioninduced dissociation (CID) technique.The MS-MS spectra of target NACs [M-H] − ions were acquired to deduce structural information.Similar to bulk carbon and light absorption measurements, NACs were primarily determined for CS-and HS-phase samples with substantial OC loadings.
Due to the limited availability of authentic standards, many of the NACs identified in cookstove combustion samples were quantified using surrogate compounds with similar MW or structures.An internal standard method with a ninepoint calibration curve (∼ 0.01-2 ng µL −1 ) was applied for quantification of concentrations.The compounds represented by each identified NAC formula were quantified individually and combined to calculate the mass ratio of total NACs to OC (µg m −3 ) ×100 % (tNAC OC %).Presently, the ratio of organic matter (OM) to OC was not measured or estimated for cookstove combustion emissions, so tNAC OC % could be up to 2 times greater than the contributions of NACs to OM (Reff et al., 2009;Turpin and Lim, 2001).Table S3 lists the chemical formulas, proposed structures, and standard assignments for the NACs identified here.The quality assurance and control (QA-QC) procedures for filter extraction and instrumental analysis were the same as Xie et al. (2017aXie et al. ( , 2019)).NACs were not detected in field blank and background samples.The average recoveries of NAC standards on prespiked blank filters ranged from 75.1 % to 116 %, and the method detection limit had a range of 0.70-17.6pg.

Data analysis
In Xie et al. (2017a), the DAD measurement directly identified the chemical compounds in chamber SOA responsible for light absorption in the near-UV and visible-light ranges.However, no light absorption from individual NACs was detected in the DAD chromatograms from open BB (Xie et al., 2019) and cookstove emissions (this work).So the contributions of individual NACs to light absorption coefficient (Abs λ ; Mm −1 ) for each sample extract at 365 nm (Abs 365,iNAC %) were calculated using the method described in Xie et al. (2017aXie et al. ( , 2019)): where C iNAC is the mass concentration (ng m −3 ) of individual NACs, and MAC 365,iNAC is the mass absorption coefficient (MAC λ , m 2 g −1 ) of individual NACs at 365 nm.Abs 365 is the light absorption coefficient (Mm −1 ) of each sample extract at 365 nm and has been widely used to represent BrC absorption (Chen and Bond, 2010;Hecobian et al., 2010;Liu et al., 2013).Each NAC compound was assumed to absorb as a standard (Table S3); of these NAC compounds, the MAC 365,iNAC value was obtained from Xie et al. (2017aXie et al. ( , 2019) ) and listed in Table S4.In this work, Student's t test was used to determine if the means of two sets of data are significantly different from each other, and a p value less than 0.05 indicates significant difference.
3 Results and discussion

Summary of total NAC concentration from cookstove emissions
Table 1 summarizes the average concentrations of total NACs and average tNAC OC % for Q f and Q b by fuel type and WBT phase.The EFs of total NACs shown in Table S5 were obtained by multiplying the EFs of OC and tNAC OC %.Filter samples of emissions from burning red oakwood had significantly (p < 0.05) higher average total NAC concentrations and tNAC OC % than the charcoal-burning samples.Wood burning generates more volatile aromatic compounds (e.g., phenols, PAHs) than charcoal burning (Kim Oanh et al., 1999), and NACs can form when aromatic compounds and reactive nitrogen (e.g., NO x ) are present during solidfuel combustion (Lin et al., 2016(Lin et al., , 2017)).While burning red oak, emissions from the CS and HS phases show similar average NAC concentrations, tNAC OC %, and NAC EFs (Tables 1 and S5).Additionally, burning low-moisture red oak in the Jiko Poa stove had higher tNAC OC % than burning highmoisture red oak (Tables S6 and S7), but the difference was not significant (p > 0.05).Thus, the NAC emissions from red oak burning are less likely influenced by the WBT phase, and the effect of fuel moisture content needs further investigation.For charcoal fuel samples, compared with the CS phase, the HS phase shows significantly higher (p < 0.05) average NAC concentrations and EFs.This is likely due to the increase in OC with the HS phase (Tables 1 and S5) as the average tNAC OC % values are much closer for the CS (0.40 ± 0.25 %) and HS phases (0.31 ± 0.21 %).
Several studies have placed a quartz-fiber filter behind a PTFE filter to evaluate the positive adsorption artifact -adsorption of gas-phase compounds onto particle filter media ("blow-on" effect; Peters et al., 2000;Subramanian et al., 2004;Watson et al., 2009;Xie et al., 2014).This method is expected to provide a consistent estimate irrespective of sampling time but may overcorrect the positive artifact by 16 %-20 % due to volatilization of OC off the upstream PTFE filter (negative artifact, "blow-off" effect; Subramanian et al., 2004).A denuder upstream of the filter for gas sampling was used to avoid positive artifacts in several studies (Ding et al., 2002;Ahrens et al., 2012).This approach can generate large negative artifacts by altering the gas-particle equilibrium after the denuder, and a denuder efficiency of 100 % might not be guaranteed (Kirchstetter et al., 2001;Subramanian et al., 2004).The present study is the first to consider sampling artifacts when measuring semivolatile NACs.This concept merits consideration as quantification of particlephase NACs may be subject to large uncertainty.Table 1 shows that the average concentrations of total NACs on Q b (0.37 ± 0.31-1.79± 0.77 µg m −3 ) are greater than 50 % and 80 % for those on Q f (0.51±0.43-3.91±2.06µg m −3 ) for red oak and charcoal burning, respectively.The average Q b -to-Q f ratio in percentage using OC concentrations is 2-3 times lower (14.8±3.87 %-38.8±18.9%).Hence, the NACs identified in this work are present in the relatively volatile bulk OC fraction emitted from cookstoves, and the NACs in the Q f samples may also be present in the gas phase in the atmosphere.Charcoal-burning emissions show even higher (p < 0.05) Q b -to-Q f total NAC mass ratios (CS: 84.1 ± 38.0 %; HS: 140 ± 52.9 %) than red oak burning (CS: 50.8 ± 13.4 %; HS: 53.4 ± 26.2 %), which is largely due to the higher OC loads on Q f from red oak burning.Along with modified combustion efficiency (MCE), the EC/OC and BC/OA (organic aerosol) ratios were used previously as indicators of biomass-burning conditions (McMeeking et al., 2014;Pokhrel et al., 2016).Here the burn condition indicates general flame intensity or combustion temperature (Chen and Bond, 2010;Saleh et al., 2014) and is parameterized to investigate combustion processes (e.g., pyrolysis).The MCE and the EC/OC and BC/OA ratios are key to understanding particulate OC absorptivity (Saleh et al., 2014;Lu et al., 2015) and NAC formation from open BB (Xie et al., 2019).Presently, the relationships of tNAC OC % versus EC/OC for Q f samples are shown in Fig. S1c and f by fuel type.Because no significant difference was observed for average total NAC concentrations, tNAC OC %, and EC/OC ratios when testing CS versus HS phases during red oak fuel burning, the CS and HS phases were pooled for a regression analysis.The tNAC OC % values of Q f samples positively correlate (r = 0.83, p < 0.05) with EC/OC for red oak burning (Fig. S1c), as observed in Xie et al. (2019) for open BB, which suggests that burn conditions influence NAC formation during BB.Note that the NAC concentrations on Q f were possibly adsorbed while in a gaseous state, while EC is particle-phase.
In Table S1, the MCE values of charcoal burning indicate that the HS phase is more smoldering than the CS phase.However, the average tNAC OC % values showed no significant difference (p = 0.29) between HS and CS phases.Like MAC 365 and Å abs in Q f samples for charcoal burning (Xie et al., 2018), tNAC OC % derived from the same samples did not correlate with EC/OC ratios in this work (Fig. S1f).Xie et al. (2018) found that the HS phase for charcoal burning had average OC EFs 5-10 times higher than the CS phase, while the EC EFs decreased by more than 90 % from the CS to HS phase.Furthermore, no correlation has been observed between MCE and EC/OC for charcoal burning in the HS phase.So, the EC/OC for charcoal burning tends to depend https://doi.org/10.5194/acp-20-14077-2020Atmos.Chem.Phys., 20, 14077-14090, 2020 more on the initial temperature in the cookstove than MCE variations and cannot be used to predict burn conditions, BrC absorption, or NAC formation.

Composition of NACs in Q f and Q b
During solid-fuel combustion, NACs may form from aromatic compounds (e.g., substituted phenols) and reactive nitrogen species (e.g., NH 3 , NO x , and HONO) in both the gas and particle phases (Harrison et al., 2005;Kwamena and Abbatt, 2008;Lu et al., 2011;Lin et al., 2016Lin et al., , 2017)).Aromatic hydrocarbons are produced during fuel pyrolysis (Simoneit et al., 1993;Simoneit, 2002;Kaal et al., 2009).Oxidation of fuel-derived nitrogen rather than molecular nitrogen in air is the major formation pathway of reactive nitrogen species (Glarborg et al., 2003).Presently, 17 chemical formulas were identified as NACs in cookstove emissions, several of which are widely observed in ambient air and open BB particles (e.g., C 6 H 5 NO 3 , C 6 H 5 NO 4 ; Claeys et al., 2012;Zhang et al., 2013;Lin et al., 2016Lin et al., , 2017;;Xie et al., 2019).Figure 1 shows the average concentrations (ng m −3 ) of individual NACs in Q f and Q b samples by fuel type and WBT phase.The corresponding average mass ratios of individual NACs to OC × 100 % (iNAC OC %) are exhibited in Fig. S2.Details of the NAC composition expressed in iNAC OC % for each fuelcookstove experiment are given in Tables S6-S9.
Generally, the CS and HS phases have consistent NAC profiles for red oak combustion (Figs.1a, b and S2a, b).C 10 H 7 NO 3 (CS-Q f : 1003 ± 803 ng m −3 ; HS-Q f : 1149 ± 1053 ng m −3 ) and C 8 H 5 NO 2 (CS-Q f : 712±921 ng m −3 ; HS-Q f : 1185 ± 1761 ng m −3 ) have the highest average concentrations on Q f , followed by C 11 H 9 NO 3 , C 10 H 11 NO 5 , and C 11 H 13 NO 5 .However, C 8 H 5 NO 2 was only detected in emission samples of Jiko Poa among the five wood stoves (Tables S6 and S7).Not considering C 8 H 5 NO 2 , Q b samples of red oak combustion emissions have similar NAC profiles and characteristic species (e.g., C 10 H 7 NO 3 , C 11 H 9 NO 3 ) as Q f samples, and the Q b -to-Q f mass ratios of individual NACs are similar between the CS and HS phases (Fig. 1a, b).It appears that the formation of NACs from red oak burning in cookstoves depends largely on burn conditions reflected by EC/OC ratios (Fig. S1c) rather than WBT phases.Among the 17 identified NACs from red oak burning, C 8 H 5 NO 2 and C 11 H 13 NO 6 have the lowest Q b -to-Q f ratios (2.03 %-9.80 %; Fig. 1a, b), indicating their low volatility.The low volatility of C 11 H 13 NO 6 might be due to its relatively high MW, while C 8 H 5 NO 2 has the second-lowest MW, and its structure likely contains functional groups that decrease vapor pressure (e.g., carboxyl group; Donahue et al., 2011).
As mentioned in Sect.3.1, using a Q b has been widely applied to evaluate the positive sampling artifact for OC and semivolatile organic compounds.This method might only work for bulk PM, OC, and organic compounds with low volatility, of which the concentrations in Q b samples are much lower than Q f samples and usually presumed to be due to positive adsorption artifacts only (Subramanian et al., 2004;Watson et al., 2009).In this work, the average Q bto-Q f mass ratios of the 17 individual NACs ranged from 50.8±13.4% to 140±52.9 %, comparable to n-alkanes with carbon number ≤ 21 (e.g., henicosane: 26.3 %-163 %) and PAHs with benzene ring number ≤ 4 (e.g., fluoranthene: 46.3 %-134 %) in the ambient of urban Denver (Xie et al., 2014).Xie et al. (2014) found that the gas-phase concentrations of n-alkanes and PAHs with vapor pressure greater than henicosane and fluoranthene were comparable to or higher than their particle-phase concentrations.The vapor pressure of five NAC standards at 25 • C (p o, * L ) was predicted using the US EPA Toxicity Estimation Software Tool (TEST) and listed in Table S10.Their p o, * L values are mostly higher than henicosane and fluoranthene (∼ 10 −8 atm; Xie et al., 2013Xie et al., , 2014)).Then the identified NACs in this study may have substantial fractions remaining in the gas phase.As the evaporation of NACs from the upstream filter (negative artifact) is unknown, the particle-phase NAC concentrations cannot be calculated by simply subtracting Q b measurements from those of Q f .Considering that most of the Q f and Q b samples were collected near ambient temperature (Table S2; ∼ 25 • C), the composition of NACs derived from Q f measurements alone can be biased due to the lack of gas-phase measurements.Future work is needed to evaluate the composition of NACs from emission sources in both the particle and gas phases.comparison, the MS-MS spectra of standard compounds used in this work are obtained from Xie et al. (2017aXie et al. ( , 2019) ) and shown in Fig. S5.Among all identified NAC formulas, C 10 H 7 NO 3 was detected in each fuel-cookstove experiment (Tables S6-S9) and showed the highest concentrations in emissions from burning red oak (Fig. 1a,  b).The MS-MS spectrum of C 10 H 7 NO 3 (Fig. S4l) is like 2-nitro-1-phenol (Fig. S5g) but shows a ∼ 1 min difference in retention time (Fig. S3i; 10.9 min, 2-nitro-1-phenol 11.8 min).C 10 H 7 NO 3 is presumed to be an isomer of 2nitro-1-phenol with a nitronaphthol structure.C 11 H 9 NO 3 has a degree of unsaturation and a fragmentation pattern (Fig. S4q) like C 10 H 7 NO 3 and is likely a structural isomer of methyl nitronaphthol.C 6 H 5 NO 3 , C 7 H 7 NO 3 , C 6 H 5 NO 4 , and C 7 H 7 NO 4 are commonly detected in combustion emissions (Lin et al., 2016(Lin et al., , 2017;;Xie et al., 2019) and atmospheric particles (Claeys et al., 2012;Zhang et al., 2013).

Identification of NAC structures
C 6 H 5 NO 3 and C 6 H 5 NO 4 are identified as 4-nitrophenol and 4-nitrocatechol using authentic standards (Figs.S4a, d and  S5a, c).C 7 H 7 NO 3 has two isomers (Fig. S3b), and the compound eluting at 9.98 min has the same retention time and MS-MS spectrum (Fig. S4c) as 2-methyl-4-nitrophenol (Fig. S5b).In ambient PM and chamber SOA, C 7 H 7 NO 4 was identified using standard compounds as a series of methylnitrocatechol isomers (4-methyl-5-nitrocatechol, 3-methyl-5-nitrocatechol, and 3-methyl-6-nitrocatechol; Iinuma et al., 2010).According to the HPLC-Q-ToF-MS data for C 7 H 7 NO 4 identified in Iinuma et al. (2010) and our previous studies (Xie et al., 2017a(Xie et al., , 2019)), the two C 7 H 7 NO 4 isomers in Fig. S3d are likely 4-methyl-5-nitrocatechol and 3-methyl-6-nitrocatechol, respectively.Here we cannot rule out the presence of 3-methyl-5-nitrocatechol, which may coelute with 4-methyl-5-nitrocatechol (Iinuma et al., 2010).In Fig. S4k, o, and p, the MS-MS spectra of https://doi.org/10.5194/acp-20-14077-2020Atmos.Chem.Phys., 20, 14077-14090, 2020 C 7 H 7 NO 5 , C 8 H 7 NO 5 , and C 8 H 9 NO 5 all show a loss of CH 3 + NO(orNO 2 ) + CO.The loss of CH 3 is typically due to a methoxy group in NAC molecules, and NO (or NO 2 ) and CO loss is commonly observed for NACs with more than one phenoxy group (Xie et al., 2019).So methoxy nitrophenol is the proposed skeleton for C 7 H 7 NO 5 , C 8 H 7 NO 5 , and C 8 H 9 NO 5 .Other functional groups were estimated using their chemical formulas and degree of unsaturation as a basis (Table S3).The present study quantifies C 8 H 7 NO 4 and C 9 H 9 NO 4 using 2-methyl-5-benzoic acid (C 8 H 7 NO 4 ) and 2,5-dimethyl-4-nitrobenzoic acid (C 9 H 9 NO 4 ), respectively.The fragmentation patterns of C 8 H 7 NO 4 (Fig. S4g, h) and C 9 H 9 NO 4 compounds (Fig. S4m, n) are different from their corresponding surrogates (Fig. S5f, h), and loss of CO 2 is not observed, so C 8 H 7 NO 4 and C 9 H 9 NO 4 compound structures do not include a carboxyl group.The MS-MS spectra of C 8 H 7 NO 4 eluting at 8.14 min (Fig. S3e) and C 9 H 9 NO 4 eluting at 9.22 min (Fig. S3j) indicate the loss of an OCN benzisoxazole structure or the presence of cyanate (−O−C≡N) or isocyanate (−O=C=N) groups.Mass spectra of selected standard compounds (Fig. S5i-n) in our previous work (Xie et al., 2019) show that the loss of an OCN group only happens during the fragmentation of phenyl cyanate.Thus, the C 8 H 7 NO 4 and C 9 H 9 NO 4 isomers containing OCN indicate a phenyl cyanate feature.However, the fragmentation mechanism related to the loss of a single nitrogen for the second C 8 H 7 NO 4 isomer (Figs.S3e and S4h) is unknown and requires further study.The MS-MS spectrum of the second C 9 H 9 NO 4 isomer had dominant ions at m/z 194 ([M-H] − ), 164 (loss of NO), and 149 (loss of NO + CH 3 ).Compared with the MS-MS spectra of 4-nitrophenol and 2-methyl-4nitrophenol (Fig. S5a, b), the second C 9 H 9 NO 4 isomer is likely a methoxy nitrophenol with an extra ethyl group.
The EIC signal of C 8 H 9 NO 4 in Fig. S3f comprises at least 3-4 isomers, and the MS-MS spectra are always dominated by ions at m/z 182 ([M-H] − ), 152 (loss of NO), and 137 (loss of NO + CH 3 ), with some changes in relative abundance.The fragmentation mechanism of C 8 H 9 NO 4 represented by the MS-MS spectrum in Fig. S4i is consistent with that of the second C 9 H 9 NO 4 isomer (Fig. S4n), so the C 8 H 9 NO 4 might also have a methoxy nitrophenol skeleton.The MS-MS spectrum of C 8 H 5 NO 2 is characterized by CO 2 loss (Fig. S4j), indicative of a carboxyl group.Considering the degree of unsaturation of the C 8 H 5 NO 2 molecule and the cyano group feature in BB tracers (e.g., hydrogen cyanide, benzonitrile; Schneider et al., 1997;Li et al., 2000;Gilman et al., 2015), C 8 H 5 NO 2 was identified as 4-cyanobenzoic acid using an authentic standard (Fig. S5o).The C 10 H 11 NO 4 , C 10 H 11 NO 5 , C 11 H 13 NO 5 , and C 11 H 13 NO 6 detected here are also observed in other BB experiments (Xie et al., 2019).Their MS-MS spectra are characterized by the loss of at least one CH 3 and/or OCN (Fig. S4r-u), suggestive of methoxy or cyanate groups.Without authentic standards, fragmentation patterns (Fig. S4r-u) were used to determine the molecu-lar structures of C 10 H 11 NO 4 , C 10 H 11 NO 5 , C 11 H 13 NO 5 , and C 11 H 13 NO 6 (Table S3).
Nearly all NAC formulas identified in this work were observed previously (Lin et al., 2016(Lin et al., , 2017;;Xie et al., 2017a;Fleming et al., 2018;Xie et al., 2019).Few studies attempt to retrieve structural information for NACs using MS-MS spectra of authentic standards.Although multiple NACs may be generated from BB and photo-oxidation of aromatics in the presence of NO x , NAC structures may differ across emission sources.Xie et al. (2019) found that fragmentation patterns of C 7 H 7 NO 5 and C 8 H 9 NO 5 from BB and photochemical reactions are distinct, and the methoxy and cyanate groups are featured only in BB NACs.Thus, knowing the NAC structure is useful to emission source identification.In this work, the chemical and structural information obtained for NACs sampled during red oak and charcoal burning is similar, presumably because the charcoal fuel used is produced by the slow pyrolysis of wood.However, NACs in red-oak-and charcoalburning emissions can be differentiated compositionally.As shown in Figs. 1 and S2, the NAC emissions from red oak burning in cookstoves are characterized by C 10 H 7 NO 3 and C 11 H 9 NO 3 .In addition to these two species, charcoal burning in cookstoves also generates high fractions of C 8 H 9 NO 5 (Fig. S2c, d).This difference among NACs may help with source apportionment using receptor models, which are commonly used and assume that the ambient pollutants measured in the field are linear combinations from a number of timevariant sources or factors.(Jaeckels et al., 2007;Shrivastava et al., 2007;Xie et al., 2013).
Figure 2 compares NAC composition from cookstove emissions (not including C 8 H 5 NO 2 ), open BB (Xie et al., 2019), and SOA chamber experiments (Xie et al., 2017a).Since previous source emission studies ignored Q b measurements and normalized individual NAC contributions to OM, only Q f measurements in this work are compared (Fig. 2a,  b) with their iNAC OC % values multiplied by 1.7 (proposed OM/OC ratio; Reff et al., 2009).The three open BB tests (Fig. 2c) were conducted with two fuel types under different ambient temperatures (10-29 • C) and relative humidity (49 %-83 %; Xie et al., 2019).But they consistently emit C 6 H 5 NO 4 , C 7 H 7 NO 4 , and C 9 H 9 NO 4 , which is compositionally distinct from cookstove emissions (Fig. 2a, b).Moreover, the average mass contribution of total NACs to OM for open BB (0.12 ± 0.051 %) was 4-14 times lower than that for cookstove emissions.This result is likely due to the hightemperature flaming combustion produced in the cookstoves (Shen et al., 2012;Xie et al., 2018).In Fig. 2d and e, the NAC profiles for photochemical reactions appear to have aromatic precursors.When using field measurement data of NACs for receptor modeling, the resulting factors can be linked with specific emission sources by comparing with the NAC patterns shown in Fig. 2. Further studies are also warranted to unveil NAC patterns of other potential sources (e.g., motor vehicle emissions).Therefore, the source of NACs can be identified by combining their characteristic structures and composition.The filter-based NACs reported for the experiments shown in Fig. 2 were all measured using an identical method and HPLC-Q-ToF-MS instrument, reducing any potential methodological bias.However, total gas-phase NAC concentrations need to be properly sampled and measured to account for the impact of gas-to-particle partitioning on their distribution.

Contributions of NACs to Abs 365
The average Abs 365,iNAC % values of Q f and Q b samples are presented by fuel type and WBT phase in the Fig. 3 stack plots, and experimental data for each fuel-cookstove experiment are provided in Tables S11-S14.The average contributions of total NACs to Abs 365 (Abs 365,tNAC %) of the sample extracts (Q f : 1.10 %-2.57%; Q b : 10.7 %-21.0 %) are up to 10 times greater than their average tNAC OC % (Q f : 0.31 %-1.01 %; Q b : 1.08 %-3.31 %; Table 1).Considering that some NACs are not light-absorbing (Table S4), and the OM/OC ratio is typically greater than unity, most NACs that contribute to Abs 365 are strong BrC chromophores.Like the mass composition of NACs (Fig. 1), C 10 H 7 NO 3 (CS: 0.24 %; HS: 0.43 %) and C 8 H 9 NO 5 (CS: 1.22 %; HS: 0.55 %) were the major contributors to Abs 365 for the Q f samples collected during red oak and charcoal burning, respectively (Fig. 3a).The average Abs 365,tNAC % values of Q b samples are 7.53 to 11.3 times higher than those of Q f samples.Unlike the Q f samples from red oak burning, C 10 H 11 NO 5 (CS: 2.77 %; HS: 3.09 %) has the highest average contribution to Abs 365 for Q b samples, followed by C 10 H 7 NO 3 (CS: 1.96 %; HS: 1.32 %) and C 8 H 9 NO 5 (CS: 1.32 %; HS: 1.44 %).C 8 H 9 NO 5 dominated the contribution (CS: 8.78 %; HS: 5.82 %) to https://doi.org/10.5194/acp-20-14077-2020Atmos.Chem.Phys., 20, 14077-14090, 2020 Abs 365 for the Q b samples from charcoal burning (Fig. 3b).All identified NACs explained 1.10 %-2.58 % (Fig. S3) of Q f extract absorption.Even if the NACs on Q b were totally derived from upstream filter evaporation, the adjusted average contributions of total NACs (Q f + Q b ) to Abs 365 of Q f extracts were still lower than 5 % (1.59 %-4.01 %).Due to the lack of authentic standards, the quantification of NAC concentrations and their contributions to Abs 365 of Q f extracts might be subject to uncertainties.However, growing evidence has shown that BrC absorption was majorly contributed by large molecules with an MW > 500-1000 Da (Di Lorenzo and Young, 2016;Di Lorenzo et al., 2017).Large molecules of NACs may be generated from flaming combustions in cookstoves, and their structures and light absorption are worth future investigations.In previous studies on ambient and biomass-burning particles, most identified NACs had an MW lower than 300-500 Da, and their total contributions to bulk BrC absorption were estimated to be less than 10 % (Mohr et al., 2013;Zhang et al., 2013;Teich et al., 2017;Xie et al., 2019).Similar results were also obtained in the current work.Therefore, further studies are needed to identify large BrC molecules (including high-MW NACs) in ambient and source particles.

Conclusion
This study investigated the composition, chemical formulas, and structures of NACs in PM 2.5 emitted from burning red oak and charcoal in a variety of cookstoves.Total NAC mass and compositional differences between Q f and Q b samples suggest that the identified NACs might have substantial gasphase concentrations.By comparing the MS-MS spectra of identified NACs to standard compound spectra, the structures of NACs featuring methoxy and cyanate groups in cookstove emissions are confirmed.The source identification of NACs would be less ambiguous if both the structures and composition of NACs are known as different emission sources have distinct NAC characteristics.However, the compositional information of NACs based on Q f measurements only is biased due to the lack of gas-phase data, and further studies are warranted to investigate the gas-to-particle distribution of NACs in the ambient and source emissions.Similar to previous work, the average contribution of total NACs to Abs 365 of Q f samples is less than 5 % (1.10 %-2.57%), suggesting the need to shift our focus from low-MW NACs (MW < 300 Da) to the chemical and optical properties of large molecules (e.g., MW > 500 Da) in particles.
Xie et al. (2018) assumed previously that the Q b -adsorbed OC represented the positive sampling artifact only and adjusted the light-absorbing properties of OC on Q f by subtracting Abs 365 and OC of Q b samples directly.In this study, the high Q b -to-Q f ratios of total NACs indicate that the volatilization of NACs from the upstream PTFE filter cannot be neglected, but the relative contributions of positive and negative artifacts to Q b measurements are unknown.Therefore, the measurement results of NACs in Q f and Q b samples were provided separately, and no correction was conducted for Q f measurements in this work.Since the gaseous NACs adsorbed in Q b samples depends on Q f loadings, tNAC OC % and total NAC concentrations in each Q f -Q b pair from matching tests are significantly correlated (p < 0.05; Fig.S1a, b, d, and e in the Supplement).

Figures
Figures S3 and S4exhibited extracted ion chromatograms (EICs) and MS-MS spectra of the 17 identified NACs.For

Figure 1 .
Figure 1.Average concentrations of individual NACs in Q f and Q b samples for (a) red oak burning in the CS phase, (b) red oak burning in the HS phase, (c) charcoal burning in the CS phase, and (d) charcoal burning in the HS phase.The blue scatters in each plot are Q b -to-Q f mass ratios of individual NACs ×100 %.

Figure 2 .
Figure 2. Average mass ratios (%) of individual NACs to organic matter from (a) red oak burning in cookstoves; (b) charcoal burning in cookstoves; (c) open BB experiments(Xie et al., 2019); and photochemical reactions of (d) toluene and ethylbenzene and (e) benzene, naphthalene, and m-cresol with NO x(Xie et al., 2017a).

Figure 3 .
Figure 3. Average contributions (%) of individual NACs to Abs 365 of bulk (a) Q f and (b) Q b extracts from burning red oak and charcoal in cookstoves in CS and HS phases.

Table 1 .
Average concentrations of total NACs and tNAC OC % in Q f and Q b samples by fuel type and WBT phase.