Comprehensive organic emission profiles, secondary organic aerosol production potential, and OH reactivity of domestic fuel combustion in Delhi, India

104 | Environ. Sci.: Atmos., 2021, 1, 10 rganic emission profiles, secondary organic aerosol production potential, and OH reactivity of domestic fuel combustion in Delhi, India† Gareth J. Stewart, a Beth S. Nelson, a W. Joe F. Acton,‡ Adam R. Vaughan, James R. Hopkins, Siti S. M. Yunus, d C. Nicholas Hewitt, b Eiko Nemitz, e Tuhin K. Mandal, Ranu Gadi, Lokesh. K. Sahu, Andrew R. Rickard, ac James D. Lee ac and Jacqueline F. Hamilton*


Introduction
Around 3 billion people globally use solid fuels to meet their daily cooking energy requirements. 1 Emissions from residential solid fuel combustion are signicant and have been shown to cause indoor air pollution which resulted in 2.8-3.9 million premature deaths globally, 1-3 with around 25% of ambient particulate matter (PM) in South Asia related to cooking emissions. 4 Approximately a quarter of worldwide residential solid fuel use is in India, 5 where cooking domestically over biomass remains popular because biomass fuel is cheaper than liqueed petroleum gas (LPG) and meals cooked with traditional methods perceived to be tastier. 6 Recent studies have shown that 16% 7 of non-methane hydrocarbons and 27% 8 of non-methane volatile organic compounds (NMVOCs) by mixing ratio at different urban sites in Delhi were from solid fuel combustion sources. Furthermore, Aerosol Mass Spectrometer measurements found that crop residue burning and solid fuel combustion jointly accounted for 24% (35.8 mg m À3 ) of the concentration of PM with a diameter <1 mm (PM 1 ) during the post-monsoon in Delhi, with likely additional contributions to the SOA. 9 Studies focussed on organic emissions from both open biomass burning and domestic solid fuel combustion have shown that organic components are released over a range of volatilities. [10][11][12][13][14] These include non-methane volatile organic compounds (NMVOCs, effective saturation concentration, C*, 3 Â 10 6 to 10 11 mg m À3 ), intermediate-volatility organic compounds (IVOCs, C* ¼ 300 to 3 Â 10 6 mg m À3 ), semi-volatile organic compounds (SVOCs, C* ¼ 0.3-300 mg m À3 ) as well as low-and extremely low-volatility organic compounds (L/ELVOCs, where LVOC C* # 0.3 mg m À3 ). 15 As a result, I/SVOCs from domestic solid fuel combustion potentially represent a large global source of SOA, however, the effect of I/SVOCs on OH reactivity, aging and SOA formation remains poorly understood. [16][17][18] The factors controlling SOA formation are complex. These include the oxidation of NMVOCs to less volatile products which partition into the particle phase, the heterogeneous oxidation of particle-phase SVOCs, and plume dilution with subsequent SVOC evaporation followed by further gas-phase oxidation. 19 Of 17 studies examining the enhancement factor of organic aerosol (OA) to CO from the aging of open biomass burning emissions, 10 found no increase in SOA, 4 found an increase and 3 reported a decrease. 20 Despite varied results, a recent lab study has shown SOA formation from combustion of fuels relevant to open biomass burning to be signicant.  showed a carbon yield of SOA from NMVOCs emitted from the combustion of western U.S. fuels of 24 AE 4% when exposed to atmospheric aging equivalent to 6 hours, which increased to 56 AE 9% aer aging equivalent to 6 days. 19 Formation of SOA from open biomass burning has been examined as part of several recent studies. Hatch et al. (2015) estimated that 8-15% of SOA from the combustion of black spruce, cut grass, Indonesian peat and ponderosa pine was because of furanic compounds. The contribution of furanic compounds to SOA was estimated to be greater still (28-50%) from rice straw and wiregrass. 11 Gilman et al. (2015) examined the relative contributors to SOA from the combustion of U.S. fuels and found the main contributors to be polyunsaturated oxygenated NMVOCs. 21 High SOA formation potential was driven by benzene diols, benzaldehyde, and phenols. Ahern et al. (2019) showed that for the combustion of coniferous fuels, which were dominated by the burning of biomass needles, biogenic NMVOCs were the most important class of SOA precursor. 22 Akherati et al. (2020) reported that oxygenated aromatic compounds resulted in just under 60% of the SOA from the combustion of western U.S. fuels. 23 These studies have also shown that reactive chemical species such as furanics, oxygenated aromatics and aliphatics are important drivers of the OH reactivity of open biomass burning emissions. 10,21,24 Recent model simulations by  focussed on modelling the OH radical chemistry in emissions from the combustion of fuels from the western U.S. showed that up to 10% of O 3 in the rst 4 h aer emission was a result of the oxidation of furanic compounds. 25 Few studies have examined SOA formation from fuels used for domestic solid fuel combustion, with little known about the impact of the species released on the reactivity of emissions. Bruns et al. (2016) examined SOA formation from the combustion of beech fuel wood and demonstrated that the main contributors were 22 compounds, and in some cases up to 80% of the SOA produced was estimated to be formed from phenol, naphthalene and benzene. 26 A further study suggested that furanic and phenolic compounds were important precursors to SOA as a result of spruce combustion. 24 The concentration of primary organic aerosol (POA) is determined by dynamic gas-to-particle partitioning of an extremely complex mixture of organics over a wide range of volatilities. Understanding the gas-to-particle partitioning represents one of the main difficulties in accurately characterising SOA formation, as measurements of organic emissions using multiple measurement techniques are required. As a result, gas-phase emissions are traditionally considered up to C 12 (saturation vapour concentration, C*, $10 6 mg m À3 ) and POA as non-volatile. [27][28][29][30] Consequently, many models neglect the importance of I/SVOCs as SOA precursors. The effect is a signicant underestimation of SOA production and an overestimation of POA in chemical transport models. 20,31,32 The concentration of organic aerosol (OA) is determined by the volatility of species and ambient conditions, with many source tests occurring at unrealistically high OA concentrations. Laboratory-based source studies typically enhance the POA emission factor relative to more dilute ambient conditions. 28,33 The inclusion of I/SVOCs leads to better agreement between modelled and measured values. 32,[34][35][36] A range of studies have been conducted to comprehensively characterise organic emissions from mobile sources 30,37-39 and aircra engines, 40,41 however, a need has been highlighted to develop source proles for both open and domestic biomass burning. 30 These have the potential to result in a better understanding of the SOA formed from the I/SVOCs released. Comprehensive source proles are far better suited to predicting SOA formation than traditional separated gas-and particle-phase emission factors developed at the point of emission. Comprehensive proles can be adjusted to real-world dilutions, aerosol concentrations and temperatures. These parameters are all likely to have a large inuence on the mass of SOA present.
This study develops comprehensive, model-ready organic emission proles for solid fuels routinely burnt in the Delhi area of India. These proles account for the full range of volatilities of organic emissions to better constrain the impact of domestic biomass burning on SOA formation. This study also compares the relative impacts of different solid fuel combustion sources to SOA production potential and OH reactivity and examines the most important chemical contributors.

Datasets
This paper combines previously published data 13,14 of organic emissions from a detailed eld campaign designed to measure organic emissions from the combustion of solid fuels widely used in Delhi, India. The data has been used to create comprehensive organic emission proles and examine the largest contributors to secondary organic aerosol production potential and the OH reactivity of emissions.
Fuels were collected from across Delhi in a manner designed to reect the range and variability of solid fuels used across the region. Sample collection included a range of fuel woods (Melia azedarach, Prosopis spp, Eucalyptus spp, Azadirachta indica, Mangifera indica, Morus spp, Pithecellobium spp, Shorea spp, Ficus religosa, Syzgium spp, Ficus spp, Vachellia spp, Dalbergai sissoo, Ricinus spp, Holopetlea spp, Saraca indica and plywood), cow dung cake, municipal solid waste (MSW, collected from 3 landll sites: Ghazipur, Bhalswa and Okhla), crop residues (Brassica spp, Solanum melongena and Cocos nucifera) and individual samples of charcoal, sawdust and LPG. Samples were stored in a manner akin to local storage practices prior to combustion. Combustion experiments were conducted using expert local judgement to ensure that the experimental setup replicated real-world burning conditions, using a combustiondilution chamber at the CSIR-National Physical Laboratory (NPL), New Delhi, that has been well described previously (see the ESI1 † for a schematic of the combustion chamber used). [42][43][44][45] Fuel (200 g) was rapidly heated to spontaneous ignition using an electric heater, with emissions convectively driven into a hood and up a ue to allow enough dilution, cooling and residence time to achieve the quenching typical of indoor environments. These conditions have been previously optimised so that the chamber setup did not alter combustion conditions. 42 The methodology was designed to replicate the convection-driven conditions of real-world combustion and was adapted from the VITA water-boiling test. Temperature measurements were made directly above the fuel being combusted (see the ESI2 † for combustion temperatures by fuel type) and ambient temperatures during experiments ranged from 25-38 C. All instruments made measurements over a 30 minute period, to ensure that emissions were collected from both highand low-temperature combustion conditions. All measurements used the same procedures to characterise emissions to create a self-consistent dataset of speciated organic emissions spanning a large range of volatilities.
NMVOCs were sampled from the top of the ue down a 1 4 00 PFA sample line, which was subsampled by three separate online gas-phase instruments designed to target a wide range of NMVOCs of different functionality and volatility. A dualchannel gas chromatograph with ame ionisation detection (DC-GC-FID) was used to sample alkanes from ethane to nhexane and a range of small alkenes from ethene-isoprene (see the ESI3 † for DC-GC-FID method). A two-dimensional gas chromatograph with ame ionisation detection (GCÂGC-FID) was used to sample alkanes from n-heptane to n-dodecane, aromatic species from benzene to monoaromatics with up to 5 carbon substituents and up to 12 monoterpenes (see the ESI4 † for GCÂGC-FID method). GC instruments were calibrated using 4 ppbv gas standards containing a range of alkanes, alkenes, and aromatics purchased from the British National Physical Laboratory. Blank measurements were made at the beginning, middle and end of the day, with mean values subtracted from measured emission factors for both GC instruments.
A proton-transfer-reaction time-of-ight mass spectrometer (PTR-ToF-MS, PTR 8000; Ionicon Analytik, Innsbruck) was used to sample a range of small oxygenates, oxygenated aromatics, alkenes, furanic species and nitrogen-containing volatile organic compounds (see the ESI5 † for PTR-ToF-MS method). Calibrations were performed twice a week with a gas calibration unit (Ionicon Analytik, Innsbruck) using a calibration gas that contained 18 compounds (Apel-Riemer Environmental Inc., Miami). Before each burn, ambient air was sampled to provide a background for the measurement. Measurements were made by the PTR-ToF-MS at 1 second, which were averaged to the 30 minute sample window of the GC instruments to create the selfconsistent dataset of gas-phase organic emissions in this study.

Mapping organics to volatility basis data set
The volatility-basis dataset (VBS) is designed to simulate the emission and evolution of I/SVOCs into the atmosphere and places NMVOCs into logarithmically spaced bins of C* at 298 K. 47 Emissions from fuel wood, cow dung cake, MSW and LPG were mapped onto a VBS to visualise and compare emissions across the entire range of volatilities measured using data collected by the DC-GC-FID, GCÂGC-FID, PTR-ToF-MS, SPE-GCÂGC-ToF-MS (>C 12 ) and PTFE-GCÂGC-ToF-MS. C* values were calculated for individual NMVOCs measured online using the DC-GC-FID, GCÂGC-FID and PTR-ToF-MS instruments. For SPE-GCÂGC-ToF-MS and PTFE-GCÂGC-ToF-MS analyses, organics were lumped into groups of unspeciated compounds. These were spaced between n-alkanes, with the volatility assigned as the mean volatility of the alkanes either side of the bin. For NMVOCs where insufficient data was available for a calculation of C*, the volatility was assigned as the C* of the nalkane with the nearest boiling point. Despite calibrating C* of organic material to n-alkanes potentially leading to an overestimation of the volatilities of polycyclic aromatic hydrocarbons (PAHs) and aromatic oxygenates, 48 this unspeciated material represented a small fraction of the total organic matter. The C * i for each NMVOC i was calculated using (1): The constant 760 was used to convert between units of atm and Torr where 1 atm ¼ 760 Torr. P o L,i values were taken from the EPA Estimation Programme Interface Suite data. 49 Grouped regions of organics from SPE disks and PTFE lters were calibrated to allow semi-quantication based on the mean total ion current (TIC) chromatogram of the two n-alkanes either side of the bin close to a concentration of $1 mg mL À1 . The approach was uncertain and suggestions for better quantication of this complex organic material are provided elsewhere. 13 Experimental and/or predicted vapour pressures of species, especially the n-alkanes used for assigning volatility bins, remained uncertain. We adopted a similar approach to Lu et al. (2018), with the factor of 10 spacing of volatility bins to minimise the chance of volatility misassignment. 30 IVOCs were in the n-alkane range $ C 12 to C 22 , SVOCs from C 23 to C 32 , and L/ ELVOCs from C 33 to C 40 . Care was taken to avoid double counting of species measured using multiple techniques (see the ESI7 † for the hierarchy of instruments used to measure different species). Gas-phase species, which were possible to measure using either of the GC instruments or the PTR-ToF-MS, were counted once only. In summary C 2 -C 6 non-methane volatile organic compounds (alkanes/alkenes) were measured using the DC-GC-FID. C 7 -C 12 non-methane volatile organic compounds (alkanes and benzene to C 3 substituted monoaromatics) were measured using the GCÂGC-FID. Remaining NMVOCs and gas-phase I/SVOCs were measured using the PTR-ToF-MS (C 4 -C 5 substituted monoaromatics, phenolics, furanics, oxygenated aromatics, oxygenated aliphatics and nitrogen containing volatile organic compounds). The uniden-tied gaseous I/SVOC fraction was estimated using SPE-GCÂGC-ToF-MS. The organic aerosol fraction was measured using the PTFE-GCÂGC-ToF-MS. To allow incorporation of I/ SVOCs species from SPE disks, species and their isomers measured using the PTR-ToF-MS were removed from the SPE and PTFE analyses.

Comparison of EPA and fuel wood source proles
Fuel wood source proles were compared to those from the EPA SPECIATE 5.0 (2019) database. Notably proles from the EPA for burning sources were split into either gas-or particle-phase measurements. All available proles for residential combustion were considered from sets of experiments including replace wood combustion (4640-4642), residential combustion using wood and pellet stoves (95 129-95 138), residential wood stove combustion (95 156-95 159) and residential wood combustion (G95467-G95470). EPA proles 95 156-95 159 for residential wood stove combustion were not directly compared due to the low number of organic species measured (n ¼ 37). Comparison was made to Pinus ponderosa (G95467), Eucalyptus spp (4640) and a wood stove (95133). This placed into context the VBS developed in this work, because the multiple and complementary techniques used here allowed simultaneous measurement of organics in both gas and aerosol phases.
EPA G95467 was a source prole derived from measurement of 179 organic species from combustion of Pinus ponderosa. 50 C 2 -C 12 compounds were collected into canisters and analysed by GC-FID/-MS, C 8 -C 20 compounds were collected onto Tenax tubes and analysed by GC-FID/-MS, carbonyls were collected onto 2,4-dinitrophenylhydrazine cartridges and analysed by HPLC and ne particles and SVOCs were collected onto lter/ PUF/XAD/PUF cartridges and analysed by GC-MS. EPA 4640 was a source prole derived from measurement of 85 organic species from Eucalyptus spp. This prole was chosen as Eucalyptus spp was also measured as part of this study. Gas phase semi-volatile species were collected onto PUF cartridges, particles collected onto lters and carbonyls onto C 18 cartridges impregnated with dinitrophenylhydrazine. Samples were then extracted and analysed by GC-MS. 51 EPA 95133 was developed by sampling VOCs into Tedlar gas sampling bags followed by GC-FID analysis and semi-volatile PAHs were collected onto PUF plugs, extracted and analysed by GC-MS. 52

Estimation of the SOA formation potential
The overall yield of SOA, g SOA , from gas-phase emissions from domestic biomass burning samples (mass of SOA produced/ mass of NMOG emissions) was estimated using: 30 where f gas,i ¼ the mass fraction of SOA precursor as a proportion of total mass of gas-phase emissions and g i ¼ yield of SOA precursor i at a concentration of OA ¼ 10 mg m À3 . SOA yields were calculated from literature (see the ESI8 † with OA mass loadings as close to 10 mg m À3 under both high and low NO x conditions, where supporting information from relevant literature was available), with gas-phase SVOCs assumed to have SOA mass yields of 1. 30 The rate of reaction of chemical species with OH was not included in (2). IVOCs usually react faster than NMVOCs with OH, and so IVOCs and NMVOCs contribute differently to SOA with respect to time. 39 As a result, the approach here estimates a lower-bound contribution of the ultimate yield of IVOCs to SOA. 30 It also does not include species which may form SOA heterogeneously, because the traditional SOA yield values have not been assigned, yet.

Estimation of OH reactivity
The OH reactivity of emissions from different fuel types was examined to understand the largest contributors. The mean concentrations of NMVOCs from the DC-GC-FID, GCÂGC-FID and PTR-ToF-MS were used to calculate OH reactivity, s À1 , using: where: ½M ¼ pressure ðmbarÞ Â 10 À4 ð8:314 Â ð273:15 þ temperatureÞÞ Â 6:023 Â 10 23 (4) and the rate constants for reaction with OH, k OH , used in this study are given in the ESI9. †

Estimation of PAH toxicity
Toxicity equivalence factors (TEFs) were used to assess the relative toxicity of PAH emissions per kg of fuel burnt. TEFs indicate the relative toxicity of a PAH to benzo[a]pyrene (BaP), one of the most carcinogenic PAHs. 53 The toxicity of a PAH was expressed in BaP equivalents ([BaP] eq ), which was calculated in (5) by multiplying the concentration of PAH i , in nanograms per cubic metre (ng m À3 ), by the corresponding TEF for i, TEF i , given in the ESI10. † 54-56  Fig. 1 displays a comprehensive characterisation of organic emissions and emphasised the importance of using multiple techniques to measure organic emissions, covering a volatility range of over 13 orders of magnitude. Cow dung cake combustion released more SVOCs and L/ELVOCs than fuel wood. Fig. 1 illustrates the particle fraction, X p , which was calculated according to the method in Lu et al. (2018), assuming all the organic emissions formed a quasi-ideal solution when diluted to ambient conditions. 30 The particle fraction demonstrated the gas-to-particle partitioning of organics at typical atmospheric conditions (T ¼ 298 K and OA concentration ¼ 10 mg m À3 ). IVOCs were predominantly found in the gas phase and SVOCs were present in both phases. The predicted particle fraction suggested that there should have been more gas-phase contributions in the I/SVOC range. The amount of organic material in the gas and particle phase is dependent on multiple factors such as temperature and concentration of OA. It is likely that at the high concentrations (OA > 10 mg m À3 ) during source testing, a larger fraction of I/SVOCs partitioned into the particle phase.  cow dung cake, (C) MSW, (D) crop residue, (E) charcoal, (F) sawdust and (G) LPG. The largest mass fraction at emission for all sources, expect LPG, was in the range C* $ 10 7 -10 9 mg m À3 and a result of small oxygenated species. LPG emission was dominated by fugitive emissions of propane and butane from C* 10 9 -10 10 mg m À3 . Fig. 2 highlights how changes in the type of source inuenced emissions of I/S/L/ELVOCs. All sources, except LPG, had signicant emissions of IVOCs. Fig. 2A also shows comparison to EPA source proles G95467 for Pinus ponderosa (red circles), 95 133 for a wood stove (green triangles) and 4640 for Eucalyptus spp (blue squares). These proles highlighted the difficulties in using current source proles to predict SOA from domestic biomass burning, due to signicantly different predictions in the range C* $ 10 2 -10 6 mg m À3 .

Chemical composition distribution
EPA 95133 reported essentially no IVOCs, EPA G95467 showed some IVOCs in the range C* $ 10 5 -10 6 mg m À3 and EPA 4640 showed considerably higher IVOC emissions. EPA 95133 did not measure important I/SVOC species released from domestic biomass burning such as phenolics and furanics and therefore no organic matter was represented for C* < 5 Â 10 6 mg m À3 . EPA G95467 was one of the best current source proles, however, no organic matter was present in this prole for C* < 10 5 mg m À3 . This may be due to lack of simultaneous gas-and particle-phase measurements of all organic species present. As a result, gas-phase organic species may have partitioned into the particle phase because of high organic aerosol concentrations during source testing and are therefore were not represented. EPA 4640 measured more phenolic and furanic compounds. Despite this, the measurement of only 85 organic species in EPA 4640 overemphasised the importance of I/SVOCs as a mass fraction. This therefore still posed signicant problems when using EPA 4640 to model SOA formation. These issues demonstrated the benet of the VBS developed here, as simultaneous measurement of organics in both gas and aerosol phases should alleviate these problems. For some sources, such as the combustion of MSW, cow dung cake, crop residues and sawdust a greater mass fraction of I/SVOCs was released. The use of a VBS for these sources is likely even more important due the presence of large amounts of I/SVOC material.
The comparison results suggested that the prole presented for Indian domestic fuel wood from this study was signicantly lower than the data in EPA G95467 and 95 133 in the two most volatile bins (C* $ 10 10 -10 11 mg m À3 ). However, this was an artefact due to the lack of measurements with the DC-GC-FID, which targeted the most volatile species. However, the underestimation of emissions from C* $ 10 10 -10 11 mg m À3 is unlikely to be signicant when calculating SOA formation using the measured species and VBS presented. Table 1 shows the mass fraction of organic material presented in Fig. 2 from the 7 different sources studied here, presented in volatility bins spanning over 13 orders of magnitude. For certain sources, such as LPG and charcoal, only one sample was taken. The lack of repeat measurements signicantly increased the uncertainty associated with the VBS presented. Despite this, multiple gas-phase NMVOC measurements were made. These showed similar results and therefore these VBS were included. The results in Table 1 should be used to better characterise SOA formation in chemical-transport models from domestic biomass combustion sources as the volatility distribution of organic emissions presented can be accurately adjusted to atmospheric dilutions, aerosol concentrations and temperatures. Fig. 3 shows that the mean mass fractions of IVOCs emitted increased from MSW (0.12 AE 0.02) to fuel wood (0.15 AE 0.04) to crop residue (0.16 AE 0.04) to cow dung cake (0.18 AE 0.02). SVOC emissions for fuel wood and crop residue were the lowest mass  [57][58][59][60] It is likely that the higher emissions of OA from cow dung cake, MSW and sawdust were a result of the lower combustion efficiency of these samples. Fig. 4 shows the sum of the mass fraction of NMVOCs released from domestic fuel burning in this study. Only the species identied as SOA precursors and assigned with SOA yields in the ESI8 † were included. The mass fraction of SOA precursors from fuel wood, crop residue and cow dung cake samples were from 0.3-0.5. Compared to sources calculated using the same method, 30 the mass fraction which resulted in SOA was less, with the exception of MSW burning, compared to gasoline ($0.65) and diesel ($0.7) engines. This was principally due to the large emission of smaller oxygenated species from burning samples. MSW burning samples released the largest mass fraction of SOA precursors (0.4-0.65). Fig. 5A and B show the estimated SOA yields from burning samples under high and low NO x conditions. These were intended to represent idealised systems for photo-oxidation of SOA precursors. Under high NO x conditions, RO 2 radicals react with NO x and under low NO x conditions RO 2 radicals react with HO 2 . 61 Consideration was given to both cases since domestic biomass burning in India impacts both urban high NO x regions and rural lower NO x regions. Under high NO x conditions, SOA yields were lower and IVOCs represented a larger proportion of the total SOA produced. Under low NO x conditions, SOA yields were greater, and NMVOCs resulted in a greater proportion of the total SOA due to higher estimated SOA yields from aromatic and furanic species. Other studies examining emissions from burning have traditionally considered yields from only one of these regimes, but greater SOA production under low NO x conditions has been well described previously. 62,63 Fig . 5C shows that high NO x SOA yields from sawdust, charcoal, cow dung cake, fuel wood and crop residue were likely  dominated by phenolics (light blue, 21-70%), with a signicant contribution from furanics (orange, 9-33%) due to high emission factors of these species and high SOA yields. Other important SOA contributions were from aromatics (2-8%), oxygenated aromatics (2-8%), oxygenated aliphatic species (2-9%), monoterpenes (0-7%) and PAHs (2-16%). A larger proportion of SOA (40%) from MSW samples under high NO x conditions was from aromatics due to a high emission factor of styrene from these samples. Fig. 5D shows that for sawdust, charcoal, cow dung cake and fuel wood samples, furanic species (17-58%) and aromatics (4-16%) were likely to provide a greater proportion of total SOA under low NO x conditions. The contribution of phenolic compounds was less (10-43%) due to larger aromatic and estimated furanic SOA yields under these conditions. Contributions remained small from oxygenated aromatics (3-11%), aliphatic species (0-2%), oxygenated aliphatics (0-2%), nitrogen containing NMVOCs (0.5-3%), monoterpenes (0-2%) and PAHs (5-15%). The contribution of aromatics to SOA from MSW remained high (43%). Bruns et al. (2016) showed that around 26% of SOA formed from the combustion of beech fuel wood was from phenolics. This was notably higher than the 5-9% contribution of phenol from Picea abies (spruce) reported by Hartikainen et al. (2018), who reported that 12-14% of the total SOA was from phenolic compounds. The results of this study appear more like that of Bruns et al. (2016), with between 10-70% of the total SOA from biomass combustion a result of phenolic compounds.

SOA formation potential
The yields under high NO x conditions for aromatics in this study (2-8%) were similarly low to those reported by Hartikainen et al. (2018) of 1.9-2.6% of the SOA from benzene and 1.9-3.3% from naphthalene, with low NO x conditions in this study suggesting aromatics could result in greater SOA yields. This study found relatively low SOA yields from monoterpenes from biomass sources (0-7%), which was like Hartikainen et al.
(2018) who found that monoterpenes contributed $1-3% to SOA. 24 This contrasted with Hatch et al. (2015) who showed that monoterpenes could result in 42-58% of SOA from black spruce and Ponderosa pine, however the fuel woods studied by Hatch et al. (2015) were likely larger emitters of monoterpenes. 11 It remained difficult to accurately characterise SOA yields from furanic species, as there is a lack of chamber simulation studies. This study suggested that furanic compounds could act as a large SOA precursor source, similar to several other studies. 11,24 In this study, the SOA yields of 2-methanol furanone, 2-(3H)-furanone, 5-hydroxymethyl-2[3H]-furanone, furfurals and methyl furfurals were estimated using the toluene yield, as a previous study by Gilman et al. (2015) indicated they had similar secondary organic aerosol formation potentials (SOAP). 21 This resulted in two different cases. Under high NO x conditions, the SOA yield in this study of furanics was 0.08, which was similar to that of by Hatch et al. (2015) 11 who used 0.10 based on the chemistry of 3-methyl furan measured from a previous study. 64 The low NO x yield used in this study was 0.33, which was like Bruns et al. (2016), who used a furfural yield of 0.32 based on the average SOAP of all assigned $C 6 compounds. The true SOA yields from furanic species from domestic biomass burning samples remained uncertain and requires further chamber studies. This issue was previously highlighted. 65 While following a different approach, this study arrived at similar estimated yields of furanic compounds as those used previously. It highlighted that SOA formation from domestic biomass burning smoke from solid fuels collected in India was predominantly driven by phenolic and furanic compounds as well as aromatics. Table 2 shows the mass fraction of NMVOCs released which were identied as SOA precursors from yield data, and the mass fraction of NMVOCs which resulted in SOA under high and low NO x conditions. These were presented as a mass fraction of total organic compounds measured during this study. Some sources, such as cow dung cake and MSW, released signicantly more NMVOCs per kg of fuel burnt than fuel wood (MSW $ 88 g kg À1 , cow dung cake $ 62 g kg À1 and fuel wood $ 19 g kg À1 ). 14 Multiplying the emission factor by the mass fraction of NMVOC which would result in SOA highlighted interesting differences in SOA production between different source types. Table 2 shows this result, with the mass of SOA which would result per kg of fuel burnt under high (SOA h , g kg À1 ) and low (SOA l , g kg À1 ) NO x conditions. The amount of SOA produced by each source was considered relative to fuel wood, due to difficulties establishing SOA precursor from the chamber background for LPG. Emissions from cow dung cake and MSW resulted in $3-4and $6-7-times greater SOA per kg of fuel burnt than fuel wood, respectively. The result for SOA formation from cow dung cake combustion was consistent with that calculated by Fleming et al. (2018), who reported a factor of 3 times higher SOA from dung-chulha stoves than brushwood-chulha for samples collected from Haryana, India. 66 It is also noteworthy that SOA estimated from chamber yield data and that observed experimentally during burning experiments have been shown to agree within a factor of 2. 22 The estimates of SOA formation should be considered relative to the heat output of specic fuels. Energy densities have been reported for LPG (45 837 kJ kg À1 ), charcoal (25 715 kJ kg À1 ), acacia fuel wood (15 099 kJ kg À1 ), Eucalyptus spp fuel wood (15 333 kJ kg À1 ), rice straw (13 027 kJ kg À1 ), Brassica spp (11 763 kJ kg À1 ) and dung cakes (11 763 kJ kg À1 ). 67 This highlights that whilst all sources were likely to result in SOA production, the burning of fuels such as cow dung cake is inadvisable due to the low caloric value and high emission factor. The result is that more fuel is required to be burnt to achieve the same heat output, which will lead to greater levels of NMVOC emission. These will subsequently degrade local and regional air quality through the formation of a greater quantity of secondary pollutants. Fig. 6A shows that LPG OH reactivity was principally driven by alkanes ($75%). The contributions of other species were small and may have arisen from difficulties in background correction for this low emission fuel. For charcoal, the reactivity with OH was principally caused by furanics (33%), phenolics (19%) and oxygenates (15%). The reactivity of fuel wood emissions with OH was driven by furanics (34%), oxygenates (27%), phenolics (13%) and alkenes (12%). Emissions from cow dung cake with OH were due to by furanics (32%), oxygenates (21%), alkenes (16%), phenolics (12%) and nitrogen containing NMVOCs (11%). The OH reactivity from crop residue was from furanics (38%), oxygenates (23%), phenolics (14%) and alkenes (11%). For sawdust, reactivity with OH was a result of furanics (34%), oxygenates (24%), phenolics (15%) and monoterpenes (9%). However, for charcoal and sawdust only 2 samples were measured. The OH reactivity from MSW samples was different and a result of aromatics (30%), followed by oxygenates (22%), furanics (19%) and phenolics (5%).

OH reactivity
This study identied the species with the largest reactivity with the OH radical from Indian solid fuels. Ozone production from emissions when these fuels are combusted will be more complex and ultimately depend on NMVOC/NO x ratios, meteorology and solar radiation. 25 Whilst the phenolic compounds here show relatively large contributions to OH reactivity (5-19%), their emission may both decrease or increase ambient O 3 . In NMVOC limited environments in India, a reduction in NO x could lead to an increase in O 3 . 68 In NO x limited environments this may result in negative O 3 formation due to the formation of nitrophenols, which reduces the amount of NO 2 available for NMVOC oxidation. 69 Gilman et al. (2015) calculated the relative contribution of different functionalities to the OH reactivity of fuel types from the U.S. 21 The fuel types studied by Gilman et al. (2015) showed that alkenes contributed 25-29% of the OH reactivity, which was larger than found in this study (7-16%) for Indian fuels. The contribution to OH reactivity of OVOC for U.S. fuels (41-54%) was less than found in this study (45-76%). The contributions of monoterpenes for fuels from the U.S. were slightly larger (4-14%) than for those from India (0-7%). This was likely due to a greater contribution of monoterpene emitting fuel woods, such as pine, to fuels from the U.S. studied by Gilman et al. (2015). Both studies found a small contribution of aromatics (<5%) and nitrogen containing NMVOCs (<11%) to OH reactivity. Fig. 6B shows the OH reactivity of each source at the top of the combustion chamber relative to LPG. This was calculated by multiplying the mean OH reactivity of ue gases by the volume of air sampled and normalising to the total reactivity of LPG. The OH reactivity of LPG was the lowest. Emissions from charcoal, fuel wood, crop residue, cow dung cake and sawdust were respectively $8, 30, 90, 120 and 150 times more reactive with OH than those from LPG. The OH reactivity of emissions from MSW was the greatest and approximately 230 times greater than from LPG. Fuel wood, cow dung cake and MSW burning are large NMVOC sources in India. 70,71 The signicantly greater OH reactivity of emissions from these sources is likely to substantially deteriorate local and regional air quality compared to users cooking over LPG.

PAH toxicity
A need has been identied to better understand the impact of PAHs from combustion sources in cities such as Delhi, where concentrations have been shown to be high and suggested to be enhanced by emissions from burning sources. 72 Fig. 7 shows that when comparing the toxicity of 21 PAHs released, fuel wood, crop residue, cow dung cake and MSW were respectively 20, 60, 130 and 220 times more toxic than LPG per kg of fuel burnt. Toxic emissions from these 21 PAHs released from LPG were small, and were principally driven by naphthalene (43%), uoranthene (24%) and methylnaphthalenes (11%). The largest drivers for fuel wood/crop residue toxicity were benzo[a]pyrene (38%/48%), naphthalene (14%/11%) and benzo[b]uoranthene (8%/8%), respectively. The contribution of naphthalene to the toxicity of cow dung cake and MSW was lower, with their toxicities driven by benzo[a]pyrene (49%/42%), dibenz[a,h] anthracene (13%/16%) and benzo[b]uoranthene (8%/13%). The real-world effect of this toxicity would be signicantly enhanced for fuel wood and cow dung cake, by around a further factor of 10. This is because signicantly more fuel wood and cow dung cake fuel is used per user than LPG, due to the higher energy density of LPG and more efficient burning of this fuel. 73,74 These results reinforced ndings of other studies assessing the health benets of LPG vs. solid fuels which suggested that to signicantly reduce the impacts of combustion, a shi to cleaner cooking technologies was required. [75][76][77]

Conclusions
This study compiled comprehensive measurements of organic emissions from the combustion of a range of domestic fuels common to India. A range of detailed and complementary techniques allowed a VBS to be generated across a wide range of C* values. This highlighted that IVOC emissions should be better represented in models for an improved understanding of SOA production from emissions caused by domestic solid fuel combustion.
The results estimated that phenolics and furanics were important to both the SOA production potential and the OH reactivity, respectively accounting for 10-70% and 9-58% of the SOA production potential and 5-22% and 9-48% of the OH reactivity of domestic biomass burning emissions. The contribution of smaller oxygenated species to the OH reactivity was also signicant at 15-42%. Different combustion sources were compared, which showed that emissions from fuel wood, crop residue, cow dung cake and MSW burning were 30, 90, 120 and 230 times more reactive with the OH radical and that PAH emissions were 20, 60, 130 and 220 times more toxic than LPG, respectively. This also showed that NMVOCs released from the combustion of cow dung cake and MSW samples in this study resulted in $ 3-4 and 6-7 times more SOA production per kg burnt than fuel wood, respectively. This demonstrated that reduction of emissions from these sources is important to improve local and regional air quality across India.
Limited measurements were made from MSW, cow dung cake, crop residue and LPG samples in this study. Future studies are needed to improve the emission proles of these domestic fuels to better understand the impact of burning emissions. This study also highlights that the C* of many species measured, including alkanes, remain uncertain and future studies are required to better understand the C* of these species. In addition, there have only been a limited number of chamber studies to determine the SOA formed during the oxidation of furanic species under high and low NO x conditions. More studies of the oxidation, and subsequent SOA formation, of these important domestic biomass burning emissions are required to better understand the impact of domestic solid fuel use on the atmosphere.

Author contributions
GJS developed volatility basis dataset, analysed I/SVOC fractions, conducted inventory comparison, developed SOA and OH reactivity models and lead paper. BSN, JRH and WJFA assisted with data analysis. SSMY contributed code to map of organics onto VBS. EN, ARR, JDL, TKM, RG and LKS contributed to data interpretation. JFH provided overall guidance for analysis and writing of paper. All authors contributed to the discussion, writing and editing of the manuscript.

Conflicts of interest
The authors declare that they have no conict of interest. drawn in the paper are based on interpretation of results by the authors and in no way reect the viewpoint of the funding agencies. GJS and BSN acknowledge the NERC SPHERES doctoral training programme for studentships. TKM is thankful to Director, CSIR-National Physical Laboratory, New Delhi for allowing to carry out this work. LKS acknowledges Physical Research Laboratory (PRL), Ahmedabad, India for the support and permission to deploy PTR-ToF-MS during the experimental campaign.