Interactive comment on “ The nitrogen budget of laboratory-simulated western U . S . wildfires during the FIREX 2016 FireLab study ”

The authors present an overview of the N-budget during the FIREX campaign, bringing together datasets from a variety of instruments and techniques (outlines in Table 1) in already peer-reviewed articles of FIREX work. This analysis is certainly a great contribution to the field and specifically to the science of reactive nitrogen budgets from forest fires. The methods used are state-of-the-art, including the instrumentation used for acquiring the data during the campaign in addition to the data analysis, including PMF. I command the authors on putting this data set together and anticipate this contribution, upon revision, to be a go-to paper for the atmospheric organic-N and inorganic-N fate

Abstract. Reactive nitrogen (N r , defined as all nitrogencontaining compounds except for N 2 and N 2 O) is one of the most important classes of compounds emitted from wildfire, as N r impacts both atmospheric oxidation processes and particle formation chemistry. In addition, several N r compounds can contribute to health impacts from wildfires. Understanding the impacts of wildfire on the atmosphere requires a thorough description of N r emissions. Total reactive nitrogen was measured by catalytic conversion to NO and detection by NO-O 3 chemiluminescence together with individual N r species during a series of laboratory fires of fuels characteristic of western US wildfires, conducted as part of the FIREX Fire Lab 2016 study. Data from 75 stack fires were analyzed to examine the systematics of nitrogen emissions. The measured N r / total-carbon ratios averaged 0.37 % for fuels characteristic of western North America, and these gas-phase emissions were compared with fuel and residue N/C ratios and mass to estimate that a mean (±SD) of 0.68 (±0.14) of fuel nitrogen was emitted as N 2 and N 2 O. The N r detected as speciated individual compounds included the following: nitric oxide (NO), nitrogen dioxide (NO 2 ), nitrous acid (HONO), isocyanic acid (HNCO), hydrogen cyanide (HCN), ammonia (NH 3 ), and 44 nitrogen-containing volatile organic compounds (NVOCs). The sum of these measured individual N r compounds averaged 84.8 (±9.8) % relative to the total N r , and much of the 15.2 % "unaccounted" N r is expected to be particle-bound species, not included in this analysis.
A number of key species, e.g., HNCO, HCN, and HONO, were confirmed not to correlate with only flaming or with only smoldering combustion when using modified combustion efficiency, MCE = CO 2 /(CO + CO 2 ), as a rough indicator. However, the systematic variations in the abundance of these species relative to other nitrogen-containing species were successfully modeled using positive matrix factorization (PMF). Three distinct factors were found for the emissions from combined coniferous fuels: a combustion factor (Comb-N) (800-1200 • C) with emissions of the inorganic compounds NO, NO 2 , and HONO, and a minor contribution from organic nitro compounds (R-NO 2 ); a high-temperature pyrolysis factor (HT-N) (500-800 • C) with emissions of HNCO, HCN, and nitriles; and a low-temperature pyrolysis factor (LT-N) (< 500 • C) with mostly ammonia and NVOCs. The temperature ranges specified are based on known com-

Introduction
Wildfires have severe impacts on the chemistry of the atmosphere from local to global scales (Crutzen and Andreae, 1990). A warmer, drier climate in western North America, coupled with policies that have allowed buildup of fuels in forest ecosystems has led to increases in frequency and severity of wildfires in this region (Abatzoglou and Williams, 2016;Westerling et al., 2006). The new strategy for management of wildfire in the US is to allow fire where possible and to fight fire where needed (D. C. . The science behind making these decisions and understanding their consequences involves, in part, a better understanding of the emissions from wildfires. The NOAA FIREX (Fire Influence on Regional and Global Environments Experiment) Fire Lab experiment was conducted in the fall of 2016, at the US Forest Service Fire Sciences Laboratory in Missoula, Montana, to acquire detailed measurements of particle and gas-phase emissions from fires involving fuels characteristic of the western US (NOAA, 2020a). Several aspects of these measurements dealing with VOC species and individual reactive nitrogen species (N r , defined as all nitrogen compounds except for N 2 and N 2 O) have already been published Manfred et al., 2018;Sekimoto et al., 2018;Selimovic et al., 2018;Zarzana et al., 2018), including emissions factors for many of the N r species .
The N r compounds emitted by natural-convection biomass burning (BB) arise solely from the N in the fuels, since the combustion temperatures are not high enough (< 1200 • C) to produce oxides of nitrogen (NO x ) from N 2 and O 2 (the so-called Zeldovich or thermal nitrogen cycle) (Lobert and Warnatz, 1993;Taylor et al., 2004;Wotton et al., 2012). The fuel nitrogen cycles that pertain to BB flaming combustion are shown schematically in Fig. 1 (Glarborg et al., 2018;Lobert and Warnatz, 1993;Lucassen et al., 2012;Manion et al., 2015). Note that the equations shown in Fig. 1 are meant to indicate the general flow of the chemistry and do not always convey the mechanistic subtleties of the reactions, which are fully covered in specialized references (Glarborg et al., 2018;Manion et al., 2015). N r compounds are emitted as small molecules, hydrogen cyanide (HCN), isocyanic acid (HNCO), and ammonia (NH 3 ) resulting from pyrolysis of the fuel, with minor contributions from larger N-containing organic species, especially at lower temperatures. Flame chemistry converts those species to N 2 , N 2 O, nitric oxide (NO), nitrogen dioxide (NO 2 ), and nitrous acid (HONO) as a result of radical chemistry. It has been recognized for some time that a significant amount of denitrification (conversion of N r compounds to N 2 ) occurs due to reactions of NO with NH i (where i = 1, 2, or 3) or N atoms, as confirmed experimentally . While N atoms are also intermediates in the thermal NO x cycle, and the reaction N + O 2 → NO + O figures in to both the fuel and thermal NO x cycles, the second reaction of the thermal NO x cycle, O + N 2 → NO + N, is too slow at BB flame temperatures to result in NO x production (Manion et al., 2015). In addition to the small molecules shown in Fig. 1, numerous N r compounds are emitted in roughly the following categories: amides, amines, heterocyclic compounds, nitriles, isocyanates, and nitro compounds (Andreae, 2019;Andreae and Merlet, 2001;Koss et al., 2018;Lobert et al., 1991Lobert et al., , 1990Lobert and Warnatz, 1993;Stockwell et al., 2015). These compounds are produced at much lower abundance from fuel pyrolysis and partial reactions with the radical species in Fig. 1.
The emissions of N compounds from BB and wildfires in general have been the subject of considerable research. Early studies by Lobert et al. (1990Lobert et al. ( , 1991 and Lobert and Warnatz (1993) measured a wide range of N r compounds in laboratory fires and suggested that considerable denitrification (conversion of fuel nitrogen to N 2 ) was taking place, a process later confirmed in experiments described by Kuhlbusch et al. (1991). Subsequent work on laboratory fires has betterdefined particle-phase nitrogen emissions (McMeeking et al., 2009) and led to the recognition of the importance of several inorganic N r species, such as HONO and HNCO Roberts et al., 2011;Veres et al., 2010), and the presence of a wider variety of organic N r species (Gilman et al., 2015;Stockwell et al., 2015;Warneke et al., 2011). A number of studies have sought to summarize both realworld and laboratory emissions of N r compounds (Akagi et al., 2011;Andreae, 2019;Andreae and Merlet, 2001;Coggon et al., 2016;Yokelson et al., 2013bYokelson et al., , 2009. The known N compounds range in oxidation state from NH 3 to HNO 3 and include N 2 and N 2 O. Among the more prominent and important N r species are the following: NO x (NO and NO 2 ), which is a key player in the atmospheric oxidant cycle; NH 3 , which has a major role in particle formation; HONO, which can be an important radical source; HCN and acetonitrile (CH 3 CN), which are toxic at high concentrations and represent valuable Figure 1. Schematic of the combustion chemistry of the small molecules that are emitted from BB and represent sources and sinks of reactive nitrogen (N r ), adapted from Glarborg et al. (2018), Lobert and Warnatz (1993), Lucassen et al. (2012), and Manion et al. (2015). H 2 O (surf) denotes the combination of H 2 O and a surface to facilitate the reaction. Red color indicates the highest temperature (combustion) processes; blue indicates intermediate temperature processes, and green indicates the lowest temperature processes. The species that are measured in this work are shown in bold and slightly larger text. tracers for following fire transport; and isocyanates, HNCO and methyl isocyanate (CH 3 NCO), which have unique health impacts . In addition, nitro (R-NO 2 ) or nitrogen heterocyclic compounds may contribute to so-called brown carbon, aerosol organic compounds exhibiting optical absorption in the near-UV or blue wavelength regions. Wildfire N emissions also have very minor contributions from gas-phase nitric acid (HNO 3 ). Nitric acid is either not efficiently produced by BB or readily incorporated into aerosol if it is produced in fresh wildfire plumes, as is clear from the absence of HNO 3 enhancements in several studies of BB plumes (Liu et al., 2016;Yokelson et al., 2009;Alvarado et al., 2010); however nitrate (NO − 3 ) has been shown to contribute to aerosol mass particularly for inefficient combustion (May et al., 2014). Flame chemistry is inefficient in forming N 2 O, relative to the pathways that form N 2 (Andreae, 2019;Andreae and Merlet, 2001;Griffith et al., 1991;Hao et al., 1991). The modeling of the emissions of these N compounds on a large scale could benefit from a better understanding of the total budget of these species as a function of fuel nitrogen content and the dependence of the individual species on fuel type and combustion conditions. The construction of N r budgets in this work is made possible by the inclusion of a total reactive nitrogen measurement (termed N r herein), a method by which all nitrogen compounds besides N 2 and N 2 O are converted to NO and detected by NO-O 3 chemiluminescence. This technology has been developed by a number of groups, typically using precious metal or NiCr catalysts that have been shown to convert all N r compounds to NO (and to some extent NO 2 ) at high temperatures (750-825 • C) (Hardy and Knarr, 1982;Kashihira et al., 1982;Marx et al., 2012;Roberts et al., 1988). There are also commercial instruments that incor-porate this technology (see for example Thermo Scientific Model 17i). This technique has been applied to gas-phase atmospheric measurements, principally to measure NH 3 by difference techniques (Saylor et al., 2010;Schwab et al., 2007), and has also been used to observe wildfire plumes that have impacted ambient air measurements (Benedict et al., 2017;Prenni et al., 2014). We have recently developed a platinum/molybdenum oxide N r catalyst system and confirmed that it quantitatively converts N r compounds including all particle-bound nitrogen compounds . To our knowledge this technique has not been applied directly to BB emissions before. This paper describes the total reactive nitrogen and individual N r compound measurements made during the Fire Lab 2016 experiment, with the intent of providing information that can be used for analysis and modeling of the impact of wildfire emissions on the atmosphere. The total N r measurements are combined with CO 2 , CO, and VOC measurements and fuel, residue, and ash C and N content to estimate the amount of N lost to N 2 and N 2 O. In addition systematics of the ratio N r / total carbon are examined for simplifying relationships. Fire-integrated N r is then compared to fireintegrated measurements of individual compounds to determine the fraction of unaccounted-for N r . The systematic behavior of individual N r species and their fractional contribution to N r are examined with respect to fuel type, N content, and combustion processes. A positive matrix factorization (PMF) technique is used to examine commonalities between fires of different fuels under different conditions and compared to the PMF analysis of the VOC emissions published by Sekimoto et al. (2018). The results are used to arrive at suggested guidelines that can be used estimate N r emissions profiles for fires representative of western North America.    a The OP-FTIR has a 10 ppbv detection for gas-phase HNO 3 , but HNO 3 was not observed above detection limit. b Ethylamine, methanimine, propeneamine, sulfinylmethanamine, trimethylamine, buteneamines. c Formamide, acetamide, methylmaleimide. d C 2 -pyrroles, dihydropyridine, ethynylpyrrole, methylpyridine, methylpyrrole, pyridinealdehyde, 4-pyrindinol, vinylpyridine. e Acetonitrile, acrylonitrile, benzonitrile, butanenitrile, butynenitrile, benzoacetonitrile, C 7 acrylonitrile, C 8 -nitriles, heptylnitrile, furancarbonitrile, methylbenzoacetonitrile, pentylnitriles, propanenitrile, propynenitrile, butenenitrile, methylisocyanate. f Butenenitrates, nitrobenzene, nitroethane, nitroethene, nitrofuran, nitromethane, nitropropanes, nitrotoluene.

Methodology
The Fire Lab 2016 study involved laboratory burns of fuels mostly characteristic of western North American wildfires such as coniferous fuels and chaparral fuels (important in central to southern California and the southwestern US). We also measured some that have global significance such as Indonesian peat and yak dung (important in areas above timberline or where wood is scarce, such as India, Nepal, and Tibet). The procedures and associated details of the study have been described previously by Selimovic et al. (2018) and will be only briefly summarized here. The detailed data on fuel types, amounts, and composition can be found in Table S1 in the Supplement and in the Supplement of Selimovic et al. (2018). The laboratory burns involved fuel samples, ranging in mass from 0.26 to 6.02 kg spread out on a fuel bed roughly 1 m × 1 m square. Fires were started without the addition of any contaminants, using an electric igniter (a series of NiCr heating elements that were flash-heated electrically), and typically lasted from approximately 5 to 30 min.
Seventy-five fires were conducted in the "stack" burn configuration where the smoke was directed up the central stack of the facility where it could be sampled simultaneously by all the instruments that measured gas-phase species and some of the particle-phase measurements. The sampling platform was about 15 m above the fire, and the sampling took place in well-mixed smoke approximately 5 s after emission (Christian et al., 2003). Thirty-one additional fires were conducted as "room" burns on most of the same fuels, when the stack was closed and the room was allowed to fill with smoke, permitting sampling to be done over the course of several hours. The following analyses will focus on the stack burns, as those measurements had little or no interferences from surfaces, where room burns are known to be compromised by the loss of materials, such as NH 3 , to the room walls at long sample times (Stockwell et al., 2014). Ash analyses were performed only on the residues from the room burns, and those values will be used for the N and C budget calculations, with the assumption that stack and room burns left similar ash, considering the combustion conditions were the same for each type of fire. Table 1 lists the compounds and associated techniques used to measure them during the Fire Lab 2016 study and describes the grouping of NVOCs measured by a proton transfer reaction time of flight mass spectrometer (PTR-ToF-MS) into common categories, e.g., amines, nitriles, etc. We specifically note that the OP-FTIR is capable of measuring gas-phase HNO 3 with comparatively good sensitivity (10 ppbv detection limit in fires where N r can be 5 ppmv or more), but HNO 3 was not observed above detection limit in any of the fires.

N r and NO measurements by chemiluminescence
Total reactive N (N r ) was measured by catalytic conversion to NO, followed by O 3 chemiluminescence using an instrument described previously (Williams et al., 1998). N r and NO were sampled from inlets inserted adjacent to the inlet-less openpath Fourier transform infrared spectrometer (OP-FTIR) instrument path during the stack burns  and from a platform approximately 4 m off the floor in the middle of the room during the room burns. The catalyst used for the N r channel, described in detail by Stockwell et al. (2018), consisted of a 11 mm i.d. quartz tube, packed with 36 platinum screens, heated to 750 • C. This tube was wrapped with high-temperature heating tape and insulated inside a 7 cm o.d. stainless steel tube that was fitted to a bulkhead placed through the wall of the stack. The N r channel was diluted by a factor of 5 : 1 (±3 %) using a flow of zero air added immediately downstream of the Pt catalyst assembly. NO was sampled through a 6.3 mm o.d. stainless steel inlet tube, which was placed through the bulkhead directly into the free air stream of the stack and connected to a 50 mm Teflon filter holder immediately outside the stack. The transfer lines for the N r and NO measurements consisted of 6.35 mm o.d., 1 mm wall thickness PFA tubing of approximately 20 m in length. N r and NO data were acquired at 1 s frequency, but the flow rate through each inlet was 1 SL min −1 , resulting in residence time in each inlet of 14 s. This time delay was corrected in the data analysis. Any chemical effects of the inlet on the sampled air stream were negligible since the analytes consisted of only NO and NO 2 , and those are known to be transmitted by PFA Teflon tubing with essentially no surface effects. However, there were possible effects of the inlets on the temporal features of the measurement through diffusion or turbulent mixing. Those effects were examined through comparison of the temporal variations in the NO signal with the NO measured by the OP-FTIR and comparison of the N r signal under smoldering conditions with the NH 3 measured by the OP-FTIR. Both of these comparisons showed that the NO and N r inlets had effective time constants of 4 s, somewhat slower than the diffusive relaxation time assuming solely laminar flow. Examples of the estimate of diffusion and dispersion on NO and N r signals and the estimate of the effective time constant of these measurements are presented in the Supplement.
The inlet streams were sampled by the NO instrument either directly (NO channel) or after passing through a second catalyst of molybdenum oxide (MoOx) to convert remaining NO 2 to NO. The MoOx catalyst consisted of a molybdenum tube at 350 • C to which a small flow of H 2 (0.8 % v/v) was added to control the redox state of the surface. Both channels of the instrument were "detuned" to keep raw photon count rates below 4 MHz, by turning down the O 3 flows and photo-multiplier tube (PMT) voltages. Calibrations were performed with both a NO standard in N 2 (Scott-Marrin) and 10.1 ppmv standard of HCN in nitrogen (Gasco). The Pt catalyst was dismounted from the stack (or room) every few days and checked for conversion efficiency by the addition of the HCN standard to the inlet. Conversion efficiencies were found to be consistently high (> 98 %) throughout the entire sampling period (5 October-12 November 2016). There were slight background signals (a few tens of ppbv) for both NO and N r in both the stack and room air prior to and after the burns, and those were subtracted from the fires' signals prior to reporting the data. The overall uncertainties in the NO and N r data were ±10 % for each measurement.

Other measurements
Measurements of individual species during the 2016 Fire Lab study have been presented in several previous publications. The OP-FTIR measurements were discussed by , and the PTR-ToF measurements were discussed by Koss et al. (2018). In addition, some of the calibration methods and GC separation and identifications rely on additional analytical work presented by Sekimoto et al. (2017) and Gilman et al. (2015). We measured the mass and elemental content of the initial fuel and the mass of unburned fuel for all the fires, and we measured the mass and the elemental

PMF analysis
Trace gas measurements from multiple instruments involved in the Fire Lab study were combined and analyzed using positive matrix factorization (PMF). PMF is a numerical method that was used in this case to partition the compounds involved in a time-varying mixture of chemicals into a few groups, or "factors", where a compound can appear in more than one factor. A factor represents a consistent profile of compounds that is representative of one of the sources contributing to the total signal. The sum of all the factors then ideally describes the total composition of the measurements, which in this case is the emissions of N r compounds. By its nature, PMF assumes that the total signal is a linear combination of individual sources that have a consistent composition, the relative contribution of which is represented by the amount of each compound or category found in each factor (Paatero and Tapper, 1994;Ulbrich et al., 2009). We hypothesize that species with dominant fractions in the same factor are related to each other via the same formation processes. With knowledge of factor composition and the amount of each factor at any given time, the original emissions measurements can be reconstructed, and this approach provides an alternate source of profiles for fire emissions. PMF has also been used by a number of groups to explore how much various source profiles contribute to complex ambient measurements (see for example Ulbrich et al., 2009)  The application of PMF to this data set is different than the instances where it is applied to data from a single instrument in which compound abundances are inherently scaled properly and error estimates are well defined and self-consistent. For example, when applied to mass spectral data from a sin-gle instrument, errors can be expected to scale as the square root of ion counts based on fundamental counting statistics . In this work we are including nitrogen measurements from several instruments, and thus we chose to use mixing ratios as the unit of comparison. The error estimates required by the PMF analysis were taken from the reported combined uncertainties: the sum of the detection limit plus the estimated random error of the measured value. For example, the uncertainty in a NO mixing ratio of 500 ppbv was ±51 ppbv. The variables that were used in this PMF analysis and their units and corresponding errors are listed in Table 2. Where compound categories are specified (e.g., nitriles), the values were the sum of the measured compounds in that category as listed in the footnotes of Table 1. The data were further adjusted by subtracting the ambient air background before and after the fires, which was a relatively minor adjustment for most compounds and categories. Any negative numbers that resulted were very small compared to the fire emissions and were set to zero. In addition to the PMF analysis for the species listed in Table 2, several exploratory runs were tried with CO 2 , CO added (in units of ppmv), and total N r (in units of ppbv) added to the list in Table 2. No significant differences were observed in the results for individual N r compounds and classes, so CO 2 , CO, and total N r were not included in this analysis.
We applied PMF to single-fire data as well as extended time series that included all fires of a particular fuel type, in line with the approach laid out by Sekimoto et al. (2018). By consolidating fuels from a particular vegetation type, the fire-to-fire variability largely driven by differences in the fuel (e.g., moisture content, structure, quantity) is constrained and the most representative fire conditions are captured. Two fuel groups were analyzed in this way: the western US coniferous ecosystem fuels, which included ponderosa pine, lodgepole pine, Douglas fir, Engelmann spruce, and subalpine fir, and the southwestern US chaparral ecosystem, which was represented by chamise and manzanita. The consolidated time series for the coniferous ecosystems included realistic mixtures, canopy only, and litter only, while duff and rotten logs were analyzed separately and not included in the time series.

Results and discussion
The measurements of total N r can be combined with N and C measurements of fuel and ash to estimate N lost to N 2 and N 2 O. The total N r emitted from laboratory fires combined with individual N compound measurements allow us to construct a budget for N r species that define what the dominant forms of N are and how those emissions depend on other fire parameters or different temperature combustion processes. The systematics of N emissions found by PMF are compared to other fire indicators, and PMF analyses previously conducted on VOCs allow the formulation of simplifying relationships that can be used in atmospheric models of wildfires.

Comparison of N r and total carbon in fire emissions
The total N r and total carbon emissions were measured for 75 stack fires in order to place the N emissions in the context of total carbon, which has been widely estimated for wildfires. Example time series of NO, N r , CO, and CO 2 (CO and CO 2 corrected for their backgrounds) are shown in Fig. 2, for a fire burning a sample of ponderosa pine realistic mix (Fire 004). In addition to the chemical species, the modified combustion efficiency (MCE) was also plotted. MCE is defined as where CO 2 and CO are the CO 2 and CO levels above the ambient. MCE has traditionally been used to indicate the relative amount of flaming and smoldering combustion in a fire. The time series for Fire 004 (Fig. 2) shows a short initial smoldering/distillation phase (MCE 0.7 to 0.8) as heat pyrolyzes the fresh fuel and releases VOCs from existing pools in the fuel, followed after ignition by a relatively efficient mix of flaming and smoldering combustion (MCE 0.95 to 0.98) and then finally a subsequent period of essentially pure smol-dering (MCE ∼ 0.80). The N r and NO timelines had many features in common because NO is often the most abundant N r compound (see below). As a result, it is useful to compare the quantities N r −NO and (N r −NO) / N r to the other measures of chemical species or combustion efficiency. As expected, (N r −NO) / N r , in Fig. 2c is anticorrelated with MCE since N r is primarily NO at high MCE. In addition to the anticorrelation, this non-NO fraction, like its approximate carbon analog CO/CO 2 , has a wider dynamic range than MCE and will often suffer less from background variability than carbon-based indices (Yokelson et al., 2013a). The concentration profiles of the background-corrected measurements of N r , CO 2 , CO, and all the carbon-containing species measured by the FTIR  during the stack burns were integrated over the entire time of the burn to obtain total carbon, termed TC here, and total N r . The additional carbon species included methane and a number of other gas phase VOCs as well as organic-and black-carbon aerosol. Altogether these carbon species should account for ≥ 98 % of emitted carbon (McMeeking et al., 2009). Total N r is plotted in Fig. 3, versus TC (Fig. 3a) and versus nitrogen burned, which is calculated from the %N in the fuel times the mass of fuel consumed (Fig. 3b). The points in Fig. 3 are colored by the fuel N/C mole % obtained from the elemental analysis of each fuel. The relationship between N r and TC in Fig. 3a clusters around the 0.37 % line, and those points are from fuels most characteristic of the North American biomes impacted by wildfire. There are clear outliers in the correlation of N r and TC; for example, yak dung and two samples of duff (Engelmann spruce and subalpine fir) were high due to the fact that they either have high fuel N/C ratios (dung; see Table S1 in the Supplement) or burned with minimal flaming (whole-fire MCEs 0.86-0.89, duff and dung), hence experienced less denitrification. The fuels that were low in N r /TC in Fig. 3a, ponderosa pine rotten log, subalpine fir, and excelsior, had low fuel N/C, so when plotted versus nitrogen burned in Fig. 3b, they cluster with the main group of characteristic fuels, i.e., they are no longer "outliers" in the distribution.

Estimates of denitrification
The removal of N to forms that are inactive in the troposphere, N 2 and N 2 O, has importance in the biogeochemistry of forest ecosystems and also determines how much N takes part in wildfire plume chemistry. The points in Fig. 3a are all lower than the corresponding fuel N/C mole ratio, due to the denitrification chemistry, shown in Fig. 1 and verified in lab studies described by Kuhlbusch et al. (1991), and the production of N 2 O, which is also not measured by the N r technique. The sum of N 2 and N 2 O produced in the fires can be estimated from the difference between the fuel N/C and the N r /total C emitted and the data on C and N content remaining in the ash. The mass balance equations used for these calculations are detailed in the Supplement.  The distribution of the N lost to N 2 and N 2 O is shown in Fig. 4. Chemical analyses were not done for all fuels during the stack burns, and the analysis above assumes that the ash residues and ash / burned fuel ratios from the stack burns were well represented by those for the same fuels used in the room burns, for which mass yields and chemical analyses were done. Data are missing for fuels that did not have a corresponding residue analysis. The median fraction of N lost to N 2 and N 2 O for ash-corrected fires was 0.70, and the mean (± standard deviation) was 0.68 (±0.14). This fuelbased estimate is uncertain by approximately 25 % because of the above assumptions concerning the applicability of the residue analyses from the room burns and because fuel moisture corrections were assumed to apply to all of the materials burned, foliage vs. woody biomass (see the Supplement for details). The emission of N 2 O relative to N 2 is approximately 10 % or less for a wide range of fuels (Andreae, 2019). Assuming the N remainder in our work is at least 90 % N 2 gives values that are somewhat higher than the N 2 values reported by Kuhlbusch et al. (1991) where N 2 accounted for 36 % of fuel N burned in flaming-stage fires. A closer inspection of Kuhlbusch et al. (1991) showed a range of N 2 yields of 40 %-54 % at highest MCEs of 0.94-0.97. Possible reasons for these differences are that the  fires were limited to grasses, hay, and pine needles, and the fires were confined to a closed container and so may not have experienced the convection and turbulence of typical biomass fires. In addition, the fires analyzed in our work were somewhat weighted towards the full canopy and higher temperature burning fuels, since ash analyses were not done for peat, dung, and many of the "litter" samples, all of which tend to burn less efficiently. Goode et al. (1999) estimated an N 2 emission of 45 ± 5 % for MCE values of 0.95 in grass and surface fuels. The range of values determined in our work overlap with these literature values but are on average higher. It should be noted that such loss of reactive nitrogen can have implications for ecosystem N budgets, as discussed by Kuhlbusch et al. (1991).

The budget of N r and individual N-containing species
The composition of the N that does not get converted to N 2 or N 2 O is of intense importance in determining atmospheric impacts of fires, since those compounds are involved in oxidation capacity (NO x ), radical production (HONO), and particle formation (NH 3 ). Emission factors (EF, defined as mass of compound emitted per mass of fuel burn) for all the individual N r compounds identified in our work have been com-piled and reported in previous publications Selimovic et al., 2018), so this paper will focus on the N r budget. The balance of N r budget for Fire 047, subalpine fir realistic mix, is shown in Fig. 5, in which the timelines of N r , NO, N r −NO, sumN, and NVOC are plotted along with MCE and (N r −NO) / N r . The quantity sumN is the sum of all other non-NO compounds, and NVOC is the subset of sumN that are organic nitrogen compounds measured by the PTR-ToF, as listed in Table 1. This fire had a mixture of flaming and smoldering combustion throughout the fire as indicated by MCE and nitrogen profiles (Fig. 5d). The comparison of N r −NO with the sumN in Fig. 5b shows that much of the N is accounted for. The major contributors to sumN for this fire were HNCO, HCN, HONO, NO 2 , and NH 3 , while NVOC was a very small contributor to sumN (Fig. 5b). Note that while HNO 3 is measurable by FTIR with good sensitivity, no HNO 3 signals were observed above detection limit, which was 10 ppbv. Fig. 5c shows the residual left after NO and sumN are subtracted from N r , corresponding to an integrated amount of 15.6 ± 8 % of N r . This residual is reasonable, considering typical published particle N r measurements (Akagi et al., 2012(Akagi et al., , 2011Liu et al., 2017;May et al., 2014), and consistent with there being some particle N r from flaming, which is most likely organic nitrates or nitro-organics, and particle ammonium from smoldering with potassium or ammonium nitrate, potentially accounting for substantial N r . Several fuels had much lower NO emissions and higher unaccounted-for N r . Yak dung was one such fuel, the emissions of which stand in contrast to the fire shown above. The nitrogen emissions from Fire 050, yak dung, are shown in Fig. S2 in the Supplement. This fuel produced mostly smoldering emissions as exemplified by the low NO levels relative to N r (panel a) and the low MCEs observed (panel d).
The sum of N r species was somewhat correlated with the quantity N r −NO but was substantially lower, and the residual N r unaccounted for by the gas-phase measurements was 33.9±16 % of N r (panel c). The majority of sumN was represented by HCN and NH 3 , with acetonitrile (CH 3 CN) higher than any of the other inorganics, HNCO, NO 2 , or HONO. The NVOCs were also a larger fraction of N r −NO than in the case of Fire 047 shown above, a feature that implies that more semivolatile organic compounds (SVOCs) survive these types of fires and could make a proportionally higher contribution to the N r budget in this fire relative to Fire 047. Fire Lab results of particle organic carbon measurements (Jen et al., 2019) and field measurements in environments with a lot of dung burning (Jayarathne et al., 2018;Stockwell et al., 2016a) are consistent with a higher EF for particle organic carbon and by extension particle NVOC compounds. The quantity (N r −NO) / N r was relatively high and had less dynamic range than for fires with more flaming combustion like Fire 047.
An overall budget of N r can be constructed for all of the stack fires through integrating the time profile of all the compounds and compound classes. The fire-integrated measure-  ments of inorganic and NVOC species are listed in the Supplement as ratios to N r for each stack fire (Table S1). The summary of all the fire integrated X i /N r fractions (where X i is the N r species or quantity) is given in Table 3 for all the fires for which we have a complete set of measurements (43 fires). In general, NO was the major species, followed by NH 3 , and the other inorganic N r species, NO 2 , HNCO, HONO, and HCN had individual contributions of 4.3 % to 9.4 %. NVOC species were less than 5 % of N r on average. The unaccounted-for N r , defined as (N r − NO − sumN)/N r had a median value of 14.3 % and a mean (± SD) of 15 (±10) %. Overall, 85 % of N r was accounted for by the gasphase measurements. The distribution of whole-fire N r residuals is plotted as a histogram in Fig. 6. We expect the residual N r was composed of either semi-or low-volatility compounds or particle-bound N r compounds, which we know are converted efficiently by the N r catalyst  but not detected by the instruments included in this analysis. Along these lines, there is some indication that the residual has a systematic variation with whole-fire MCE, with higher residuals (up to 30 %) observed at lower MCEs and higher (N r −NO) / N r (see Fig. S3a, b), which would be consistent with higher EF for SVOC at low MCE (Jen et al., 2019) and particle N r having a higher contribution from NO − 3 (May et al., 2014) and, perhaps, particle ammonium or reduced-N r compounds. In general, there is more particulate organic material emitted from fires at low MCE (Jen et al., 2019), so we would expect more particle N at low MCE to go along with that.

Systematic dependences of N r composition on combustion processes
The features noted in fires shown above, as well as the anticorrelation of MCE and (N r −NO) / N r lead to the question of whether there are systematic dependences in N r -compound composition on fire stage that can be used to formally classify and/or potentially predict the relative emissions of N r compounds. MCE has been used as a rough indicator of the relative amounts of flaming and smoldering combustion in a fire, with high MCE (∼ 99 %) being "pure" flaming, low MCE (∼ 80 %) being "pure smoldering," and an MCE of ∼ 0.9 being roughly equal amounts of both (Sect. 2.1.1 in Akagi et al., 2011). It should be understood that "smoldering" in this framework is a lumped term that includes all nonflame processes such as pyrolysis, glowing, and distillation, which are the processes that produce gaseous fuel to support flaming (Yokelson et al., 1996). In addition, pure flaming is essentially the efficient oxidation of smoldering products before they enter the atmosphere. However, for MCE to predict flaming and smoldering N r species well, the variable fuel N must be considered. For instance, NO x is clearly produced by flaming based on its temporal profile, but fire-integrated EF NO x may not correlate with MCE due to variable fuel N. In these cases, EF NO x / fuel N or NH 3 / NO x may still correlate (or anticorrelate) well with MCE (e.g., Fig. 4 in Burling et al., 2010, or Yokelson et al., 1996. Finally, the flame chemistry involving NH 3 , HNCO, and HCN both produces and destroys NO in a fashion that does not conserve N r . This chemistry is explored in Fig. 7, in which NO x , NH 3 , HNCO, HCN, HONO, and CH 3 CN ratios to N r are plotted vs. real-time MCE for Fire 047 as a typical example for fires that have a substantial range of MCEs (e.g., from 0.8 to above 0.98). The relationship between NH 3 /N r and MCE confirms that NH 3 is primarily a smoldering emission, and NO x /N r increases with increasing MCE in a nonlinear fashion that confirms it is primarily a flaming compound. Such a nonlinear dependence has also been seen for other flamingrelated quantities such as elemental carbon/TC or EF HCl (Christian et al., 2003;Stockwell et al., 2014). Most importantly, the variations in HNCO/N r , HCN/N r , HONO/N r , and CH 3 CN/N r versus MCE do not arise dominantly from either regime, as these are species that are likely produced by multiple pathways (e.g., "incomplete flaming", pyrolysis, or possibly glowing). By incomplete flame chemistry we mean the production of incompletely oxidized products in flames such as the complex system of reactions shown in Fig. 1. These reactions involving HNCO, HCN, and NH 3 both produce and destroy NO, while HONO is produced from reactions of NO and NO 2 that are faster at slightly lower temperatures, for example the three-body association reaction of NO with OH radical (Manion et al., 2015). Variable turbulence in the turbulent diffusion flames that are characteristic of open BB likely contributes to varying temperatures and, therefore, varying amounts of incomplete oxidation of the fuel N (Shaddix et al., 1994).

The PMF analysis of coniferous fuels
The complexity of the dependence of N r speciation on combustion chemistry suggests that MCE is an insufficient model to use for applying lab results to real-world fire emissions Figure 8. (a) The measured N r signal for Fire 063 (lodgepole pine) (blue line), the sum of the signal reconstructed by the PMF (purple points), and the three PMF factors -combustion (grey), high temperature (green), and low temperature (red) -plotted in a stacked fashion (i.e., added on top of one another). (b) The "residual" of the PMF fit consisting of the measured N r signal minus the N r signal reconstructed by the PMF, as a percentage of the N r signal. (Stockwell et al., 2016a;Yokelson et al., 2013b). Accordingly, we employed the positive matrix factorization (PMF) method (see Methodology section) that has been used by a number of groups to probe the sources contributing to complex mixtures (see for example Ulbrich et al., 2009;Sekimoto et al., 2018). Our PMF results showed several general features, irrespective of the inclusion or exclusion of CO 2 , CO, and N r . The emissions were best fit by three factors (with approximate descriptive names justified below and prime species): (1) a combustion (flaming) factor (abbreviated Comb-N), (2) a high-temperature pyrolysis factor (HT-N), and (3) a low-temperature pyrolysis factor (LT-N). We use these terms in part to harmonize our discussion with the VOC results discussed by Sekimoto et al. (2018). An example time series for the PMF analysis of a coniferous fuel with just the N r species included is shown in Fig. 8 for a realistic mix of lodgepole pine (Fire 063). In this case, several different F peak values were tried (−1, 0, +1) and runs with 100 different seeds (initial factor profiles) were performed. The results of those analyses (Fig. S4) show that a three-factor PMF result is robust. A PMF analysis was performed on the consolidated time series of all coniferous fuels fit using just the N r species, as shown in Fig. S5. In this case F peak = 0 was used, and the Q/Q expected showed an inflection for the three-factor solution at a value of 5.3. The three factors successfully describe the majority of the N r emissions where the difference between the measured and calculated mass is on average 5.1 % for coniferous fuels and 4.6 % for chaparrals as indicated in Table 4. Several metrics of the PMF analysis quantify how the compounds or compound classes contribute to each factor. The "loadings" of the three different factors, i.e., the contribution of compounds to each factor, for coniferous fuels are shown in Fig. 9a, and the distribution of a given compound or compound class amongst the three factors is shown in Fig. 9b as normalized fraction. Normalized fraction is equal to the PMF-determined contribution of a compound to a factor, divided by the sum of the contribution of the compound to all three factors. The Comb-N factor contained NO, NO 2 , and HONO; the HT-N factor had mostly HCN, HNCO, and nitriles, with contributions from NO 2 and nitro compounds; and the LT-N factor contained NH 3 , amines, amides, and heterocyclics. Within the Comb-N factor there is some evidence that the relative amounts of HONO and NO x depend on fuel moisture. For example, the ratio HONO/NO x for whole fires shows some correlation with needle moisture in coniferous fires that were canopy fuels (foliage and small woody biomass), as shown in Fig. S6. This may be due to flame processes that interconvert NO x and HONO in the presence of water vapor or OH (see Fig. 1).
Literature values from studies where flame temperature was measured are typically in the range of 1100-1200 • C (Taylor et al., 2004;Wotton et al., 2012), so we would assume that would constitute the upper range of our Comb-N factor. The radical chemistry involving HCN, HNCO, and NH 3 starts to shut down below about 800-900 • C, according to the modeling of Glarborg et al. (2018), so we set 800 • C as a lower limit for the Comb-N factor. The HT-N factor species are known to be produced by the intense pyrolysis of fuel N r compounds (Hansson et al., 2004;Liu et al., 2018;Ren et al., 2010), which for these compounds becomes important at temperatures around 500-600 • C. Accordingly, we estimate the temperature range for the HT-N factor at 500-800 • C. The remaining LT-N factor results from mild pyrolysis and pertains to fire conditions of roughly 500 • C and below, and it was dominated by NH 3 , amines, amides, and some of the more complex organics . The names and temperature ranges are approximate and likely include processes that occur inside flames as part of the flame proper, as turbulent diffusive flames are highly variable in space and time.

The comparison of N-PMF factors to other fire parameters and VOC emission factors
It is useful to explore the correlation of N-PMF factors with other fire indicators to determine relationships for parameterizing N r emissions together with carbon and VOC emissions, in order to simplify how emissions might be parameterized in models. The Comb-N factor for coniferous fuels, which consisted of NO x and HONO, would be expected to correlate with CO 2 but not as well with MCE since the latter includes an indicator of incomplete combustion. The time series of Comb-N along with CO 2 and with MCE for Fire 037 (ponderosa pine) are plotted in Fig. 10. As expected they show an excellent correlation of Comb-N with CO 2 (R 2 = 0.942), since all the species are flaming compounds, but nonlinear correlation of Comb-N with MCE (R 2 = 0.363), since the latter factors in a smoldering compound (CO), similar to the NO x /N r vs. MCE plot for Fire 047 in Fig. 7. The excellent correlation of Comb-N with CO 2 is a broadly applicable result; the R 2 parameters for all the fires shown in Fig. S5 had an average of 0.898 and ranged from 0.806 to 0.966. As a consequence, we can conclude that CO 2 would be the best tracer for Comb-N in many western US ecosystems where conifers predominate, provided ambient CO 2 backgrounds can be properly accounted for as described by Yokelson et al. (2013a).
Our Comb-N factor did not correspond to the hightemperature VOC factor (HT-VOC) found by  in their pyrolysis study, because our broader study includes flaming combustion, which produces NO x and HONO, and almost none of the compounds classified as VOCs survive flaming conditions. However, the HT-N and HT-VOC factors are both linked to pyrolysis and were well correlated for many fires. An example of this is shown in Fig. 11 for Fire 037, a sample that was broadly representative of ponderosa pine (i.e., canopy and litter). This result can be rationalized by the fact that while HT-VOC factors have large contributions from many more compounds than the N compounds measured here, they also have large contributions (> 85 %) from HCN, HNCO, and HONO (in other words > 85 % of HCN, HNCO, and HONO are found in the HT-VOC factor). Since the HT-N factors are also heavily weighted by HCN and HNCO, it is reassuring that both of these PMF analyses have independently identified these species as important contributors to the HT fire regime. The R 2 correlation coefficients between the HT-N and HT-VOC factors for the coniferous fires shown in Fig. S5 averaged 0.866 and ranged from 0.419 to 0.959. As a consequence of this correlation, we can conclude that HCN is the best marker for the HT-N and HT-VOC factors in most western US wildfires, since HCN is essentially inert on the timescales of fire plumes (Li et al., 2000). It should be noted that other nitriles, particularly acetonitrile, also show up in the HT-N factor, and acetonitrile has also been used as a tracer of biomass combustion. However, it has been shown that this acetonitrile signal can be obscured in urban or industrial areas by solvent usage or can be quite small in woodstove emissions due to low N in the fuel (Coggon et al., 2016).
The correlations of LT-N and LT-VOC factors were not particularly high for most of the coniferous fires shown in Fig. S5. The average R 2 was 0.427 with a range of between 0.072 and 0.827. The reasons for this lack of correlation are not clear, as NH 3 , amines, and amides appear predominantly in both LT factors, and the absolute concentrations of NH 3 are usually quite high in these fires relative to VOCs . However, the LT-VOC factor includes many more compounds with a variety of functional groups not found in the LT-N factor, so it appears that the VOC and N compounds have sufficiently different pyrolysis chemistry that the LT factors do not show much correlation. We conclude that NH 3 (and particle NH + 4 ) will be the best marker for the LT-N factor in western US coniferous wildfires, but the LT-VOC chemistry might not be captured reliably by this marker.

PMF analysis of chaparral fuels
Chaparral is an important ecosystem of concern in wildfires that occur in central and southern California and other areas of the southwestern US. The emissions from burning chaparral fuels (manzanita and chamise) collected at two sites in California were also analyzed as a group and yielded three separate factors in a fashion similar to the coniferous fuels (see Fig. S7 for the PMF timeline). As with the coniferous fuels, there was essentially no change in the three-factor solution with F peak , so F peak = 0 was used, and the Q/Q exp was 3.8. The chaparral factors had slightly different composition (Fig. S8); the combustion factor was mostly NO, with small amounts of HNCO, HONO, and NH 3 ; the high-temperature factor was dominated by NO 2 and included HONO, HCN, and HNCO; and the low-temperature factor was mostly NH 3 with a slight amount of NO contributing. The NVOC species were found in both the medium-and low-temperature factors.
There was less similarity between the Comb-N factor and CO 2 emissions for chaparral fuels compared to those found for coniferous fuels, with an average correlation coefficient (R 2 ) of 0.689, with a range from 0.244 and 0.950. As a result, there may not be a simple conserved tracer for the combustion factor of these fuel types; however total odd nitrogen (NO y ) which is NO x and all the compounds that are produced from NO x in the troposphere, may be useful as it is a reasonably conserved tracer in the absence of wet or dry deposition of particles. Correlation coefficients between the HT-N and HT-VOC factors were on average R 2 = 0.551, with a range 0.047-0.911. The correlations between LT-N and LT-VOC factors were in the same range for chaparral fuels as for coniferous, with average R 2 = 0.447, range 0.028-0.827.
There were some fuels that do not sustain flaming combustion well, specifically duff, Yak dung, and Indonesian peat. These fires exhibited little or no NO emission commensurate with minimal flaming combustion. Instead the emissions were mostly the pyrolysis products NH 3 , (0.22-0.53 N r fraction), and HCN (up to 0.32 N r fraction for peat). It was also apparent that these fires also had unaccounted-for N r , close to or just over 0.30 (Table S1). The distribution of N r compounds in the one peat fire that we measured (Fire 055) is in line with those reported for fires measured in situ, which showed relatively high EFs for HCN and NH 3 (Stockwell et al., 2016b(Stockwell et al., , 2015.

Application to real-world fires
The application of our N r emissions results to real-world fires will depend somewhat on the nature of the information available on a particular fire or fire complex. As a good starting point, or in the absence of detailed N and C analyses of fuels, a N r /C ratio of 0.37 % appears to capture most of the fires studied in this work. The N r could be apportioned according to the results summarized in Table 3. Adjustments to those fractions can be made by scaling slightly by average MCE, with the knowledge that intermediate species (HT-N pyrolysis species) such as HCN and HNCO do not scale in the simple manner that NH 3 and NO x do. If measurements of marker compounds are available then CO 2 , HCN, and the sum NH 3 + NH + 4 can be used for the combustion, Figure 11. Details of the PMF factors for Fire 037 (ponderosa pine). Panel (a) shows the total N r signal (magenta) and the Comb-N factor (black); panel (b) shows the HT-N factor (green) and HT-VOC factor (blue), and panel (c) shows the LT-N factor (red) and LT-VOC factor (orange). The insets (d) and (e) show the correlation of the two HT factors and the correlation between the two LT factors, respectively.
high-temperature pyrolysis, and low-temperature pyrolysis factors, respectively.

Conclusions
Seventy-five stack fire experiments were conducted during the FIREX Fire Lab experiments in Fall 2016. A range of fuels characteristic of the western US was burned under conditions and in mixtures meant to represent authentic wildfire conditions, as closely as is possible in the laboratory. Total reactive nitrogen (N r : all N-containing compounds except N 2 and N 2 O) was measured along with a suite of N-containing compounds in order to obtain a budget for N r emissions and to examine relationships between fuels, combustion conditions, and emissions chemistry. Natural convection wildfires do not burn hot enough to produce NO x from N 2 and O 2 , so all N r emissions come from the fuel N. Almost all of the fires representative of North American ecosystems had emissions that clustered around a N r /C ratio of 0.37 %, which can serve as a starting point for scaling emissions from these ecosystems. Comparing total N r and total carbon emissions with the N/C ratios of both the original fuel and remaining ash allowed us to estimate that an average of 68 % (±14 %) of the fuel nitrogen ends up as N 2 and N 2 O. This loss of nitrogen can be used to estimate how much fuel nitrogen ends up as N r . Of the remaining N emitted as N r , approximately 85 % (±10 %) was accounted for by individually measured gas-phase species, while the rest was most likely particle-bound NH + 4 and NO − 3 , with a smaller contribution from low-volatility species or other species such as cyanogen (Lobert and Warnatz, 1993), that were not quantified by the instruments for individual measurements we used in this study. The speciation and modeling of N r we present promotes accurate modeling of fire plume chemistry since the photochemistry of many fire plumes is NO x -limited, and NH 3 is an important contributor to particle chemistry.
The individual N r species composition normalized to total N r , to account for fuel N variability, correlated monotonically with flaming versus vs. smoldering combustion as indicated by modified combustion efficiency (MCE) for some species (e.g., NH 3 , NO x ). Other species, such as HCN and HNCO, peaked at intermediate MCE values. Positive matrix factorization (PMF) showed that all the measured N r emissions from the main two categories of fuels, conifers, and chaparral, grouped into three mixtures (factors), roughly attributed to temperature: combustion (NO x , HONO), high temperature (HNCO, HCN, nitriles), and low-temperature (NH 3 , amines, amides). Chemical kinetic and pyrolysis considerations set the temperature ranges for these regimes at approximately 800-1200, 500-800, and < 500 • C, respectively.
This paper connects mechanistic aspects of N combustion chemistry to the budget of N r emissions from biomass burning. The emission composition measurements detailed here give useful information concerning what the initial conditions will be in actual fire plumes. These results suggest that for coniferous fuels characteristic of the western US, CO 2 is the best marker for flaming combustion, HCN is the best marker for high-temperature pyrolysis processes, and NH 3 /NH + 4 is the best marker for low-temperature pyrolysis processes. The HT-N and HT-VOC pyrolysis factors showed a high degree of correlation especially for coniferous fuels, which can simplify how these different classes of emissions can be estimated. Results from less comprehensive field experiments can be combined with this emissions information to improve the representation of N r chemistry in the modeling frameworks needed to predict fire plume chemistry and impacts.
Author contributions. JMR, RJY, CW, and JdG designed the research. The measurements were conducted by JMR, CES, CW, RJY, JdG, YL, VS, ARK, KS, MMC, BY, KJZ, SSB, CS, and SHD. All authors contributed to the discussion and interpretation of the results and writing the paper.
Competing interests. Joost de Gouw worked as a consultant for Aerodyne Research during part of the preparation phase of this paper.
Disclaimer. Any mention of brand names or manufacturers is for information purposes only and does not constitute an endorsement.